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

Private Types

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

Private Member Functions

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

Private Attributes

SG::ReadHandleKeyArray< SpacePointContainer > m_spacePointKeys
SG::ReadHandleKeyArray< xAOD::MuonSegmentContainer > m_inSegmentKeys {this, "SegmentKey", {"MsSegments"}}
SG::ReadDecorHandleKeyArray< xAOD::MuonSegmentContainer > m_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< IAthRNGSvc > m_rndmSvc {this, "RndmSvc", "AthRNGSvc", ""}
ToolHandle< MuonML::IGraphInferenceTool > m_inferenceTool
 Inference tool for ML bucket scoring (optional).
ToolHandle< MuonValR4::IPatternVisualizationTool > m_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 199 of file AthHistogramming.h.

◆ GraphMap_t

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

Typedef for convenience.

Definition at line 213 of file AthHistogramming.h.

◆ hash_t

typedef uint32_t AthHistogramming::hash_t
privateinherited

typedef for the internal hash

Definition at line 171 of file AthHistogramming.h.

◆ HistMap_t

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

Typedef for convenience.

Definition at line 192 of file AthHistogramming.h.

◆ StoreGateSvc_t

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

Definition at line 376 of file AthCommonDataStore.h.

◆ TreeMap_t

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

Typedef for convenience.

Definition at line 206 of file AthHistogramming.h.

Constructor & Destructor Documentation

◆ ~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,
std::string_view tDir = "",
std::string_view stream = "" )
inlineprotectedinherited

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

Definition at line 329 of file AthHistogramming.h.

330{
331 // We need to create a non-const clone
332 TEfficiency* effClone = dynamic_cast< TEfficiency* >( eff.Clone() );
333 if ( !effClone ) {
334 m_msg << MSG::ERROR << "Couldn't create a TEfficiency clone" << endmsg;
335 return StatusCode::FAILURE;
336 }
337 return this->book( *effClone, tDir, stream );
338}
#define endmsg
StatusCode book(const TH1 &hist, std::string_view tDir="", std::string_view stream="")
Simplify the booking and registering (into THistSvc) of histograms.
MsgStream m_msg
Cached Message Stream.

◆ book() [2/8]

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

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

◆ book() [3/8]

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

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

Definition at line 297 of file AthHistogramming.h.

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

◆ book() [4/8]

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

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

Definition at line 397 of file AthHistogramming.h.

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

◆ book() [5/8]

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

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

Definition at line 350 of file AthHistogramming.h.

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

◆ book() [6/8]

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

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

Definition at line 340 of file AthHistogramming.h.

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

◆ book() [7/8]

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

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

Definition at line 318 of file AthHistogramming.h.

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

◆ book() [8/8]

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

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

Definition at line 308 of file AthHistogramming.h.

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

◆ bookGetPointer() [1/8]

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

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

Definition at line 275 of file AthHistogramming.h.

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

◆ bookGetPointer() [2/8]

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

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

Definition at line 374 of file AthHistogramming.cxx.

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

◆ bookGetPointer() [3/8]

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

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

Definition at line 254 of file AthHistogramming.h.

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

◆ bookGetPointer() [4/8]

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

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

Definition at line 268 of file AthHistogramming.cxx.

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

◆ bookGetPointer() [5/8]

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

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

Definition at line 139 of file AthHistogramming.cxx.

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

◆ bookGetPointer() [6/8]

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

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

Definition at line 287 of file AthHistogramming.h.

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

◆ bookGetPointer() [7/8]

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

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

Definition at line 97 of file AthHistogramming.cxx.

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

◆ bookGetPointer() [8/8]

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

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

Definition at line 266 of file AthHistogramming.h.

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

◆ buildBookingString()

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

Method to build individual booking string.

Definition at line 500 of file AthHistogramming.cxx.

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

◆ 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_DEBUG(x,...)
#define ATH_MSG_ERROR(x,...)
std::vector< std::vector< float > > scores
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 to fill the internal configuration.

Definition at line 64 of file AthHistogramming.cxx.

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

◆ declareGaudiProperty()

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

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

Definition at line 156 of file AthCommonDataStore.h.

158 {
160 hndl.value(),
161 hndl.documentation());
162
163 }

◆ declareProperty()

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

Definition at line 145 of file AthCommonDataStore.h.

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

◆ detStore()

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

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

Definition at line 95 of file AthCommonDataStore.h.

◆ 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_WARNING(x,...)
#define ATH_MSG_VERBOSE(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.
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 ( std::string_view effName,
std::string_view tDir = "",
std::string_view stream = "" )
protectedinherited

Simplify the retrieval of registered TEfficiency.

Definition at line 214 of file AthHistogramming.cxx.

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

◆ evtStore()

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

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

Definition at line 85 of file AthCommonDataStore.h.

◆ execute()

StatusCode 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 ( std::string_view graphName,
std::string_view tDir = "",
std::string_view stream = "" )
protectedinherited

Simplify the retrieval of registered TGraphs.

Definition at line 438 of file AthHistogramming.cxx.

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

◆ hash()

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

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

Definition at line 422 of file AthHistogramming.h.

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

◆ hist()

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

Simplify the retrieval of registered histograms of any type.

Definition at line 182 of file AthHistogramming.cxx.

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

◆ hist2d()

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

Simplify the retrieval of registered 2-d histograms.

Definition at line 363 of file AthHistogramming.h.

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

◆ hist3d()

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

Simplify the retrieval of registered 3-d histograms.

Definition at line 380 of file AthHistogramming.h.

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

◆ histSvc()

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

