ATLAS Offline Software
Loading...
Searching...
No Matches
MuonValR4::MsTrackTester Class Reference

#include <MsTrackTester.h>

Inheritance diagram for MuonValR4::MsTrackTester:

Public Member Functions

StatusCode initialize () override final
StatusCode execute (const EventContext &ctx) override final
 Execute method.
StatusCode finalize () override final
 AthHistogramAlgorithm (const std::string &name, ISvcLocator *pSvcLocator)
 Constructor with parameters:
virtual StatusCode sysInitialize ()
 Initialization method invoked by the framework.
const ServiceHandle< ITHistSvc > & histSvc () const
 The standard THistSvc (for writing histograms and TTrees and more to a root file) Returns (kind of) a pointer to the THistSvc.
virtual bool isClonable () const override
 Specify if the algorithm is clonable.
virtual StatusCode sysExecute (const EventContext &ctx) override
 Execute an algorithm.
virtual const DataObjIDColl & extraOutputDeps () const override
 Return the list of extra output dependencies.
virtual bool filterPassed (const EventContext &ctx) const
 Get filter decision:
virtual void setFilterPassed (bool state, const EventContext &ctx) const
 Set filter decision:
ServiceHandle< StoreGateSvc > & evtStore ()
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
const ServiceHandle< StoreGateSvc > & detStore () const
 The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
virtual StatusCode sysStart () override
 Handle START transition.
virtual std::vector< Gaudi::DataHandle * > inputHandles () const override
 Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * > outputHandles () const override
 Return this algorithm's output handles.
Gaudi::Details::PropertyBase & declareProperty (Gaudi::Property< T, V, H > &t)
void updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream & msg () const
bool msgLvl (const MSG::Level lvl) const

Protected Member Functions

virtual bool isReEntrant () const override final
 Legacy algorithms are not thread-safe.
void renounceArray (SG::VarHandleKeyArray &handlesArray)
 remove all handles from I/O resolution
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > renounce (T &h)
void extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
 Add StoreName to extra input/output deps as needed.
StatusCode configAthHistogramming (const ServiceHandle< ITHistSvc > &histSvc, const std::string &prefix, const std::string &rootDir, const std::string &histNamePrefix, const std::string &histNamePostfix, const std::string &histTitlePrefix, const std::string &histTitlePostfix)
 To be called by the derived classes to fill the internal configuration.
