ATLAS Offline Software
Loading...
Searching...
No Matches
TruthHitAnalysis Class Reference

#include <TruthHitAnalysis.h>

Inheritance diagram for TruthHitAnalysis:
Collaboration diagram for TruthHitAnalysis:

Public Member Functions

 ~TruthHitAnalysis ()=default
virtual StatusCode execute (const EventContext &ctx) override final
 Execute method.
virtual StatusCode initialize () override final
 AthHistogramAlgorithm (const std::string &name, ISvcLocator *pSvcLocator)
 Constructor with parameters:
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 bool isClonable () const override
 Specify if the algorithm is clonable.
virtual StatusCode sysExecute (const EventContext &ctx) override
 Execute an algorithm.
virtual const DataObjIDColl & extraOutputDeps () const override
 Return the list of extra output dependencies.
virtual bool filterPassed (const EventContext &ctx) const
 Get filter decision:
virtual void setFilterPassed (bool state, const EventContext &ctx) const
 Set filter decision:
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

virtual bool isReEntrant () const override final
 Legacy algorithms are not thread-safe.
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 to fill the internal configuration.
TH1 * bookGetPointer (const TH1 &hist, std::string_view tDir="", std::string_view stream="")
 Simplify the booking and registering (into THistSvc) of histograms.
TH1 * bookGetPointer (TH1 *hist, std::string_view tDir="", std::string_view stream="")
 Simplify the booking and registering (into THistSvc) of histograms.
TH1 * bookGetPointer (TH1 &histRef, std::string_view tDir="", std::string_view stream="")
 Simplify the booking and registering (into THistSvc) of histograms.
TTree * bookGetPointer (const TTree &treeRef, std::string_view tDir="", std::string_view stream="")
 Simplify the booking and registering (into THistSvc) of TTrees.
TGraph * bookGetPointer (const TGraph &graphRef, std::string_view tDir="", std::string_view stream="")
 Simplify the booking and registering (into THistSvc) of TGraphs.
