ATLAS Offline Software
Loading...
Searching...
No Matches
ActsTrk::ActsMuonTrackingGeometryTest Class Reference

Extrapolation test for the ActsMuonTrackingGeometry for gen3. More...

#include <ActsMuonTrackingGeometryTest.h>

Inheritance diagram for ActsTrk::ActsMuonTrackingGeometryTest:

Public Member Functions

 ~ActsMuonTrackingGeometryTest ()=default
StatusCode initialize () override
StatusCode execute () override
StatusCode finalize () override
 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 const DataObjIDColl & extraOutputDeps () const override
 Return the list of extra output dependencies.
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

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

Private Types

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

Private Member Functions

Amg::Transform3D toLocalTrf (const ActsTrk::GeometryContext &gctx, const Identifier &hitId) const
Amg::Transform3D toGlobalTrf (const ActsTrk::GeometryContext &gctx, const Identifier &hitId) const
IdentifierHash layerHash (const Identifier &id) const
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey>
void buildBookingString (std::string &bookingString, std::string &histName, std::string &tDir, std::string &stream, bool usePrefixPostfix=false)
 Method to build individual booking string.
void myReplace (std::string &str, const std::string &oldStr, const std::string &newStr)
 Helper method to replace sub-string.
hash_t hash (const std::string &histName) const
 Method to calculate a 32-bit hash from a string.

Private Attributes

const MuonGMR4::MuonDetectorManagerm_r4DetMgr {nullptr}
ServiceHandle< Muon::IMuonIdHelperSvcm_idHelperSvc {this, "IdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}
ServiceHandle< IAthRNGSvcm_rndmGenSvc {this, "AthRNGSvc", "AthRNGSvc"}
PublicToolHandle< ActsTrk::ITrackingGeometryToolm_trackingGeometryTool {this, "TrackingGeometryTool", "ActsTrackingGeometryTool"}
SG::ReadCondHandleKey< MuonGM::MuonDetectorManagerm_detMgrKey
SG::ReadCondHandleKey< AtlasFieldCacheCondObjm_fieldCacheCondObjInputKey
SG::ReadHandleKey< ActsTrk::GeometryContextm_geoCtxKey {this, "AlignmentKey", "ActsAlignment", "cond handle key"}
SG::ReadHandleKey< xAOD::TruthParticleContainerm_truthParticleKey {this, "TruthKey", "MuonTruthParticles"}
SG::ReadDecorHandleKey< xAOD::TruthParticleContainerm_truthSegLinkKey {this, "TruthKeyToSeg", m_truthParticleKey, "truthSegmentLinks", "TruthParticle to TruthSegment link"}
UnsignedIntegerProperty m_nEvents {this, "nEvents", 100}
DoubleProperty m_pathLimit {this, "PathLimit", 30.*Gaudi::Units::m}
DoubleProperty m_maxStepSize {this, "MaxStepSize", 1*Gaudi::Units::m}
UnsignedIntegerProperty m_maxSteps {this, "MaxSteps", 100000}
UnsignedIntegerProperty m_maxTargetSkipping {this, "MaxTargetSkipping", 10000}
DoubleProperty m_stepTolerance {this, "StepTolerance", 1e-10}
Gaudi::Property< bool > m_startFromFirstHit {this, "StartFromFirstHit", false, "Start from first hit"}
Gaudi::Property< std::vector< double > > m_etaRange {this, "EtaRange", {-3, 3}, "The eta range for particles"}
MuonVal::MuonTesterTree m_tree {"MuonNavigationTestGen3R4", "MuonNavigationTestGen3R4"}
MuonVal::MuonIdentifierBranch m_detId {m_tree, "detId"}
MuonVal::VectorBranch< unsigned short > & m_techIdx {m_tree.newVector<unsigned short>("detId_techIdx")}
MuonVal::VectorBranch< unsigned short > & m_gasGapId {m_tree.newVector<unsigned short>("detId_gasGap")}
MuonVal::ThreeVectorBranch m_truthLoc {m_tree, "truthLoc"}
MuonVal::ThreeVectorBranch m_truthGlob {m_tree, "truthGlob"}
MuonVal::ThreeVectorBranch m_truthDir {m_tree, "truthDir"}
MuonVal::ScalarBranch< float > & m_truthPt {m_tree.newScalar<float>("truthPt")}
MuonVal::ScalarBranch< float > & m_truthP {m_tree.newScalar<float>("truthP")}
MuonVal::ThreeVectorBranch m_startGlob {m_tree, "startGlob"}
MuonVal::VectorBranch< float > & m_actsPropabsMomentum {m_tree.newVector<float>("actsPropabsMomentum")}
MuonVal::ThreeVectorBranch m_actsPropLoc {m_tree, "actsPropLoc"}
MuonVal::ThreeVectorBranch m_actsPropGlob {m_tree, "actsPropGlob"}
MuonVal::ThreeVectorBranch m_actsPropDir {m_tree, "actsPropDir"}
MuonVal::VectorBranch< float > & m_actsHitWireDist {m_tree.newVector<float>("actsHitWireDist")}
MuonVal::VectorBranch< float > & m_actsStepSize {m_tree.newVector<float>("actsStepSize")}
MuonVal::ScalarBranch< float > & m_propLength {m_tree.newScalar<float>("propLength")}
MuonVal::VectorBranch< unsigned short > & m_isPropagated {m_tree.newVector<unsigned short>("isPropagated")}
MuonVal::ScalarBranch< unsigned int > & m_propSteps {m_tree.newScalar<unsigned int>("propSteps")}
MuonVal::ScalarBranch< unsigned int > & m_event {m_tree.newScalar<unsigned int>("event")}
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
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

Extrapolation test for the ActsMuonTrackingGeometry for gen3.

Definition at line 37 of file ActsMuonTrackingGeometryTest.h.

Member Typedef Documentation

◆ EffMap_t

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

Typedef for convenience.

Definition at line 205 of file AthHistogramming.h.

◆ GraphMap_t

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

Typedef for convenience.

Definition at line 219 of file AthHistogramming.h.

◆ hash_t

typedef uint32_t AthHistogramming::hash_t
privateinherited

typedef for the internal hash

Definition at line 169 of file AthHistogramming.h.

◆ HistMap_t

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

Typedef for convenience.

Definition at line 198 of file AthHistogramming.h.

◆ StoreGateSvc_t

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

Definition at line 388 of file AthCommonDataStore.h.

◆ TreeMap_t

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

Typedef for convenience.

Definition at line 212 of file AthHistogramming.h.

Constructor & Destructor Documentation

◆ ~ActsMuonTrackingGeometryTest()

ActsTrk::ActsMuonTrackingGeometryTest::~ActsMuonTrackingGeometryTest ( )
default

Member Function Documentation

◆ AthHistogramAlgorithm()

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

Constructor with parameters:

Definition at line 40 of file AthHistogramAlgorithm.cxx.

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

◆ book() [1/8]

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

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

Definition at line 335 of file AthHistogramming.h.

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

◆ book() [2/8]

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

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

◆ book() [3/8]

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

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

Definition at line 303 of file AthHistogramming.h.

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

◆ book() [4/8]

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

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

Definition at line 403 of file AthHistogramming.h.

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

◆ book() [5/8]

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

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

Definition at line 356 of file AthHistogramming.h.

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

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

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

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

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

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

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

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

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

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

428{
429 // Get a pointer
430 const TGraph* graphPointer = &graphRef;
431
432 // Check that we got a valid pointer
433 if ( !graphPointer )
434 {
435 m_msg << MSG::WARNING
436 << "We got an invalid TGraph pointer in the BookGetPointer(TGraph*) method of the class" << m_name
437 << "!" << endmsg;
438 return NULL;
439 }
440
441 // Modify the name and title according to the prefixes of this classes instance
442 std::string graphName = graphPointer->GetName();
443 const std::string graphTitle = graphPointer->GetTitle();
444
445 // Check if the hash for this graphName already exists, i.e., if we have a hash collision
446 const hash_t graphHash = this->hash(graphName);
447 GraphMap_t::const_iterator it = m_graphMap.find( graphHash );
448 if ( it != m_graphMap.end() ) // It does exist!
449 {
450 m_msg << MSG::WARNING
451 << "Detected a hash collision. The hash for the TGraph with name=" << graphName
452 << " already exists and points to a TGraph with name=" << it->second->GetName()
453 << " NOT going to book the new histogram and returning a NULL pointer!" << endmsg;
454 return NULL;
455 }
456
457 // Create a clone that has the new name
458 TGraph* graphClone = dynamic_cast< TGraph* >( graphPointer->Clone((m_histNamePrefix+graphName+m_histNamePostfix).c_str()) );
459 if( !graphClone )
460 {
461 m_msg << MSG::WARNING
462 << "We couldn't clone the TGraph in the BookGetPointer(TGraph&) method of the class" << m_name
463 << "!" << endmsg;
464 return NULL;
465 }
466 graphClone->SetTitle ((m_histTitlePrefix+graphTitle+m_histTitlePostfix).c_str());
467
468 // Massage the final string to book things
469 std::string bookingString("");
470 this->buildBookingString( bookingString, graphName, tDir, stream );
471
472 // Register the TGraph into the THistSvc
473 if ( !((histSvc()->regGraph(bookingString, graphClone)).isSuccess()) )
474 {
475 m_msg << MSG::WARNING
476 << "Problem registering TGraph with name " << graphName
477 << ", title " << graphTitle
478 << " in " << m_name << "!" << endmsg;
479 return NULL;
480 }
481
482 // Also register it in the local map of string to pointer
483 m_graphMap.insert( m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphClone ) );
484
485 return graphClone;
486}
const ServiceHandle< ITHistSvc > & histSvc() const
The standard THistSvc (for writing histograms and TTrees and more to a root file) Returns (kind of) a...
std::string m_histNamePostfix
The postfix for the histogram THx name.
hash_t hash(const std::string &histName) const
Method to calculate a 32-bit hash from a string.
uint32_t hash_t
typedef for the internal hash
std::string m_histTitlePostfix
The postfix for the histogram THx title.
std::string m_histTitlePrefix
The prefix for the histogram THx title.
std::string m_name
Instance name.
void buildBookingString(std::string &bookingString, std::string &histName, std::string &tDir, std::string &stream, bool usePrefixPostfix=false)
Method to build individual booking string.
std::string m_histNamePrefix
The prefix for the histogram THx name.
GraphMap_t m_graphMap
The map of TGraph names to their pointers.

◆ bookGetPointer() [3/8]

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

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

Definition at line 260 of file AthHistogramming.h.

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

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

313{
314 // Get a pointer
315 const TTree* treePointer = &treeRef;
316
317 // Check that we got a valid pointer
318 if ( !treePointer )
319 {
320 m_msg << MSG::WARNING
321 << "We got an invalid TTree pointer in the BookGetPointer(TTree*) method of the class" << m_name
322 << "!" << endmsg;
323 return NULL;
324 }
325
326 // Modify the name and title according to the prefixes of this classes instance
327 std::string treeName = treePointer->GetName();
328 const std::string treeTitle = treePointer->GetTitle();
329
330 // Check if the hash for this treeName already exists, i.e., if we have a hash collision
331 const hash_t treeHash = this->hash(treeName);
332 TreeMap_t::const_iterator it = m_treeMap.find( treeHash );
333 if ( it != m_treeMap.end() ) // It does exist!
334 {
335 m_msg << MSG::WARNING
336 << "Detected a hash collision. The hash for the TTree with name=" << treeName
337 << " already exists and points to a TTree with name=" << it->second->GetName()
338 << " NOT going to book the new histogram and returning a NULL pointer!" << endmsg;
339 return NULL;
340 }
341
342 // Create a clone that has the new name
343 TTree* treeClone = dynamic_cast< TTree* >( treePointer->Clone(treeName.c_str()) );
344 if( !treeClone )
345 {
346 m_msg << MSG::WARNING
347 << "We couldn't clone the TTree in the BookGetPointer(TTree&) method of the class" << m_name
348 << "!" << endmsg;
349 return NULL;
350 }
351 treeClone->SetTitle (treeTitle.c_str());
352
353 // Massage the final string to book things
354 std::string bookingString("");
355 this->buildBookingString( bookingString, treeName, tDir, stream );
356
357 // Register the TTree into the THistSvc
358 if ( !((histSvc()->regTree(bookingString, treeClone)).isSuccess()) )
359 {
360 m_msg << MSG::WARNING
361 << "Problem registering TTree with name " << treeName
362 << ", title " << treeTitle
363 << " in " << m_name << "!" << endmsg;
364 return NULL;
365 }
366
367 // Also register it in the local map of string to pointer
368 m_treeMap.insert( m_treeMap.end(), std::pair< const hash_t, TTree* >( treeHash, treeClone ) );
369
370 return treeClone;
371}
TreeMap_t m_treeMap
The map of TTree names to their pointers.

◆ bookGetPointer() [5/8]

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

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

Definition at line 146 of file AthHistogramming.cxx.

147{
148 // Modify the name and title according to the prefixes of this classes instance
149 std::string effName(effRef.GetName());
150 const std::string effTitle(effRef.GetTitle());
151 std::string bookingString("");
152
153 this->buildBookingString( bookingString, effName, tDir, stream );
154 effRef.SetTitle((m_histTitlePrefix+effTitle+m_histTitlePostfix).c_str() );
155 effRef.SetName(effName.c_str());
156
157 // Check if the hash for this effName already exists, i.e., if we have a hash collision
158 const hash_t effHash = this->hash(effName);
159 EffMap_t::const_iterator it = m_effMap.find( effHash );
160 if ( it != m_effMap.end() ) // It does exist!
161 {
162 m_msg << MSG::WARNING
163 << "Detected a hash collision. The hash for the TEfficiency with name=" << effName
164 << " already exists and points to a TEfficiency with name=" << it->second->GetName()
165 << " NOT going to book the new TEfficiency and returning a NULL pointer!" << endmsg;
166 return NULL;
167 }
168
169 // Set the new name and title for the TEfficiency, based on the prefixes that the user set for this class instance
170 // Create a clone that has the new name
171
172 // Massage the final string to book things
173
174 // Register the TEfficiency into the THistSvc
175 if ( !((histSvc()->regEfficiency(bookingString, &effRef)).isSuccess()) )
176 {
177 m_msg << MSG::WARNING
178 << "Problem registering TEfficiency with name " << effName
179 << ", name prefix " << m_histNamePrefix
180 << ", title " << effTitle
181 << ", tile prefix " << m_histTitlePrefix
182 << ", and tile postfix " << m_histTitlePostfix
183 << " in " << m_name << "!" << endmsg;
184 return NULL;
185 }
186
187 // Also register it in the local map of string to pointer
188 m_effMap.insert( m_effMap.end(), std::pair< const hash_t, TEfficiency* >( effHash, &effRef ) );
189
190 return &effRef;
191}
EffMap_t m_effMap
The map of histogram names to their pointers.

◆ bookGetPointer() [6/8]

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

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

Definition at line 293 of file AthHistogramming.h.

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

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

99{
100 // Modify the name and title according to the prefixes of this classes instance
101 std::string histName(histRef.GetName());
102 const std::string histTitle(histRef.GetTitle());
103 std::string bookingString("");
104
105 this->buildBookingString( bookingString, histName, tDir, stream );
106 histRef.SetTitle((m_histTitlePrefix+histTitle+m_histTitlePostfix).c_str() );
107 histRef.SetName(histName.c_str());
108
109 // Check if the hash for this histName already exists, i.e., if we have a hash collision
110 const hash_t histHash = this->hash(histName);
111 HistMap_t::const_iterator it = m_histMap.find( histHash );
112 if ( it != m_histMap.end() ) // It does exist!
113 {
114 m_msg << MSG::WARNING
115 << "Detected a hash collision. The hash for the histogram with name=" << histName
116 << " already exists and points to a histogram with name=" << it->second->GetName()
117 << " NOT going to book the new histogram and returning a NULL pointer!" << endmsg;
118 return NULL;
119 }
120
121 // Set the new name and title for the histogram, based on the prefixes that the user set for this class instance
122 // Create a clone that has the new name
123
124 // Massage the final string to book things
125
126 // Register the histogram into the THistSvc
127 if ( !((histSvc()->regHist(bookingString, &histRef)).isSuccess()) )
128 {
129 m_msg << MSG::WARNING
130 << "Problem registering histogram with name " << histName
131 << ", name prefix " << m_histNamePrefix
132 << ", title " << histTitle
133 << ", tile prefix " << m_histTitlePrefix
134 << ", and tile postfix " << m_histTitlePostfix
135 << " in " << m_name << "!" << endmsg;
136 return NULL;
137 }
138
139 // Also register it in the local map of string to pointer
140 m_histMap.insert( m_histMap.end(), std::pair< const hash_t, TH1* >( histHash, &histRef ) );
141
142 return &histRef;
143}
HistMap_t m_histMap
The map of histogram names to their pointers.

◆ bookGetPointer() [8/8]

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

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

Definition at line 272 of file AthHistogramming.h.

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

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

565{
566 // Massage the final string to book things
567 if(tDir.empty()) tDir = m_rootDir;
568 size_t pos = histName.rfind('/');
569 if(pos != std::string::npos){
570 tDir+='/';
571 tDir.append(histName, 0,pos);
572 histName.erase(0,pos+1);
573 };
574 if(stream.empty()) stream = m_streamName;
575
576 if(usePrefixPostfix){
577 bookingString = "/"+stream+"/"+tDir+"/"+m_histNamePrefix+histName+m_histNamePostfix;
578 } else {
579 bookingString = "/"+stream+"/"+tDir+"/"+histName;
580 }
581 while(bookingString.find("//") != std::string::npos){
582 this->myReplace(bookingString,"//","/");
583 }
584
585 return;
586}
std::string m_rootDir
Name of the ROOT directory.
std::string m_streamName
Name of the ROOT output stream (file)
void myReplace(std::string &str, const std::string &oldStr, const std::string &newStr)
Helper method to replace sub-string.

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

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

◆ declareGaudiProperty()

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

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

Definition at line 156 of file AthCommonDataStore.h.

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

◆ declareProperty()

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

Definition at line 145 of file AthCommonDataStore.h.

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

◆ detStore()

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

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

Definition at line 95 of file AthCommonDataStore.h.

◆ efficiency()

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

Simplify the retrieval of registered TEfficiency.

Definition at line 250 of file AthHistogramming.cxx.

251{
252 // Build a 32 bit hash out of the name
253 const hash_t effHash = this->hash(effName);
254
255 // See if this entry exists in the map
256 EffMap_t::const_iterator it = m_effMap.find( effHash );
257 if ( it == m_effMap.end() ) // It doesn't exist!
258 { // Let's see into the THistSvc if somebody else has registered the TEfficiency...
259
260 // Need to copy the strings as we will massage them from here on
261 std::string effNameCopy = effName;
262 std::string tDirCopy = tDir;
263 std::string streamCopy = stream;
264
265 // Massage the final string to book things
266 std::string bookingString("");
267 this->buildBookingString( bookingString, effNameCopy, tDirCopy, streamCopy ,false);
268
269 TEfficiency* effPointer(NULL);
270 if ( !((histSvc()->getEfficiency(bookingString, effPointer)).isSuccess()) )
271 {
272 // Massage the final string to book things
273 std::string bookingString("");
274 this->buildBookingString( bookingString, effNameCopy, tDirCopy, streamCopy, true );
275
276 if ( !((histSvc()->getEfficiency(bookingString, effPointer)).isSuccess()) )
277 {
278 m_msg << MSG::WARNING
279 << "Problem retrieving the TEfficiency with name (including pre- and post-fixes) "
280 << m_histNamePrefix + effNameCopy + m_histNamePostfix
281 << " or with name " << effNameCopy
282 << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
283 << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
284 return NULL;
285 }
286 // If we get to here, we actually found the TEfficiency in the THistSvc.
287 // So let's add it to the local cache map and return its pointer
288 m_effMap.insert( m_effMap.end(), std::pair< const hash_t, TEfficiency* >( effHash, effPointer ) );
289 return effPointer;
290 }
291 // If we get to here, we actually found the TEfficiency in the THistSvc.
292 // So let's add it to the local cache map and return its pointer
293 m_effMap.insert( m_effMap.end(), std::pair< const hash_t, TEfficiency* >( effHash, effPointer ) );
294 return effPointer;
295 }
296
297 // Return the pointer to the TEfficiency that we got from the local cache map
298 return it->second;
299}

◆ evtStore()

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

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

Definition at line 85 of file AthCommonDataStore.h.

◆ execute()

StatusCode ActsTrk::ActsMuonTrackingGeometryTest::execute ( )
override

find the closest propagated hit to the truth hit

Definition at line 107 of file ActsMuonTrackingGeometryTest.cxx.

