ATLAS Offline Software
Loading...
Searching...
No Matches
MuonR4::BucketDumperAlg Class Reference

#include <BucketDumperAlg.h>

Inheritance diagram for MuonR4::BucketDumperAlg:

Public Member Functions

 ~BucketDumperAlg ()=default
virtual StatusCode initialize () override final
virtual StatusCode finalize () override final
virtual StatusCode execute (const EventContext &ctx) override final
 Execute method.
 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 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

StatusCode dumpContainer (const EventContext &ctx, const SG::ReadHandleKey< SpacePointContainer > &spacePointKey, const SG::ReadHandleKey< xAOD::MuonSegmentContainer > &segmentKey)
 Dumps the space point container with the associated muon segment container.
StatusCode computeAllBucketScores (const EventContext &ctx, const SpacePointContainer *spContainer, std::unordered_map< const SpacePointBucket *, std::vector< float > > &bucketScoreMap) const
 Computes ML bucket scores using the inference tool for a specific container.
CLHEP::HepRandomEngine * getRandomEngine (const EventContext &ctx) const
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::ReadHandleKeyArray< SpacePointContainerm_spacePointKeys
SG::ReadHandleKeyArray< xAOD::MuonSegmentContainerm_inSegmentKeys {this, "SegmentKey", {"MuonSegmentsFromR4"}}
SG::ReadDecorHandleKeyArray< xAOD::MuonSegmentContainerm_truthDecorKeys {this, "TruthDecorLinks", {}}
ActsTrk::GeoContextReadKey_t m_geoCtxKey {this, "AlignmentKey", "ActsAlignment", "cond handle key"}
Gaudi::Property< bool > m_isMC {this, "isMC", true}
Gaudi::Property< double > m_fracToKeep {this,"dataFracToKeep", 1.}
Gaudi::Property< std::string > m_streamName {this, "StreamName", ""}
Gaudi::Property< bool > m_doMLBucketScore
Gaudi::Property< bool > m_doMLBucketFilter
ServiceHandle< IAthRNGSvcm_rndmSvc {this, "RndmSvc", "AthRNGSvc", ""}
ToolHandle< MuonML::IGraphInferenceToolm_inferenceTool
 Inference tool for ML bucket scoring (optional).
ToolHandle< MuonValR4::IPatternVisualizationToolm_visionTool {this, "VisualizationTool", ""}
 Pattern visualization tool.
MuonVal::MuonTesterTree m_tree {"MuonBucketDump","MuonBucketDump"}
MuonVal::ScalarBranch< float > & m_bucket_min {m_tree.newScalar<float>("bucket_min", -1)}
MuonVal::ScalarBranch< float > & m_bucket_max {m_tree.newScalar<float>("bucket_max", -1)}
MuonVal::ScalarBranch< uint16_t > & m_bucket_spacePoints {m_tree.newScalar<uint16_t>("bucket_spacePoints", 0)}
MuonVal::ScalarBranch< uint16_t > & m_bucket_segments {m_tree.newScalar<uint16_t>("bucket_segments", 0)}
MuonVal::ScalarBranch< uint16_t > & m_bucket_layers {m_tree.newScalar<uint16_t>("bucket_layers", 0)}
MuonVal::ScalarBranch< uint8_t > & m_bucket_truthHit {m_tree.newScalar<uint8_t>("bucket_hasTruth", 0)}
MuonVal::ScalarBranch< float > & m_bucket_posX {m_tree.newScalar<float>("bucket_positionX")}
MuonVal::ScalarBranch< float > & m_bucket_posY {m_tree.newScalar<float>("bucket_positionY")}
MuonVal::ScalarBranch< float > & m_bucket_posZ {m_tree.newScalar<float>("bucket_positionZ")}
MuonVal::ScalarBranch< uint8_t > & m_bucket_sector {m_tree.newScalar<uint8_t>("bucket_sector")}
MuonVal::ScalarBranch< uint8_t > & m_bucket_chamberIdx {m_tree.newScalar<uint8_t>("bucket_chamberIndex")}
MuonVal::ScalarBranch< Char_t > & m_bucket_side {m_tree.newScalar<Char_t>("bucket_side")}
MuonVal::VectorBranch< float > & m_bucket_ml_score_class0 {m_tree.newVector<float>("bucket_ml_score_class0", false)}
 ML bucket filter scores (3 classes for the filter model).
MuonVal::VectorBranch< float > & m_bucket_ml_score_class1 {m_tree.newVector<float>("bucket_ml_score_class1", false)}
MuonVal::VectorBranch< float > & m_bucket_ml_score_class2 {m_tree.newVector<float>("bucket_ml_score_class2", false)}
MuonVal::ThreeVectorBranch m_spoint_localPosition {m_tree, "localPosition"}
MuonVal::ThreeVectorBranch m_spoint_globalPosition {m_tree, "globalPosition"}
MuonVal::MuonIdentifierBranch m_spoint_id {m_tree, "id"}
MuonVal::VectorBranch< uint16_t > & m_spoint_layer {m_tree.newVector<uint16_t>("Layer")}
MuonVal::VectorBranch< unsigned short > & m_spoint_isStrip {m_tree.newVector<unsigned short>("isStrip", false)}
MuonVal::VectorBranch< unsigned short > & m_spoint_isMdt {m_tree.newVector<unsigned short>("isMdt", false)}
MuonVal::VectorBranch< uint16_t > & m_spoint_adc {m_tree.newVector<uint16_t>("adc")}
MuonVal::VectorBranch< uint16_t > & m_spoint_tdc {m_tree.newVector<uint16_t>("tdc")}
MuonVal::VectorBranch< float > & m_spoint_covX {m_tree.newVector<float>("covX")}
MuonVal::VectorBranch< float > & m_spoint_covY {m_tree.newVector<float>("covY")}
MuonVal::VectorBranch< float > & m_spoint_driftR {m_tree.newVector<float>("driftR")}
MuonVal::VectorBranch< unsigned short > & m_spoint_measuresEta {m_tree.newVector<unsigned short>("measuresEta")}
MuonVal::VectorBranch< unsigned short > & m_spoint_measuresPhi {m_tree.newVector<unsigned short>("measuresPhi")}
MuonVal::VectorBranch< unsigned short > & m_spoint_trueLabel {m_tree.newVector<unsigned short>("trueLabel")}
MuonVal::VectorBranch< unsigned int > & m_spoint_nEtaInstances {m_tree.newVector<unsigned int>("nEtaInUse")}
MuonVal::VectorBranch< unsigned int > & m_spoint_nPhiInstances {m_tree.newVector<unsigned int>("nPhiInUse")}
MuonVal::VectorBranch< unsigned int > & m_spoint_dimension {m_tree.newVector<unsigned int>("dimension")}
MuonVal::VectorBranch< uint16_t > & m_spoint_nSegments {m_tree.newVector<uint16_t>("nSegments")}
MuonVal::MatrixBranch< int16_t > & m_spoint_mat {m_tree.newMatrix<int16_t>("sp_seg_matching",-1)}
MuonVal::MatrixBranch< int16_t > & m_spoint_trueSeg {m_tree.newMatrix<int16_t>("spacePoint_truthSegLink",-1)}
MuonVal::VectorBranch< float > & m_segmentLocX {m_tree.newVector<float>("segmentLocalX")}
MuonVal::VectorBranch< float > & m_segmentLocY {m_tree.newVector<float>("segmentLocalY")}
MuonVal::VectorBranch< float > & m_segmentLocTheta {m_tree.newVector<float>("segmentLocalTheta")}
MuonVal::VectorBranch< float > & m_segmentLocPhi {m_tree.newVector<float>("segmentLocalPhi")}
MuonVal::VectorBranch< uint16_t > & m_segmentTruthIdx {m_tree.newVector<uint16_t>("segmentTruthPart")}
MuonVal::ThreeVectorBranch m_segmentPos {m_tree, "segmentPosition"}
MuonVal::ThreeVectorBranch m_segmentDir {m_tree, "segmentDirection"}
MuonVal::VectorBranch< float > & m_segment_chiSquared {m_tree.newVector<float>("segment_chiSquared")}
MuonVal::VectorBranch< float > & m_segment_numberDoF {m_tree.newVector<float>("segment_numberDoF")}
MuonVal::VectorBranch< float > & m_truthSegLocX {m_tree.newVector<float>("truthSegmentLocalX")}
MuonVal::VectorBranch< float > & m_truthSegLocY {m_tree.newVector<float>("truthSegmentLocalY")}
MuonVal::VectorBranch< float > & m_truthSegLocTheta {m_tree.newVector<float>("truthSegmentLocalTheta")}
MuonVal::VectorBranch< float > & m_truthSegLocPhi {m_tree.newVector<float>("truthSegmentLocalPhi")}
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_name
 Instance name.
MsgStream m_msg
 Cached Message Stream.

Detailed Description

Definition at line 35 of file BucketDumperAlg.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< Gaudi::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