TH1 * bookGetPointer (const TH1 &hist, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of histograms.
TH1 * bookGetPointer (TH1 *hist, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of histograms.
TH1 * bookGetPointer (TH1 &histRef, std::string tDir="", std::string stream="")
 Simplify the booking and registering (into THistSvc) of histograms.
TTree * bookGetPointer (const TTree &treeRef, std::string tDir="", std::string stream="")
 Simplify the booking and registering (into THistSvc) of TTrees.
TGraph * bookGetPointer (const TGraph &graphRef, std::string tDir="", std::string stream="")
 Simplify the booking and registering (into THistSvc) of TGraphs.
TEfficiency * bookGetPointer (const TEfficiency &eff, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency.
TEfficiency * bookGetPointer (TEfficiency *eff, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency.
TEfficiency * bookGetPointer (TEfficiency &effRef, std::string tDir="", std::string stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency.
StatusCode book (const TH1 &hist, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of histograms.
StatusCode book (TH1 *hist, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of histograms.
StatusCode book (TH1 &histRef, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of histograms.
StatusCode book (const TTree &treeRef, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of TTrees.
StatusCode book (const TGraph &graphRef, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of TGraphs.
StatusCode book (const TEfficiency &eff, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency.
StatusCode book (TEfficiency *eff, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency.
StatusCode book (TEfficiency &effRef, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency.
TH1 * hist (const std::string &histName, const std::string &tDir="", const std::string &stream="")
 Simplify the retrieval of registered histograms of any type.
TH2 * hist2d (const std::string &histName, const std::string &tDir="", const std::string &stream="")
 Simplify the retrieval of registered 2-d histograms.
TH3 * hist3d (const std::string &histName, const std::string &tDir="", const std::string &stream="")
 Simplify the retrieval of registered 3-d histograms.
TTree * tree (const std::string &treeName, const std::string &tDir="", const std::string &stream="")
 Simplify the retrieval of registered TTrees.
TGraph * graph (const std::string &graphName, const std::string &tDir="", const std::string &stream="")
 Simplify the retrieval of registered TGraphs.
TEfficiency * efficiency (const std::string &effName, const std::string &tDir="", const std::string &stream="")
 Simplify the retrieval of registered TEfficiency.

Private Types

using Location = MuonR4::MsTrackSeed::Location
using SectorProjector = MuonR4::ExpandedSector::SectorProjector
using TruthHitCol = std::unordered_set<const xAOD::MuonSimHit*>
using SegmentKey_t = SG::ReadHandleKey<xAOD::MuonSegmentContainer>
 Abrivate the ReadHandleKey_t for the segment container.
using TrackKey_t = SG::ReadHandleKey<xAOD::TrackParticleContainer>
 Abrivate the key type for the track particle container.
using MuonKey_t = SG::ReadHandleKey<xAOD::MuonContainer>
 Abrivate the key type for the muon container.
using ParticleBranchPtr_t = std::shared_ptr<MuonVal::IParticleFourMomBranch>
using SegmentBranchPtr_t = std::shared_ptr<MuonPRDTest::SegmentVariables>
typedef ServiceHandle< StoreGateSvcStoreGateSvc_t
typedef uint32_t hash_t
 typedef for the internal hash
typedef std::map< const hash_t, TH1 * > HistMap_t
 Typedef for convenience.
typedef std::map< const hash_t, TEfficiency * > EffMap_t
 Typedef for convenience.
typedef std::map< const hash_t, TTree * > TreeMap_t
 Typedef for convenience.
typedef std::map< const hash_t, TGraph * > GraphMap_t
 Typedef for convenience.

Private Member Functions

std::optional< MuonR4::MsTrackSeedmakeSeedFromTruth (const ActsTrk::GeometryContext &gctx, const xAOD::TruthParticle &truthMuon) const
 Construct MS track seed from the truth associated segments.
std::pair< double, double > calcSeedLength (const ActsTrk::GeometryContext &gctx, const MuonR4::MsTrackSeed &seed) const
 Calculate the length of the seed and the theta deflection angle The length is defined as the spread of the seed's segments in the cylinder coordinate.
StatusCode dumpLegacyTracks (const EventContext &ctx)
 Dumps the legacy containers to the TTree.
StatusCode dumpTruthContent (const EventContext &ctx)
 Dump truth information.
StatusCode dumpRecoContent (const EventContext &ctx)
 Dump the reconstructed information.
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey>
void buildBookingString (std::string &bookingString, std::string &histName, std::string &tDir, std::string &stream, bool usePrefixPostfix=false)
 Method to build individual booking string.
void myReplace (std::string &str, const std::string &oldStr, const std::string &newStr)
 Helper method to replace sub-string.
hash_t hash (const std::string &histName) const
 Method to calculate a 32-bit hash from a string.

Private Attributes

MuonVal::MuonTesterTree m_tree {"MsTrackValidTest", "MuonTrackTester"}
const MuonGMR4::MuonDetectorManagerm_detMgr {nullptr}
Gaudi::Property< bool > m_isMC {this, "isMC", false}
SegmentKey_t m_truthSegmentKey {this, "TruthSegmentKey", "MuonTruthSegments"}
 Segment from the truth hits.
SegmentKey_t m_recoSegmentKey {this, "SegmentKey", "MuonSegmentsFromR4"}
 Primary segment container.
SegmentKey_t m_legacySegmentKey {this, "LegacySegmentKey", "MuonSegments"}
 Legacy segment container.
SG::ReadHandleKey< xAOD::TruthParticleContainerm_truthKey {this, "TruthKey", "MuonTruthParticles"}
 Key to the truth particle collection.
SG::ReadDecorHandleKeyArray< SG::AuxVectorBasem_trkTruthLinks {this, "TruthTrackLinks", {}}
 Decoration dependency to the MS truth track links.
SG::ReadHandleKey< MuonR4::MsTrackSeedContainerm_msTrkSeedKey {this, "MsTrkSeedKey", "MsTrackSeeds"}
 Temporary container write handle to push the seeds to store gate for later efficiency analysis.
ActsTrk::GeoContextReadKey_t m_geoCtxKey {this, "AlignmentKey", "ActsAlignment", "cond handle key"}
 Dependency on the geometry alignment.
MuonKey_t m_muonKey {this, "MuonKey", "MuonsR4"}
 Dependency on the R4 muon container.
ToolHandle< MuonR4::ITrackSummaryToolm_summaryTool {this, "SummaryTool" ,""}
 Hit summary tool.
ToolHandle< MuonR4::ITrackSeedingDiagnosticsToolm_seedingTool {this, "SeedingTool", ""}
 The track seeding tool to construct the seed candidates and to estimate the initial parameters.
ToolHandle< MuonR4::ISegmentSelectionToolm_segSelector {this, "SegmentSelectionTool" , "" }
 Selection tool to quantify the segment candidate quality.
TrackKey_t m_legacyTrackKey {this,"LegacyTrackKey", "MuonSpectrometerTrackParticles"}
 Legacy track reconstruction chain.
SG::ReadCondHandleKey< AtlasFieldCacheCondObjm_fieldCacheKey {this, "MagFieldKey", "fieldCondObj", "Name of the Magnetic Field conditions object key"}
 Dependency on the magnetic field.
MuonKey_t m_legacyMuonKey {this,"LegacyMuonKey", "Muons"}
 Legacy muons.
ParticleBranchPtr_t m_truthTrks {}
ParticleBranchPtr_t m_muonTrks {}
 Stored muon information from the Acts muon reco chain.
SegmentBranchPtr_t m_truthSegs {}
SegmentBranchPtr_t m_recoSegs {}
SegmentBranchPtr_t m_legacyRecoSegs {}
MuonVal::ThreeVectorBranch m_seedPos {m_tree, "MsTrkSeed_position"}
 Simple seed information.
MuonVal::UnitThreeVectorBranch m_seedDir {m_tree, "MsTrkSeed_direction"}
 Seed direction vector.
MuonVal::VectorBranch< char > & m_seedType {m_tree.newVector<char>("MsTrkSeed_type")}
 Is the seed in the encap or in the barrel chambers.
MuonVal::VectorBranch< int > & m_seedSector {m_tree.newVector<int>("MsTrkSeed_sector")}
 Sector of the seed, even center, odd overlap regions, for details see:
MuonVal::VectorBranch< float > & m_seedLength {m_tree.newVector<float>("MsTrkSeed_length")}
 Maximum separation between the segments on the reference plane.
MuonVal::VectorBranch< float > & m_seedThetaCone {m_tree.newVector<float>("MsTrkSeed_thetaCone")}
 Maximum angular difference between the segments part of the seed.
MuonVal::VectorBranch< char > & m_seedGood {m_tree.newVector<char>("MsTrkSeed_goodSeed")}
 Does the seeding tool construct valid parameters from the seed.
MuonVal::VectorBranch< float > & m_seedQP {m_tree.newVector<float>("MsTrkSeed_qTimesP")}
 Estimated momentum times charge from the track seed.
MuonVal::VectorBranch< unsigned short > & m_seedTruthLink {m_tree.newVector<unsigned short>("MsTrkSeed_truthLink", -1)}
 Link to the truth muon.
std::shared_ptr< TrackSummaryModulem_seedSummary {}
 Hit summary on the track seed.
MuonVal::MatrixBranch< unsigned short > & m_seedRecoSegMatch {m_tree.newMatrix<unsigned short>("MsTrkSeed_segmentLinks")}
 Link of the track seed to the building segment.
MuonVal::MatrixBranch< unsigned short > & m_truthSegToRecoLink {m_tree.newMatrix<unsigned short>("TruthSegments_recoSegLinks",-1)}
 Link of the truth segments to the matchin reco segments.
MuonVal::VectorBranch< unsigned short > & m_legacySegToTrkLinks {m_tree.newVector<unsigned short>("LegacyRecoSegments_trkLinks", -1)}
 Link of the legacy track to the legacy segment.
MuonVal::MatrixBranch< unsigned short > & m_truthMuToSeedIdx {m_tree.newMatrix<unsigned short>("TruthMuons_seedLinks", -1)}
 Links to all MsTrkSeeds that could be matched to the truthMuon, i.e.
MuonVal::MatrixBranch< unsigned short > & m_truthMuToSeedCounter {m_tree.newMatrix<unsigned short>("TruthMuons_seedNSeg")}
 Number of matched segments in the seed.
MuonVal::MatrixBranch< unsigned short > & m_truthMuRecoSegLinks {m_tree.newMatrix<unsigned short>("TruthMuons_recoSegLinks")}
 Links from the truth muon to the segments.
ParticleBranchPtr_t m_legacyTrks {}
 Output branches of the legacy MS tracks.
ServiceHandle< ITHistSvc > m_histSvc
 Default constructor: AthHistogramAlgorithm();.
std::string m_prefix
 Name of the ROOT output stream (file).
std::string m_rootDir
 Name of the ROOT directory.
std::string m_histNamePrefix
 The prefix for the histogram THx name.
std::string m_histNamePostfix
 The postfix for the histogram THx name.
std::string m_histTitlePrefix
 The prefix for the histogram THx title.
std::string m_histTitlePostfix
 The postfix for the histogram THx title.
DataObjIDColl m_extendedExtraObjects
 Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default).
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default).
std::vector< SG::VarHandleKeyArray * > m_vhka
bool m_varHandleArraysDeclared
HistMap_t m_histMap
 The map of histogram names to their pointers.
EffMap_t m_effMap
 The map of histogram names to their pointers.
TreeMap_t m_treeMap
 The map of TTree names to their pointers.
GraphMap_t m_graphMap
 The map of TGraph names to their pointers.
std::string m_streamName
 Name of the ROOT output stream (file).
std::string m_name
 Instance name.
MsgStream m_msg
 Cached Message Stream.

Detailed Description

Definition at line 33 of file MsTrackTester.h.

Member Typedef Documentation

◆ EffMap_t

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

Typedef for convenience.

Definition at line 207 of file AthHistogramming.h.

◆ GraphMap_t

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

Typedef for convenience.

Definition at line 221 of file AthHistogramming.h.

◆ hash_t

typedef uint32_t AthHistogramming::hash_t
privateinherited

typedef for the internal hash

Definition at line 171 of file AthHistogramming.h.

◆ HistMap_t

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

Typedef for convenience.

Definition at line 200 of file AthHistogramming.h.

◆ Location

◆ MuonKey_t

Abrivate the key type for the muon container.

Definition at line 79 of file MsTrackTester.h.

◆ ParticleBranchPtr_t

Definition at line 111 of file MsTrackTester.h.

◆ SectorProjector

◆ SegmentBranchPtr_t

Definition at line 116 of file MsTrackTester.h.

◆ SegmentKey_t

Abrivate the ReadHandleKey_t for the segment container.

Definition at line 75 of file MsTrackTester.h.

◆ StoreGateSvc_t

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

Definition at line 388 of file AthCommonDataStore.h.

◆ TrackKey_t

Abrivate the key type for the track particle container.

Definition at line 77 of file MsTrackTester.h.

◆ TreeMap_t

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

Typedef for convenience.

Definition at line 214 of file AthHistogramming.h.

◆ TruthHitCol

using MuonValR4::MsTrackTester::TruthHitCol = std::unordered_set<const xAOD::MuonSimHit*>
private

Definition at line 73 of file MsTrackTester.h.

Member Function Documentation

◆ AthHistogramAlgorithm()

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

Constructor with parameters:

Definition at line 40 of file AthHistogramAlgorithm.cxx.

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

◆ book() [1/8]

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

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

Definition at line 337 of file AthHistogramming.h.

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

◆ book() [2/8]

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

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

◆ book() [3/8]

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

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

Definition at line 305 of file AthHistogramming.h.

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

◆ book() [4/8]

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

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

Definition at line 405 of file AthHistogramming.h.

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

◆ book() [5/8]

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

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

Definition at line 358 of file AthHistogramming.h.

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

◆ book() [6/8]

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

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

Definition at line 348 of file AthHistogramming.h.

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

◆ book() [7/8]

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

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

Definition at line 326 of file AthHistogramming.h.

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

◆ book() [8/8]

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

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

Definition at line 316 of file AthHistogramming.h.

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

◆ bookGetPointer() [1/8]

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

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

Definition at line 283 of file AthHistogramming.h.

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

◆ bookGetPointer() [2/8]

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

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

Definition at line 388 of file AthHistogramming.cxx.

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

◆ bookGetPointer() [3/8]

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

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

Definition at line 262 of file AthHistogramming.h.

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

◆ bookGetPointer() [4/8]

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

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

Definition at line 273 of file AthHistogramming.cxx.

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

◆ bookGetPointer() [5/8]

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

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

Definition at line 123 of file AthHistogramming.cxx.

124{
125 std::string effName(effRef.GetName());
126 const std::string effTitle(effRef.GetTitle());
127 std::string bookingString;
128 this->buildBookingString(bookingString, effName, tDir, stream);
129 const std::string finalEffName = m_histNamePrefix + effName + m_histNamePostfix;
130 effRef.SetTitle((m_histTitlePrefix + effTitle + m_histTitlePostfix).c_str());
131 effRef.SetName(finalEffName.c_str());
132
133 const hash_t effHash = this->hash(effName);
134 EffMap_t::const_iterator it = m_effMap.find(effHash);
135 if (it != m_effMap.end()) {
136 m_msg << MSG::WARNING
137 << "Detected a hash collision. The hash for the TEfficiency with name=" << effName
138 << " already exists and points to a TEfficiency with name=" << it->second->GetName()
139 << " NOT going to book the new TEfficiency and returning a NULL pointer!" << endmsg;
140 return nullptr;
141 }
142
143 if (!histSvc()->regEfficiency(bookingString, &effRef).isSuccess()) {
144 m_msg << MSG::WARNING
145 << "Problem registering TEfficiency with name " << effName
146 << ", name prefix " << m_histNamePrefix
147 << ", title " << effTitle
148 << ", title prefix " << m_histTitlePrefix
149 << ", and title postfix " << m_histTitlePostfix
150 << " in " << m_name << "!" << endmsg;
151 return nullptr;
152 }
153
154 m_effMap.insert(m_effMap.end(), std::pair<const hash_t, TEfficiency*>(effHash, &effRef));
155
156 return &effRef;
157}
EffMap_t m_effMap
The map of histogram names to their pointers.

◆ bookGetPointer() [6/8]

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

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

Definition at line 295 of file AthHistogramming.h.

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

◆ bookGetPointer() [7/8]

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

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

Definition at line 83 of file AthHistogramming.cxx.

84{
85 std::string histName(histRef.GetName());
86 const std::string histTitle(histRef.GetTitle());
87 std::string bookingString;
88
89 this->buildBookingString(bookingString, histName, tDir, stream);
90
91 const std::string finalHistName = m_histNamePrefix + histName + m_histNamePostfix;
92
93 histRef.SetTitle((m_histTitlePrefix + histTitle + m_histTitlePostfix).c_str());
94 histRef.SetName(finalHistName.c_str());
95
96 const hash_t histHash = this->hash(histName);
97 HistMap_t::const_iterator it = m_histMap.find(histHash);
98 if (it != m_histMap.end()) {
99 m_msg << MSG::WARNING
100 << "Detected a hash collision. The hash for the histogram with name=" << histName
101 << " already exists and points to a histogram with name=" << it->second->GetName()
102 << " NOT going to book the new histogram and returning a NULL pointer!" << endmsg;
103 return nullptr;
104 }
105
106 if (!histSvc()->regHist(bookingString, &histRef).isSuccess()) {
107 m_msg << MSG::WARNING
108 << "Problem registering histogram with name " << histName
109 << ", name prefix " << m_histNamePrefix
110 << ", title " << histTitle
111 << ", title prefix " << m_histTitlePrefix
112 << ", and title postfix " << m_histTitlePostfix
113 << " in " << m_name << "!" << endmsg;
114 return nullptr;
115 }
116
117 m_histMap.insert(m_histMap.end(), std::pair<const hash_t, TH1*>(histHash, &histRef));
118
119 return &histRef;
120}
HistMap_t m_histMap
The map of histogram names to their pointers.

◆ bookGetPointer() [8/8]

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

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

Definition at line 274 of file AthHistogramming.h.

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

◆ buildBookingString()

void AthHistogramming::buildBookingString ( std::string & bookingString,
std::string & histName,
std::string & tDir,
std::string & stream,
bool usePrefixPostfix = false )
privateinherited

Method to build individual booking string.

Definition at line 521 of file AthHistogramming.cxx.

526{
527 // Massage the final string to book things
528 if(tDir.empty()) tDir = m_rootDir;
529 size_t pos = histName.rfind('/');
530 if(pos != std::string::npos){
531 tDir+='/';
532 tDir.append(histName, 0,pos);
533 histName.erase(0,pos+1);
534 };
535 if(stream.empty()) stream = m_streamName;
536
537 if(usePrefixPostfix){
538 bookingString = "/"+stream+"/"+tDir+"/"+m_histNamePrefix+histName+m_histNamePostfix;
539 } else {
540 bookingString = "/"+stream+"/"+tDir+"/"+histName;
541 }
542 while(bookingString.find("//") != std::string::npos){
543 this->myReplace(bookingString,"//","/");
544 }
545
546 return;
547}
std::string m_rootDir
Name of the ROOT directory.
std::string m_streamName
Name of the ROOT output stream (file).
void myReplace(std::string &str, const std::string &oldStr, const std::string &newStr)
Helper method to replace sub-string.

◆ calcSeedLength()

std::pair< double, double > MuonValR4::MsTrackTester::calcSeedLength ( const ActsTrk::GeometryContext & gctx,
const MuonR4::MsTrackSeed & seed ) const
private

Calculate the length of the seed and the theta deflection angle The length is defined as the spread of the seed's segments in the cylinder coordinate.

The deflection angle is calculates as the spread of the theta angles of the individual segments

Parameters
gctxThe geometry context to align the segment within ATLAS
seedThe seed with the contributing segments

Definition at line 67 of file MsTrackTester.cxx.

68 {
69 double maxL{-1.*Gaudi::Units::km}, minL{1.*Gaudi::Units::km},
70 maxTheta{-181.}, minTheta{181};
71 for (const xAOD::MuonSegment* seg : seed.segments()) {
72 const Amg::Vector2D projPos{m_seedingTool->expressOnCylinder(gctx.context(),*seg, seed.location(), seed.sector())};
73 if (!m_seedingTool->withinBounds(projPos, seed.location())) {
74 continue;
75 }
76 const double projected = projPos[seed.location()==Location::Barrel];
77 const double theta = seg->direction().theta();
78 minL = std::min(minL, projected);
79 maxL = std::max(maxL, projected);
80 minTheta = std::min(minTheta, toDeg(theta));
81 maxTheta = std::max(maxTheta, toDeg(theta));
82 }
83 return std::make_pair(maxL - minL, maxTheta - minTheta);
84 }
Scalar theta() const
theta method
Acts::GeometryContext context() const
ToolHandle< MuonR4::ITrackSeedingDiagnosticsTool > m_seedingTool
The track seeding tool to construct the seed candidates and to estimate the initial parameters.
Eigen::Matrix< double, 2, 1 > Vector2D
MuonSegment_v1 MuonSegment
Reference the current persistent version:

◆ configAthHistogramming()

StatusCode AthHistogramming::configAthHistogramming ( const ServiceHandle< ITHistSvc > & histSvc,
const std::string & prefix,
const std::string & rootDir,
const std::string & histNamePrefix,
const std::string & histNamePostfix,
const std::string & histTitlePrefix,
const std::string & histTitlePostfix )
protectedinherited

To be called by the derived classes to fill the internal configuration.

Definition at line 51 of file AthHistogramming.cxx.

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

◆ declareGaudiProperty()

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

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

Definition at line 156 of file AthCommonDataStore.h.

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

◆ declareProperty()

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

Definition at line 145 of file AthCommonDataStore.h.

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

◆ detStore()

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

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

Definition at line 95 of file AthCommonDataStore.h.

◆ dumpLegacyTracks()

StatusCode MuonValR4::MsTrackTester::dumpLegacyTracks ( const EventContext & ctx)
private

Dumps the legacy containers to the TTree.

Definition at line 338 of file MsTrackTester.cxx.

338 {
339
340 //This for now is to be able to retrieve the matching between the legacy segments and tracks ...
341 const xAOD::MuonContainer* legacyMuons{nullptr};
342 ATH_CHECK(SG::get(legacyMuons, m_legacyMuonKey, ctx));
343
344 if (!legacyMuons) {
345 return StatusCode::SUCCESS;
346 }
347
348 for (const xAOD::Muon* muon : *legacyMuons) {
349 ATH_MSG_VERBOSE("iMuon " << muon->index() << " pT: "<<(muon->pt() *MeVtoGeV)<<" [GeV], eta: "
350 <<muon->eta() <<", phi: "<<toDeg(muon->phi())<<", q: "<<muon->charge()
351 <<", nSegments: "<<muon->nMuonSegments());
352 const xAOD::TrackParticle* track = muon->trackParticle(xAOD::Muon::TrackParticleType::MuonSpectrometerTrackParticle);
353 if (!track) {
354 continue;
355 }
356 m_summaryTool->copySummary(m_summaryTool->makeSummary(ctx, *track->track()), *track);
357 m_legacyTrks->push_back(track);
358 auto trkIdx = m_legacyTrks->find(track);
359 for (size_t s = 0; s < muon->nMuonSegments(); ++s) {
360 const xAOD::MuonSegment* segment = muon->muonSegment(s);
361 auto segIdx = m_legacyRecoSegs->push_back(*segment);
362 ATH_MSG_VERBOSE(std::format( "Legacy muon-segment link: segment {:} @{:}, eta: {:.2f}, phi {:.2f}",
363 printID(*segment), Amg::toString(segment->position()),
364 segment->direction().eta(), toDeg(segment->direction().phi())));
365 m_legacySegToTrkLinks[segIdx] = trkIdx;
366 }
367 }
368
369 //Dump also legacy segments
370 const xAOD::MuonSegmentContainer* legacyRecoSegs{nullptr};
371 ATH_CHECK(SG::get(legacyRecoSegs, m_legacySegmentKey, ctx));
372
373 if (legacyRecoSegs->size()) {
374 m_legacySegToTrkLinks[legacyRecoSegs->size()-1];
375 for (const xAOD::MuonSegment* seg : *legacyRecoSegs) {
376 //Store all segments
377 m_legacyRecoSegs->push_back(*seg);
378 }
379 }
380
381 const xAOD::TrackParticleContainer* legacyTrks{nullptr};
382 ATH_CHECK(SG::get(legacyTrks, m_legacyTrackKey, ctx));
383
384 for (const xAOD::TrackParticle* track : *legacyTrks) {
385 m_summaryTool->copySummary(m_summaryTool->makeSummary(ctx, *track->track()), *track);
386 m_legacyTrks->push_back(track);
387 }
388
389 return StatusCode::SUCCESS;
390 }
#define ATH_CHECK
Evaluate an expression and check for errors.
#define ATH_MSG_VERBOSE(x)
size_type size() const noexcept
Returns the number of elements in the collection.
ToolHandle< MuonR4::ITrackSummaryTool > m_summaryTool
Hit summary tool.
SegmentKey_t m_legacySegmentKey
Legacy segment container.
MuonVal::VectorBranch< unsigned short > & m_legacySegToTrkLinks
Link of the legacy track to the legacy segment.
ParticleBranchPtr_t m_legacyTrks
Output branches of the legacy MS tracks.
SegmentBranchPtr_t m_legacyRecoSegs
TrackKey_t m_legacyTrackKey
Legacy track reconstruction chain.
MuonKey_t m_legacyMuonKey
Legacy muons.
Amg::Vector3D direction() const
Returns the direction as Amg::Vector.
Amg::Vector3D position() const
Returns the position as Amg::Vector.
std::string toString(const Translation3D &translation, int precision=4)
GeoPrimitvesToStringConverter.
std::string printID(const xAOD::MuonSegment &seg)
Print the chamber ID of a segment, e.g.
const T * get(const ReadCondHandleKey< T > &key, const EventContext &ctx)
Convenience function to retrieve an object given a ReadCondHandleKey.
MuonSegmentContainer_v1 MuonSegmentContainer
Definition of the current "MuonSegment container version".
TrackParticle_v1 TrackParticle
Reference the current persistent version:
TrackParticleContainer_v1 TrackParticleContainer
Definition of the current "TrackParticle container version".
Muon_v1 Muon
Reference the current persistent version:
MuonContainer_v1 MuonContainer
Definition of the current "Muon container version".

◆ dumpRecoContent()

StatusCode MuonValR4::MsTrackTester::dumpRecoContent ( const EventContext & ctx)
private

Dump the reconstructed information.

Fetch the containers from store gate

Link the truth muons to the seeds

Link the seeds by the number of matched segments

Definition at line 428 of file MsTrackTester.cxx.

428 {
429 const xAOD::MuonContainer* muons{nullptr};
430 ATH_CHECK(SG::get(muons, m_muonKey, ctx));
432 const xAOD::MuonSegmentContainer* recoSegments{nullptr};
433 ATH_CHECK(SG::get(recoSegments, m_recoSegmentKey, ctx));
434
435 const MuonR4::MsTrackSeedContainer* trkSeeds{nullptr};
436 ATH_CHECK(SG::get(trkSeeds, m_msTrkSeedKey, ctx));
437
438 const ActsTrk::GeometryContext* gctx{nullptr};
439 ATH_CHECK(SG::get(gctx, m_geoCtxKey, ctx));
440
441
442 const AtlasFieldCacheCondObj* magCache{nullptr};
443 ATH_CHECK(SG::get(magCache, m_fieldCacheKey, ctx));
444
445 MagField::AtlasFieldCache magField{};
446 magCache->getInitializedCache(magField);
447 std::unordered_map<const xAOD::TruthParticle*,
448 std::vector<unsigned>> truthToSeedMatchCounter{};
449
450 if (!trkSeeds->empty()) {
451 m_seedTruthLink[trkSeeds->size() -1];
452 }
453
454 for (const MuonR4::MsTrackSeed& seed : *trkSeeds) {
455 unsigned int seedIdx = m_seedPos.size();
456 m_seedPos += seed.position();
457 m_seedType+= Acts::toUnderlying(seed.location());
458 m_seedSector += seed.sector().sector();
459 m_seedSummary->push_back(ctx, seed);
460
461 auto startPars = m_seedingTool->estimateStartParameters(ctx, seed);
462 if (startPars.ok()) {
463 m_seedDir += (*startPars).direction();
464 } else {
465 m_seedDir += Amg::Vector3D::UnitZ();
466 }
467 // m_seedDir
468 ATH_MSG_VERBOSE(" Dump new seed: "<<seed);
469 for (const xAOD::MuonSegment* seg : seed.segments()){
470 m_seedRecoSegMatch[seedIdx].push_back(m_recoSegs->push_back(*seg));
471 if (const xAOD::MuonSegment* truthSeg = MuonR4::getMatchedTruthSegment(*seg);
472 truthSeg != nullptr) {
473 std::vector<unsigned>& matchCounter = truthToSeedMatchCounter[MuonR4::getTruthMatchedParticle(*truthSeg)];
474 if (seedIdx >= matchCounter.size()) {
475 matchCounter.resize(seedIdx +1);
476 }
477 ++matchCounter[seedIdx];
478 }
479 }
480 const auto[seedLength, theta] = calcSeedLength(*gctx, seed);
481 m_seedLength+= seedLength;
483 m_seedQP += m_seedingTool->estimateQtimesP(gctx->context(), seed, magField) / Gaudi::Units::GeV;
484 m_seedGood += startPars.ok();
485 }
487 for (auto& [truthMuon, matches] : truthToSeedMatchCounter) {
488 const unsigned tIndex = m_truthTrks->find(truthMuon);
489 if (tIndex >= m_truthTrks->size()) {
490 continue;
491 }
493 while (std::count_if(matches.begin(), matches.end(),
494 [](const unsigned nMatched){
495 return nMatched > 0;
496 })) {
497 auto bestMatch = std::ranges::max_element(matches);
498 const std::size_t seedIdx = std::distance(matches.begin(), bestMatch);
499 m_truthMuToSeedIdx[tIndex].push_back(seedIdx);
500
501 m_seedTruthLink[seedIdx] = tIndex;
502 m_truthMuToSeedCounter[tIndex].push_back(*bestMatch);
503 (*bestMatch) = 0;
504 }
505 }
506
507 for (const xAOD::Muon* muon : *muons) {
508 m_muonTrks->push_back(muon);
509 }
510
511 for (const xAOD::MuonSegment* seg : *recoSegments) {
512 m_recoSegs->push_back(*seg);
513 }
514
515 return StatusCode::SUCCESS;
516 }
void getInitializedCache(MagField::AtlasFieldCache &cache) const
get B field cache for evaluation as a function of 2-d or 3-d position.
ParticleBranchPtr_t m_truthTrks
MuonVal::VectorBranch< int > & m_seedSector
Sector of the seed, even center, odd overlap regions, for details see:
MuonVal::VectorBranch< float > & m_seedThetaCone
Maximum angular difference between the segments part of the seed.
MuonVal::MatrixBranch< unsigned short > & m_truthMuToSeedIdx
Links to all MsTrkSeeds that could be matched to the truthMuon, i.e.
ActsTrk::GeoContextReadKey_t m_geoCtxKey
Dependency on the geometry alignment.
MuonVal::VectorBranch< char > & m_seedType
Is the seed in the encap or in the barrel chambers.
ParticleBranchPtr_t m_muonTrks
Stored muon information from the Acts muon reco chain.
SG::ReadCondHandleKey< AtlasFieldCacheCondObj > m_fieldCacheKey
Dependency on the magnetic field.
SegmentBranchPtr_t m_recoSegs
SG::ReadHandleKey< MuonR4::MsTrackSeedContainer > m_msTrkSeedKey
Temporary container write handle to push the seeds to store gate for later efficiency analysis.
MuonVal::MatrixBranch< unsigned short > & m_seedRecoSegMatch
Link of the track seed to the building segment.
MuonVal::VectorBranch< float > & m_seedQP
Estimated momentum times charge from the track seed.
MuonVal::MatrixBranch< unsigned short > & m_truthMuToSeedCounter
Number of matched segments in the seed.
MuonVal::ThreeVectorBranch m_seedPos
Simple seed information.
MuonVal::VectorBranch< float > & m_seedLength
Maximum separation between the segments on the reference plane.
MuonVal::UnitThreeVectorBranch m_seedDir
Seed direction vector.
SegmentKey_t m_recoSegmentKey
Primary segment container.
MuonVal::VectorBranch< char > & m_seedGood
Does the seeding tool construct valid parameters from the seed.
std::shared_ptr< TrackSummaryModule > m_seedSummary
Hit summary on the track seed.
MuonVal::VectorBranch< unsigned short > & m_seedTruthLink
Link to the truth muon.
std::pair< double, double > calcSeedLength(const ActsTrk::GeometryContext &gctx, const MuonR4::MsTrackSeed &seed) const
Calculate the length of the seed and the theta deflection angle The length is defined as the spread o...
MuonKey_t m_muonKey
Dependency on the R4 muon container.
const xAOD::TruthParticle * getTruthMatchedParticle(const xAOD::MuonSegment &segment)
Returns the particle truth-matched to the segment.
std::vector< MsTrackSeed > MsTrackSeedContainer
Definition MsTrackSeed.h:69
const xAOD::MuonSegment * getMatchedTruthSegment(const xAOD::MuonSegment &segment)
Returns the truth-matched segment.
TruthParticle_v1 TruthParticle
Typedef to implementation.

◆ dumpTruthContent()

StatusCode MuonValR4::MsTrackTester::dumpTruthContent ( const EventContext & ctx)
private

Dump truth information.

Allocate the memory

Definition at line 391 of file MsTrackTester.cxx.

391 {
392 if (!m_isMC) {
393 return StatusCode::SUCCESS;
394 }
395 const xAOD::MuonSegmentContainer* truthSegs{nullptr};
396 ATH_CHECK(SG::get(truthSegs, m_truthSegmentKey, ctx));
397
398
399 for (const xAOD::MuonSegment* seg : *truthSegs) {
400 ATH_MSG_VERBOSE("Dump truth segment "<<printID(*seg)<<" @"<<
401 Amg::toString(seg->position())<<", eta: "<<seg->direction().eta()
402 <<", phi: "<<toDeg(seg->direction().phi()));
403 m_truthSegs->push_back(*seg);
404 }
405 if (truthSegs->size()) {
406 m_truthSegToRecoLink[truthSegs->size()-1];
407 }
408
409
410 const xAOD::TruthParticleContainer* truthMuons{nullptr};
411 ATH_CHECK(SG::get(truthMuons, m_truthKey, ctx));
412 const std::size_t nT = truthMuons->size() - 1;
413 if (!truthMuons->empty()) {
418 }
419
420 for (const xAOD::TruthParticle* truth : *truthMuons) {
421 ATH_MSG_DEBUG("Truth muon: pT: "<<(truth->pt() *MeVtoGeV)
422 <<", eta: "<<truth->eta()<<", phi: "<<toDeg(truth->phi())<<", q: "<<truth->charge());
423 m_truthTrks->push_back(*truth);
424 }
425
426 return StatusCode::SUCCESS;
427 }
#define ATH_MSG_DEBUG(x)
bool empty() const noexcept
Returns true if the collection is empty.
SegmentBranchPtr_t m_truthSegs
MuonVal::MatrixBranch< unsigned short > & m_truthMuRecoSegLinks
Links from the truth muon to the segments.
SG::ReadHandleKey< xAOD::TruthParticleContainer > m_truthKey
Key to the truth particle collection.
Gaudi::Property< bool > m_isMC
MuonVal::MatrixBranch< unsigned short > & m_truthSegToRecoLink
Link of the truth segments to the matchin reco segments.
SegmentKey_t m_truthSegmentKey
Segment from the truth hits.
TruthParticleContainer_v1 TruthParticleContainer
Declare the latest version of the truth particle container.

◆ efficiency()

TEfficiency * AthHistogramming::efficiency ( const std::string & effName,
const std::string & tDir = "",
const std::string & stream = "" )
protectedinherited

Simplify the retrieval of registered TEfficiency.

Definition at line 211 of file AthHistogramming.cxx.

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

◆ evtStore()

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

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

Definition at line 85 of file AthCommonDataStore.h.

◆ execute()

StatusCode MuonValR4::MsTrackTester::execute ( const EventContext & ctx)
finaloverridevirtual

Execute method.

Implements AthAlgorithm.

Definition at line 517 of file MsTrackTester.cxx.

517 {
521
522 ATH_CHECK(m_tree.fill(ctx));
523 return StatusCode::SUCCESS;
524 }
StatusCode dumpLegacyTracks(const EventContext &ctx)
Dumps the legacy containers to the TTree.
StatusCode dumpRecoContent(const EventContext &ctx)
Dump the reconstructed information.
StatusCode dumpTruthContent(const EventContext &ctx)
Dump truth information.
MuonVal::MuonTesterTree m_tree

◆ 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 MuonValR4::MsTrackTester::finalize ( )
finaloverride

Definition at line 525 of file MsTrackTester.cxx.

525 {
526 ATH_CHECK(m_tree.write());
527 return StatusCode::SUCCESS;
528 }

◆ graph()

TGraph * AthHistogramming::graph ( const std::string & graphName,
const std::string & tDir = "",
const std::string & stream = "" )
protectedinherited

Simplify the retrieval of registered TGraphs.

Definition at line 453 of file AthHistogramming.cxx.

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

◆ hash()

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

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

Definition at line 430 of file AthHistogramming.h.

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

◆ hist()

TH1 * AthHistogramming::hist ( const std::string & histName,
const std::string & tDir = "",
const std::string & stream = "" )
protectedinherited

Simplify the retrieval of registered histograms of any type.

Definition at line 164 of file AthHistogramming.cxx.

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

◆ hist2d()

TH2 * AthHistogramming::hist2d ( const std::string & histName,
const std::string & tDir = "",
const std::string & stream = "" )
inlineprotectedinherited

Simplify the retrieval of registered 2-d histograms.

Definition at line 371 of file AthHistogramming.h.

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

◆ hist3d()

TH3 * AthHistogramming::hist3d ( const std::string & histName,
const std::string & tDir = "",
const std::string & stream = "" )
inlineprotectedinherited

Simplify the retrieval of registered 3-d histograms.

Definition at line 388 of file AthHistogramming.h.

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

◆ histSvc()

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

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

Definition at line 113 of file AthHistogramAlgorithm.h.

114{
115 return m_histSvc;
116}

◆ initialize()

StatusCode MuonValR4::MsTrackTester::initialize ( )
finaloverride

Link the truth segment to the reconstructed segment...

Link the truth segments to the truth partcle

Count the number of matched truth segments

Calculate the truth seed length

Link the reconstructed segments to the muon

Link the associated seed

Definition at line 85 of file MsTrackTester.cxx.

85 {
87 ATH_CHECK(m_truthKey.initialize(m_isMC));
88
89 ATH_CHECK(m_muonKey.initialize());
90 ATH_CHECK(m_msTrkSeedKey.initialize());
91 ATH_CHECK(m_recoSegmentKey.initialize());
92
93 ATH_CHECK(m_summaryTool.retrieve());
94 ATH_CHECK(m_seedingTool.retrieve());
95 ATH_CHECK(m_geoCtxKey.initialize());
96 ATH_CHECK(m_fieldCacheKey.initialize());
97
98 ATH_CHECK(m_legacyMuonKey.initialize(!m_legacyMuonKey.empty()));
99 ATH_CHECK(m_legacyTrackKey.initialize(!m_legacyMuonKey.empty()));
100 ATH_CHECK(m_legacySegmentKey.initialize(!m_legacyMuonKey.empty()));
101
102 ATH_CHECK(detStore()->retrieve(m_detMgr));
103
104 int evOpts{0};
105
106 m_recoSegs = std::make_unique<SegmentVariables>(m_tree, m_recoSegmentKey.key(), "Segments", msgLevel());
107 m_recoSegs->addVariable(std::make_unique<MuonVal::GenericAuxDecorationBranch<unsigned short>>(m_tree,
108 "Segments_passSeedQual",[this](const SG::AuxElement* aux){
109 const auto* seg = static_cast<const xAOD::MuonSegment*>(aux);
110 return m_segSelector->passSeedingQuality(Gaudi::Hive::currentContext(),
111 *seg); }));
112 m_recoSegs->addVariable(std::make_unique<MuonVal::GenericAuxDecorationBranch<unsigned short>>(m_tree,
113 "Segments_passTrackQual",[this](const SG::AuxElement* aux){
114 const auto* seg = static_cast<const xAOD::MuonSegment*>(aux);
115 return m_segSelector->passTrackQuality(Gaudi::Hive::currentContext(),
116 *seg);
117 }));
118
119
120 if (m_isMC) {
121 evOpts |= EventInfoBranch::isMC;
122 m_recoSegs->addVariable(std::make_unique<MuonVal::GenericAuxDecorationBranch<unsigned short>>(m_tree,
123 "Segments_truthSegLink",[this](const SG::AuxElement* aux){
124 const auto* seg = static_cast<const xAOD::MuonSegment*>(aux);
126 const unsigned linkIdx = truthS ? m_truthSegs->push_back(*truthS) : -1;
128 if (truthS) {
129 m_truthSegToRecoLink.push_back(linkIdx, m_recoSegs->push_back(*seg));
131 if (truthMuon){
132 m_truthTrks->push_back(truthMuon);
133 m_truthMuRecoSegLinks[m_truthTrks->find(truthMuon)].push_back(m_recoSegs->push_back(*seg));
134 }
135 }
136 return linkIdx;
137 }));
138
139 m_truthSegs = std::make_unique<SegmentVariables>(m_tree, m_truthSegmentKey.key(), "TruthSegments", msgLevel());
140
141 m_truthSegs->addVariable(std::make_unique<MuonVal::GenericAuxDecorationBranch<unsigned short>>(m_tree,
142 "TruthSegments_truthLink",[this](const SG::AuxElement* aux){
143 const auto* seg = static_cast<const xAOD::MuonSegment*>(aux);
145 m_truthTrks->push_back(truthP);
146 unsigned short linkIdx = m_truthTrks->find(truthP);
147 return linkIdx;
148 }));
149 for (auto loc : {Location::Barrel, Location::Endcap}) {
150 m_truthSegs->addVariable(std::make_unique<MuonVal::GenericAuxDecorationBranch<unsigned short>>(m_tree,
151 std::format("TruthSegments_has{}Proj", loc), [loc, this](const SG::AuxElement* aux) -> unsigned short {
152
153 const auto* seg = static_cast<const xAOD::MuonSegment*>(aux);
154 const ActsTrk::GeometryContext* gctx{};
155 (void) SG::get(gctx, m_geoCtxKey, Gaudi::Hive::currentContext()).isSuccess();
156 ExpandedSector sector{seg->position().phi()};
157 const Amg::Vector2D projPos{m_seedingTool->expressOnCylinder(gctx->context(), *seg, loc, sector)};
158 return m_seedingTool->withinBounds(projPos, loc);
159 }));
160 }
161
162 m_tree.addBranch(m_truthSegs);
163
164 m_truthTrks = std::make_unique<IParticleFourMomBranch>(m_tree, "TruthMuons");
165 m_truthTrks->addVariable<int>(-1, "truthOrigin");
166 m_truthTrks->addVariable<int>(-1, "truthType");
168 m_truthTrks->addVariable(
169 std::make_unique<GenericPartDecorBranch<xAOD::TruthParticle,
170 std::vector<unsigned short>>>(m_tree,
171 std::format("{:}_truthSegLinks", m_truthTrks->name()), [&] (const xAOD::TruthParticle& p){
172 std::vector<unsigned short> idx{};
173 for (const xAOD::MuonSegment* truthSeg: getTruthSegments(p)){
174 idx.push_back(m_truthSegs->push_back(*truthSeg));
175 }
176 return idx;
177 }));
179 m_truthTrks->addVariable(
180 std::make_unique<GenericPartDecorBranch<xAOD::TruthParticle, unsigned short>>(m_tree,
181 std::format("{:}_nTruthSegments", m_truthTrks->name()), [&] (const xAOD::TruthParticle& p) -> unsigned short {
182 return getTruthSegments(p).size();
183 }));
185 auto cone = std::make_shared<VectorBranch<float>>(m_tree,
186 std::format("{}_seedThetaCone", m_truthTrks->name()));
187 auto qTimesP = std::make_shared<VectorBranch<float>>(m_tree,
188 std::format("{}_qTimesP", m_truthTrks->name()));
189
190 m_tree.addBranch(cone);
191 m_tree.addBranch(qTimesP);
192 m_truthTrks->addVariable(
193 std::make_unique<GenericPartDecorBranch<xAOD::TruthParticle, float>>(m_tree,
194 std::format("{:}_seedLength", m_truthTrks->name()), [cone, qTimesP, this] (const xAOD::TruthParticle& p) -> float {
195 const ActsTrk::GeometryContext* gctx{};
196 const AtlasFieldCacheCondObj* magCache{nullptr};
197 const EventContext& ctx{Gaudi::Hive::currentContext()};
198 (void) SG::get(gctx, m_geoCtxKey, ctx).isSuccess();
199 (void) SG::get(magCache, m_fieldCacheKey, ctx).isSuccess();
200 MagField::AtlasFieldCache magField{};
201 magCache->getInitializedCache(magField);
202 auto truthSeed = makeSeedFromTruth(*gctx, p);
203 if (!truthSeed) {
204 cone->push_back(-1);
205 qTimesP->push_back(0);
206 return -1.;
207 }
208 auto [length, theta] = calcSeedLength(*gctx, *truthSeed);
209 cone->push_back(theta);
210 qTimesP->push_back(m_seedingTool->estimateQtimesP(gctx->context(), *truthSeed, magField) / Gaudi::Units::GeV);
211
212 return length;
213 }));
214
215 m_tree.addBranch(m_truthTrks);
216 m_trkTruthLinks.emplace_back(m_truthSegmentKey, "truthParticleLink");
217 m_trkTruthLinks.emplace_back(m_truthKey, "truthSegmentLinks");
218 m_trkTruthLinks.emplace_back(m_recoSegmentKey, "truthSegmentLink");
219 m_trkTruthLinks.emplace_back(m_muonKey, "truthParticleLink");
220 }
221
222 m_tree.addBranch(m_recoSegs);
223 m_tree.addBranch(std::make_unique<EventInfoBranch>(m_tree, evOpts));
224
225 static const std::vector<std::string> trackSummaries{
226 // Inner
227 "innerSmallHits", "innerLargeHits", "innerSmallHoles", "innerLargeHoles",
228 "innerClosePrecisionHits",
229 // Middle
230 "middleSmallHits", "middleLargeHits", "middleSmallHoles",
231 "middleLargeHoles", "middleClosePrecisionHits",
232 // Outer
233 "outerSmallHits", "outerLargeHits", "outerSmallHoles", "outerLargeHoles",
234 "outerClosePrecisionHits",
235 // Extended
236 "extendedSmallHits", "extendedLargeHits", "extendedSmallHoles",
237 "extendedLargeHoles", "extendedClosePrecisionHits",
238 "innerTriggerEtaHits", "innerTriggerPhiHits",
239 "middleTriggerEtaHits", "middleTriggerPhiHits",
240 "outerTriggerEtaHits", "outerTriggerPhiHits",
241
242 "innerTriggerEtaHoles", "innerTriggerPhiHoles",
243 "middleTriggerEtaHoles", "middleTriggerPhiHoles",
244 "outerTriggerEtaHoles", "outerTriggerPhiHoles",
245 };
246 if(!m_legacyTrackKey.empty()) {
247 m_legacyTrks = std::make_unique<IParticleFourMomBranch>(m_tree, "LegacyMSTrks");
248 m_legacyTrks->addVariable(std::make_unique<TrackChi2Branch>(*m_legacyTrks));
249 if (m_isMC) {
251 return xAOD::TruthHelpers::getTruthParticle(*trk); }, "truth", "LegacyMS");
252 }
253 for (const auto& summary : trackSummaries) {
254 m_legacyTrks->addVariable<uint8_t>(-1, summary);
255 }
256 m_tree.addBranch(m_legacyTrks);
257 m_legacyRecoSegs = std::make_unique<SegmentVariables>(m_tree, m_legacySegmentKey.key(), "LegacyRecoSegments", msgLevel());
258 m_tree.addBranch(m_legacyRecoSegs);
259 }
260
261 m_seedSummary = std::make_shared<TrackSummaryModule>(m_tree, "MsTrkSeed", m_summaryTool.get());
262 m_muonTrks = std::make_shared<IParticleFourMomBranch>(m_tree, "ActsMuons");
263 m_muonTrks->addVariable(std::make_unique<TrackChi2Branch>(*m_muonTrks));
264 m_muonTrks->addVariable(std::make_unique<TrackFitIterBranch>(*m_muonTrks));
265
266 using TrkType = xAOD::Muon::TrackParticleType;
267 auto dumpTrack = [&](const std::string& trkName,
268 TrkType type) {
269 auto trkColl = std::make_shared<IParticleFourMomBranch>(m_tree, std::format("Acts{:}", trkName));
270 trkColl->addVariable(std::make_unique<TrackChi2Branch>(*trkColl));
271 trkColl->addVariable(std::make_unique<TrackFitIterBranch>(*trkColl));
272 trkColl->addVariable(std::make_unique<MaterialRecorderBranch>(*trkColl));
273 trkColl->addVariable(std::make_unique<EnergyLossBranch>(*trkColl));
274
275 trkColl->addVariable<float>("d0");
276 trkColl->addVariable<float>("z0");
277 for (const auto& summary : trackSummaries) {
278 trkColl->addVariable<uint8_t>(-1, summary);
279 }
280 m_muonTrks->addVariable(std::make_unique<MuonVal::LinkerBranch>(*m_muonTrks, trkColl,
281 [type](const xAOD::IParticle* muonP) -> const xAOD::IParticle* {
282 const auto* muon = dynamic_cast<const xAOD::Muon*>(muonP);
283 if (!muon) {
284 return nullptr;
285 }
286 return muon->trackParticle(type);
287 }, trkName));
288 };
289
290 dumpTrack("MsTrk", TrkType::MuonSpectrometerTrackParticle);
291 dumpTrack("MeTrk", TrkType::ExtrapolatedMuonSpectrometerTrackParticle);
292
293 m_muonTrks->addVariable<uint16_t>("allAuthors");
294 m_muonTrks->addVariable<uint16_t>("author");
296 m_muonTrks->addVariable(std::make_unique<GenericPartDecorBranch<xAOD::Muon,
297 std::vector<unsigned short>>>(m_tree,
298 std::format("{:}_segmentLinks", m_muonTrks->name()), [&] (const xAOD::Muon& p){
299 std::vector<unsigned short> idx{};
300 for (unsigned seg = 0 ; seg < p.nMuonSegments(); ++seg) {
301 idx.push_back(m_recoSegs->push_back(*p.muonSegment(seg)));
302 }
303 return idx;
304 }));
306 m_muonTrks->addVariable(std::make_unique<GenericPartDecorBranch<xAOD::Muon, unsigned short>>(m_tree,
307 std::format("{:}_seedLink", m_muonTrks->name()), [&] (const xAOD::Muon& p) -> unsigned short {
308 using enum xAOD::Muon::TrackParticleType;
309 const xAOD::TrackParticle* msTrack = p.trackParticle(MuonSpectrometerTrackParticle);
310 if (!msTrack) {
311 return -1;
312 }
313 auto actsTrk = ActsTrk::getActsTrack(*msTrack);
314 if (!actsTrk) {
315 THROW_EXCEPTION("Cannot find the associated ms track from the primary track");
316 }
317 return actsTrk->component<std::size_t>("parentSeed");
318 }));
319
320 if (m_isMC) {
322 [](const xAOD::IParticle* trk) -> const xAOD::TruthParticle* {
324 }, "truth", "ActsMuon");
325 }
326
327 for (const auto& summary : trackSummaries) {
328 m_muonTrks->addVariable<uint8_t>(summary);
329 }
330 m_tree.addBranch(m_muonTrks);
331
332 m_tree.addBranch(m_seedSummary);
333
334 ATH_CHECK(m_trkTruthLinks.initialize());
335 ATH_CHECK(m_tree.init(this));
336 return StatusCode::SUCCESS;
337 }
double length(const pvec &v)
if(pathvar)
const ServiceHandle< StoreGateSvc > & detStore() const
const MuonGMR4::MuonDetectorManager * m_detMgr
std::optional< MuonR4::MsTrackSeed > makeSeedFromTruth(const ActsTrk::GeometryContext &gctx, const xAOD::TruthParticle &truthMuon) const
Construct MS track seed from the truth associated segments.
SG::ReadDecorHandleKeyArray< SG::AuxVectorBase > m_trkTruthLinks
Decoration dependency to the MS truth track links.
ToolHandle< MuonR4::ISegmentSelectionTool > m_segSelector
Selection tool to quantify the segment candidate quality.
static bool connectCollections(ParticleBranch_ptr primColl, ParticleBranch_ptr secondColl, Linker_t fromPrimToSec, const std::string &altPrimName="", const std::string &altSecName="")
@ isMC
Flag determining whether the branch is simulation.
std::optional< ActsTrk::TrackContainer::ConstTrackProxy > getActsTrack(const xAOD::TrackParticle &trkPart)
Return the proxy to the Acts track from which the track particle was made frome.
Definition Decoration.cxx:9
for(size_t i=0;i< m_blockFillers.size();i++)
Fill one block.
std::vector< const xAOD::MuonSegment * > getTruthSegments(const xAOD::TruthParticle &truthMuon)
Returns the segments associated to the truth muon.
AuxElement(SG::AuxVectorData *container, size_t index)
Base class for elements of a container that can have aux data.
const xAOD::TruthParticle * getTruthParticle(const xAOD::IParticle &p)
Return the truthParticle associated to the given IParticle (if any).
setWord1 uint16_t
#define THROW_EXCEPTION(MESSAGE)
Definition throwExcept.h:10

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

◆ makeSeedFromTruth()

std::optional< MsTrackSeed > MuonValR4::MsTrackTester::makeSeedFromTruth ( const ActsTrk::GeometryContext & gctx,
const xAOD::TruthParticle & truthMuon ) const
private

Construct MS track seed from the truth associated segments.

A nullopt is returned if either no segment is matched to the particle or no valid seed could be constructed

Parameters
gctxThe geometry context to align the segment within ATLAS
truthMuonReference to the truth muon for which a seed should be constructed

Construct the 2 candidate seeds

Definition at line 35 of file MsTrackTester.cxx.

36 {
37 std::vector<const xAOD::MuonSegment*> matchedSegs = MuonR4::getTruthSegments(truthMuon);
38 if (matchedSegs.empty() || toLayerIndex(matchedSegs.front()->chamberIndex()) == toLayerIndex(matchedSegs.back()->chamberIndex())) {
39 return std::nullopt;
40 }
41 ExpandedSector sector{matchedSegs[0]->position().phi()};
43 MsTrackSeed barrelSeed{Location::Barrel, sector};
44 MsTrackSeed endcapSeed{Location::Endcap, sector};
45 for (const xAOD::MuonSegment* seg : matchedSegs) {
46 barrelSeed.addSegment(seg);
47 endcapSeed.addSegment(seg);
48 }
49 barrelSeed.setPosition(matchedSegs[0]->position());
50 endcapSeed.setPosition(matchedSegs[0]->position());
51 const auto [barrelLength, barrelTheta] = calcSeedLength(gctx, barrelSeed);
52 const auto [endcapLength, endcapTheta] = calcSeedLength(gctx, endcapSeed);
53 ATH_MSG_VERBOSE(__func__<<"() "<<__LINE__<<" - Constructed new seed from truth muon wih pT:"
54 <<(truthMuon.pt()/ Gaudi::Units::GeV)<<" [GeV], eta: "<<truthMuon.eta()
55 <<", phi: "<<toDeg(truthMuon.phi())<<", q: "<<truthMuon.charge()
56 <<", matchedSeg: "<<matchedSegs.size()<< " barrel (L/theta): "<<barrelLength
57 <<"/"<<barrelTheta<<" - endcap (L/theta): "
58 <<endcapLength<<"/"<<endcapTheta<<"\n"<<barrelSeed);
59 if (barrelLength < 0 && endcapLength < 0) {
60 ATH_MSG_WARNING(__func__<<"() "<<__LINE__<<" - Invalid seed");
61 return std::nullopt;
62 }
63 return barrelLength < 0 || std::abs(endcapLength) < barrelLength
64 ? endcapSeed : barrelSeed;
65 }
#define ATH_MSG_WARNING(x)
void addSegment(const xAOD::MuonSegment *seg)
Append a segment to the seed.
void setPosition(Amg::Vector3D &&pos)
set the seed's position
LayerIndex toLayerIndex(ChIndex index)
convert ChIndex into LayerIndex
const Amg::Vector3D & position() const
Method to retrieve the position of the Intersection.

◆ msg()

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

Definition at line 24 of file AthCommonMsg.h.

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

◆ msgLvl()

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

Definition at line 30 of file AthCommonMsg.h.

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

◆ myReplace()

void AthHistogramming::myReplace ( std::string & str,
const std::string & oldStr,
const std::string & newStr )
privateinherited

Helper method to replace sub-string.

Definition at line 551 of file AthHistogramming.cxx.

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

◆ outputHandles()

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

Return this algorithm's output handles.

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

◆ renounce()

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

Definition at line 380 of file AthCommonDataStore.h.

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

◆ renounceArray()

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

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364 {
366 }

◆ setFilterPassed()

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

Set filter decision:

Reimplemented in AthFilterAlgorithm.

Definition at line 99 of file AthCommonAlgorithm.h.

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

◆ sysExecute()

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

Execute an algorithm.

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

Reimplemented in AthAnalysisAlgorithm.

Definition at line 80 of file AthCommonAlgorithm.cxx.

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

◆ sysInitialize()

StatusCode AthHistogramAlgorithm::sysInitialize ( )
virtualinherited

Initialization method invoked by the framework.

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

Reimplemented from AthCommonAlgorithm< Gaudi::Algorithm >.

Reimplemented in AthAnalysisAlgorithm.

Definition at line 75 of file AthHistogramAlgorithm.cxx.

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

◆ sysStart()

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

Handle START transition.

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

◆ tree()

TTree * AthHistogramming::tree ( const std::string & treeName,
const std::string & tDir = "",
const std::string & stream = "" )
protectedinherited

Simplify the retrieval of registered TTrees.

Definition at line 339 of file AthHistogramming.cxx.

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

◆ updateVHKA()

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

Definition at line 308 of file AthCommonDataStore.h.

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

Member Data Documentation

◆ m_detMgr

const MuonGMR4::MuonDetectorManager* MuonValR4::MsTrackTester::m_detMgr {nullptr}
private

Definition at line 68 of file MsTrackTester.h.

68{nullptr};

◆ m_detStore

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

Pointer to StoreGate (detector store by default).

Definition at line 393 of file AthCommonDataStore.h.

◆ m_effMap

EffMap_t AthHistogramming::m_effMap
privateinherited

The map of histogram names to their pointers.

Definition at line 210 of file AthHistogramming.h.

◆ m_evtStore

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

Pointer to StoreGate (event store by default).

Definition at line 390 of file AthCommonDataStore.h.

◆ m_extendedExtraObjects

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

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

Empty if no symlinks were found.

Definition at line 108 of file AthCommonAlgorithm.h.

◆ m_fieldCacheKey

SG::ReadCondHandleKey<AtlasFieldCacheCondObj> MuonValR4::MsTrackTester::m_fieldCacheKey {this, "MagFieldKey", "fieldCondObj", "Name of the Magnetic Field conditions object key"}
private

Dependency on the magnetic field.

Definition at line 105 of file MsTrackTester.h.

105{this, "MagFieldKey", "fieldCondObj", "Name of the Magnetic Field conditions object key"};

◆ m_geoCtxKey

ActsTrk::GeoContextReadKey_t MuonValR4::MsTrackTester::m_geoCtxKey {this, "AlignmentKey", "ActsAlignment", "cond handle key"}
private

Dependency on the geometry alignment.

Definition at line 93 of file MsTrackTester.h.

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

◆ m_graphMap

GraphMap_t AthHistogramming::m_graphMap
privateinherited

The map of TGraph names to their pointers.

Definition at line 224 of file AthHistogramming.h.

◆ m_histMap

HistMap_t AthHistogramming::m_histMap
privateinherited

The map of histogram names to their pointers.

Definition at line 203 of file AthHistogramming.h.

◆ m_histNamePostfix

std::string AthHistogramAlgorithm::m_histNamePostfix
privateinherited

The postfix for the histogram THx name.

Definition at line 97 of file AthHistogramAlgorithm.h.

◆ m_histNamePrefix

std::string AthHistogramAlgorithm::m_histNamePrefix
privateinherited

The prefix for the histogram THx name.

Definition at line 94 of file AthHistogramAlgorithm.h.

◆ m_histSvc

ServiceHandle<ITHistSvc> AthHistogramAlgorithm::m_histSvc
privateinherited

Default constructor: AthHistogramAlgorithm();.

a handle on the Hist/TTree registration service

Definition at line 83 of file AthHistogramAlgorithm.h.

◆ m_histTitlePostfix

std::string AthHistogramAlgorithm::m_histTitlePostfix
privateinherited

The postfix for the histogram THx title.

Definition at line 103 of file AthHistogramAlgorithm.h.

◆ m_histTitlePrefix

std::string AthHistogramAlgorithm::m_histTitlePrefix
privateinherited

The prefix for the histogram THx title.

Definition at line 100 of file AthHistogramAlgorithm.h.

◆ m_isMC

Gaudi::Property<bool> MuonValR4::MsTrackTester::m_isMC {this, "isMC", false}
private

Definition at line 70 of file MsTrackTester.h.

70{this, "isMC", false};

◆ m_legacyMuonKey

MuonKey_t MuonValR4::MsTrackTester::m_legacyMuonKey {this,"LegacyMuonKey", "Muons"}
private

Legacy muons.

Definition at line 109 of file MsTrackTester.h.

109{this,"LegacyMuonKey", "Muons"};

◆ m_legacyRecoSegs

SegmentBranchPtr_t MuonValR4::MsTrackTester::m_legacyRecoSegs {}
private

Definition at line 119 of file MsTrackTester.h.

119{};

◆ m_legacySegmentKey

SegmentKey_t MuonValR4::MsTrackTester::m_legacySegmentKey {this, "LegacySegmentKey", "MuonSegments"}
private

Legacy segment container.

Definition at line 85 of file MsTrackTester.h.

85{this, "LegacySegmentKey", "MuonSegments"};

◆ m_legacySegToTrkLinks

MuonVal::VectorBranch<unsigned short>& MuonValR4::MsTrackTester::m_legacySegToTrkLinks {m_tree.newVector<unsigned short>("LegacyRecoSegments_trkLinks", -1)}
private

Link of the legacy track to the legacy segment.

Definition at line 147 of file MsTrackTester.h.

147{m_tree.newVector<unsigned short>("LegacyRecoSegments_trkLinks", -1)};

◆ m_legacyTrackKey

TrackKey_t MuonValR4::MsTrackTester::m_legacyTrackKey {this,"LegacyTrackKey", "MuonSpectrometerTrackParticles"}
private

Legacy track reconstruction chain.

Definition at line 103 of file MsTrackTester.h.

103{this,"LegacyTrackKey", "MuonSpectrometerTrackParticles"};

◆ m_legacyTrks

ParticleBranchPtr_t MuonValR4::MsTrackTester::m_legacyTrks {}
private

Output branches of the legacy MS tracks.

Definition at line 157 of file MsTrackTester.h.

157{};

◆ m_msg

MsgStream AthHistogramming::m_msg
privateinherited

Cached Message Stream.

Definition at line 250 of file AthHistogramming.h.

◆ m_msTrkSeedKey

SG::ReadHandleKey<MuonR4::MsTrackSeedContainer> MuonValR4::MsTrackTester::m_msTrkSeedKey {this, "MsTrkSeedKey", "MsTrackSeeds"}
private

Temporary container write handle to push the seeds to store gate for later efficiency analysis.

Definition at line 91 of file MsTrackTester.h.

91{this, "MsTrkSeedKey", "MsTrackSeeds"};

◆ m_muonKey

MuonKey_t MuonValR4::MsTrackTester::m_muonKey {this, "MuonKey", "MuonsR4"}
private

Dependency on the R4 muon container.

Definition at line 95 of file MsTrackTester.h.

95{this, "MuonKey", "MuonsR4"};

◆ m_muonTrks

ParticleBranchPtr_t MuonValR4::MsTrackTester::m_muonTrks {}
private

Stored muon information from the Acts muon reco chain.

Definition at line 114 of file MsTrackTester.h.

114{};

◆ m_name

std::string AthHistogramming::m_name
privateinherited

Instance name.

Definition at line 247 of file AthHistogramming.h.

◆ m_prefix

std::string AthHistogramAlgorithm::m_prefix
privateinherited

Name of the ROOT output stream (file).

Definition at line 88 of file AthHistogramAlgorithm.h.

◆ m_recoSegmentKey

SegmentKey_t MuonValR4::MsTrackTester::m_recoSegmentKey {this, "SegmentKey", "MuonSegmentsFromR4"}
private

Primary segment container.

Definition at line 83 of file MsTrackTester.h.

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

◆ m_recoSegs

SegmentBranchPtr_t MuonValR4::MsTrackTester::m_recoSegs {}
private

Definition at line 118 of file MsTrackTester.h.

118{};

◆ m_rootDir

std::string AthHistogramAlgorithm::m_rootDir
privateinherited

Name of the ROOT directory.

Definition at line 91 of file AthHistogramAlgorithm.h.

◆ m_seedDir

MuonVal::UnitThreeVectorBranch MuonValR4::MsTrackTester::m_seedDir {m_tree, "MsTrkSeed_direction"}
private

Seed direction vector.

Definition at line 124 of file MsTrackTester.h.

124{m_tree, "MsTrkSeed_direction"};

◆ m_seedGood

MuonVal::VectorBranch<char>& MuonValR4::MsTrackTester::m_seedGood {m_tree.newVector<char>("MsTrkSeed_goodSeed")}
private

Does the seeding tool construct valid parameters from the seed.

Definition at line 134 of file MsTrackTester.h.

134{m_tree.newVector<char>("MsTrkSeed_goodSeed")};

◆ m_seedingTool

ToolHandle<MuonR4::ITrackSeedingDiagnosticsTool> MuonValR4::MsTrackTester::m_seedingTool {this, "SeedingTool", ""}
private

The track seeding tool to construct the seed candidates and to estimate the initial parameters.

Definition at line 99 of file MsTrackTester.h.

99{this, "SeedingTool", ""};

◆ m_seedLength

MuonVal::VectorBranch<float>& MuonValR4::MsTrackTester::m_seedLength {m_tree.newVector<float>("MsTrkSeed_length")}
private

Maximum separation between the segments on the reference plane.

Definition at line 130 of file MsTrackTester.h.

130{m_tree.newVector<float>("MsTrkSeed_length")};

◆ m_seedPos

MuonVal::ThreeVectorBranch MuonValR4::MsTrackTester::m_seedPos {m_tree, "MsTrkSeed_position"}
private

Simple seed information.

Definition at line 122 of file MsTrackTester.h.

122{m_tree, "MsTrkSeed_position"};

◆ m_seedQP

MuonVal::VectorBranch<float>& MuonValR4::MsTrackTester::m_seedQP {m_tree.newVector<float>("MsTrkSeed_qTimesP")}
private

Estimated momentum times charge from the track seed.

Definition at line 136 of file MsTrackTester.h.

136{m_tree.newVector<float>("MsTrkSeed_qTimesP")};

◆ m_seedRecoSegMatch

MuonVal::MatrixBranch<unsigned short>& MuonValR4::MsTrackTester::m_seedRecoSegMatch {m_tree.newMatrix<unsigned short>("MsTrkSeed_segmentLinks")}
private

Link of the track seed to the building segment.

Definition at line 142 of file MsTrackTester.h.

142{m_tree.newMatrix<unsigned short>("MsTrkSeed_segmentLinks")};

◆ m_seedSector

MuonVal::VectorBranch<int>& MuonValR4::MsTrackTester::m_seedSector {m_tree.newVector<int>("MsTrkSeed_sector")}
private

Sector of the seed, even center, odd overlap regions, for details see:

Definition at line 128 of file MsTrackTester.h.

128{m_tree.newVector<int>("MsTrkSeed_sector")};

◆ m_seedSummary

std::shared_ptr<TrackSummaryModule> MuonValR4::MsTrackTester::m_seedSummary {}
private

Hit summary on the track seed.

Definition at line 140 of file MsTrackTester.h.

140{};

◆ m_seedThetaCone

MuonVal::VectorBranch<float>& MuonValR4::MsTrackTester::m_seedThetaCone {m_tree.newVector<float>("MsTrkSeed_thetaCone")}
private

Maximum angular difference between the segments part of the seed.

Definition at line 132 of file MsTrackTester.h.

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

◆ m_seedTruthLink

MuonVal::VectorBranch<unsigned short>& MuonValR4::MsTrackTester::m_seedTruthLink {m_tree.newVector<unsigned short>("MsTrkSeed_truthLink", -1)}
private

Link to the truth muon.

Definition at line 138 of file MsTrackTester.h.

138{m_tree.newVector<unsigned short>("MsTrkSeed_truthLink", -1)};

◆ m_seedType

MuonVal::VectorBranch<char>& MuonValR4::MsTrackTester::m_seedType {m_tree.newVector<char>("MsTrkSeed_type")}
private

Is the seed in the encap or in the barrel chambers.

Definition at line 126 of file MsTrackTester.h.

126{m_tree.newVector<char>("MsTrkSeed_type")};

◆ m_segSelector

ToolHandle<MuonR4::ISegmentSelectionTool> MuonValR4::MsTrackTester::m_segSelector {this, "SegmentSelectionTool" , "" }
private

Selection tool to quantify the segment candidate quality.

Definition at line 101 of file MsTrackTester.h.

101{this, "SegmentSelectionTool" , "" };

◆ m_streamName

std::string AthHistogramming::m_streamName
privateinherited

Name of the ROOT output stream (file).

Definition at line 228 of file AthHistogramming.h.

◆ m_summaryTool

ToolHandle<MuonR4::ITrackSummaryTool> MuonValR4::MsTrackTester::m_summaryTool {this, "SummaryTool" ,""}
private

Hit summary tool.

Definition at line 97 of file MsTrackTester.h.

97{this, "SummaryTool" ,""};

◆ m_tree

MuonVal::MuonTesterTree MuonValR4::MsTrackTester::m_tree {"MsTrackValidTest", "MuonTrackTester"}
private

Definition at line 66 of file MsTrackTester.h.

66{"MsTrackValidTest", "MuonTrackTester"};

◆ m_treeMap

TreeMap_t AthHistogramming::m_treeMap
privateinherited

The map of TTree names to their pointers.

Definition at line 217 of file AthHistogramming.h.

◆ m_trkTruthLinks

SG::ReadDecorHandleKeyArray<SG::AuxVectorBase> MuonValR4::MsTrackTester::m_trkTruthLinks {this, "TruthTrackLinks", {}}
private

Decoration dependency to the MS truth track links.

Definition at line 89 of file MsTrackTester.h.

89{this, "TruthTrackLinks", {}};

◆ m_truthKey

SG::ReadHandleKey<xAOD::TruthParticleContainer> MuonValR4::MsTrackTester::m_truthKey {this, "TruthKey", "MuonTruthParticles"}
private

Key to the truth particle collection.

Definition at line 87 of file MsTrackTester.h.

87{this, "TruthKey", "MuonTruthParticles"};

◆ m_truthMuRecoSegLinks

MuonVal::MatrixBranch<unsigned short>& MuonValR4::MsTrackTester::m_truthMuRecoSegLinks {m_tree.newMatrix<unsigned short>("TruthMuons_recoSegLinks")}
private

Links from the truth muon to the segments.

Definition at line 154 of file MsTrackTester.h.

154{m_tree.newMatrix<unsigned short>("TruthMuons_recoSegLinks")};

◆ m_truthMuToSeedCounter

MuonVal::MatrixBranch<unsigned short>& MuonValR4::MsTrackTester::m_truthMuToSeedCounter {m_tree.newMatrix<unsigned short>("TruthMuons_seedNSeg")}
private

Number of matched segments in the seed.

Definition at line 152 of file MsTrackTester.h.

152{m_tree.newMatrix<unsigned short>("TruthMuons_seedNSeg")};

◆ m_truthMuToSeedIdx

MuonVal::MatrixBranch<unsigned short>& MuonValR4::MsTrackTester::m_truthMuToSeedIdx {m_tree.newMatrix<unsigned short>("TruthMuons_seedLinks", -1)}
private

Links to all MsTrkSeeds that could be matched to the truthMuon, i.e.

>= 1 segment

Definition at line 150 of file MsTrackTester.h.

150{m_tree.newMatrix<unsigned short>("TruthMuons_seedLinks", -1)};

◆ m_truthSegmentKey

SegmentKey_t MuonValR4::MsTrackTester::m_truthSegmentKey {this, "TruthSegmentKey", "MuonTruthSegments"}
private

Segment from the truth hits.

Definition at line 81 of file MsTrackTester.h.

81{this, "TruthSegmentKey", "MuonTruthSegments"};

◆ m_truthSegs

SegmentBranchPtr_t MuonValR4::MsTrackTester::m_truthSegs {}
private

Definition at line 117 of file MsTrackTester.h.

117{};

◆ m_truthSegToRecoLink

MuonVal::MatrixBranch<unsigned short>& MuonValR4::MsTrackTester::m_truthSegToRecoLink {m_tree.newMatrix<unsigned short>("TruthSegments_recoSegLinks",-1)}
private

Link of the truth segments to the matchin reco segments.

Definition at line 144 of file MsTrackTester.h.

144{m_tree.newMatrix<unsigned short>("TruthSegments_recoSegLinks",-1)};

◆ m_truthTrks

ParticleBranchPtr_t MuonValR4::MsTrackTester::m_truthTrks {}
private

Definition at line 112 of file MsTrackTester.h.

112{};

◆ m_varHandleArraysDeclared

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

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

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

Definition at line 398 of file AthCommonDataStore.h.


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