107 {
108
109 const EventContext& ctx = Gaudi::Hive::currentContext();
110
111 const ActsTrk::GeometryContext* gctx{nullptr};
112 const AtlasFieldCacheCondObj* fieldCondObj{nullptr};
113 const MuonGM::MuonDetectorManager* detMgr{nullptr};
114 const xAOD::TruthParticleContainer* truthParticles{nullptr};
115
116 ATH_CHECK(SG::get(fieldCondObj, m_fieldCacheCondObjInputKey, ctx));
117 ATH_CHECK(SG::get(detMgr, m_detMgrKey, ctx));
118 ATH_CHECK(SG::get(truthParticles, m_truthParticleKey, ctx));
119 ATH_CHECK(SG::get(gctx, m_geoCtxKey, ctx));
120
121
122 const Acts::MagneticFieldContext mfContext = Acts::MagneticFieldContext(fieldCondObj);
123
124 auto anygctx = gctx->context();
125
126 //Get the tracking geometry
127 auto trackingGeometry = m_trackingGeometryTool->trackingGeometry();
128
129 if (!trackingGeometry) {
130 ATH_MSG_ERROR("Failed to retrieve the tracking geometry");
131 return StatusCode::FAILURE;
132 }
133
134 //Configure the ACTS propagator with the navigator and the stepper
135 using Stepper = Acts::EigenStepper<>;
136 using Navigator = Acts::Navigator;
137 using Propagator = Acts::Propagator<Stepper,Navigator>;
138 using ActorList = Acts::ActorList<Acts::detail::SteppingLogger, Acts::MaterialInteractor, Acts::EndOfWorldReached>;
139 using PropagatorOptions = Propagator::Options<ActorList>;
140
141 Navigator::Config navCfg;
142 navCfg.trackingGeometry = trackingGeometry;
143 navCfg.resolveSensitive = true;
144 navCfg.resolveMaterial= false;
145 navCfg.resolvePassive = true;
146
147 Navigator navigator(navCfg, Acts::getDefaultLogger("Navigator", ActsTrk::actsLevelVector(msgLevel())));
148
149 auto bField = std::make_shared<ATLASMagneticFieldWrapper>();
150 auto bfieldCache = bField->makeCache(mfContext);
151
152 auto stepper = Stepper(bField);
153
154 PropagatorOptions options(anygctx, mfContext);
155 options.pathLimit = m_pathLimit;
156 options.stepping.maxStepSize = m_maxStepSize;
157 options.stepping.stepTolerance = m_stepTolerance;
158 options.maxSteps = m_maxSteps;
159 options.maxTargetSkipping = m_maxTargetSkipping;
160
161
162 //switch off material interactions
163 auto& materialInteractor = options.actorList.get<Acts::MaterialInteractor>();
164 materialInteractor.energyLoss = false;
165 materialInteractor.multipleScattering = false;
166 materialInteractor.recordInteractions = false;
167
168 Propagator propagator(std::move(stepper), std::move(navigator),
169 Acts::getDefaultLogger("Propagator", ActsTrk::actsLevelVector(msgLevel())));
170
171 for(const xAOD::TruthParticle* truthParticle : *truthParticles) {
172
173 ATH_MSG_DEBUG("Processing truth particle with PDG ID: " << truthParticle->pdgId() << " ,pT: "
174 << truthParticle->pt() << " , p: " << truthParticle->p4().P() << ", eta: " << truthParticle->eta() << " , phi: " << truthParticle->phi());
175
176 //select only muons or geantinos particles
177 if(std::abs(truthParticle->pdgId()) != 13 && std::abs(truthParticle->pdgId()) != 998) {
178 ATH_MSG_VERBOSE("Skipping truth particle with PDG ID: " << truthParticle->pdgId()<<" only muons or charged geantinos are being processed");
179 continue;
180 }
181
182 Acts::ParticleHypothesis actsParticleHypothesis = truthParticle->pdgId() == 998 ?
183 Acts::ParticleHypothesis::chargedGeantino() : Acts::ParticleHypothesis::muon();
184
185
186 //take the truth segments witht he sim hits
187 std::vector<std::pair<const xAOD::MuonSegment*, std::vector<const xAOD::MuonSimHit*>>> muonSegmentWithSimHits;
188 const SegLink_t& segLink = segAcc(*truthParticle);
189 if (segLink.empty()) {
190 ATH_MSG_WARNING("No segment link found for truth particle with PDG ID: " << truthParticle->pdgId());
191 continue;
192 }
193 for(const auto& truthSegLink: segLink) {
194 const xAOD::MuonSegment* seg{*truthSegLink};
195 if(!seg){
196 continue;
197 }
198 auto unordedHits = MuonR4::getMatchingSimHits(*seg);
199 std::vector<const xAOD::MuonSimHit*> muonSimHits{unordedHits.begin(), unordedHits.end()};
200 ATH_MSG_VERBOSE("Segment at : "<<Amg::toString(seg->position())<<" with Sim Hits: "<<unordedHits.size());
201 muonSegmentWithSimHits.emplace_back(seg,muonSimHits);
202
203 }
204
205 const auto& particle = truthParticle->p4();
206 m_truthPt = particle.Pt();
207 m_truthP = particle.P();
208
209 const xAOD::TruthVertex* prodVertex = truthParticle->prodVtx();
210 Amg::Vector3D prodPos = prodVertex ? Amg::Vector3D(prodVertex->x(), prodVertex->y(), prodVertex->z())
211 : Amg::Vector3D::Zero();
212 Amg::Vector3D prodDir{Amg::dirFromAngles(particle.Phi(), particle.Theta())};
213 ATH_MSG_VERBOSE("Truth particle produced at " << Amg::toString(prodPos) << " with direction " << Amg::toString(prodDir));
214
215 //position, direction and momentum at the start of the propagation
216 Amg::Vector3D startPropPos{prodPos};
217 Amg::Vector3D startPropDir{prodDir};
218 double startPropP{energyToActs(particle.P())};
219
221
222 if(muonSegmentWithSimHits.empty()) {
223 ATH_MSG_DEBUG("No segments found for truth particle with PDG ID: " << truthParticle->pdgId());
224 continue;
225 }
226 // sort the segments by their radial distance
227 std::ranges::sort(muonSegmentWithSimHits,
228 [](const auto& a, const auto& b) {
229 return a.first->position().perp() < b.first->position().perp();
230 });
231
232
233 //sort the sim hits of the first segment to get the starting propagation position as the first one
234 std::vector<const xAOD::MuonSimHit*>& muonSimHits = muonSegmentWithSimHits.front().second;
235 std::ranges::sort(muonSimHits,
236 [this, &gctx](const xAOD::MuonSimHit* hit1, const xAOD::MuonSimHit* hit2){
237 const Amg::Vector3D globalPos1 = toGlobalTrf(*gctx,hit1->identify())*xAOD::toEigen(hit1->localPosition());
238 const Amg::Vector3D globalPos2 = toGlobalTrf(*gctx,hit2->identify())*xAOD::toEigen(hit2->localPosition());
239 return globalPos1.norm()<globalPos2.norm();
240 });
241
242 ATH_MSG_VERBOSE("After sorting..");
243 for(const auto& simHit: muonSimHits){
244 ATH_MSG_VERBOSE("Sim Hit of first segment at position: "<<Amg::toString(toGlobalTrf(*gctx, simHit->identify())*xAOD::toEigen(simHit->localPosition())));
245
246 }
247
248
249 startPropPos = toGlobalTrf(*gctx, muonSimHits.front()->identify())*xAOD::toEigen(muonSimHits.front()->localPosition());
250 startPropDir = toGlobalTrf(*gctx, muonSimHits.front()->identify()).linear()*xAOD::toEigen(muonSimHits.front()->localDirection());
251 startPropP = energyToActs(muonSimHits.front()->kineticEnergy());
252
253
254 ATH_MSG_VERBOSE("Kinetic Energy from the simHit: "<<startPropP / Gaudi::Units::GeV<<" and mass: "<<muonSimHits.front()->mass()<<" and energy deposit: "<<muonSimHits.front()->energyDeposit() / Gaudi::Units::eV);
255
256
257 }
258
259 ATH_MSG_VERBOSE("Starting propagation from"<<Amg::toString(startPropPos)<<" with direction"
260 <<Amg::toString(startPropDir)<<" and momentum "<<startPropP);
261
262
263
264 Acts::BoundTrackParameters start = Acts::BoundTrackParameters::createCurvilinear(
265 Acts::VectorHelpers::makeVector4(startPropPos, 0.), startPropDir,
266 truthParticle->charge() / startPropP,
267 std::nullopt,
268 actsParticleHypothesis);
269
270 ATH_MSG_DEBUG("start propagating here");
271 auto result = propagator.propagate(start, options);
272 const Acts::detail::SteppingLogger::result_type state = result.value().get<Acts::detail::SteppingLogger::result_type>();
273 const Acts::MaterialInteractor::result_type material = result.value().get<Acts::MaterialInteractor::result_type>();
274
275 m_propSteps = state.steps.size();
276 ATH_MSG_DEBUG("Number of propagated steps : " << m_propSteps);
277 m_propLength = result.value().pathLength;
278 std::vector<PropagatorRecorder> propagatedHits;
279
280 for(const auto& step : state.steps) {
281 if(!step.surface) {
282 continue;
283 }
284
285 const SurfaceCache* sCache = dynamic_cast<const SurfaceCache *>(step.surface->associatedDetectorElement());
286 if(!sCache) {
287 ATH_MSG_VERBOSE("Surface found but it's a portal, continuing..");
288 continue;
289 }
290 const Identifier ID = sCache->identify();
291 const Amg::Transform3D toGap{toLocalTrf(*gctx, ID)};
292 ATH_MSG_VERBOSE("Identify propagated hit " << m_idHelperSvc->toString(ID) << " with hit at local position " << Amg::toString(toGap*step.position)<<" and global direction "<<Amg::toString(step.momentum.unit()));
293
294
295 PropagatorRecorder newRecord{};
296 newRecord.id = ID;
297 newRecord.actsPropPos = toGap*step.position;
298 newRecord.actsGlobalPos = step.position;
299 newRecord.actsPropDir = toGap.linear()*step.momentum.unit();
300 newRecord.actsPropabsMomentum = step.momentum.norm();
301 newRecord.actsStepSize = step.stepSize.value();
302 //calculate the distance of the propagated hit to the wire in case of MDTs
303
304 if(m_idHelperSvc->isMdt(ID)){
305 //get the surface
306 const auto* tubeSurf = dynamic_cast<const Acts::StrawSurface*>(&sCache->surface());
307 if(tubeSurf){
308 Amg::Vector3D wirePos = toGap*tubeSurf->center(anygctx);
309 Amg::Vector3D wireDir = toGap.linear()*tubeSurf->lineDirection(anygctx);
310 Amg::Vector3D localPropHit = newRecord.actsPropPos;
311 Amg::Vector3D dirPropHit = newRecord.actsPropDir;
312 double distToWire = Amg::lineDistance(wirePos, wireDir, localPropHit, dirPropHit);
313 newRecord.actsHitWireDist = distToWire;
314 }
315
316 }
317
318 propagatedHits.emplace_back(newRecord);
319
320 }
321
322
323 for (const auto&[segment, simHits] : muonSegmentWithSimHits) {
324 for(const xAOD::MuonSimHit* simHit : simHits){
325
326 const Identifier ID = simHit->identify();
327
328 const Amg::Vector3D localPos = xAOD::toEigen(simHit->localPosition());
329
330 const Amg::Vector3D globalPos = toGlobalTrf(*gctx, simHit->identify())*xAOD::toEigen(simHit->localPosition());
331 const Amg::Vector3D localDir = xAOD::toEigen(simHit->localDirection());
332
333 if(m_r4DetMgr->getReadoutElement(ID)->detectorType() == ActsTrk::DetectorType::sTgc && localPos.z() != 0.0){
334 continue;
335
336 }
337
338 m_detId.push_back(ID);
339 m_techIdx.push_back(toInt(m_idHelperSvc->technologyIndex(ID)));
340 m_gasGapId.push_back(layerHash(ID));
341
342 m_truthLoc.push_back(localPos);
343 m_truthDir.push_back(localDir);
344 m_truthGlob.push_back(globalPos);
345 m_startGlob.push_back(startPropPos);
346 ATH_MSG_DEBUG("Truth hit with ID: " << m_idHelperSvc->toString(ID)
347 << " at local position: " << Amg::toString(localPos)
348 << " and global position: " << Amg::toString(globalPos)
349 << " and direction: " << Amg::toString(localDir));
350
351 //get the matching propagated hits
352 auto it_begin = std::ranges::find_if(propagatedHits,
353 [this, ID](const auto& propagatedHit) {
354 return m_idHelperSvc->detElId(ID) == m_idHelperSvc->detElId(propagatedHit.id) &&
355 layerHash(ID) == layerHash(propagatedHit.id);
356 });
357 m_isPropagated.push_back(it_begin != propagatedHits.end());
358
359 if(it_begin == propagatedHits.end()){
360 m_actsPropLoc.push_back(dummyVec);
361 Amg::Vector3D zero = Amg::Vector3D::Zero();
362 m_actsPropDir.push_back(zero);
363 m_actsPropGlob.push_back(dummyVec);
364 m_actsPropabsMomentum.push_back(0.);
365 m_actsStepSize.push_back(0.);
366 m_actsHitWireDist.push_back(-1.);
367 continue;
368 }
369
370 //find - if any - propagated hits on the same layer
371 auto it_end = std::find_if(it_begin, propagatedHits.end(),
372 [this, ID](const auto& propagatedHit) {
373 return m_idHelperSvc->detElId(ID) != m_idHelperSvc->detElId(propagatedHit.id) ||
374 layerHash(ID) != layerHash(propagatedHit.id);
375 });
376
378 auto it = std::min_element(it_begin, it_end,
379 [&localPos](const PropagatorRecorder& a,
380 const PropagatorRecorder& b){
381 return (localPos - a.actsPropPos).mag() < (localPos - b.actsPropPos).mag();
382 });
383
384 ATH_MSG_DEBUG("Found propagated hit with ID: " << m_idHelperSvc->toString(it->id)
385 << " at local position: " << Amg::toString(it->actsPropPos)
386 << " and global position: " << Amg::toString(it->actsGlobalPos)
387 << " and direction: " << Amg::toString(it->actsPropDir));
388
389 m_actsPropLoc.push_back(it->actsPropPos);
390 m_actsPropDir.push_back(it->actsPropDir);
391 m_actsPropGlob.push_back(it->actsGlobalPos);
392 m_actsPropabsMomentum.push_back(it->actsPropabsMomentum);
393 m_actsStepSize.push_back(it->actsStepSize);
394 m_actsHitWireDist.push_back(it->actsHitWireDist);
395
396
397 }
398 }
399
400 m_event = ctx.eventID().event_number();
401 m_tree.fill(ctx);
402
403
404 }
405
406
407 return StatusCode::SUCCESS;
408 }
Scalar mag() const
mag method
#define ATH_CHECK
Evaluate an expression and check for errors.
#define ATH_MSG_ERROR(x)
#define ATH_MSG_VERBOSE(x)
#define ATH_MSG_WARNING(x)
#define ATH_MSG_DEBUG(x)
std::vector< Identifier > ID
static Double_t a
Eigen::Matrix< double, 3, 1 > Vector3D
MuonVal::ScalarBranch< unsigned int > & m_propSteps
MuonVal::VectorBranch< unsigned short > & m_isPropagated
SG::ReadHandleKey< ActsTrk::GeometryContext > m_geoCtxKey
MuonVal::VectorBranch< unsigned short > & m_gasGapId
MuonVal::VectorBranch< float > & m_actsHitWireDist
SG::ReadHandleKey< xAOD::TruthParticleContainer > m_truthParticleKey
const MuonGMR4::MuonDetectorManager * m_r4DetMgr
ServiceHandle< Muon::IMuonIdHelperSvc > m_idHelperSvc
Amg::Transform3D toGlobalTrf(const ActsTrk::GeometryContext &gctx, const Identifier &hitId) const
MuonVal::VectorBranch< float > & m_actsPropabsMomentum
SG::ReadCondHandleKey< AtlasFieldCacheCondObj > m_fieldCacheCondObjInputKey
IdentifierHash layerHash(const Identifier &id) const
MuonVal::ScalarBranch< unsigned int > & m_event
Amg::Transform3D toLocalTrf(const ActsTrk::GeometryContext &gctx, const Identifier &hitId) const
MuonVal::VectorBranch< unsigned short > & m_techIdx
SG::ReadCondHandleKey< MuonGM::MuonDetectorManager > m_detMgrKey
MuonVal::VectorBranch< float > & m_actsStepSize
PublicToolHandle< ActsTrk::ITrackingGeometryTool > m_trackingGeometryTool
Acts::GeometryContext context() const
Amg::Vector3D position() const
Returns the position as Amg::Vector.
Identifier identify() const
Returns the global ATLAS identifier of the SimHit.
ConstVectorMap< 3 > localPosition() const
Returns the local postion of the traversing particle.
float z() const
Vertex longitudinal distance along the beam line form the origin.
float y() const
Vertex y displacement.
float x() const
Vertex x displacement.
void zero(TH2 *h)
zero the contents of a 2d histogram
Acts::SympyStepper Stepper
Adapted from Acts Examples/Algorithms/TrackFinding/src/TrackFindingAlgorithmFunction....
constexpr double energyToActs(const double athenaE)
Converts an energy scalar from Athena to Acts units.
Acts::Logging::Level actsLevelVector(MSG::Level lvl)
@ sTgc
Micromegas (NSW)
std::string toString(const Translation3D &translation, int precision=4)
GeoPrimitvesToStringConverter.
double lineDistance(const AmgVector(N)&posA, const AmgVector(N)&dirA, const AmgVector(N)&posB, const AmgVector(N)&dirB)
: Calculates the shortest distance between two lines
Amg::Vector3D dirFromAngles(const double phi, const double theta)
Constructs a direction vector from the azimuthal & polar angles.
Eigen::Affine3d Transform3D
Eigen::Matrix< double, 3, 1 > Vector3D
std::unordered_set< const xAOD::MuonSimHit * > getMatchingSimHits(const xAOD::MuonSegment &segment)
: Returns all sim hits matched to a xAOD::MuonSegment
ElementLink< MuonR4::SegmentContainer > SegLink_t
Abrivation of the link to the reco segment container.
constexpr int toInt(const EnumType enumVal)
const T * get(const ReadCondHandleKey< T > &key, const EventContext &ctx)
Convenience function to retrieve an object given a ReadCondHandleKey.
constexpr ParticleHypothesis particle[PARTICLEHYPOTHESES]
the array of masses
void Zero(TH1D *hin)
Definition generate.cxx:32
TruthVertex_v1 TruthVertex
Typedef to implementation.
Definition TruthVertex.h:15
MuonSimHit_v1 MuonSimHit
Defined the version of the MuonSimHit.
Definition MuonSimHit.h:12
TruthParticle_v1 TruthParticle
Typedef to implementation.
MuonSegment_v1 MuonSegment
Reference the current persistent version:
TruthParticleContainer_v1 TruthParticleContainer
Declare the latest version of the truth particle container.

