MuonValR4::MsTrackTester Class Reference

#include <MsTrackTester.h>

Inheritance diagram for MuonValR4::MsTrackTester:

Classes
struct	SegmentTagVariables
	Searched the MuTagIMO tag. More...

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	MuonTagKey_t = SG::ReadHandleKey<MuonR4::MuonTagContainer>
	Abrivate the muon tag container.
using	ParticleBranchPtr_t = std::shared_ptr<MuonVal::IParticleFourMomBranch>
using	SegmentBranchPtr_t = std::shared_ptr<MuonPRDTest::SegmentVariables>
	Abrivation of the branches containing sement information.
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t
typedef uint32_t	hash_t
	typedef for the internal hash
typedef std::map< const hash_t, TH1 * >	HistMap_t
	Typedef for convenience.
typedef std::map< const hash_t, TEfficiency * >	EffMap_t
	Typedef for convenience.
typedef std::map< const hash_t, TTree * >	TreeMap_t
	Typedef for convenience.
typedef std::map< const hash_t, TGraph * >	GraphMap_t
	Typedef for convenience.

Private Member Functions
std::optional< MuonR4::MsTrackSeed >	makeSeedFromTruth (const Acts::GeometryContext &tgContext, const xAOD::TruthParticle &truthMuon) const
	Construct MS track seed from the truth associated segments.