The standard THistSvc (for writing histograms and TTrees and more to a root file) Returns (kind of) a pointer to the THistSvc.

Definition at line 113 of file AthHistogramAlgorithm.h.

114{
115 return m_histSvc;
116}

◆ initialize()

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

◆ outputHandles()

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

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ renounce()

std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::renounce ( T & h)
inlineprotectedinherited

Definition at line 368 of file AthCommonDataStore.h.

369 {
370 h.renounce();
372 }
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > renounce(T &h)

◆ renounceArray()

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

remove all handles from I/O resolution

Definition at line 352 of file AthCommonDataStore.h.

352 {
354 }

◆ setFilterPassed()

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

Set filter decision:

Definition at line 99 of file AthCommonAlgorithm.h.

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

◆ sysExecute()

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

Execute an algorithm.

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

Reimplemented in AthAnalysisAlgorithm.

Definition at line 80 of file AthCommonAlgorithm.cxx.

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

◆ sysInitialize()

StatusCode AthHistogramAlgorithm::sysInitialize ( )
virtualinherited

Initialization method invoked by the framework.

This method is responsible for any bookkeeping of initialization required by the framework itself. It will in turn invoke the initialize() method of the derived algorithm, and of any sub-algorithms which it creates.

Reimplemented from AthCommonAlgorithm< Gaudi::Algorithm >.

Reimplemented in AthAnalysisAlgorithm.

Definition at line 75 of file AthHistogramAlgorithm.cxx.

76{
77 // ---- stolen from GaudiKernel/Algorithm::sysInitialize -------
78 // Bypass the initialization if the algorithm
79 // has already been initialized.
80 if ( Gaudi::StateMachine::INITIALIZED <= FSMState() ) return StatusCode::SUCCESS;
81
82 // Set the Algorithm's properties
83 bindPropertiesTo( serviceLocator()->getOptsSvc() );
84
85 // Bypass the initialization if the algorithm is disabled.
86 // Need to do this after bindPropertiesTo.
87 if ( !isEnabled( ) ) return StatusCode::SUCCESS;
88
89 // ---- stolen from GaudiKernel/Algorithm::sysInitialize ------- END ---
90
91
92 // Get the THistSvc
93 ATH_CHECK ( histSvc().retrieve() );
94
95 // Configure the underlying AthHistogramming helper
100
101 // Print some setup information into the log file
102 ATH_MSG_DEBUG ("Initializing {}...", name());
103 ATH_MSG_DEBUG (" using THistService = " << m_histSvc );
104 ATH_MSG_DEBUG (" using RootStreamName = {}", m_prefix );
105 ATH_MSG_DEBUG (" using RootDirName = {}", m_rootDir );
106 ATH_MSG_DEBUG (" using HistNamePrefix = {}", m_histNamePrefix );
107 ATH_MSG_DEBUG (" using HistNamePostfix = {}", m_histNamePostfix );
108 ATH_MSG_DEBUG (" using HistTitlePrefix = {}", m_histTitlePrefix );
109 ATH_MSG_DEBUG (" using HistTitlePostfix = {}", m_histTitlePostfix );
110
111
112 // re-direct to base class...
114}
virtual StatusCode sysInitialize() override
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 to fill the internal configuration.

◆ sysStart()

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

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

◆ tree()

TTree * AthHistogramming::tree ( std::string_view treeName,
std::string_view tDir = "",
std::string_view stream = "" )
protectedinherited

Simplify the retrieval of registered TTrees.

Definition at line 333 of file AthHistogramming.cxx.

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

◆ updateVHKA()

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

Definition at line 298 of file AthCommonDataStore.h.

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

Member Data Documentation

◆ m_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 381 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 202 of file AthHistogramming.h.

◆ m_evtStore

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

Pointer to StoreGate (event store by default).

Definition at line 378 of file AthCommonDataStore.h.

◆ m_extendedExtraObjects

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

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

Empty if no symlinks were found.

Definition at line 108 of file AthCommonAlgorithm.h.

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

◆ m_histMap

HistMap_t AthHistogramming::m_histMap
privateinherited

The map of histogram names to their pointers.

Definition at line 195 of file AthHistogramming.h.

◆ m_histNamePostfix

std::string AthHistogramAlgorithm::m_histNamePostfix
privateinherited

The postfix for the histogram THx name.

Definition at line 97 of file AthHistogramAlgorithm.h.

◆ m_histNamePrefix

std::string AthHistogramAlgorithm::m_histNamePrefix
privateinherited

The prefix for the histogram THx name.

Definition at line 94 of file AthHistogramAlgorithm.h.

◆ m_histSvc

ServiceHandle<ITHistSvc> AthHistogramAlgorithm::m_histSvc
privateinherited

Default constructor: AthHistogramAlgorithm();.

a handle on the Hist/TTree registration service

Definition at line 83 of file AthHistogramAlgorithm.h.

◆ m_histTitlePostfix

std::string AthHistogramAlgorithm::m_histTitlePostfix
privateinherited

The postfix for the histogram THx title.

Definition at line 103 of file AthHistogramAlgorithm.h.

◆ m_histTitlePrefix

std::string AthHistogramAlgorithm::m_histTitlePrefix
privateinherited

The prefix for the histogram THx title.

Definition at line 100 of file AthHistogramAlgorithm.h.

◆ m_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", {"MsSegments"}}
private

Definition at line 61 of file BucketDumperAlg.h.

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

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

◆ m_name

std::string AthHistogramming::m_name
privateinherited

Instance name.

Definition at line 239 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 209 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 387 of file AthCommonDataStore.h.

◆ m_vhka

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

Definition at line 386 of file AthCommonDataStore.h.

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