◆ extraDeps_update_handler()

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

Return the list of extra output dependencies.

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

Definition at line 50 of file AthAlgorithm.cxx.

51{
52 // If we didn't find any symlinks to add, just return the collection
53 // from the base class. Otherwise, return the extended collection.
54 if (!m_extendedExtraObjects.empty()) {
56 }
57 return Algorithm::extraOutputDeps();
58}
DataObjIDColl m_extendedExtraObjects

◆ finalize()

StatusCode ActsTrk::ActsMuonTrackingGeometryTest::finalize ( )
override

Definition at line 102 of file ActsMuonTrackingGeometryTest.cxx.

102 {
103 ATH_CHECK(m_tree.write());
104 return StatusCode::SUCCESS;
105 }

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

493{
494 // Build a 32 bit hash out of the name
495 const hash_t graphHash = this->hash(graphName);
496
497 // See if this entry exists in the map
498 GraphMap_t::const_iterator it = m_graphMap.find( graphHash );
499 if ( it == m_graphMap.end() ) // It doesn't exist!
500 { // Let's see into the THistSvc if somebody else has registered the TGraph...
501
502 // Need to copy the strings as we will massage them from here on
503 std::string graphNameCopy = graphName;
504 std::string tDirCopy = tDir;
505 std::string streamCopy = stream;
506
507 // Massage the final string to book things
508 std::string bookingString("");
509 this->buildBookingString( bookingString, graphNameCopy, tDirCopy, streamCopy, true);
510
511 TGraph* graphPointer(NULL);
512 if ( !((histSvc()->getGraph(bookingString, graphPointer)).isSuccess()) )
513 {
514 // Massage the final string to book things
515 std::string bookingString("");
516 this->buildBookingString( bookingString, graphNameCopy, tDirCopy, streamCopy, false );
517
518 if ( !((histSvc()->getGraph(bookingString, graphPointer)).isSuccess()) )
519 {
520 m_msg << MSG::WARNING
521 << "Problem retrieving the TGraph with name (including pre- and post-fixes) "
522 << m_histNamePrefix + graphNameCopy + m_histNamePostfix
523 << " or with name " << graphNameCopy
524 << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
525 << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
526 return NULL;
527 }
528 // If we get to here, we actually found the TGraph in the THistSvc.
529 // So let's add it to the local cache map and return its pointer
530 m_graphMap.insert( m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphPointer ) );
531 return graphPointer;
532 }
533 // If we get to here, we actually found the TGraph in the THistSvc.
534 // So let's add it to the local cache map and return its pointer
535 m_graphMap.insert( m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphPointer ) );
536 return graphPointer;
537 }
538
539
540 // Return the pointer to the TGraph that we got from the local cache map
541 return it->second;
542}

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

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

◆ hist()

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

Simplify the retrieval of registered histograms of any type.

Definition at line 198 of file AthHistogramming.cxx.

199{
200 // Build a 32 bit hash out of the name
201 const hash_t histHash = this->hash(histName);
202
203 // See if this entry exists in the map
204 HistMap_t::const_iterator it = m_histMap.find( histHash );
205 if ( it == m_histMap.end() ) // It doesn't exist!
206 { // Let's see into the THistSvc if somebody else has registered the histogram...
207
208 // Need to copy the strings as we will massage them from here on
209 std::string histNameCopy = histName;
210 std::string tDirCopy = tDir;
211 std::string streamCopy = stream;
212
213 // Massage the final string to book things
214 std::string bookingString("");
215 this->buildBookingString( bookingString, histNameCopy, tDirCopy, streamCopy ,false);
216
217 TH1* histPointer(NULL);
218 if ( !((histSvc()->getHist(bookingString, histPointer)).isSuccess()) )
219 {
220 // Massage the final string to book things
221 std::string bookingString("");
222 this->buildBookingString( bookingString, histNameCopy, tDirCopy, streamCopy, true );
223
224 if ( !((histSvc()->getHist(bookingString, histPointer)).isSuccess()) )
225 {
226 m_msg << MSG::WARNING
227 << "Problem retrieving the histogram with name (including pre- and post-fixes) "
228 << m_histNamePrefix + histNameCopy + m_histNamePostfix
229 << " or with name " << histNameCopy
230 << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
231 << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
232 return NULL;
233 }
234 // If we get to here, we actually found the histogram in the THistSvc.
235 // So let's add it to the local cache map and return its pointer
236 m_histMap.insert( m_histMap.end(), std::pair< const hash_t, TH1* >( histHash, histPointer ) );
237 return histPointer;
238 }
239 // If we get to here, we actually found the histogram in the THistSvc.
240 // So let's add it to the local cache map and return its pointer
241 m_histMap.insert( m_histMap.end(), std::pair< const hash_t, TH1* >( histHash, histPointer ) );
242 return histPointer;
243 }
244
245
246 // Return the pointer to the histogram that we got from the local cache map
247 return it->second;
248}

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

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

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

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

◆ histSvc()

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

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

Definition at line 113 of file AthHistogramAlgorithm.h.

114{
115 return m_histSvc;
116}

◆ initialize()

StatusCode ActsTrk::ActsMuonTrackingGeometryTest::initialize ( )
override

Definition at line 85 of file ActsMuonTrackingGeometryTest.cxx.

85 {
86 ATH_CHECK(AthHistogramAlgorithm::initialize());
87 ATH_CHECK(m_idHelperSvc.retrieve());
89 ATH_CHECK(m_detMgrKey.initialize());
90 ATH_CHECK(m_rndmGenSvc.retrieve());
91 ATH_CHECK(m_geoCtxKey.initialize());
93 ATH_CHECK(m_truthParticleKey.initialize());
94 ATH_CHECK(m_truthSegLinkKey.initialize());
96
97 ATH_CHECK(m_tree.init(this));
98 ATH_MSG_INFO("ActsMuonTrackingGeometryTest successfully initialized");
99 return StatusCode::SUCCESS;
100 }
#define ATH_MSG_INFO(x)
SG::ReadDecorHandleKey< xAOD::TruthParticleContainer > m_truthSegLinkKey
const ServiceHandle< StoreGateSvc > & detStore() const
retrieve(aClass, aKey=None)
Definition PyKernel.py:110

◆ inputHandles()

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

◆ layerHash()

IdentifierHash ActsTrk::ActsMuonTrackingGeometryTest::layerHash ( const Identifier & id) const
private

Definition at line 80 of file ActsMuonTrackingGeometryTest.cxx.

80 {
81 return m_r4DetMgr->getReadoutElement(id)->layerHash(id);
82 }

◆ msg()

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

Definition at line 24 of file AthCommonMsg.h.

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

◆ msgLvl()

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

Definition at line 30 of file AthCommonMsg.h.

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

◆ myReplace()

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

Helper method to replace sub-string.

Definition at line 590 of file AthHistogramming.cxx.

593{
594 size_t pos = 0;
595 while((pos = str.find(oldStr, pos)) != std::string::npos)
596 {
597 str.replace(pos, oldStr.length(), newStr);
598 pos += newStr.length();
599 }
600}

◆ outputHandles()

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

Definition at line 380 of file AthCommonDataStore.h.

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

◆ renounceArray()

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

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364 {
366 }

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

Reimplemented in AthAnalysisAlgorithm.

Definition at line 75 of file AthHistogramAlgorithm.cxx.

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

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< 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.

◆ toGlobalTrf()

Amg::Transform3D ActsTrk::ActsMuonTrackingGeometryTest::toGlobalTrf ( const ActsTrk::GeometryContext & gctx,
const Identifier & hitId ) const
private

Definition at line 72 of file ActsMuonTrackingGeometryTest.cxx.

72 {
73 const MuonGMR4::MuonReadoutElement* reElement = m_r4DetMgr->getReadoutElement(hitId);
74 const IdentifierHash trfHash = reElement->detectorType() == ActsTrk::DetectorType::Mdt ?
75 reElement->measurementHash(hitId) : reElement->layerHash(hitId);
76 return reElement->localToGlobalTrans(gctx, trfHash);
77 }
virtual DetectorType detectorType() const =0
Returns the detector element type.
const Amg::Transform3D & localToGlobalTrans(const ActsTrk::GeometryContext &ctx) const
Returns the local to global transformation into the ATLAS coordinate system.
virtual IdentifierHash measurementHash(const Identifier &measId) const =0
Constructs the identifier hash from the full measurement Identifier.
virtual IdentifierHash layerHash(const Identifier &measId) const =0
@ Mdt
MuonSpectrometer.

◆ toLocalTrf()

Amg::Transform3D ActsTrk::ActsMuonTrackingGeometryTest::toLocalTrf ( const ActsTrk::GeometryContext & gctx,
const Identifier & hitId ) const
private

Definition at line 65 of file ActsMuonTrackingGeometryTest.cxx.

65 {
66 const MuonGMR4::MuonReadoutElement* reElement = m_r4DetMgr->getReadoutElement(hitId);
67 const IdentifierHash trfHash = reElement->detectorType() == ActsTrk::DetectorType::Mdt ?
68 reElement->measurementHash(hitId) : reElement->layerHash(hitId);
69 return reElement->globalToLocalTrans(gctx, trfHash);
70 }
Amg::Transform3D globalToLocalTrans(const ActsTrk::GeometryContext &ctx) const
Transformations to translate between local <-> global coordinates.

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

379{
380 // Build a 32 bit hash out of the name
381 const hash_t treeHash = this->hash(treeName);
382
383 // See if this entry exists in the map
384 TreeMap_t::const_iterator it = m_treeMap.find( treeHash );
385 if ( it == m_treeMap.end() ) // It doesn't exist!
386 { // Let's see into the THistSvc if somebody else has registered the TTree...
387
388 // Need to copy the strings as we will massage them from here on
389 std::string treeNameCopy = treeName;
390 std::string tDirCopy = tDir;
391 std::string streamCopy = stream;
392
393 // Massage the final string to book things
394 std::string bookingString("");
395 this->buildBookingString( bookingString, treeNameCopy, tDirCopy, streamCopy );
396
397 TTree* treePointer(NULL);
398 if ( !((histSvc()->getTree(bookingString, treePointer)).isSuccess()) )
399 {
400 m_msg << MSG::WARNING
401 << "Problem retrieving the TTree with name " << treeNameCopy
402 << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
403 << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
404 return NULL;
405 }
406 // If we get to here, we actually found the TTree in the THistSvc.
407 // So let's add it to the local cache map and return its pointer
408 m_treeMap.insert( m_treeMap.end(), std::pair< const hash_t, TTree* >( treeHash, treePointer ) );
409 return treePointer;
410 }
411
412 // Return the pointer to the TTree that we got from the local cache map
413 return it->second;
414}

◆ updateVHKA()

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

Definition at line 308 of file AthCommonDataStore.h.

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

Member Data Documentation

◆ m_actsHitWireDist

MuonVal::VectorBranch<float>& ActsTrk::ActsMuonTrackingGeometryTest::m_actsHitWireDist {m_tree.newVector<float>("actsHitWireDist")}
private

Definition at line 106 of file ActsMuonTrackingGeometryTest.h.

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

◆ m_actsPropabsMomentum

MuonVal::VectorBranch<float>& ActsTrk::ActsMuonTrackingGeometryTest::m_actsPropabsMomentum {m_tree.newVector<float>("actsPropabsMomentum")}
private

Definition at line 102 of file ActsMuonTrackingGeometryTest.h.

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

◆ m_actsPropDir

MuonVal::ThreeVectorBranch ActsTrk::ActsMuonTrackingGeometryTest::m_actsPropDir {m_tree, "actsPropDir"}
private

Definition at line 105 of file ActsMuonTrackingGeometryTest.h.

105{m_tree, "actsPropDir"};

◆ m_actsPropGlob

MuonVal::ThreeVectorBranch ActsTrk::ActsMuonTrackingGeometryTest::m_actsPropGlob {m_tree, "actsPropGlob"}
private

Definition at line 104 of file ActsMuonTrackingGeometryTest.h.

104{m_tree, "actsPropGlob"};

◆ m_actsPropLoc

MuonVal::ThreeVectorBranch ActsTrk::ActsMuonTrackingGeometryTest::m_actsPropLoc {m_tree, "actsPropLoc"}
private

Definition at line 103 of file ActsMuonTrackingGeometryTest.h.

103{m_tree, "actsPropLoc"};

◆ m_actsStepSize

MuonVal::VectorBranch<float>& ActsTrk::ActsMuonTrackingGeometryTest::m_actsStepSize {m_tree.newVector<float>("actsStepSize")}
private

Definition at line 107 of file ActsMuonTrackingGeometryTest.h.

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

◆ m_detId

MuonVal::MuonIdentifierBranch ActsTrk::ActsMuonTrackingGeometryTest::m_detId {m_tree, "detId"}
private

Definition at line 93 of file ActsMuonTrackingGeometryTest.h.

93{m_tree, "detId"};

◆ m_detMgrKey

SG::ReadCondHandleKey<MuonGM::MuonDetectorManager> ActsTrk::ActsMuonTrackingGeometryTest::m_detMgrKey
private
Initial value:
{this, "MuonManagerKey",
"MuonDetectorManager", "MuonManager ReadKey for IOV Range intersection"}

Definition at line 58 of file ActsMuonTrackingGeometryTest.h.

58 {this, "MuonManagerKey",
59 "MuonDetectorManager", "MuonManager ReadKey for IOV Range intersection"};

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< 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 208 of file AthHistogramming.h.

◆ m_etaRange

Gaudi::Property<std::vector<double> > ActsTrk::ActsMuonTrackingGeometryTest::m_etaRange {this, "EtaRange", {-3, 3}, "The eta range for particles"}
private

Definition at line 90 of file ActsMuonTrackingGeometryTest.h.

90{this, "EtaRange", {-3, 3}, "The eta range for particles"};

◆ m_event

MuonVal::ScalarBranch<unsigned int>& ActsTrk::ActsMuonTrackingGeometryTest::m_event {m_tree.newScalar<unsigned int>("event")}
private

Definition at line 111 of file ActsMuonTrackingGeometryTest.h.

111{m_tree.newScalar<unsigned int>("event")};

◆ m_evtStore

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

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_extendedExtraObjects

DataObjIDColl AthAlgorithm::m_extendedExtraObjects
privateinherited

Definition at line 79 of file AthAlgorithm.h.

◆ m_fieldCacheCondObjInputKey

SG::ReadCondHandleKey<AtlasFieldCacheCondObj> ActsTrk::ActsMuonTrackingGeometryTest::m_fieldCacheCondObjInputKey
private
Initial value:
{this, "AtlasFieldCacheCondObj",
"fieldCondObj", "Name of the Magnetic Field conditions object key"}

Definition at line 61 of file ActsMuonTrackingGeometryTest.h.

61 {this, "AtlasFieldCacheCondObj",
62 "fieldCondObj", "Name of the Magnetic Field conditions object key"};

◆ m_gasGapId

MuonVal::VectorBranch<unsigned short>& ActsTrk::ActsMuonTrackingGeometryTest::m_gasGapId {m_tree.newVector<unsigned short>("detId_gasGap")}
private

Definition at line 95 of file ActsMuonTrackingGeometryTest.h.

95{m_tree.newVector<unsigned short>("detId_gasGap")};

◆ m_geoCtxKey

SG::ReadHandleKey<ActsTrk::GeometryContext> ActsTrk::ActsMuonTrackingGeometryTest::m_geoCtxKey {this, "AlignmentKey", "ActsAlignment", "cond handle key"}
private

Definition at line 64 of file ActsMuonTrackingGeometryTest.h.

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

◆ m_graphMap

GraphMap_t AthHistogramming::m_graphMap
privateinherited

The map of TGraph names to their pointers.

Definition at line 222 of file AthHistogramming.h.

◆ m_histMap

HistMap_t AthHistogramming::m_histMap
privateinherited

The map of histogram names to their pointers.

Definition at line 201 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_idHelperSvc