std::pair< double, double >	calcSeedLength (const Acts::GeometryContext &tgContext, 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.
const xAOD::TruthParticle *	truthTreeParticle (const xAOD::IParticle &part) const
	Returns the matched truth particle that is dumped in the tree.
std::vector< const xAOD::MuonSegment * >	getAssociatedSegments (const xAOD::TrackParticle &idTrack, const EventContext &ctx) const
	Returns the collection of reconstructed segments that can be associated to the id track If the track is connected with a truth muon, the list of reconstructed segments matched to the truth muon is returned.
const MuonR4::MuonTag *	findBaseIdTag (const xAOD::TrackParticle &idTrack, const EventContext &ctx) const
	Searches the Id track in the basline collection of all ID tracks selected for the combined muon reconstructiion (STACO/MuidCo/MuTagIMO/Calo).
const MuonR4::MuonTag *	findMuTagIMO (const xAOD::TrackParticle &idTrack, const EventContext &ctx) const
const MuonGMR4::SpectrometerSector *	getEnvelope (const xAOD::MuonSegment &segment) const
SegmentTagVariables	calcSegTagVariables (const xAOD::TrackParticle &idTrack, const EventContext &ctx)
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.
const Acts::Logger &	logger () const
	Return the reference to the Acts logger.
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"}
	The output tree object.
const MuonGMR4::MuonDetectorManager *	m_detMgr {nullptr}
	The detector manager to fetch the segment surfaces.
Gaudi::Property< bool >	m_isMC {this, "isMC", false}
	Toggle whether the job runs on MC or not.
Gaudi::Property< bool >	m_storeID {this, "storeIdTrks", true}
	Toggle whether to process ID tracks or not.
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::TruthParticleContainer >	m_truthKey {this, "TruthKey", "MuonTruthParticles"}
	Key to the truth particle collection.
SG::ReadDecorHandleKeyArray< SG::AuxVectorBase >	m_trkTruthLinks {this, "TruthTrackLinks", {}}
	Decoration dependency to the MS truth track links.
SG::ReadHandleKey< MuonR4::MsTrackSeedContainer >	m_msTrkSeedKey {this, "MsTrkSeedKey", "MsTrackSeeds"}
	Temporary container write handle to push the seeds to store gate for later efficiency analysis.
ActsTrk::ContextUtility	m_ctxProvider {this}
	Context provider for geometry, magnetic field and calibration contexts.
MuonKey_t	m_muonKey {this, "MuonKey", "MuonsR4"}
	Dependency on the R4 muon container.
ToolHandle< MuonR4::ITrackSummaryTool >	m_summaryTool {this, "SummaryTool" ,""}
	Hit summary tool.
ToolHandle< MuonR4::ITrackSeedingDiagnosticsTool >	m_seedingTool {this, "SeedingTool", ""}
	The track seeding tool to construct the seed candidates and to estimate the initial parameters.
ToolHandle< MuonR4::ISegmentSelectionTool >	m_segSelector {this, "SegmentSelectionTool" , "" }
	Selection tool to quantify the segment candidate quality.
ToolHandle< ActsTrk::IExtrapolationTool >	m_extrapolationTool {this, "ExtrapolationTool" ,"" }
	Track extrapolation tool.
TrackKey_t	m_legacyTrackKey {this,"LegacyTrackKey", "MuonSpectrometerTrackParticles"}
	Legacy track reconstruction chain.
TrackKey_t	m_idTrackKey {this, "IdTrackKey", "InDetTrackParticles"}
	The collection of ID tracks associated with the truth particle.
MuonTagKey_t	m_idTagKey {this, "IdTagKey", "MuonInDetCandidates"}
	The collection of all ID tracks that are selected for combined reconstruction.
MuonTagKey_t	m_segTagKey {this, "SegTagKey", "SegmentTags"}
	The collection of the original segment tags.
Gaudi::Property< double >	m_toleranceX0 {this, "toleranceX0", 20.*Gaudi::Units::cm}
	Extra tolerance applied on the non-bending intercept when calculating the matching score.
Gaudi::Property< double >	m_toleranceY0 {this, "toleranceY0", 5.*Gaudi::Units::cm}
	Extra tolerance applied on the bending intercept when calculating the matching score.
Gaudi::Property< double >	m_toleranceTheta {this, "toleranceTheta", 1.*Gaudi::Units::deg}
	Extra tolerance applied on the bending direction when calculating the matching score.
Gaudi::Property< double >	m_tolerancePhi {this, "tolerancePhi", 2.*Gaudi::Units::deg}
	Extra tolerance applied on the bending direction when calculating the matching score.
MuonKey_t	m_legacyMuonKey {this,"LegacyMuonKey", "Muons"}
	Legacy muons.
std::unique_ptr< const Acts::Logger >	m_logger {}
	Instance to the Acts logger.
ParticleBranchPtr_t	m_truthTrks {}
ParticleBranchPtr_t	m_muonTrks {}
	Stored muon information from the Acts muon reco chain.
ParticleBranchPtr_t	m_idTracks {}
	Stored ID track information.
SegmentBranchPtr_t	m_truthSegs {}
SegmentBranchPtr_t	m_recoSegs {}
SegmentBranchPtr_t	m_legacyRecoSegs {}
MuonVal::ScalarBranch< std::uint16_t > &	m_nIdTracks {m_tree.newScalar<std::uint16_t>("nIdTracks", 0)}
	The number of ID tracks in the event.
MuonVal::ScalarBranch< std::uint16_t > &	m_nIdTags {m_tree.newScalar<std::uint16_t>("nIdTags", 0)}
	The number of selected ID tracks in the event.
MuonVal::ThreeVectorBranch	m_seedPos {m_tree, "TrkSeed_position"}
	Simple seed information.
MuonVal::UnitThreeVectorBranch	m_seedDir {m_tree, "TrkSeed_direction"}
	Seed direction vector.
MuonVal::VectorBranch< char > &	m_seedType {m_tree.newVector<char>("TrkSeed_type")}
	Is the seed in the encap or in the barrel chambers.
MuonVal::VectorBranch< int > &	m_seedSector {m_tree.newVector<int>("TrkSeed_sector")}
	Sector of the seed, even center, odd overlap regions, for details see:
MuonVal::VectorBranch< float > &	m_seedLength {m_tree.newVector<float>("TrkSeed_length")}
	Maximum separation between the segments on the reference plane.
MuonVal::VectorBranch< float > &	m_seedThetaCone {m_tree.newVector<float>("TrkSeed_thetaCone")}
	Maximum angular difference between the segments part of the seed.
MuonVal::VectorBranch< char > &	m_seedGood {m_tree.newVector<char>("TrkSeed_goodSeed")}
	Does the seeding tool construct valid parameters from the seed.
MuonVal::VectorBranch< float > &	m_seedQP {m_tree.newVector<float>("TrkSeed_qTimesP")}
	Estimated momentum times charge from the track seed.
MuonVal::VectorBranch< unsigned short > &	m_seedTruthLink {m_tree.newVector<unsigned short>("TrkSeed_truthLink", -1)}
	Link to the truth muon.
std::shared_ptr< TrackSummaryModule >	m_seedSummary {}
	Hit summary on the track seed.
MuonVal::MatrixBranch< unsigned short > &	m_seedRecoSegMatch {m_tree.newMatrix<unsigned short>("TrkSeed_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 38 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

using MuonValR4::MsTrackTester::Location = MuonR4::MsTrackSeed::Location

private

Definition at line 46 of file MsTrackTester.h.

MuonKey_t

using MuonValR4::MsTrackTester::MuonKey_t = SG::ReadHandleKey<xAOD::MuonContainer>

private

Abrivate the key type for the muon container.

Definition at line 138 of file MsTrackTester.h.

MuonTagKey_t

using MuonValR4::MsTrackTester::MuonTagKey_t = SG::ReadHandleKey<MuonR4::MuonTagContainer>

private

Abrivate the muon tag container.

Definition at line 140 of file MsTrackTester.h.

ParticleBranchPtr_t

using MuonValR4::MsTrackTester::ParticleBranchPtr_t = std::shared_ptr<MuonVal::IParticleFourMomBranch>

private

Definition at line 197 of file MsTrackTester.h.

SectorProjector

using MuonValR4::MsTrackTester::SectorProjector = MuonR4::ExpandedSector::SectorProjector

private

Definition at line 47 of file MsTrackTester.h.

SegmentBranchPtr_t

using MuonValR4::MsTrackTester::SegmentBranchPtr_t = std::shared_ptr<MuonPRDTest::SegmentVariables>

private

Abrivation of the branches containing sement information.

Definition at line 205 of file MsTrackTester.h.

SegmentKey_t

using MuonValR4::MsTrackTester::SegmentKey_t = SG::ReadHandleKey<xAOD::MuonSegmentContainer>

private

Abrivate the ReadHandleKey_t for the segment container.

Definition at line 134 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

using MuonValR4::MsTrackTester::TrackKey_t = SG::ReadHandleKey<xAOD::TrackParticleContainer>

private

Abrivate the key type for the track particle container.

Definition at line 136 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 132 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.

                                                                         : 
  ::AthAlgorithm( name, pSvcLocator ),
  ::AthHistogramming( name ),
  m_histSvc ( "THistSvc/THistSvc", name )
{
  //
  // Property declaration
  // 
  // The histogram service
  declareProperty("THistSvc",
                  m_histSvc = ServiceHandle<ITHistSvc>("THistSvc/THistSvc", name),
                  "Handle to a THistSvc instance: it will be used to write "
                  "ROOT objects to ROOT files" );
  //  declareProperty("THistService",      m_histSvc, "The THistSvc" );
  
  declareProperty("RootStreamName",    m_prefix  = "/ANALYSIS", "Name of the output ROOT stream (file) that the THistSvc uses");
  declareProperty("RootDirName",       m_rootDir = "",
                  "Name of the ROOT directory inside the ROOT file where the histograms will go");
 
  declareProperty( "HistNamePrefix",   m_histNamePrefix  = "", "The prefix for the histogram THx name" );
  declareProperty( "HistNamePostfix",  m_histNamePostfix = "",      "The postfix for the histogram THx name" );
 
  declareProperty( "HistTitlePrefix",  m_histTitlePrefix  = "", "The prefix for the histogram THx title" );
  declareProperty( "HistTitlePostfix", m_histTitlePostfix = "", "The postfix for the histogram THx title" );
}

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.

{
  // We need to create a non-const clone
  TEfficiency* effClone = dynamic_cast< TEfficiency* >( eff.Clone() );
  if ( !effClone ) {
    m_msg << MSG::ERROR << "Couldn't create a TEfficiency clone" << endmsg;
    return StatusCode::FAILURE;
  }
  return this->book( *effClone, tDir, 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.

{
  // We need to create a non-const clone
  TH1* histClone = dynamic_cast< TH1* >( hist.Clone() );
  if ( !histClone ) {
    m_msg << MSG::ERROR << "Couldn't create a TH1 clone" << endmsg;
    return StatusCode::FAILURE;
  }
  return this->book( *histClone, tDir, stream );
}

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.

{
  // Call the other Book method and see if it returns a valid pointer
  TTree* treePointer = this->bookGetPointer( treeRef, tDir, stream );
  if ( treePointer )
    {
      return StatusCode::SUCCESS;
    }
  else
    {
      return StatusCode::FAILURE;
    }
}

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.

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

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.

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

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.

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

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.

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

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.

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

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.

{
  // Get a pointer
  const TGraph* graphPointer = &graphRef;
 
  // Check that we got a valid pointer
  if ( !graphPointer )
    {
      m_msg << MSG::WARNING
            << "We got an invalid TGraph pointer in the BookGetPointer(TGraph*) method of the class" << m_name
            << "!" << endmsg;
      return NULL;
    }
 
  // Modify the name and title according to the prefixes of this classes instance
  std::string graphName  = graphPointer->GetName();
  const std::string graphTitle = graphPointer->GetTitle();
 
  // Check if the hash for this graphName already exists, i.e., if we have a hash collision
  const hash_t graphHash = this->hash(graphName);
  GraphMap_t::const_iterator it = m_graphMap.find( graphHash );
  if ( it != m_graphMap.end() ) // It does exist!
    {
      m_msg << MSG::WARNING
            << "Detected a hash collision. The hash for the TGraph with name=" << graphName
            << " already exists and points to a TGraph with name=" << it->second->GetName()
            << " NOT going to book the new histogram and returning a NULL pointer!" << endmsg;
      return NULL;
    }
 
  // Create a clone that has the new name
  TGraph* graphClone = dynamic_cast< TGraph* >( graphPointer->Clone((m_histNamePrefix+graphName+m_histNamePostfix).c_str()) );
  if( !graphClone )
    {
      m_msg << MSG::WARNING
            << "We couldn't clone the TGraph in the BookGetPointer(TGraph&) method of the class" << m_name
            << "!" << endmsg;
      return NULL;
    }
  graphClone->SetTitle ((m_histTitlePrefix+graphTitle+m_histTitlePostfix).c_str());
 
  // Massage the final string to book things
  std::string bookingString("");
  this->buildBookingString( bookingString, graphName, tDir, stream );
 
  // Register the TGraph into the THistSvc
  if ( !((histSvc()->regGraph(bookingString, graphClone)).isSuccess()) )
    {
      m_msg << MSG::WARNING
            << "Problem registering TGraph with name " << graphName
            << ", title " << graphTitle
            << " in " << m_name << "!" << endmsg;
      return NULL;
    }
 
  // Also register it in the local map of string to pointer
  m_graphMap.insert( m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphClone ) );
 
  return graphClone;
}

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.

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

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.

{
  // Get a pointer
  const TTree* treePointer = &treeRef;
 
  // Check that we got a valid pointer
  if ( !treePointer )
    {
      m_msg << MSG::WARNING
            << "We got an invalid TTree pointer in the BookGetPointer(TTree*) method of the class" << m_name
            << "!" << endmsg;
      return NULL;
    }
 
  // Modify the name and title according to the prefixes of this classes instance
  std::string treeName  = treePointer->GetName();
  const std::string treeTitle = treePointer->GetTitle();
 
  // Check if the hash for this treeName already exists, i.e., if we have a hash collision
  const hash_t treeHash = this->hash(treeName);
  TreeMap_t::const_iterator it = m_treeMap.find( treeHash );
  if ( it != m_treeMap.end() ) // It does exist!
    {
      m_msg << MSG::WARNING
            << "Detected a hash collision. The hash for the TTree with name=" << treeName
            << " already exists and points to a TTree with name=" << it->second->GetName()
            << " NOT going to book the new histogram and returning a NULL pointer!" << endmsg;
      return NULL;
    }
 
  // Create a clone that has the new name
  TTree* treeClone = dynamic_cast< TTree* >( treePointer->Clone(treeName.c_str()) );
  if( !treeClone )
    {
      m_msg << MSG::WARNING
            << "We couldn't clone the TTree in the BookGetPointer(TTree&) method of the class" << m_name
            << "!" << endmsg;
      return NULL;
    }
  treeClone->SetTitle (treeTitle.c_str());
 
  // Massage the final string to book things
  std::string bookingString("");
  this->buildBookingString( bookingString, treeName, tDir, stream );
 
  // Register the TTree into the THistSvc
  if ( !((histSvc()->regTree(bookingString, treeClone)).isSuccess()) )
    {
      m_msg << MSG::WARNING
            << "Problem registering TTree with name " << treeName
            << ", title " << treeTitle
            << " in " << m_name << "!" << endmsg;
      return NULL;
    }
 
  // Also register it in the local map of string to pointer
  m_treeMap.insert( m_treeMap.end(), std::pair< const hash_t, TTree* >( treeHash, treeClone ) );
 
  return treeClone;
}

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.

{
  std::string effName(effRef.GetName());
  const std::string effTitle(effRef.GetTitle());
  std::string bookingString;
  this->buildBookingString(bookingString, effName, tDir, stream);
  const std::string finalEffName = m_histNamePrefix + effName + m_histNamePostfix;
  effRef.SetTitle((m_histTitlePrefix + effTitle + m_histTitlePostfix).c_str());
  effRef.SetName(finalEffName.c_str());
 
  const hash_t effHash = this->hash(effName);
  EffMap_t::const_iterator it = m_effMap.find(effHash);
  if (it != m_effMap.end()) {
    m_msg << MSG::WARNING
          << "Detected a hash collision. The hash for the TEfficiency with name=" << effName
          << " already exists and points to a TEfficiency with name=" << it->second->GetName()
          << " NOT going to book the new TEfficiency and returning a NULL pointer!" << endmsg;
    return nullptr;
  }
 
  if (!histSvc()->regEfficiency(bookingString, &effRef).isSuccess()) {
    m_msg << MSG::WARNING
          << "Problem registering TEfficiency with name " << effName
          << ", name prefix " << m_histNamePrefix
          << ", title " << effTitle
          << ", title prefix " << m_histTitlePrefix
          << ", and title postfix " << m_histTitlePostfix
          << " in " << m_name << "!" << endmsg;
    return nullptr;
  }
 
  m_effMap.insert(m_effMap.end(), std::pair<const hash_t, TEfficiency*>(effHash, &effRef));
 
  return &effRef;
}

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.

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

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.

{
  std::string histName(histRef.GetName());
  const std::string histTitle(histRef.GetTitle());
  std::string bookingString;
 
  this->buildBookingString(bookingString, histName, tDir, stream);
 
  const std::string finalHistName = m_histNamePrefix + histName + m_histNamePostfix;
 
  histRef.SetTitle((m_histTitlePrefix + histTitle + m_histTitlePostfix).c_str());
  histRef.SetName(finalHistName.c_str());
 
  const hash_t histHash = this->hash(histName);
  HistMap_t::const_iterator it = m_histMap.find(histHash);
  if (it != m_histMap.end()) {
    m_msg << MSG::WARNING
          << "Detected a hash collision. The hash for the histogram with name=" << histName
          << " already exists and points to a histogram with name=" << it->second->GetName()
          << " NOT going to book the new histogram and returning a NULL pointer!" << endmsg;
    return nullptr;
  }
 
  if (!histSvc()->regHist(bookingString, &histRef).isSuccess()) {
    m_msg << MSG::WARNING
          << "Problem registering histogram with name " << histName
          << ", name prefix " << m_histNamePrefix
          << ", title " << histTitle
          << ", title prefix " << m_histTitlePrefix
          << ", and title postfix " << m_histTitlePostfix
          << " in " << m_name << "!" << endmsg;
    return nullptr;
  }
 
  m_histMap.insert(m_histMap.end(), std::pair<const hash_t, TH1*>(histHash, &histRef));
 
  return &histRef;
}

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.

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

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.

{
  // Massage the final string to book things
  if(tDir.empty()) tDir = m_rootDir;
  size_t pos = histName.rfind('/');
  if(pos != std::string::npos){
    tDir+='/';
    tDir.append(histName, 0,pos);
    histName.erase(0,pos+1);
  };
  if(stream.empty()) stream = m_streamName;
 
  if(usePrefixPostfix){
    bookingString = "/"+stream+"/"+tDir+"/"+m_histNamePrefix+histName+m_histNamePostfix;
  } else {
    bookingString = "/"+stream+"/"+tDir+"/"+histName;
  }
  while(bookingString.find("//") != std::string::npos){
    this->myReplace(bookingString,"//","/");
  }
 
  return;
}

calcSeedLength()

std::pair< double, double > MuonValR4::MsTrackTester::calcSeedLength ( const Acts::GeometryContext & tgContext, 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

gctx	The geometry context to align the segment within ATLAS
seed	The seed with the contributing segments

Definition at line 100 of file MsTrackTester.cxx.

                                                                                                 {
        double maxL{-1.*Gaudi::Units::km}, minL{1.*Gaudi::Units::km},
               maxTheta{-181.}, minTheta{181};
        for (const xAOD::MuonSegment* seg : seed.segments()) {
            const Amg::Vector2D projPos{m_seedingTool->expressOnCylinder(tgContext, *seg, seed.location(), seed.sector())};
            if (!m_seedingTool->withinBounds(projPos, seed.location())) {
                continue;
            }
            const double projected = projPos[seed.location()==Location::Barrel];
            const double theta = seg->direction().theta();
            minL = std::min(minL, projected);
            maxL = std::max(maxL, projected);
            minTheta = std::min(minTheta, toDeg(theta));
            maxTheta = std::max(maxTheta, toDeg(theta));
        }
        return std::make_pair(maxL - minL, maxTheta - minTheta);
    }

calcSegTagVariables()

MsTrackTester::SegmentTagVariables MuonValR4::MsTrackTester::calcSegTagVariables ( const xAOD::TrackParticle & idTrack, const EventContext & ctx )

private

No associated segments found

Calo extension does not work

In case that there is already a

Store in the tree that the segment cannot be associated due to extrapolation failure

If there is a best segment dump the best segment to the tree

The segment does not measure phi. Reset anything in loc0 and non-precision direction

Definition at line 538 of file MsTrackTester.cxx.

                                                                     {
 
        SegmentTagVariables result{};
        const std::vector<const xAOD::MuonSegment*> compatibleSegs = getAssociatedSegments(idTrack, ctx);
        if (compatibleSegs.empty()) {
            return result;
        } 
        const ActsTrk::CaloExtension* caloExt = ActsTrk::getCaloExtension(idTrack);
        if (!caloExt) {
            return result; 
        }
        std::optional<Acts::BoundTrackParameters> startPars = caloExt->lastParameters();
        if (!startPars) {
            return result;
        }
 
        const MuonR4::MuonTag* muTag = findMuTagIMO(idTrack, ctx);
    
        const Acts::GeometryContext tgContext = m_ctxProvider.getGeometryContext(ctx);
 
        auto nextParameters = [&](const xAOD::MuonSegment& segment) {
            const MuonGMR4::SpectrometerSector* msSector = getEnvelope(segment);
            const Acts::Surface& target{msSector->surface()};
            if (target.geometryId() == startPars->referenceSurface().geometryId()) {
                return true;
            }
            if (muTag) {
                auto surfPars = muTag->extrapolatedParsID(Acts::toUnderlying(segment.chamberIndex()));
                if (surfPars) {
                    startPars = surfPars;
                    ATH_MSG_VERBOSE(__func__<<"() "<<__LINE__<<" - Recycle cached parameters from muTag\n"
                                    <<(*startPars));
                    return true;
                }
                ATH_MSG_DEBUG(__func__<<"() "<<__LINE__<<" - Extrapolation to "<<printID(segment)<<" failed somehow");
            }
 
            auto extpPars = m_extrapolationTool->propagate(ctx, *startPars, target);
            if (!extpPars.ok()) {
                ATH_MSG_DEBUG(__func__<<"() "<<__LINE__<<" - Extrapolation keeps failing");
                return false;
            }
            ATH_MSG_DEBUG(__func__<<"() "<<__LINE__<<" - Updated parameters to \n"<<(*extpPars));
            startPars = *extpPars;
            return true;
        };
 
        const xAOD::MuonSegment* bestSeg{nullptr};
        float bestMatchScore{std::numeric_limits<float>::max()};
        Acts::BoundVector bestPars{Acts::BoundVector::Zero()};
 
        auto dumpBestSeg = [&]() {
            if (!bestSeg) {
                return;
            }
            result.recoSegs.push_back(m_recoSegs->push_back(*bestSeg));
            result.matchScores.push_back(bestMatchScore);
            result.deltaPhi.push_back(toDeg(bestPars[Acts::eBoundPhi]));
            result.deltaTheta.push_back(toDeg(bestPars[Acts::eBoundTheta]));
            result.deltaY0.push_back(bestPars[Acts::eBoundLoc1]);
            result.deltaX0.push_back(bestPars[Acts::eBoundLoc0]);
            result.goodExtp.push_back(1);
            result.taggedSeg.push_back(muTag && Acts::rangeContainsValue(muTag->segments(), bestSeg));
            bestSeg = nullptr;
            bestMatchScore = std::numeric_limits<float>::max();
            bestPars = Acts::BoundVector::Zero();
        };
 
        for (auto segItr = compatibleSegs.begin(); segItr != compatibleSegs.end(); ) {
            const xAOD::MuonSegment* testMe{*segItr};
            if (!nextParameters(*testMe)) {
                result.recoSegs.push_back(m_recoSegs->push_back(*testMe));
                result.matchScores.push_back(-1.f);
                result.deltaX0.push_back(1._m);
                result.deltaY0.push_back(1._m);
                result.deltaTheta.push_back(180.);
                result.deltaPhi.push_back(180.);
                result.goodExtp.push_back(0);
                result.taggedSeg.push_back(0);
                segItr = std::find_if(std::next(segItr),compatibleSegs.end(), [&](const xAOD::MuonSegment* assocSeg){
                    return getEnvelope(*testMe) != getEnvelope(*assocSeg);
                });
            } else {
                ++segItr;
            } 
            if (bestSeg && getEnvelope(*testMe) != getEnvelope(*bestSeg)) {
                dumpBestSeg();
            }
            const Acts::BoundTrackParameters segmentPars{MuonR4::SegmentFit::boundSegmentPars(tgContext, *m_detMgr, *testMe)};
            Acts::BoundVector dPars = segmentPars.parameters() - startPars->parameters();
            if (!testMe->nPhiLayers()) {
                dPars[Acts::eBoundPhi] = dPars[Acts::eBoundLoc0] = 0.;
            } else {
                dPars[Acts::eBoundPhi] = P4Helpers::deltaPhi(dPars[Acts::eBoundPhi], 0.);
            }
            
            dPars[Acts::eBoundQOverP] = dPars[Acts::eBoundTime] = 0.;
            Acts::BoundMatrix covariance{Acts::BoundMatrix::Identity()};
            covariance(Acts::eBoundLoc0, Acts::eBoundLoc0) = Acts::square(m_toleranceX0.value());
            covariance(Acts::eBoundLoc1, Acts::eBoundLoc1) = Acts::square(m_toleranceY0.value());
            covariance(Acts::eBoundTheta, Acts::eBoundTheta) = Acts::square(m_toleranceTheta.value());
            covariance(Acts::eBoundPhi, Acts::eBoundPhi) = Acts::square(m_tolerancePhi.value());
            if (startPars->covariance()) {
                covariance += (*startPars->covariance());
            }
            if (segmentPars.covariance()) {
                covariance += (*segmentPars.covariance());
            }
            const float chi2 = dPars.dot(covariance.inverse()*dPars);
            ATH_MSG_VERBOSE(__func__<<"() "<<__LINE__<<" - Difference: "<<Acts::toString(dPars)
                        <<", covariance: \n"<<Acts::toString(covariance)<<",\nchi2: "<<chi2);
            if (chi2 < bestMatchScore) {
                ATH_MSG_DEBUG(__func__<<"() "<<__LINE__<<" - Best parameters on surface thus far");
                bestMatchScore = chi2;
                bestSeg = testMe;
                bestPars = std::move(dPars);
            }
        }
        dumpBestSeg();
        return result;
    }

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.

{
  m_histSvc          = histSvc;
  m_streamName       = prefix;
  m_rootDir          = rootDir;
  m_histNamePrefix   = histNamePrefix;
  m_histNamePostfix  = histNamePostfix;
  m_histTitlePrefix  = histTitlePrefix;
  m_histTitlePostfix = histTitlePostfix;
 
  return StatusCode::SUCCESS;
}

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.

  {
    return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                       hndl.value(), 
                                                       hndl.documentation());
 
  }

declareProperty()

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

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

                                                                     {
    typedef typename SG::HandleClassifier<T>::type htype;
    return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
  }

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.

{ return m_detStore; }

dumpLegacyTracks()

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

private

Dumps the legacy containers to the TTree.

Definition at line 668 of file MsTrackTester.cxx.

                                                                      {
 
        //This for now is to be able to retrieve the matching between the legacy segments and tracks ...
        const xAOD::MuonContainer* legacyMuons{nullptr};
        ATH_CHECK(SG::get(legacyMuons, m_legacyMuonKey, ctx));
 
        if (!legacyMuons) {
            return StatusCode::SUCCESS;
        }
 
        for (const xAOD::Muon* muon : *legacyMuons) {
            ATH_MSG_VERBOSE("iMuon " << muon->index() << " pT: "<<(muon->pt() *MeVtoGeV)<<" [GeV], eta: "
                    <<muon->eta() <<", phi: "<<toDeg(muon->phi())<<", q: "<<muon->charge()  
                    <<", nSegments: "<<muon->nMuonSegments());
            const xAOD::TrackParticle* track = muon->trackParticle(xAOD::Muon::TrackParticleType::MuonSpectrometerTrackParticle);
            if (!track) {
                continue;
            }
            m_summaryTool->copySummary(m_summaryTool->makeSummary(ctx, *track->track()), *track);
            m_legacyTrks->push_back(track);
            auto trkIdx = m_legacyTrks->find(track);
            for (size_t s = 0; s < muon->nMuonSegments(); ++s) {
                const xAOD::MuonSegment* segment = muon->muonSegment(s);
                auto segIdx = m_legacyRecoSegs->push_back(*segment);
                ATH_MSG_VERBOSE(std::format( "Legacy muon-segment link: segment {:}  @{:}, eta: {:.2f}, phi {:.2f}", 
                                            printID(*segment), Amg::toString(segment->position()),
                                    segment->direction().eta(), toDeg(segment->direction().phi())));
                m_legacySegToTrkLinks[segIdx] = trkIdx;
            }
        }
 
        //Dump also legacy segments
        const xAOD::MuonSegmentContainer* legacyRecoSegs{nullptr};
        ATH_CHECK(SG::get(legacyRecoSegs, m_legacySegmentKey, ctx));
 
        if (legacyRecoSegs->size()) {
            m_legacySegToTrkLinks[legacyRecoSegs->size()-1];
            for (const xAOD::MuonSegment* seg : *legacyRecoSegs) {
                //Store all segments
                m_legacyRecoSegs->push_back(*seg);
            }
        }
 
        const xAOD::TrackParticleContainer* legacyTrks{nullptr};
        ATH_CHECK(SG::get(legacyTrks, m_legacyTrackKey, ctx));
 
        for (const xAOD::TrackParticle* track : *legacyTrks) {
            m_summaryTool->copySummary(m_summaryTool->makeSummary(ctx, *track->track()), *track);
            m_legacyTrks->push_back(track);
        }
 
        return StatusCode::SUCCESS;
    }

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

Collect only the ID tracks that belong to a truth muon of interest.

Select only the id tags that are extrapolated to the calo exit

Definition at line 758 of file MsTrackTester.cxx.

                                                                     {
        const xAOD::MuonContainer* muons{nullptr};
        ATH_CHECK(SG::get(muons, m_muonKey, ctx));
        const xAOD::MuonSegmentContainer* recoSegments{nullptr};
        ATH_CHECK(SG::get(recoSegments, m_recoSegmentKey, ctx));
 
        const MuonR4::MsTrackSeedContainer* trkSeeds{nullptr};
        ATH_CHECK(SG::get(trkSeeds, m_msTrkSeedKey, ctx));
 
        const xAOD::TrackParticleContainer* idTracks{nullptr};
        ATH_CHECK(SG::get(idTracks, m_idTrackKey, ctx));
 
        const MuonR4::MuonTagContainer* idTags{nullptr};
        ATH_CHECK(SG::get(idTags, m_idTagKey, ctx));
 
        const Acts::GeometryContext tgContext = m_ctxProvider.getGeometryContext(ctx);
        const Acts::MagneticFieldContext mfContext = m_ctxProvider.getMagneticFieldContext(ctx);
 
        MagField::AtlasFieldCache magField{};
        mfContext.get<const AtlasFieldCacheCondObj*>()->getInitializedCache(magField);
 
        std::unordered_map<const xAOD::TruthParticle*, 
                           std::vector<unsigned>> truthToSeedMatchCounter{};
 
        if (!trkSeeds->empty()) {
            m_seedTruthLink[trkSeeds->size() -1];
        }
 
        for (const MuonR4::MsTrackSeed& seed : *trkSeeds) {
            unsigned int seedIdx = m_seedPos.size();
            m_seedPos += seed.position();
            m_seedType+= Acts::toUnderlying(seed.location());
            m_seedSector += seed.sector().sector();
            m_seedSummary->push_back(ctx, seed);
 
            auto startPars = m_seedingTool->estimateStartParameters(ctx, seed);
            if (startPars.ok()) {
                m_seedDir += (*startPars).direction();
            } else {
                m_seedDir += Amg::Vector3D::UnitZ();
            }
            // m_seedDir
            ATH_MSG_VERBOSE(" Dump new seed: "<<seed);
            for (const xAOD::MuonSegment* seg : seed.segments()){
                m_seedRecoSegMatch[seedIdx].push_back(m_recoSegs->push_back(*seg));
                if (const xAOD::MuonSegment* truthSeg =  MuonR4::getMatchedTruthSegment(*seg);
                    truthSeg != nullptr) {
                    std::vector<unsigned>& matchCounter = truthToSeedMatchCounter[MuonR4::getTruthMatchedParticle(*truthSeg)];
                    if (seedIdx >= matchCounter.size()) {
                        matchCounter.resize(seedIdx +1);
                    }
                    ++matchCounter[seedIdx];
                }
            }
            const auto[seedLength, theta] = calcSeedLength(tgContext, seed);
            m_seedLength+= seedLength;
            m_seedThetaCone+=theta;
            m_seedQP += m_seedingTool->estimateQtimesP(tgContext, seed, magField) / Gaudi::Units::GeV; 
            m_seedGood += startPars.ok();
        }
        for (auto& [truthMuon, matches] : truthToSeedMatchCounter) {
            const unsigned tIndex = m_truthTrks->find(truthMuon);
            if (tIndex >= m_truthTrks->size()) {
                continue;
            }
            while (std::count_if(matches.begin(), matches.end(), 
                   [](const unsigned nMatched){
                       return nMatched > 0;
                   })) {
                auto bestMatch = std::ranges::max_element(matches);
                const std::size_t seedIdx = std::distance(matches.begin(), bestMatch);
                m_truthMuToSeedIdx[tIndex].push_back(seedIdx);
 
                m_seedTruthLink[seedIdx] = tIndex;
                m_truthMuToSeedCounter[tIndex].push_back(*bestMatch);
                (*bestMatch) = 0;
            }
        }
 
        for (const xAOD::Muon* muon : *muons) {
            m_muonTrks->push_back(muon);
        }
 
        for (const xAOD::MuonSegment* seg : *recoSegments) {
            m_recoSegs->push_back(*seg);
        }

        if (idTracks) {
            m_nIdTracks = idTracks->size();
            for (const xAOD::TrackParticle* idTrk : *idTracks) {
                if (m_truthTrks->find([idTrk](const xAOD::IParticle* p){
                    return xAOD::TruthHelpers::getTruthParticle(*idTrk) == 
                          xAOD::TruthHelpers::getTruthParticle(*p);
                }) < m_truthTrks->size()) {
                    m_idTracks->push_back(idTrk);
                }
            }
        }
        if (idTags) {
            m_nIdTags = idTags->size();
            for (const MuonR4::MuonTag* idTag : *idTags) {
                if (idTag->extrapolatedParsID(Acts::hashString("@CaloExit"))){
                    m_idTracks->push_back(idTag->idTrack());
                }
            }
        }
        return StatusCode::SUCCESS;
    }

dumpTruthContent()

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

private

Dump truth information.

Allocate the memory

Definition at line 721 of file MsTrackTester.cxx.

                                                                      {
        if (!m_isMC) {
            return StatusCode::SUCCESS;
        }
        const xAOD::MuonSegmentContainer* truthSegs{nullptr};
        ATH_CHECK(SG::get(truthSegs, m_truthSegmentKey, ctx));
 
 
        for (const xAOD::MuonSegment* seg : *truthSegs) {
            ATH_MSG_VERBOSE("Dump truth segment "<<printID(*seg)<<" @"<<
            Amg::toString(seg->position())<<", eta: "<<seg->direction().eta()
            <<", phi: "<<toDeg(seg->direction().phi()));
            m_truthSegs->push_back(*seg);
        }
        if (truthSegs->size()) {
            m_truthSegToRecoLink[truthSegs->size()-1];
        }
 
 
        const xAOD::TruthParticleContainer* truthMuons{nullptr};
        ATH_CHECK(SG::get(truthMuons, m_truthKey, ctx));
        const std::size_t nT = truthMuons->size() - 1;
        if (!truthMuons->empty()) {
            m_truthMuToSeedIdx[nT];
            m_truthMuToSeedCounter[nT];
            m_truthMuRecoSegLinks[nT];
        }
 
        for (const xAOD::TruthParticle* truth : *truthMuons) {
            ATH_MSG_DEBUG("Truth muon: pT: "<<(truth->pt() *MeVtoGeV)
                <<", eta: "<<truth->eta()<<", phi: "<<toDeg(truth->phi())<<", q: "<<truth->charge());
            m_truthTrks->push_back(*truth);
        }
        
        return StatusCode::SUCCESS;
    }

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.

{
  // Build a 32 bit hash out of the name
  const hash_t effHash = this->hash(effName);
 
  // See if this entry exists in the map
  EffMap_t::const_iterator it = m_effMap.find( effHash );
  if ( it == m_effMap.end() ) // It doesn't exist!
    { // Let's see into the THistSvc if somebody else has registered the TEfficiency...
 
      // Need to copy the strings as we will massage them from here on
      std::string effNameCopy = effName;
      std::string tDirCopy     = tDir;
      std::string streamCopy   = stream;
 
      // Massage the final string to book things
      std::string bookingString("");
      this->buildBookingString( bookingString, effNameCopy, tDirCopy, streamCopy ,false);
 
      TEfficiency* effPointer(NULL);
      if ( !((histSvc()->getEfficiency(bookingString, effPointer)).isSuccess()) )
        {
          // Massage the final string to book things
          std::string bookingString("");
          this->buildBookingString( bookingString, effNameCopy, tDirCopy, streamCopy, true );
 
          if ( !((histSvc()->getEfficiency(bookingString, effPointer)).isSuccess()) )
            {
              m_msg << MSG::WARNING
                    << "Problem retrieving the TEfficiency with name (including pre- and post-fixes) "
                    << m_histNamePrefix + effNameCopy + m_histNamePostfix
                    << " or with name " << effNameCopy
                    << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
                    << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
              return NULL;
            }
          // If we get to here, we actually found the TEfficiency in the THistSvc.
          // So let's add it to the local cache map and return its pointer
          m_effMap.insert( m_effMap.end(), std::pair< const hash_t, TEfficiency* >( effHash, effPointer ) );
          return effPointer;
        }
      // If we get to here, we actually found the TEfficiency in the THistSvc.
      // So let's add it to the local cache map and return its pointer
      m_effMap.insert( m_effMap.end(), std::pair< const hash_t, TEfficiency* >( effHash, effPointer ) );
      return effPointer;
    }
 
  // Return the pointer to the TEfficiency that we got from the local cache map
  return it->second;
}

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.

{ return m_evtStore; }

execute()

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

finaloverridevirtual

Execute method.

Implements AthAlgorithm.

Definition at line 871 of file MsTrackTester.cxx.

                                                             {
        ATH_CHECK(dumpTruthContent(ctx));
        ATH_CHECK(dumpLegacyTracks(ctx));
        ATH_CHECK(dumpRecoContent(ctx));
 
        ATH_CHECK(m_tree.fill(ctx));
        return StatusCode::SUCCESS;
    }

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.

{
  // If we didn't find any symlinks to add, just return the collection
  // from the base class.  Otherwise, return the extended collection.
  if (!m_extendedExtraObjects.empty()) {
    return m_extendedExtraObjects;
  }
  return BaseAlg::extraOutputDeps();
}

filterPassed()

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

inlinevirtualinherited

Get filter decision:

Definition at line 93 of file AthCommonAlgorithm.h.

                                                           {
    return execState( ctx ).filterPassed();
  }

finalize()

StatusCode MuonValR4::MsTrackTester::finalize ( )

finaloverride

Definition at line 879 of file MsTrackTester.cxx.

                                       {
        ATH_CHECK(m_tree.write());
        return StatusCode::SUCCESS;
    }

findBaseIdTag()

const MuonR4::MuonTag * MuonValR4::MsTrackTester::findBaseIdTag ( const xAOD::TrackParticle & idTrack, const EventContext & ctx ) const

private

Searches the Id track in the basline collection of all ID tracks selected for the combined muon reconstructiion (STACO/MuidCo/MuTagIMO/Calo).

Parameters

idTrack	The inner detector track of interest
ctx	The current event context to access the segment tag collection from store gate

Definition at line 494 of file MsTrackTester.cxx.

                                                                                       {
        return findTagMatchingToId(idTrack, m_idTagKey, ctx);
    }

findMuTagIMO()

const MuonR4::MuonTag * MuonValR4::MsTrackTester::findMuTagIMO ( const xAOD::TrackParticle & idTrack, const EventContext & ctx ) const

private

Definition at line 498 of file MsTrackTester.cxx.

                                                                                       {
        return findTagMatchingToId(idTrack, m_segTagKey, ctx);
    }

getAssociatedSegments()

std::vector< const xAOD::MuonSegment * > MuonValR4::MsTrackTester::getAssociatedSegments ( const xAOD::TrackParticle & idTrack, const EventContext & ctx ) const

private

Returns the collection of reconstructed segments that can be associated to the id track If the track is connected with a truth muon, the list of reconstructed segments matched to the truth muon is returned.

Otherwise, the list of segments associated with the MuTagIMO tag is returned.

Parameters

idTrack	The inner detector track of interest
ctx	The current event context to access the segment tag collection from store gate

Definition at line 504 of file MsTrackTester.cxx.

                                                                            {
        std::vector<const xAOD::MuonSegment*> compatibleSegs{};
        if (const xAOD::TruthParticle* truthPart = truthTreeParticle(idTrack); truthPart != nullptr) {
            std::vector<const xAOD::MuonSegment*> truthSegs = getTruthSegments(*truthPart);
            const xAOD::MuonSegmentContainer* recoSegs{nullptr};
            SG::get(recoSegs, m_recoSegmentKey, ctx).ignore();
            std::copy_if(recoSegs->begin(), recoSegs->end(), std::back_inserter(compatibleSegs), 
                        [&truthSegs](const xAOD::MuonSegment* recoSeg){
                             return Acts::rangeContainsValue(truthSegs, getMatchedTruthSegment(*recoSeg));
                        });
            if (!compatibleSegs.empty()) {
                return compatibleSegs;
            }
        }
        
        const xAOD::MuonContainer* muons{nullptr};
        SG::get(muons,m_muonKey, ctx).ignore();
        for (const xAOD::Muon* muon : *muons) {
            if (muon->trackParticle(xAOD::Muon::TrackParticleType::InnerDetectorTrackParticle) == &idTrack) {
                for (std::size_t s = 0; s < muon->nMuonSegments() ; ++s) {
                    compatibleSegs.push_back(muon->muonSegment(s));
                }
                break;
            }
        }
        return compatibleSegs;
    }

getEnvelope()

const MuonGMR4::SpectrometerSector * MuonValR4::MsTrackTester::getEnvelope ( const xAOD::MuonSegment & segment ) const

private

Definition at line 533 of file MsTrackTester.cxx.

                                                                                                     {
        return m_detMgr->getSectorEnvelope(segment.chamberIndex(), segment.sector(), segment.etaIndex()); 
    }

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.

{
  // Build a 32 bit hash out of the name
  const hash_t graphHash = this->hash(graphName);
 
  // See if this entry exists in the map
  GraphMap_t::const_iterator it = m_graphMap.find( graphHash );
  if ( it == m_graphMap.end() ) // It doesn't exist!
    { // Let's see into the THistSvc if somebody else has registered the TGraph...
 
      // Need to copy the strings as we will massage them from here on
      std::string graphNameCopy = graphName;
      std::string tDirCopy      = tDir;
      std::string streamCopy    = stream;
 
      // Massage the final string to book things
      std::string bookingString("");
      this->buildBookingString( bookingString, graphNameCopy, tDirCopy, streamCopy, false);
 
      TGraph* graphPointer(nullptr);
      if ( !((histSvc()->getGraph(bookingString, graphPointer)).isSuccess()) )
        {
          // Massage the final string to book things
          std::string bookingString("");
          this->buildBookingString( bookingString, graphNameCopy, tDirCopy, streamCopy, true );
 
          if ( !((histSvc()->getGraph(bookingString, graphPointer)).isSuccess()) )
            {
              m_msg << MSG::WARNING
                    << "Problem retrieving the TGraph with name (including pre- and post-fixes) "
                    << m_histNamePrefix + graphNameCopy + m_histNamePostfix
                    << " or with name " << graphNameCopy
                    << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
                    << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
              return nullptr;
            }
          // If we get to here, we actually found the TGraph in the THistSvc.
          // So let's add it to the local cache map and return its pointer
          m_graphMap.insert( m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphPointer ) );
          return graphPointer;
        }
      // If we get to here, we actually found the TGraph in the THistSvc.
      // So let's add it to the local cache map and return its pointer
      m_graphMap.insert( m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphPointer ) );
      return graphPointer;
    }
 
 
  // Return the pointer to the TGraph that we got from the local cache map
  return it->second;
}

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.

{
  const uint64_t hash64 = CxxUtils::crc64( histName );
  return (hash_t)(hash64 & 0xFFFFFFFF);
}

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.

{
  std::string histNameCopy = histName;
  std::string tDirCopy     = tDir;
  std::string streamCopy   = stream;
 
  std::string bookingString;
  this->buildBookingString(bookingString, histNameCopy, tDirCopy, streamCopy, false);
 
  const hash_t histHash = this->hash(histNameCopy);
 
  HistMap_t::const_iterator it = m_histMap.find(histHash);
  if (it == m_histMap.end()) {
    TH1* histPointer(nullptr);
 
    if (!histSvc()->getHist(bookingString, histPointer).isSuccess()) {
      std::string prefixedHistNameCopy = histName;
      std::string prefixedTDirCopy     = tDir;
      std::string prefixedStreamCopy   = stream;
 
      std::string prefixedBookingString;
      this->buildBookingString(prefixedBookingString,
                               prefixedHistNameCopy,
                               prefixedTDirCopy,
                               prefixedStreamCopy,
                               true);
 
      if (!histSvc()->getHist(prefixedBookingString, histPointer).isSuccess()) {
        m_msg << MSG::WARNING
              << "Problem retrieving the histogram with name (including pre- and post-fixes) "
              << m_histNamePrefix + histNameCopy + m_histNamePostfix
              << " or with name " << histNameCopy
              << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
              << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
        return nullptr;
      }
    }
 
    m_histMap.insert(m_histMap.end(), std::pair<const hash_t, TH1*>(histHash, histPointer));
    return histPointer;
  }
 
  return it->second;
}

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.

{
  // Get the TH1 pointer
  TH1* th1Pointer = this->hist(histName, tDir, stream);
  if ( !th1Pointer )
    {
      m_msg << MSG::ERROR
            << "Cannot get a 2-d histogram with name " << histName
            << "... will probably seg-fault!" << endmsg;
      return NULL;
    }
  // If the TH1 pointer is valid, simply return the dynamic_cast
  return dynamic_cast<TH2*>( th1Pointer );
}

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.

{
  // Get the TH1 pointer
  TH1* th1Pointer = this->hist(histName, tDir, stream);
  if ( !th1Pointer )
    {
      m_msg << MSG::ERROR
            << "Cannot get a 3-d histogram with name " << histName
            << "... will probably seg-fault!" << endmsg;
      return NULL;
    }
  // If the TH1 pointer is valid, simply return the dynamic_cast
  return dynamic_cast<TH3*>( th1Pointer );
}

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.

{
  return m_histSvc;
}

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

Dump whether the id track is selected and dump the parameters at calo exit

Link the reconstructed segments to the muon

Link the associated seed

Definition at line 118 of file MsTrackTester.cxx.

                                         {
        ATH_CHECK(m_truthSegmentKey.initialize(m_isMC));
        ATH_CHECK(m_truthKey.initialize(m_isMC));
 
        ATH_CHECK(m_muonKey.initialize());
        ATH_CHECK(m_msTrkSeedKey.initialize());
        ATH_CHECK(m_recoSegmentKey.initialize());
 
        ATH_CHECK(m_summaryTool.retrieve());
        ATH_CHECK(m_seedingTool.retrieve());
 
        ATH_CHECK(m_legacyMuonKey.initialize(!m_legacyMuonKey.empty()));
        ATH_CHECK(m_legacyTrackKey.initialize(!m_legacyMuonKey.empty()));
        ATH_CHECK(m_legacySegmentKey.initialize(!m_legacyMuonKey.empty()));
        ATH_CHECK(m_idTrackKey.initialize(m_isMC && m_storeID));
        ATH_CHECK(m_idTagKey.initialize(!m_idTrackKey.empty()));
        ATH_CHECK(m_segTagKey.initialize(!m_idTrackKey.empty()));
        
        ATH_CHECK(m_ctxProvider.initialize());
        ATH_CHECK(m_extrapolationTool.retrieve(EnableTool{m_storeID}));
        ATH_CHECK(detStore()->retrieve(m_detMgr));
        m_logger = makeActsAthenaLogger(this, name());
 
        int evOpts{0};
 
        m_recoSegs = std::make_unique<SegmentVariables>(m_tree, m_recoSegmentKey.key(), "Segments", msgLevel());
        m_recoSegs->addVariable(std::make_unique<MuonVal::GenericAuxDecorationBranch<unsigned short>>(m_tree, 
            "Segments_passSeedQual",[this](const SG::AuxElement* aux){
            const auto* seg = static_cast<const xAOD::MuonSegment*>(aux);
            return m_segSelector->passSeedingQuality(Gaudi::Hive::currentContext(),
                                                     *seg); }));
        m_recoSegs->addVariable(std::make_unique<MuonVal::GenericAuxDecorationBranch<unsigned short>>(m_tree, 
                                "Segments_passTrackQual",[this](const SG::AuxElement* aux){
            const auto* seg = static_cast<const xAOD::MuonSegment*>(aux);
            return m_segSelector->passTrackQuality(Gaudi::Hive::currentContext(),
                                                   *seg); 
        }));
 
 
        if (m_isMC) {
            evOpts |= EventInfoBranch::isMC;
            m_recoSegs->addVariable(std::make_unique<MuonVal::GenericAuxDecorationBranch<unsigned short>>(m_tree, 
                "Segments_truthSegLink",[this](const SG::AuxElement* aux){
                const auto* seg = static_cast<const xAOD::MuonSegment*>(aux);
                const xAOD::MuonSegment* truthS = MuonR4::getMatchedTruthSegment(*seg);
                const unsigned linkIdx = truthS ? m_truthSegs->push_back(*truthS) : -1;
                if (truthS) {
                    m_truthSegToRecoLink.push_back(linkIdx, m_recoSegs->push_back(*seg));
                    const xAOD::TruthParticle* truthMuon = MuonR4::getTruthMatchedParticle(*truthS);
                    if (truthMuon){
                        m_truthTrks->push_back(truthMuon);
                        m_truthMuRecoSegLinks[m_truthTrks->find(truthMuon)].push_back(m_recoSegs->push_back(*seg));
                    }
                }
                return linkIdx;
            }));
 
            m_truthSegs = std::make_unique<SegmentVariables>(m_tree, m_truthSegmentKey.key(), "TruthSegments", msgLevel());
            
            m_truthSegs->addVariable(std::make_unique<MuonVal::GenericAuxDecorationBranch<unsigned short>>(m_tree, 
                    "TruthSegments_truthLink",[this](const SG::AuxElement* aux){
                    const auto* seg = static_cast<const xAOD::MuonSegment*>(aux);
                    const xAOD::TruthParticle* truthP = MuonR4::getTruthMatchedParticle(*seg);
                    m_truthTrks->push_back(truthP);
                    unsigned short linkIdx = m_truthTrks->find(truthP);                    
                    return linkIdx;
                }));
            for (auto loc : {Location::Barrel, Location::Endcap}) {
                m_truthSegs->addVariable(std::make_unique<MuonVal::GenericAuxDecorationBranch<unsigned short>>(m_tree,
                    std::format("TruthSegments_has{}Proj", loc), [loc, this](const SG::AuxElement* aux) -> unsigned short {
                    
                    const auto* seg = static_cast<const xAOD::MuonSegment*>(aux);
                    const Acts::GeometryContext tgContext = m_ctxProvider.getGeometryContext(Gaudi::Hive::currentContext());
                    ExpandedSector sector{seg->position().phi()};
                    const Amg::Vector2D projPos{m_seedingTool->expressOnCylinder(tgContext, *seg, loc, sector)};
                    return m_seedingTool->withinBounds(projPos, loc);
                }));
            }
 
            m_tree.addBranch(m_truthSegs);
 
            m_truthTrks = std::make_unique<IParticleFourMomBranch>(m_tree, "TruthMuons");
            m_truthTrks->addVariable<int>(-1, "truthOrigin");
            m_truthTrks->addVariable<int>(-1, "truthType");
            m_truthTrks->addVariable(
                std::make_unique<GenericPartDecorBranch<xAOD::TruthParticle,
                                                        std::vector<unsigned short>>>(m_tree, 
                std::format("{:}_truthSegLinks", m_truthTrks->name()), [&] (const xAOD::TruthParticle& p){
                    std::vector<unsigned short> idx{};
                    for (const xAOD::MuonSegment* truthSeg: getTruthSegments(p)){
                        idx.push_back(m_truthSegs->push_back(*truthSeg));
                    }
                    return idx;
                }));
            m_truthTrks->addVariable(
                std::make_unique<GenericPartDecorBranch<xAOD::TruthParticle, unsigned short>>(m_tree, 
                std::format("{:}_nTruthSegments", m_truthTrks->name()), [&] (const xAOD::TruthParticle& p) -> unsigned short {
                    return getTruthSegments(p).size();
                }));
 
            m_truthTrks->addVariable(
                std::make_unique<GenericPartDecorBranch<xAOD::TruthParticle, float>>(m_tree,
                std::format("{:}_eLossInMS", m_truthTrks->name()), [] (const xAOD::TruthParticle& p) -> float {
                    std::vector<const xAOD::MuonSimHit*> allHits{};
                    for (const xAOD::MuonSegment* seg : getTruthSegments(p)){
                        auto hits = getMatchingSimHits(*seg);
                        allHits.insert(allHits.end(), hits.begin(), hits.end());
                    }
                    if (allHits.empty()) {
                        return 0.;
                    }
                    auto [min, max] = std::ranges::minmax(allHits, [](const xAOD::MuonSimHit* a, const xAOD::MuonSimHit* b){
                        return a->kineticEnergy() < b->kineticEnergy();
                    });
                    return (min->kineticEnergy() - max->kineticEnergy()) / Gaudi::Units::GeV;
                }));
            auto cone = std::make_shared<VectorBranch<float>>(m_tree,
                                        std::format("{}_seedThetaCone", m_truthTrks->name()));
            auto qTimesP = std::make_shared<VectorBranch<float>>(m_tree,
                                        std::format("{}_qTimesP", m_truthTrks->name()));
 
            m_tree.addBranch(cone);
            m_tree.addBranch(qTimesP);
            m_truthTrks->addVariable(
                std::make_unique<GenericPartDecorBranch<xAOD::TruthParticle, float>>(m_tree, 
                std::format("{:}_seedLength", m_truthTrks->name()), [cone, qTimesP, this] (const xAOD::TruthParticle& p) -> float {
                    
                    const EventContext& ctx{Gaudi::Hive::currentContext()};
                    const Acts::GeometryContext tgContext = m_ctxProvider.getGeometryContext(ctx);
                    const Acts::MagneticFieldContext mfContext = m_ctxProvider.getMagneticFieldContext(ctx);
                    MagField::AtlasFieldCache magField{};
                    mfContext.get<const AtlasFieldCacheCondObj*>()->getInitializedCache(magField);
                    auto truthSeed = makeSeedFromTruth(tgContext, p);
                    if (!truthSeed) {
                        cone->push_back(-1);
                        qTimesP->push_back(0);
                        return -1.;
                    }
                    auto [length, theta] = calcSeedLength(tgContext, *truthSeed);
                    cone->push_back(theta);
                    qTimesP->push_back(m_seedingTool->estimateQtimesP(tgContext, *truthSeed, magField) / Gaudi::Units::GeV);
                    return length;
                }));
 
            m_tree.addBranch(m_truthTrks);
            m_trkTruthLinks.emplace_back(m_truthSegmentKey, "truthParticleLink");
            m_trkTruthLinks.emplace_back(m_truthKey, "truthSegmentLinks");
            m_trkTruthLinks.emplace_back(m_recoSegmentKey, "truthSegmentLink");
            m_trkTruthLinks.emplace_back(m_muonKey, "truthParticleLink");
        }
 
        m_tree.addBranch(m_recoSegs);
        m_tree.addBranch(std::make_unique<EventInfoBranch>(m_tree, evOpts));
 
        static const std::vector<std::string> trackSummaries{
                 // Inner
                "innerSmallHits", "innerLargeHits", "innerSmallHoles", "innerLargeHoles",
                "innerClosePrecisionHits",
                // Middle
                "middleSmallHits", "middleLargeHits", "middleSmallHoles",
                "middleLargeHoles", "middleClosePrecisionHits",
                // Outer
                "outerSmallHits", "outerLargeHits", "outerSmallHoles", "outerLargeHoles",
                "outerClosePrecisionHits",
                // Extended
                "extendedSmallHits", "extendedLargeHits", "extendedSmallHoles",
                "extendedLargeHoles", "extendedClosePrecisionHits",
                "innerTriggerEtaHits", "innerTriggerPhiHits", 
                "middleTriggerEtaHits", "middleTriggerPhiHits",
                "outerTriggerEtaHits", "outerTriggerPhiHits", 
 
                "innerTriggerEtaHoles", "innerTriggerPhiHoles",
                "middleTriggerEtaHoles", "middleTriggerPhiHoles", 
                "outerTriggerEtaHoles", "outerTriggerPhiHoles",
                };
        if(!m_legacyTrackKey.empty()) {
            m_legacyTrks = std::make_unique<IParticleFourMomBranch>(m_tree, "LegacyMSTrks");
            m_legacyTrks->addVariable(std::make_unique<TrackChi2Branch>(*m_legacyTrks));
            if (m_isMC) {
                BilateralLinkerBranch::connectCollections(m_legacyTrks, m_truthTrks, 
                                                          [](const xAOD::IParticle* trk){ 
                                                            return xAOD::TruthHelpers::getTruthParticle(*trk); 
                                                          }, "truth", "LegacyMS");
            }
            for (const auto& summary : trackSummaries) {
                m_legacyTrks->addVariable<uint8_t>(-1, summary); 
            }
            m_tree.addBranch(m_legacyTrks);
            m_legacyRecoSegs = std::make_unique<SegmentVariables>(m_tree, m_legacySegmentKey.key(), "LegacyRecoSegments", msgLevel());
            m_tree.addBranch(m_legacyRecoSegs);
        }
 
        m_seedSummary = std::make_shared<TrackSummaryModule>(m_tree, "TrkSeed", m_summaryTool.get());
        m_muonTrks = std::make_shared<IParticleFourMomBranch>(m_tree, "ActsMuons");
        m_muonTrks->addVariable(std::make_unique<TrackChi2Branch>(*m_muonTrks));
        m_muonTrks->addVariable<float>(4.5, "segmentDeltaEta");
        m_muonTrks->addVariable<float>(std::numbers::pi, "segmentDeltaPhi");
        m_muonTrks->addVariable<float>(-1., "segmentChi2OverDoF");
 
        using TrkType = xAOD::Muon::TrackParticleType;
        auto dumpTrack = [&](const std::string& trkName, TrkType type) {
            auto trkColl = std::make_shared<IParticleFourMomBranch>(m_tree, trkName);
            trkColl->addVariable(std::make_unique<TrackChi2Branch>(*trkColl));
            if (type != TrkType::InnerDetectorTrackParticle) {
                trkColl->addVariable(std::make_unique<TrackFitIterBranch>(*trkColl));
            }
            trkColl->addVariable(std::make_unique<MaterialRecorderBranch>(*trkColl));
            trkColl->addVariable(std::make_unique<EnergyLossBranch>(*trkColl));
            
            trkColl->addVariable<float>("d0");
            trkColl->addVariable<float>("z0");
            if (type != TrkType::InnerDetectorTrackParticle) {
                for (const auto& summary : trackSummaries) {
                    trkColl->addVariable<uint8_t>(-1, summary); 
                }
            }
            BilateralLinkerBranch::connectCollections(m_muonTrks, trkColl, 
                    [type](const xAOD::IParticle* muonP) -> const xAOD::IParticle* {
                        const auto* muon = dynamic_cast<const xAOD::Muon*>(muonP);
                        if (!muon) {
                            return nullptr;
                        }
                        return muon->trackParticle(type);
                    }, trkName, "muon");
 
            return trkColl;
        };
 
        auto msTracks = dumpTrack("MsTrk", TrkType::MuonSpectrometerTrackParticle);
        dumpTrack("MeTrk", TrkType::ExtrapolatedMuonSpectrometerTrackParticle);
        m_idTracks = dumpTrack("IdTrk", TrkType::InnerDetectorTrackParticle);

        {
            auto& l0{m_tree.newVector<float>(std::format("{:}_caloExitLoc0", m_idTracks->name()))};
            auto& l1{m_tree.newVector<float>(std::format("{:}_caloExitLoc1", m_idTracks->name()))};
            auto& theta{m_tree.newVector<float>(std::format("{:}_caloExitTheta", m_idTracks->name()))};
            auto& phi{m_tree.newVector<float>(std::format("{:}_caloExitPhi", m_idTracks->name()))};
            auto& mom{m_tree.newVector<float>(std::format("{:}_caloExitMom", m_idTracks->name()))};
            
            m_idTracks->addVariable(std::make_unique<GenericPartDecorBranch<xAOD::TrackParticle, std::uint8_t>>(m_tree, 
                std::format("{:}_isCandidate", m_idTracks->name()), 
                [this, &l0, &l1, &theta, &phi, &mom] (const xAOD::TrackParticle& idTrack) -> std::uint8_t {
                    const MuonR4::MuonTag* tag = findBaseIdTag(idTrack, Gaudi::Hive::currentContext());
                    std::optional<Acts::BoundTrackParameters> trkPars = tag ? tag->extrapolatedParsID(Acts::hashString("@CaloExit"))
                                                                            : std::nullopt;
                    const std::size_t idx = l0.size();
                    using namespace MuonCombinedR4;
                    l0[idx] = trkPars ? longitudinalParam(*trkPars, logger()) : 1._km;
                    l1[idx] = trkPars ? toDeg(localPolarAngle(*trkPars, logger())) : 360.;
                    
                    theta[idx] = trkPars ? toDeg(trkPars->get<Acts::eBoundTheta>()) : 360.;
                    phi[idx] = trkPars ? toDeg(trkPars->get<Acts::eBoundPhi>()) : 360.;
                    mom[idx] = trkPars ? trkPars->absoluteMomentum() : -1._GeV;
                    return trkPars != std::nullopt;
                }));
        }
        {
            auto& dX0{m_tree.newMatrix<float>(std::format("{:}_segTagDx0", m_idTracks->name()))};
            auto& dY0{m_tree.newMatrix<float>(std::format("{:}_segTagDy0", m_idTracks->name()))};
            auto& dTheta{m_tree.newMatrix<float>(std::format("{:}_segTagDtheta", m_idTracks->name()))};
            auto& dPhi{m_tree.newMatrix<float>(std::format("{:}_segTagDphi", m_idTracks->name()))};
            auto& segments{m_tree.newMatrix<std::uint8_t>(std::format("{:}_segTagSegments", m_idTracks->name()))};
            auto& goodExtp{m_tree.newMatrix<std::uint8_t>(std::format("{:}_segTagExtpGood",m_idTracks->name()))};
            auto& taggedSeg{m_tree.newMatrix<std::uint8_t>(std::format("{:}_segTagOnTag", m_idTracks->name()))};
            m_idTracks->addVariable(
            std::make_unique<GenericPartDecorBranch<xAOD::TrackParticle, std::vector<float>>>(m_tree, 
                std::format("{:}_segTagScore", m_idTracks->name()), 
                    [&dX0, &dY0, &dTheta, &dPhi, &segments, 
                     &goodExtp, &taggedSeg,this] (const xAOD::TrackParticle& idTrack) -> std::vector<float> {
                    auto result = calcSegTagVariables(idTrack, Gaudi::Hive::currentContext());
 
                    segments.push_back(std::move(result.recoSegs));
                    dX0.push_back(std::move(result.deltaX0));
                    dY0.push_back(std::move(result.deltaY0));
                    dTheta.push_back(std::move(result.deltaTheta));
                    dPhi.push_back(std::move(result.deltaPhi));
                    goodExtp.push_back(std::move(result.goodExtp));
                    taggedSeg.push_back(std::move(result.taggedSeg));
                    return std::move(result.matchScores);
            }));
        }
        if (m_isMC) {
            BilateralLinkerBranch::connectCollections(m_idTracks, m_truthTrks, 
                [this](const xAOD::IParticle* trk) { 
                    return truthTreeParticle(*trk); 
                }, "truth", "IdTrack");
 
        }
        
        m_muonTrks->addVariable<uint16_t>("allAuthors");
        m_muonTrks->addVariable<uint16_t>("author");
        m_muonTrks->addVariable(std::make_unique<GenericPartDecorBranch<xAOD::Muon,
                                                        std::vector<unsigned short>>>(m_tree, 
                    std::format("{:}_segmentLinks", m_muonTrks->name()), [&] (const xAOD::Muon& p){
                    std::vector<unsigned short> idx{};
                    for (unsigned seg = 0 ; seg < p.nMuonSegments(); ++seg) {
                        idx.push_back(m_recoSegs->push_back(*p.muonSegment(seg)));
                    }
                    return idx;
                }));
        msTracks->addVariable(std::make_unique<GenericPartDecorBranch<xAOD::TrackParticle, unsigned short>>(m_tree, 
                    std::format("{:}_seedLink", msTracks->name()), [&] (const xAOD::TrackParticle& p) -> unsigned short {
                        auto actsTrk = ActsTrk::getActsTrack(p);
                        if (!actsTrk) {
                            THROW_EXCEPTION("Cannot find the associated ms track from the primary track");
                        }
                        return actsTrk->component<std::size_t>("parentSeed");
                }));
 
        msTracks->addVariable(std::make_unique<GenericPartDecorBranch<xAOD::TrackParticle, float>>(m_tree, 
                    std::format("{:}_loc0", msTracks->name()), [&] (const xAOD::TrackParticle& p){
                        auto actsTrk = ActsTrk::getActsTrack(p);
                        if (!actsTrk) {
                            THROW_EXCEPTION("Cannot find the associated ms track from the primary track");
                        }
                        auto refPars = actsTrk->createParametersAtReference();
                        return isOnCaloExit(refPars) ? MuonCombinedR4::longitudinalParam(refPars, logger())
                                                     : refPars.get<Acts::eBoundLoc0>();
                }));
        msTracks->addVariable(std::make_unique<GenericPartDecorBranch<xAOD::TrackParticle, float>>(m_tree, 
                    std::format("{:}_loc1", msTracks->name()), [&] (const xAOD::TrackParticle& p){
                        auto actsTrk = ActsTrk::getActsTrack(p);
                        if (!actsTrk) {
                            THROW_EXCEPTION("Cannot find the associated ms track from the primary track");
                        }
                        auto refPars = actsTrk->createParametersAtReference();
                        return isOnCaloExit(refPars) ? toDeg(MuonCombinedR4::localPolarAngle(refPars, logger()))
                                                     : refPars.get<Acts::eBoundLoc0>(); 
                }));
        msTracks->addVariable(std::make_unique<GenericPartDecorBranch<xAOD::TrackParticle, std::uint8_t>>(m_tree, 
                    std::format("{:}_isOnCaloExit", msTracks->name()), [&] (const xAOD::TrackParticle& p){
                        auto actsTrk = ActsTrk::getActsTrack(p);
                        if (!actsTrk) {
                            THROW_EXCEPTION("Cannot find the associated ms track from the primary track");
                        }
                        return isOnCaloExit(actsTrk->referenceSurface());
                }));
        if (m_isMC) {
            BilateralLinkerBranch::connectCollections(m_muonTrks, m_truthTrks, 
                [](const xAOD::IParticle* trk) -> const xAOD::TruthParticle* { 
                    return xAOD::TruthHelpers::getTruthParticle(*trk); 
                }, "truth", "ActsMuon");
        }
 
        for (const auto& summary : trackSummaries) {
             m_muonTrks->addVariable<uint8_t>(summary); 
        }
        m_tree.addBranch(m_muonTrks);
 
        m_tree.addBranch(m_seedSummary);
        
        ATH_CHECK(m_trkTruthLinks.initialize());
        ATH_CHECK(m_tree.init(this));
        return StatusCode::SUCCESS;
    }

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.

                                           {
    return true;
  }

isReEntrant()

virtual bool AthAlgorithm::isReEntrant ( ) const

inlinefinaloverrideprotectedvirtualinherited

Legacy algorithms are not thread-safe.

Definition at line 47 of file AthAlgorithm.h.

{ return false; }

logger()

const Acts::Logger & MuonValR4::MsTrackTester::logger ( ) const

inlineprivate

Return the reference to the Acts logger.

Definition at line 194 of file MsTrackTester.h.

{ return *m_logger; }

makeSeedFromTruth()

std::optional< MsTrackSeed > MuonValR4::MsTrackTester::makeSeedFromTruth ( const Acts::GeometryContext & tgContext, 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

gctx	The geometry context to align the segment within ATLAS
truthMuon	Reference to the truth muon for which a seed should be constructed

Construct the 2 candidate seeds

Definition at line 68 of file MsTrackTester.cxx.

                                                                                                          {
        std::vector<const xAOD::MuonSegment*> matchedSegs = MuonR4::getTruthSegments(truthMuon);
        if (matchedSegs.empty() || toLayerIndex(matchedSegs.front()->chamberIndex()) == toLayerIndex(matchedSegs.back()->chamberIndex())) {
            return std::nullopt;
        }
        ExpandedSector sector{matchedSegs[0]->position().phi()};
        MsTrackSeed barrelSeed{Location::Barrel, sector};
        MsTrackSeed endcapSeed{Location::Endcap, sector};
        for (const xAOD::MuonSegment* seg : matchedSegs) {
            barrelSeed.addSegment(seg);
            endcapSeed.addSegment(seg);
        }
        barrelSeed.setPosition(matchedSegs[0]->position());
        endcapSeed.setPosition(matchedSegs[0]->position());
        const auto [barrelLength, barrelTheta] = calcSeedLength(tgContext, barrelSeed);
        const auto [endcapLength, endcapTheta] = calcSeedLength(tgContext, endcapSeed);
        ATH_MSG_VERBOSE(__func__<<"() "<<__LINE__<<" - Constructed new seed from truth muon wih pT:"
                <<(truthMuon.pt()/ Gaudi::Units::GeV)<<" [GeV], eta: "<<truthMuon.eta()
                <<", phi: "<<toDeg(truthMuon.phi())<<", q: "<<truthMuon.charge()
                <<", matchedSeg: "<<matchedSegs.size()<< " barrel (L/theta): "<<barrelLength
                <<"/"<<barrelTheta<<" - endcap (L/theta): "
                <<endcapLength<<"/"<<endcapTheta<<"\n"<<barrelSeed);
        if (barrelLength < 0 && endcapLength < 0) {
            ATH_MSG_WARNING(__func__<<"() "<<__LINE__<<" - Invalid seed");
            return std::nullopt;
        }
        return barrelLength < 0 || std::abs(endcapLength) < barrelLength 
              ? endcapSeed : barrelSeed;
    }

msg()

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

msgLvl()

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

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

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.

{
  size_t pos = 0;
  while((pos = str.find(oldStr, pos)) != std::string::npos)
    {
      str.replace(pos, oldStr.length(), newStr);
      pos += newStr.length();
    }
}

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.

  {
    h.renounce();
    PBASE::renounce (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.

                                                          {
    handlesArray.renounce();
  }

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.

                                                                            {
    execState( ctx ).setFilterPassed( state );
  }

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.

{
  return BaseAlg::sysExecute (ctx);
}

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.

{
  // ---- stolen from GaudiKernel/Algorithm::sysInitialize -------
  // Bypass the initialization if the algorithm
  // has already been initialized.
  if ( Gaudi::StateMachine::INITIALIZED <= FSMState() ) return StatusCode::SUCCESS;
 
  // Set the Algorithm's properties
  bindPropertiesTo( serviceLocator()->getOptsSvc() );
 
  // Bypass the initialization if the algorithm is disabled.
  // Need to do this after bindPropertiesTo.
  if ( !isEnabled( ) ) return StatusCode::SUCCESS;
  
  // ---- stolen from GaudiKernel/Algorithm::sysInitialize ------- END ---
 
 
  // Get the THistSvc
  ATH_CHECK ( histSvc().retrieve() );
 
  // Configure the underlying AthHistogramming helper
  ATH_CHECK ( AthHistogramming::configAthHistogramming( m_histSvc,
                                                        m_prefix,          m_rootDir, 
                                                        m_histNamePrefix,  m_histNamePostfix,
                                                        m_histTitlePrefix, m_histTitlePostfix ) );
  
  // Print some setup information into the log file
  ATH_MSG_DEBUG ("Initializing " << name() << "...");
  ATH_MSG_DEBUG (" using THistService     = " << m_histSvc );
  ATH_MSG_DEBUG (" using RootStreamName   = " << m_prefix );
  ATH_MSG_DEBUG (" using RootDirName      = " << m_rootDir );
  ATH_MSG_DEBUG (" using HistNamePrefix   = " << m_histNamePrefix );
  ATH_MSG_DEBUG (" using HistNamePostfix  = " << m_histNamePostfix );
  ATH_MSG_DEBUG (" using HistTitlePrefix  = " << m_histTitlePrefix );
  ATH_MSG_DEBUG (" using HistTitlePostfix = " << m_histTitlePostfix );
 
 
  // re-direct to base class...
  return AthAlgorithm::sysInitialize();
}

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.

{
  // Build a 32 bit hash out of the name
  const hash_t treeHash = this->hash(treeName);
 
  // See if this entry exists in the map
  TreeMap_t::const_iterator it = m_treeMap.find( treeHash );
  if ( it == m_treeMap.end() ) // It doesn't exist!
    { // Let's see into the THistSvc if somebody else has registered the TTree...
 
      // Need to copy the strings as we will massage them from here on
      std::string treeNameCopy = treeName;
      std::string tDirCopy     = tDir;
      std::string streamCopy   = stream;
 
      // Massage the final string to book things
      std::string bookingString("");
      this->buildBookingString( bookingString, treeNameCopy, tDirCopy, streamCopy );
 
      TTree* treePointer(NULL);
      if ( !((histSvc()->getTree(bookingString, treePointer)).isSuccess()) )
        {
          m_msg << MSG::WARNING
                << "Problem retrieving the TTree with name " << treeNameCopy
                << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
                << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
          return NULL;
        }
      // If we get to here, we actually found the TTree in the THistSvc.
      // So let's add it to the local cache map and return its pointer
      m_treeMap.insert( m_treeMap.end(), std::pair< const hash_t, TTree* >( treeHash, treePointer ) );
      return treePointer;
    }
 
  // Return the pointer to the TTree that we got from the local cache map
  return it->second;
}

truthTreeParticle()

const xAOD::TruthParticle * MuonValR4::MsTrackTester::truthTreeParticle ( const xAOD::IParticle & part ) const

private

Returns the matched truth particle that is dumped in the tree.

This method only works downstream the dumpTruthContent call

Parameters

part	Reference to the particle of interest

Definition at line 481 of file MsTrackTester.cxx.

                                                                                             {
        const xAOD::TruthParticle* truthTrk = xAOD::TruthHelpers::getTruthParticle(part);
        if (!truthTrk || !m_truthTrks) {
            return nullptr;
        }
        for (const xAOD::IParticle* thisPart : m_truthTrks->getCached()) {
            if (thisPart == truthTrk || truthTrk ==  xAOD::TruthHelpers::getTruthParticle(*thisPart)) {
                return dynamic_cast<const xAOD::TruthParticle*>(thisPart);
            }
        }
        return nullptr;
    }

updateVHKA()

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

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

                                               {
    // debug() << "updateVHKA for property " << p.name() << " " << p.toString() 
    //         << "  size: " << m_vhka.size() << endmsg;
    for (auto &a : m_vhka) {
      std::vector<SG::VarHandleKey*> keys = a->keys();
      for (auto k : keys) {
        k->setOwner(this);
      }
    }
  }

Member Data Documentation

m_ctxProvider

ActsTrk::ContextUtility MuonValR4::MsTrackTester::m_ctxProvider {this}

private

Context provider for geometry, magnetic field and calibration contexts.

Definition at line 155 of file MsTrackTester.h.

{this};

m_detMgr

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

private

The detector manager to fetch the segment surfaces.

Definition at line 126 of file MsTrackTester.h.

{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_extrapolationTool

ToolHandle<ActsTrk::IExtrapolationTool> MuonValR4::MsTrackTester::m_extrapolationTool {this, "ExtrapolationTool" ,"" }

private

Track extrapolation tool.

Definition at line 165 of file MsTrackTester.h.

{this, "ExtrapolationTool" ,"" };

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_idTagKey

MuonTagKey_t MuonValR4::MsTrackTester::m_idTagKey {this, "IdTagKey", "MuonInDetCandidates"}

private

The collection of all ID tracks that are selected for combined reconstruction.

Definition at line 171 of file MsTrackTester.h.

{this, "IdTagKey", "MuonInDetCandidates"};

m_idTrackKey

TrackKey_t MuonValR4::MsTrackTester::m_idTrackKey {this, "IdTrackKey", "InDetTrackParticles"}

private

The collection of ID tracks associated with the truth particle.

Definition at line 169 of file MsTrackTester.h.

{this, "IdTrackKey", "InDetTrackParticles"};

m_idTracks

ParticleBranchPtr_t MuonValR4::MsTrackTester::m_idTracks {}

private

Stored ID track information.

Definition at line 202 of file MsTrackTester.h.

{};

m_isMC

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

private

Toggle whether the job runs on MC or not.

Definition at line 128 of file MsTrackTester.h.

{this, "isMC", false};

m_legacyMuonKey

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

private

Legacy muons.

Definition at line 189 of file MsTrackTester.h.

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

m_legacyRecoSegs

SegmentBranchPtr_t MuonValR4::MsTrackTester::m_legacyRecoSegs {}

private

Definition at line 208 of file MsTrackTester.h.

{};

m_legacySegmentKey

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

private

Legacy segment container.

Definition at line 146 of file MsTrackTester.h.

{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 239 of file MsTrackTester.h.

{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 167 of file MsTrackTester.h.

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

m_legacyTrks

ParticleBranchPtr_t MuonValR4::MsTrackTester::m_legacyTrks {}

private

Output branches of the legacy MS tracks.

Definition at line 249 of file MsTrackTester.h.

{};

m_logger

std::unique_ptr<const Acts::Logger> MuonValR4::MsTrackTester::m_logger {}

private

Instance to the Acts logger.

Definition at line 192 of file MsTrackTester.h.

{};

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 152 of file MsTrackTester.h.

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

m_muonKey

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

private

Dependency on the R4 muon container.

Definition at line 157 of file MsTrackTester.h.

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

m_muonTrks

ParticleBranchPtr_t MuonValR4::MsTrackTester::m_muonTrks {}

private

Stored muon information from the Acts muon reco chain.

Definition at line 200 of file MsTrackTester.h.

{};

m_name

std::string AthHistogramming::m_name

privateinherited

Instance name.

Definition at line 247 of file AthHistogramming.h.

m_nIdTags

MuonVal::ScalarBranch<std::uint16_t>& MuonValR4::MsTrackTester::m_nIdTags {m_tree.newScalar<std::uint16_t>("nIdTags", 0)}

private

The number of selected ID tracks in the event.

Definition at line 212 of file MsTrackTester.h.

{m_tree.newScalar<std::uint16_t>("nIdTags", 0)};

m_nIdTracks

MuonVal::ScalarBranch<std::uint16_t>& MuonValR4::MsTrackTester::m_nIdTracks {m_tree.newScalar<std::uint16_t>("nIdTracks", 0)}

private

The number of ID tracks in the event.

Definition at line 210 of file MsTrackTester.h.

{m_tree.newScalar<std::uint16_t>("nIdTracks", 0)};

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 144 of file MsTrackTester.h.

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

m_recoSegs

SegmentBranchPtr_t MuonValR4::MsTrackTester::m_recoSegs {}

private

Definition at line 207 of file MsTrackTester.h.

{};

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, "TrkSeed_direction"}

private

Seed direction vector.

Definition at line 216 of file MsTrackTester.h.

{m_tree, "TrkSeed_direction"};

m_seedGood

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

private

Does the seeding tool construct valid parameters from the seed.

Definition at line 226 of file MsTrackTester.h.

{m_tree.newVector<char>("TrkSeed_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 161 of file MsTrackTester.h.

{this, "SeedingTool", ""};

m_seedLength

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

private

Maximum separation between the segments on the reference plane.

Definition at line 222 of file MsTrackTester.h.

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

m_seedPos

MuonVal::ThreeVectorBranch MuonValR4::MsTrackTester::m_seedPos {m_tree, "TrkSeed_position"}

private

Simple seed information.

Definition at line 214 of file MsTrackTester.h.

{m_tree, "TrkSeed_position"};

m_seedQP

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

private

Estimated momentum times charge from the track seed.

Definition at line 228 of file MsTrackTester.h.

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

m_seedRecoSegMatch

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

private

Link of the track seed to the building segment.

Definition at line 234 of file MsTrackTester.h.

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

m_seedSector

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

private

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

Definition at line 220 of file MsTrackTester.h.

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

m_seedSummary

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

private

Hit summary on the track seed.

Definition at line 232 of file MsTrackTester.h.

{};

m_seedThetaCone

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

private

Maximum angular difference between the segments part of the seed.

Definition at line 224 of file MsTrackTester.h.

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

m_seedTruthLink

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

private

Link to the truth muon.

Definition at line 230 of file MsTrackTester.h.

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

m_seedType

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

private

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

Definition at line 218 of file MsTrackTester.h.

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

m_segSelector

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

private

Selection tool to quantify the segment candidate quality.

Definition at line 163 of file MsTrackTester.h.

{this, "SegmentSelectionTool" , "" };

m_segTagKey

MuonTagKey_t MuonValR4::MsTrackTester::m_segTagKey {this, "SegTagKey", "SegmentTags"}

private

The collection of the original segment tags.

Definition at line 173 of file MsTrackTester.h.

{this, "SegTagKey", "SegmentTags"};

m_storeID

Gaudi::Property<bool> MuonValR4::MsTrackTester::m_storeID {this, "storeIdTrks", true}

private

Toggle whether to process ID tracks or not.

Definition at line 130 of file MsTrackTester.h.

{this, "storeIdTrks", true};

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 159 of file MsTrackTester.h.

{this, "SummaryTool" ,""};

m_tolerancePhi

Gaudi::Property<double> MuonValR4::MsTrackTester::m_tolerancePhi {this, "tolerancePhi", 2.*Gaudi::Units::deg}

private

Extra tolerance applied on the bending direction when calculating the matching score.

Definition at line 186 of file MsTrackTester.h.

{this, "tolerancePhi", 2.*Gaudi::Units::deg};

m_toleranceTheta

Gaudi::Property<double> MuonValR4::MsTrackTester::m_toleranceTheta {this, "toleranceTheta", 1.*Gaudi::Units::deg}

private

Extra tolerance applied on the bending direction when calculating the matching score.

Definition at line 183 of file MsTrackTester.h.

{this, "toleranceTheta", 1.*Gaudi::Units::deg};

m_toleranceX0

Gaudi::Property<double> MuonValR4::MsTrackTester::m_toleranceX0 {this, "toleranceX0", 20.*Gaudi::Units::cm}

private

Extra tolerance applied on the non-bending intercept when calculating the matching score.

Definition at line 177 of file MsTrackTester.h.

{this, "toleranceX0", 20.*Gaudi::Units::cm};

m_toleranceY0

Gaudi::Property<double> MuonValR4::MsTrackTester::m_toleranceY0 {this, "toleranceY0", 5.*Gaudi::Units::cm}

private

Extra tolerance applied on the bending intercept when calculating the matching score.

Definition at line 180 of file MsTrackTester.h.

{this, "toleranceY0", 5.*Gaudi::Units::cm};

m_tree

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

private

The output tree object.

Definition at line 124 of file MsTrackTester.h.

{"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 150 of file MsTrackTester.h.

{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 148 of file MsTrackTester.h.

{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 246 of file MsTrackTester.h.

{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 244 of file MsTrackTester.h.

{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 242 of file MsTrackTester.h.

{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 142 of file MsTrackTester.h.

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

m_truthSegs

SegmentBranchPtr_t MuonValR4::MsTrackTester::m_truthSegs {}

private

Definition at line 206 of file MsTrackTester.h.

{};

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 236 of file MsTrackTester.h.

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

m_truthTrks

ParticleBranchPtr_t MuonValR4::MsTrackTester::m_truthTrks {}

private

Definition at line 198 of file MsTrackTester.h.

{};

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:

• MsTrackTester.h
• MsTrackTester.cxx