◆ ~BucketDumperAlg()

MuonR4::BucketDumperAlg::~BucketDumperAlg ( )
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,
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 388 of file AthHistogramming.cxx.

389{
390 // Get a pointer
391 const TGraph* graphPointer = &graphRef;
392
393 // Check that we got a valid pointer
394 if ( !graphPointer )
395 {
396 m_msg << MSG::WARNING
397 << "We got an invalid TGraph pointer in the BookGetPointer(TGraph*) method of the class" << m_name
398 << "!" << endmsg;
399 return NULL;
400 }
401
402 // Modify the name and title according to the prefixes of this classes instance
403 std::string graphName = graphPointer->GetName();
404 const std::string graphTitle = graphPointer->GetTitle();
405
406 // Check if the hash for this graphName already exists, i.e., if we have a hash collision
407 const hash_t graphHash = this->hash(graphName);
408 GraphMap_t::const_iterator it = m_graphMap.find( graphHash );
409 if ( it != m_graphMap.end() ) // It does exist!
410 {
411 m_msg << MSG::WARNING
412 << "Detected a hash collision. The hash for the TGraph with name=" << graphName
413 << " already exists and points to a TGraph with name=" << it->second->GetName()
414 << " NOT going to book the new histogram and returning a NULL pointer!" << endmsg;
415 return NULL;
416 }
417
418 // Create a clone that has the new name
419 TGraph* graphClone = dynamic_cast< TGraph* >( graphPointer->Clone((m_histNamePrefix+graphName+m_histNamePostfix).c_str()) );
420 if( !graphClone )
421 {
422 m_msg << MSG::WARNING
423 << "We couldn't clone the TGraph in the BookGetPointer(TGraph&) method of the class" << m_name
424 << "!" << endmsg;
425 return NULL;
426 }
427 graphClone->SetTitle ((m_histTitlePrefix+graphTitle+m_histTitlePostfix).c_str());
428
429 // Massage the final string to book things
430 std::string bookingString("");
431 this->buildBookingString( bookingString, graphName, tDir, stream );
432
433 // Register the TGraph into the THistSvc
434 if ( !((histSvc()->regGraph(bookingString, graphClone)).isSuccess()) )
435 {
436 m_msg << MSG::WARNING
437 << "Problem registering TGraph with name " << graphName
438 << ", title " << graphTitle
439 << " in " << m_name << "!" << endmsg;
440 return NULL;
441 }
442
443 // Also register it in the local map of string to pointer
444 m_graphMap.insert( m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphClone ) );
445
446 return graphClone;
447}
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 273 of file AthHistogramming.cxx.

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

124{
125 std::string effName(effRef.GetName());
126 const std::string effTitle(effRef.GetTitle());
127 std::string bookingString;
128 this->buildBookingString(bookingString, effName, tDir, stream);
129 const std::string finalEffName = m_histNamePrefix + effName + m_histNamePostfix;
130 effRef.SetTitle((m_histTitlePrefix + effTitle + m_histTitlePostfix).c_str());
131 effRef.SetName(finalEffName.c_str());
132
133 const hash_t effHash = this->hash(effName);
134 EffMap_t::const_iterator it = m_effMap.find(effHash);
135 if (it != m_effMap.end()) {
136 m_msg << MSG::WARNING
137 << "Detected a hash collision. The hash for the TEfficiency with name=" << effName
138 << " already exists and points to a TEfficiency with name=" << it->second->GetName()
139 << " NOT going to book the new TEfficiency and returning a NULL pointer!" << endmsg;
140 return nullptr;
141 }
142
143 if (!histSvc()->regEfficiency(bookingString, &effRef).isSuccess()) {
144 m_msg << MSG::WARNING
145 << "Problem registering TEfficiency with name " << effName
146 << ", name prefix " << m_histNamePrefix
147 << ", title " << effTitle
148 << ", title prefix " << m_histTitlePrefix
149 << ", and title postfix " << m_histTitlePostfix
150 << " in " << m_name << "!" << endmsg;
151 return nullptr;
152 }
153
154 m_effMap.insert(m_effMap.end(), std::pair<const hash_t, TEfficiency*>(effHash, &effRef));
155
156 return &effRef;
157}
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 std::string histName(histRef.GetName());
86 const std::string histTitle(histRef.GetTitle());
87 std::string bookingString;
88
89 this->buildBookingString(bookingString, histName, tDir, stream);
90
91 const std::string finalHistName = m_histNamePrefix + histName + m_histNamePostfix;
92
93 histRef.SetTitle((m_histTitlePrefix + histTitle + m_histTitlePostfix).c_str());
94 histRef.SetName(finalHistName.c_str());
95
96 const hash_t histHash = this->hash(histName);
97 HistMap_t::const_iterator it = m_histMap.find(histHash);
98 if (it != m_histMap.end()) {
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 nullptr;
104 }
105
106 if (!histSvc()->regHist(bookingString, &histRef).isSuccess()) {
107 m_msg << MSG::WARNING
108 << "Problem registering histogram with name " << histName
109 << ", name prefix " << m_histNamePrefix
110 << ", title " << histTitle
111 << ", title prefix " << m_histTitlePrefix
112 << ", and title postfix " << m_histTitlePostfix
113 << " in " << m_name << "!" << endmsg;
114 return nullptr;
115 }
116
117 m_histMap.insert(m_histMap.end(), std::pair<const hash_t, TH1*>(histHash, &histRef));
118
119 return &histRef;
120}
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 521 of file AthHistogramming.cxx.

526{
527 // Massage the final string to book things
528 if(tDir.empty()) tDir = m_rootDir;
529 size_t pos = histName.rfind('/');
530 if(pos != std::string::npos){
531 tDir+='/';
532 tDir.append(histName, 0,pos);
533 histName.erase(0,pos+1);
534 };
535 if(stream.empty()) stream = m_streamName;
536
537 if(usePrefixPostfix){
538 bookingString = "/"+stream+"/"+tDir+"/"+m_histNamePrefix+histName+m_histNamePostfix;
539 } else {
540 bookingString = "/"+stream+"/"+tDir+"/"+histName;
541 }
542 while(bookingString.find("//") != std::string::npos){
543 this->myReplace(bookingString,"//","/");
544 }
545
546 return;
547}
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.

◆ computeAllBucketScores()

StatusCode MuonR4::BucketDumperAlg::computeAllBucketScores ( const EventContext & ctx,
const SpacePointContainer * spContainer,
std::unordered_map< const SpacePointBucket *, std::vector< float > > & bucketScoreMap ) const
private

Computes ML bucket scores using the inference tool for a specific container.

Definition at line 338 of file BucketDumperAlg.cxx.

340 {
341 bucketScoreMap.clear();
342
343 ATH_MSG_DEBUG("computeAllBucketScores: doMLBucketScore=" << m_doMLBucketScore
344 << ", tool empty=" << m_inferenceTool.empty()
345 << ", spContainer=" << (spContainer ? "valid" : "null"));
346
347 if (!m_doMLBucketScore || m_inferenceTool.empty() || !spContainer) {
348 return StatusCode::SUCCESS;
349 }
350
351 // Count non-empty buckets
352 size_t nonEmptyBuckets = 0;
353 for (const MuonR4::SpacePointBucket* bucket : *spContainer) {
354 if (bucket && !bucket->empty()) {
355 ++nonEmptyBuckets;
356 }
357 }
358
359 // If no non-empty buckets, skip inference
360 if (nonEmptyBuckets == 0) {
361 ATH_MSG_DEBUG("Container has no non-empty buckets, skipping ML scoring");
362 return StatusCode::SUCCESS;
363 }
364
365 // Run the inference tool
366 ATH_MSG_DEBUG("Calling runGraphInference on tool: " << m_inferenceTool.name());
367 MuonML::GraphRawData graphData{};
368 ATH_CHECK(m_inferenceTool->runGraphInference(ctx, graphData));
369 ATH_MSG_DEBUG("runGraphInference completed successfully");
370
371 // Extract logits from the inference output
372 if (!graphData.graph || graphData.graph->dataTensor.size() <= 2) {
373 ATH_MSG_ERROR("Missing output logits tensor at index 2");
374 return StatusCode::FAILURE;
375 }
376
377 const Ort::Value& outTensor = graphData.graph->dataTensor[2];
378 const auto& info = outTensor.GetTensorTypeAndShapeInfo();
379 std::vector<int64_t> outShape = info.GetShape();
380
381 if (outShape.size() != 2 || outShape[1] != 3) {
382 ATH_MSG_ERROR("Unexpected ONNX output tensor shape");
383 return StatusCode::FAILURE;
384 }
385
386 const float* logitsPtr = outTensor.GetTensorData<float>();
387 if (!logitsPtr) {
388 ATH_MSG_ERROR("Failed to get logits data pointer");
389 return StatusCode::FAILURE;
390 }
391
392 // Map logits to buckets - only considering non-empty buckets
393 size_t totalNumPredictions = static_cast<size_t>(outShape[0]);
394 size_t predIdx = 0;
395
396 for (const MuonR4::SpacePointBucket* bucket : *spContainer) {
397 if (!bucket || bucket->empty()) {
398 continue; // Skip empty buckets
399 }
400
401 if (predIdx >= totalNumPredictions) {
402 ATH_MSG_ERROR("More non-empty buckets than predictions from model");
403 return StatusCode::FAILURE;
404 }
405
406 std::vector<float> scores(3);
407 scores[0] = logitsPtr[3 * predIdx + 0];
408 scores[1] = logitsPtr[3 * predIdx + 1];
409 scores[2] = logitsPtr[3 * predIdx + 2];
410
411 bucketScoreMap[bucket] = scores;
412 ++predIdx;
413 }
414
415 if (predIdx != totalNumPredictions) {
416 ATH_MSG_ERROR("Number of non-empty buckets (" << predIdx << ") does not match predictions ("
417 << totalNumPredictions << ")");
418 return StatusCode::FAILURE;
419 }
420
421 ATH_MSG_DEBUG("Successfully computed ML bucket scores for " << bucketScoreMap.size() << " buckets");
422 return StatusCode::SUCCESS;
423 }
#define ATH_CHECK
Evaluate an expression and check for errors.
#define ATH_MSG_ERROR(x)
#define ATH_MSG_DEBUG(x)
ToolHandle< MuonML::IGraphInferenceTool > m_inferenceTool
Inference tool for ML bucket scoring (optional).
Gaudi::Property< bool > m_doMLBucketScore
std::unique_ptr< InferenceGraph > graph
Pointer to the graph to be parsed to ONNX.
Definition GraphData.h:46

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

