ATLAS Offline Software
Loading...
Searching...
No Matches
InDetGNNHardScatterSelection::VertexDecoratorAlg Class Referencefinal

#include <VertexDecoratorAlg.h>

Inheritance diagram for InDetGNNHardScatterSelection::VertexDecoratorAlg:

Public Member Functions

 VertexDecoratorAlg (const std::string &name, ISvcLocator *pSvcLocator)
StatusCode initialize () override
StatusCode execute () override
virtual StatusCode sysInitialize ()
 Initialization method invoked by the framework.
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.
virtual const DataObjIDColl & extraOutputDeps () const override
 Return the list of extra output dependencies.
ServiceHandle< StoreGateSvc > & evtStore ()
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
const ServiceHandle< StoreGateSvc > & detStore () const
 The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
virtual StatusCode sysStart () override
 Handle START transition.
virtual std::vector< Gaudi::DataHandle * > inputHandles () const override
 Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * > outputHandles () const override
 Return this algorithm's output handles.
Gaudi::Details::PropertyBase & declareProperty (Gaudi::Property< T, V, H > &t)
void updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream & msg () const
bool msgLvl (const MSG::Level lvl) const

Protected Member Functions

void renounceArray (SG::VarHandleKeyArray &handlesArray)
 remove all handles from I/O resolution
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.
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.
TH1 * bookGetPointer (const TH1 &hist, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of histograms.
TH1 * bookGetPointer (TH1 *hist, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of histograms.
TH1 * bookGetPointer (TH1 &histRef, std::string tDir="", std::string stream="")
 Simplify the booking and registering (into THistSvc) of histograms.
TTree * bookGetPointer (const TTree &treeRef, std::string tDir="", std::string stream="")
 Simplify the booking and registering (into THistSvc) of TTrees.
TGraph * bookGetPointer (const TGraph &graphRef, std::string tDir="", std::string stream="")
 Simplify the booking and registering (into THistSvc) of TGraphs.
TEfficiency * bookGetPointer (const TEfficiency &eff, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency.
TEfficiency * bookGetPointer (TEfficiency *eff, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency.
TEfficiency * bookGetPointer (TEfficiency &effRef, std::string tDir="", std::string stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency.
StatusCode book (const TH1 &hist, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of histograms.
StatusCode book (TH1 *hist, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of histograms.
StatusCode book (TH1 &histRef, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of histograms.
StatusCode book (const TTree &treeRef, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of TTrees.
StatusCode book (const TGraph &graphRef, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of TGraphs.
StatusCode book (const TEfficiency &eff, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency.
StatusCode book (TEfficiency *eff, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency.
StatusCode book (TEfficiency &effRef, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency.
TH1 * hist (const std::string &histName, const std::string &tDir="", const std::string &stream="")
 Simplify the retrieval of registered histograms of any type.
TH2 * hist2d (const std::string &histName, const std::string &tDir="", const std::string &stream="")
 Simplify the retrieval of registered 2-d histograms.
TH3 * hist3d (const std::string &histName, const std::string &tDir="", const std::string &stream="")
 Simplify the retrieval of registered 3-d histograms.
TTree * tree (const std::string &treeName, const std::string &tDir="", const std::string &stream="")
 Simplify the retrieval of registered TTrees.
TGraph * graph (const std::string &graphName, const std::string &tDir="", const std::string &stream="")
 Simplify the retrieval of registered TGraphs.
TEfficiency * efficiency (const std::string &effName, const std::string &tDir="", const std::string &stream="")
 Simplify the retrieval of registered TEfficiency.

Private Types

typedef ServiceHandle< StoreGateSvcStoreGateSvc_t
typedef uint32_t hash_t
 typedef for the internal hash
typedef std::map< const hash_t, TH1 * > HistMap_t
 Typedef for convenience.
typedef std::map< const hash_t, TEfficiency * > EffMap_t
 Typedef for convenience.
typedef std::map< const hash_t, TTree * > TreeMap_t
 Typedef for convenience.
typedef std::map< const hash_t, TGraph * > GraphMap_t
 Typedef for convenience.

Private Member Functions

Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey>
void buildBookingString (std::string &bookingString, std::string &histName, std::string &tDir, std::string &stream, bool usePrefixPostfix=false)
 Method to build individual booking string.
void myReplace (std::string &str, const std::string &oldStr, const std::string &newStr)
 Helper method to replace sub-string.
hash_t hash (const std::string &histName) const
 Method to calculate a 32-bit hash from a string.

Private Attributes

SG::ReadHandleKey< xAOD::VertexContainerm_vertexInKey
SG::ReadHandleKey< xAOD::EventInfom_eventInKey
SG::ReadHandleKey< xAOD::ElectronContainerm_electronsInKey
SG::ReadHandleKey< xAOD::MuonContainerm_muonsInKey
SG::ReadHandleKey< xAOD::JetContainerm_jetsInKey
SG::ReadHandleKey< xAOD::PhotonContainerm_photonsInKey
SG::WriteDecorHandleKey< xAOD::VertexContainerm_mDecor_ntrk
SG::WriteDecorHandleKey< xAOD::VertexContainerm_mDecor_sumPt
SG::WriteDecorHandleKey< xAOD::VertexContainerm_mDecor_chi2Over_ndf
SG::WriteDecorHandleKey< xAOD::VertexContainerm_mDecor_z_asym
SG::WriteDecorHandleKey< xAOD::VertexContainerm_mDecor_weighted_z_asym
SG::WriteDecorHandleKey< xAOD::VertexContainerm_mDecor_weighted_z_kurt
SG::WriteDecorHandleKey< xAOD::VertexContainerm_mDecor_z_skew
SG::WriteDecorHandleKey< xAOD::VertexContainerm_mDecor_photon_deltaz
SG::WriteDecorHandleKey< xAOD::VertexContainerm_mDecor_photon_deltaPhi
SG::WriteDecorHandleKey< xAOD::VertexContainerm_mDecor_actualInterPerXing
SG::WriteDecorHandleKey< xAOD::VertexContainerm_photonLinksKey
SG::WriteDecorHandleKey< xAOD::VertexContainerm_jetLinksKey
SG::WriteDecorHandleKey< xAOD::VertexContainerm_electronLinksKey
SG::WriteDecorHandleKey< xAOD::VertexContainerm_muonLinksKey
SG::WriteHandleKey< xAOD::CompositeParticleContainerm_multiPhotonsOutKey
SG::WriteDecorHandleKey< xAOD::VertexContainerm_multiPhotonLinksKey
SG::ReadDecorHandleKey< xAOD::VertexContainerm_deltaZKey
SG::ReadDecorHandleKey< xAOD::VertexContainerm_deltaPhiKey
SG::ReadDecorHandleKey< xAOD::PhotonContainerm_caloPointingZKey
SG::ReadDecorHandleKey< xAOD::PhotonContainerm_zCommonKey
SG::ReadDecorHandleKey< xAOD::PhotonContainerm_zCommonErrorKey
ToolHandle< CP::TrackVertexAssociationToolm_trkVtxAssociationTool
ToolHandle< GNNToolm_gnnTool
ServiceHandle< ITHistSvc > m_histSvc
 Default constructor: AthHistogramAlgorithm();.
std::string m_prefix
 Name of the ROOT output stream (file).
std::string m_rootDir
 Name of the ROOT directory.
std::string m_histNamePrefix
 The prefix for the histogram THx name.
std::string m_histNamePostfix
 The postfix for the histogram THx name.
std::string m_histTitlePrefix
 The prefix for the histogram THx title.
std::string m_histTitlePostfix
 The postfix for the histogram THx title.
DataObjIDColl m_extendedExtraObjects
StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default).
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default).
std::vector< SG::VarHandleKeyArray * > m_vhka
bool m_varHandleArraysDeclared
HistMap_t m_histMap
 The map of histogram names to their pointers.
EffMap_t m_effMap
 The map of histogram names to their pointers.
TreeMap_t m_treeMap
 The map of TTree names to their pointers.
GraphMap_t m_graphMap
 The map of TGraph names to their pointers.
std::string m_streamName
 Name of the ROOT output stream (file).
std::string m_name
 Instance name.
MsgStream m_msg
 Cached Message Stream.

Detailed Description

Definition at line 40 of file VertexDecoratorAlg.h.

Member Typedef Documentation

◆ EffMap_t

typedef std::map< const hash_t, TEfficiency* > AthHistogramming::EffMap_t
privateinherited

Typedef for convenience.

Definition at line 207 of file AthHistogramming.h.

◆ GraphMap_t

typedef std::map< const hash_t, TGraph* > AthHistogramming::GraphMap_t
privateinherited

Typedef for convenience.

Definition at line 221 of file AthHistogramming.h.

◆ hash_t

typedef uint32_t AthHistogramming::hash_t
privateinherited

typedef for the internal hash

Definition at line 171 of file AthHistogramming.h.

◆ HistMap_t

typedef std::map< const hash_t, TH1* > AthHistogramming::HistMap_t
privateinherited

Typedef for convenience.

Definition at line 200 of file AthHistogramming.h.

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< Algorithm > >::StoreGateSvc_t
privateinherited

Definition at line 388 of file AthCommonDataStore.h.

◆ TreeMap_t

typedef std::map< const hash_t, TTree* > AthHistogramming::TreeMap_t
privateinherited

Typedef for convenience.

Definition at line 214 of file AthHistogramming.h.

Constructor & Destructor Documentation

◆ VertexDecoratorAlg()

VertexDecoratorAlg::VertexDecoratorAlg ( const std::string & name,
ISvcLocator * pSvcLocator )

Definition at line 19 of file VertexDecoratorAlg.cxx.

21 : AthHistogramAlgorithm (name, pSvcLocator)
22 {
23 }
AthHistogramAlgorithm(const std::string &name, ISvcLocator *pSvcLocator)
Constructor with parameters:

Member Function Documentation

◆ 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 337 of file AthHistogramming.h.

338{
339 // We need to create a non-const clone
340 TEfficiency* effClone = dynamic_cast< TEfficiency* >( eff.Clone() );
341 if ( !effClone ) {
342 m_msg << MSG::ERROR << "Couldn't create a TEfficiency clone" << endmsg;
343 return StatusCode::FAILURE;
344 }
345 return this->book( *effClone, tDir, stream );
346}
#define endmsg
StatusCode book(const TH1 &hist, const std::string &tDir="", const std::string &stream="")
Simplify the booking and registering (into THistSvc) of histograms.
MsgStream m_msg
Cached Message Stream.

◆ 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 305 of file AthHistogramming.h.

306{
307 // We need to create a non-const clone
308 TH1* histClone = dynamic_cast< TH1* >( hist.Clone() );
309 if ( !histClone ) {
310 m_msg << MSG::ERROR << "Couldn't create a TH1 clone" << endmsg;
311 return StatusCode::FAILURE;
312 }
313 return this->book( *histClone, tDir, stream );
314}
TH1 * hist(const std::string &histName, const std::string &tDir="", const std::string &stream="")
Simplify the retrieval of registered histograms of any type.

◆ 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 405 of file AthHistogramming.h.

406{
407 // Call the other Book method and see if it returns a valid pointer
408 TTree* treePointer = this->bookGetPointer( treeRef, tDir, stream );
409 if ( treePointer )
410 {
411 return StatusCode::SUCCESS;
412 }
413 else
414 {
415 return StatusCode::FAILURE;
416 }
417}
TH1 * bookGetPointer(const TH1 &hist, const std::string &tDir="", const std::string &stream="")
Simplify the booking and registering (into THistSvc) of histograms.

◆ 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 358 of file AthHistogramming.h.

359{
360 // Call the other Book method and see if it returns a valid pointer
361 TEfficiency* effPointer = this->bookGetPointer( effRef, tDir, stream );
362 if ( !effPointer ) {
363 m_msg << MSG::ERROR << "Couldn't book a TEfficiency" << endmsg;
364 return StatusCode::FAILURE;
365 }
366 return StatusCode::SUCCESS;
367}

◆ 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 348 of file AthHistogramming.h.

349{
350 if ( !eff ) {
351 m_msg << MSG::ERROR << "Got a zero pointer to a TEfficiency" << endmsg;
352 return StatusCode::FAILURE;
353 }
354 return this->book( *eff, tDir, stream );
355}

◆ 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 326 of file AthHistogramming.h.

327{
328 // Call the other Book method and see if it returns a valid pointer
329 TH1* histPointer = this->bookGetPointer( histRef, tDir, stream );
330 if ( !histPointer ) {
331 m_msg << MSG::ERROR << "Couldn't book a TH1" << endmsg;
332 return StatusCode::FAILURE;
333 }
334 return StatusCode::SUCCESS;
335}

◆ 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 316 of file AthHistogramming.h.

317{
318 if ( !hist ) {
319 m_msg << MSG::ERROR << "Got a zero pointer to a TH1" << endmsg;
320 return StatusCode::FAILURE;
321 }
322 return this->book( *hist, tDir, stream );
323}

◆ 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 283 of file AthHistogramming.h.

284{
285 // We need to create a non-const clone
286 TEfficiency* histClone = dynamic_cast< TEfficiency* >( hist.Clone() );
287 if ( !histClone ) {
288 m_msg << MSG::ERROR << "Couldn't create a TEfficiency clone in bookGetPointer" << endmsg;
289 return 0;
290 }
291 return this->bookGetPointer( *histClone, tDir, stream );
292
293}

◆ 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 412 of file AthHistogramming.cxx.

413{
414 // Get a pointer
415 const TGraph* graphPointer = &graphRef;
416
417 // Check that we got a valid pointer
418 if ( !graphPointer )
419 {
420 m_msg << MSG::WARNING
421 << "We got an invalid TGraph pointer in the BookGetPointer(TGraph*) method of the class" << m_name
422 << "!" << endmsg;
423 return NULL;
424 }
425
426 // Modify the name and title according to the prefixes of this classes instance
427 std::string graphName = graphPointer->GetName();
428 const std::string graphTitle = graphPointer->GetTitle();
429
430 // Check if the hash for this graphName already exists, i.e., if we have a hash collision
431 const hash_t graphHash = this->hash(graphName);
432 GraphMap_t::const_iterator it = m_graphMap.find( graphHash );
433 if ( it != m_graphMap.end() ) // It does exist!
434 {
435 m_msg << MSG::WARNING
436 << "Detected a hash collision. The hash for the TGraph with name=" << graphName
437 << " already exists and points to a TGraph with name=" << it->second->GetName()
438 << " NOT going to book the new histogram and returning a NULL pointer!" << endmsg;
439 return NULL;
440 }
441
442 // Create a clone that has the new name
443 TGraph* graphClone = dynamic_cast< TGraph* >( graphPointer->Clone((m_histNamePrefix+graphName+m_histNamePostfix).c_str()) );
444 if( !graphClone )
445 {
446 m_msg << MSG::WARNING
447 << "We couldn't clone the TGraph in the BookGetPointer(TGraph&) method of the class" << m_name
448 << "!" << endmsg;
449 return NULL;
450 }
451 graphClone->SetTitle ((m_histTitlePrefix+graphTitle+m_histTitlePostfix).c_str());
452
453 // Massage the final string to book things
454 std::string bookingString("");
455 this->buildBookingString( bookingString, graphName, tDir, stream );
456
457 // Register the TGraph into the THistSvc
458 if ( !((histSvc()->regGraph(bookingString, graphClone)).isSuccess()) )
459 {
460 m_msg << MSG::WARNING
461 << "Problem registering TGraph with name " << graphName
462 << ", title " << graphTitle
463 << " in " << m_name << "!" << endmsg;
464 return NULL;
465 }
466
467 // Also register it in the local map of string to pointer
468 m_graphMap.insert( m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphClone ) );
469
470 return graphClone;
471}
const ServiceHandle< ITHistSvc > & histSvc() const
The standard THistSvc (for writing histograms and TTrees and more to a root file) Returns (kind of) a...
std::string m_histNamePostfix
The postfix for the histogram THx name.
hash_t hash(const std::string &histName) const
Method to calculate a 32-bit hash from a string.
uint32_t hash_t
typedef for the internal hash
std::string m_histTitlePostfix
The postfix for the histogram THx title.
std::string m_histTitlePrefix
The prefix for the histogram THx title.
std::string m_name
Instance name.
void buildBookingString(std::string &bookingString, std::string &histName, std::string &tDir, std::string &stream, bool usePrefixPostfix=false)
Method to build individual booking string.
std::string m_histNamePrefix
The prefix for the histogram THx name.
GraphMap_t m_graphMap
The map of TGraph names to their pointers.

◆ 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 262 of file AthHistogramming.h.

263{
264 // We need to create a non-const clone
265 TH1* histClone = dynamic_cast< TH1* >( hist.Clone() );
266 if ( !histClone ) {
267 m_msg << MSG::ERROR << "Couldn't create a TH1 clone in bookGetPointer" << endmsg;
268 return 0;
269 }
270 return this->bookGetPointer( *histClone, tDir, stream );
271
272}

◆ 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 297 of file AthHistogramming.cxx.

298{
299 // Get a pointer
300 const TTree* treePointer = &treeRef;
301
302 // Check that we got a valid pointer
303 if ( !treePointer )
304 {
305 m_msg << MSG::WARNING
306 << "We got an invalid TTree pointer in the BookGetPointer(TTree*) method of the class" << m_name
307 << "!" << endmsg;
308 return NULL;
309 }
310
311 // Modify the name and title according to the prefixes of this classes instance
312 std::string treeName = treePointer->GetName();
313 const std::string treeTitle = treePointer->GetTitle();
314
315 // Check if the hash for this treeName already exists, i.e., if we have a hash collision
316 const hash_t treeHash = this->hash(treeName);
317 TreeMap_t::const_iterator it = m_treeMap.find( treeHash );
318 if ( it != m_treeMap.end() ) // It does exist!
319 {
320 m_msg << MSG::WARNING
321 << "Detected a hash collision. The hash for the TTree with name=" << treeName
322 << " already exists and points to a TTree with name=" << it->second->GetName()
323 << " NOT going to book the new histogram and returning a NULL pointer!" << endmsg;
324 return NULL;
325 }
326
327 // Create a clone that has the new name
328 TTree* treeClone = dynamic_cast< TTree* >( treePointer->Clone(treeName.c_str()) );
329 if( !treeClone )
330 {
331 m_msg << MSG::WARNING
332 << "We couldn't clone the TTree in the BookGetPointer(TTree&) method of the class" << m_name
333 << "!" << endmsg;
334 return NULL;
335 }
336 treeClone->SetTitle (treeTitle.c_str());
337
338 // Massage the final string to book things
339 std::string bookingString("");
340 this->buildBookingString( bookingString, treeName, tDir, stream );
341
342 // Register the TTree into the THistSvc
343 if ( !((histSvc()->regTree(bookingString, treeClone)).isSuccess()) )
344 {
345 m_msg << MSG::WARNING
346 << "Problem registering TTree with name " << treeName
347 << ", title " << treeTitle
348 << " in " << m_name << "!" << endmsg;
349 return NULL;
350 }
351
352 // Also register it in the local map of string to pointer
353 m_treeMap.insert( m_treeMap.end(), std::pair< const hash_t, TTree* >( treeHash, treeClone ) );
354
355 return treeClone;
356}
TreeMap_t m_treeMap
The map of TTree names to their pointers.

◆ 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 131 of file AthHistogramming.cxx.

132{
133 // Modify the name and title according to the prefixes of this classes instance
134 std::string effName(effRef.GetName());
135 const std::string effTitle(effRef.GetTitle());
136 std::string bookingString("");
137
138 this->buildBookingString( bookingString, effName, tDir, stream );
139 effRef.SetTitle((m_histTitlePrefix+effTitle+m_histTitlePostfix).c_str() );
140 effRef.SetName(effName.c_str());
141
142 // Check if the hash for this effName already exists, i.e., if we have a hash collision
143 const hash_t effHash = this->hash(effName);
144 EffMap_t::const_iterator it = m_effMap.find( effHash );
145 if ( it != m_effMap.end() ) // It does exist!
146 {
147 m_msg << MSG::WARNING
148 << "Detected a hash collision. The hash for the TEfficiency with name=" << effName
149 << " already exists and points to a TEfficiency with name=" << it->second->GetName()
150 << " NOT going to book the new TEfficiency and returning a NULL pointer!" << endmsg;
151 return NULL;
152 }
153
154 // Set the new name and title for the TEfficiency, based on the prefixes that the user set for this class instance
155 // Create a clone that has the new name
156
157 // Massage the final string to book things
158
159 // Register the TEfficiency into the THistSvc
160 if ( !((histSvc()->regEfficiency(bookingString, &effRef)).isSuccess()) )
161 {
162 m_msg << MSG::WARNING
163 << "Problem registering TEfficiency with name " << effName
164 << ", name prefix " << m_histNamePrefix
165 << ", title " << effTitle
166 << ", tile prefix " << m_histTitlePrefix
167 << ", and tile postfix " << m_histTitlePostfix
168 << " in " << m_name << "!" << endmsg;
169 return NULL;
170 }
171
172 // Also register it in the local map of string to pointer
173 m_effMap.insert( m_effMap.end(), std::pair< const hash_t, TEfficiency* >( effHash, &effRef ) );
174
175 return &effRef;
176}
EffMap_t m_effMap
The map of histogram names to their pointers.

◆ 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 295 of file AthHistogramming.h.

296{
297 if ( !hist ) {
298 m_msg << MSG::ERROR << "Got a zero pointer to a TEfficiency in bookGetPointer" << endmsg;
299 return 0;
300 }
301 return this->bookGetPointer( *hist, tDir, stream );
302}

◆ 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 83 of file AthHistogramming.cxx.

84{
85 // Modify the name and title according to the prefixes of this classes instance
86 std::string histName(histRef.GetName());
87 const std::string histTitle(histRef.GetTitle());
88 std::string bookingString("");
89
90 this->buildBookingString( bookingString, histName, tDir, stream );
91 histRef.SetTitle((m_histTitlePrefix+histTitle+m_histTitlePostfix).c_str() );
92 histRef.SetName(histName.c_str());
93
94 // Check if the hash for this histName already exists, i.e., if we have a hash collision
95 const hash_t histHash = this->hash(histName);
96 HistMap_t::const_iterator it = m_histMap.find( histHash );
97 if ( it != m_histMap.end() ) // It does exist!
98 {
99 m_msg << MSG::WARNING
100 << "Detected a hash collision. The hash for the histogram with name=" << histName
101 << " already exists and points to a histogram with name=" << it->second->GetName()
102 << " NOT going to book the new histogram and returning a NULL pointer!" << endmsg;
103 return NULL;
104 }
105
106 // Set the new name and title for the histogram, based on the prefixes that the user set for this class instance
107 // Create a clone that has the new name
108
109 // Massage the final string to book things
110
111 // Register the histogram into the THistSvc
112 if ( !((histSvc()->regHist(bookingString, &histRef)).isSuccess()) )
113 {
114 m_msg << MSG::WARNING
115 << "Problem registering histogram with name " << histName
116 << ", name prefix " << m_histNamePrefix
117 << ", title " << histTitle
118 << ", tile prefix " << m_histTitlePrefix
119 << ", and tile postfix " << m_histTitlePostfix
120 << " in " << m_name << "!" << endmsg;
121 return NULL;
122 }
123
124 // Also register it in the local map of string to pointer
125 m_histMap.insert( m_histMap.end(), std::pair< const hash_t, TH1* >( histHash, &histRef ) );
126
127 return &histRef;
128}
HistMap_t m_histMap
The map of histogram names to their pointers.

◆ 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 274 of file AthHistogramming.h.

275{
276 if ( !hist ) {
277 m_msg << MSG::ERROR << "Got a zero pointer to a TH1 in bookGetPointer" << endmsg;
278 return 0;
279 }
280 return this->bookGetPointer( *hist, tDir, stream );
281}

◆ 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 545 of file AthHistogramming.cxx.

550{
551 // Massage the final string to book things
552 if(tDir.empty()) tDir = m_rootDir;
553 size_t pos = histName.rfind('/');
554 if(pos != std::string::npos){
555 tDir+='/';
556 tDir.append(histName, 0,pos);
557 histName.erase(0,pos+1);
558 };
559 if(stream.empty()) stream = m_streamName;
560
561 if(usePrefixPostfix){
562 bookingString = "/"+stream+"/"+tDir+"/"+m_histNamePrefix+histName+m_histNamePostfix;
563 } else {
564 bookingString = "/"+stream+"/"+tDir+"/"+histName;
565 }
566 while(bookingString.find("//") != std::string::npos){
567 this->myReplace(bookingString,"//","/");
568 }
569
570 return;
571}
std::string m_rootDir
Name of the ROOT directory.
std::string m_streamName
Name of the ROOT output stream (file).
void myReplace(std::string &str, const std::string &oldStr, const std::string &newStr)
Helper method to replace sub-string.

◆ 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 51 of file AthHistogramming.cxx.

55{
58 m_rootDir = rootDir;
59 m_histNamePrefix = histNamePrefix;
60 m_histNamePostfix = histNamePostfix;
61 m_histTitlePrefix = histTitlePrefix;
62 m_histTitlePostfix = histTitlePostfix;
63
64 return StatusCode::SUCCESS;
65}
ServiceHandle< ITHistSvc > m_histSvc
Pointer to the THistSvc (event store by default).

◆ declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty ( Gaudi::Property< T, V, H > & hndl,
const SG::VarHandleKeyType &  )
inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

158 {
160 hndl.value(),
161 hndl.documentation());
162
163 }
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)

◆ declareProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( Gaudi::Property< T, V, H > & t)
inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

145 {
146 typedef typename SG::HandleClassifier<T>::type htype;
148 }
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>

◆ detStore()

const ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< Algorithm > >::detStore ( ) const
inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

◆ 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 235 of file AthHistogramming.cxx.

236{
237 // Build a 32 bit hash out of the name
238 const hash_t effHash = this->hash(effName);
239
240 // See if this entry exists in the map
241 EffMap_t::const_iterator it = m_effMap.find( effHash );
242 if ( it == m_effMap.end() ) // It doesn't exist!
243 { // Let's see into the THistSvc if somebody else has registered the TEfficiency...
244
245 // Need to copy the strings as we will massage them from here on
246 std::string effNameCopy = effName;
247 std::string tDirCopy = tDir;
248 std::string streamCopy = stream;
249
250 // Massage the final string to book things
251 std::string bookingString("");
252 this->buildBookingString( bookingString, effNameCopy, tDirCopy, streamCopy ,false);
253
254 TEfficiency* effPointer(NULL);
255 if ( !((histSvc()->getEfficiency(bookingString, effPointer)).isSuccess()) )
256 {
257 // Massage the final string to book things
258 std::string bookingString("");
259 this->buildBookingString( bookingString, effNameCopy, tDirCopy, streamCopy, true );
260
261 if ( !((histSvc()->getEfficiency(bookingString, effPointer)).isSuccess()) )
262 {
263 m_msg << MSG::WARNING
264 << "Problem retrieving the TEfficiency with name (including pre- and post-fixes) "
265 << m_histNamePrefix + effNameCopy + m_histNamePostfix
266 << " or with name " << effNameCopy
267 << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
268 << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
269 return NULL;
270 }
271 // If we get to here, we actually found the TEfficiency in the THistSvc.
272 // So let's add it to the local cache map and return its pointer
273 m_effMap.insert( m_effMap.end(), std::pair< const hash_t, TEfficiency* >( effHash, effPointer ) );
274 return effPointer;
275 }
276 // If we get to here, we actually found the TEfficiency in the THistSvc.
277 // So let's add it to the local cache map and return its pointer
278 m_effMap.insert( m_effMap.end(), std::pair< const hash_t, TEfficiency* >( effHash, effPointer ) );
279 return effPointer;
280 }
281
282 // Return the pointer to the TEfficiency that we got from the local cache map
283 return it->second;
284}
std::pair< StatusCode, TEfficiency * > getEfficiency(ITHistSvc &svc, const std::string &name)

◆ evtStore()

ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< Algorithm > >::evtStore ( )
inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

◆ execute()

StatusCode VertexDecoratorAlg::execute ( )
override

Definition at line 94 of file VertexDecoratorAlg.cxx.

95 {
96 const EventContext &ctx = Gaudi::Hive::currentContext();
97 SG::ReadHandle<xAOD::VertexContainer> vertices(m_vertexInKey, ctx);
98 ATH_CHECK(vertices.isValid());
99
100 SG::ReadHandle<xAOD::PhotonContainer> photonsIn(m_photonsInKey, ctx);
101 ATH_CHECK(photonsIn.isValid());
102
103 SG::ReadHandle<xAOD::ElectronContainer> electronsIn(m_electronsInKey, ctx);
104 ATH_CHECK(electronsIn.isValid());
105 SG::ReadHandle<xAOD::MuonContainer> muonsIn(m_muonsInKey, ctx);
106 ATH_CHECK(muonsIn.isValid());
107 SG::ReadHandle<xAOD::JetContainer> jetsIn(m_jetsInKey, ctx);
108 ATH_CHECK(jetsIn.isValid());
109 SG::ReadHandle<xAOD::EventInfo> eventInfo(m_eventInKey,ctx);
110 ATH_CHECK(eventInfo.isValid());
111
112 SG::WriteDecorHandle<xAOD::VertexContainer, std::vector<ElementLink<xAOD::PhotonContainer>>> dec_photonLinks(m_photonLinksKey, ctx);
113 SG::WriteDecorHandle<xAOD::VertexContainer, std::vector<ElementLink<xAOD::JetContainer>>> dec_jetLinks(m_jetLinksKey, ctx);
114 SG::WriteDecorHandle<xAOD::VertexContainer, std::vector<ElementLink<xAOD::ElectronContainer>>> dec_electronLinks(m_electronLinksKey, ctx);
115 SG::WriteDecorHandle<xAOD::VertexContainer, std::vector<ElementLink<xAOD::MuonContainer>>> dec_muonLinks(m_muonLinksKey, ctx);
116 SG::WriteDecorHandle<xAOD::VertexContainer, std::vector<ElementLink<xAOD::CompositeParticleContainer>>> dec_multiPhotonLinks(m_multiPhotonLinksKey, ctx);
117
118 SG::ReadDecorHandle<xAOD::VertexContainer, float> acc_deltaZ(m_deltaZKey, ctx);
119 SG::ReadDecorHandle<xAOD::VertexContainer, float> acc_deltaPhi(m_deltaPhiKey, ctx);
120 SG::ReadDecorHandle<xAOD::PhotonContainer, float> acc_caloPointingZ(m_caloPointingZKey, ctx);
121 SG::ReadDecorHandle<xAOD::PhotonContainer, float> acc_zCommon(m_zCommonKey, ctx);
122 SG::ReadDecorHandle<xAOD::PhotonContainer, float> acc_zCommonError(m_zCommonErrorKey, ctx);
123
124 // Decorations needed by the GNNTool
125 SG::WriteDecorHandle<xAOD::VertexContainer, int> dec_ntrk(m_mDecor_ntrk, ctx);
126 SG::WriteDecorHandle<xAOD::VertexContainer, float> dec_sumPt(m_mDecor_sumPt,ctx);
127 SG::WriteDecorHandle<xAOD::VertexContainer, float> dec_chi2Over_ndf(m_mDecor_chi2Over_ndf,ctx);
128 SG::WriteDecorHandle<xAOD::VertexContainer, float> dec_z_asym(m_mDecor_z_asym,ctx);
129 SG::WriteDecorHandle<xAOD::VertexContainer, float> dec_weighted_z_asym(m_mDecor_weighted_z_asym,ctx);
130 SG::WriteDecorHandle<xAOD::VertexContainer, float> dec_z_kurt(m_mDecor_weighted_z_kurt,ctx);
131 SG::WriteDecorHandle<xAOD::VertexContainer, float> dec_z_skew(m_mDecor_z_skew,ctx);
132 SG::WriteDecorHandle<xAOD::VertexContainer, float> dec_photon_deltaz(m_mDecor_photon_deltaz,ctx);
133 SG::WriteDecorHandle<xAOD::VertexContainer, float> dec_photon_deltaPhi(m_mDecor_photon_deltaPhi,ctx);
134 SG::WriteDecorHandle<xAOD::VertexContainer, float> dec_actualInterPerXing(m_mDecor_actualInterPerXing,ctx);
135 SG::WriteHandle<xAOD::CompositeParticleContainer> mpHandle(m_multiPhotonsOutKey, ctx);
136
137 auto mpCont = std::make_unique<xAOD::CompositeParticleContainer>();
138 auto mpAux = std::make_unique<xAOD::CompositeParticleAuxContainer>();
139 mpCont->setStore(mpAux.get());
140
141 const xAOD::Vertex* bestVtxForMP = nullptr;
142 ElementLink<xAOD::CompositeParticleContainer> mpLink; // will be valid only if created
143 bool haveMP = false;
144
145 // For MT safety, decoration keys must be initialized and available upfront.
146 const bool haveZCommonDecor = acc_zCommon.isAvailable() && acc_zCommonError.isAvailable();
147 const bool haveCaloPointingDecor = acc_caloPointingZ.isAvailable();
148
149 if (photonsIn->size() > 1 && (!haveZCommonDecor || !haveCaloPointingDecor))
150 {
151 ATH_MSG_ERROR("photonsIn has size > 1 but required decorations are missing: "
152 << m_zCommonKey.key() << " and/or " << m_zCommonErrorKey.key()
153 << " and/or " << m_caloPointingZKey.key());
154 return StatusCode::FAILURE;
155 }
156
157 if (photonsIn->size() > 1)
158 {
159 // Combined pointing Z (weighted)
160 double sumW = 0.0;
161 double sumWZ = 0.0;
162
163 TLorentzVector p4sum;
164 int nUsed = 0;
165 for (const xAOD::Photon* ph : *photonsIn)
166 {
167 const float zCommon = acc_zCommon(*ph);
168 const float zCaloPointing = acc_caloPointingZ(*ph);
169 // Value-based selection only: prefer zCommon, then fallback to caloPointing value.
170 const bool useZCommonValue = std::isfinite(zCommon);
171 const float z = useZCommonValue ? zCommon : zCaloPointing;
172 if (!std::isfinite(z)) {
173 ATH_MSG_ERROR("Non-finite photon pointing values found in both "
174 << m_zCommonKey.key() << " and " << m_caloPointingZKey.key());
175 return StatusCode::FAILURE;
176 }
177
178 // Keep weighting consistent with the z source:
179 // use zCommonError only when z comes from zCommon.
180 double w = 1.0;
181 if (useZCommonValue) {
182 const float s = acc_zCommonError(*ph);
183 if (std::isfinite(s) && s > 0.0f) w = 1.0 / (double(s) * double(s));
184 }
185
186 sumW += w;
187 sumWZ += w * z;
188 p4sum += ph->p4();
189 ++nUsed;
190 }
191
192 if (nUsed > 1 && sumW > 0.0)
193 {
194 const float zPoint = sumWZ / sumW;
195
196 // Choose vertex closest in z to combined pointing
197 float bestAbsDZ = 1e30;
198 for (const xAOD::Vertex* vtx : *vertices)
199 {
200 if (vtx->vertexType() == xAOD::VxType::NoVtx) continue;
201
202 const float dz = zPoint - vtx->z();
203 const float adz = std::abs(dz);
204 if (adz < bestAbsDZ)
205 {
206 bestAbsDZ = adz;
207 bestVtxForMP = vtx;
208 }
209 }
210
211 if (bestVtxForMP)
212 {
213 // Create exactly one CompositeParticle node
214 auto* mp = new xAOD::CompositeParticle();
215 mpCont->push_back(mp);
216
217 mp->setP4(p4sum);
218
219 SG::AuxElement::Decorator<float> dec_mp_deltaZ("deltaZ");
220 SG::AuxElement::Decorator<float> dec_mp_deltaPhi("deltaPhi");
221 SG::AuxElement::Decorator<int> dec_mp_nPhotons("nPhotons");
222 SG::AuxElement::Decorator<float> dec_mp_zPointing("zPointing");
223
224 const float deltaZ = zPoint - bestVtxForMP->z();
225 dec_mp_deltaZ(*mp) = deltaZ;
226 dec_mp_zPointing(*mp) = zPoint;
227 dec_mp_nPhotons(*mp) = static_cast<int>(photonsIn->size());
228
229 float dphi = 0.0f;
230 if (acc_deltaPhi.isAvailable())
231 {
232 dphi = acc_deltaPhi(*bestVtxForMP);
233 if (!std::isfinite(dphi)) dphi = 0.0f;
234 }
235 dec_mp_deltaPhi(*mp) = dphi;
236 mpLink.setElement(mp);
237 mpLink.setStorableObject(*mpCont, true);
238
239 haveMP = true;
240 }
241 }
242 }
243
244 // Record the MultiPhotons container (even if empty)
245 ATH_CHECK(mpHandle.record(std::move(mpCont), std::move(mpAux)));
246
247 std::map< const xAOD::Vertex*, std::vector<ElementLink<xAOD::JetContainer>> > jetsInVertex;
248 std::map< const xAOD::Jet*, std::map< const xAOD::Vertex*, int> > jetVertexPt;
249
250 // initialize jet-vertex map
251 for(const xAOD::Vertex *vertex : *vertices){
252 jetsInVertex[vertex] = {};
253 for(const xAOD::Jet *jet : *jetsIn){
254 jetVertexPt[jet][vertex] = 0;
255 }
256 }
257
258 // pre-fill the jet-vertex map
259 for(const xAOD::Jet* jet : *jetsIn){
260
261 // for each jet, calculate the track pT associated to each vertex
262 std::vector<const xAOD::TrackParticle*> ghostTracks = jet->getAssociatedObjects<xAOD::TrackParticle >(xAOD::JetAttribute::GhostTrack);
263 for(const xAOD::TrackParticle* jtrk : ghostTracks){
264 if( !jtrk ) continue;
265 auto jetTrackVertex = m_trkVtxAssociationTool->getUniqueMatchVertexLink(*jtrk, *vertices);
266 if(jetTrackVertex) jetVertexPt[jet][*jetTrackVertex] += jtrk->pt();
267 }
268
269 // find vertex with the largest fraction of jet track pt
270 float maxPtFrac = -1;
271 const xAOD::Vertex* uniqueVertexAddress = nullptr;
272 for (const xAOD::Vertex *vertex : *vertices) {
273 if (vertex->vertexType() == xAOD::VxType::NoVtx) continue;
274 if(jetVertexPt[jet][vertex] > maxPtFrac){
275 maxPtFrac = jetVertexPt[jet][vertex];
276 uniqueVertexAddress = vertex;
277 }
278 }
279
280 // add jet to that vertex's vector of links
281 ElementLink<xAOD::JetContainer> jetLink;
282 jetLink.setElement(jet);
283 jetLink.setStorableObject(*jetsIn.ptr(), true);
284 jetsInVertex[uniqueVertexAddress].push_back(jetLink);
285 }
286
287 for (const xAOD::Vertex *vertex : *vertices)
288 {
289 if (vertex->vertexType() == xAOD::VxType::NoVtx)
290 continue;
291
292 if (vertex->nTrackParticles() < 2)
293 continue;
294
295 dec_actualInterPerXing(*vertex) = eventInfo->actualInteractionsPerCrossing();
296
297 // variables for calculation of delta Z asymmetry and delta d asymmetry
298 float z_asym = 0;
299 float sumDZ = 0;
300 float deltaZ = 0;
301 float modsumDZ = 0;
302 float weighted_sumDZ = 0;
303 float weighted_deltaZ = 0;
304 float weighted_modsumDZ = 0;
305 float weighted_z_asym = 0;
306
307 // make vector
308 std::vector<float> track_deltaZ;
309
310 for (size_t i = 0; i < vertex->nTrackParticles(); i++) {
311
312 const xAOD::TrackParticle *trackTmp = vertex->trackParticle(i);
313
314 if(!trackTmp) continue;
315
316 deltaZ = trackTmp->z0() + trackTmp->vz() - vertex->z();
317 track_deltaZ.push_back(deltaZ);
318 // get the track weight for each track to get the deltaZ/trk_weight
319 float trk_weight = vertex->trackWeight(i);
320 weighted_deltaZ = deltaZ * trk_weight;
321 // sum of delta z
322 sumDZ += deltaZ;
323 modsumDZ += std::abs(deltaZ);
324 weighted_sumDZ += weighted_deltaZ;
325 weighted_modsumDZ += std::abs(weighted_deltaZ);
326 } // end loop over tracks
327
328 if (modsumDZ > 0) {
329 z_asym = sumDZ / modsumDZ;
330 }
331 if (weighted_modsumDZ > 0) {
332 weighted_z_asym = weighted_sumDZ / weighted_modsumDZ;
333 }
334
335 const float number_tracks = track_deltaZ.size(); // get number of tracks
336 const float mean_Dz =
337 number_tracks > 0 ? sumDZ / number_tracks : 0.F; // calculate average
338
339 float z_skew = 0; // skewness of DeltaZ asymmetry
340 float z_kurt = 0; // Kurtosis of DeltaZ asymmetry
341 float z_var = 0; // variance of DeltaZ
342
343 for (auto i : track_deltaZ)
344 {
345 float z_zbar = (i - mean_Dz);
346 z_var += std::pow(z_zbar, 2);
347 z_skew += std::pow(z_zbar, 3);
348 z_kurt += std::pow(z_zbar, 4);
349 }
350 if (number_tracks > 1 && z_var > 0) {
351 z_var /= (number_tracks - 1);
352 float z_sd = std::sqrt(z_var);
353 const float skew_denom = (number_tracks - 1) * std::pow(z_sd, 3);
354 const float kurt_denom = (number_tracks - 1) * std::pow(z_sd, 4);
355 if (std::isfinite(skew_denom) && skew_denom != 0.F) {
356 z_skew /= skew_denom;
357 } else {
358 z_skew = 0.F;
359 }
360 if (std::isfinite(kurt_denom) && kurt_denom != 0.F) {
361 z_kurt /= kurt_denom;
362 } else {
363 z_kurt = 0.F;
364 }
365 }
366 else
367 {
368 ATH_MSG_WARNING("z momenta are NaN: setting to zero");
369 z_skew = 0.;
370 z_kurt = 0.;
371 }
372
373 dec_ntrk(*vertex) = number_tracks;
374
375 if(!dec_sumPt.isAvailable()){
376 dec_sumPt(*vertex) = xAOD::PVHelpers::getVertexSumPt(vertex, 1, false);
377 }
378 const float numberDoF = vertex->numberDoF();
379 dec_chi2Over_ndf(*vertex) =
380 numberDoF > 0.F ? vertex->chiSquared() / numberDoF : 0.F;
381 dec_z_asym(*vertex) = z_asym;
382 dec_weighted_z_asym(*vertex) = weighted_z_asym;
383 dec_z_kurt(*vertex) = z_kurt;
384 dec_z_skew(*vertex) = z_skew;
385
386 if (acc_deltaZ.isAvailable() && photonsIn->size() > 0 ) {
387 //protect against rare NaNs before assigning decorator: setting to 0 (-999 cause NaNs)
388 if (!std::isfinite(acc_deltaZ(*vertex))) {
389 ATH_MSG_WARNING("photon deltaPhi is NaN: setting to -1!");
390 dec_photon_deltaz(*vertex) = -1;
391 }
392 else{
393 dec_photon_deltaz(*vertex) = acc_deltaZ(*vertex);
394 }
395 }
396 else{
397 dec_photon_deltaz(*vertex) = -1;
398 }
399 if (acc_deltaPhi.isAvailable() && photonsIn->size() > 0) {
400
401 if (!std::isfinite(acc_deltaPhi(*vertex))) {
402 ATH_MSG_WARNING("photon deltaPhi is NaN: setting to 0!");
403 dec_photon_deltaPhi(*vertex) = -1;
404 }
405 else{
406 dec_photon_deltaPhi(*vertex) = acc_deltaPhi(*vertex);
407 }
408 }
409 else{
410 dec_photon_deltaPhi(*vertex) = -1;
411 }
412
413 // associate objects to vertices
414 std::vector<ElementLink<xAOD::ElectronContainer>> electronLinks;
415 for(const xAOD::Electron* electron : *electronsIn){
416 const auto *id_trk = xAOD::EgammaHelpers::getOriginalTrackParticle(electron);
417 if(!id_trk) continue;
418 auto eleVertex = m_trkVtxAssociationTool->getUniqueMatchVertexLink(*id_trk, *vertices);
419 if(!eleVertex) continue;
420 ElementLink<xAOD::ElectronContainer> elLink;
421 if(*eleVertex == vertex){
422 elLink.setElement(electron);
423 elLink.setStorableObject(*electronsIn.ptr(), true);
424 electronLinks.push_back(elLink);
425 }
426 }
427 dec_electronLinks(*vertex) = electronLinks;
428
429 std::vector<ElementLink<xAOD::PhotonContainer>> photonLinks;
430 for(const xAOD::Photon* photon : *photonsIn){
431 ElementLink<xAOD::PhotonContainer> phLink;
432 phLink.setElement(photon);
433 phLink.setStorableObject(*photonsIn.ptr(), true);
434 photonLinks.push_back(phLink);
435 }
436 dec_photonLinks(*vertex) = photonLinks;
437
438 // multi-photon link
439 std::vector<ElementLink<xAOD::CompositeParticleContainer>> mpLinks;
440 if (haveMP && vertex == bestVtxForMP)
441 {
442 mpLinks.push_back(mpLink);
443 }
444 dec_multiPhotonLinks(*vertex) = mpLinks;
445
446 // for jets, use prefilled map
447 dec_jetLinks(*vertex) = jetsInVertex[vertex];
448
449 std::vector<ElementLink<xAOD::MuonContainer>> muonLinks;
450 for(const xAOD::Muon* muon : *muonsIn){
451 const auto *tp = muon->trackParticle(xAOD::Muon::InnerDetectorTrackParticle);
452 if(!tp) continue;
453 try{
454 auto muonVertex = m_trkVtxAssociationTool->getUniqueMatchVertexLink(*tp, *vertices);
455 if(!muonVertex) continue;
456 ElementLink<xAOD::MuonContainer> muonLink;
457 if(*muonVertex == vertex){
458 muonLink.setElement(muon);
459 muonLink.setStorableObject(*muonsIn.ptr(), true);
460 muonLinks.push_back(muonLink);
461 }
462 }catch(...) {
463 ATH_MSG_DEBUG("Skipping muon as the track is not associated to any PV ");
464 ATH_MSG_DEBUG("Muon pT, eta = " << muon->pt() << " " << muon->eta());
465 }
466
467 }
468 dec_muonLinks(*vertex) = muonLinks;
469
470 // Finally, decorate the vertices with the GNN score
471 m_gnnTool->decorate(*vertex);
472 }
473
474 return StatusCode::SUCCESS;
475 }
#define ATH_CHECK
Evaluate an expression and check for errors.
#define ATH_MSG_ERROR(x)
#define ATH_MSG_WARNING(x)
#define ATH_MSG_DEBUG(x)
#define z
SG::WriteDecorHandleKey< xAOD::VertexContainer > m_mDecor_photon_deltaPhi
SG::WriteDecorHandleKey< xAOD::VertexContainer > m_mDecor_z_skew
SG::WriteDecorHandleKey< xAOD::VertexContainer > m_muonLinksKey
SG::WriteDecorHandleKey< xAOD::VertexContainer > m_mDecor_sumPt
SG::WriteHandleKey< xAOD::CompositeParticleContainer > m_multiPhotonsOutKey
SG::WriteDecorHandleKey< xAOD::VertexContainer > m_mDecor_chi2Over_ndf
SG::ReadDecorHandleKey< xAOD::VertexContainer > m_deltaZKey
SG::WriteDecorHandleKey< xAOD::VertexContainer > m_photonLinksKey
SG::ReadHandleKey< xAOD::MuonContainer > m_muonsInKey
SG::WriteDecorHandleKey< xAOD::VertexContainer > m_mDecor_weighted_z_asym
SG::ReadDecorHandleKey< xAOD::PhotonContainer > m_zCommonKey
SG::WriteDecorHandleKey< xAOD::VertexContainer > m_electronLinksKey
SG::ReadHandleKey< xAOD::EventInfo > m_eventInKey
SG::WriteDecorHandleKey< xAOD::VertexContainer > m_mDecor_ntrk
SG::WriteDecorHandleKey< xAOD::VertexContainer > m_mDecor_actualInterPerXing
ToolHandle< CP::TrackVertexAssociationTool > m_trkVtxAssociationTool
SG::ReadDecorHandleKey< xAOD::VertexContainer > m_deltaPhiKey
SG::WriteDecorHandleKey< xAOD::VertexContainer > m_mDecor_weighted_z_kurt
SG::ReadDecorHandleKey< xAOD::PhotonContainer > m_zCommonErrorKey
SG::WriteDecorHandleKey< xAOD::VertexContainer > m_mDecor_photon_deltaz
SG::ReadDecorHandleKey< xAOD::PhotonContainer > m_caloPointingZKey
SG::WriteDecorHandleKey< xAOD::VertexContainer > m_multiPhotonLinksKey
SG::WriteDecorHandleKey< xAOD::VertexContainer > m_mDecor_z_asym
SG::ReadHandleKey< xAOD::VertexContainer > m_vertexInKey
SG::ReadHandleKey< xAOD::PhotonContainer > m_photonsInKey
SG::WriteDecorHandleKey< xAOD::VertexContainer > m_jetLinksKey
SG::ReadHandleKey< xAOD::ElectronContainer > m_electronsInKey
SG::ReadHandleKey< xAOD::JetContainer > m_jetsInKey
float z0() const
Returns the parameter.
float vz() const
The z origin for the parameters.
float z() const
Returns the z position.
const xAOD::TrackParticle * getOriginalTrackParticle(const xAOD::Electron *el)
Helper function for getting the "Original" Track Particle (i.e before GSF) via the electron.
float getVertexSumPt(const xAOD::Vertex *vertex, int power=1, bool useAux=true)
Loop over track particles associated with vertex and return scalar sum of pT^power in GeV (from auxda...
@ NoVtx
Dummy vertex. TrackParticle was not used in vertex fit.
Jet_v1 Jet
Definition of the current "jet version".
TrackParticle_v1 TrackParticle
Reference the current persistent version:
Vertex_v1 Vertex
Define the latest version of the vertex class.
CompositeParticle_v1 CompositeParticle
Define the latest version of the composite particle class.
Muon_v1 Muon
Reference the current persistent version:
Photon_v1 Photon
Definition of the current "egamma version".
Electron_v1 Electron
Definition of the current "egamma version".

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Algorithm > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase & ExtraDeps)
protectedinherited

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.

51{
52 // If we didn't find any symlinks to add, just return the collection
53 // from the base class. Otherwise, return the extended collection.
54 if (!m_extendedExtraObjects.empty()) {
56 }
57 return Algorithm::extraOutputDeps();
58}
DataObjIDColl m_extendedExtraObjects

◆ 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 477 of file AthHistogramming.cxx.

478{
479 // Build a 32 bit hash out of the name
480 const hash_t graphHash = this->hash(graphName);
481
482 // See if this entry exists in the map
483 GraphMap_t::const_iterator it = m_graphMap.find( graphHash );
484 if ( it == m_graphMap.end() ) // It doesn't exist!
485 { // Let's see into the THistSvc if somebody else has registered the TGraph...
486
487 // Need to copy the strings as we will massage them from here on
488 std::string graphNameCopy = graphName;
489 std::string tDirCopy = tDir;
490 std::string streamCopy = stream;
491
492 // Massage the final string to book things
493 std::string bookingString("");
494 this->buildBookingString( bookingString, graphNameCopy, tDirCopy, streamCopy, true);
495
496 TGraph* graphPointer(NULL);
497 if ( !((histSvc()->getGraph(bookingString, graphPointer)).isSuccess()) )
498 {
499 // Massage the final string to book things
500 std::string bookingString("");
501 this->buildBookingString( bookingString, graphNameCopy, tDirCopy, streamCopy, false );
502
503 if ( !((histSvc()->getGraph(bookingString, graphPointer)).isSuccess()) )
504 {
505 m_msg << MSG::WARNING
506 << "Problem retrieving the TGraph with name (including pre- and post-fixes) "
507 << m_histNamePrefix + graphNameCopy + m_histNamePostfix
508 << " or with name " << graphNameCopy
509 << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
510 << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
511 return NULL;
512 }
513 // If we get to here, we actually found the TGraph in the THistSvc.
514 // So let's add it to the local cache map and return its pointer
515 m_graphMap.insert( m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphPointer ) );
516 return graphPointer;
517 }
518 // If we get to here, we actually found the TGraph in the THistSvc.
519 // So let's add it to the local cache map and return its pointer
520 m_graphMap.insert( m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphPointer ) );
521 return graphPointer;
522 }
523
524
525 // Return the pointer to the TGraph that we got from the local cache map
526 return it->second;
527}
std::pair< StatusCode, TGraph * > getGraph(ITHistSvc &svc, const std::string &name)

◆ hash()

AthHistogramming::hash_t AthHistogramming::hash ( const std::string & histName) const
inlineprivateinherited

Method to calculate a 32-bit hash from a string.

Definition at line 430 of file AthHistogramming.h.

431{
432 const uint64_t hash64 = CxxUtils::crc64( histName );
433 return (hash_t)(hash64 & 0xFFFFFFFF);
434}
uint64_t crc64(const CRCTable &table, const char *data, size_t data_len)
Find the CRC-64 of a string,.
Definition crc64.cxx:696
std::uint64_t hash64(const void *data, std::size_t size)
Passthrough to XXH3_64bits.
Definition XXH.cxx:9

◆ 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 183 of file AthHistogramming.cxx.

184{
185 // Build a 32 bit hash out of the name
186 const hash_t histHash = this->hash(histName);
187
188 // See if this entry exists in the map
189 HistMap_t::const_iterator it = m_histMap.find( histHash );
190 if ( it == m_histMap.end() ) // It doesn't exist!
191 { // Let's see into the THistSvc if somebody else has registered the histogram...
192
193 // Need to copy the strings as we will massage them from here on
194 std::string histNameCopy = histName;
195 std::string tDirCopy = tDir;
196 std::string streamCopy = stream;
197
198 // Massage the final string to book things
199 std::string bookingString("");
200 this->buildBookingString( bookingString, histNameCopy, tDirCopy, streamCopy ,false);
201
202 TH1* histPointer(NULL);
203 if ( !((histSvc()->getHist(bookingString, histPointer)).isSuccess()) )
204 {
205 // Massage the final string to book things
206 std::string bookingString("");
207 this->buildBookingString( bookingString, histNameCopy, tDirCopy, streamCopy, true );
208
209 if ( !((histSvc()->getHist(bookingString, histPointer)).isSuccess()) )
210 {
211 m_msg << MSG::WARNING
212 << "Problem retrieving the histogram with name (including pre- and post-fixes) "
213 << m_histNamePrefix + histNameCopy + m_histNamePostfix
214 << " or with name " << histNameCopy
215 << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
216 << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
217 return NULL;
218 }
219 // If we get to here, we actually found the histogram in the THistSvc.
220 // So let's add it to the local cache map and return its pointer
221 m_histMap.insert( m_histMap.end(), std::pair< const hash_t, TH1* >( histHash, histPointer ) );
222 return histPointer;
223 }
224 // If we get to here, we actually found the histogram in the THistSvc.
225 // So let's add it to the local cache map and return its pointer
226 m_histMap.insert( m_histMap.end(), std::pair< const hash_t, TH1* >( histHash, histPointer ) );
227 return histPointer;
228 }
229
230
231 // Return the pointer to the histogram that we got from the local cache map
232 return it->second;
233}
std::pair< StatusCode, TH1 * > getHist(ITHistSvc &svc, const std::string &name, size_t index=0)

◆ 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 371 of file AthHistogramming.h.

372{
373 // Get the TH1 pointer
374 TH1* th1Pointer = this->hist(histName, tDir, stream);
375 if ( !th1Pointer )
376 {
377 m_msg << MSG::ERROR
378 << "Cannot get a 2-d histogram with name " << histName
379 << "... will probably seg-fault!" << endmsg;
380 return NULL;
381 }
382 // If the TH1 pointer is valid, simply return the dynamic_cast
383 return dynamic_cast<TH2*>( th1Pointer );
384}

◆ 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 388 of file AthHistogramming.h.

389{
390 // Get the TH1 pointer
391 TH1* th1Pointer = this->hist(histName, tDir, stream);
392 if ( !th1Pointer )
393 {
394 m_msg << MSG::ERROR
395 << "Cannot get a 3-d histogram with name " << histName
396 << "... will probably seg-fault!" << endmsg;
397 return NULL;
398 }
399 // If the TH1 pointer is valid, simply return the dynamic_cast
400 return dynamic_cast<TH3*>( th1Pointer );
401}

◆ histSvc()

const ServiceHandle< ITHistSvc > & AthHistogramAlgorithm::histSvc ( ) const
inlineinherited

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.

114{
115 return m_histSvc;
116}
ServiceHandle< ITHistSvc > m_histSvc
Default constructor: AthHistogramAlgorithm();.

◆ initialize()

StatusCode VertexDecoratorAlg::initialize ( )
override

Definition at line 25 of file VertexDecoratorAlg.cxx.

26 {
27 ATH_CHECK(m_vertexInKey.initialize());
28 ATH_CHECK(m_eventInKey.initialize());
29 ATH_CHECK(m_photonsInKey.initialize());
30 ATH_CHECK(m_electronsInKey.initialize());
31 ATH_CHECK(m_muonsInKey.initialize());
32 ATH_CHECK(m_jetsInKey.initialize());
33 ATH_CHECK(m_multiPhotonsOutKey.initialize());
34
35 const std::string baseName = m_vertexInKey.key();
36
37 // WriteHandleKeys
38 m_photonLinksKey = baseName + ".photonLinks";
39 m_jetLinksKey = baseName + ".jetLinks";
40 m_electronLinksKey = baseName + ".electronLinks";
41 m_muonLinksKey = baseName + ".muonLinks";
42 m_mDecor_ntrk = baseName + "." + m_mDecor_ntrk.key();
43 m_mDecor_sumPt = baseName + "." + m_mDecor_sumPt.key();
44 m_mDecor_chi2Over_ndf = baseName + "." + m_mDecor_chi2Over_ndf.key();
45 m_mDecor_z_asym = baseName + "." + m_mDecor_z_asym.key();
48 m_mDecor_z_skew = baseName + "." + m_mDecor_z_skew.key();
49 m_mDecor_photon_deltaz = baseName + "." + m_mDecor_photon_deltaz.key();
52 m_multiPhotonLinksKey = baseName + ".multiPhotonLinks";
53
54 ATH_CHECK(m_multiPhotonLinksKey.initialize());
55 ATH_CHECK(m_photonLinksKey.initialize());
56 ATH_CHECK(m_jetLinksKey.initialize());
57 ATH_CHECK(m_electronLinksKey.initialize());
58 ATH_CHECK(m_muonLinksKey.initialize());
59 ATH_CHECK(m_mDecor_ntrk.initialize());
60 ATH_CHECK(m_mDecor_sumPt.initialize());
61 ATH_CHECK(m_mDecor_chi2Over_ndf.initialize());
62 ATH_CHECK(m_mDecor_z_asym.initialize());
65 ATH_CHECK(m_mDecor_z_skew.initialize());
69
70 // ReadHandleKeys
71 m_deltaZKey = baseName + ".deltaZ";
72 m_deltaPhiKey = baseName + ".deltaPhi";
73 ATH_CHECK(m_deltaZKey.initialize());
74 ATH_CHECK(m_deltaPhiKey.initialize());
75
76 // additional ReadHandleKeys to declare dependencies to the scheduler
77 std::string photonBaseName = m_photonsInKey.key();
78
79 m_caloPointingZKey = photonBaseName + ".caloPointingZ";
80 m_zCommonKey = photonBaseName + ".zCommon";
81 m_zCommonErrorKey = photonBaseName + ".zCommonError";
82
83 ATH_CHECK(m_caloPointingZKey.initialize());
84 ATH_CHECK(m_zCommonKey.initialize());
85 ATH_CHECK(m_zCommonErrorKey.initialize());
86
87 // Tools
88 ATH_CHECK(m_gnnTool.retrieve());
90
91 return StatusCode::SUCCESS;
92 }

◆ inputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< Algorithm > >::inputHandles ( ) const
overridevirtualinherited

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.

◆ msg()

MsgStream & AthCommonMsg< Algorithm >::msg ( ) const
inlineinherited

Definition at line 24 of file AthCommonMsg.h.

24 {
25 return this->msgStream();
26 }

◆ msgLvl()

bool AthCommonMsg< Algorithm >::msgLvl ( const MSG::Level lvl) const
inlineinherited

Definition at line 30 of file AthCommonMsg.h.

30 {
31 return this->msgLevel(lvl);
32 }

◆ 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 575 of file AthHistogramming.cxx.

578{
579 size_t pos = 0;
580 while((pos = str.find(oldStr, pos)) != std::string::npos)
581 {
582 str.replace(pos, oldStr.length(), newStr);
583 pos += newStr.length();
584 }
585}

◆ outputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< Algorithm > >::outputHandles ( ) const
overridevirtualinherited

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.

◆ renounce()

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 > AthCommonDataStore< AthCommonMsg< Algorithm > >::renounce ( T & h)
inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

381 {
382 h.renounce();
384 }
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)

◆ renounceArray()

void AthCommonDataStore< AthCommonMsg< Algorithm > >::renounceArray ( SG::VarHandleKeyArray & handlesArray)
inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364 {
366 }

◆ 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.

76{
77 // ---- stolen from GaudiKernel/Algorithm::sysInitialize -------
78 // Bypass the initialization if the algorithm
79 // has already been initialized.
80 if ( Gaudi::StateMachine::INITIALIZED <= FSMState() ) return StatusCode::SUCCESS;
81
82 // Set the Algorithm's properties
83 bindPropertiesTo( serviceLocator()->getOptsSvc() );
84
85 // Bypass the initialization if the algorithm is disabled.
86 // Need to do this after bindPropertiesTo.
87 if ( !isEnabled( ) ) return StatusCode::SUCCESS;
88
89 // ---- stolen from GaudiKernel/Algorithm::sysInitialize ------- END ---
90
91
92 // Get the THistSvc
93 ATH_CHECK ( histSvc().retrieve() );
94
95 // Configure the underlying AthHistogramming helper
100
101 // Print some setup information into the log file
102 ATH_MSG_DEBUG ("Initializing " << name() << "...");
103 ATH_MSG_DEBUG (" using THistService = " << m_histSvc );
104 ATH_MSG_DEBUG (" using RootStreamName = " << m_prefix );
105 ATH_MSG_DEBUG (" using RootDirName = " << m_rootDir );
106 ATH_MSG_DEBUG (" using HistNamePrefix = " << m_histNamePrefix );
107 ATH_MSG_DEBUG (" using HistNamePostfix = " << m_histNamePostfix );
108 ATH_MSG_DEBUG (" using HistTitlePrefix = " << m_histTitlePrefix );
109 ATH_MSG_DEBUG (" using HistTitlePostfix = " << m_histTitlePostfix );
110
111
112 // re-direct to base class...
114}
virtual StatusCode sysInitialize() override
Override sysInitialize.
const ServiceHandle< ITHistSvc > & histSvc() const
The standard THistSvc (for writing histograms and TTrees and more to a root file) Returns (kind of) a...
std::string m_histNamePrefix
The prefix for the histogram THx name.
std::string m_prefix
Name of the ROOT output stream (file).
std::string m_histNamePostfix
The postfix for the histogram THx name.
std::string m_histTitlePostfix
The postfix for the histogram THx title.
std::string m_rootDir
Name of the ROOT directory.
std::string m_histTitlePrefix
The prefix for the histogram THx title.
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.

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Algorithm > >::sysStart ( )
overridevirtualinherited

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 363 of file AthHistogramming.cxx.

364{
365 // Build a 32 bit hash out of the name
366 const hash_t treeHash = this->hash(treeName);
367
368 // See if this entry exists in the map
369 TreeMap_t::const_iterator it = m_treeMap.find( treeHash );
370 if ( it == m_treeMap.end() ) // It doesn't exist!
371 { // Let's see into the THistSvc if somebody else has registered the TTree...
372
373 // Need to copy the strings as we will massage them from here on
374 std::string treeNameCopy = treeName;
375 std::string tDirCopy = tDir;
376 std::string streamCopy = stream;
377
378 // Massage the final string to book things
379 std::string bookingString("");
380 this->buildBookingString( bookingString, treeNameCopy, tDirCopy, streamCopy );
381
382 TTree* treePointer(NULL);
383 if ( !((histSvc()->getTree(bookingString, treePointer)).isSuccess()) )
384 {
385 m_msg << MSG::WARNING
386 << "Problem retrieving the TTree with name " << treeNameCopy
387 << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
388 << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
389 return NULL;
390 }
391 // If we get to here, we actually found the TTree in the THistSvc.
392 // So let's add it to the local cache map and return its pointer
393 m_treeMap.insert( m_treeMap.end(), std::pair< const hash_t, TTree* >( treeHash, treePointer ) );
394 return treePointer;
395 }
396
397 // Return the pointer to the TTree that we got from the local cache map
398 return it->second;
399}
std::pair< StatusCode, TTree * > getTree(ITHistSvc &svc, const std::string &name)

◆ updateVHKA()

void AthCommonDataStore< AthCommonMsg< Algorithm > >::updateVHKA ( Gaudi::Details::PropertyBase & )
inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

308 {
309 // debug() << "updateVHKA for property " << p.name() << " " << p.toString()
310 // << " size: " << m_vhka.size() << endmsg;
311 for (auto &a : m_vhka) {
313 for (auto k : keys) {
314 k->setOwner(this);
315 }
316 }
317 }
std::vector< SG::VarHandleKeyArray * > m_vhka

Member Data Documentation

◆ m_caloPointingZKey

SG::ReadDecorHandleKey<xAOD::PhotonContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_caloPointingZKey
private
Initial value:
{
this, "caloPointingZKey", "", "" }

Definition at line 105 of file VertexDecoratorAlg.h.

105 {
106 this, "caloPointingZKey", "", "" };

◆ m_deltaPhiKey

SG::ReadDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_deltaPhiKey
private
Initial value:
{
this, "deltaPhiKey", "", "" }

Definition at line 100 of file VertexDecoratorAlg.h.

100 {
101 this, "deltaPhiKey", "", "" };

◆ m_deltaZKey

SG::ReadDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_deltaZKey
private
Initial value:
{
this, "deltaZKey", "", "" }

Definition at line 98 of file VertexDecoratorAlg.h.

98 {
99 this, "deltaZKey", "", "" };

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Algorithm > >::m_detStore
privateinherited

Pointer to StoreGate (detector store by default).

Definition at line 393 of file AthCommonDataStore.h.

◆ m_effMap

EffMap_t AthHistogramming::m_effMap
privateinherited

The map of histogram names to their pointers.

Definition at line 210 of file AthHistogramming.h.

◆ m_electronLinksKey

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_electronLinksKey
private
Initial value:
{
this, "electronLinks", "", "" }

Definition at line 87 of file VertexDecoratorAlg.h.

87 {
88 this, "electronLinks", "", "" };

◆ m_electronsInKey

SG::ReadHandleKey<xAOD::ElectronContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_electronsInKey
private
Initial value:
{
this, "electronsIn", "Electrons", "containerName to read"}

Definition at line 53 of file VertexDecoratorAlg.h.

53 {
54 this, "electronsIn", "Electrons", "containerName to read"};

◆ m_eventInKey

SG::ReadHandleKey<xAOD::EventInfo> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_eventInKey
private
Initial value:
{
this, "eventIn", "EventInfo", "containerName to read"}

Definition at line 51 of file VertexDecoratorAlg.h.

51 {
52 this, "eventIn", "EventInfo", "containerName to read"};

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Algorithm > >::m_evtStore
privateinherited

Pointer to StoreGate (event store by default).

Definition at line 390 of file AthCommonDataStore.h.

◆ m_extendedExtraObjects

DataObjIDColl AthAlgorithm::m_extendedExtraObjects
privateinherited

Definition at line 79 of file AthAlgorithm.h.

◆ m_gnnTool

ToolHandle<GNNTool> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_gnnTool
private
Initial value:
{
this, "gnnTool",
"InDetGNNHardScatterSelection/GNNTool",
"GNN tool to use"}

Definition at line 118 of file VertexDecoratorAlg.h.

118 {
119 this, "gnnTool",
120 "InDetGNNHardScatterSelection/GNNTool",
121 "GNN tool to use"};

◆ m_graphMap

GraphMap_t AthHistogramming::m_graphMap
privateinherited

The map of TGraph names to their pointers.

Definition at line 224 of file AthHistogramming.h.

◆ m_histMap

HistMap_t AthHistogramming::m_histMap
privateinherited

The map of histogram names to their pointers.

Definition at line 203 of file AthHistogramming.h.

◆ m_histNamePostfix

std::string AthHistogramAlgorithm::m_histNamePostfix
privateinherited

The postfix for the histogram THx name.

Definition at line 97 of file AthHistogramAlgorithm.h.

◆ m_histNamePrefix

std::string AthHistogramAlgorithm::m_histNamePrefix
privateinherited

The prefix for the histogram THx name.

Definition at line 94 of file AthHistogramAlgorithm.h.

◆ m_histSvc

ServiceHandle<ITHistSvc> AthHistogramAlgorithm::m_histSvc
privateinherited

Default constructor: AthHistogramAlgorithm();.

a handle on the Hist/TTree registration service

Definition at line 83 of file AthHistogramAlgorithm.h.

◆ m_histTitlePostfix

std::string AthHistogramAlgorithm::m_histTitlePostfix
privateinherited

The postfix for the histogram THx title.

Definition at line 103 of file AthHistogramAlgorithm.h.

◆ m_histTitlePrefix

std::string AthHistogramAlgorithm::m_histTitlePrefix
privateinherited

The prefix for the histogram THx title.

Definition at line 100 of file AthHistogramAlgorithm.h.

◆ m_jetLinksKey

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_jetLinksKey
private
Initial value:
{
this, "jetLinks", "", "" }

Definition at line 85 of file VertexDecoratorAlg.h.

85 {
86 this, "jetLinks", "", "" };

◆ m_jetsInKey

SG::ReadHandleKey<xAOD::JetContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_jetsInKey
private
Initial value:
{
this, "jetsIn", "AntiKt4EMTopoJets", "containerName to read"}

Definition at line 57 of file VertexDecoratorAlg.h.

57 {
58 this, "jetsIn", "AntiKt4EMTopoJets", "containerName to read"};

◆ m_mDecor_actualInterPerXing

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_mDecor_actualInterPerXing
private
Initial value:
{
this, "decor_actualInterPerXing", "actualIntPerXing", ""}

Definition at line 80 of file VertexDecoratorAlg.h.

80 {
81 this, "decor_actualInterPerXing", "actualIntPerXing", ""};

◆ m_mDecor_chi2Over_ndf

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_mDecor_chi2Over_ndf
private
Initial value:
{
this, "decor_chi2Over_ndfKey", "chi2Over_ndf", ""}

Definition at line 66 of file VertexDecoratorAlg.h.

66 {
67 this, "decor_chi2Over_ndfKey", "chi2Over_ndf", ""};

◆ m_mDecor_ntrk

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_mDecor_ntrk
private
Initial value:
{
this, "decor_ntrkKey", "ntrk", ""}

Definition at line 62 of file VertexDecoratorAlg.h.

62 {
63 this, "decor_ntrkKey", "ntrk", ""};

◆ m_mDecor_photon_deltaPhi

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_mDecor_photon_deltaPhi
private
Initial value:
{
this, "decor_photon_deltaPhi", "photon_deltaPhi", ""}

Definition at line 78 of file VertexDecoratorAlg.h.

78 {
79 this, "decor_photon_deltaPhi", "photon_deltaPhi", ""};

◆ m_mDecor_photon_deltaz

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_mDecor_photon_deltaz
private
Initial value:
{
this, "decor_photon_deltaz", "photon_deltaz", ""}

Definition at line 76 of file VertexDecoratorAlg.h.

76 {
77 this, "decor_photon_deltaz", "photon_deltaz", ""};

◆ m_mDecor_sumPt

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_mDecor_sumPt
private
Initial value:
{
this, "decor_sumPtKey", "sumPt", ""}

Definition at line 64 of file VertexDecoratorAlg.h.

64 {
65 this, "decor_sumPtKey", "sumPt", ""};

◆ m_mDecor_weighted_z_asym

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_mDecor_weighted_z_asym
private
Initial value:
{
this, "decor_weighted_z_asym", "weighted_z_asymmetry", ""}

Definition at line 70 of file VertexDecoratorAlg.h.

70 {
71 this, "decor_weighted_z_asym", "weighted_z_asymmetry", ""};

◆ m_mDecor_weighted_z_kurt

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_mDecor_weighted_z_kurt
private
Initial value:
{
this, "decor_weighted_z_kurt", "z_kurtosis", ""}

Definition at line 72 of file VertexDecoratorAlg.h.

72 {
73 this, "decor_weighted_z_kurt", "z_kurtosis", ""};

◆ m_mDecor_z_asym

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_mDecor_z_asym
private
Initial value:
{
this, "decor_z_asymmetryKey", "z_asymmetry", ""}

Definition at line 68 of file VertexDecoratorAlg.h.

68 {
69 this, "decor_z_asymmetryKey", "z_asymmetry", ""};

◆ m_mDecor_z_skew

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_mDecor_z_skew
private
Initial value:
{
this, "decor_z_skew", "z_skewness", ""}

Definition at line 74 of file VertexDecoratorAlg.h.

74 {
75 this, "decor_z_skew", "z_skewness", ""};

◆ m_msg

MsgStream AthHistogramming::m_msg
privateinherited

Cached Message Stream.

Definition at line 250 of file AthHistogramming.h.

◆ m_multiPhotonLinksKey

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_multiPhotonLinksKey
private
Initial value:
{
this, "multiPhotonLinks", "multiPhotonLinks", "Per-vertex link(s) to multi-photon node"}

Definition at line 94 of file VertexDecoratorAlg.h.

94 {
95 this, "multiPhotonLinks", "multiPhotonLinks", "Per-vertex link(s) to multi-photon node"};

◆ m_multiPhotonsOutKey

SG::WriteHandleKey<xAOD::CompositeParticleContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_multiPhotonsOutKey
private
Initial value:
{
this, "multiPhotonsOut", "MultiPhotons", "Output container for multi-photon node"}

Definition at line 92 of file VertexDecoratorAlg.h.

92 {
93 this, "multiPhotonsOut", "MultiPhotons", "Output container for multi-photon node"};

◆ m_muonLinksKey

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_muonLinksKey
private
Initial value:
{
this, "muonLinks", "", "" }

Definition at line 89 of file VertexDecoratorAlg.h.

89 {
90 this, "muonLinks", "", "" };

◆ m_muonsInKey

SG::ReadHandleKey<xAOD::MuonContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_muonsInKey
private
Initial value:
{
this, "muonsIn", "Muons", "containerName to read"}

Definition at line 55 of file VertexDecoratorAlg.h.

55 {
56 this, "muonsIn", "Muons", "containerName to read"};

◆ m_name

std::string AthHistogramming::m_name
privateinherited

Instance name.

Definition at line 247 of file AthHistogramming.h.

◆ m_photonLinksKey

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_photonLinksKey
private
Initial value:
{
this, "photonLinks", "", "" }

Definition at line 83 of file VertexDecoratorAlg.h.

83 {
84 this, "photonLinks", "", "" };

◆ m_photonsInKey

SG::ReadHandleKey<xAOD::PhotonContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_photonsInKey
private
Initial value:
{
this, "photonsIn", "Photons", "containerName to read"}

Definition at line 59 of file VertexDecoratorAlg.h.

59 {
60 this, "photonsIn", "Photons", "containerName to read"};

◆ m_prefix

std::string AthHistogramAlgorithm::m_prefix
privateinherited

Name of the ROOT output stream (file).

Definition at line 88 of file AthHistogramAlgorithm.h.

◆ m_rootDir

std::string AthHistogramAlgorithm::m_rootDir
privateinherited

Name of the ROOT directory.

Definition at line 91 of file AthHistogramAlgorithm.h.

◆ m_streamName

std::string AthHistogramming::m_streamName
privateinherited

Name of the ROOT output stream (file).

Definition at line 228 of file AthHistogramming.h.

◆ m_treeMap

TreeMap_t AthHistogramming::m_treeMap
privateinherited

The map of TTree names to their pointers.

Definition at line 217 of file AthHistogramming.h.

◆ m_trkVtxAssociationTool

ToolHandle<CP::TrackVertexAssociationTool> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_trkVtxAssociationTool
private
Initial value:
{
this, "TrackVertexAssociationTool",
"CP::TrackVertexAssociationTool/TrackVertexAssociationTool",
"Track vertex association tool to use"}

Definition at line 113 of file VertexDecoratorAlg.h.

113 {
114 this, "TrackVertexAssociationTool",
115 "CP::TrackVertexAssociationTool/TrackVertexAssociationTool",
116 "Track vertex association tool to use"};

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Algorithm > >::m_varHandleArraysDeclared
privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vertexInKey

SG::ReadHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_vertexInKey
private
Initial value:
{
this, "vertexIn", "PrimaryVertices", "containerName to read"}

Definition at line 49 of file VertexDecoratorAlg.h.

49 {
50 this, "vertexIn", "PrimaryVertices", "containerName to read"};

◆ m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< Algorithm > >::m_vhka
privateinherited

Definition at line 398 of file AthCommonDataStore.h.

◆ m_zCommonErrorKey

SG::ReadDecorHandleKey<xAOD::PhotonContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_zCommonErrorKey
private
Initial value:
{
this, "zCommonErrorKey", "", "" }

Definition at line 109 of file VertexDecoratorAlg.h.

109 {
110 this, "zCommonErrorKey", "", "" };

◆ m_zCommonKey

SG::ReadDecorHandleKey<xAOD::PhotonContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_zCommonKey
private
Initial value:
{
this, "zCommonKey", "", "" }

Definition at line 107 of file VertexDecoratorAlg.h.

107 {
108 this, "zCommonKey", "", "" };

The documentation for this class was generated from the following files: