|
ATLAS Offline Software
|
Algorithm to perform truth matching on secondary vertices.
More...
#include <SecVertexTruthMatchAlg.h>
|
| SecVertexTruthMatchAlg (const std::string &name, ISvcLocator *svcLoc) |
| Regular Algorithm constructor. More...
|
|
virtual StatusCode | initialize () override |
|
virtual StatusCode | execute () override |
|
::StatusCode | requestFileExecute () |
| register this algorithm to have an implementation of fileexecute More...
|
|
::StatusCode | requestBeginInputFile () |
| register this algorithm to have an implementation of beginInputFile More...
|
|
::StatusCode | requestEndInputFile () |
| register this algorithm to have an implementation of endInputFile More...
|
|
void | handle (const Incident &inc) |
| receive the given incident More...
|
|
virtual StatusCode | sysInitialize () |
| Initialization method invoked by the framework. More...
|
|
const ServiceHandle< ITHistSvc > & | histSvc () const |
| The standard THistSvc (for writing histograms and TTrees and more to a root file) Returns (kind of) a pointer to the THistSvc . More...
|
|
virtual const DataObjIDColl & | extraOutputDeps () const override |
| Return the list of extra output dependencies. More...
|
|
ServiceHandle< StoreGateSvc > & | evtStore () |
| The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc . More...
|
|
const ServiceHandle< StoreGateSvc > & | evtStore () const |
| The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc . More...
|
|
const ServiceHandle< StoreGateSvc > & | detStore () const |
| The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc . More...
|
|
virtual StatusCode | sysStart () override |
| Handle START transition. More...
|
|
virtual std::vector< Gaudi::DataHandle * > | inputHandles () const override |
| Return this algorithm's input handles. More...
|
|
virtual std::vector< Gaudi::DataHandle * > | outputHandles () const override |
| Return this algorithm's output handles. More...
|
|
Gaudi::Details::PropertyBase & | declareProperty (Gaudi::Property< T > &t) |
|
Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &) |
| Declare a new Gaudi property. More...
|
|
Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &) |
| Declare a new Gaudi property. More...
|
|
Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &) |
|
Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &) |
| Declare a new Gaudi property. More...
|
|
Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, T &property, const std::string &doc="none") |
| Declare a new Gaudi property. More...
|
|
void | updateVHKA (Gaudi::Details::PropertyBase &) |
|
MsgStream & | msg () const |
|
MsgStream & | msg (const MSG::Level lvl) const |
|
bool | msgLvl (const MSG::Level lvl) const |
|
|
virtual ::StatusCode | finalize () |
| finalize this algorithm More...
|
|
virtual void | print () const |
| print the state of the algorithm More...
|
|
virtual ::StatusCode | fileExecute () |
| perform the action exactly once for each file in the dataset More...
|
|
virtual ::StatusCode | beginInputFile () |
| perform the action for the beginning of an input file More...
|
|
virtual ::StatusCode | endInputFile () |
| perform the action for the end of an input file More...
|
|
void | renounceArray (SG::VarHandleKeyArray &handlesArray) |
| remove all handles from I/O resolution More...
|
|
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > | renounce (T &h) |
|
void | extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps) |
| Add StoreName to extra input/output deps as needed. More...
|
|
StatusCode | configAthHistogramming (const ServiceHandle< ITHistSvc > &histSvc, const std::string &prefix, const std::string &rootDir, const std::string &histNamePrefix, const std::string &histNamePostfix, const std::string &histTitlePrefix, const std::string &histTitlePostfix) |
| To be called by the derived classes to fill the internal configuration. More...
|
|
TH1 * | bookGetPointer (const TH1 &hist, const std::string &tDir="", const std::string &stream="") |
| Simplify the booking and registering (into THistSvc) of histograms. More...
|
|
TH1 * | bookGetPointer (TH1 *hist, const std::string &tDir="", const std::string &stream="") |
| Simplify the booking and registering (into THistSvc) of histograms. More...
|
|
TH1 * | bookGetPointer (TH1 &histRef, std::string tDir="", std::string stream="") |
| Simplify the booking and registering (into THistSvc) of histograms. More...
|
|
TTree * | bookGetPointer (const TTree &treeRef, std::string tDir="", std::string stream="") |
| Simplify the booking and registering (into THistSvc) of TTrees. More...
|
|
TGraph * | bookGetPointer (const TGraph &graphRef, std::string tDir="", std::string stream="") |
| Simplify the booking and registering (into THistSvc) of TGraphs. More...
|
|
TEfficiency * | bookGetPointer (const TEfficiency &eff, const std::string &tDir="", const std::string &stream="") |
| Simplify the booking and registering (into THistSvc) of TEfficiency. More...
|
|
TEfficiency * | bookGetPointer (TEfficiency *eff, const std::string &tDir="", const std::string &stream="") |
| Simplify the booking and registering (into THistSvc) of TEfficiency. More...
|
|
TEfficiency * | bookGetPointer (TEfficiency &effRef, std::string tDir="", std::string stream="") |
| Simplify the booking and registering (into THistSvc) of TEfficiency. More...
|
|
StatusCode | book (const TH1 &hist, const std::string &tDir="", const std::string &stream="") |
| Simplify the booking and registering (into THistSvc) of histograms. More...
|
|
StatusCode | book (TH1 *hist, const std::string &tDir="", const std::string &stream="") |
| Simplify the booking and registering (into THistSvc) of histograms. More...
|
|
StatusCode | book (TH1 &histRef, const std::string &tDir="", const std::string &stream="") |
| Simplify the booking and registering (into THistSvc) of histograms. More...
|
|
StatusCode | book (const TTree &treeRef, const std::string &tDir="", const std::string &stream="") |
| Simplify the booking and registering (into THistSvc) of TTrees. More...
|
|
StatusCode | book (const TGraph &graphRef, const std::string &tDir="", const std::string &stream="") |
| Simplify the booking and registering (into THistSvc) of TGraphs. More...
|
|
StatusCode | book (const TEfficiency &eff, const std::string &tDir="", const std::string &stream="") |
| Simplify the booking and registering (into THistSvc) of TEfficiency. More...
|
|
StatusCode | book (TEfficiency *eff, const std::string &tDir="", const std::string &stream="") |
| Simplify the booking and registering (into THistSvc) of TEfficiency. More...
|
|
StatusCode | book (TEfficiency &effRef, const std::string &tDir="", const std::string &stream="") |
| Simplify the booking and registering (into THistSvc) of TEfficiency. More...
|
|
TH1 * | hist (const std::string &histName, const std::string &tDir="", const std::string &stream="") |
| Simplify the retrieval of registered histograms of any type. More...
|
|
TH2 * | hist2d (const std::string &histName, const std::string &tDir="", const std::string &stream="") |
| Simplify the retrieval of registered 2-d histograms. More...
|
|
TH3 * | hist3d (const std::string &histName, const std::string &tDir="", const std::string &stream="") |
| Simplify the retrieval of registered 3-d histograms. More...
|
|
TTree * | tree (const std::string &treeName, const std::string &tDir="", const std::string &stream="") |
| Simplify the retrieval of registered TTrees. More...
|
|
TGraph * | graph (const std::string &graphName, const std::string &tDir="", const std::string &stream="") |
| Simplify the retrieval of registered TGraphs. More...
|
|
TEfficiency * | efficiency (const std::string &effName, const std::string &tDir="", const std::string &stream="") |
| Simplify the retrieval of registered TEfficiency. More...
|
|
|
void | fillRecoHistograms (const xAOD::Vertex *secVtx, const std::string &matchType) |
|
void | fillTruthHistograms (const xAOD::TruthVertex *truthVtx, const std::string &truthType) |
|
Gaudi::Details::PropertyBase & | declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &) |
| specialization for handling Gaudi::Property<SG::VarHandleKey> More...
|
|
Gaudi::Details::PropertyBase & | declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &) |
| specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
|
|
Gaudi::Details::PropertyBase & | declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &) |
| specialization for handling Gaudi::Property<SG::VarHandleBase> More...
|
|
Gaudi::Details::PropertyBase & | declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &) |
| specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...
|
|
void | buildBookingString (std::string &bookingString, std::string &histName, std::string &tDir, std::string &stream, bool usePrefixPostfix=false) |
| Method to build individual booking string. More...
|
|
void | myReplace (std::string &str, const std::string &oldStr, const std::string &newStr) |
| Helper method to replace sub-string. More...
|
|
hash_t | hash (const std::string &histName) const |
| Method to calculate a 32-bit hash from a string. More...
|
|
Algorithm to perform truth matching on secondary vertices.
- Author
- Jackson Burzynski jacks.nosp@m.on.c.nosp@m.arl.b.nosp@m.urzy.nosp@m.nski@.nosp@m.cern.nosp@m..ch
Definition at line 31 of file SecVertexTruthMatchAlg.h.
◆ ConstMetaStorePtr_t
◆ EffMap_t
◆ GraphMap_t
◆ hash_t
◆ HistMap_t
◆ MetaStore_t
◆ MetaStorePtr_t
Type of the metadata store pointer in standalone mode.
Definition at line 110 of file AnaAlgorithm.h.
◆ StoreGateSvc_t
◆ TreeMap_t
◆ SecVertexTruthMatchAlg()
CP::SecVertexTruthMatchAlg::SecVertexTruthMatchAlg |
( |
const std::string & |
name, |
|
|
ISvcLocator * |
svcLoc |
|
) |
| |
◆ beginInputFile()
StatusCode EL::AnaAlgorithm::beginInputFile |
( |
| ) |
|
|
protectedinherited |
perform the action for the beginning of an input file
Ideally you don't use this, but instead rely on meta-data tools instead. However, there are enough people asking for it that I decided to implement it anyways.
\warn To use this you have to call requestBeginInputFile to use this.
\warn If a file is split across multiple jobs this will be called more than once. This only happens for specific batch drivers and/or if it is explicitly configured by the user. With PROOF it could even happen multiple times within the same job, and while PROOF is no longer supported that behavior may come back if support for a similar framework is added in the future. As such, this method should not be used for accounting that relies to be called exactly once per file, take a look at fileExecute if you want something that is guaranteed to be executed exactly once per input file.
\warn The execution order of beginInputFile and fileExecute is currently unspecified.
Definition at line 350 of file AnaAlgorithm.cxx.
353 return StatusCode::SUCCESS;
◆ book() [1/8]
StatusCode AthHistogramming::book |
( |
const TEfficiency & |
eff, |
|
|
const std::string & |
tDir = "" , |
|
|
const std::string & |
stream = "" |
|
) |
| |
|
inlineprotectedinherited |
Simplify the booking and registering (into THistSvc) of TEfficiency.
Definition at line 335 of file AthHistogramming.h.
338 TEfficiency* effClone =
dynamic_cast< TEfficiency*
>(
eff.Clone() );
340 m_msg << MSG::ERROR <<
"Couldn't create a TEfficiency clone" <<
endmsg;
341 return StatusCode::FAILURE;
◆ book() [2/8]
StatusCode AthHistogramming::book |
( |
const TGraph & |
graphRef, |
|
|
const std::string & |
tDir = "" , |
|
|
const std::string & |
stream = "" |
|
) |
| |
|
inlineprotectedinherited |
Simplify the booking and registering (into THistSvc) of TGraphs.
◆ book() [3/8]
StatusCode AthHistogramming::book |
( |
const TH1 & |
hist, |
|
|
const std::string & |
tDir = "" , |
|
|
const std::string & |
stream = "" |
|
) |
| |
|
inlineprotectedinherited |
Simplify the booking and registering (into THistSvc) of histograms.
Definition at line 303 of file AthHistogramming.h.
306 TH1* histClone =
dynamic_cast< TH1*
>(
hist.Clone() );
308 m_msg << MSG::ERROR <<
"Couldn't create a TH1 clone" <<
endmsg;
309 return StatusCode::FAILURE;
311 return this->
book( *histClone, tDir,
stream );
◆ book() [4/8]
StatusCode AthHistogramming::book |
( |
const TTree & |
treeRef, |
|
|
const std::string & |
tDir = "" , |
|
|
const std::string & |
stream = "" |
|
) |
| |
|
inlineprotectedinherited |
Simplify the booking and registering (into THistSvc) of TTrees.
Definition at line 403 of file AthHistogramming.h.
409 return StatusCode::SUCCESS;
413 return StatusCode::FAILURE;
◆ book() [5/8]
StatusCode AthHistogramming::book |
( |
TEfficiency & |
effRef, |
|
|
const std::string & |
tDir = "" , |
|
|
const std::string & |
stream = "" |
|
) |
| |
|
inlineprotectedinherited |
Simplify the booking and registering (into THistSvc) of TEfficiency.
Definition at line 356 of file AthHistogramming.h.
361 m_msg << MSG::ERROR <<
"Couldn't book a TEfficiency" <<
endmsg;
362 return StatusCode::FAILURE;
364 return StatusCode::SUCCESS;
◆ book() [6/8]
StatusCode AthHistogramming::book |
( |
TEfficiency * |
eff, |
|
|
const std::string & |
tDir = "" , |
|
|
const std::string & |
stream = "" |
|
) |
| |
|
inlineprotectedinherited |
Simplify the booking and registering (into THistSvc) of TEfficiency.
Definition at line 346 of file AthHistogramming.h.
349 m_msg << MSG::ERROR <<
"Got a zero pointer to a TEfficiency" <<
endmsg;
350 return StatusCode::FAILURE;
◆ book() [7/8]
StatusCode AthHistogramming::book |
( |
TH1 & |
histRef, |
|
|
const std::string & |
tDir = "" , |
|
|
const std::string & |
stream = "" |
|
) |
| |
|
inlineprotectedinherited |
Simplify the booking and registering (into THistSvc) of histograms.
Definition at line 324 of file AthHistogramming.h.
328 if ( !histPointer ) {
329 m_msg << MSG::ERROR <<
"Couldn't book a TH1" <<
endmsg;
330 return StatusCode::FAILURE;
332 return StatusCode::SUCCESS;
◆ book() [8/8]
StatusCode AthHistogramming::book |
( |
TH1 * |
hist, |
|
|
const std::string & |
tDir = "" , |
|
|
const std::string & |
stream = "" |
|
) |
| |
|
inlineprotectedinherited |
Simplify the booking and registering (into THistSvc) of histograms.
Definition at line 314 of file AthHistogramming.h.
317 m_msg << MSG::ERROR <<
"Got a zero pointer to a TH1" <<
endmsg;
318 return StatusCode::FAILURE;
◆ bookGetPointer() [1/8]
TEfficiency * AthHistogramming::bookGetPointer |
( |
const TEfficiency & |
eff, |
|
|
const std::string & |
tDir = "" , |
|
|
const std::string & |
stream = "" |
|
) |
| |
|
inlineprotectedinherited |
Simplify the booking and registering (into THistSvc) of TEfficiency.
Definition at line 281 of file AthHistogramming.h.
284 TEfficiency* histClone =
dynamic_cast< TEfficiency*
>(
hist.Clone() );
286 m_msg << MSG::ERROR <<
"Couldn't create a TEfficiency clone in bookGetPointer" <<
endmsg;
◆ bookGetPointer() [2/8]
TGraph * AthHistogramming::bookGetPointer |
( |
const TGraph & |
graphRef, |
|
|
std::string |
tDir = "" , |
|
|
std::string |
stream = "" |
|
) |
| |
|
protectedinherited |
Simplify the booking and registering (into THistSvc) of TGraphs.
Definition at line 427 of file AthHistogramming.cxx.
430 const TGraph* graphPointer = &graphRef;
435 m_msg << MSG::WARNING
436 <<
"We got an invalid TGraph pointer in the BookGetPointer(TGraph*) method of the class" <<
m_name
442 std::string graphName = graphPointer->GetName();
443 const std::string graphTitle = graphPointer->GetTitle();
446 const hash_t graphHash = this->
hash(graphName);
447 GraphMap_t::const_iterator
it =
m_graphMap.find( graphHash );
450 m_msg << MSG::WARNING
451 <<
"Detected a hash collision. The hash for the TGraph with name=" << graphName
452 <<
" already exists and points to a TGraph with name=" <<
it->second->GetName()
453 <<
" NOT going to book the new histogram and returning a NULL pointer!" <<
endmsg;
461 m_msg << MSG::WARNING
462 <<
"We couldn't clone the TGraph in the BookGetPointer(TGraph&) method of the class" <<
m_name
469 std::string bookingString(
"");
473 if ( !((
histSvc()->regGraph(bookingString, graphClone)).isSuccess()) )
475 m_msg << MSG::WARNING
476 <<
"Problem registering TGraph with name " << graphName
477 <<
", title " << graphTitle
483 m_graphMap.insert(
m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphClone ) );
◆ bookGetPointer() [3/8]
TH1 * AthHistogramming::bookGetPointer |
( |
const TH1 & |
hist, |
|
|
const std::string & |
tDir = "" , |
|
|
const std::string & |
stream = "" |
|
) |
| |
|
inlineprotectedinherited |
Simplify the booking and registering (into THistSvc) of histograms.
Definition at line 260 of file AthHistogramming.h.
263 TH1* histClone =
dynamic_cast< TH1*
>(
hist.Clone() );
265 m_msg << MSG::ERROR <<
"Couldn't create a TH1 clone in bookGetPointer" <<
endmsg;
◆ bookGetPointer() [4/8]
TTree * AthHistogramming::bookGetPointer |
( |
const TTree & |
treeRef, |
|
|
std::string |
tDir = "" , |
|
|
std::string |
stream = "" |
|
) |
| |
|
protectedinherited |
Simplify the booking and registering (into THistSvc) of TTrees.
Definition at line 312 of file AthHistogramming.cxx.
315 const TTree* treePointer = &treeRef;
320 m_msg << MSG::WARNING
321 <<
"We got an invalid TTree pointer in the BookGetPointer(TTree*) method of the class" <<
m_name
327 std::string
treeName = treePointer->GetName();
328 const std::string treeTitle = treePointer->GetTitle();
332 TreeMap_t::const_iterator
it =
m_treeMap.find( treeHash );
335 m_msg << MSG::WARNING
336 <<
"Detected a hash collision. The hash for the TTree with name=" <<
treeName
337 <<
" already exists and points to a TTree with name=" <<
it->second->GetName()
338 <<
" NOT going to book the new histogram and returning a NULL pointer!" <<
endmsg;
343 TTree* treeClone =
dynamic_cast< TTree*
>( treePointer->Clone(
treeName.c_str()) );
346 m_msg << MSG::WARNING
347 <<
"We couldn't clone the TTree in the BookGetPointer(TTree&) method of the class" <<
m_name
351 treeClone->SetTitle (treeTitle.c_str());
354 std::string bookingString(
"");
358 if ( !((
histSvc()->regTree(bookingString, treeClone)).isSuccess()) )
360 m_msg << MSG::WARNING
361 <<
"Problem registering TTree with name " <<
treeName
362 <<
", title " << treeTitle
368 m_treeMap.insert(
m_treeMap.end(), std::pair< const hash_t, TTree* >( treeHash, treeClone ) );
◆ bookGetPointer() [5/8]
TEfficiency * AthHistogramming::bookGetPointer |
( |
TEfficiency & |
effRef, |
|
|
std::string |
tDir = "" , |
|
|
std::string |
stream = "" |
|
) |
| |
|
protectedinherited |
Simplify the booking and registering (into THistSvc) of TEfficiency.
Definition at line 146 of file AthHistogramming.cxx.
149 std::string effName(effRef.GetName());
150 const std::string effTitle(effRef.GetTitle());
151 std::string bookingString(
"");
155 effRef.SetName(effName.c_str());
159 EffMap_t::const_iterator
it =
m_effMap.find( effHash );
162 m_msg << MSG::WARNING
163 <<
"Detected a hash collision. The hash for the TEfficiency with name=" << effName
164 <<
" already exists and points to a TEfficiency with name=" <<
it->second->GetName()
165 <<
" NOT going to book the new TEfficiency and returning a NULL pointer!" <<
endmsg;
175 if ( !((
histSvc()->regEfficiency(bookingString, &effRef)).isSuccess()) )
177 m_msg << MSG::WARNING
178 <<
"Problem registering TEfficiency with name " << effName
180 <<
", title " << effTitle
188 m_effMap.insert(
m_effMap.end(), std::pair< const hash_t, TEfficiency* >( effHash, &effRef ) );
◆ bookGetPointer() [6/8]
TEfficiency * AthHistogramming::bookGetPointer |
( |
TEfficiency * |
eff, |
|
|
const std::string & |
tDir = "" , |
|
|
const std::string & |
stream = "" |
|
) |
| |
|
inlineprotectedinherited |
Simplify the booking and registering (into THistSvc) of TEfficiency.
Definition at line 293 of file AthHistogramming.h.
296 m_msg << MSG::ERROR <<
"Got a zero pointer to a TEfficiency in bookGetPointer" <<
endmsg;
◆ bookGetPointer() [7/8]
TH1 * AthHistogramming::bookGetPointer |
( |
TH1 & |
histRef, |
|
|
std::string |
tDir = "" , |
|
|
std::string |
stream = "" |
|
) |
| |
|
protectedinherited |
Simplify the booking and registering (into THistSvc) of histograms.
Definition at line 98 of file AthHistogramming.cxx.
101 std::string
histName(histRef.GetName());
102 const std::string histTitle(histRef.GetTitle());
103 std::string bookingString(
"");
111 HistMap_t::const_iterator
it =
m_histMap.find( histHash );
114 m_msg << MSG::WARNING
115 <<
"Detected a hash collision. The hash for the histogram with name=" <<
histName
116 <<
" already exists and points to a histogram with name=" <<
it->second->GetName()
117 <<
" NOT going to book the new histogram and returning a NULL pointer!" <<
endmsg;
127 if ( !((
histSvc()->regHist(bookingString, &histRef)).isSuccess()) )
129 m_msg << MSG::WARNING
130 <<
"Problem registering histogram with name " <<
histName
132 <<
", title " << histTitle
140 m_histMap.insert(
m_histMap.end(), std::pair< const hash_t, TH1* >( histHash, &histRef ) );
◆ bookGetPointer() [8/8]
TH1 * AthHistogramming::bookGetPointer |
( |
TH1 * |
hist, |
|
|
const std::string & |
tDir = "" , |
|
|
const std::string & |
stream = "" |
|
) |
| |
|
inlineprotectedinherited |
Simplify the booking and registering (into THistSvc) of histograms.
Definition at line 272 of file AthHistogramming.h.
275 m_msg << MSG::ERROR <<
"Got a zero pointer to a TH1 in bookGetPointer" <<
endmsg;
◆ buildBookingString()
void AthHistogramming::buildBookingString |
( |
std::string & |
bookingString, |
|
|
std::string & |
histName, |
|
|
std::string & |
tDir, |
|
|
std::string & |
stream, |
|
|
bool |
usePrefixPostfix = false |
|
) |
| |
|
privateinherited |
Method to build individual booking string.
Definition at line 560 of file AthHistogramming.cxx.
569 if(
pos != std::string::npos){
576 if(usePrefixPostfix){
581 while(bookingString.find(
"//") != std::string::npos){
◆ configAthHistogramming()
StatusCode AthHistogramming::configAthHistogramming |
( |
const ServiceHandle< ITHistSvc > & |
histSvc, |
|
|
const std::string & |
prefix, |
|
|
const std::string & |
rootDir, |
|
|
const std::string & |
histNamePrefix, |
|
|
const std::string & |
histNamePostfix, |
|
|
const std::string & |
histTitlePrefix, |
|
|
const std::string & |
histTitlePostfix |
|
) |
| |
|
protectedinherited |
To be called by the derived classes to fill the internal configuration.
Definition at line 66 of file AthHistogramming.cxx.
79 return StatusCode::SUCCESS;
◆ declareGaudiProperty() [1/4]
specialization for handling Gaudi::Property<SG::VarHandleKeyArray>
Definition at line 170 of file AthCommonDataStore.h.
175 hndl.documentation());
◆ declareGaudiProperty() [2/4]
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition at line 156 of file AthCommonDataStore.h.
161 hndl.documentation());
◆ declareGaudiProperty() [3/4]
specialization for handling Gaudi::Property<SG::VarHandleBase>
Definition at line 184 of file AthCommonDataStore.h.
189 hndl.documentation());
◆ declareGaudiProperty() [4/4]
◆ declareProperty() [1/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
hndl | Object holding the property value. |
doc | Documentation string for the property. |
This is the version for types that derive from SG::VarHandleBase
. The property value object is put on the input and output lists as appropriate; then we forward to the base class.
Definition at line 245 of file AthCommonDataStore.h.
250 this->declare(hndl.
vhKey());
251 hndl.
vhKey().setOwner(
this);
253 return PBASE::declareProperty(
name,hndl,
doc);
◆ declareProperty() [2/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
hndl | Object holding the property value. |
doc | Documentation string for the property. |
This is the version for types that derive from SG::VarHandleKey
. The property value object is put on the input and output lists as appropriate; then we forward to the base class.
Definition at line 221 of file AthCommonDataStore.h.
229 return PBASE::declareProperty(
name,hndl,
doc);
◆ declareProperty() [3/6]
◆ declareProperty() [4/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
property | Object holding the property value. |
doc | Documentation string for the property. |
This is the generic version, for types that do not derive from SG::VarHandleKey
. It just forwards to the base class version of declareProperty
.
Definition at line 333 of file AthCommonDataStore.h.
338 return PBASE::declareProperty(
name, property,
doc);
◆ declareProperty() [5/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
property | Object holding the property value. |
doc | Documentation string for the property. |
This dispatches to either the generic declareProperty
or the one for VarHandle/Key/KeyArray.
Definition at line 352 of file AthCommonDataStore.h.
◆ declareProperty() [6/6]
◆ detStore()
◆ efficiency()
TEfficiency * AthHistogramming::efficiency |
( |
const std::string & |
effName, |
|
|
const std::string & |
tDir = "" , |
|
|
const std::string & |
stream = "" |
|
) |
| |
|
protectedinherited |
Simplify the retrieval of registered TEfficiency.
Definition at line 250 of file AthHistogramming.cxx.
256 EffMap_t::const_iterator
it =
m_effMap.find( effHash );
261 std::string effNameCopy = effName;
262 std::string tDirCopy = tDir;
263 std::string streamCopy =
stream;
266 std::string bookingString(
"");
269 TEfficiency* effPointer(NULL);
270 if ( !((
histSvc()->getEfficiency(bookingString, effPointer)).isSuccess()) )
273 std::string bookingString(
"");
276 if ( !((
histSvc()->getEfficiency(bookingString, effPointer)).isSuccess()) )
278 m_msg << MSG::WARNING
279 <<
"Problem retrieving the TEfficiency with name (including pre- and post-fixes) "
281 <<
" or with name " << effNameCopy
282 <<
" in " <<
m_name <<
"... it doesn't exist, neither in the cached map nor in the THistSvc!"
283 <<
" Will return an NULL pointer... you have to handle it correctly!" <<
endmsg;
288 m_effMap.insert(
m_effMap.end(), std::pair< const hash_t, TEfficiency* >( effHash, effPointer ) );
293 m_effMap.insert(
m_effMap.end(), std::pair< const hash_t, TEfficiency* >( effHash, effPointer ) );
◆ endInputFile()
StatusCode EL::AnaAlgorithm::endInputFile |
( |
| ) |
|
|
protectedinherited |
perform the action for the end of an input file
Ideally you don't use this, but instead rely on meta-data tools instead. However, there are enough people asking for it that I decided to implement it anyways.
\warn To use this you have to call requestEndInputFile to use this.
\warn If a file is split across multiple jobs this will be called more than once. This only happens for specific batch drivers and/or if it is explicitly configured by the user. With PROOF it could even happen multiple times within the same job, and while PROOF is no longer supported that behavior may come back if support for a similar framework is added in the future. As such, this method should not be used for accounting that relies to be called exactly once per file, take a look at fileExecute if you want something that is guaranteed to be executed exactly once per input file.
\warn The execution order of endInputFile and fileExecute is currently unspecified.
Definition at line 358 of file AnaAlgorithm.cxx.
361 return StatusCode::SUCCESS;
◆ evtStore() [1/2]
◆ evtStore() [2/2]
◆ execute()
StatusCode CP::SecVertexTruthMatchAlg::execute |
( |
| ) |
|
|
overridevirtual |
Definition at line 121 of file SecVertexTruthMatchAlg.cxx.
128 std::vector<const xAOD::Vertex*> recoVerticesToMatch;
129 std::vector<const xAOD::TruthVertex*> truthVerticesToMatch;
131 for(
const auto recoVertex : *recoVertexContainer) {
138 recoVerticesToMatch.push_back(recoVertex);
141 for(
const auto truthVertex : *truthVertexContainer) {
142 if(truthVertex->nIncomingParticles() != 1) {
152 if(truthVertex->nOutgoingParticles() < 2) {
155 truthVerticesToMatch.push_back(truthVertex);
159 ATH_CHECK(
m_matchTool->matchVertices( recoVerticesToMatch, truthVerticesToMatch, trackParticleContainer.cptr() ) );
163 for(
const auto& secVtx : recoVerticesToMatch) {
164 int matchTypeBitset = matchTypeDecor(*secVtx);
165 hist(
"RecoVertex/matchType")->Fill(matchTypeBitset);
186 for(
const auto& truthVtx : truthVerticesToMatch) {
187 int truthTypeBitset = truthTypeDecor(*truthVtx);
214 return StatusCode::SUCCESS;
◆ extraDeps_update_handler()
Add StoreName to extra input/output deps as needed.
use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given
◆ extraOutputDeps()
const DataObjIDColl & AthAlgorithm::extraOutputDeps |
( |
| ) |
const |
|
overridevirtualinherited |
Return the list of extra output dependencies.
This list is extended to include symlinks implied by inheritance relations.
Definition at line 50 of file AthAlgorithm.cxx.
57 return Algorithm::extraOutputDeps();
◆ fileExecute()
StatusCode EL::AnaAlgorithm::fileExecute |
( |
| ) |
|
|
protectedinherited |
perform the action exactly once for each file in the dataset
Ideally you don't use this, but instead rely on meta-data tools instead. However, there are enough people asking for it that I decided to implement it anyways.
\warn To use this you have to call requestFileExecute to use this.
\warn The user should not expect this to be called at any particular point in execution. If a file is split between multiple jobs this will be called in only one of these jobs, and not the others. It usually gets called before the first event in a file, but that is not guaranteed and relying on this is a bug.
\warn The execution order of beginInputFile and fileExecute is currently unspecified.
\warn fileExecute does not work with sub-file splitting in Athena, i.e. processing half the events of a file in one job the other half in another job. this should not normally happen, unless you do crazy things like run AthenaMP or explicitly select sub-file splitting in panda. in that case you are on your own.
Definition at line 342 of file AnaAlgorithm.cxx.
345 return StatusCode::SUCCESS;
◆ fillRecoHistograms()
void CP::SecVertexTruthMatchAlg::fillRecoHistograms |
( |
const xAOD::Vertex * |
secVtx, |
|
|
const std::string & |
matchType |
|
) |
| |
|
private |
Definition at line 217 of file SecVertexTruthMatchAlg.cxx.
226 TVector3 reco_pos(secVtx->
x(), secVtx->
y(), secVtx->
z());
227 float Lxy = reco_pos.Perp();
231 ntracks = trkParts.size();
233 TLorentzVector sumP4(0,0,0,0);
237 double minOpAng = -1.0* 1.e10;
238 double maxOpAng = 1.0* 1.e10;
239 double minD0 = 1.0* 1.e10;
246 for(
size_t t = 0;
t < ntracks;
t++){
256 if(trk_d0 < minD0){ minD0 = trk_d0; }
257 if(trk_d0 > maxD0){ maxD0 = trk_d0; }
261 vv.SetPtEtaPhiM(trk.
pt(),trk.
eta(), trk.
phi0(), trk.
m());
263 H +=
vv.Vect().Mag();
266 TLorentzVector v_minus_iv(0,0,0,0);
267 for(
size_t j = 0; j < ntracks; j++){
268 if (j ==
t){
continue; }
278 tmp.SetPtEtaPhiM(trk_2.
pt(),trk_2.
eta(), trk_2.
phi0(), trk_2.
m());
282 double tm =
vv *
tmp / (
vv.Mag() *
tmp.Mag() );
283 if( minOpAng < tm ) minOpAng = tm;
284 if( maxOpAng > tm ) maxOpAng = tm;
287 double DR =
vv.DeltaR(v_minus_iv);
288 if( DR > maxDR ){ maxDR = DR;}
295 if ( Trk_Chi2.isAvailable(trk) && Trk_Chi2(trk) && Trk_nDoF.isAvailable(trk) && Trk_nDoF(trk) ) {
296 hist(
"RecoVertex/" + matchType +
"_Trk_Chi2")->Fill(Trk_Chi2(trk) / Trk_nDoF(trk));
297 hist(
"RecoVertex/" + matchType +
"_Trk_nDoF")->Fill(Trk_nDoF(trk));
299 hist(
"RecoVertex/" + matchType +
"_Trk_D0")->Fill(trk_d0);
300 hist(
"RecoVertex/" + matchType +
"_Trk_Z0")->Fill(trk_z0);
303 hist(
"RecoVertex/" + matchType +
"_Trk_Eta")->Fill(trk.
eta());
304 hist(
"RecoVertex/" + matchType +
"_Trk_Phi")->Fill(trk.
phi0());
305 hist(
"RecoVertex/" + matchType +
"_Trk_E")->Fill(trk.
e() /
GeV);
306 hist(
"RecoVertex/" + matchType +
"_Trk_M")->Fill(trk.
m() /
GeV);
307 hist(
"RecoVertex/" + matchType +
"_Trk_Pt")->Fill(trk.
pt() /
GeV);
308 hist(
"RecoVertex/" + matchType +
"_Trk_Px")->Fill(trk.
p4().Px() /
GeV);
309 hist(
"RecoVertex/" + matchType +
"_Trk_Py")->Fill(trk.
p4().Py() /
GeV);
310 hist(
"RecoVertex/" + matchType +
"_Trk_Pz")->Fill(trk.
p4().Pz() /
GeV);
311 hist(
"RecoVertex/" + matchType +
"_Trk_charge")->Fill(trk.
charge());
316 const double dir = sumP4.Vect().Dot( reco_pos ) / sumP4.Vect().Mag() / reco_pos.Mag();
321 hist(
"RecoVertex/" + matchType +
"_x")->Fill(secVtx->
x());
322 hist(
"RecoVertex/" + matchType +
"_y")->Fill(secVtx->
y());
323 hist(
"RecoVertex/" + matchType +
"_z")->Fill(secVtx->
z());
324 hist(
"RecoVertex/" + matchType +
"_Lxy")->Fill(Lxy);
325 hist(
"RecoVertex/" + matchType +
"_ntrk")->Fill(ntracks);
326 hist(
"RecoVertex/" + matchType +
"_pT")->Fill(sumP4.Pt() /
GeV);
327 hist(
"RecoVertex/" + matchType +
"_eta")->Fill(sumP4.Eta());
328 hist(
"RecoVertex/" + matchType +
"_phi")->Fill(sumP4.Phi());
329 hist(
"RecoVertex/" + matchType +
"_mass")->Fill(sumP4.M() /
GeV);
330 hist(
"RecoVertex/" + matchType +
"_mu")->Fill(sumP4.M()/maxDR /
GeV);
331 hist(
"RecoVertex/" + matchType +
"_chi2")->Fill(Chi2(*secVtx)/nDoF(*secVtx));
332 hist(
"RecoVertex/" + matchType +
"_dir")->Fill(
dir);
333 hist(
"RecoVertex/" + matchType +
"_charge")->Fill(
charge);
334 hist(
"RecoVertex/" + matchType +
"_H")->Fill(H /
GeV);
335 hist(
"RecoVertex/" + matchType +
"_HT")->Fill(HT /
GeV);
336 hist(
"RecoVertex/" + matchType +
"_minOpAng")->Fill(minOpAng);
337 hist(
"RecoVertex/" + matchType +
"_maxOpAng")->Fill(maxOpAng);
338 hist(
"RecoVertex/" + matchType +
"_mind0")->Fill(minD0);
339 hist(
"RecoVertex/" + matchType +
"_maxd0")->Fill(maxD0);
340 hist(
"RecoVertex/" + matchType +
"_maxdR")->Fill(maxDR);
342 std::vector<InDetSecVtxTruthMatchUtils::VertexTruthMatchInfo> truthmatchinfo;
343 truthmatchinfo = matchInfoDecor(*secVtx);
346 if (matchType !=
"All" and matchType !=
"Fake") {
347 if(not truthmatchinfo.empty()){
348 float matchScore_weight = std::get<1>(truthmatchinfo.at(0));
349 float matchScore_pt = std::get<2>(truthmatchinfo.at(0));
351 ATH_MSG_DEBUG(
"Match Score and probability: " << matchScore_weight <<
" " << matchScore_pt/0.01);
356 hist(
"RecoVertex/" + matchType +
"_positionRes_R")->Fill(Lxy - truthVtx.
perp());
357 hist(
"RecoVertex/" + matchType +
"_positionRes_Z")->Fill(secVtx->
z() - truthVtx.
z());
358 hist(
"RecoVertex/" + matchType +
"_matchScore_weight")->Fill(matchScore_weight);
359 hist(
"RecoVertex/" + matchType +
"_matchScore_pt")->Fill(matchScore_pt);
◆ fillTruthHistograms()
◆ finalize()
StatusCode EL::AnaAlgorithm::finalize |
( |
| ) |
|
|
protectedinherited |
finalize this algorithm
This gets called after event processing has finished. The last event may no longer be in memory, and the code should not try to access it.
Definition at line 328 of file AnaAlgorithm.cxx.
331 return StatusCode::SUCCESS;
◆ graph()
TGraph * AthHistogramming::graph |
( |
const std::string & |
graphName, |
|
|
const std::string & |
tDir = "" , |
|
|
const std::string & |
stream = "" |
|
) |
| |
|
protectedinherited |
Simplify the retrieval of registered TGraphs.
Definition at line 492 of file AthHistogramming.cxx.
495 const hash_t graphHash = this->
hash(graphName);
498 GraphMap_t::const_iterator
it =
m_graphMap.find( graphHash );
503 std::string graphNameCopy = graphName;
504 std::string tDirCopy = tDir;
505 std::string streamCopy =
stream;
508 std::string bookingString(
"");
511 TGraph* graphPointer(NULL);
512 if ( !((
histSvc()->getGraph(bookingString, graphPointer)).isSuccess()) )
515 std::string bookingString(
"");
516 this->
buildBookingString( bookingString, graphNameCopy, tDirCopy, streamCopy,
false );
518 if ( !((
histSvc()->getGraph(bookingString, graphPointer)).isSuccess()) )
520 m_msg << MSG::WARNING
521 <<
"Problem retrieving the TGraph with name (including pre- and post-fixes) "
523 <<
" or with name " << graphNameCopy
524 <<
" in " <<
m_name <<
"... it doesn't exist, neither in the cached map nor in the THistSvc!"
525 <<
" Will return an NULL pointer... you have to handle it correctly!" <<
endmsg;
530 m_graphMap.insert(
m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphPointer ) );
535 m_graphMap.insert(
m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphPointer ) );
◆ handle()
void EL::AnaAlgorithm::handle |
( |
const Incident & |
inc | ) |
|
|
inherited |
receive the given incident
- Guarantee
- basic
- Failures
- incident handling errors
Definition at line 520 of file AnaAlgorithm.cxx.
523 if (inc.type() == IncidentType::BeginInputFile)
529 }
else if (inc.type() == IncidentType::EndInputFile)
◆ hash()
◆ hist()
TH1 * AthHistogramming::hist |
( |
const std::string & |
histName, |
|
|
const std::string & |
tDir = "" , |
|
|
const std::string & |
stream = "" |
|
) |
| |
|
protectedinherited |
Simplify the retrieval of registered histograms of any type.
Definition at line 198 of file AthHistogramming.cxx.
204 HistMap_t::const_iterator
it =
m_histMap.find( histHash );
209 std::string histNameCopy =
histName;
210 std::string tDirCopy = tDir;
211 std::string streamCopy =
stream;
214 std::string bookingString(
"");
217 TH1* histPointer(NULL);
218 if ( !((
histSvc()->
getHist(bookingString, histPointer)).isSuccess()) )
221 std::string bookingString(
"");
224 if ( !((
histSvc()->
getHist(bookingString, histPointer)).isSuccess()) )
226 m_msg << MSG::WARNING
227 <<
"Problem retrieving the histogram with name (including pre- and post-fixes) "
229 <<
" or with name " << histNameCopy
230 <<
" in " <<
m_name <<
"... it doesn't exist, neither in the cached map nor in the THistSvc!"
231 <<
" Will return an NULL pointer... you have to handle it correctly!" <<
endmsg;
236 m_histMap.insert(
m_histMap.end(), std::pair< const hash_t, TH1* >( histHash, histPointer ) );
241 m_histMap.insert(
m_histMap.end(), std::pair< const hash_t, TH1* >( histHash, histPointer ) );
◆ hist2d()
TH2 * AthHistogramming::hist2d |
( |
const std::string & |
histName, |
|
|
const std::string & |
tDir = "" , |
|
|
const std::string & |
stream = "" |
|
) |
| |
|
inlineprotectedinherited |
Simplify the retrieval of registered 2-d histograms.
Definition at line 369 of file AthHistogramming.h.
376 <<
"Cannot get a 2-d histogram with name " <<
histName
377 <<
"... will probably seg-fault!" <<
endmsg;
381 return dynamic_cast<TH2*
>( th1Pointer );
◆ hist3d()
TH3 * AthHistogramming::hist3d |
( |
const std::string & |
histName, |
|
|
const std::string & |
tDir = "" , |
|
|
const std::string & |
stream = "" |
|
) |
| |
|
inlineprotectedinherited |
Simplify the retrieval of registered 3-d histograms.
Definition at line 386 of file AthHistogramming.h.
393 <<
"Cannot get a 3-d histogram with name " <<
histName
394 <<
"... will probably seg-fault!" <<
endmsg;
398 return dynamic_cast<TH3*
>( th1Pointer );
◆ histSvc()
The standard THistSvc
(for writing histograms and TTrees and more to a root file) Returns (kind of) a pointer to the THistSvc
.
Definition at line 113 of file AthHistogramAlgorithm.h.
◆ initialize()
StatusCode CP::SecVertexTruthMatchAlg::initialize |
( |
| ) |
|
|
overridevirtual |
Definition at line 20 of file SecVertexTruthMatchAlg.cxx.
31 std::vector<std::string> recoTypes{
"All",
"Matched",
"Merged",
"Fake",
"Split",
"Other"};
32 std::vector<std::string> truthTypes{
"Inclusive",
"Reconstructable",
"Accepted",
"Seeded",
"Reconstructed",
"ReconstructedSplit"};
34 ANA_CHECK (
book(
TH1F(
"RecoVertex/matchType",
"Vertex Match Type", 65, -0.5, 64.5)));
36 for(
const auto& recoType : recoTypes) {
37 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_x").c_str(),
"Reco vertex x [mm]", 1000, -500, 500)));
38 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_y").c_str(),
"Reco vertex y [mm]", 1000, -500, 500)));
39 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_z").c_str(),
"Reco vertex z [mm]", 1000, -500, 500)));
40 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_Lxy").c_str(),
"Reco vertex L_{xy} [mm]", 500, 0, 500)));
41 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_pT").c_str(),
"Reco vertex p_{T} [GeV]", 100, 0, 100)));
42 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_eta").c_str(),
"Reco vertex #eta", 100, -5, 5)));
43 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_phi").c_str(),
"Reco vertex #phi", 100, -TMath::Pi(), TMath::Pi())));
44 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_mass").c_str(),
"Reco vertex mass [GeV]", 500, 0, 100)));
45 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_mu").c_str(),
"Reco vertex Red. Mass [GeV]", 500, 0, 100)));
46 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_chi2").c_str(),
"Reco vertex recoChi2", 100, 0, 10)));
47 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_dir").c_str(),
"Reco vertex recoDirection", 100, -1, 1)));
48 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_charge").c_str(),
"Reco vertex recoCharge", 20, -10, 10)));
49 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_H").c_str(),
"Reco vertex H [GeV]", 100, 0, 100)));
50 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_HT").c_str(),
"Reco vertex Mass [GeV]", 100, 0, 100)));
51 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_minOpAng").c_str(),
"Reco vertex minOpAng", 100, -1, 1)));
52 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_maxOpAng").c_str(),
"Reco vertex maxOpAng", 100, -1, 1)));
53 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_maxdR").c_str(),
"Reco vertex maxDR", 100, 0, 10)));
54 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_mind0").c_str(),
"Reco vertex min d0 [mm]", 100, 0, 100)));
55 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_maxd0").c_str(),
"Reco vertex max d0 [mm]", 100, 0, 100)));
56 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_ntrk").c_str(),
"Reco vertex n tracks", 30, 0, 30)));
59 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_Trk_qOverP").c_str(),
"Reco track qOverP ", 100, 0, .01)));
60 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_Trk_theta").c_str(),
"Reco track theta ", 64, 0, 3.2)));
61 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_Trk_E").c_str(),
"Reco track E ", 100, 0, 100)));
62 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_Trk_M").c_str(),
"Reco track M ", 100, 0, 10)));
63 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_Trk_Pt").c_str(),
"Reco track Pt ", 100, 0, 100)));
64 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_Trk_Px").c_str(),
"Reco track Px ", 100, 0, 100)));
65 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_Trk_Py").c_str(),
"Reco track Py ", 100, 0, 100)));
66 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_Trk_Pz").c_str(),
"Reco track Pz ", 100, 0, 100)));
67 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_Trk_Eta").c_str(),
"Reco track Eta ", 100, -5, 5)));
68 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_Trk_Phi").c_str(),
"Reco track Phi ", 63, -3.2, 3.2)));
69 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_Trk_D0").c_str(),
"Reco track D0 ", 300, -300, 300)));
70 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_Trk_Z0").c_str(),
"Reco track Z0 ", 500, -500, 500)));
71 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_Trk_errD0").c_str(),
"Reco track errD0 ", 300, 0, 30)));
72 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_Trk_errZ0").c_str(),
"Reco track errZ0 ", 500, 0, 50)));
73 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_Trk_Chi2").c_str(),
"Reco track Chi2 ", 100, 0, 10)));
74 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_Trk_nDoF").c_str(),
"Reco track nDoF ", 100, 0, 100)));
75 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_Trk_charge").c_str(),
"Reco track charge ", 3, -1.5, 1.5)));
78 if ( recoType !=
"All" and recoType !=
"Fake" ) {
79 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_positionRes_R").c_str(),
"Position resolution for vertices matched to truth decays", 400, -20, 20)));
80 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_positionRes_Z").c_str(),
"Position resolution for vertices matched to truth decays", 400, -20, 20)));
81 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_matchScore_weight").c_str(),
"Vertex Match Score (weight)", 101, 0, 1.01)));
82 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_matchScore_pt").c_str(),
"Vertex Match Score (pT)", 101, 0, 1.01)));
83 ANA_CHECK (
book(
TH1F((
"RecoVertex/" + recoType +
"_matchedTruthID").c_str(),
"Vertex Truth Match ID", 100, 0, 100)));
101 ANA_CHECK (
book(
TH1F((
"TruthVertex/" +
truthType +
"_ParentProdX").c_str(),
"truthParentProd vertex x [mm]", 500, -500, 500)));
102 ANA_CHECK (
book(
TH1F((
"TruthVertex/" +
truthType +
"_ParentProdY").c_str(),
"truthParentProd vertex y [mm]", 500, -500, 500)));
103 ANA_CHECK (
book(
TH1F((
"TruthVertex/" +
truthType +
"_ParentProdZ").c_str(),
"truthParentProd vertex z [mm]", 500, -500, 500)));
105 ANA_CHECK (
book(
TH1F((
"TruthVertex/" +
truthType +
"_ParentProddistFromPV").c_str(),
"truthParentProd vertex distFromPV [mm]", 500, 0, 500)));
110 Double_t
bins[] = {0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 15, 20, 25, 30, 35, 40, 50, 60, 70, 80, 90, 100, 125, 150, 200, 300, 500};
118 return StatusCode::SUCCESS;
◆ inputHandles()
Return this algorithm's input handles.
We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.
◆ inputMetaStore() [1/2]
AnaAlgorithm::MetaStorePtr_t EL::AnaAlgorithm::inputMetaStore |
( |
| ) |
|
|
inherited |
Definition at line 81 of file AnaAlgorithm.cxx.
83 #ifdef XAOD_STANDALONE
87 #endif // XAOD_STANDALONE
◆ inputMetaStore() [2/2]
AnaAlgorithm::ConstMetaStorePtr_t EL::AnaAlgorithm::inputMetaStore |
( |
| ) |
const |
|
inherited |
Accessor for the input metadata store
Definition at line 72 of file AnaAlgorithm.cxx.
74 #ifdef XAOD_STANDALONE
78 #endif // XAOD_STANDALONE
◆ msg() [1/2]
◆ msg() [2/2]
◆ msgLvl()
◆ myReplace()
void AthHistogramming::myReplace |
( |
std::string & |
str, |
|
|
const std::string & |
oldStr, |
|
|
const std::string & |
newStr |
|
) |
| |
|
privateinherited |
Helper method to replace sub-string.
Definition at line 590 of file AthHistogramming.cxx.
595 while((
pos =
str.find(oldStr,
pos)) != std::string::npos)
597 str.replace(
pos, oldStr.length(), newStr);
598 pos += newStr.length();
◆ outputHandles()
Return this algorithm's output handles.
We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.
◆ outputMetaStore() [1/2]
AnaAlgorithm::MetaStorePtr_t EL::AnaAlgorithm::outputMetaStore |
( |
| ) |
|
|
inherited |
Definition at line 101 of file AnaAlgorithm.cxx.
103 #ifdef XAOD_STANDALONE
107 #endif // XAOD_STANDALONE
◆ outputMetaStore() [2/2]
AnaAlgorithm::ConstMetaStorePtr_t EL::AnaAlgorithm::outputMetaStore |
( |
| ) |
const |
|
inherited |
Accessor for the output metadata store
Definition at line 92 of file AnaAlgorithm.cxx.
94 #ifdef XAOD_STANDALONE
98 #endif // XAOD_STANDALONE
◆ print()
void EL::AnaAlgorithm::print |
( |
| ) |
const |
|
protectedvirtualinherited |
print the state of the algorithm
This is mostly to allow algorithms to add a little debugging information if they feel like it.
Definition at line 336 of file AnaAlgorithm.cxx.
◆ renounce()
◆ renounceArray()
◆ requestBeginInputFile()
StatusCode EL::AnaAlgorithm::requestBeginInputFile |
( |
| ) |
|
|
inherited |
register this algorithm to have an implementation of beginInputFile
- Guarantee
- strong
- Failures
- beginInputFile not supported
Definition at line 270 of file AnaAlgorithm.cxx.
275 #ifndef XAOD_STANDALONE
282 incSvc->removeListener(
this, IncidentType::BeginInputFile );
283 incSvc->addListener(
this, IncidentType::BeginInputFile, 0,
true );
286 return StatusCode::SUCCESS;
◆ requestEndInputFile()
StatusCode EL::AnaAlgorithm::requestEndInputFile |
( |
| ) |
|
|
inherited |
register this algorithm to have an implementation of endInputFile
- Guarantee
- strong
- Failures
- endInputFile not supported
Definition at line 291 of file AnaAlgorithm.cxx.
296 #ifndef XAOD_STANDALONE
303 incSvc->removeListener(
this, IncidentType::EndInputFile );
304 incSvc->addListener(
this, IncidentType::EndInputFile, 0,
true );
307 return StatusCode::SUCCESS;
◆ requestFileExecute()
StatusCode EL::AnaAlgorithm::requestFileExecute |
( |
| ) |
|
|
inherited |
register this algorithm to have an implementation of fileexecute
- Guarantee
- strong
- Failures
- fileExecute not supported
Definition at line 249 of file AnaAlgorithm.cxx.
254 #ifndef XAOD_STANDALONE
261 incSvc->removeListener(
this, IncidentType::BeginInputFile );
262 incSvc->addListener(
this, IncidentType::BeginInputFile, 0,
true );
265 return StatusCode::SUCCESS;
◆ sysInitialize()
StatusCode AthHistogramAlgorithm::sysInitialize |
( |
| ) |
|
|
virtualinherited |
Initialization method invoked by the framework.
This method is responsible for any bookkeeping of initialization required by the framework itself. It will in turn invoke the initialize() method of the derived algorithm, and of any sub-algorithms which it creates.
Reimplemented from AthAlgorithm.
Reimplemented in AthAnalysisAlgorithm.
Definition at line 75 of file AthHistogramAlgorithm.cxx.
80 if ( Gaudi::StateMachine::INITIALIZED <= FSMState() )
return StatusCode::SUCCESS;
83 bindPropertiesTo( serviceLocator()->getOptsSvc() );
87 if ( !isEnabled( ) )
return StatusCode::SUCCESS;
◆ sysStart()
Handle START transition.
We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.
◆ tree()
TTree * AthHistogramming::tree |
( |
const std::string & |
treeName, |
|
|
const std::string & |
tDir = "" , |
|
|
const std::string & |
stream = "" |
|
) |
| |
|
protectedinherited |
Simplify the retrieval of registered TTrees.
Definition at line 378 of file AthHistogramming.cxx.
384 TreeMap_t::const_iterator
it =
m_treeMap.find( treeHash );
389 std::string treeNameCopy =
treeName;
390 std::string tDirCopy = tDir;
391 std::string streamCopy =
stream;
394 std::string bookingString(
"");
397 TTree* treePointer(NULL);
398 if ( !((
histSvc()->
getTree(bookingString, treePointer)).isSuccess()) )
400 m_msg << MSG::WARNING
401 <<
"Problem retrieving the TTree with name " << treeNameCopy
402 <<
" in " <<
m_name <<
"... it doesn't exist, neither in the cached map nor in the THistSvc!"
403 <<
" Will return an NULL pointer... you have to handle it correctly!" <<
endmsg;
408 m_treeMap.insert(
m_treeMap.end(), std::pair< const hash_t, TTree* >( treeHash, treePointer ) );
◆ updateVHKA()
◆ m_detStore
◆ m_effMap
◆ m_evtStore
◆ m_extendedExtraObjects
DataObjIDColl AthAlgorithm::m_extendedExtraObjects |
|
privateinherited |
◆ m_graphMap
◆ m_hasBeginInputFile
bool EL::AnaAlgorithm::m_hasBeginInputFile {false} |
|
privateinherited |
◆ m_hasEndInputFile
bool EL::AnaAlgorithm::m_hasEndInputFile {false} |
|
privateinherited |
◆ m_hasFileExecute
bool EL::AnaAlgorithm::m_hasFileExecute {false} |
|
privateinherited |
◆ m_histMap
◆ m_histNamePostfix
std::string AthHistogramAlgorithm::m_histNamePostfix |
|
privateinherited |
◆ m_histNamePrefix
std::string AthHistogramAlgorithm::m_histNamePrefix |
|
privateinherited |
◆ m_histSvc
◆ m_histTitlePostfix
std::string AthHistogramAlgorithm::m_histTitlePostfix |
|
privateinherited |
◆ m_histTitlePrefix
std::string AthHistogramAlgorithm::m_histTitlePrefix |
|
privateinherited |
◆ m_inputMetaStore
Object accessing the input metadata store.
Definition at line 551 of file AnaAlgorithm.h.
◆ m_matchTool
◆ m_msg
MsgStream AthHistogramming::m_msg |
|
privateinherited |
◆ m_name
std::string AthHistogramming::m_name |
|
privateinherited |
◆ m_outputMetaStore
Object accessing the output metadata store.
Definition at line 555 of file AnaAlgorithm.h.
◆ m_prefix
std::string AthHistogramAlgorithm::m_prefix |
|
privateinherited |
◆ m_rootDir
std::string AthHistogramAlgorithm::m_rootDir |
|
privateinherited |
◆ m_secVtxContainerKey
Initial value:{this, "SecondaryVertexContainer", "VrtSecInclusive_SecondaryVertices",
"Secondary vertex container"}
Definition at line 41 of file SecVertexTruthMatchAlg.h.
◆ m_streamName
std::string AthHistogramming::m_streamName |
|
privateinherited |
◆ m_targetPDGIDs
Gaudi::Property<std::vector<int> > CP::SecVertexTruthMatchAlg::m_targetPDGIDs {this, "TargetPDGIDs", {}, "List of PDGIDs to select for matching"} |
|
private |
◆ m_trackParticleContainerKey
Initial value:{this, "TrackParticleContainer", "InDetTrackParticles",
"Track container"}
Definition at line 47 of file SecVertexTruthMatchAlg.h.
◆ m_treeMap
◆ m_truthVtxContainerKey
Initial value:{this, "TruthVertexContainer", "TruthVertices",
"Truth vertex container"}
Definition at line 44 of file SecVertexTruthMatchAlg.h.
◆ m_varHandleArraysDeclared
◆ m_vhka
◆ m_writeHistograms
Gaudi::Property<bool> CP::SecVertexTruthMatchAlg::m_writeHistograms {this, "WriteHistograms", true, "Write histograms"} |
|
private |
The documentation for this class was generated from the following files:
std::string m_rootDir
Name of the ROOT directory.
def retrieve(aClass, aKey=None)
const ServiceHandle< ITHistSvc > & histSvc() const
The standard THistSvc (for writing histograms and TTrees and more to a root file) Returns (kind of) a...
const ServiceHandle< ITHistSvc > & histSvc() const
The standard THistSvc (for writing histograms and TTrees and more to a root file) Returns (kind of) a...
size_t nOutgoingParticles() const
Get the number of outgoing particles.
virtual double pt() const override final
The transverse momentum ( ) of the particle.
bool m_hasBeginInputFile
the value of hasBeginInputFile
float x() const
Returns the x position.
virtual ::StatusCode endInputFile()
perform the action for the end of an input file
virtual double m() const override final
The invariant mass of the particle..
void myReplace(std::string &str, const std::string &oldStr, const std::string &newStr)
Helper method to replace sub-string.
MsgStream m_msg
Cached Message Stream.
HistMap_t m_histMap
The map of histogram names to their pointers.
float phi() const
Vertex azimuthal angle.
ServiceHandle< ITHistSvc > m_histSvc
Default constructor: AthHistogramAlgorithm();.
std::string find(const std::string &s)
return a remapped string
StatusCode book(const TH1 &hist, const std::string &tDir="", const std::string &stream="")
Simplify the booking and registering (into THistSvc) of histograms.
bool isFake(int matchInfo)
float charge() const
Returns the charge.
bool isReconstructable(int matchInfo)
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
virtual double eta() const override final
The pseudorapidity ( ) of the particle.
TEfficiency * efficiency(const std::string &effName, const std::string &tDir="", const std::string &stream="")
Simplify the retrieval of registered TEfficiency.
MetaStore_t m_inputMetaStore
Object accessing the input metadata store.
MetaStore_t m_outputMetaStore
Object accessing the output metadata store.
std::string m_prefix
Name of the ROOT output stream (file)
virtual ::StatusCode beginInputFile()
perform the action for the beginning of an input file
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
std::vector< SG::VarHandleKeyArray * > m_vhka
bool isMerged(int matchInfo)
std::string m_histNamePrefix
The prefix for the histogram THx name.
bool isMatched(int matchInfo)
std::string m_histTitlePrefix
The prefix for the histogram THx title.
Helper class to provide constant type-safe access to aux data.
bool isSplit(int matchInfo)
bool isValid(const T &p)
Av: we implement here an ATLAS-sepcific convention: all particles which are 99xxxxx are fine.
std::string m_histNamePostfix
The postfix for the histogram THx name.
virtual ::StatusCode fileExecute()
perform the action exactly once for each file in the dataset
std::string getTree(const std::string &file_name)
float y() const
Vertex y displacement.
int truthType(const U &p)
virtual void setOwner(IDataHandleHolder *o)=0
ToolHandle< IInDetSecVtxTruthMatchTool > m_matchTool
virtual FourMom_t p4() const override final
The full 4-momentum of the particle.
void fillTruthHistograms(const xAOD::TruthVertex *truthVtx, const std::string &truthType)
Gaudi::Property< bool > m_writeHistograms
the (new) base class for EventLoop algorithms
virtual StatusCode sysInitialize() override
Override sysInitialize.
bool isOther(int matchInfo)
bool isReconstructed(int matchInfo)
virtual double e() const override final
The total energy of the particle.
Helper class to provide type-safe access to aux data.
bool isAccepted(int matchInfo)
@ SecVtx
Secondary vertex.
Class describing a truth particle in the MC record.
SG::ReadHandleKey< xAOD::VertexContainer > m_secVtxContainerKey
Gaudi::Property< std::vector< int > > m_targetPDGIDs
float perp() const
Vertex transverse distance from the beam line.
bool hasProdVtx() const
Check for a production vertex on this particle.
bool m_hasFileExecute
the value of hasFileExecute
std::vector< ElementLink< xAOD::TrackParticleContainer > > TrackParticleLinks_t
Type for the associated track particles.
GraphMap_t m_graphMap
The map of TGraph names to their pointers.
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
float z() const
Returns the z position.
const TruthParticle_v1 * incomingParticle(size_t index) const
Get one of the incoming particles.
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
std::string m_streamName
Name of the ROOT output stream (file)
StatusCode configAthHistogramming(const ServiceHandle< ITHistSvc > &histSvc, const std::string &prefix, const std::string &rootDir, const std::string &histNamePrefix, const std::string &histNamePostfix, const std::string &histTitlePrefix, const std::string &histTitlePostfix)
To be called by the derived classes to fill the internal configuration.
uint32_t hash_t
typedef for the internal hash
uint64_t crc64(const CRCTable &table, const char *data, size_t data_len)
Find the CRC-64 of a string,.
std::string m_histTitlePostfix
The postfix for the histogram THx title.
virtual void renounce()=0
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
ElementLink implementation for ROOT usage.
const TruthVertex_v1 * prodVtx() const
The production vertex of this particle.
Class describing a truth vertex in the MC record.
float phi0() const
Returns the parameter, which has range to .
const ParametersCovMatrix_t definingParametersCovMatrix() const
Returns the 5x5 symmetric matrix containing the defining parameters covariance matrix.
void fillRecoHistograms(const xAOD::Vertex *secVtx, const std::string &matchType)
double charge(const T &p)
hash_t hash(const std::string &histName) const
Method to calculate a 32-bit hash from a string.
def getHist(name, makehist)
virtual double e() const override final
The total energy of the particle.
virtual double eta() const override final
The pseudorapidity ( ) of the particle.
DefiningParameters_t definingParameters() const
Returns a SVector of the Perigee track parameters.
float x() const
Vertex x displacement.
std::uint64_t hash64(const void *data, std::size_t size)
Passthrough to XXH3_64bits.
bool isReconstructedSplit(int matchInfo)
DataObjIDColl m_extendedExtraObjects
EffMap_t m_effMap
The map of histogram names to their pointers.
std::string m_histTitlePrefix
The prefix for the histogram THx title.
float eta() const
Vertex pseudorapidity.
virtual double phi() const override final
The azimuthal angle ( ) of the particle.
#define ATH_MSG_WARNING(x)
std::string m_histNamePostfix
The postfix for the histogram THx name.
float z() const
Vertex longitudinal distance along the beam line form the origin.
TreeMap_t m_treeMap
The map of TTree names to their pointers.
float y() const
Returns the y position.
std::string m_histNamePrefix
The prefix for the histogram THx name.
bool m_hasEndInputFile
the value of hasEndInputFile
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
virtual double pt() const override final
The transverse momentum ( ) of the particle.
std::string m_histTitlePostfix
The postfix for the histogram THx title.
TH1 * hist(const std::string &histName, const std::string &tDir="", const std::string &stream="")
Simplify the retrieval of registered histograms of any type.
ServiceHandle< ITHistSvc > m_histSvc
Pointer to the THistSvc (event store by default)
def TH1F(name, title, nxbins, bins_par2, bins_par3=None, path='', **kwargs)
SG::ReadHandleKey< xAOD::TruthVertexContainer > m_truthVtxContainerKey
Class describing a TrackParticle.
bool isSeeded(int matchInfo)
std::string m_name
Instance name.
std::string m_rootDir
Name of the ROOT directory.
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
SG::ReadHandleKey< xAOD::TrackParticleContainer > m_trackParticleContainerKey
virtual double m() const override final
The mass of the particle.
double charge() const
Physical charge.
void buildBookingString(std::string &bookingString, std::string &histName, std::string &tDir, std::string &stream, bool usePrefixPostfix=false)
Method to build individual booking string.
TH1 * bookGetPointer(const TH1 &hist, const std::string &tDir="", const std::string &stream="")
Simplify the booking and registering (into THistSvc) of histograms.