◆ dumpContainer()

StatusCode MuonR4::BucketDumperAlg::dumpContainer ( const EventContext & ctx,
const SG::ReadHandleKey< SpacePointContainer > & spacePointKey,
const SG::ReadHandleKey< xAOD::MuonSegmentContainer > & segmentKey )
private

Dumps the space point container with the associated muon segment container.

Parameters
spacePointKeyKey to the space point key (Legacy /Nsw)
segmentKeyKey to the fitted segments (Legacy / Nsw)

Filter random noise

Bucket identifier

Global bucket position

Flag whether the bucket contains at least one good hit

Number of reconstructed segments in the bucket

Associate the Space points to the reconstructed segments

For the moment this only works on MC

Truth segment parameters

Calculate the layer of the space point

Definition at line 95 of file BucketDumperAlg.cxx.

97 {
98
99 // Define the segment mapping type
100 using SegmentsPerBucket_t = std::unordered_map<const SpacePointBucket*,
101 std::set<const xAOD::MuonSegment*, LocalSegSorter>>;
102 SegmentsPerBucket_t segmentMap{};
103
104 const xAOD::MuonSegmentContainer* readSegment{nullptr};
105 ATH_CHECK(SG::get(readSegment, segmentKey, ctx));
106 if (readSegment) {
107 for (const xAOD::MuonSegment* segment : *readSegment) {
108 segmentMap[detailedSegment(*segment)->parent()->parentBucket()].insert(segment);
109 }
110 }
111
112 const ActsTrk::GeometryContext* gctx{nullptr};
113 ATH_CHECK(SG::get(gctx, m_geoCtxKey, ctx));
114
115 const SpacePointContainer* spContainer{nullptr};
116 ATH_CHECK(SG::get(spContainer, spacePointKey, ctx));
117
118 // Compute ML bucket scores if enabled - only for MuonSpacePoints container
119 // The tool is configured to read "MuonSpacePoints" from StoreGate, so we only run it for that container
120 std::unordered_map<const SpacePointBucket*, std::vector<float>> bucketScores;
121 const std::string& containerName = spacePointKey.key();
122 if (m_doMLBucketScore && containerName == "MuonSpacePoints") {
123 ATH_MSG_DEBUG("Computing ML bucket scores for container: " << containerName);
124 ATH_CHECK(computeAllBucketScores(ctx, spContainer, bucketScores));
125 ATH_MSG_DEBUG("ML scoring completed, got scores for " << bucketScores.size() << " buckets");
126 } else if (m_doMLBucketScore && containerName != "MuonSpacePoints") {
127 ATH_MSG_DEBUG("Skipping ML scoring for container: " << containerName
128 << " (only MuonSpacePoints is supported)");
129 }
130
131 CLHEP::HepRandomEngine* rndEngine = getRandomEngine(ctx);
132
133 const SpacePointPerLayerSorter layerSorter{};
134
135 for(const SpacePointBucket* bucket : *spContainer) {
137 if (!m_isMC && segmentMap[bucket].empty() && m_fracToKeep < 1. &&
138 CLHEP::RandFlat::shoot(rndEngine,0.,1.) > m_fracToKeep) {
139 ATH_MSG_VERBOSE("Skipping bucket without segment");
140 continue;
141 }
142
143 if (m_doMLBucketFilter && !bucketScores.empty()) {
144 auto scoreIt = bucketScores.find(bucket);
145 if (scoreIt != bucketScores.end() && scoreIt->second.size() >= 3) {
146 const auto& logits = scoreIt->second;
147 int predictedClass = 0;
148 float maxLogit = logits[0];
149 if (logits[1] > maxLogit) { maxLogit = logits[1]; predictedClass = 1; }
150 if (logits[2] > maxLogit) { maxLogit = logits[2]; predictedClass = 2; }
151 if (predictedClass == 0) continue;
152 }
153 }
154
156 m_bucket_sector = bucket->msSector()->sector();
157 m_bucket_chamberIdx = static_cast<uint8_t>(bucket->msSector()->chamberIndex());
158 m_bucket_side = bucket->msSector()->side();
160 m_bucket_min = bucket->coveredMin();
161 m_bucket_max = bucket->coveredMax();
162
164 const Amg::Vector3D bucketPos = bucket->msSector()->localToGlobalTransform(*gctx) *
165 (0.5*(bucket->coveredMin() + bucket->coveredMax()) * Amg::Vector3D::UnitY());
166 m_bucket_posX = bucketPos.x();
167 m_bucket_posY = bucketPos.y();
168 m_bucket_posZ = bucketPos.z();
169
171 m_bucket_truthHit = std::ranges::any_of(*bucket,[this](const SpacePointBucket::value_type & sp){
172 return m_visionTool->isLabeled(*sp);
173 });
174
176 m_bucket_segments = segmentMap[bucket].size();
177
178
179
180 std::unordered_map<const SpacePoint*, std::vector<int16_t>> spacePointToSegment{};
181 std::set<const xAOD::MuonSegment*, LocalSegSorter> truthSegments{};
183 auto match_itr = segmentMap.find(bucket);
184 if (match_itr != segmentMap.end()) {
185 for (const xAOD::MuonSegment* segment : match_itr->second) {
186 for (const auto& meas : detailedSegment(*segment)->measurements()) {
187 if (meas->fitState() == CalibratedSpacePoint::State::Valid) {
188 spacePointToSegment[meas->spacePoint()].push_back(segment->index());
189 }
190 }
191 using namespace SegmentFit;
192 const auto pars = localSegmentPars(*segment);
193 m_segmentLocX += pars[Acts::toUnderlying(ParamDefs::x0)];
194 m_segmentLocY += pars[Acts::toUnderlying(ParamDefs::y0)];
195 m_segmentLocTheta += pars[Acts::toUnderlying(ParamDefs::theta)];
196 m_segmentLocPhi += pars[Acts::toUnderlying(ParamDefs::phi)];
197
199 unsigned truthLink = -1;
200 if (const xAOD::TruthParticle* truthPart = getTruthMatchedParticle(*segment)) {
201 truthLink = truthPart->index();
202 }
204 if (const xAOD::MuonSegment* truthSeg = getMatchedTruthSegment(*segment)) {
205 truthSegments.insert(truthSeg);
206 }
207 m_segmentTruthIdx+=truthLink;
208 m_segmentPos.push_back(segment->position());
209 m_segmentDir.push_back(segment->direction());
210 m_segment_chiSquared.push_back(segment->chiSquared());
211 m_segment_numberDoF.push_back(segment->numberDoF());
212 }
213 }
214
215 std::vector<unsigned int> layNumbers{};
216 std::unordered_map<const SpacePoint*, std::vector<const xAOD::MuonSegment*>> spToTrueSeg{};
217 if (m_isMC) {
218 using SegLinkVec_t = std::vector<ElementLink<xAOD::MuonSegmentContainer>>;
219 static const SG::ConstAccessor<SegLinkVec_t> segAcc{"truthSegmentLinks"};
220 for (const auto& sp : *bucket){
221 for (const auto& link : segAcc(*sp->primaryMeasurement())) {
222 spToTrueSeg[sp.get()].push_back(*link);
223 truthSegments.insert(*link);
224 }
225 }
226 }
227
228 for(const SpacePointBucket::value_type& sp : *bucket) {
230 const unsigned int layNum = layerSorter.sectorLayerNum(*sp);
231 if (std::find(layNumbers.begin(), layNumbers.end(), layNum) == layNumbers.end()) {
232 layNumbers.push_back(layNum);
233 }
234 const unsigned layer = layNumbers.size()-1;
235
236 const Identifier& id = sp->identify();
237
239 const auto* dc = static_cast<const xAOD::MdtDriftCircle*>(sp->primaryMeasurement());
240 if (dc->status() != Muon::MdtDriftCircleStatus::MdtStatusDriftTime){
241 continue;
242 }
243 m_spoint_isMdt.push_back(true);
244 m_spoint_isStrip.push_back(false);
245 m_spoint_adc.push_back(dc->adc());
246 m_spoint_tdc.push_back(dc->tdc());
247 } else {
248 // check the technology to fill channel, adc and tdc... pushing back 0 for now
249 m_spoint_adc.push_back(0);
250 m_spoint_tdc.push_back(0);
251 m_spoint_isMdt.push_back(false);
252 m_spoint_isStrip.push_back(true);
253 }
254
255 m_spoint_id.push_back(id);
256
257
258 const std::vector<int16_t>& segIdxs = spacePointToSegment[sp.get()];
259 m_spoint_mat[m_spoint_mat.size()] = segIdxs;
260 auto& trueSegLinks = m_spoint_trueSeg[m_spoint_trueSeg.size()];
261 for (const xAOD::MuonSegment* matchedSeg : spToTrueSeg[sp.get()]) {
262 trueSegLinks.push_back(std::distance(truthSegments.begin(), truthSegments.find(matchedSeg)));
263 }
264 m_spoint_nSegments.push_back(segIdxs.size());
265
266 m_bucket_spacePoints = bucket->size();
267 m_spoint_localPosition.push_back(sp->localPosition());
268 using CovIdx = SpacePoint::CovIdx;
269 m_spoint_covX.push_back(sp->covariance()[Acts::toUnderlying(CovIdx::phiCov)]);
270 m_spoint_covY.push_back(sp->covariance()[Acts::toUnderlying(CovIdx::etaCov)]);
271 m_spoint_driftR.push_back(sp->driftRadius());
272 m_spoint_measuresEta.push_back(sp->measuresEta());
273 m_spoint_measuresPhi.push_back(sp->measuresPhi());
274 m_spoint_nEtaInstances.push_back(sp->nEtaInstanceCounts());
275 m_spoint_nPhiInstances.push_back(sp->nPhiInstanceCounts());
276 m_spoint_dimension.push_back(sp->dimension());
277 m_spoint_layer.push_back(layer);
278 if (m_visionTool.isEnabled()) {
279 m_spoint_trueLabel.push_back(m_visionTool->isLabeled(*sp));
280 }
281
282 Amg::Vector3D globalPos = sp->msSector()->localToGlobalTransform(*gctx) * sp->localPosition();
283 m_spoint_globalPosition.push_back( globalPos );
284 }
285
286 for (const xAOD::MuonSegment* truthSeg: truthSegments) {
287 using namespace SegmentFit;
288 const auto truthPars = localSegmentPars(*truthSeg);
289 m_truthSegLocX += truthPars[Acts::toUnderlying(ParamDefs::x0)];
290 m_truthSegLocY += truthPars[Acts::toUnderlying(ParamDefs::y0)];
291 m_truthSegLocTheta += truthPars[Acts::toUnderlying(ParamDefs::theta)];
292 m_truthSegLocPhi += truthPars[Acts::toUnderlying(ParamDefs::phi)];
293 }
294
295 m_bucket_layers = layNumbers.size();
296
297 // Add ML bucket filter scores if available
298 // Only add scores if we actually computed them for this container
299 if (m_doMLBucketScore && !bucketScores.empty()) {
300 auto scoreIt = bucketScores.find(bucket);
301 if (scoreIt != bucketScores.end() && scoreIt->second.size() >= 3) {
302 m_bucket_ml_score_class0.push_back(scoreIt->second[0]);
303 m_bucket_ml_score_class1.push_back(scoreIt->second[1]);
304 m_bucket_ml_score_class2.push_back(scoreIt->second[2]);
305 } else {
306 // Empty buckets are skipped during inference, so missing scores is expected
307 // Only warn if a non-empty bucket is missing scores (indicates a real problem)
308 if (bucket && !bucket->empty()) {
309 ATH_MSG_WARNING("Non-empty bucket from scored container not found in ML scores map");
310 }
311 m_bucket_ml_score_class0.push_back(0.f);
312 m_bucket_ml_score_class1.push_back(0.f);
313 m_bucket_ml_score_class2.push_back(0.f);
314 }
315 } else if (m_doMLBucketScore) {
316 // No scores computed for this container (e.g., NswSpacePoints) - fill with zeros
317 m_bucket_ml_score_class0.push_back(0.f);
318 m_bucket_ml_score_class1.push_back(0.f);
319 m_bucket_ml_score_class2.push_back(0.f);
320 }
321
322 if (!m_tree.fill(ctx)) {
323 return StatusCode::FAILURE;
324 }
325 }
326
327 return StatusCode::SUCCESS;
328
329 }
#define ATH_MSG_VERBOSE(x)
#define ATH_MSG_WARNING(x)
static Double_t sp
static const Attributes_t empty
ToolHandle< MuonValR4::IPatternVisualizationTool > m_visionTool
Pattern visualization tool.
MuonVal::ScalarBranch< Char_t > & m_bucket_side
MuonVal::ThreeVectorBranch m_spoint_localPosition
MuonVal::VectorBranch< unsigned short > & m_spoint_measuresEta
MuonVal::ScalarBranch< uint8_t > & m_bucket_sector
MuonVal::VectorBranch< float > & m_truthSegLocTheta
MuonVal::VectorBranch< float > & m_segment_numberDoF
Gaudi::Property< bool > m_doMLBucketFilter
MuonVal::ScalarBranch< float > & m_bucket_min
MuonVal::VectorBranch< uint16_t > & m_spoint_layer
MuonVal::ThreeVectorBranch m_segmentPos
MuonVal::VectorBranch< float > & m_spoint_covX
MuonVal::VectorBranch< float > & m_segmentLocTheta
MuonVal::VectorBranch< float > & m_truthSegLocY
MuonVal::VectorBranch< float > & m_segment_chiSquared
CLHEP::HepRandomEngine * getRandomEngine(const EventContext &ctx) const
MuonVal::MatrixBranch< int16_t > & m_spoint_mat
MuonVal::ScalarBranch< float > & m_bucket_posZ
MuonVal::ScalarBranch< uint16_t > & m_bucket_layers
MuonVal::VectorBranch< float > & m_truthSegLocPhi
ActsTrk::GeoContextReadKey_t m_geoCtxKey
MuonVal::VectorBranch< float > & m_segmentLocY
MuonVal::VectorBranch< uint16_t > & m_segmentTruthIdx
MuonVal::VectorBranch< float > & m_bucket_ml_score_class1
MuonVal::MuonTesterTree m_tree
MuonVal::VectorBranch< float > & m_segmentLocX
Gaudi::Property< double > m_fracToKeep
MuonVal::ScalarBranch< float > & m_bucket_max
MuonVal::VectorBranch< uint16_t > & m_spoint_nSegments
Gaudi::Property< bool > m_isMC
MuonVal::VectorBranch< float > & m_bucket_ml_score_class0
ML bucket filter scores (3 classes for the filter model).
MuonVal::VectorBranch< unsigned int > & m_spoint_nPhiInstances
StatusCode computeAllBucketScores(const EventContext &ctx, const SpacePointContainer *spContainer, std::unordered_map< const SpacePointBucket *, std::vector< float > > &bucketScoreMap) const
Computes ML bucket scores using the inference tool for a specific container.
MuonVal::ScalarBranch< uint8_t > & m_bucket_truthHit
MuonVal::ThreeVectorBranch m_segmentDir
MuonVal::VectorBranch< uint16_t > & m_spoint_tdc
MuonVal::VectorBranch< float > & m_spoint_driftR
MuonVal::ScalarBranch< float > & m_bucket_posY
MuonVal::VectorBranch< float > & m_segmentLocPhi
MuonVal::ScalarBranch< float > & m_bucket_posX
MuonVal::VectorBranch< uint16_t > & m_spoint_adc
MuonVal::MuonIdentifierBranch m_spoint_id
MuonVal::ScalarBranch< uint16_t > & m_bucket_segments
MuonVal::VectorBranch< unsigned int > & m_spoint_dimension
MuonVal::ScalarBranch< uint8_t > & m_bucket_chamberIdx
MuonVal::VectorBranch< float > & m_spoint_covY
MuonVal::MatrixBranch< int16_t > & m_spoint_trueSeg
MuonVal::VectorBranch< unsigned short > & m_spoint_isStrip
MuonVal::VectorBranch< unsigned int > & m_spoint_nEtaInstances
MuonVal::VectorBranch< float > & m_truthSegLocX
MuonVal::VectorBranch< unsigned short > & m_spoint_trueLabel
MuonVal::ScalarBranch< uint16_t > & m_bucket_spacePoints
MuonVal::VectorBranch< float > & m_bucket_ml_score_class2
MuonVal::VectorBranch< unsigned short > & m_spoint_isMdt
MuonVal::ThreeVectorBranch m_spoint_globalPosition
MuonVal::VectorBranch< unsigned short > & m_spoint_measuresPhi
const SpacePointBucket * parentBucket() const
Returns the bucket out of which the seed was formed.
const SegmentSeed * parent() const
Returns the seed out of which the segment was built.
const MeasVec & measurements() const
Returns the associated measurements.
const std::string & key() const
Return the StoreGate ID for the referenced object.
float numberDoF() const
Returns the numberDoF.
Amg::Vector3D direction() const
Returns the direction as Amg::Vector.
float chiSquared() const
Amg::Vector3D position() const
Returns the position as Amg::Vector.
Eigen::Matrix< double, 3, 1 > Vector3D
Parameters localSegmentPars(const xAOD::MuonSegment &seg)
Returns the localSegPars decoration from a xAODMuon::Segment.
const xAOD::TruthParticle * getTruthMatchedParticle(const xAOD::MuonSegment &segment)
Returns the particle truth-matched to the segment.
const xAOD::MuonSegment * getMatchedTruthSegment(const xAOD::MuonSegment &segment)
Returns the truth-matched segment.
std::vector< SegLink_t > SegLinkVec_t
DataVector< SpacePointBucket > SpacePointContainer
Abrivation of the space point container type.
const Segment * detailedSegment(const xAOD::MuonSegment &seg)
Helper function to navigate from the xAOD::MuonSegment to the MuonR4::Segment.
SG::ConstAccessor< T, ALLOC > ConstAccessor
Definition AuxElement.h:570
const T * get(const ReadCondHandleKey< T > &key, const EventContext &ctx)
Convenience function to retrieve an object given a ReadCondHandleKey.
@ layer
Definition HitInfo.h:79
MdtDriftCircle_v1 MdtDriftCircle
MuonSegmentContainer_v1 MuonSegmentContainer
Definition of the current "MuonSegment container version".
TruthParticle_v1 TruthParticle
Typedef to implementation.
MuonSegment_v1 MuonSegment
Reference the current persistent version:

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