ServiceHandle<Muon::IMuonIdHelperSvc> ActsTrk::ActsMuonTrackingGeometryTest::m_idHelperSvc {this, "IdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}
private

Definition at line 52 of file ActsMuonTrackingGeometryTest.h.

52{this, "IdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"};

◆ m_isPropagated

MuonVal::VectorBranch<unsigned short>& ActsTrk::ActsMuonTrackingGeometryTest::m_isPropagated {m_tree.newVector<unsigned short>("isPropagated")}
private

Definition at line 109 of file ActsMuonTrackingGeometryTest.h.

109{m_tree.newVector<unsigned short>("isPropagated")};

◆ m_maxSteps

UnsignedIntegerProperty ActsTrk::ActsMuonTrackingGeometryTest::m_maxSteps {this, "MaxSteps", 100000}
private

Definition at line 76 of file ActsMuonTrackingGeometryTest.h.

76{this, "MaxSteps", 100000};

◆ m_maxStepSize

DoubleProperty ActsTrk::ActsMuonTrackingGeometryTest::m_maxStepSize {this, "MaxStepSize", 1*Gaudi::Units::m}
private

Definition at line 74 of file ActsMuonTrackingGeometryTest.h.

74{this, "MaxStepSize", 1*Gaudi::Units::m};

◆ m_maxTargetSkipping

UnsignedIntegerProperty ActsTrk::ActsMuonTrackingGeometryTest::m_maxTargetSkipping {this, "MaxTargetSkipping", 10000}
private

Definition at line 78 of file ActsMuonTrackingGeometryTest.h.

78{this, "MaxTargetSkipping", 10000};

◆ m_msg

MsgStream AthHistogramming::m_msg
privateinherited

Cached Message Stream.

Definition at line 248 of file AthHistogramming.h.

◆ m_name

std::string AthHistogramming::m_name
privateinherited

Instance name.

Definition at line 245 of file AthHistogramming.h.

◆ m_nEvents

UnsignedIntegerProperty ActsTrk::ActsMuonTrackingGeometryTest::m_nEvents {this, "nEvents", 100}
private

Definition at line 70 of file ActsMuonTrackingGeometryTest.h.

70{this, "nEvents", 100};

◆ m_pathLimit

DoubleProperty ActsTrk::ActsMuonTrackingGeometryTest::m_pathLimit {this, "PathLimit", 30.*Gaudi::Units::m}
private

Definition at line 72 of file ActsMuonTrackingGeometryTest.h.

72{this, "PathLimit", 30.*Gaudi::Units::m};

◆ m_prefix

std::string AthHistogramAlgorithm::m_prefix
privateinherited

Name of the ROOT output stream (file)

Definition at line 88 of file AthHistogramAlgorithm.h.

◆ m_propLength

MuonVal::ScalarBranch<float>& ActsTrk::ActsMuonTrackingGeometryTest::m_propLength {m_tree.newScalar<float>("propLength")}
private

Definition at line 108 of file ActsMuonTrackingGeometryTest.h.

108{m_tree.newScalar<float>("propLength")};

◆ m_propSteps

MuonVal::ScalarBranch<unsigned int>& ActsTrk::ActsMuonTrackingGeometryTest::m_propSteps {m_tree.newScalar<unsigned int>("propSteps")}
private

Definition at line 110 of file ActsMuonTrackingGeometryTest.h.

110{m_tree.newScalar<unsigned int>("propSteps")};

◆ m_r4DetMgr

const MuonGMR4::MuonDetectorManager* ActsTrk::ActsMuonTrackingGeometryTest::m_r4DetMgr {nullptr}
private

Definition at line 50 of file ActsMuonTrackingGeometryTest.h.

50{nullptr};

◆ m_rndmGenSvc

ServiceHandle<IAthRNGSvc> ActsTrk::ActsMuonTrackingGeometryTest::m_rndmGenSvc {this, "AthRNGSvc", "AthRNGSvc"}
private

Definition at line 54 of file ActsMuonTrackingGeometryTest.h.

54{this, "AthRNGSvc", "AthRNGSvc"};

◆ m_rootDir

std::string AthHistogramAlgorithm::m_rootDir
privateinherited

Name of the ROOT directory.

Definition at line 91 of file AthHistogramAlgorithm.h.

◆ m_startFromFirstHit

Gaudi::Property<bool> ActsTrk::ActsMuonTrackingGeometryTest::m_startFromFirstHit {this, "StartFromFirstHit", false, "Start from first hit"}
private

Definition at line 82 of file ActsMuonTrackingGeometryTest.h.

82{this, "StartFromFirstHit", false, "Start from first hit"};

◆ m_startGlob

MuonVal::ThreeVectorBranch ActsTrk::ActsMuonTrackingGeometryTest::m_startGlob {m_tree, "startGlob"}
private

Definition at line 101 of file ActsMuonTrackingGeometryTest.h.

101{m_tree, "startGlob"};

◆ m_stepTolerance

DoubleProperty ActsTrk::ActsMuonTrackingGeometryTest::m_stepTolerance {this, "StepTolerance", 1e-10}
private

Definition at line 80 of file ActsMuonTrackingGeometryTest.h.

80{this, "StepTolerance", 1e-10};

◆ m_streamName

std::string AthHistogramming::m_streamName
privateinherited

Name of the ROOT output stream (file)

Definition at line 226 of file AthHistogramming.h.

◆ m_techIdx

MuonVal::VectorBranch<unsigned short>& ActsTrk::ActsMuonTrackingGeometryTest::m_techIdx {m_tree.newVector<unsigned short>("detId_techIdx")}
private

Definition at line 94 of file ActsMuonTrackingGeometryTest.h.

94{m_tree.newVector<unsigned short>("detId_techIdx")};

◆ m_trackingGeometryTool

PublicToolHandle<ActsTrk::ITrackingGeometryTool> ActsTrk::ActsMuonTrackingGeometryTest::m_trackingGeometryTool {this, "TrackingGeometryTool", "ActsTrackingGeometryTool"}
private

Definition at line 56 of file ActsMuonTrackingGeometryTest.h.

56{this, "TrackingGeometryTool", "ActsTrackingGeometryTool"};

◆ m_tree

MuonVal::MuonTesterTree ActsTrk::ActsMuonTrackingGeometryTest::m_tree {"MuonNavigationTestGen3R4", "MuonNavigationTestGen3R4"}
private

Definition at line 92 of file ActsMuonTrackingGeometryTest.h.

92{"MuonNavigationTestGen3R4", "MuonNavigationTestGen3R4"};

◆ m_treeMap

TreeMap_t AthHistogramming::m_treeMap
privateinherited

The map of TTree names to their pointers.

Definition at line 215 of file AthHistogramming.h.

◆ m_truthDir

MuonVal::ThreeVectorBranch ActsTrk::ActsMuonTrackingGeometryTest::m_truthDir {m_tree, "truthDir"}
private

Definition at line 98 of file ActsMuonTrackingGeometryTest.h.

98{m_tree, "truthDir"};

◆ m_truthGlob

MuonVal::ThreeVectorBranch ActsTrk::ActsMuonTrackingGeometryTest::m_truthGlob {m_tree, "truthGlob"}
private

Definition at line 97 of file ActsMuonTrackingGeometryTest.h.

97{m_tree, "truthGlob"};

◆ m_truthLoc

MuonVal::ThreeVectorBranch ActsTrk::ActsMuonTrackingGeometryTest::m_truthLoc {m_tree, "truthLoc"}
private

Definition at line 96 of file ActsMuonTrackingGeometryTest.h.

96{m_tree, "truthLoc"};

◆ m_truthP

MuonVal::ScalarBranch<float>& ActsTrk::ActsMuonTrackingGeometryTest::m_truthP {m_tree.newScalar<float>("truthP")}
private

Definition at line 100 of file ActsMuonTrackingGeometryTest.h.

100{m_tree.newScalar<float>("truthP")};

◆ m_truthParticleKey

SG::ReadHandleKey<xAOD::TruthParticleContainer> ActsTrk::ActsMuonTrackingGeometryTest::m_truthParticleKey {this, "TruthKey", "MuonTruthParticles"}
private

Definition at line 66 of file ActsMuonTrackingGeometryTest.h.

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

◆ m_truthPt

MuonVal::ScalarBranch<float>& ActsTrk::ActsMuonTrackingGeometryTest::m_truthPt {m_tree.newScalar<float>("truthPt")}
private

Definition at line 99 of file ActsMuonTrackingGeometryTest.h.

99{m_tree.newScalar<float>("truthPt")};

◆ m_truthSegLinkKey

SG::ReadDecorHandleKey<xAOD::TruthParticleContainer> ActsTrk::ActsMuonTrackingGeometryTest::m_truthSegLinkKey {this, "TruthKeyToSeg", m_truthParticleKey, "truthSegmentLinks", "TruthParticle to TruthSegment link"}
private

Definition at line 68 of file ActsMuonTrackingGeometryTest.h.

68{this, "TruthKeyToSeg", m_truthParticleKey, "truthSegmentLinks", "TruthParticle to TruthSegment link"};

◆ m_varHandleArraysDeclared

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

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

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

Definition at line 398 of file AthCommonDataStore.h.


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