TEfficiency * bookGetPointer (const TEfficiency &eff, std::string_view tDir="", std::string_view stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency.
TEfficiency * bookGetPointer (TEfficiency *eff, std::string_view tDir="", std::string_view stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency.
TEfficiency * bookGetPointer (TEfficiency &effRef, std::string_view tDir="", std::string_view stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency.
StatusCode book (const TH1 &hist, std::string_view tDir="", std::string_view stream="")
 Simplify the booking and registering (into THistSvc) of histograms.
StatusCode book (TH1 *hist, std::string_view tDir="", std::string_view stream="")
 Simplify the booking and registering (into THistSvc) of histograms.
StatusCode book (TH1 &histRef, std::string_view tDir="", std::string_view stream="")
 Simplify the booking and registering (into THistSvc) of histograms.
StatusCode book (const TTree &treeRef, std::string_view tDir="", std::string_view stream="")
 Simplify the booking and registering (into THistSvc) of TTrees.
StatusCode book (const TGraph &graphRef, std::string_view tDir="", std::string_view stream="")
 Simplify the booking and registering (into THistSvc) of TGraphs.
StatusCode book (const TEfficiency &eff, std::string_view tDir="", std::string_view stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency.
StatusCode book (TEfficiency *eff, std::string_view tDir="", std::string_view stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency.
StatusCode book (TEfficiency &effRef, std::string_view tDir="", std::string_view stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency.
TH1 * hist (std::string_view histName, std::string_view tDir="", std::string_view stream="")
 Simplify the retrieval of registered histograms of any type.
TH2 * hist2d (std::string_view histName, std::string_view tDir="", std::string_view stream="")
 Simplify the retrieval of registered 2-d histograms.
TH3 * hist3d (std::string_view histName, std::string_view tDir="", std::string_view stream="")
 Simplify the retrieval of registered 3-d histograms.
TTree * tree (std::string_view treeName, std::string_view tDir="", std::string_view stream="")
 Simplify the retrieval of registered TTrees.
TGraph * graph (std::string_view graphName, std::string_view tDir="", std::string_view stream="")
 Simplify the retrieval of registered TGraphs.
TEfficiency * efficiency (std::string_view effName, std::string_view tDir="", std::string_view stream="")
 Simplify the retrieval of registered TEfficiency.

Private Types

typedef ServiceHandle< StoreGateSvc > StoreGateSvc_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>
std::string buildBookingString (std::string_view histName, std::string_view tDir, std::string_view stream, bool usePrefixPostfix=false) const
 Method to build individual booking string.
hash_t hash (std::string_view histName) const
 Method to calculate a 32-bit hash from a string.

Private Attributes

TH1 * m_h_n_vert {nullptr}
 Some variables.
TH1 * m_h_n_part {nullptr}
TH1 * m_h_n_vert_prim {nullptr}
TH1 * m_h_n_part_prim {nullptr}
TH1 * m_h_n_vert_sec {nullptr}
TH1 * m_h_n_part_sec {nullptr}
TH1 * m_h_vtx_x {nullptr}
TH1 * m_h_vtx_y {nullptr}
TH1 * m_h_vtx_z {nullptr}
TH1 * m_h_vtx_r {nullptr}
TH1 * m_h_vtx_proc_x {nullptr}
TH1 * m_h_vtx_proc_y {nullptr}
TH1 * m_h_vtx_proc_z {nullptr}
TH1 * m_h_vtx_proc_r {nullptr}
TH1 * m_h_vtx_prim_x {nullptr}
TH1 * m_h_vtx_prim_y {nullptr}
TH1 * m_h_vtx_prim_z {nullptr}
TH1 * m_h_vtx_prim_r {nullptr}
TH2 * m_h_vtx_prim_xy {nullptr}
TH2 * m_h_vtx_prim_zr {nullptr}
TH2 * m_h_vtx_sec_xy {nullptr}
TH2 * m_h_vtx_sec_zr {nullptr}
TH1 * m_h_n_generations {nullptr}
TH1 * m_h_truth_px {nullptr}
TH1 * m_h_truth_py {nullptr}
TH1 * m_h_truth_pz {nullptr}
TH1 * m_h_truth_pt {nullptr}
TH1 * m_h_truth_eta {nullptr}
TH1 * m_h_truth_phi {nullptr}
TH1 * m_h_barcode {nullptr}
TH1 * m_h_part_status {nullptr}
TH1 * m_h_part_pdgid {nullptr}
TH1 * m_h_part_pdgid_sec {nullptr}
TH1 * m_h_part_eta {nullptr}
TH1 * m_h_part_phi {nullptr}
TH1 * m_h_part_p {nullptr}
std::vector< float > * m_vtx_x {nullptr}
std::vector< float > * m_vtx_y {nullptr}
std::vector< float > * m_vtx_z {nullptr}
std::vector< float > * m_vtx_r {nullptr}
std::vector< float > * m_vtx_barcode {nullptr}
std::vector< float > * m_vtx_proc_x {nullptr}
std::vector< float > * m_vtx_proc_y {nullptr}
std::vector< float > * m_vtx_proc_z {nullptr}
std::vector< float > * m_vtx_proc_r {nullptr}
std::vector< float > * m_vtx_proc_barcode {nullptr}
std::vector< float > * m_vtx_prim_x {nullptr}
std::vector< float > * m_vtx_prim_y {nullptr}
std::vector< float > * m_vtx_prim_z {nullptr}
std::vector< float > * m_vtx_prim_r {nullptr}
std::vector< float > * m_vtx_prim_barcode {nullptr}
std::vector< float > * m_truth_px {nullptr}
std::vector< float > * m_truth_py {nullptr}
std::vector< float > * m_truth_pz {nullptr}
std::vector< float > * m_truth_pt {nullptr}
std::vector< float > * m_truth_eta {nullptr}
std::vector< float > * m_truth_phi {nullptr}
std::vector< float > * m_barcode {nullptr}
std::vector< float > * m_status {nullptr}
std::vector< float > * m_pdgid {nullptr}
TTree * m_tree {nullptr}
Gaudi::Property< std::string > m_path {this, "HistPath","/TruthHitAnalysis/"}
Gaudi::Property< std::string > m_ntupleFileName {this, "NtupleFileName","/TruthHitAnalysis/"}
SG::ReadHandleKey< McEventCollection > m_readKey {this, "InputKey", "TruthEvent"}
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
 Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
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 14 of file TruthHitAnalysis.h.

Member Typedef Documentation

◆ EffMap_t

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

Typedef for convenience.

Definition at line 199 of file AthHistogramming.h.

◆ GraphMap_t

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

Typedef for convenience.

Definition at line 213 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 192 of file AthHistogramming.h.

◆ StoreGateSvc_t

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

Definition at line 376 of file AthCommonDataStore.h.

◆ TreeMap_t

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

Typedef for convenience.

Definition at line 206 of file AthHistogramming.h.

Constructor & Destructor Documentation

◆ ~TruthHitAnalysis()

TruthHitAnalysis::~TruthHitAnalysis ( )
default

Member Function Documentation

◆ AthHistogramAlgorithm()

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

Constructor with parameters:

Definition at line 40 of file AthHistogramAlgorithm.cxx.

32 :
33 ::AthAlgorithm( name, pSvcLocator ),
34 ::AthHistogramming( name ),
35 m_histSvc ( "THistSvc/THistSvc", name )
36{
37 //
38 // Property declaration
39 //
40 // The histogram service
41 declareProperty("THistSvc",
42 m_histSvc = ServiceHandle<ITHistSvc>("THistSvc/THistSvc", name),
43 "Handle to a THistSvc instance: it will be used to write "
44 "ROOT objects to ROOT files" );
45 // declareProperty("THistService", m_histSvc, "The THistSvc" );
46
47 declareProperty("RootStreamName", m_prefix = "/ANALYSIS", "Name of the output ROOT stream (file) that the THistSvc uses");
48 declareProperty("RootDirName", m_rootDir = "",
49 "Name of the ROOT directory inside the ROOT file where the histograms will go");
50
51 declareProperty( "HistNamePrefix", m_histNamePrefix = "", "The prefix for the histogram THx name" );
52 declareProperty( "HistNamePostfix", m_histNamePostfix = "", "The postfix for the histogram THx name" );
53
54 declareProperty( "HistTitlePrefix", m_histTitlePrefix = "", "The prefix for the histogram THx title" );
55 declareProperty( "HistTitlePostfix", m_histTitlePostfix = "", "The postfix for the histogram THx title" );
56}
AthAlgorithm(const std::string &name, ISvcLocator *pSvcLocator)
Constructor.
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)
std::string m_histNamePrefix
The prefix for the histogram THx name.
ServiceHandle< ITHistSvc > m_histSvc
Default constructor: AthHistogramAlgorithm();.
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.
AthHistogramming(const std::string &name)
Constructor with parameters:

◆ book() [1/8]

StatusCode AthHistogramming::book ( const TEfficiency & eff,
std::string_view tDir = "",
std::string_view stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TEfficiency.

Definition at line 329 of file AthHistogramming.h.

330{
331 // We need to create a non-const clone
332 TEfficiency* effClone = dynamic_cast< TEfficiency* >( eff.Clone() );
333 if ( !effClone ) {
334 m_msg << MSG::ERROR << "Couldn't create a TEfficiency clone" << endmsg;
335 return StatusCode::FAILURE;
336 }
337 return this->book( *effClone, tDir, stream );
338}
#define endmsg
StatusCode book(const TH1 &hist, std::string_view tDir="", std::string_view 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,
std::string_view tDir = "",
std::string_view stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TGraphs.

◆ book() [3/8]

StatusCode AthHistogramming::book ( const TH1 & hist,
std::string_view tDir = "",
std::string_view stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of histograms.

Definition at line 297 of file AthHistogramming.h.

298{
299 // We need to create a non-const clone
300 TH1* histClone = dynamic_cast< TH1* >( hist.Clone() );
301 if ( !histClone ) {
302 m_msg << MSG::ERROR << "Couldn't create a TH1 clone" << endmsg;
303 return StatusCode::FAILURE;
304 }
305 return this->book( *histClone, tDir, stream );
306}
TH1 * hist(std::string_view histName, std::string_view tDir="", std::string_view stream="")
Simplify the retrieval of registered histograms of any type.

◆ book() [4/8]

StatusCode AthHistogramming::book ( const TTree & treeRef,
std::string_view tDir = "",
std::string_view stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TTrees.

Definition at line 397 of file AthHistogramming.h.

398{
399 // Call the other Book method and see if it returns a valid pointer
400 TTree* treePointer = this->bookGetPointer( treeRef, tDir, stream );
401 if ( treePointer )
402 {
403 return StatusCode::SUCCESS;
404 }
405 else
406 {
407 return StatusCode::FAILURE;
408 }
409}
TH1 * bookGetPointer(const TH1 &hist, std::string_view tDir="", std::string_view stream="")
Simplify the booking and registering (into THistSvc) of histograms.

◆ book() [5/8]

StatusCode AthHistogramming::book ( TEfficiency & effRef,
std::string_view tDir = "",
std::string_view stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TEfficiency.

Definition at line 350 of file AthHistogramming.h.

351{
352 // Call the other Book method and see if it returns a valid pointer
353 TEfficiency* effPointer = this->bookGetPointer( effRef, tDir, stream );
354 if ( !effPointer ) {
355 m_msg << MSG::ERROR << "Couldn't book a TEfficiency" << endmsg;
356 return StatusCode::FAILURE;
357 }
358 return StatusCode::SUCCESS;
359}

◆ book() [6/8]

StatusCode AthHistogramming::book ( TEfficiency * eff,
std::string_view tDir = "",
std::string_view stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TEfficiency.

Definition at line 340 of file AthHistogramming.h.

341{
342 if ( !eff ) {
343 m_msg << MSG::ERROR << "Got a zero pointer to a TEfficiency" << endmsg;
344 return StatusCode::FAILURE;
345 }
346 return this->book( *eff, tDir, stream );
347}

◆ book() [7/8]

StatusCode AthHistogramming::book ( TH1 & histRef,
std::string_view tDir = "",
std::string_view stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of histograms.

Definition at line 318 of file AthHistogramming.h.

319{
320 // Call the other Book method and see if it returns a valid pointer
321 TH1* histPointer = this->bookGetPointer( histRef, tDir, stream );
322 if ( !histPointer ) {
323 m_msg << MSG::ERROR << "Couldn't book a TH1" << endmsg;
324 return StatusCode::FAILURE;
325 }
326 return StatusCode::SUCCESS;
327}

◆ book() [8/8]

StatusCode AthHistogramming::book ( TH1 * hist,
std::string_view tDir = "",
std::string_view stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of histograms.

Definition at line 308 of file AthHistogramming.h.

309{
310 if ( !hist ) {
311 m_msg << MSG::ERROR << "Got a zero pointer to a TH1" << endmsg;
312 return StatusCode::FAILURE;
313 }
314 return this->book( *hist, tDir, stream );
315}

◆ bookGetPointer() [1/8]

TEfficiency * AthHistogramming::bookGetPointer ( const TEfficiency & eff,
std::string_view tDir = "",
std::string_view stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TEfficiency.

Definition at line 275 of file AthHistogramming.h.

276{
277 // We need to create a non-const clone
278 TEfficiency* histClone = dynamic_cast< TEfficiency* >( hist.Clone() );
279 if ( !histClone ) {
280 m_msg << MSG::ERROR << "Couldn't create a TEfficiency clone in bookGetPointer" << endmsg;
281 return 0;
282 }
283 return this->bookGetPointer( *histClone, tDir, stream );
284
285}

◆ bookGetPointer() [2/8]

TGraph * AthHistogramming::bookGetPointer ( const TGraph & graphRef,
std::string_view tDir = "",
std::string_view stream = "" )
protectedinherited

Simplify the booking and registering (into THistSvc) of TGraphs.

Definition at line 374 of file AthHistogramming.cxx.

375{
376 // Get a pointer
377 const TGraph* graphPointer = &graphRef;
378
379 // Check that we got a valid pointer
380 if ( !graphPointer )
381 {
382 m_msg << MSG::WARNING
383 << "We got an invalid TGraph pointer in the BookGetPointer(TGraph*) method of the class" << m_name
384 << "!" << endmsg;
385 return NULL;
386 }
387
388 // Modify the name and title according to the prefixes of this classes instance
389 std::string graphName = graphPointer->GetName();
390 const std::string graphTitle = graphPointer->GetTitle();
391
392 // Check if the hash for this graphName already exists, i.e., if we have a hash collision
393 const hash_t graphHash = this->hash(graphName);
394 GraphMap_t::const_iterator it = m_graphMap.find( graphHash );
395 if ( it != m_graphMap.end() ) // It does exist!
396 {
397 m_msg << MSG::WARNING
398 << "Detected a hash collision. The hash for the TGraph with name=" << graphName
399 << " already exists and points to a TGraph with name=" << it->second->GetName()
400 << " NOT going to book the new histogram and returning a NULL pointer!" << endmsg;
401 return NULL;
402 }
403
404 // Create a clone that has the new name
405 TGraph* graphClone = dynamic_cast< TGraph* >( graphPointer->Clone((m_histNamePrefix+graphName+m_histNamePostfix).c_str()) );
406 if( !graphClone )
407 {
408 m_msg << MSG::WARNING
409 << "We couldn't clone the TGraph in the BookGetPointer(TGraph&) method of the class" << m_name
410 << "!" << endmsg;
411 return NULL;
412 }
413 graphClone->SetTitle ((m_histTitlePrefix+graphTitle+m_histTitlePostfix).c_str());
414
415 // Massage the final string to book things
416 std::string bookingString = buildBookingString( graphName, tDir, stream );
417
418 // Register the TGraph into the THistSvc
419 if ( !((histSvc()->regGraph(bookingString, graphClone)).isSuccess()) )
420 {
421 m_msg << MSG::WARNING
422 << "Problem registering TGraph with name " << graphName
423 << ", title " << graphTitle
424 << " in " << m_name << "!" << endmsg;
425 return NULL;
426 }
427
428 // Also register it in the local map of string to pointer
429 m_graphMap.insert( m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphClone ) );
430
431 return graphClone;
432}
const ServiceHandle< ITHistSvc > & histSvc() const
The standard THistSvc (for writing histograms and TTrees and more to a root file) Returns (kind of) a...
hash_t hash(std::string_view histName) const
Method to calculate a 32-bit hash from a string.
std::string buildBookingString(std::string_view histName, std::string_view tDir, std::string_view stream, bool usePrefixPostfix=false) const
Method to build individual booking string.
std::string m_histNamePostfix
The postfix for the histogram THx name.
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.
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,
std::string_view tDir = "",
std::string_view stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of histograms.

Definition at line 254 of file AthHistogramming.h.

255{
256 // We need to create a non-const clone
257 TH1* histClone = dynamic_cast< TH1* >( hist.Clone() );
258 if ( !histClone ) {
259 m_msg << MSG::ERROR << "Couldn't create a TH1 clone in bookGetPointer" << endmsg;
260 return 0;
261 }
262 return this->bookGetPointer( *histClone, tDir, stream );
263
264}

◆ bookGetPointer() [4/8]

TTree * AthHistogramming::bookGetPointer ( const TTree & treeRef,
std::string_view tDir = "",
std::string_view stream = "" )
protectedinherited

Simplify the booking and registering (into THistSvc) of TTrees.

Definition at line 268 of file AthHistogramming.cxx.

269{
270 // Get a pointer
271 const TTree* treePointer = &treeRef;
272
273 // Check that we got a valid pointer
274 if ( !treePointer )
275 {
276 m_msg << MSG::WARNING
277 << "We got an invalid TTree pointer in the BookGetPointer(TTree*) method of the class" << m_name
278 << "!" << endmsg;
279 return NULL;
280 }
281
282 // Modify the name and title according to the prefixes of this classes instance
283 std::string treeName = treePointer->GetName();
284 const std::string treeTitle = treePointer->GetTitle();
285
286 // Check if the hash for this treeName already exists, i.e., if we have a hash collision
287 const hash_t treeHash = this->hash(treeName);
288 TreeMap_t::const_iterator it = m_treeMap.find( treeHash );
289 if ( it != m_treeMap.end() ) // It does exist!
290 {
291 m_msg << MSG::WARNING
292 << "Detected a hash collision. The hash for the TTree with name=" << treeName
293 << " already exists and points to a TTree with name=" << it->second->GetName()
294 << " NOT going to book the new histogram and returning a NULL pointer!" << endmsg;
295 return NULL;
296 }
297
298 // Create a clone that has the new name
299 TTree* treeClone = dynamic_cast< TTree* >( treePointer->Clone(treeName.c_str()) );
300 if( !treeClone )
301 {
302 m_msg << MSG::WARNING
303 << "We couldn't clone the TTree in the BookGetPointer(TTree&) method of the class" << m_name
304 << "!" << endmsg;
305 return NULL;
306 }
307 treeClone->SetTitle (treeTitle.c_str());
308
309 // Massage the final string to book things
310 std::string bookingString = buildBookingString( treeName, tDir, stream );
311
312 // Register the TTree into the THistSvc
313 if ( !((histSvc()->regTree(bookingString, treeClone)).isSuccess()) )
314 {
315 m_msg << MSG::WARNING
316 << "Problem registering TTree with name " << treeName
317 << ", title " << treeTitle
318 << " in " << m_name << "!" << endmsg;
319 return NULL;
320 }
321
322 // Also register it in the local map of string to pointer
323 m_treeMap.insert( m_treeMap.end(), std::pair< const hash_t, TTree* >( treeHash, treeClone ) );
324
325 return treeClone;
326}
TreeMap_t m_treeMap
The map of TTree names to their pointers.

◆ bookGetPointer() [5/8]

TEfficiency * AthHistogramming::bookGetPointer ( TEfficiency & effRef,
std::string_view tDir = "",
std::string_view stream = "" )
protectedinherited

Simplify the booking and registering (into THistSvc) of TEfficiency.

Definition at line 139 of file AthHistogramming.cxx.

140{
141 std::string originalName{effRef.GetName()};
142 const auto effName = splitName(originalName).second;
143 const std::string effTitle(effRef.GetTitle());
144 std::string bookingString = buildBookingString(originalName, tDir, stream);
145 std::string finalEffName{m_histNamePrefix};
146 finalEffName += effName;
147 finalEffName += m_histNamePostfix;
148 effRef.SetTitle((m_histTitlePrefix + effTitle + m_histTitlePostfix).c_str());
149 effRef.SetName(finalEffName.c_str());
150
151 const hash_t effHash = this->hash(effName);
152 EffMap_t::const_iterator it = m_effMap.find(effHash);
153 if (it != m_effMap.end()) {
154 m_msg << MSG::WARNING
155 << "Detected a hash collision. The hash for the TEfficiency with name=" << effName
156 << " already exists and points to a TEfficiency with name=" << it->second->GetName()
157 << " NOT going to book the new TEfficiency and returning a NULL pointer!" << endmsg;
158 return nullptr;
159 }
160
161 if (!histSvc()->regEfficiency(bookingString, &effRef).isSuccess()) {
162 m_msg << MSG::WARNING
163 << "Problem registering TEfficiency with name " << effName
164 << ", name prefix " << m_histNamePrefix
165 << ", title " << effTitle
166 << ", title prefix " << m_histTitlePrefix
167 << ", and title postfix " << m_histTitlePostfix
168 << " in " << m_name << "!" << endmsg;
169 return nullptr;
170 }
171
172 m_effMap.insert(m_effMap.end(), std::pair<const hash_t, TEfficiency*>(effHash, &effRef));
173
174 return &effRef;
175}
EffMap_t m_effMap
The map of histogram names to their pointers.

◆ bookGetPointer() [6/8]

TEfficiency * AthHistogramming::bookGetPointer ( TEfficiency * eff,
std::string_view tDir = "",
std::string_view stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TEfficiency.

Definition at line 287 of file AthHistogramming.h.

288{
289 if ( !hist ) {
290 m_msg << MSG::ERROR << "Got a zero pointer to a TEfficiency in bookGetPointer" << endmsg;
291 return 0;
292 }
293 return this->bookGetPointer( *hist, tDir, stream );
294}

◆ bookGetPointer() [7/8]

TH1 * AthHistogramming::bookGetPointer ( TH1 & histRef,
std::string_view tDir = "",
std::string_view stream = "" )
protectedinherited

Simplify the booking and registering (into THistSvc) of histograms.

Definition at line 97 of file AthHistogramming.cxx.

97 {
98
99 const std::string originalName{histRef.GetName()};
100 const auto histName = splitName(originalName).second;
101
102 const std::string histTitle(histRef.GetTitle());
103 std::string bookingString = buildBookingString(originalName, tDir, stream);
104
105 std::string finalHistName{m_histNamePrefix};
106 finalHistName += histName;
107 finalHistName += m_histNamePostfix;
108
109 histRef.SetTitle((m_histTitlePrefix + histTitle + m_histTitlePostfix).c_str());
110 histRef.SetName(finalHistName.c_str());
111
112 const hash_t histHash = this->hash(histName);
113 HistMap_t::const_iterator it = m_histMap.find(histHash);
114 if (it != m_histMap.end()) {
115 m_msg << MSG::WARNING
116 << "Detected a hash collision. The hash for the histogram with name=" << histName
117 << " already exists and points to a histogram with name=" << it->second->GetName()
118 << " NOT going to book the new histogram and returning a NULL pointer!" << endmsg;
119 return nullptr;
120 }
121
122 if (!histSvc()->regHist(bookingString, &histRef).isSuccess()) {
123 m_msg << MSG::WARNING
124 << "Problem registering histogram with name " << histName
125 << ", name prefix " << m_histNamePrefix
126 << ", title " << histTitle
127 << ", title prefix " << m_histTitlePrefix
128 << ", and title postfix " << m_histTitlePostfix
129 << " in " << m_name << "!" << endmsg;
130 return nullptr;
131 }
132
133 m_histMap.insert(m_histMap.end(), std::pair<const hash_t, TH1*>(histHash, &histRef));
134
135 return &histRef;
136}
HistMap_t m_histMap
The map of histogram names to their pointers.

◆ bookGetPointer() [8/8]

TH1 * AthHistogramming::bookGetPointer ( TH1 * hist,
std::string_view tDir = "",
std::string_view stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of histograms.

Definition at line 266 of file AthHistogramming.h.

267{
268 if ( !hist ) {
269 m_msg << MSG::ERROR << "Got a zero pointer to a TH1 in bookGetPointer" << endmsg;
270 return 0;
271 }
272 return this->bookGetPointer( *hist, tDir, stream );
273}

◆ buildBookingString()

std::string AthHistogramming::buildBookingString ( std::string_view histName,
std::string_view tDir,
std::string_view stream,
bool usePrefixPostfix = false ) const
privateinherited

Method to build individual booking string.

Definition at line 500 of file AthHistogramming.cxx.

501 {
502 if (tDir.empty()) {
503 tDir = m_rootDir;
504 }
505 if (stream.empty()) {
507 }
508 const auto [histDir, baseName] = splitName(histName);
509 std::string bookingString{"/"};
510 bookingString += stream;
511 bookingString += '/';
512 bookingString += tDir;
513 if (!histDir.empty()) {
514 bookingString += '/';
515 bookingString += histDir;
516 }
517 bookingString += '/';
518 if (usePrefixPostfix) {
519 bookingString += m_histNamePrefix;
520 }
521 bookingString += baseName;
522 if (usePrefixPostfix) {
523 bookingString += m_histNamePostfix;
524 }
525 const auto tail = std::ranges::unique(
526 bookingString, [](char lhs, char rhs) { return lhs == '/' && rhs == '/'; });
527 bookingString.erase(tail.begin(), tail.end());
528 return bookingString;
529}
std::string m_rootDir
Name of the ROOT directory.
std::string m_streamName
Name of the ROOT output stream (file).
std::string tail(std::string s, const std::string &pattern)
tail of a 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 to fill the internal configuration.

Definition at line 64 of file AthHistogramming.cxx.

68{
71 m_rootDir = rootDir;
72 m_histNamePrefix = histNamePrefix;
73 m_histNamePostfix = histNamePostfix;
74 m_histTitlePrefix = histTitlePrefix;
75 m_histTitlePostfix = histTitlePostfix;
76
77 return StatusCode::SUCCESS;
78}
ServiceHandle< ITHistSvc > m_histSvc
Pointer to the THistSvc (event store by default).

◆ declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Gaudi::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 }

◆ declareProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Gaudi::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< Gaudi::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 ( std::string_view effName,
std::string_view tDir = "",
std::string_view stream = "" )
protectedinherited

Simplify the retrieval of registered TEfficiency.

Definition at line 214 of file AthHistogramming.cxx.

215{
216 // Build a 32 bit hash out of the name
217 const hash_t effHash = this->hash(effName);
218
219 // See if this entry exists in the map
220 EffMap_t::const_iterator it = m_effMap.find( effHash );
221 if ( it == m_effMap.end() ) // It doesn't exist!
222 { // Let's see into the THistSvc if somebody else has registered the TEfficiency...
223 // Massage the final string to book things
224 std::string bookingString = buildBookingString(effName, tDir, stream ,false);
225
226 TEfficiency* effPointer{};
227 if ( !((histSvc()->getEfficiency(bookingString, effPointer)).isSuccess()) )
228 {
229 // Book things
230 std::string bookingString = buildBookingString( effName, tDir, stream, true );
231
232 if ( !((histSvc()->getEfficiency(bookingString, effPointer)).isSuccess()) )
233 {
234 m_msg << MSG::WARNING
235 << "Problem retrieving the TEfficiency with name (including pre- and post-fixes) "
236 << bookingString
237 << " or with name " << effName
238 << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
239 << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
240 return NULL;
241 }
242 // If we get to here, we actually found the TEfficiency in the THistSvc.
243 // So let's add it to the local cache map and return its pointer
244 m_effMap.insert( m_effMap.end(), std::pair< const hash_t, TEfficiency* >( effHash, effPointer ) );
245 return effPointer;
246 }
247 // If we get to here, we actually found the TEfficiency in the THistSvc.
248 // So let's add it to the local cache map and return its pointer
249 m_effMap.insert( m_effMap.end(), std::pair< const hash_t, TEfficiency* >( effHash, effPointer ) );
250 return effPointer;
251 }
252
253 // Return the pointer to the TEfficiency that we got from the local cache map
254 return it->second;
255}
std::pair< StatusCode, TEfficiency * > getEfficiency(ITHistSvc &svc, const std::string &name)

◆ evtStore()

ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< Gaudi::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 TruthHitAnalysis::execute ( const EventContext & ctx)
finaloverridevirtual

Execute method.

Implements AthAlgorithm.

Definition at line 198 of file TruthHitAnalysis.cxx.

198 {
199 ATH_MSG_DEBUG( "In TruthHitAnalysis::execute()" );
200
201 m_vtx_x->clear();
202 m_vtx_y->clear();
203 m_vtx_z->clear();
204 m_vtx_r->clear();
205 m_vtx_barcode->clear();
206 m_vtx_proc_x->clear();
207 m_vtx_proc_y->clear();
208 m_vtx_proc_z->clear();
209 m_vtx_proc_r->clear();
210 m_vtx_proc_barcode->clear();
211 m_vtx_prim_x->clear();
212 m_vtx_prim_y->clear();
213 m_vtx_prim_z->clear();
214 m_vtx_prim_r->clear();
215 m_vtx_prim_barcode->clear();
216 m_truth_px->clear();
217 m_truth_py->clear();
218 m_truth_pz->clear();
219 m_truth_pt->clear();
220 m_truth_eta->clear();
221 m_truth_phi->clear();
222 m_barcode->clear();
223 m_status->clear();
224 m_pdgid->clear();
225
226 const McEventCollection* mcCollection{nullptr};
227 ATH_CHECK(SG::get(mcCollection, m_readKey, ctx));
228
229
230 McEventCollection::const_iterator currentGenEventIter = mcCollection->begin();
231 if (currentGenEventIter != mcCollection->end()) {
232 int nvtx = 0;
233 int nvtx_sec=0;
234
235 const auto &procVtx = HepMC::signal_process_vertex(*currentGenEventIter);
236 const auto &barcodes = (*currentGenEventIter)->attribute<HepMC::GenEventBarcodes> (HepMCStr::barcodes);
237 std::map<int,int> id_to_barcode_map;
238 if (barcodes) id_to_barcode_map = barcodes->id_to_barcode_map();
239 for (const auto& vtx: (*currentGenEventIter)->vertices()) {
240 int bcode = id_to_barcode_map[vtx->id()];
241 double x = vtx->position().x();
242 double y = vtx->position().y();
243 double z = vtx->position().z();
244 double r = std::sqrt(x*x+y*y);
245 m_h_vtx_x->Fill(x);
246 m_h_vtx_y->Fill(y);
247 m_h_vtx_r->Fill(r);
248 m_h_vtx_z->Fill(z);
249
250 m_vtx_x->push_back(x);
251 m_vtx_y->push_back(y);
252 m_vtx_r->push_back(r);
253 m_vtx_z->push_back(z);
254 m_vtx_barcode->push_back(bcode);
255
256 if (vtx == procVtx) {
257 m_h_vtx_proc_x->Fill(x);
258 m_h_vtx_proc_y->Fill(y);
259 m_h_vtx_proc_r->Fill(r);
260 m_h_vtx_proc_z->Fill(z);
261
262 m_vtx_proc_x->push_back(x);
263 m_vtx_proc_y->push_back(y);
264 m_vtx_proc_r->push_back(r);
265 m_vtx_proc_z->push_back(z);
266 m_vtx_proc_barcode->push_back(bcode);
267 }
268
269 if (!HepMC::is_simulation_vertex(vtx)) {
270 m_h_vtx_prim_x->Fill(x);
271 m_h_vtx_prim_y->Fill(y);
272 m_h_vtx_prim_r->Fill(r);
273 m_h_vtx_prim_z->Fill(z);
274 m_h_vtx_prim_xy->Fill(x,y);
275 m_h_vtx_prim_zr->Fill(z,r);
276
277 m_vtx_prim_x->push_back(x);
278 m_vtx_prim_y->push_back(y);
279 m_vtx_prim_r->push_back(r);
280 m_vtx_prim_z->push_back(z);
281 m_vtx_prim_barcode->push_back(bcode);
282 ++nvtx;
283 }
284 else {
285 m_h_vtx_sec_xy->Fill(x,y);
286 m_h_vtx_sec_zr->Fill(z,r);
287 ++nvtx_sec;
288 }
289 } //End iteration over vertices
290
291 m_h_n_vert->Fill(nvtx+nvtx_sec);
292 m_h_n_vert_prim->Fill(nvtx);
293 m_h_n_vert_sec->Fill(nvtx_sec);
294
295 int npart_prim=0;
296 int npart_sec=0;
297
298 for (auto currentGenParticle: *(*currentGenEventIter)) {
299 const HepMC::FourVector mom = currentGenParticle->momentum();
300 int currentGenParticlebarcode = id_to_barcode_map[currentGenParticle->id()];
301 m_h_truth_px->Fill(mom.x());
302 m_h_truth_py->Fill(mom.y());
303 m_h_truth_pz->Fill(mom.z());
304 m_h_truth_pt->Fill(mom.perp());
305 m_h_truth_eta->Fill(mom.eta());
306 m_h_truth_phi->Fill(mom.phi());
307 m_h_barcode->Fill(currentGenParticlebarcode);
308 m_h_part_status->Fill(currentGenParticle->status());
309 m_truth_px->push_back(mom.x());
310 m_truth_py->push_back(mom.y());
311 m_truth_pz->push_back(mom.z());
312 m_truth_pt->push_back(mom.perp());
313 m_truth_eta->push_back(mom.eta());
314 m_truth_phi->push_back(mom.phi());
315 m_barcode->push_back(currentGenParticlebarcode);
316 m_status->push_back(currentGenParticle->status());
317
318 int pdg = currentGenParticle->pdg_id();
319 m_pdgid->push_back(pdg);
320
321 if (!HepMC::is_simulation_particle(currentGenParticle)) {
322 m_h_part_pdgid->Fill(pdg);
323 m_h_part_p->Fill(std::sqrt(mom.x()*mom.x()+mom.y()*mom.y()+mom.z()*mom.z()));
324 m_h_part_eta->Fill(mom.eta());
325 m_h_part_phi->Fill(mom.phi());
326 ++npart_prim;
327 m_h_n_generations->Fill(1);
328/* The conditions don't make any sense.
329 if (currentGenParticlebarcode < HepMC::PHOTOSMIN) {
330 m_h_n_generations->Fill(0);
331 }
332 else {
333 m_h_n_generations->Fill(1);
334 }
335*/
336 } //End simulation particle
337 else {
338 m_h_part_pdgid_sec->Fill(pdg);
339 ++npart_sec;
340 const int gen = HepMC::generations(currentGenParticle) + 2;
341 m_h_n_generations->Fill(gen);
342 }
343 } // End iteration over particles
344
345 m_h_n_part_prim->Fill(npart_prim);
346 m_h_n_part_sec->Fill(npart_sec);
347 m_h_n_part->Fill(npart_prim+npart_sec);
348 } // End mcCollection
349
350 m_tree->Fill();
351
352 return StatusCode::SUCCESS;
353}
#define ATH_CHECK
Evaluate an expression and check for errors.
#define ATH_MSG_DEBUG(x,...)
#define y
#define x
#define z
DataModel_detail::const_iterator< DataVector > const_iterator
Definition DataVector.h:861
const_iterator end() const noexcept
Return a const_iterator pointing past the end of the collection.
const_iterator begin() const noexcept
Return a const_iterator pointing at the beginning of the collection.
std::vector< float > * m_vtx_z
std::vector< float > * m_truth_pt
std::vector< float > * m_truth_py
std::vector< float > * m_vtx_prim_barcode
std::vector< float > * m_vtx_proc_r
std::vector< float > * m_pdgid
std::vector< float > * m_vtx_y
std::vector< float > * m_vtx_prim_z
TH1 * m_h_n_vert
Some variables.
std::vector< float > * m_vtx_proc_y
std::vector< float > * m_vtx_prim_x
std::vector< float > * m_vtx_x
SG::ReadHandleKey< McEventCollection > m_readKey
std::vector< float > * m_vtx_proc_x
std::vector< float > * m_vtx_barcode
std::vector< float > * m_vtx_proc_z
std::vector< float > * m_vtx_prim_r
std::vector< float > * m_vtx_prim_y
std::vector< float > * m_truth_pz
std::vector< float > * m_truth_eta
std::vector< float > * m_barcode
std::vector< float > * m_truth_phi
std::vector< float > * m_status
std::vector< float > * m_vtx_r
std::vector< float > * m_truth_px
std::vector< float > * m_vtx_proc_barcode
int r
Definition globals.cxx:22
bool is_simulation_vertex(const T &v)
Method to establish if the vertex was created during simulation (TODO migrate to be based on status).
HepMC3::FourVector FourVector
ConstGenVertexPtr signal_process_vertex(const GenEvent *e)
Definition GenEvent.h:597
bool is_simulation_particle(const T &p)
Method to establish if a particle (or barcode) was created during the simulation (TODO update to be s...
int generations(const T &p)
Method to return how many interactions a particle has undergone during simulation (TODO migrate to be...
const T * get(const ReadCondHandleKey< T > &key, const EventContext &ctx)
Convenience function to retrieve an object given a ReadCondHandleKey.
barcodes(beg, end, sz)
Definition Dumpers.py:2322

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Gaudi::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 & AthCommonAlgorithm< Gaudi::Algorithm >::extraOutputDeps ( ) const
overridevirtualinherited

Return the list of extra output dependencies.

This list is extended to include symlinks implied by inheritance relations.

Definition at line 89 of file AthCommonAlgorithm.cxx.

54{
55 // If we didn't find any symlinks to add, just return the collection
56 // from the base class. Otherwise, return the extended collection.
57 if (!m_extendedExtraObjects.empty()) {
59 }
61}
Common base class for algorithms.

◆ filterPassed()

virtual bool AthCommonAlgorithm< Gaudi::Algorithm >::filterPassed ( const EventContext & ctx) const
inlinevirtualinherited

Get filter decision:

Definition at line 93 of file AthCommonAlgorithm.h.

93 {
94 return execState( ctx ).filterPassed();
95 }
virtual bool filterPassed(const EventContext &ctx) const
Get filter decision:

◆ graph()

TGraph * AthHistogramming::graph ( std::string_view graphName,
std::string_view tDir = "",
std::string_view stream = "" )
protectedinherited

Simplify the retrieval of registered TGraphs.

Definition at line 438 of file AthHistogramming.cxx.

439{
440 // Build a 32 bit hash out of the name
441 const hash_t graphHash = this->hash(graphName);
442
443 // See if this entry exists in the map
444 GraphMap_t::const_iterator it = m_graphMap.find( graphHash );
445 if ( it == m_graphMap.end() ) // It doesn't exist!
446 { // Let's see into the THistSvc if somebody else has registered the TGraph...
447
448 // Massage the final string to book things
449 std::string bookingString = buildBookingString( graphName, tDir, stream, false);
450
451 TGraph* graphPointer(nullptr);
452 if ( !((histSvc()->getGraph(bookingString, graphPointer)).isSuccess()) )
453 {
454 // Massage the final string to book things
455 std::string bookingString = buildBookingString( graphName, tDir, stream, true );
456
457 if ( !((histSvc()->getGraph(bookingString, graphPointer)).isSuccess()) )
458 {
459 m_msg << MSG::WARNING
460 << "Problem retrieving the TGraph with name (including pre- and post-fixes) "
461 << bookingString
462 << " or with name " << graphName
463 << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
464 << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
465 return nullptr;
466 }
467 // If we get to here, we actually found the TGraph in the THistSvc.
468 // So let's add it to the local cache map and return its pointer
469 m_graphMap.insert( m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphPointer ) );
470 return graphPointer;
471 }
472 // If we get to here, we actually found the TGraph in the THistSvc.
473 // So let's add it to the local cache map and return its pointer
474 m_graphMap.insert( m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphPointer ) );
475 return graphPointer;
476 }
477
478
479 // Return the pointer to the TGraph that we got from the local cache map
480 return it->second;
481}
std::pair< StatusCode, TGraph * > getGraph(ITHistSvc &svc, const std::string &name)

◆ hash()

AthHistogramming::hash_t AthHistogramming::hash ( std::string_view histName) const
inlineprivateinherited

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

Definition at line 422 of file AthHistogramming.h.

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

◆ hist()

TH1 * AthHistogramming::hist ( std::string_view histName,
std::string_view tDir = "",
std::string_view stream = "" )
protectedinherited

Simplify the retrieval of registered histograms of any type.

Definition at line 182 of file AthHistogramming.cxx.

183 {
184
185 const auto histBaseName = splitName(histName).second;
186 std::string bookingString = buildBookingString(histName, tDir, stream, false);
187
188 const hash_t histHash = this->hash(histBaseName);
189
190 HistMap_t::const_iterator it = m_histMap.find(histHash);
191 if (it == m_histMap.end()) {
192 TH1* histPointer(nullptr);
193
194 if (!histSvc()->getHist(bookingString, histPointer).isSuccess()) {
195 std::string prefixedBookingString = buildBookingString(histName, tDir, stream, true);
196 if (!histSvc()->getHist(prefixedBookingString, histPointer).isSuccess()) {
197 m_msg << MSG::WARNING
198 << "Problem retrieving the histogram with name (including pre- and post-fixes) "
199 << prefixedBookingString
200 << " or with name " << histName
201 << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
202 << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
203 return nullptr;
204 }
205 }
206
207 m_histMap.insert(m_histMap.end(), std::pair<const hash_t, TH1*>(histHash, histPointer));
208 return histPointer;
209 }
210
211 return it->second;
212}
std::pair< StatusCode, TH1 * > getHist(ITHistSvc &svc, const std::string &name, size_t index=0)

◆ hist2d()

TH2 * AthHistogramming::hist2d ( std::string_view histName,
std::string_view tDir = "",
std::string_view stream = "" )
inlineprotectedinherited

Simplify the retrieval of registered 2-d histograms.

Definition at line 363 of file AthHistogramming.h.

364{
365 // Get the TH1 pointer
366 TH1* th1Pointer = this->hist(histName, tDir, stream);
367 if ( !th1Pointer )
368 {
369 m_msg << MSG::ERROR
370 << "Cannot get a 2-d histogram with name " << histName
371 << "... will probably seg-fault!" << endmsg;
372 return NULL;
373 }
374 // If the TH1 pointer is valid, simply return the dynamic_cast
375 return dynamic_cast<TH2*>( th1Pointer );
376}

◆ hist3d()

TH3 * AthHistogramming::hist3d ( std::string_view histName,
std::string_view tDir = "",
std::string_view stream = "" )
inlineprotectedinherited

Simplify the retrieval of registered 3-d histograms.

Definition at line 380 of file AthHistogramming.h.

381{
382 // Get the TH1 pointer
383 TH1* th1Pointer = this->hist(histName, tDir, stream);
384 if ( !th1Pointer )
385 {
386 m_msg << MSG::ERROR
387 << "Cannot get a 3-d histogram with name " << histName
388 << "... will probably seg-fault!" << endmsg;
389 return NULL;
390 }
391 // If the TH1 pointer is valid, simply return the dynamic_cast
392 return dynamic_cast<TH3*>( th1Pointer );
393}

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

◆ initialize()

StatusCode TruthHitAnalysis::initialize ( )
finaloverridevirtual

histograms declaration

now add branches and leaves to the tree

Definition at line 13 of file TruthHitAnalysis.cxx.

13 {
14 ATH_MSG_DEBUG( "Initializing TruthHitAnalysis" );
15
16 // Grab the Ntuple and histogramming service for the tree
17 ATH_CHECK(m_readKey.initialize());
18
20 m_h_n_vert = new TH1D("h_n_vert","n_vert", 100,200, 1500);
21 m_h_n_vert->StatOverflows();
22 ATH_CHECK(histSvc()->regHist(m_path + m_h_n_vert->GetName(), m_h_n_vert));
23
24 m_h_n_part = new TH1D("h_n_part","n_part", 100,1000, 10000);
25 m_h_n_part->StatOverflows();
26 ATH_CHECK(histSvc()->regHist(m_path + m_h_n_part->GetName(), m_h_n_part));
27
28 m_h_n_vert_prim = new TH1D("h_n_vert_prim","n_vert prim", 100,0, 1000);
29 m_h_n_vert_prim->StatOverflows();
30 ATH_CHECK(histSvc()->regHist(m_path + m_h_n_vert_prim->GetName(), m_h_n_vert_prim));
31
32 m_h_n_part_prim = new TH1D("h_n_part_prim","n_part prim", 100,200, 1500);
33 m_h_n_part_prim->StatOverflows();
34 ATH_CHECK(histSvc()->regHist(m_path + m_h_n_part_prim->GetName(), m_h_n_part_prim));
35
36 m_h_n_vert_sec = new TH1D("h_n_vert_sec","n_vert sec", 100,0, 1000);
37 m_h_n_vert_sec->StatOverflows();
38 ATH_CHECK(histSvc()->regHist(m_path + m_h_n_vert_sec->GetName(), m_h_n_vert_sec));
39
40 m_h_n_part_sec = new TH1D("h_n_part_sec","n_part sec", 100,0, 5000);
41 m_h_n_part_sec->StatOverflows();
42 ATH_CHECK(histSvc()->regHist(m_path + m_h_n_part_sec->GetName(), m_h_n_part_sec));
43
44 m_h_vtx_x = new TH1D("h_vtx_x","vtx_x", 100,-1300, 1300);
45 m_h_vtx_x->StatOverflows();
46 ATH_CHECK(histSvc()->regHist(m_path + m_h_vtx_x->GetName(), m_h_vtx_x));
47
48 m_h_vtx_y = new TH1D("h_vtx_y","vtx_y", 100,-1200, 1200);
49 m_h_vtx_y->StatOverflows();
50 ATH_CHECK(histSvc()->regHist(m_path + m_h_vtx_y->GetName(), m_h_vtx_y));
51
52 m_h_vtx_z = new TH1D("h_vtx_z","vtx_z", 100,-5000, 5000);
53 m_h_vtx_z->StatOverflows();
54 ATH_CHECK(histSvc()->regHist(m_path + m_h_vtx_z->GetName(), m_h_vtx_z));
55
56 m_h_vtx_r = new TH1D("h_vtx_r","vtx_r", 100,0, 1160);
57 m_h_vtx_r->StatOverflows();
58 ATH_CHECK(histSvc()->regHist(m_path + m_h_vtx_r->GetName(), m_h_vtx_r));
59
60 m_h_vtx_proc_x = new TH1D("h_vtx_proc_x","vtx_proc_x", 100,-1300, 1300);
61 m_h_vtx_proc_x->StatOverflows();
62 ATH_CHECK(histSvc()->regHist(m_path + m_h_vtx_proc_x->GetName(), m_h_vtx_proc_x));
63
64 m_h_vtx_proc_y = new TH1D("h_vtx_proc_y","vtx_proc_y", 100,-1200, 1200);
65 m_h_vtx_proc_y->StatOverflows();
66 ATH_CHECK(histSvc()->regHist(m_path + m_h_vtx_proc_y->GetName(), m_h_vtx_proc_y));
67
68 m_h_vtx_proc_z = new TH1D("h_vtx_proc_z","vtx_proc_z", 100,-5000, 5000);
69 m_h_vtx_proc_z->StatOverflows();
70 ATH_CHECK(histSvc()->regHist(m_path + m_h_vtx_proc_z->GetName(), m_h_vtx_proc_z));
71
72 m_h_vtx_proc_r = new TH1D("h_vtx_proc_r","vtx_proc_r", 100,0, 1160);
73 m_h_vtx_proc_r->StatOverflows();
74 ATH_CHECK(histSvc()->regHist(m_path + m_h_vtx_proc_r->GetName(), m_h_vtx_proc_r));
75
76 m_h_vtx_prim_x = new TH1D("h_vtx_prim_x","vtx_prim_x", 100,-1300, 1300);
77 m_h_vtx_prim_x->StatOverflows();
78 ATH_CHECK(histSvc()->regHist(m_path + m_h_vtx_prim_x->GetName(), m_h_vtx_prim_x));
79
80 m_h_vtx_prim_y = new TH1D("h_vtx_prim_y","vtx_prim_y", 100,-1200, 1200);
81 m_h_vtx_prim_y->StatOverflows();
82 ATH_CHECK(histSvc()->regHist(m_path + m_h_vtx_prim_y->GetName(), m_h_vtx_prim_y));
83
84 m_h_vtx_prim_z = new TH1D("h_vtx_prim_z","vtx_prim_z", 100,-5000, 5000);
85 m_h_vtx_prim_z->StatOverflows();
86 ATH_CHECK(histSvc()->regHist(m_path + m_h_vtx_prim_z->GetName(), m_h_vtx_prim_z));
87
88 m_h_vtx_prim_r = new TH1D("h_vtx_prim_r","vtx_prim_r", 100,0, 1160);
89 m_h_vtx_prim_r->StatOverflows();
90 ATH_CHECK(histSvc()->regHist(m_path + m_h_vtx_prim_r->GetName(), m_h_vtx_prim_r));
91
92 m_h_vtx_prim_xy = new TH2D("h_vtx_prim_xy","vtx_prim_xy", 100,-100, 100, 100,-100, 100);
93 m_h_vtx_prim_xy->StatOverflows();
94 ATH_CHECK(histSvc()->regHist(m_path + m_h_vtx_prim_xy->GetName(), m_h_vtx_prim_xy));
95
96 m_h_vtx_prim_zr = new TH2D("h_vtx_prim_zr","vtx_prim_zr", 100,-1500, 1500, 100,0, 150);
97 m_h_vtx_prim_zr->StatOverflows();
98 ATH_CHECK(histSvc()->regHist(m_path + m_h_vtx_prim_zr->GetName(), m_h_vtx_prim_zr));
99
100 m_h_vtx_sec_xy = new TH2D("h_vtx_sec_xy","vtx_sec_xy", 100,-1200, 1200, 100,-1200, 1200);
101 m_h_vtx_sec_xy->StatOverflows();
102 ATH_CHECK(histSvc()->regHist(m_path + m_h_vtx_sec_xy->GetName(), m_h_vtx_sec_xy));
103
104 m_h_vtx_sec_zr = new TH2D("h_vtx_sec_zr","vtx_sec_zr", 100,-6000, 6000, 100,0, 1160);
105 m_h_vtx_sec_zr->StatOverflows();
106 ATH_CHECK(histSvc()->regHist(m_path + m_h_vtx_sec_zr->GetName(), m_h_vtx_sec_zr));
107
108 m_h_n_generations = new TH1D("h_n_generations","h_generations", 100,0, 25);
109 m_h_n_generations->StatOverflows();
111
112 m_h_truth_px = new TH1D("h_truth_px","truth_px", 100,0, 4000);
113 m_h_truth_px->StatOverflows();
114 ATH_CHECK(histSvc()->regHist(m_path + m_h_truth_px->GetName(), m_h_truth_px));
115
116 m_h_truth_py = new TH1D("h_truth_py","truth_py", 100,0, 4000);
117 m_h_truth_py->StatOverflows();
118 ATH_CHECK(histSvc()->regHist(m_path + m_h_truth_py->GetName(), m_h_truth_py));
119
120 m_h_truth_pz = new TH1D("h_truth_pz","truth_pz", 100,0, 4000);
121 m_h_truth_pz->StatOverflows();
122 ATH_CHECK(histSvc()->regHist(m_path + m_h_truth_pz->GetName(), m_h_truth_pz));
123
124 m_h_truth_pt = new TH1D("h_truth_pt","truth_pt", 100,0, 4000);
125 m_h_truth_pt->StatOverflows();
126 ATH_CHECK(histSvc()->regHist(m_path + m_h_truth_pt->GetName(), m_h_truth_pt));
127
128 m_h_truth_eta = new TH1D("h_truth_eta","truth_eta", 50,-10, 10);
129 m_h_truth_eta->StatOverflows();
130 ATH_CHECK(histSvc()->regHist(m_path + m_h_truth_eta->GetName(), m_h_truth_eta));
131
132 m_h_truth_phi = new TH1D("h_truth_phi","truth_phi", 25,-3.1416, 3.1416);
133 m_h_truth_phi->StatOverflows();
134 ATH_CHECK(histSvc()->regHist(m_path + m_h_truth_phi->GetName(), m_h_truth_phi));
135
136 m_h_barcode = new TH1D("h_truth_barcode","truth_barcode", 100,0, 300000);
137 m_h_barcode->StatOverflows();
138 ATH_CHECK(histSvc()->regHist(m_path + m_h_barcode->GetName(), m_h_barcode));
139
140 m_h_part_status = new TH1D("h_part_status","part status", 100,0,50);
141 m_h_part_status->StatOverflows();
142 ATH_CHECK(histSvc()->regHist(m_path + m_h_part_status->GetName(), m_h_part_status));
143
144 m_h_part_pdgid = new TH1D("h_part_pdgid","part pdgid", 100,-5000, 5000);
145 m_h_part_pdgid->StatOverflows();
146 ATH_CHECK(histSvc()->regHist(m_path + m_h_part_pdgid->GetName(), m_h_part_pdgid));
147
148 m_h_part_pdgid_sec = new TH1D("h_part_pdgid_sec","part pdgid sec", 100,-5000, 5000);
149 m_h_part_pdgid_sec->StatOverflows();
151
152 m_h_part_eta = new TH1D("h_part_eta","part eta", 100,-10, 10);
153 m_h_part_eta->StatOverflows();
154 ATH_CHECK(histSvc()->regHist(m_path + m_h_part_eta->GetName(), m_h_part_eta));
155
156 m_h_part_phi = new TH1D("h_part_phi","part phi", 100,-3.2, 3.2);
157 m_h_part_phi->StatOverflows();
158 ATH_CHECK(histSvc()->regHist(m_path + m_h_part_phi->GetName(), m_h_part_phi));
159
160 m_h_part_p = new TH1D("h_part_p","part p", 100,0, 5000);
161 m_h_part_p->StatOverflows();
162 ATH_CHECK(histSvc()->regHist(m_path + m_h_part_p->GetName(), m_h_part_p));
163
165 m_tree = new TTree("Truth", "Truth");
166 std::string fullNtupleName = "/" + m_ntupleFileName + "/";
167 ATH_CHECK(histSvc()->regTree(fullNtupleName,m_tree));
168
169 m_tree->Branch("vtx_x", &m_vtx_x);
170 m_tree->Branch("vtx_y", &m_vtx_y);
171 m_tree->Branch("vtx_z", &m_vtx_z);
172 m_tree->Branch("vtx_r", &m_vtx_r);
173 m_tree->Branch("vtx_barcode", &m_vtx_barcode);
174 m_tree->Branch("vtx_proc_x", &m_vtx_proc_x);
175 m_tree->Branch("vtx_proc_y", &m_vtx_proc_y);
176 m_tree->Branch("vtx_proc_z", &m_vtx_proc_z);
177 m_tree->Branch("vtx_proc_r", &m_vtx_proc_r);
178 m_tree->Branch("vtx_proc_barcode", &m_vtx_proc_barcode);
179 m_tree->Branch("vtx_prim_x", &m_vtx_prim_x);
180 m_tree->Branch("vtx_prim_y", &m_vtx_prim_y);
181 m_tree->Branch("vtx_prim_z", &m_vtx_prim_z);
182 m_tree->Branch("vtx_prim_r", &m_vtx_prim_r);
183 m_tree->Branch("vtx_prim_barcode", &m_vtx_prim_barcode);
184 m_tree->Branch("truth_px", &m_truth_px);
185 m_tree->Branch("truth_py", &m_truth_py);
186 m_tree->Branch("truth_pz", &m_truth_pz);
187 m_tree->Branch("truth_pt", &m_truth_pt);
188 m_tree->Branch("truth_eta", &m_truth_eta);
189 m_tree->Branch("truth_phi", &m_truth_phi);
190 m_tree->Branch("barcode", &m_barcode);
191 m_tree->Branch("status", &m_status);
192 m_tree->Branch("pdg_id", &m_pdgid);
193
194 return StatusCode::SUCCESS;
195}
const ServiceHandle< ITHistSvc > & histSvc() const
The standard THistSvc (for writing histograms and TTrees and more to a root file) Returns (kind of) a...
Gaudi::Property< std::string > m_path
Gaudi::Property< std::string > m_ntupleFileName

◆ inputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< Gaudi::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.

◆ isClonable()

virtual bool AthCommonAlgorithm< Gaudi::Algorithm >::isClonable ( ) const
inlineoverridevirtualinherited

Specify if the algorithm is clonable.

Only relevant for non-reentrant algorithms. Actual number of clones needs to be set via the "Cardinality" property.

Reimplemented in AFP_DigiTop, AlgB, AlgT, BCM_Digitization, CscDigitBuilder, CscDigitToCscRDO, G4AtlasAlg, G4RunAlg, HGTD_Digitization, HiveAlgBase, InDet::GNNSeedingTrackMaker, InDet::SCT_Clusterization, InDet::SiSPGNNTrackMaker, InDet::SiSPSeededTrackFinder, InDet::SiTrackerSpacePointFinder, ISF::SimKernelMT, ITk::StripDigitization, ITkPixelCablingAlg, ITkStripCablingAlg, LArHitEMapMaker, LArTTL1Maker, LUCID_DigiTop, LVL1::L1TopoSimulation, MergeCalibHits, MergeGenericMuonSimHitColl, MergeHijingPars, MergeMcEventCollection, MergeTrackRecordCollection, MergeTruthJets, MergeTruthParticles, MuonDigitizer, PileUpMTAlg, PixelDigitization, RoIBResultToxAOD, SCT_ByteStreamErrorsTestAlg, SCT_CablingCondAlgFromCoraCool, SCT_CablingCondAlgFromText, SCT_ConditionsParameterTestAlg, SCT_ConditionsSummaryTestAlg, SCT_ConfigurationConditionsTestAlg, SCT_Digitization, SCT_FlaggedConditionTestAlg, SCT_LinkMaskingTestAlg, SCT_MajorityConditionsTestAlg, SCT_ModuleVetoTestAlg, SCT_MonitorConditionsTestAlg, SCT_PrepDataToxAOD, SCT_RawDataToxAOD, SCT_ReadCalibChipDataTestAlg, SCT_ReadCalibDataTestAlg, SCT_RODVetoTestAlg, SCT_SensorsTestAlg, SCT_SiliconConditionsTestAlg, SCT_StripVetoTestAlg, SCT_TdaqEnabledTestAlg, SCT_TestCablingAlg, SCTEventFlagWriter, SCTRawDataProvider, SCTSiLorentzAngleTestAlg, SCTSiPropertiesTestAlg, SGInputLoader, Simulation::BeamEffectsAlg, TileHitVecToCnt, TileMuonFitter, TilePulseForTileMuonReceiver, TileRawChannelMaker, TRTDigitization, and ZDC_DigiTop.

Definition at line 68 of file AthCommonAlgorithm.h.

68 {
69 return true;
70 }

◆ isReEntrant()

virtual bool AthAlgorithm::isReEntrant ( ) const
inlinefinaloverrideprotectedvirtualinherited

Legacy algorithms are not thread-safe.

Definition at line 47 of file AthAlgorithm.h.

47{ return false; }

◆ msg()

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

Definition at line 24 of file AthCommonMsg.h.

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

◆ msgLvl()

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

Definition at line 30 of file AthCommonMsg.h.

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

◆ outputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< Gaudi::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< Gaudi::Algorithm > >::renounce ( T & h)
inlineprotectedinherited

Definition at line 368 of file AthCommonDataStore.h.

369 {
370 h.renounce();
372 }
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< Gaudi::Algorithm > >::renounceArray ( SG::VarHandleKeyArray & handlesArray)
inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 352 of file AthCommonDataStore.h.

352 {
354 }

◆ setFilterPassed()

virtual void AthCommonAlgorithm< Gaudi::Algorithm >::setFilterPassed ( bool state,
const EventContext & ctx ) const
inlinevirtualinherited

Set filter decision:

Definition at line 99 of file AthCommonAlgorithm.h.

99 {
101 }
virtual void setFilterPassed(bool state, const EventContext &ctx) const
Set filter decision:

◆ sysExecute()

StatusCode AthCommonAlgorithm< Gaudi::Algorithm >::sysExecute ( const EventContext & ctx)
overridevirtualinherited

Execute an algorithm.

We override this in order to work around an issue with the Algorithm base class storing the event context in a member variable that can cause crashes in MT jobs.

Reimplemented in AthAnalysisAlgorithm.

Definition at line 80 of file AthCommonAlgorithm.cxx.

41{
42 return BaseAlg::sysExecute (ctx);
43}

◆ 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 AthCommonAlgorithm< Gaudi::Algorithm >.

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
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 to fill the internal configuration.

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Gaudi::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 ( std::string_view treeName,
std::string_view tDir = "",
std::string_view stream = "" )
protectedinherited

Simplify the retrieval of registered TTrees.

Definition at line 333 of file AthHistogramming.cxx.

334{
335 // Build a 32 bit hash out of the name
336 const hash_t treeHash = this->hash(treeName);
337 // See if this entry exists in the map
338 TreeMap_t::const_iterator it = m_treeMap.find( treeHash );
339 if ( it == m_treeMap.end() ) // It doesn't exist!
340 { // Let's see into the THistSvc if somebody else has registered the TTree...
341 // Massage the final string to book things
342 std::string bookingString = buildBookingString( treeName, tDir, stream);
343
344 TTree* treePointer(NULL);
345 if ( !((histSvc()->getTree(bookingString, treePointer)).isSuccess()) )
346 {
347 m_msg << MSG::WARNING
348 << "Problem retrieving the TTree with name " << treeName
349 << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
350 << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
351 return NULL;
352 }
353 // If we get to here, we actually found the TTree in the THistSvc.
354 // So let's add it to the local cache map and return its pointer
355 m_treeMap.insert( m_treeMap.end(), std::pair< const hash_t, TTree* >( treeHash, treePointer ) );
356 return treePointer;
357 }
358
359 // Return the pointer to the TTree that we got from the local cache map
360 return it->second;
361}
std::pair< StatusCode, TTree * > getTree(ITHistSvc &svc, const std::string &name)

◆ updateVHKA()

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

Definition at line 298 of file AthCommonDataStore.h.

298 {
299 for (auto &a : m_vhka) {
301 for (auto k : keys) {
302 k->setOwner(this);
303 }
304 }
305 }

Member Data Documentation

◆ m_barcode

std::vector<float>* TruthHitAnalysis::m_barcode {nullptr}
private

Definition at line 84 of file TruthHitAnalysis.h.

84{nullptr};

◆ m_detStore

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

Pointer to StoreGate (detector store by default).

Definition at line 381 of file AthCommonDataStore.h.

◆ m_effMap

EffMap_t AthHistogramming::m_effMap
privateinherited

The map of histogram names to their pointers.

Definition at line 202 of file AthHistogramming.h.

◆ m_evtStore

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

Pointer to StoreGate (event store by default).

Definition at line 378 of file AthCommonDataStore.h.

◆ m_extendedExtraObjects

DataObjIDColl AthCommonAlgorithm< Gaudi::Algorithm >::m_extendedExtraObjects
privateinherited

Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.

Empty if no symlinks were found.

Definition at line 108 of file AthCommonAlgorithm.h.

◆ m_graphMap

GraphMap_t AthHistogramming::m_graphMap
privateinherited

The map of TGraph names to their pointers.

Definition at line 216 of file AthHistogramming.h.

◆ m_h_barcode

TH1* TruthHitAnalysis::m_h_barcode {nullptr}
private

Definition at line 55 of file TruthHitAnalysis.h.

55{nullptr};

◆ m_h_n_generations

TH1* TruthHitAnalysis::m_h_n_generations {nullptr}
private

Definition at line 48 of file TruthHitAnalysis.h.

48{nullptr};

◆ m_h_n_part

TH1* TruthHitAnalysis::m_h_n_part {nullptr}
private

Definition at line 27 of file TruthHitAnalysis.h.

27{nullptr};

◆ m_h_n_part_prim

TH1* TruthHitAnalysis::m_h_n_part_prim {nullptr}
private

Definition at line 29 of file TruthHitAnalysis.h.

29{nullptr};

◆ m_h_n_part_sec

TH1* TruthHitAnalysis::m_h_n_part_sec {nullptr}
private

Definition at line 31 of file TruthHitAnalysis.h.

31{nullptr};

◆ m_h_n_vert

TH1* TruthHitAnalysis::m_h_n_vert {nullptr}
private

Some variables.

Definition at line 26 of file TruthHitAnalysis.h.

26{nullptr};

◆ m_h_n_vert_prim

TH1* TruthHitAnalysis::m_h_n_vert_prim {nullptr}
private

Definition at line 28 of file TruthHitAnalysis.h.

28{nullptr};

◆ m_h_n_vert_sec

TH1* TruthHitAnalysis::m_h_n_vert_sec {nullptr}
private

Definition at line 30 of file TruthHitAnalysis.h.

30{nullptr};

◆ m_h_part_eta

TH1* TruthHitAnalysis::m_h_part_eta {nullptr}
private

Definition at line 59 of file TruthHitAnalysis.h.

59{nullptr};

◆ m_h_part_p

TH1* TruthHitAnalysis::m_h_part_p {nullptr}
private

Definition at line 61 of file TruthHitAnalysis.h.

61{nullptr};

◆ m_h_part_pdgid

TH1* TruthHitAnalysis::m_h_part_pdgid {nullptr}
private

Definition at line 57 of file TruthHitAnalysis.h.

57{nullptr};

◆ m_h_part_pdgid_sec

TH1* TruthHitAnalysis::m_h_part_pdgid_sec {nullptr}
private

Definition at line 58 of file TruthHitAnalysis.h.

58{nullptr};

◆ m_h_part_phi

TH1* TruthHitAnalysis::m_h_part_phi {nullptr}
private

Definition at line 60 of file TruthHitAnalysis.h.

60{nullptr};

◆ m_h_part_status

TH1* TruthHitAnalysis::m_h_part_status {nullptr}
private

Definition at line 56 of file TruthHitAnalysis.h.

56{nullptr};

◆ m_h_truth_eta

TH1* TruthHitAnalysis::m_h_truth_eta {nullptr}
private

Definition at line 53 of file TruthHitAnalysis.h.

53{nullptr};

◆ m_h_truth_phi

TH1* TruthHitAnalysis::m_h_truth_phi {nullptr}
private

Definition at line 54 of file TruthHitAnalysis.h.

54{nullptr};

◆ m_h_truth_pt

TH1* TruthHitAnalysis::m_h_truth_pt {nullptr}
private

Definition at line 52 of file TruthHitAnalysis.h.

52{nullptr};

◆ m_h_truth_px

TH1* TruthHitAnalysis::m_h_truth_px {nullptr}
private

Definition at line 49 of file TruthHitAnalysis.h.

49{nullptr};

◆ m_h_truth_py

TH1* TruthHitAnalysis::m_h_truth_py {nullptr}
private

Definition at line 50 of file TruthHitAnalysis.h.

50{nullptr};

◆ m_h_truth_pz

TH1* TruthHitAnalysis::m_h_truth_pz {nullptr}
private

Definition at line 51 of file TruthHitAnalysis.h.

51{nullptr};

◆ m_h_vtx_prim_r

TH1* TruthHitAnalysis::m_h_vtx_prim_r {nullptr}
private

Definition at line 43 of file TruthHitAnalysis.h.

43{nullptr};

◆ m_h_vtx_prim_x

TH1* TruthHitAnalysis::m_h_vtx_prim_x {nullptr}
private

Definition at line 40 of file TruthHitAnalysis.h.

40{nullptr};

◆ m_h_vtx_prim_xy

TH2* TruthHitAnalysis::m_h_vtx_prim_xy {nullptr}
private

Definition at line 44 of file TruthHitAnalysis.h.

44{nullptr};

◆ m_h_vtx_prim_y

TH1* TruthHitAnalysis::m_h_vtx_prim_y {nullptr}
private

Definition at line 41 of file TruthHitAnalysis.h.

41{nullptr};

◆ m_h_vtx_prim_z

TH1* TruthHitAnalysis::m_h_vtx_prim_z {nullptr}
private

Definition at line 42 of file TruthHitAnalysis.h.

42{nullptr};

◆ m_h_vtx_prim_zr

TH2* TruthHitAnalysis::m_h_vtx_prim_zr {nullptr}
private

Definition at line 45 of file TruthHitAnalysis.h.

45{nullptr};

◆ m_h_vtx_proc_r

TH1* TruthHitAnalysis::m_h_vtx_proc_r {nullptr}
private

Definition at line 39 of file TruthHitAnalysis.h.

39{nullptr};

◆ m_h_vtx_proc_x

TH1* TruthHitAnalysis::m_h_vtx_proc_x {nullptr}
private

Definition at line 36 of file TruthHitAnalysis.h.

36{nullptr};

◆ m_h_vtx_proc_y

TH1* TruthHitAnalysis::m_h_vtx_proc_y {nullptr}
private

Definition at line 37 of file TruthHitAnalysis.h.

37{nullptr};

◆ m_h_vtx_proc_z

TH1* TruthHitAnalysis::m_h_vtx_proc_z {nullptr}
private

Definition at line 38 of file TruthHitAnalysis.h.

38{nullptr};

◆ m_h_vtx_r

TH1* TruthHitAnalysis::m_h_vtx_r {nullptr}
private

Definition at line 35 of file TruthHitAnalysis.h.

35{nullptr};

◆ m_h_vtx_sec_xy

TH2* TruthHitAnalysis::m_h_vtx_sec_xy {nullptr}
private

Definition at line 46 of file TruthHitAnalysis.h.

46{nullptr};

◆ m_h_vtx_sec_zr

TH2* TruthHitAnalysis::m_h_vtx_sec_zr {nullptr}
private

Definition at line 47 of file TruthHitAnalysis.h.

47{nullptr};

◆ m_h_vtx_x

TH1* TruthHitAnalysis::m_h_vtx_x {nullptr}
private

Definition at line 32 of file TruthHitAnalysis.h.

32{nullptr};

◆ m_h_vtx_y

TH1* TruthHitAnalysis::m_h_vtx_y {nullptr}
private

Definition at line 33 of file TruthHitAnalysis.h.

33{nullptr};

◆ m_h_vtx_z

TH1* TruthHitAnalysis::m_h_vtx_z {nullptr}
private

Definition at line 34 of file TruthHitAnalysis.h.

34{nullptr};

◆ m_histMap

HistMap_t AthHistogramming::m_histMap
privateinherited

The map of histogram names to their pointers.

Definition at line 195 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_msg

MsgStream AthHistogramming::m_msg
privateinherited

Cached Message Stream.

Definition at line 242 of file AthHistogramming.h.

◆ m_name

std::string AthHistogramming::m_name
privateinherited

Instance name.

Definition at line 239 of file AthHistogramming.h.

◆ m_ntupleFileName

Gaudi::Property<std::string> TruthHitAnalysis::m_ntupleFileName {this, "NtupleFileName","/TruthHitAnalysis/"}
private

Definition at line 91 of file TruthHitAnalysis.h.

91{this, "NtupleFileName","/TruthHitAnalysis/"};

◆ m_path

Gaudi::Property<std::string> TruthHitAnalysis::m_path {this, "HistPath","/TruthHitAnalysis/"}
private

Definition at line 90 of file TruthHitAnalysis.h.

90{this, "HistPath","/TruthHitAnalysis/"};

◆ m_pdgid

std::vector<float>* TruthHitAnalysis::m_pdgid {nullptr}
private

Definition at line 86 of file TruthHitAnalysis.h.

86{nullptr};

◆ m_prefix

std::string AthHistogramAlgorithm::m_prefix
privateinherited

Name of the ROOT output stream (file).

Definition at line 88 of file AthHistogramAlgorithm.h.

◆ m_readKey

SG::ReadHandleKey<McEventCollection> TruthHitAnalysis::m_readKey {this, "InputKey", "TruthEvent"}
private

Definition at line 92 of file TruthHitAnalysis.h.

92{this, "InputKey", "TruthEvent"};

◆ m_rootDir

std::string AthHistogramAlgorithm::m_rootDir
privateinherited

Name of the ROOT directory.

Definition at line 91 of file AthHistogramAlgorithm.h.

◆ m_status

std::vector<float>* TruthHitAnalysis::m_status {nullptr}
private

Definition at line 85 of file TruthHitAnalysis.h.

85{nullptr};

◆ m_streamName

std::string AthHistogramming::m_streamName
privateinherited

Name of the ROOT output stream (file).

Definition at line 220 of file AthHistogramming.h.

◆ m_tree

TTree* TruthHitAnalysis::m_tree {nullptr}
private

Definition at line 88 of file TruthHitAnalysis.h.

88{nullptr};

◆ m_treeMap

TreeMap_t AthHistogramming::m_treeMap
privateinherited

The map of TTree names to their pointers.

Definition at line 209 of file AthHistogramming.h.

◆ m_truth_eta

std::vector<float>* TruthHitAnalysis::m_truth_eta {nullptr}
private

Definition at line 82 of file TruthHitAnalysis.h.

82{nullptr};

◆ m_truth_phi

std::vector<float>* TruthHitAnalysis::m_truth_phi {nullptr}
private

Definition at line 83 of file TruthHitAnalysis.h.

83{nullptr};

◆ m_truth_pt

std::vector<float>* TruthHitAnalysis::m_truth_pt {nullptr}
private

Definition at line 81 of file TruthHitAnalysis.h.

81{nullptr};

◆ m_truth_px

std::vector<float>* TruthHitAnalysis::m_truth_px {nullptr}
private

Definition at line 78 of file TruthHitAnalysis.h.

78{nullptr};

◆ m_truth_py

std::vector<float>* TruthHitAnalysis::m_truth_py {nullptr}
private

Definition at line 79 of file TruthHitAnalysis.h.

79{nullptr};

◆ m_truth_pz

std::vector<float>* TruthHitAnalysis::m_truth_pz {nullptr}
private

Definition at line 80 of file TruthHitAnalysis.h.

80{nullptr};

◆ m_varHandleArraysDeclared

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

Definition at line 387 of file AthCommonDataStore.h.

◆ m_vhka

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

Definition at line 386 of file AthCommonDataStore.h.

◆ m_vtx_barcode

std::vector<float>* TruthHitAnalysis::m_vtx_barcode {nullptr}
private

Definition at line 67 of file TruthHitAnalysis.h.

67{nullptr};

◆ m_vtx_prim_barcode

std::vector<float>* TruthHitAnalysis::m_vtx_prim_barcode {nullptr}
private

Definition at line 77 of file TruthHitAnalysis.h.

77{nullptr};

◆ m_vtx_prim_r

std::vector<float>* TruthHitAnalysis::m_vtx_prim_r {nullptr}
private

Definition at line 76 of file TruthHitAnalysis.h.

76{nullptr};

◆ m_vtx_prim_x

std::vector<float>* TruthHitAnalysis::m_vtx_prim_x {nullptr}
private

Definition at line 73 of file TruthHitAnalysis.h.

73{nullptr};

◆ m_vtx_prim_y

std::vector<float>* TruthHitAnalysis::m_vtx_prim_y {nullptr}
private

Definition at line 74 of file TruthHitAnalysis.h.

74{nullptr};

◆ m_vtx_prim_z

std::vector<float>* TruthHitAnalysis::m_vtx_prim_z {nullptr}
private

Definition at line 75 of file TruthHitAnalysis.h.

75{nullptr};

◆ m_vtx_proc_barcode

std::vector<float>* TruthHitAnalysis::m_vtx_proc_barcode {nullptr}
private

Definition at line 72 of file TruthHitAnalysis.h.

72{nullptr};

◆ m_vtx_proc_r

std::vector<float>* TruthHitAnalysis::m_vtx_proc_r {nullptr}
private

Definition at line 71 of file TruthHitAnalysis.h.

71{nullptr};

◆ m_vtx_proc_x

std::vector<float>* TruthHitAnalysis::m_vtx_proc_x {nullptr}
private

Definition at line 68 of file TruthHitAnalysis.h.

68{nullptr};

◆ m_vtx_proc_y

std::vector<float>* TruthHitAnalysis::m_vtx_proc_y {nullptr}
private

Definition at line 69 of file TruthHitAnalysis.h.

69{nullptr};

◆ m_vtx_proc_z

std::vector<float>* TruthHitAnalysis::m_vtx_proc_z {nullptr}
private

Definition at line 70 of file TruthHitAnalysis.h.

70{nullptr};

◆ m_vtx_r

std::vector<float>* TruthHitAnalysis::m_vtx_r {nullptr}
private

Definition at line 66 of file TruthHitAnalysis.h.

66{nullptr};

◆ m_vtx_x

std::vector<float>* TruthHitAnalysis::m_vtx_x {nullptr}
private

Definition at line 63 of file TruthHitAnalysis.h.

63{nullptr};

◆ m_vtx_y

std::vector<float>* TruthHitAnalysis::m_vtx_y {nullptr}
private

Definition at line 64 of file TruthHitAnalysis.h.

64{nullptr};

◆ m_vtx_z

std::vector<float>* TruthHitAnalysis::m_vtx_z {nullptr}
private

Definition at line 65 of file TruthHitAnalysis.h.

65{nullptr};

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