212{
213 // Build a 32 bit hash out of the name
214 const hash_t effHash = this->hash(effName);
215
216 // See if this entry exists in the map
217 EffMap_t::const_iterator it = m_effMap.find( effHash );
218 if ( it == m_effMap.end() ) // It doesn't exist!
219 { // Let's see into the THistSvc if somebody else has registered the TEfficiency...
220
221 // Need to copy the strings as we will massage them from here on
222 std::string effNameCopy = effName;
223 std::string tDirCopy = tDir;
224 std::string streamCopy = stream;
225
226 // Massage the final string to book things
227 std::string bookingString("");
228 this->buildBookingString( bookingString, effNameCopy, tDirCopy, streamCopy ,false);
229
230 TEfficiency* effPointer(NULL);
231 if ( !((histSvc()->getEfficiency(bookingString, effPointer)).isSuccess()) )
232 {
233 // Massage the final string to book things
234 std::string bookingString("");
235 this->buildBookingString( bookingString, effNameCopy, tDirCopy, streamCopy, true );
236
237 if ( !((histSvc()->getEfficiency(bookingString, effPointer)).isSuccess()) )
238 {
239 m_msg << MSG::WARNING
240 << "Problem retrieving the TEfficiency with name (including pre- and post-fixes) "
241 << m_histNamePrefix + effNameCopy + m_histNamePostfix
242 << " or with name " << effNameCopy
243 << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
244 << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
245 return NULL;
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 // If we get to here, we actually found the TEfficiency in the THistSvc.
253 // So let's add it to the local cache map and return its pointer
254 m_effMap.insert( m_effMap.end(), std::pair< const hash_t, TEfficiency* >( effHash, effPointer ) );
255 return effPointer;
256 }
257
258 // Return the pointer to the TEfficiency that we got from the local cache map
259 return it->second;
260}
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 MuonR4::BucketDumperAlg::execute ( const EventContext & ctx)
finaloverridevirtual

Execute method.

Implements AthAlgorithm.

Definition at line 85 of file BucketDumperAlg.cxx.

85 {
86 SG::ReadHandleKey<xAOD::MuonSegmentContainer> emptyKey{};
88
89 for (unsigned keyNum = 0 ; keyNum < m_spacePointKeys.size(); ++keyNum) {
91 keyNum < m_inSegmentKeys.size() ? m_inSegmentKeys[keyNum] : emptyKey));
92 }
93 return StatusCode::SUCCESS;
94 }
SG::ReadHandleKeyArray< SpacePointContainer > m_spacePointKeys
StatusCode dumpContainer(const EventContext &ctx, const SG::ReadHandleKey< SpacePointContainer > &spacePointKey, const SG::ReadHandleKey< xAOD::MuonSegmentContainer > &segmentKey)
Dumps the space point container with the associated muon segment container.
SG::ReadHandleKeyArray< xAOD::MuonSegmentContainer > m_inSegmentKeys
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.

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

◆ finalize()

StatusCode MuonR4::BucketDumperAlg::finalize ( )
finaloverridevirtual

Definition at line 80 of file BucketDumperAlg.cxx.

80 {
81 ATH_CHECK(m_tree.write());
82 return StatusCode::SUCCESS;
83 }

◆ getRandomEngine()

CLHEP::HepRandomEngine * MuonR4::BucketDumperAlg::getRandomEngine ( const EventContext & ctx) const
private

Definition at line 331 of file BucketDumperAlg.cxx.

331 {
332 ATHRNG::RNGWrapper* rngWrapper = m_rndmSvc->getEngine(this, m_streamName);
333 std::string rngName = m_streamName;
334 rngWrapper->setSeed(rngName, ctx);
335 return rngWrapper->getEngine(ctx);
336 }
void setSeed(const std::string &algName, const EventContext &ctx)
Set the random seed using a string (e.g.
Definition RNGWrapper.h:154
CLHEP::HepRandomEngine * getEngine(const EventContext &ctx) const
Retrieve the random engine corresponding to the provided EventContext.
Definition RNGWrapper.h:108
Gaudi::Property< std::string > m_streamName
ServiceHandle< IAthRNGSvc > m_rndmSvc

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

454{
455 // Build a 32 bit hash out of the name
456 const hash_t graphHash = this->hash(graphName);
457
458 // See if this entry exists in the map
459 GraphMap_t::const_iterator it = m_graphMap.find( graphHash );
460 if ( it == m_graphMap.end() ) // It doesn't exist!
461 { // Let's see into the THistSvc if somebody else has registered the TGraph...
462
463 // Need to copy the strings as we will massage them from here on
464 std::string graphNameCopy = graphName;
465 std::string tDirCopy = tDir;
466 std::string streamCopy = stream;
467
468 // Massage the final string to book things
469 std::string bookingString("");
470 this->buildBookingString( bookingString, graphNameCopy, tDirCopy, streamCopy, false);
471
472 TGraph* graphPointer(nullptr);
473 if ( !((histSvc()->getGraph(bookingString, graphPointer)).isSuccess()) )
474 {
475 // Massage the final string to book things
476 std::string bookingString("");
477 this->buildBookingString( bookingString, graphNameCopy, tDirCopy, streamCopy, true );
478
479 if ( !((histSvc()->getGraph(bookingString, graphPointer)).isSuccess()) )
480 {
481 m_msg << MSG::WARNING
482 << "Problem retrieving the TGraph with name (including pre- and post-fixes) "
483 << m_histNamePrefix + graphNameCopy + m_histNamePostfix
484 << " or with name " << graphNameCopy
485 << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
486 << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
487 return nullptr;
488 }
489 // If we get to here, we actually found the TGraph in the THistSvc.
490 // So let's add it to the local cache map and return its pointer
491 m_graphMap.insert( m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphPointer ) );
492 return graphPointer;
493 }
494 // If we get to here, we actually found the TGraph in the THistSvc.
495 // So let's add it to the local cache map and return its pointer
496 m_graphMap.insert( m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphPointer ) );
497 return graphPointer;
498 }
499
500
501 // Return the pointer to the TGraph that we got from the local cache map
502 return it->second;
503}
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 164 of file AthHistogramming.cxx.

167{
168 std::string histNameCopy = histName;
169 std::string tDirCopy = tDir;
170 std::string streamCopy = stream;
171
172 std::string bookingString;
173 this->buildBookingString(bookingString, histNameCopy, tDirCopy, streamCopy, false);
174
175 const hash_t histHash = this->hash(histNameCopy);
176
177 HistMap_t::const_iterator it = m_histMap.find(histHash);
178 if (it == m_histMap.end()) {
179 TH1* histPointer(nullptr);
180
181 if (!histSvc()->getHist(bookingString, histPointer).isSuccess()) {
182 std::string prefixedHistNameCopy = histName;
183 std::string prefixedTDirCopy = tDir;
184 std::string prefixedStreamCopy = stream;
185
186 std::string prefixedBookingString;
187 this->buildBookingString(prefixedBookingString,
188 prefixedHistNameCopy,
189 prefixedTDirCopy,
190 prefixedStreamCopy,
191 true);
192
193 if (!histSvc()->getHist(prefixedBookingString, histPointer).isSuccess()) {
194 m_msg << MSG::WARNING
195 << "Problem retrieving the histogram with name (including pre- and post-fixes) "
196 << m_histNamePrefix + histNameCopy + m_histNamePostfix
197 << " or with name " << histNameCopy
198 << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
199 << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
200 return nullptr;
201 }
202 }
203
204 m_histMap.insert(m_histMap.end(), std::pair<const hash_t, TH1*>(histHash, histPointer));
205 return histPointer;
206 }
207
208 return it->second;
209}
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}

◆ initialize()

StatusCode MuonR4::BucketDumperAlg::initialize ( )
finaloverridevirtual

Definition at line 45 of file BucketDumperAlg.cxx.

45 {
46 ATH_CHECK(m_spacePointKeys.initialize());
47 ATH_CHECK(m_inSegmentKeys.initialize());
48 if (m_isMC) {
49 for (const auto& key : m_inSegmentKeys) {
50 m_truthDecorKeys.emplace_back(key, "truthParticleLink");
51 }
52 }
53 ATH_CHECK(m_truthDecorKeys.initialize());
54 ATH_CHECK(m_geoCtxKey.initialize());
55 m_tree.addBranch(std::make_shared<MuonVal::EventHashBranch>(m_tree.tree()));
56 ATH_CHECK(m_visionTool.retrieve(EnableTool{!m_visionTool.empty()}));
57 if (m_visionTool.empty()) {
58 m_tree.disableBranch(m_spoint_trueLabel.name());
59 }
60
61 // Setup ML bucket score branches
63 ATH_CHECK(m_inferenceTool.retrieve());
64 ATH_MSG_INFO("ML bucket scoring enabled with tool: " << m_inferenceTool.name());
65 }
66
67 // Disable ML score branches if not enabled
68 if (!m_doMLBucketScore) {
69 m_tree.disableBranch(m_bucket_ml_score_class0.name());
70 m_tree.disableBranch(m_bucket_ml_score_class1.name());
71 m_tree.disableBranch(m_bucket_ml_score_class2.name());
72 }
73
74 ATH_CHECK(m_tree.init(this));
75 ATH_MSG_DEBUG("Successfully initialized");
76
77 return StatusCode::SUCCESS;
78 }
#define ATH_MSG_INFO(x)
SG::ReadDecorHandleKeyArray< xAOD::MuonSegmentContainer > m_truthDecorKeys

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

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

554{
555 size_t pos = 0;
556 while((pos = str.find(oldStr, pos)) != std::string::npos)
557 {
558 str.replace(pos, oldStr.length(), newStr);
559 pos += newStr.length();
560 }
561}

◆ 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 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< Gaudi::Algorithm > >::renounceArray ( SG::VarHandleKeyArray & handlesArray)
inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364 {
366 }

◆ setFilterPassed()

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

Set filter decision:

Reimplemented in AthFilterAlgorithm.

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
const ServiceHandle< ITHistSvc > & histSvc() const
The standard THistSvc (for writing histograms and TTrees and more to a root file) Returns (kind of) a...
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< 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 ( const std::string & treeName,
const std::string & tDir = "",
const std::string & stream = "" )
protectedinherited

Simplify the retrieval of registered TTrees.

Definition at line 339 of file AthHistogramming.cxx.

340{
341 // Build a 32 bit hash out of the name
342 const hash_t treeHash = this->hash(treeName);
343
344 // See if this entry exists in the map
345 TreeMap_t::const_iterator it = m_treeMap.find( treeHash );
346 if ( it == m_treeMap.end() ) // It doesn't exist!
347 { // Let's see into the THistSvc if somebody else has registered the TTree...
348
349 // Need to copy the strings as we will massage them from here on
350 std::string treeNameCopy = treeName;
351 std::string tDirCopy = tDir;
352 std::string streamCopy = stream;
353
354 // Massage the final string to book things
355 std::string bookingString("");
356 this->buildBookingString( bookingString, treeNameCopy, tDirCopy, streamCopy );
357
358 TTree* treePointer(NULL);
359 if ( !((histSvc()->getTree(bookingString, treePointer)).isSuccess()) )
360 {
361 m_msg << MSG::WARNING
362 << "Problem retrieving the TTree with name " << treeNameCopy
363 << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
364 << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
365 return NULL;
366 }
367 // If we get to here, we actually found the TTree in the THistSvc.
368 // So let's add it to the local cache map and return its pointer
369 m_treeMap.insert( m_treeMap.end(), std::pair< const hash_t, TTree* >( treeHash, treePointer ) );
370 return treePointer;
371 }
372
373 // Return the pointer to the TTree that we got from the local cache map
374 return it->second;
375}
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 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 }

Member Data Documentation

◆ m_bucket_chamberIdx

MuonVal::ScalarBranch<uint8_t>& MuonR4::BucketDumperAlg::m_bucket_chamberIdx {m_tree.newScalar<uint8_t>("bucket_chamberIndex")}
private

Definition at line 96 of file BucketDumperAlg.h.

96{m_tree.newScalar<uint8_t>("bucket_chamberIndex")};

◆ m_bucket_layers

MuonVal::ScalarBranch<uint16_t>& MuonR4::BucketDumperAlg::m_bucket_layers {m_tree.newScalar<uint16_t>("bucket_layers", 0)}
private

Definition at line 89 of file BucketDumperAlg.h.

89{m_tree.newScalar<uint16_t>("bucket_layers", 0)};
setWord1 uint16_t

◆ m_bucket_max

MuonVal::ScalarBranch<float>& MuonR4::BucketDumperAlg::m_bucket_max {m_tree.newScalar<float>("bucket_max", -1)}
private

Definition at line 86 of file BucketDumperAlg.h.

86{m_tree.newScalar<float>("bucket_max", -1)};

◆ m_bucket_min

MuonVal::ScalarBranch<float>& MuonR4::BucketDumperAlg::m_bucket_min {m_tree.newScalar<float>("bucket_min", -1)}
private

Definition at line 85 of file BucketDumperAlg.h.

85{m_tree.newScalar<float>("bucket_min", -1)};

◆ m_bucket_ml_score_class0

MuonVal::VectorBranch<float>& MuonR4::BucketDumperAlg::m_bucket_ml_score_class0 {m_tree.newVector<float>("bucket_ml_score_class0", false)}
private

ML bucket filter scores (3 classes for the filter model).

Definition at line 100 of file BucketDumperAlg.h.

100{m_tree.newVector<float>("bucket_ml_score_class0", false)};

◆ m_bucket_ml_score_class1

MuonVal::VectorBranch<float>& MuonR4::BucketDumperAlg::m_bucket_ml_score_class1 {m_tree.newVector<float>("bucket_ml_score_class1", false)}
private

Definition at line 101 of file BucketDumperAlg.h.

101{m_tree.newVector<float>("bucket_ml_score_class1", false)};

◆ m_bucket_ml_score_class2

MuonVal::VectorBranch<float>& MuonR4::BucketDumperAlg::m_bucket_ml_score_class2 {m_tree.newVector<float>("bucket_ml_score_class2", false)}
private

Definition at line 102 of file BucketDumperAlg.h.

102{m_tree.newVector<float>("bucket_ml_score_class2", false)};

◆ m_bucket_posX

MuonVal::ScalarBranch<float>& MuonR4::BucketDumperAlg::m_bucket_posX {m_tree.newScalar<float>("bucket_positionX")}
private

Definition at line 91 of file BucketDumperAlg.h.

91{m_tree.newScalar<float>("bucket_positionX")};

◆ m_bucket_posY

MuonVal::ScalarBranch<float>& MuonR4::BucketDumperAlg::m_bucket_posY {m_tree.newScalar<float>("bucket_positionY")}
private

Definition at line 92 of file BucketDumperAlg.h.

92{m_tree.newScalar<float>("bucket_positionY")};

◆ m_bucket_posZ

MuonVal::ScalarBranch<float>& MuonR4::BucketDumperAlg::m_bucket_posZ {m_tree.newScalar<float>("bucket_positionZ")}
private

Definition at line 93 of file BucketDumperAlg.h.

93{m_tree.newScalar<float>("bucket_positionZ")};

◆ m_bucket_sector

MuonVal::ScalarBranch<uint8_t>& MuonR4::BucketDumperAlg::m_bucket_sector {m_tree.newScalar<uint8_t>("bucket_sector")}
private

Definition at line 95 of file BucketDumperAlg.h.

95{m_tree.newScalar<uint8_t>("bucket_sector")};

◆ m_bucket_segments

MuonVal::ScalarBranch<uint16_t>& MuonR4::BucketDumperAlg::m_bucket_segments {m_tree.newScalar<uint16_t>("bucket_segments", 0)}
private

Definition at line 88 of file BucketDumperAlg.h.

88{m_tree.newScalar<uint16_t>("bucket_segments", 0)};

◆ m_bucket_side

MuonVal::ScalarBranch<Char_t>& MuonR4::BucketDumperAlg::m_bucket_side {m_tree.newScalar<Char_t>("bucket_side")}
private

Definition at line 97 of file BucketDumperAlg.h.

97{m_tree.newScalar<Char_t>("bucket_side")};

◆ m_bucket_spacePoints

MuonVal::ScalarBranch<uint16_t>& MuonR4::BucketDumperAlg::m_bucket_spacePoints {m_tree.newScalar<uint16_t>("bucket_spacePoints", 0)}
private

Definition at line 87 of file BucketDumperAlg.h.

87{m_tree.newScalar<uint16_t>("bucket_spacePoints", 0)};

◆ m_bucket_truthHit

MuonVal::ScalarBranch<uint8_t>& MuonR4::BucketDumperAlg::m_bucket_truthHit {m_tree.newScalar<uint8_t>("bucket_hasTruth", 0)}
private

Definition at line 90 of file BucketDumperAlg.h.

90{m_tree.newScalar<uint8_t>("bucket_hasTruth", 0)};

◆ m_detStore

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

Pointer to StoreGate (detector store by default).

Definition at line 393 of file AthCommonDataStore.h.

◆ m_doMLBucketFilter

Gaudi::Property<bool> MuonR4::BucketDumperAlg::m_doMLBucketFilter
private
Initial value:
{this, "DoMLBucketFilter", false,
"Skip dumping buckets classified as class 0 (reject) by ML model. Requires DoMLBucketScore=true."}

Definition at line 71 of file BucketDumperAlg.h.

71 {this, "DoMLBucketFilter", false,
72 "Skip dumping buckets classified as class 0 (reject) by ML model. Requires DoMLBucketScore=true."};

◆ m_doMLBucketScore

Gaudi::Property<bool> MuonR4::BucketDumperAlg::m_doMLBucketScore
private
Initial value:
{this, "DoMLBucketScore", false,
"Enable ML bucket score computation and dumping (requires InferenceTool to be configured)"}

Definition at line 69 of file BucketDumperAlg.h.

69 {this, "DoMLBucketScore", false,
70 "Enable ML bucket score computation and dumping (requires InferenceTool to be configured)"};

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

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

Pointer to StoreGate (event store by default).

Definition at line 390 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_fracToKeep

Gaudi::Property<double> MuonR4::BucketDumperAlg::m_fracToKeep {this,"dataFracToKeep", 1.}
private

Definition at line 67 of file BucketDumperAlg.h.

67{this,"dataFracToKeep", 1.}; // 0.055 to balanced dataset without MC

◆ m_geoCtxKey

ActsTrk::GeoContextReadKey_t MuonR4::BucketDumperAlg::m_geoCtxKey {this, "AlignmentKey", "ActsAlignment", "cond handle key"}
private

Definition at line 64 of file BucketDumperAlg.h.

64{this, "AlignmentKey", "ActsAlignment", "cond handle key"};

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

ToolHandle<MuonML::IGraphInferenceTool> MuonR4::BucketDumperAlg::m_inferenceTool
private
Initial value:
{this, "InferenceTool", "",
"Optional ML inference tool for computing bucket scores. If not provided, scores will be empty."}

Inference tool for ML bucket scoring (optional).

Definition at line 76 of file BucketDumperAlg.h.

76 {this, "InferenceTool", "",
77 "Optional ML inference tool for computing bucket scores. If not provided, scores will be empty."};

◆ m_inSegmentKeys

SG::ReadHandleKeyArray<xAOD::MuonSegmentContainer> MuonR4::BucketDumperAlg::m_inSegmentKeys {this, "SegmentKey", {"MuonSegmentsFromR4"}}
private

Definition at line 61 of file BucketDumperAlg.h.

61{this, "SegmentKey", {"MuonSegmentsFromR4"}};

◆ m_isMC

Gaudi::Property<bool> MuonR4::BucketDumperAlg::m_isMC {this, "isMC", true}
private

Definition at line 66 of file BucketDumperAlg.h.

66{this, "isMC", true};

◆ m_msg

MsgStream AthHistogramming::m_msg
privateinherited

Cached Message Stream.

Definition at line 250 of file AthHistogramming.h.

◆ m_name

std::string AthHistogramming::m_name
privateinherited

Instance name.

Definition at line 247 of file AthHistogramming.h.

◆ m_prefix

std::string AthHistogramAlgorithm::m_prefix
privateinherited

Name of the ROOT output stream (file).

Definition at line 88 of file AthHistogramAlgorithm.h.

◆ m_rndmSvc

ServiceHandle<IAthRNGSvc> MuonR4::BucketDumperAlg::m_rndmSvc {this, "RndmSvc", "AthRNGSvc", ""}
private

Definition at line 73 of file BucketDumperAlg.h.

73{this, "RndmSvc", "AthRNGSvc", ""};

◆ m_rootDir

std::string AthHistogramAlgorithm::m_rootDir
privateinherited

Name of the ROOT directory.

Definition at line 91 of file AthHistogramAlgorithm.h.

◆ m_segment_chiSquared

MuonVal::VectorBranch<float>& MuonR4::BucketDumperAlg::m_segment_chiSquared {m_tree.newVector<float>("segment_chiSquared")}
private

Definition at line 141 of file BucketDumperAlg.h.

141{m_tree.newVector<float>("segment_chiSquared")};

◆ m_segment_numberDoF

MuonVal::VectorBranch<float>& MuonR4::BucketDumperAlg::m_segment_numberDoF {m_tree.newVector<float>("segment_numberDoF")}
private

Definition at line 142 of file BucketDumperAlg.h.

142{m_tree.newVector<float>("segment_numberDoF")};

◆ m_segmentDir

MuonVal::ThreeVectorBranch MuonR4::BucketDumperAlg::m_segmentDir {m_tree, "segmentDirection"}
private

Definition at line 140 of file BucketDumperAlg.h.

140{m_tree, "segmentDirection"};

◆ m_segmentLocPhi

MuonVal::VectorBranch<float>& MuonR4::BucketDumperAlg::m_segmentLocPhi {m_tree.newVector<float>("segmentLocalPhi")}
private

Definition at line 135 of file BucketDumperAlg.h.

135{m_tree.newVector<float>("segmentLocalPhi")};

◆ m_segmentLocTheta

MuonVal::VectorBranch<float>& MuonR4::BucketDumperAlg::m_segmentLocTheta {m_tree.newVector<float>("segmentLocalTheta")}
private

Definition at line 134 of file BucketDumperAlg.h.

134{m_tree.newVector<float>("segmentLocalTheta")};

◆ m_segmentLocX

MuonVal::VectorBranch<float>& MuonR4::BucketDumperAlg::m_segmentLocX {m_tree.newVector<float>("segmentLocalX")}
private

Definition at line 132 of file BucketDumperAlg.h.

132{m_tree.newVector<float>("segmentLocalX")};

◆ m_segmentLocY

MuonVal::VectorBranch<float>& MuonR4::BucketDumperAlg::m_segmentLocY {m_tree.newVector<float>("segmentLocalY")}
private

Definition at line 133 of file BucketDumperAlg.h.

133{m_tree.newVector<float>("segmentLocalY")};

◆ m_segmentPos

MuonVal::ThreeVectorBranch MuonR4::BucketDumperAlg::m_segmentPos {m_tree, "segmentPosition"}
private

Definition at line 139 of file BucketDumperAlg.h.

139{m_tree, "segmentPosition"};

◆ m_segmentTruthIdx

MuonVal::VectorBranch<uint16_t>& MuonR4::BucketDumperAlg::m_segmentTruthIdx {m_tree.newVector<uint16_t>("segmentTruthPart")}
private

Definition at line 137 of file BucketDumperAlg.h.

137{m_tree.newVector<uint16_t>("segmentTruthPart")};

◆ m_spacePointKeys

SG::ReadHandleKeyArray<SpacePointContainer> MuonR4::BucketDumperAlg::m_spacePointKeys
private
Initial value:
{this, "SpacePointKeys", {"MuonSpacePoints"},
"Key to the space point container"}

Definition at line 58 of file BucketDumperAlg.h.

58 {this, "SpacePointKeys", {"MuonSpacePoints"},
59 "Key to the space point container"};

◆ m_spoint_adc

MuonVal::VectorBranch<uint16_t>& MuonR4::BucketDumperAlg::m_spoint_adc {m_tree.newVector<uint16_t>("adc")}
private

Definition at line 112 of file BucketDumperAlg.h.

112{m_tree.newVector<uint16_t>("adc")};

◆ m_spoint_covX

MuonVal::VectorBranch<float>& MuonR4::BucketDumperAlg::m_spoint_covX {m_tree.newVector<float>("covX")}
private

Definition at line 115 of file BucketDumperAlg.h.

115{m_tree.newVector<float>("covX")};

◆ m_spoint_covY

MuonVal::VectorBranch<float>& MuonR4::BucketDumperAlg::m_spoint_covY {m_tree.newVector<float>("covY")}
private

Definition at line 116 of file BucketDumperAlg.h.

116{m_tree.newVector<float>("covY")};

◆ m_spoint_dimension

MuonVal::VectorBranch<unsigned int>& MuonR4::BucketDumperAlg::m_spoint_dimension {m_tree.newVector<unsigned int>("dimension")}
private

Definition at line 124 of file BucketDumperAlg.h.

124{m_tree.newVector<unsigned int>("dimension")};

◆ m_spoint_driftR

MuonVal::VectorBranch<float>& MuonR4::BucketDumperAlg::m_spoint_driftR {m_tree.newVector<float>("driftR")}
private

Definition at line 117 of file BucketDumperAlg.h.

117{m_tree.newVector<float>("driftR")};

◆ m_spoint_globalPosition

MuonVal::ThreeVectorBranch MuonR4::BucketDumperAlg::m_spoint_globalPosition {m_tree, "globalPosition"}
private

Definition at line 105 of file BucketDumperAlg.h.

105{m_tree, "globalPosition"};

◆ m_spoint_id

MuonVal::MuonIdentifierBranch MuonR4::BucketDumperAlg::m_spoint_id {m_tree, "id"}
private

Definition at line 107 of file BucketDumperAlg.h.

107{m_tree, "id"};

◆ m_spoint_isMdt

MuonVal::VectorBranch<unsigned short>& MuonR4::BucketDumperAlg::m_spoint_isMdt {m_tree.newVector<unsigned short>("isMdt", false)}
private

Definition at line 110 of file BucketDumperAlg.h.

110{m_tree.newVector<unsigned short>("isMdt", false)};

◆ m_spoint_isStrip

MuonVal::VectorBranch<unsigned short>& MuonR4::BucketDumperAlg::m_spoint_isStrip {m_tree.newVector<unsigned short>("isStrip", false)}
private

Definition at line 109 of file BucketDumperAlg.h.

109{m_tree.newVector<unsigned short>("isStrip", false)};

◆ m_spoint_layer

MuonVal::VectorBranch<uint16_t>& MuonR4::BucketDumperAlg::m_spoint_layer {m_tree.newVector<uint16_t>("Layer")}
private

Definition at line 108 of file BucketDumperAlg.h.

108{m_tree.newVector<uint16_t>("Layer")};

◆ m_spoint_localPosition

MuonVal::ThreeVectorBranch MuonR4::BucketDumperAlg::m_spoint_localPosition {m_tree, "localPosition"}
private

Definition at line 104 of file BucketDumperAlg.h.

104{m_tree, "localPosition"};

◆ m_spoint_mat

MuonVal::MatrixBranch<int16_t>& MuonR4::BucketDumperAlg::m_spoint_mat {m_tree.newMatrix<int16_t>("sp_seg_matching",-1)}
private

Definition at line 127 of file BucketDumperAlg.h.

127{m_tree.newMatrix<int16_t>("sp_seg_matching",-1)};
setScaleOne setStatusOne setSaturated int16_t

◆ m_spoint_measuresEta

MuonVal::VectorBranch<unsigned short>& MuonR4::BucketDumperAlg::m_spoint_measuresEta {m_tree.newVector<unsigned short>("measuresEta")}
private

Definition at line 119 of file BucketDumperAlg.h.

119{m_tree.newVector<unsigned short>("measuresEta")};

◆ m_spoint_measuresPhi

MuonVal::VectorBranch<unsigned short>& MuonR4::BucketDumperAlg::m_spoint_measuresPhi {m_tree.newVector<unsigned short>("measuresPhi")}
private

Definition at line 120 of file BucketDumperAlg.h.

120{m_tree.newVector<unsigned short>("measuresPhi")};

◆ m_spoint_nEtaInstances

MuonVal::VectorBranch<unsigned int>& MuonR4::BucketDumperAlg::m_spoint_nEtaInstances {m_tree.newVector<unsigned int>("nEtaInUse")}
private

Definition at line 122 of file BucketDumperAlg.h.

122{m_tree.newVector<unsigned int>("nEtaInUse")};

◆ m_spoint_nPhiInstances

MuonVal::VectorBranch<unsigned int>& MuonR4::BucketDumperAlg::m_spoint_nPhiInstances {m_tree.newVector<unsigned int>("nPhiInUse")}
private

Definition at line 123 of file BucketDumperAlg.h.

123{m_tree.newVector<unsigned int>("nPhiInUse")};

◆ m_spoint_nSegments

MuonVal::VectorBranch<uint16_t>& MuonR4::BucketDumperAlg::m_spoint_nSegments {m_tree.newVector<uint16_t>("nSegments")}
private

Definition at line 126 of file BucketDumperAlg.h.

126{m_tree.newVector<uint16_t>("nSegments")};

◆ m_spoint_tdc

MuonVal::VectorBranch<uint16_t>& MuonR4::BucketDumperAlg::m_spoint_tdc {m_tree.newVector<uint16_t>("tdc")}
private

Definition at line 113 of file BucketDumperAlg.h.

113{m_tree.newVector<uint16_t>("tdc")};

◆ m_spoint_trueLabel

MuonVal::VectorBranch<unsigned short>& MuonR4::BucketDumperAlg::m_spoint_trueLabel {m_tree.newVector<unsigned short>("trueLabel")}
private

Definition at line 121 of file BucketDumperAlg.h.

121{m_tree.newVector<unsigned short>("trueLabel")};

◆ m_spoint_trueSeg

MuonVal::MatrixBranch<int16_t>& MuonR4::BucketDumperAlg::m_spoint_trueSeg {m_tree.newMatrix<int16_t>("spacePoint_truthSegLink",-1)}
private

Definition at line 128 of file BucketDumperAlg.h.

128{m_tree.newMatrix<int16_t>("spacePoint_truthSegLink",-1)};

◆ m_streamName

Gaudi::Property<std::string> MuonR4::BucketDumperAlg::m_streamName {this, "StreamName", ""}
private

Definition at line 68 of file BucketDumperAlg.h.

68{this, "StreamName", ""};

◆ m_tree

MuonVal::MuonTesterTree MuonR4::BucketDumperAlg::m_tree {"MuonBucketDump","MuonBucketDump"}
private

Definition at line 83 of file BucketDumperAlg.h.

83{"MuonBucketDump","MuonBucketDump"};

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

SG::ReadDecorHandleKeyArray<xAOD::MuonSegmentContainer> MuonR4::BucketDumperAlg::m_truthDecorKeys {this, "TruthDecorLinks", {}}
private

Definition at line 63 of file BucketDumperAlg.h.

63{this, "TruthDecorLinks", {}};

◆ m_truthSegLocPhi

MuonVal::VectorBranch<float>& MuonR4::BucketDumperAlg::m_truthSegLocPhi {m_tree.newVector<float>("truthSegmentLocalPhi")}
private

Definition at line 147 of file BucketDumperAlg.h.

147{m_tree.newVector<float>("truthSegmentLocalPhi")};

◆ m_truthSegLocTheta

MuonVal::VectorBranch<float>& MuonR4::BucketDumperAlg::m_truthSegLocTheta {m_tree.newVector<float>("truthSegmentLocalTheta")}
private

Definition at line 146 of file BucketDumperAlg.h.

146{m_tree.newVector<float>("truthSegmentLocalTheta")};

◆ m_truthSegLocX

MuonVal::VectorBranch<float>& MuonR4::BucketDumperAlg::m_truthSegLocX {m_tree.newVector<float>("truthSegmentLocalX")}
private

Definition at line 144 of file BucketDumperAlg.h.

144{m_tree.newVector<float>("truthSegmentLocalX")};

◆ m_truthSegLocY

MuonVal::VectorBranch<float>& MuonR4::BucketDumperAlg::m_truthSegLocY {m_tree.newVector<float>("truthSegmentLocalY")}
private

Definition at line 145 of file BucketDumperAlg.h.

145{m_tree.newVector<float>("truthSegmentLocalY")};

◆ m_varHandleArraysDeclared

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

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

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

Definition at line 398 of file AthCommonDataStore.h.

◆ m_visionTool

ToolHandle<MuonValR4::IPatternVisualizationTool> MuonR4::BucketDumperAlg::m_visionTool {this, "VisualizationTool", ""}
private

Pattern visualization tool.

Definition at line 80 of file BucketDumperAlg.h.

80{this, "VisualizationTool", ""};

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