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

#include <MuonHoughTransformTester.h>

Inheritance diagram for MuonValR4::MuonHoughTransformTester:

Classes

struct  ObjectMatching

Public Types

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

Public Member Functions

virtual ~MuonHoughTransformTester ()=default
virtual StatusCode initialize () override
virtual StatusCode execute () override
virtual 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

std::vector< ObjectMatchingmatchWithTruth (const ActsTrk::GeometryContext &gctx, const MuonR4::SegmentSeedContainer &seedContainer, const xAOD::MuonSegmentContainer &segmentContainer, const xAOD::MuonSegmentContainer *truthSegments) const
unsigned int countOnSameSide (const ActsTrk::GeometryContext &gctx, const xAOD::MuonSegment &truthSeg, const MuonR4::Segment &recoSeg) const
 Calculates how many measurements from the segment fit have the same drift sign as when evaluated with the truth parameters.
void fillChamberInfo (const MuonGMR4::SpectrometerSector *chamber)
 Fill the current chamber info into the output.
void fillTruthInfo (const ActsTrk::GeometryContext &gctx, const xAOD::MuonSegment *truthSegment)
 Fill the associated truth information into the tree.
void fillBucketInfo (const MuonR4::SpacePointBucket &bucket)
 Fill the hit summary info of the associated bucket.
void fillSeedInfo (const ObjectMatching &obj)
 Fill the info associated to the seed.
void fillSegmentInfo (const ActsTrk::GeometryContext &gctx, const ObjectMatching &obj)
 Fill the info assciated to the segment.
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 {"MuonEtaHoughTest","MuonEtaHoughTransformTest"}
SG::ReadHandleKey< xAOD::MuonSegmentContainerm_truthSegmentKey {this, "TruthSegmentKey","MuonTruthSegments", "truth segment container"}
 Key to the truth segment.
SG::ReadDecorHandleKeyArray< xAOD::MuonSegmentContainerm_truthSegLinkKeys {this, "TruthSegLinkKeys", {}}
 Declare the dependencies on the decorations.
Gaudi::Property< std::vector< std::string > > m_truthLinks {this, "TruthSegLinks", {"simHitLinks", "truthParticleLink"}}
 Name of the decorations for the truth segment.
SG::ReadHandleKey< xAOD::MuonSegmentContainerm_recoSegKey {this, "SegmentKey", "MuonSegmentsFromR4"}
 Key to the xAOD::MuonSegment container.
Gaudi::Property< std::vector< std::string > > m_recoSegLinks {this, "RecoSegLinks", {"truthSegmentLink", "truthParticleLink"}}
 name of the truth link decorations for the reco segment container
SG::ReadHandleKeyArray< MuonR4::SegmentSeedContainerm_patternSeedKeys {this, "SegmentSeedKeys", {"MuonHoughStationSegmentSeeds"}}
 List of the two segment seed containers from which the segments are buiit (Complets the pattern finding step).
SG::ReadHandleKeyArray< MuonR4::SpacePointContainerm_spKeys {this, "SpacePointKeys", {"MuonSpacePoints"}}
 List of the space point containers in the event legacy + NSW containers.
SG::ReadHandleKey< ActsTrk::GeometryContextm_geoCtxKey {this, "AlignmentKey", "ActsAlignment", "cond handle key"}
 Tracking geometry context.
ServiceHandle< Muon::IMuonIdHelperSvcm_idHelperSvc {this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}
Gaudi::Property< bool > m_isMC {this, "isMC", false, "Toggle whether the job is ran on MC or not"}
Gaudi::Property< bool > m_writeSpacePoints
MuonVal::ScalarBranch< int > & m_out_chamberIndex {m_tree.newScalar<int>("chamberIndex")}
 ====== Common block: Filled for all entries ===========
MuonVal::ScalarBranch< short > & m_out_stationSide {m_tree.newScalar<short>("stationSide")}
 +1 for A-, -1 of C-side
MuonVal::ScalarBranch< int > & m_out_stationPhi {m_tree.newScalar<int>("stationPhi")}
 phi index of the station
MuonVal::ScalarBranch< float > & m_out_bucketStart {m_tree.newScalar<float>("bucketStart", 1)}
MuonVal::ScalarBranch< float > & m_out_bucketEnd {m_tree.newScalar<float>("bucketEnd", -1)}
MuonVal::ScalarBranch< float > & m_out_bucketHitGap {m_tree.newScalar<float>("bucketHiGap", 0.)}
std::shared_ptr< SpacePointTesterModulem_spTester {}
 Branch dumping all the space points from the difference buckets.
MuonVal::MatrixBranch< unsigned char > & m_spMatchedToPattern {m_tree.newMatrix<unsigned char>("seedMatchedSp")}
 Branch indicating which space points in the tree are associated to the i-th pattern.
MuonVal::MatrixBranch< unsigned char > & m_spMatchedToSegment {m_tree.newMatrix<unsigned char>("segmentMatchedSp")}
 Branch indicating which space points in the tree are associated to the i-th segment.
MuonVal::ScalarBranch< unsigned char > & m_out_nSpacePoints {m_tree.newScalar<unsigned char>("bucketNHits",0)}
 Number of all space points in the bucket.
MuonVal::ScalarBranch< unsigned char > & m_out_nPrecSpacePoints {m_tree.newScalar<unsigned char>("bucketNPrecMeas",0)}
 Number of precision hits in the bucket.
MuonVal::ScalarBranch< unsigned char > & m_out_nPhiSpacePoints {m_tree.newScalar<unsigned char>("bucketNPhiMeass",0)}
 Number of phi hits in the bucket.
MuonVal::ScalarBranch< unsigned char > & m_out_nTrueSpacePoints {m_tree.newScalar<unsigned char>("bucketNTrueMeas",0)}
 Number of all space points in the bucket.
MuonVal::ScalarBranch< unsigned char > & m_out_nTruePrecSpacePoints {m_tree.newScalar<unsigned char>("bucketNTruePrecMeas",0)}
 Number of precision hits in the bucket.
MuonVal::ScalarBranch< unsigned char > & m_out_nTruePhiSpacePoints {m_tree.newScalar<unsigned char>("bucketNTruePhiMeass",0)}
 Number of phi hits in the bucket.
MuonVal::ScalarBranch< bool > & m_out_hasTruth {m_tree.newScalar<bool>("hasTruth",false)}
 ======= Truth block: Filled if we have a truth match. ============
MuonVal::ScalarBranch< float > & m_out_gen_Eta {m_tree.newScalar<float>("genEta",-10.)}
 global particle properties
MuonVal::ScalarBranch< float > & m_out_gen_Phi {m_tree.newScalar<float>("genPhi",-10.)}
MuonVal::ScalarBranch< float > & m_out_gen_Pt {m_tree.newScalar<float>("genPt",-10.)}
MuonVal::ScalarBranch< short > & m_out_gen_Q {m_tree.newScalar<short>("genQ", 0)}
MuonVal::ScalarBranch< float > & m_out_gen_y0 {m_tree.newScalar<float>("genY0", 0.0)}
 Truth - segment parameters.
MuonVal::ScalarBranch< float > & m_out_gen_tantheta {m_tree.newScalar<float>("genTanTheta", 0.0)}
MuonVal::ScalarBranch< float > & m_out_gen_tanphi {m_tree.newScalar<float>("genTanPhi", 0.0)}
MuonVal::ScalarBranch< float > & m_out_gen_x0 {m_tree.newScalar<float>("genX0", 0.0)}
MuonVal::ScalarBranch< float > & m_out_gen_time {m_tree.newScalar<float>("genTime", 0.0)}
MuonVal::ScalarBranch< int > & m_out_gen_truthOrigin {m_tree.newScalar<int>("genTruthOrigin", -1)}
MuonVal::ScalarBranch< int > & m_out_gen_truthType {m_tree.newScalar<int>("genTruthType", -1)}
MuonVal::ScalarBranch< unsigned short > & m_out_gen_nHits {m_tree.newScalar<unsigned short>("genNHits",0)}
 Truth - hit count summary.
MuonVal::ScalarBranch< unsigned short > & m_out_gen_nRPCHits {m_tree.newScalar<unsigned short>("genNRpcHits",0)}
MuonVal::ScalarBranch< unsigned short > & m_out_gen_nMDTHits {m_tree.newScalar<unsigned short>("genNMdtHits",0)}
MuonVal::ScalarBranch< unsigned short > & m_out_gen_nTGCHits {m_tree.newScalar<unsigned short>("genNTgcHits",0)}
MuonVal::ScalarBranch< unsigned short > & m_out_gen_nNswHits {m_tree.newScalar<unsigned short>("genNNswHits",0)}
MuonVal::ScalarBranch< unsigned short > & m_out_gen_nMmHits {m_tree.newScalar<unsigned short>("genNMmHits",0)}
MuonVal::ScalarBranch< unsigned short > & m_out_gen_nSTGCHits {m_tree.newScalar<unsigned short>("genNsTgcHits",0)}
MuonVal::ScalarBranch< float > & m_out_gen_minYhit {m_tree.newScalar<float>("genMinYhit", 1.0)}
MuonVal::ScalarBranch< float > & m_out_gen_maxYhit {m_tree.newScalar<float>("genMaxYhit", -1.0)}
MuonVal::ScalarBranch< unsigned > & m_out_seed_n {m_tree.newScalar<unsigned>("nSeeds", 0)}
 ========== Seed block: Filled when we have one or multiple seeds ============= seed count
MuonVal::VectorBranch< unsigned short > & m_out_seed_hasPhiExtension {m_tree.newVector<unsigned short>("seedHasPhiExtension", false)}
MuonVal::VectorBranch< unsigned short > & m_out_seed_nMatchedHits {m_tree.newVector<unsigned short>("seedTruthMatchedHits", false)}
MuonVal::VectorBranch< float > & m_out_seed_y0 {m_tree.newVector<float>("seedY0", 0.0)}
MuonVal::VectorBranch< float > & m_out_seed_x0 {m_tree.newVector<float>("seedX0", 0.0)}
MuonVal::VectorBranch< float > & m_out_seed_tantheta {m_tree.newVector<float>("seedTanTheta", 0.0)}
MuonVal::VectorBranch< float > & m_out_seed_tanphi {m_tree.newVector<float>("seedTanPhi", 0.0)}
MuonVal::VectorBranch< float > & m_out_seed_minYhit {m_tree.newVector<float>("seedMinYhit", 1.0)}
MuonVal::VectorBranch< float > & m_out_seed_maxYhit {m_tree.newVector<float>("seedMaxYhit", -1.0)}
MuonVal::VectorBranch< unsigned short > & m_out_seed_nHits {m_tree.newVector<unsigned short>("seedNHits", 0)}
MuonVal::VectorBranch< unsigned short > & m_out_seed_nPrecHits {m_tree.newVector<unsigned short>("seedNPrecHits", 0)}
MuonVal::VectorBranch< unsigned short > & m_out_seed_nEtaHits {m_tree.newVector<unsigned short>("seedNEtaHits", 0)}
MuonVal::VectorBranch< unsigned short > & m_out_seed_nPhiHits {m_tree.newVector<unsigned short>("seedNPhiHits", 0)}
MuonVal::VectorBranch< unsigned short > & m_out_seed_nMdt {m_tree.newVector<unsigned short>("seedNMdtHits", 0)}
MuonVal::VectorBranch< unsigned short > & m_out_seed_nRpc {m_tree.newVector<unsigned short>("seedNRpcHits", 0)}
MuonVal::VectorBranch< unsigned short > & m_out_seed_nTgc {m_tree.newVector<unsigned short>("seedNTgcHits", 0)}
MuonVal::VectorBranch< unsigned short > & m_out_seed_nsTgcStrip {m_tree.newVector<unsigned short>("seedNsTgcStripHits", 0)}
MuonVal::VectorBranch< unsigned short > & m_out_seed_nsTgcWire {m_tree.newVector<unsigned short>("seedNsTgcWireHits", 0)}
MuonVal::VectorBranch< unsigned short > & m_out_seed_nsTgcPad {m_tree.newVector<unsigned short>("seedNsTgcPadHits", 0)}
MuonVal::VectorBranch< unsigned short > & m_out_seed_nMmEta {m_tree.newVector<unsigned short>("seedNMmHitsEta", 0)}
MuonVal::VectorBranch< unsigned short > & m_out_seed_nMmStereo {m_tree.newVector<unsigned short>("seedNMmHitsStereo", 0)}
MuonVal::VectorBranch< unsigned short > & m_out_seed_nTrueHits {m_tree.newVector<unsigned short>("seedNTrueHits", 0)}
 Labelled hits from the pattern visualization tool.
MuonVal::VectorBranch< unsigned short > & m_out_seed_nTruePrecHits {m_tree.newVector<unsigned short>("seedNTruePrecHits", 0)}
MuonVal::VectorBranch< unsigned short > & m_out_seed_nTruePhiHits {m_tree.newVector<unsigned short>("seedNTruePhiHits", 0)}
MuonVal::VectorBranch< unsigned short > & m_out_seed_nTrueEtaHits {m_tree.newVector<unsigned short>("seedNTrueEtaHits", 0)}
MuonVal::VectorBranch< unsigned char > & m_out_seed_ledToSegment {m_tree.newVector<unsigned char>("seedConvertedToSegment",false)}
MuonVal::ScalarBranch< unsigned > & m_out_segment_n {m_tree.newScalar<unsigned>("nSegments", 0)}
 ========== Segment block: Filled when we have one or multiple segments =============
MuonVal::VectorBranch< float > & m_out_segment_chi2 {m_tree.newVector<float>("segmentChi2", -1.)}
MuonVal::VectorBranch< uint16_t > & m_out_segment_nDoF {m_tree.newVector<uint16_t>("segmentNdoF", 0)}
MuonVal::VectorBranch< bool > & m_out_segment_hasTimeFit {m_tree.newVector<bool>("segmentHasTimeFit", false)}
MuonVal::VectorBranch< uint16_t > & m_out_segment_fitIter {m_tree.newVector<uint16_t>("segmentFitIterations", 0)}
MuonVal::VectorBranch< bool > & m_out_segment_hasPhi {m_tree.newVector<bool>("segmentHasPhiHits", false)}
MuonVal::VectorBranch< float > & m_out_segment_y0 {m_tree.newVector<float>("segmentY0", 0.0)}
MuonVal::VectorBranch< float > & m_out_segment_x0 {m_tree.newVector<float>("segmentX0", 0.0)}
MuonVal::VectorBranch< float > & m_out_segment_tantheta {m_tree.newVector<float>("segmentTanTheta", 0.0)}
MuonVal::VectorBranch< float > & m_out_segment_tanphi {m_tree.newVector<float>("segmentTanPhi", 0.0)}
MuonVal::VectorBranch< float > & m_out_segment_time {m_tree.newVector<float>("segmentTime", 0.)}
MuonVal::VectorBranch< float > & m_out_segment_err_y0 {m_tree.newVector<float>("segmentErrY0", -1.0)}
MuonVal::VectorBranch< float > & m_out_segment_err_x0 {m_tree.newVector<float>("segmentErrX0", -1.0)}
MuonVal::VectorBranch< float > & m_out_segment_err_tantheta {m_tree.newVector<float>("segmentErrTanTheta", -1.0)}
MuonVal::VectorBranch< float > & m_out_segment_err_tanphi {m_tree.newVector<float>("segmentErrTanPhi", -1.0)}
MuonVal::VectorBranch< float > & m_out_segment_err_time {m_tree.newVector<float>("segmentErrTime", -1.0)}
MuonVal::VectorBranch< unsigned short > & m_out_segment_truthMatchedHits {m_tree.newVector<unsigned short>("segmentTruthMatchedHits", 0)}
MuonVal::VectorBranch< unsigned short > & m_out_segment_nMdtHits {m_tree.newVector<unsigned short>("segmentNMdtHits", 0)}
MuonVal::VectorBranch< unsigned short > & m_out_segment_nRpcEtaHits {m_tree.newVector<unsigned short>("segmentNRpcEtaHits", 0)}
MuonVal::VectorBranch< unsigned short > & m_out_segment_nRpcPhiHits {m_tree.newVector<unsigned short>("segmentNRpcPhiHits", 0)}
MuonVal::VectorBranch< unsigned short > & m_out_segment_nTgcEtaHits {m_tree.newVector<unsigned short>("segmentNTgcEtaHits", 0)}
MuonVal::VectorBranch< unsigned short > & m_out_segment_nTgcPhiHits {m_tree.newVector<unsigned short>("segmentNTgcPhiHits", 0)}
MuonVal::VectorBranch< unsigned short > & m_out_segment_nMmEtaHits {m_tree.newVector<unsigned short>("segmentNMmEtaHits", 0)}
MuonVal::VectorBranch< unsigned short > & m_out_segment_nMmStereoHits {m_tree.newVector<unsigned short>("segmentNMmStereoHits", 0)}
MuonVal::VectorBranch< unsigned short > & m_out_segment_nsTgcStripHits {m_tree.newVector<unsigned short>("segmentNsTgcStripHits", 0)}
MuonVal::VectorBranch< unsigned short > & m_out_segment_nsTgcWireHits {m_tree.newVector<unsigned short>("segmentNsTgcWireHits", 0)}
MuonVal::VectorBranch< unsigned short > & m_out_segment_nsTgcPadpHits {m_tree.newVector<unsigned short>("segmentNsTgcPadHits", 0)}
MuonVal::VectorBranch< unsigned short > & m_out_segment_nTrueHits {m_tree.newVector<unsigned short>("segmentNTrueHits", 0)}
 Labelled hits from the pattern visualization tool.
MuonVal::VectorBranch< unsigned short > & m_out_segment_nTruePrecHits {m_tree.newVector<unsigned short>("segmentNTruePrecHits", 0)}
MuonVal::VectorBranch< unsigned short > & m_out_segment_nTruePhiHits {m_tree.newVector<unsigned short>("segmentNTruePhiHits", 0)}
MuonVal::VectorBranch< unsigned short > & m_out_segment_nTrueEtaHits {m_tree.newVector<unsigned short>("segmentNTrueEtaHits", 0)}
MuonVal::VectorBranch< float > & m_out_segment_minYhit {m_tree.newVector<float>("segmentMinYhit", 1.0)}
MuonVal::VectorBranch< float > & m_out_segment_maxYhit {m_tree.newVector<float>("segmentMaxYhit", -1.0)}
ToolHandle< MuonValR4::IPatternVisualizationToolm_visionTool {this, "VisualizationTool", ""}
 Pattern visualization tool.
const MuonGMR4::MuonDetectorManagerm_detMgr {nullptr}
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

Definition at line 38 of file MuonHoughTransformTester.h.

Member Typedef Documentation

◆ EffMap_t

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

Typedef for convenience.

Definition at line 207 of file AthHistogramming.h.

◆ GraphMap_t

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

Typedef for convenience.

Definition at line 221 of file AthHistogramming.h.

◆ hash_t

typedef uint32_t AthHistogramming::hash_t
privateinherited

typedef for the internal hash

Definition at line 171 of file AthHistogramming.h.

◆ HistMap_t

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

Typedef for convenience.

Definition at line 200 of file AthHistogramming.h.

◆ StoreGateSvc_t

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

Definition at line 388 of file AthCommonDataStore.h.

◆ TreeMap_t

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

Typedef for convenience.

Definition at line 214 of file AthHistogramming.h.

◆ TruthHitCol

Constructor & Destructor Documentation

◆ ~MuonHoughTransformTester()

virtual MuonValR4::MuonHoughTransformTester::~MuonHoughTransformTester ( )
virtualdefault

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

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

◆ book() [2/8]

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

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

◆ book() [3/8]

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

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

Definition at line 305 of file AthHistogramming.h.

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

◆ book() [4/8]

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

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

Definition at line 405 of file AthHistogramming.h.

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

◆ book() [5/8]

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

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

Definition at line 358 of file AthHistogramming.h.

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

◆ book() [6/8]

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

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

Definition at line 348 of file AthHistogramming.h.

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

◆ book() [7/8]

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

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

Definition at line 326 of file AthHistogramming.h.

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

◆ book() [8/8]

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

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

Definition at line 316 of file AthHistogramming.h.

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

◆ bookGetPointer() [1/8]

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

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

Definition at line 283 of file AthHistogramming.h.

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

◆ bookGetPointer() [2/8]

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

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

Definition at line 412 of file AthHistogramming.cxx.

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

◆ bookGetPointer() [3/8]

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

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

Definition at line 262 of file AthHistogramming.h.

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

◆ bookGetPointer() [4/8]

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

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

Definition at line 297 of file AthHistogramming.cxx.

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

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

◆ bookGetPointer() [6/8]

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

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

Definition at line 295 of file AthHistogramming.h.

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

◆ bookGetPointer() [7/8]

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

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

Definition at line 83 of file AthHistogramming.cxx.

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

◆ bookGetPointer() [8/8]

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

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

Definition at line 274 of file AthHistogramming.h.

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

◆ buildBookingString()

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

Method to build individual booking string.

Definition at line 545 of file AthHistogramming.cxx.

550{
551 // Massage the final string to book things
552 if(tDir.empty()) tDir = m_rootDir;
553 size_t pos = histName.rfind('/');
554 if(pos != std::string::npos){
555 tDir+='/';
556 tDir.append(histName, 0,pos);
557 histName.erase(0,pos+1);
558 };
559 if(stream.empty()) stream = m_streamName;
560
561 if(usePrefixPostfix){
562 bookingString = "/"+stream+"/"+tDir+"/"+m_histNamePrefix+histName+m_histNamePostfix;
563 } else {
564 bookingString = "/"+stream+"/"+tDir+"/"+histName;
565 }
566 while(bookingString.find("//") != std::string::npos){
567 this->myReplace(bookingString,"//","/");
568 }
569
570 return;
571}
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 51 of file AthHistogramming.cxx.

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

◆ countOnSameSide()

unsigned int MuonValR4::MuonHoughTransformTester::countOnSameSide ( const ActsTrk::GeometryContext & gctx,
const xAOD::MuonSegment & truthSeg,
const MuonR4::Segment & recoSeg ) const
private

Calculates how many measurements from the segment fit have the same drift sign as when evaluated with the truth parameters.

Parameters
gctxGeometry context to fetch the alignment constants
truthSegReference to the truth segment
recoSegReference to the reco segment.

Definition at line 84 of file MuonHoughTransformTester.cxx.

86 {
87 unsigned int same{0};
88 using namespace SegmentFit;
89 const auto[truePos, trueDir] = makeLine(localSegmentPars(truthSeg));
90 const auto[recoPos, recoDir] = makeLine(localSegmentPars(gctx, recoSeg));
91 const std::vector<int> truthSigns = SeedingAux::strawSigns(truePos, trueDir, recoSeg.measurements());
92 const std::vector<int> recoSigns = SeedingAux::strawSigns(recoPos, recoDir, recoSeg.measurements());
93 for (unsigned int s = 0 ; s < truthSigns.size(); ++s) {
94 same += (truthSigns[s] != 0) && truthSigns[s] == recoSigns[s];
95 }
96 return same;
97 }
const MeasVec & measurements() const
Returns the associated measurements.
Parameters localSegmentPars(const xAOD::MuonSegment &seg)
Returns the localSegPars decoration from a xAODMuon::Segment.
std::pair< Amg::Vector3D, Amg::Vector3D > makeLine(const Parameters &pars)
Returns the parsed parameters into an Eigen line parametrization.

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

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

◆ 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 MuonValR4::MuonHoughTransformTester::execute ( )
overridevirtual

Definition at line 222 of file MuonHoughTransformTester.cxx.

222 {
223
224 const EventContext & ctx = Gaudi::Hive::currentContext();
225
226 const ActsTrk::GeometryContext* gctxPtr{nullptr};
227 ATH_CHECK(SG::get(gctxPtr, m_geoCtxKey, ctx));
228 const ActsTrk::GeometryContext& gctx{*gctxPtr};
229
230 ConstDataVector<MuonR4::SegmentSeedContainer> segmentSeeds{SG::VIEW_ELEMENTS};
231 for (const SG::ReadHandleKey<SegmentSeedContainer>& key : m_patternSeedKeys) {
232 const SegmentSeedContainer* readSegmentSeeds{nullptr};
233 ATH_CHECK(SG::get(readSegmentSeeds, key, ctx));
234 segmentSeeds.insert(segmentSeeds.end(), readSegmentSeeds->begin(), readSegmentSeeds->end());
235 }
236
237 const xAOD::MuonSegmentContainer* truthSegments{nullptr};
238 ATH_CHECK(SG::get(truthSegments, m_truthSegmentKey, ctx));
239
240 const xAOD::MuonSegmentContainer* recoSegments{nullptr};
241 ATH_CHECK(SG::get(recoSegments, m_recoSegKey, ctx));
242
243 ATH_MSG_DEBUG("Succesfully retrieved input collections. Seeds: "<<segmentSeeds.size()
244 <<", segments: "<<recoSegments->size()
245 <<", truth segments: "<<(truthSegments? truthSegments->size() : -1)
246 <<".");
247 std::vector<ObjectMatching> objects = matchWithTruth(gctx, *segmentSeeds.asDataVector(),
248 *recoSegments, truthSegments);
249 for (const ObjectMatching& obj : objects) {
250 fillChamberInfo(obj.chamber);
251 fillSeedInfo(obj);
252 fillSegmentInfo(gctx, obj);
253 if(m_isMC) fillTruthInfo(gctx, obj.truthSegment);
254 ATH_CHECK(m_tree.fill(ctx));
255 }
256 return StatusCode::SUCCESS;
257 }
#define ATH_CHECK
Evaluate an expression and check for errors.
#define ATH_MSG_DEBUG(x)
iterator end() noexcept
Return an iterator pointing past the end of the collection.
iterator insert(iterator position, value_type pElem)
Add a new element to the collection.
const DV * asDataVector() const
Return a pointer to this object, as a const DataVector.
size_type size() const noexcept
Returns the number of elements in the collection.
void fillSegmentInfo(const ActsTrk::GeometryContext &gctx, const ObjectMatching &obj)
Fill the info assciated to the segment.
void fillTruthInfo(const ActsTrk::GeometryContext &gctx, const xAOD::MuonSegment *truthSegment)
Fill the associated truth information into the tree.
SG::ReadHandleKey< xAOD::MuonSegmentContainer > m_recoSegKey
Key to the xAOD::MuonSegment container.
SG::ReadHandleKeyArray< MuonR4::SegmentSeedContainer > m_patternSeedKeys
List of the two segment seed containers from which the segments are buiit (Complets the pattern findi...
std::vector< ObjectMatching > matchWithTruth(const ActsTrk::GeometryContext &gctx, const MuonR4::SegmentSeedContainer &seedContainer, const xAOD::MuonSegmentContainer &segmentContainer, const xAOD::MuonSegmentContainer *truthSegments) const
void fillChamberInfo(const MuonGMR4::SpectrometerSector *chamber)
Fill the current chamber info into the output.
SG::ReadHandleKey< xAOD::MuonSegmentContainer > m_truthSegmentKey
Key to the truth segment.
void fillSeedInfo(const ObjectMatching &obj)
Fill the info associated to the seed.
SG::ReadHandleKey< ActsTrk::GeometryContext > m_geoCtxKey
Tracking geometry context.
DataVector< SegmentSeed > SegmentSeedContainer
@ VIEW_ELEMENTS
this data object is a view, it does not own its elmts
const T * get(const ReadCondHandleKey< T > &key, const EventContext &ctx)
Convenience function to retrieve an object given a ReadCondHandleKey.
MuonSegmentContainer_v1 MuonSegmentContainer
Definition of the current "MuonSegment container version".

◆ 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

◆ fillBucketInfo()

void MuonValR4::MuonHoughTransformTester::fillBucketInfo ( const MuonR4::SpacePointBucket & bucket)
private

Fill the hit summary info of the associated bucket.

Definition at line 328 of file MuonHoughTransformTester.cxx.

328 {
329 m_out_bucketEnd = bucket.coveredMax();
330 m_out_bucketStart = bucket.coveredMin();
331 m_out_nSpacePoints = bucket.size();
332 double maxHitGap{-1.f};
333 for (auto itr = bucket.begin(); itr != bucket.end(); ++itr){
334 if (!(*itr)->measuresEta()) {
335 continue;
336 }
337 auto itr1 = std::find_if(itr+1, bucket.end(),[](const SpacePointBucket::value_type& sp){
338 return sp->measuresEta();
339 });
340 if (itr1 == bucket.end()){
341 break;
342 }
343 maxHitGap = std::max(maxHitGap, std::abs( (*itr)->localPosition().y() -
344 (*itr1)->localPosition().y()));
345 }
346 m_out_bucketHitGap = maxHitGap;
347 m_out_nPrecSpacePoints = std::ranges::count_if(bucket, [](const SpacePointBucket::value_type& sp){
348 return isPrecisionHit(*sp);
349 });
350 m_out_nPhiSpacePoints = std::ranges::count_if(bucket, [](const SpacePointBucket::value_type& sp){
351 return sp->measuresPhi();
352 });
353 if (!m_visionTool.isEnabled()){
354 return;
355 }
356 m_out_nTrueSpacePoints = std::ranges::count_if(bucket,[this](const SpacePointBucket::value_type& sp){
357 return m_visionTool->isLabeled(*sp);
358 });
359 m_out_nTruePrecSpacePoints = std::ranges::count_if(bucket,[this](const SpacePointBucket::value_type& sp){
360 return isPrecisionHit(*sp) && m_visionTool->isLabeled(*sp);
361 });
362 m_out_nTruePhiSpacePoints = std::ranges::count_if(bucket,[this](const SpacePointBucket::value_type& sp){
363 return sp->measuresPhi() && m_visionTool->isLabeled(*sp);
364 });
365 }
static Double_t sp
MuonVal::ScalarBranch< unsigned char > & m_out_nTruePhiSpacePoints
Number of phi hits in the bucket.
ToolHandle< MuonValR4::IPatternVisualizationTool > m_visionTool
Pattern visualization tool.
MuonVal::ScalarBranch< unsigned char > & m_out_nPhiSpacePoints
Number of phi hits in the bucket.
MuonVal::ScalarBranch< float > & m_out_bucketEnd
MuonVal::ScalarBranch< unsigned char > & m_out_nTrueSpacePoints
Number of all space points in the bucket.
MuonVal::ScalarBranch< unsigned char > & m_out_nPrecSpacePoints
Number of precision hits in the bucket.
MuonVal::ScalarBranch< unsigned char > & m_out_nTruePrecSpacePoints
Number of precision hits in the bucket.
MuonVal::ScalarBranch< float > & m_out_bucketHitGap
MuonVal::ScalarBranch< unsigned char > & m_out_nSpacePoints
Number of all space points in the bucket.
MuonVal::ScalarBranch< float > & m_out_bucketStart
bool isPrecisionHit(const SpacePoint &hit)
Returns whether the uncalibrated spacepoint is a precision hit (Mdt, micromegas, stgc strips).

◆ fillChamberInfo()

void MuonValR4::MuonHoughTransformTester::fillChamberInfo ( const MuonGMR4::SpectrometerSector * chamber)
private

Fill the current chamber info into the output.

Parameters
chamberPointer to the reference sector to which the view belongs to

Definition at line 258 of file MuonHoughTransformTester.cxx.

258 {
259 m_out_chamberIndex = Acts::toUnderlying(msSector->chamberIndex());
260 m_out_stationSide = msSector->side();
261 m_out_stationPhi = msSector->stationPhi();
262 }
MuonVal::ScalarBranch< int > & m_out_stationPhi
phi index of the station
MuonVal::ScalarBranch< int > & m_out_chamberIndex
====== Common block: Filled for all entries ===========
MuonVal::ScalarBranch< short > & m_out_stationSide
+1 for A-, -1 of C-side

◆ fillSeedInfo()

void MuonValR4::MuonHoughTransformTester::fillSeedInfo ( const ObjectMatching & obj)
private

Fill the info associated to the seed.

Parameters
objPointer to the matching object connecting the seeds & segment & truth

Definition at line 366 of file MuonHoughTransformTester.cxx.

366 {
367
368 m_out_seed_n = obj.matchedSeeds.size();
369 for (const auto [iseed, seed] : Acts::enumerate(obj.matchedSeeds)){
370 if (iseed ==0) {
371 fillBucketInfo(*seed->parentBucket());
372 }
373 double minYhit = m_out_bucketEnd.getVariable();
374 double maxYhit = m_out_bucketStart.getVariable();
375 for (const SpacePoint* hit : seed->getHitsInMax()){
376 minYhit = std::min(hit->localPosition().y(),minYhit);
377 maxYhit = std::max(hit->localPosition().y(),maxYhit);
378 }
379 m_out_seed_minYhit.push_back(minYhit);
380 m_out_seed_maxYhit.push_back(maxYhit);
381
382 m_out_seed_hasPhiExtension.push_back(seed->hasPhiExtension());
383 m_out_seed_nMatchedHits.push_back(countMatched(obj.truthSegment, *seed));
384 m_out_seed_y0.push_back(seed->interceptY());
385 m_out_seed_tantheta.push_back(seed->tanBeta());
386 if (seed->hasPhiExtension()){
387 m_out_seed_x0.push_back(seed->interceptX());
388 m_out_seed_tanphi.push_back(seed->tanAlpha());
389 } else{
390 m_out_seed_x0.push_back(-999);
391 m_out_seed_tanphi.push_back(-999);
392 }
393
394 m_out_seed_nHits.push_back(seed->getHitsInMax().size());
395 unsigned nMdtSeed{0}, nRpcSeed{0}, nTgcSeed{0}, nMmEtaSeed{0}, nMmStereoSeed{0},
396 nsTgcStripSeed{0}, nsTgcWireSeed{0}, nsTgcPadSeed{0};
397 unsigned nPrecHits{0}, nEtaHits{0}, nPhiHits{0}, nTrueHits{0}, nTruePrecHits{0}, nTrueEtaHits{0}, nTruePhiHits{0};
398 std::vector<unsigned char> treeIdxs{};
399
400 for (const HoughHitType & houghSP: seed->getHitsInMax()){
402 unsigned treeIdx = m_spTester->push_back(*houghSP);
403 treeIdxs.push_back(treeIdx);
404 }
405 nPrecHits += isPrecisionHit(*houghSP);
406 nPhiHits += houghSP->measuresPhi();
407 nEtaHits += houghSP->measuresEta();
408
409 if (m_visionTool.isEnabled()) {
410 nTrueHits += m_visionTool->isLabeled(*houghSP);
411 nTruePrecHits += m_visionTool->isLabeled(*houghSP) && isPrecisionHit(*houghSP);
412 nTruePhiHits += m_visionTool->isLabeled(*houghSP) && houghSP->measuresPhi();
413 nTrueEtaHits += m_visionTool->isLabeled(*houghSP) && houghSP->measuresEta();
414 }
415 switch (houghSP->type()) {
417 ++nMdtSeed;
418 break;
420 nRpcSeed+=houghSP->measuresEta();
421 nRpcSeed+=houghSP->measuresPhi();
422 break;
424 nTgcSeed+=houghSP->measuresEta();
425 nTgcSeed+=houghSP->measuresPhi();
426 break;
428 const Identifier sTgc = houghSP->primaryMeasurement()->identify();
429 const Identifier sTgc2 = houghSP->secondaryMeasurement() ?
430 houghSP->secondaryMeasurement()->identify() : Identifier{};
431 const sTgcIdHelper& idHelper{m_idHelperSvc->stgcIdHelper()};
432 const int primType = idHelper.channelType(sTgc);
433 const int secType = idHelper.channelType(sTgc2);
434 nsTgcStripSeed += primType == sTgcIdHelper::sTgcChannelTypes::Strip;
435 nsTgcWireSeed += primType == sTgcIdHelper::sTgcChannelTypes::Wire;
436 nsTgcPadSeed += primType == sTgcIdHelper::sTgcChannelTypes::Pad;
437
438 nsTgcWireSeed += secType == sTgcIdHelper::sTgcChannelTypes::Wire;
439 nsTgcPadSeed += primType != sTgcIdHelper::sTgcChannelTypes::Pad &&
441 break;
443 if (m_idHelperSvc->mmIdHelper().isStereo(houghSP->identify())) {
444 ++nMmEtaSeed;
445 } else {
446 ++nMmStereoSeed;
447 }
448 break;
449 }default:
450 ATH_MSG_WARNING("Technology "<<houghSP->identify() <<" not yet implemented");
451 }
452 }
453 m_out_seed_nMdt.push_back(nMdtSeed);
454 m_out_seed_nRpc.push_back(nRpcSeed);
455 m_out_seed_nTgc.push_back(nTgcSeed);
456
457 m_out_seed_nMmEta.push_back(nMmEtaSeed);
458 m_out_seed_nMmStereo.push_back(nMmStereoSeed);
459
460 m_out_seed_nsTgcStrip.push_back(nsTgcStripSeed);
461 m_out_seed_nsTgcWire.push_back(nsTgcWireSeed);
462 m_out_seed_nsTgcPad.push_back(nsTgcPadSeed);
463
464 m_out_seed_nPrecHits.push_back(nPrecHits);
465 m_out_seed_nEtaHits.push_back(nEtaHits);
466 m_out_seed_nPhiHits.push_back(nPhiHits);
467
468 m_out_seed_nTrueHits.push_back(nTrueHits);
469 m_out_seed_nTruePrecHits.push_back(nTruePrecHits);
470 m_out_seed_nTrueEtaHits.push_back(nTrueEtaHits);
471 m_out_seed_nTruePhiHits.push_back(nTruePhiHits);
472
473 m_out_seed_ledToSegment.push_back(obj.matchedSeedFoundSegment.at(iseed));
474 if (m_writeSpacePoints) {
475 m_spMatchedToPattern[iseed] = std::move(treeIdxs);
476
477 }
478 }
479 }
#define ATH_MSG_WARNING(x)
MuonVal::VectorBranch< float > & m_out_seed_maxYhit
MuonVal::VectorBranch< unsigned short > & m_out_seed_nHits
MuonVal::VectorBranch< unsigned short > & m_out_seed_nMmStereo
MuonVal::VectorBranch< unsigned char > & m_out_seed_ledToSegment
MuonVal::VectorBranch< unsigned short > & m_out_seed_nTrueEtaHits
ServiceHandle< Muon::IMuonIdHelperSvc > m_idHelperSvc
MuonVal::VectorBranch< float > & m_out_seed_x0
MuonVal::VectorBranch< unsigned short > & m_out_seed_nTruePrecHits
MuonVal::VectorBranch< float > & m_out_seed_y0
MuonVal::VectorBranch< unsigned short > & m_out_seed_hasPhiExtension
MuonVal::MatrixBranch< unsigned char > & m_spMatchedToPattern
Branch indicating which space points in the tree are associated to the i-th pattern.
MuonVal::VectorBranch< unsigned short > & m_out_seed_nsTgcWire
MuonVal::VectorBranch< unsigned short > & m_out_seed_nsTgcStrip
MuonVal::VectorBranch< unsigned short > & m_out_seed_nEtaHits
void fillBucketInfo(const MuonR4::SpacePointBucket &bucket)
Fill the hit summary info of the associated bucket.
MuonVal::VectorBranch< float > & m_out_seed_tanphi
MuonVal::ScalarBranch< unsigned > & m_out_seed_n
========== Seed block: Filled when we have one or multiple seeds ============= seed count
MuonVal::VectorBranch< unsigned short > & m_out_seed_nTruePhiHits
MuonVal::VectorBranch< unsigned short > & m_out_seed_nMdt
MuonVal::VectorBranch< float > & m_out_seed_minYhit
MuonVal::VectorBranch< unsigned short > & m_out_seed_nMmEta
MuonVal::VectorBranch< unsigned short > & m_out_seed_nRpc
MuonVal::VectorBranch< unsigned short > & m_out_seed_nPrecHits
MuonVal::VectorBranch< float > & m_out_seed_tantheta
MuonVal::VectorBranch< unsigned short > & m_out_seed_nsTgcPad
MuonVal::VectorBranch< unsigned short > & m_out_seed_nTrueHits
Labelled hits from the pattern visualization tool.
MuonVal::VectorBranch< unsigned short > & m_out_seed_nPhiHits
MuonVal::VectorBranch< unsigned short > & m_out_seed_nMatchedHits
MuonVal::VectorBranch< unsigned short > & m_out_seed_nTgc
std::shared_ptr< SpacePointTesterModule > m_spTester
Branch dumping all the space points from the difference buckets.
int channelType(const Identifier &id) const
@ sTgc
Micromegas (NSW).
const SpacePoint * HoughHitType

◆ fillSegmentInfo()

void MuonValR4::MuonHoughTransformTester::fillSegmentInfo ( const ActsTrk::GeometryContext & gctx,
const ObjectMatching & obj )
private

Fill the info assciated to the segment.

Parameters
objPointer to the matching object connecting the seeds & segment & truth

Definition at line 481 of file MuonHoughTransformTester.cxx.

482 {
483 using namespace SegmentFit;
484
485 m_out_segment_n = obj.matchedSegments.size();
486 for (const Segment* segment : obj.matchedSegments) {
487 m_out_segment_hasPhi.push_back(std::ranges::any_of(segment->measurements(),
488 [](const auto& meas){ return meas->measuresPhi();}));
489 m_out_segment_fitIter.push_back(segment->nFitIterations());
490 m_out_segment_truthMatchedHits.push_back(countMatched(obj.truthSegment, *segment));
491 m_out_segment_chi2.push_back(segment->chi2());
492 m_out_segment_nDoF.push_back(segment->nDoF());
493 m_out_segment_hasTimeFit.push_back(segment->hasTimeFit());
494
495 m_out_segment_err_x0.push_back(segment->covariance()(Acts::toUnderlying(ParamDefs::x0), Acts::toUnderlying(ParamDefs::x0)));
496 m_out_segment_err_y0.push_back(segment->covariance()(Acts::toUnderlying(ParamDefs::y0), Acts::toUnderlying(ParamDefs::y0)));
497 m_out_segment_err_tantheta.push_back(segment->covariance()(Acts::toUnderlying(ParamDefs::theta), Acts::toUnderlying(ParamDefs::theta)));
498 m_out_segment_err_tanphi.push_back(segment->covariance()(Acts::toUnderlying(ParamDefs::phi), Acts::toUnderlying(ParamDefs::phi)));
499 m_out_segment_err_time.push_back(segment->covariance()(Acts::toUnderlying(ParamDefs::t0), Acts::toUnderlying(ParamDefs::t0)));
500 const auto [locPos, locDir] = makeLine(localSegmentPars(gctx, *segment));
501 m_out_segment_tanphi.push_back(houghTanAlpha(locDir));
502 m_out_segment_tantheta.push_back(houghTanBeta(locDir));
503 m_out_segment_y0.push_back(locPos.y());
504 m_out_segment_x0.push_back(locPos.x());
505 m_out_segment_time.push_back(segment->segementT0() + segment->position().mag() * c_inv);
506
507 unsigned nMdtHits{0}, nRpcEtaHits{0}, nRpcPhiHits{0}, nTgcEtaHits{0}, nTgcPhiHits{0},
508 nMmEtaHits{0}, nMmStereoHits{0}, nStgcStripHits{0},nStgcWireHits{0}, nStgcPadHits{0};
509 unsigned nTrueHits{0}, nTruePrecHits{0}, nTrueEtaHits{0}, nTruePhiHits{0};
510
511 double minYhit = std::numeric_limits<double>::max();
512 double maxYhit = -1 * std::numeric_limits<double>::max();
513
514 std::vector<unsigned char> matched;
515 for (const auto & meas : segment->measurements()){
516 // skip dummy measurement from beam spot constraint
517 ATH_MSG_VERBOSE(__func__<<"() - "<<__LINE__<<" Dump "<<(*meas));
518 if (meas->type() == xAOD::UncalibMeasType::Other) {
519 continue;
520 }
521 minYhit = std::min(meas->localPosition().y(),minYhit);
522 maxYhit = std::max(meas->localPosition().y(),maxYhit);
523 if (m_writeSpacePoints) {
524 unsigned treeIdx = m_spTester->push_back(*meas->spacePoint());
525 if (treeIdx >= matched.size()){
526 matched.resize(treeIdx +1);
527 }
528 matched[treeIdx] = true;
529 }
530 if (m_visionTool.isEnabled()) {
531 nTrueHits += m_visionTool->isLabeled(*meas->spacePoint());
532 nTruePrecHits += isPrecisionHit(*meas) && m_visionTool->isLabeled(*meas->spacePoint());
533 nTrueEtaHits += meas->measuresEta() && m_visionTool->isLabeled(*meas->spacePoint());
534 nTruePhiHits += meas->measuresPhi() && m_visionTool->isLabeled(*meas->spacePoint());
535 }
536 switch (meas->type()) {
538 ++nMdtHits;
539 break;
541 nRpcEtaHits += meas->measuresEta();
542 nRpcPhiHits += meas->measuresPhi();
543 break;
545 nTgcEtaHits += meas->measuresEta();
546 nTgcPhiHits += meas->measuresPhi();
547 break;
549 const MmIdHelper& idHelper{m_idHelperSvc->mmIdHelper()};
550 nMmEtaHits += !idHelper.isStereo(meas->spacePoint()->identify());
551 nMmStereoHits += !idHelper.isStereo(meas->spacePoint()->identify());
552 break;
554 const auto* prd = static_cast<const xAOD::sTgcMeasurement*>(meas->spacePoint()->primaryMeasurement());
555 const int primType = prd->channelType();
556 prd = dynamic_cast<const xAOD::sTgcMeasurement*>(meas->spacePoint()->secondaryMeasurement());
557 const int secType = (prd != nullptr ? prd->channelType() : -1);
558 nStgcStripHits += primType == sTgcIdHelper::sTgcChannelTypes::Strip;
559 nStgcWireHits += primType == sTgcIdHelper::sTgcChannelTypes::Wire;
560 nStgcPadHits += primType == sTgcIdHelper::sTgcChannelTypes::Pad;
561 nStgcWireHits += secType == sTgcIdHelper::sTgcChannelTypes::Wire;
562 nStgcPadHits += primType != secType && secType == sTgcIdHelper::sTgcChannelTypes::Pad;
563 break;
564 } default:
565 break;
566 }
567 }
568 m_out_segment_nMdtHits.push_back(nMdtHits);
569 m_out_segment_nRpcEtaHits.push_back(nRpcEtaHits);
570 m_out_segment_nRpcPhiHits.push_back(nRpcPhiHits);
571 m_out_segment_nTgcEtaHits.push_back(nTgcEtaHits);
572 m_out_segment_nTgcPhiHits.push_back(nTgcPhiHits);
573
574 m_out_segment_nMmEtaHits.push_back(nMmEtaHits);
575 m_out_segment_nMmStereoHits.push_back(nMmStereoHits);
576 m_out_segment_nsTgcStripHits.push_back(nStgcStripHits);
577 m_out_segment_nsTgcWireHits.push_back(nStgcWireHits);
578 m_out_segment_nsTgcPadpHits.push_back(nStgcPadHits);
579
580
581 m_out_segment_nTrueHits.push_back(nTrueHits);
582 m_out_segment_nTruePrecHits.push_back(nTruePrecHits);
583 m_out_segment_nTruePhiHits.push_back(nTruePhiHits);
584 m_out_segment_nTrueEtaHits.push_back(nTrueEtaHits);
585
586 m_out_segment_minYhit.push_back(minYhit);
587 m_out_segment_maxYhit.push_back(maxYhit);
588 if (m_writeSpacePoints) {
589 m_spMatchedToSegment.push_back(std::move(matched));
590 }
591 }
592 }
#define ATH_MSG_VERBOSE(x)
bool isStereo(const Identifier &id) const
MuonVal::VectorBranch< unsigned short > & m_out_segment_nMmStereoHits
MuonVal::VectorBranch< unsigned short > & m_out_segment_nRpcPhiHits
MuonVal::ScalarBranch< unsigned > & m_out_segment_n
========== Segment block: Filled when we have one or multiple segments =============
MuonVal::VectorBranch< float > & m_out_segment_err_y0
MuonVal::VectorBranch< float > & m_out_segment_err_tantheta
MuonVal::VectorBranch< unsigned short > & m_out_segment_truthMatchedHits
MuonVal::VectorBranch< unsigned short > & m_out_segment_nsTgcPadpHits
MuonVal::VectorBranch< unsigned short > & m_out_segment_nsTgcWireHits
MuonVal::VectorBranch< unsigned short > & m_out_segment_nsTgcStripHits
MuonVal::VectorBranch< unsigned short > & m_out_segment_nTruePrecHits
MuonVal::VectorBranch< unsigned short > & m_out_segment_nTgcPhiHits
MuonVal::VectorBranch< float > & m_out_segment_time
MuonVal::MatrixBranch< unsigned char > & m_spMatchedToSegment
Branch indicating which space points in the tree are associated to the i-th segment.
MuonVal::VectorBranch< bool > & m_out_segment_hasPhi
MuonVal::VectorBranch< float > & m_out_segment_tanphi
MuonVal::VectorBranch< float > & m_out_segment_err_tanphi
MuonVal::VectorBranch< float > & m_out_segment_chi2
MuonVal::VectorBranch< unsigned short > & m_out_segment_nTrueHits
Labelled hits from the pattern visualization tool.
MuonVal::VectorBranch< float > & m_out_segment_tantheta
MuonVal::VectorBranch< float > & m_out_segment_err_time
MuonVal::VectorBranch< unsigned short > & m_out_segment_nTruePhiHits
MuonVal::VectorBranch< float > & m_out_segment_minYhit
MuonVal::VectorBranch< float > & m_out_segment_x0
MuonVal::VectorBranch< bool > & m_out_segment_hasTimeFit
MuonVal::VectorBranch< unsigned short > & m_out_segment_nTgcEtaHits
MuonVal::VectorBranch< unsigned short > & m_out_segment_nMmEtaHits
MuonVal::VectorBranch< float > & m_out_segment_maxYhit
MuonVal::VectorBranch< unsigned short > & m_out_segment_nMdtHits
MuonVal::VectorBranch< unsigned short > & m_out_segment_nTrueEtaHits
MuonVal::VectorBranch< unsigned short > & m_out_segment_nRpcEtaHits
MuonVal::VectorBranch< float > & m_out_segment_err_x0
MuonVal::VectorBranch< uint16_t > & m_out_segment_nDoF
MuonVal::VectorBranch< uint16_t > & m_out_segment_fitIter
MuonVal::VectorBranch< float > & m_out_segment_y0
virtual sTgcChannelTypes channelType() const =0
Returns the channel type of the measurement (Pad/Wire/Strip).
double houghTanBeta(const Amg::Vector3D &v)
Returns the hough tanBeta [y] / [z].
double houghTanAlpha(const Amg::Vector3D &v)
: Returns the hough tanAlpha [x] / [z]
sTgcMeasurement_v1 sTgcMeasurement

◆ fillTruthInfo()

void MuonValR4::MuonHoughTransformTester::fillTruthInfo ( const ActsTrk::GeometryContext & gctx,
const xAOD::MuonSegment * truthSegment )
private

Fill the associated truth information into the tree.

Parameters
truthSegmentPointer to the truth parameters in form of a segment
gctxGeometry context for the alignment of the spectrometer sector

Definition at line 263 of file MuonHoughTransformTester.cxx.

264 {
265 if (!segment) {
266 return;
267 }
268 m_out_hasTruth = true;
269
270 const Amg::Vector3D segDir{segment->direction()};
271 static const SG::ConstAccessor<float> acc_pt{"pt"};
272 static const SG::ConstAccessor<float> acc_charge{"charge"};
273 // eta is interpreted as the eta-location
274 m_out_gen_Eta = segDir.eta();
275 m_out_gen_Phi = segDir.phi();
276 m_out_gen_Pt = acc_pt(*segment);
277 m_out_gen_Q = acc_charge(*segment);
278
279 const auto [chamberPos, chamberDir] = SegmentFit::makeLine(SegmentFit::localSegmentPars(*segment));
280 m_out_gen_nHits = segment->nPrecisionHits()+segment->nPhiLayers() + segment->nTrigEtaLayers();
281 using namespace Muon::MuonStationIndex;
282 m_out_gen_nMDTHits = segment->nPrecisionHits() * (segment->technology() == TechnologyIndex::MDT);
283 m_out_gen_nNswHits = segment->nPrecisionHits() * (segment->technology() != TechnologyIndex::MDT);
284 m_out_gen_nTGCHits = (segment->nPhiLayers() + segment->nTrigEtaLayers()) * !isBarrel(segment->chamberIndex());
285 m_out_gen_nRPCHits = (segment->nPhiLayers() + segment->nTrigEtaLayers()) * isBarrel(segment->chamberIndex());
286
287 unsigned nMMHits{0}, nSTGHits{0};
288 for (const xAOD::MuonSimHit* simHit : getMatchingSimHits(*segment)) {
289 const TechnologyIndex simIdx = m_idHelperSvc->technologyIndex(simHit->identify());
290 nMMHits += simIdx == TechnologyIndex::MM;
291 nSTGHits += simIdx == TechnologyIndex::STGC;
292 }
293 m_out_gen_nMmHits = nMMHits;
294 m_out_gen_nSTGCHits = nSTGHits;
295
296 m_out_gen_tantheta = houghTanBeta(chamberDir);
297 m_out_gen_tanphi = houghTanAlpha(chamberDir);
298 m_out_gen_y0 = chamberPos.y();
299 m_out_gen_x0 = chamberPos.x();
300 m_out_gen_time = segment->t0();
301
302 double minYhit = std::numeric_limits<double>::max();
303 double maxYhit = -1 * std::numeric_limits<double>::max();
304 for (const xAOD::MuonSimHit* hit : getMatchingSimHits(*segment)){
305 const Identifier hitId = hit->identify();
306 const MuonGMR4::MuonReadoutElement* RE = m_detMgr->getReadoutElement(hitId);
307 const IdentifierHash hash{m_idHelperSvc->isMdt(hitId) ? RE->measurementHash(hitId)
308 : RE->layerHash(hitId) };
309 const Amg::Transform3D localToChamber = RE->msSector()->globalToLocalTransform(gctx) * RE->localToGlobalTransform(gctx, hash);
310 const Amg::Vector3D chamberPos = localToChamber * xAOD::toEigen(hit->localPosition());
311 minYhit = std::min(chamberPos.y(), minYhit);
312 maxYhit = std::max(chamberPos.y(), maxYhit);
313 }
314 m_out_gen_minYhit = minYhit;
315 m_out_gen_maxYhit = maxYhit;
316
317 ATH_MSG_DEBUG("A true max on chamber index "<<m_out_chamberIndex.getVariable()<<" side "<<m_out_stationSide.getVariable()<<" phi "<<m_out_stationPhi.getVariable()<<" with "
318 <<m_out_gen_nMDTHits.getVariable()<<" MDT and "<<m_out_gen_nRPCHits.getVariable()+m_out_gen_nTGCHits.getVariable()<< " trigger hits is at "
319 <<m_out_gen_tantheta.getVariable()<<" and "<<m_out_gen_y0.getVariable());
320
322 if (truthMuon) {
323 using namespace xAOD::TruthHelpers;
326 }
327 }
const Amg::Transform3D & localToGlobalTransform(const ActsTrk::GeometryContext &ctx) const
Returns the transformation from the local coordinate system of the readout element into the global AT...
const SpectrometerSector * msSector() const
Returns the pointer to the envelope volume enclosing all chambers in the sector.
virtual IdentifierHash measurementHash(const Identifier &measId) const =0
The measurement hash is a continous numbering schema of all readout channels described by the specifi...
virtual IdentifierHash layerHash(const Identifier &measId) const =0
The layer hash removes the bits from the IdentifierHash corresponding to the measurement's channel nu...
Amg::Transform3D globalToLocalTransform(const ActsTrk::GeometryContext &gctx) const
Returns the global -> local transformation from the ATLAS global.
MuonVal::ScalarBranch< unsigned short > & m_out_gen_nMmHits
MuonVal::ScalarBranch< float > & m_out_gen_Pt
const MuonGMR4::MuonDetectorManager * m_detMgr
MuonVal::ScalarBranch< unsigned short > & m_out_gen_nHits
Truth - hit count summary.
MuonVal::ScalarBranch< unsigned short > & m_out_gen_nMDTHits
MuonVal::ScalarBranch< bool > & m_out_hasTruth
======= Truth block: Filled if we have a truth match. ============
MuonVal::ScalarBranch< int > & m_out_gen_truthType
MuonVal::ScalarBranch< unsigned short > & m_out_gen_nRPCHits
MuonVal::ScalarBranch< short > & m_out_gen_Q
MuonVal::ScalarBranch< float > & m_out_gen_y0
Truth - segment parameters.
MuonVal::ScalarBranch< float > & m_out_gen_tanphi
MuonVal::ScalarBranch< float > & m_out_gen_Phi
MuonVal::ScalarBranch< unsigned short > & m_out_gen_nSTGCHits
MuonVal::ScalarBranch< int > & m_out_gen_truthOrigin
MuonVal::ScalarBranch< float > & m_out_gen_maxYhit
MuonVal::ScalarBranch< float > & m_out_gen_minYhit
MuonVal::ScalarBranch< unsigned short > & m_out_gen_nNswHits
MuonVal::ScalarBranch< float > & m_out_gen_Eta
global particle properties
MuonVal::ScalarBranch< float > & m_out_gen_time
MuonVal::ScalarBranch< float > & m_out_gen_x0
MuonVal::ScalarBranch< float > & m_out_gen_tantheta
MuonVal::ScalarBranch< unsigned short > & m_out_gen_nTGCHits
Eigen::Affine3d Transform3D
Eigen::Matrix< double, 3, 1 > Vector3D
const xAOD::TruthParticle * getTruthMatchedParticle(const xAOD::MuonSegment &segment)
Returns the particle truth-matched to the segment.
std::unordered_set< const xAOD::MuonSimHit * > getMatchingSimHits(const xAOD::MuonSegment &segment)
: Returns all sim hits matched to a xAOD::MuonSegment
TechnologyIndex
enum to classify the different layers in the muon spectrometer
bool isBarrel(const ChIndex index)
Returns true if the chamber index points to a barrel chamber.
int getParticleTruthType(const xAOD::IParticle &p)
Return the particle's truth type (as defined by the MC Truth Classifier).
int getParticleTruthOrigin(const xAOD::IParticle &p)
Return the particle's truth origin (as defined by the MC Truth Classifier).
MuonSimHit_v1 MuonSimHit
Defined the version of the MuonSimHit.
Definition MuonSimHit.h:12
TruthParticle_v1 TruthParticle
Typedef to implementation.

◆ finalize()

StatusCode MuonValR4::MuonHoughTransformTester::finalize ( )
overridevirtual

Definition at line 218 of file MuonHoughTransformTester.cxx.

218 {
219 ATH_CHECK(m_tree.write());
220 return StatusCode::SUCCESS;
221 }

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

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

◆ hash()

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

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

Definition at line 430 of file AthHistogramming.h.

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

◆ hist()

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

Simplify the retrieval of registered histograms of any type.

Definition at line 183 of file AthHistogramming.cxx.

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

◆ hist2d()

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

Simplify the retrieval of registered 2-d histograms.

Definition at line 371 of file AthHistogramming.h.

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

◆ hist3d()

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

Simplify the retrieval of registered 3-d histograms.

Definition at line 388 of file AthHistogramming.h.

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

◆ histSvc()

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

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

Definition at line 113 of file AthHistogramAlgorithm.h.

114{
115 return m_histSvc;
116}

◆ initialize()

StatusCode MuonValR4::MuonHoughTransformTester::initialize ( )
overridevirtual

The collection of readHandle keys should be either 1 or 2

Definition at line 45 of file MuonHoughTransformTester.cxx.

45 {
46 ATH_CHECK(m_geoCtxKey.initialize());
47
48 {
49 int infoOpts = 0;
50 if (m_isMC) infoOpts = EventInfoBranch::isMC;
51 m_tree.addBranch(std::make_unique<EventInfoBranch>(m_tree, infoOpts));
52 }
53
54 ATH_CHECK(m_recoSegKey.initialize());
55 for (const std::string& recoLink : m_recoSegLinks) {
56 m_truthSegLinkKeys.emplace_back(m_recoSegKey, recoLink);
57 }
58 ATH_CHECK(m_truthSegmentKey.initialize(!m_truthSegmentKey.empty()));
59 for (const std::string& link: m_truthLinks) {
60 m_truthSegLinkKeys.emplace_back(m_truthSegmentKey, link);
61 }
63
65 ATH_CHECK(m_patternSeedKeys.initialize());
68 m_spTester = std::make_unique<SpacePointTesterModule>(m_tree, m_spKeys.front().key(), msgLevel());
69 m_tree.addBranch(m_spTester);
70 } else {
71 m_tree.disableBranch(m_spMatchedToPattern.name());
72 m_tree.disableBranch(m_spMatchedToSegment.name());
73 }
74 ATH_CHECK(m_tree.init(this));
75 ATH_CHECK(m_idHelperSvc.retrieve());
76 ATH_CHECK(detStore()->retrieve(m_detMgr));
77
78 ATH_CHECK(m_visionTool.retrieve(EnableTool{!m_visionTool.empty()}));
79 ATH_MSG_DEBUG("Succesfully initialised");
80 return StatusCode::SUCCESS;
81 }
const ServiceHandle< StoreGateSvc > & detStore() const
SG::ReadHandleKeyArray< MuonR4::SpacePointContainer > m_spKeys
List of the space point containers in the event legacy + NSW containers.
SG::ReadDecorHandleKeyArray< xAOD::MuonSegmentContainer > m_truthSegLinkKeys
Declare the dependencies on the decorations.
Gaudi::Property< std::vector< std::string > > m_truthLinks
Name of the decorations for the truth segment.
Gaudi::Property< std::vector< std::string > > m_recoSegLinks
name of the truth link decorations for the reco segment container
@ isMC
Flag determining whether the branch is simulation.

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

◆ matchWithTruth()

std::vector< ObjectMatching > MuonValR4::MuonHoughTransformTester::matchWithTruth ( const ActsTrk::GeometryContext & gctx,
const MuonR4::SegmentSeedContainer & seedContainer,
const xAOD::MuonSegmentContainer & segmentContainer,
const xAOD::MuonSegmentContainer * truthSegments ) const
private

Step 1 Loop over all reconstructed segments && fill them into the association table

Next loop over all seeds

Don't recycle the used seeds again

Find the best matching truth segment

find the best matching segment

Finally loop over all truth segments and push back the unused ones

Check that the segment is not used yet

Definition at line 99 of file MuonHoughTransformTester.cxx.

102 {
103 std::vector<ObjectMatching> allAssociations{};
104 std::unordered_set<const SegmentSeed*> usedSeeds{};
107 for (const xAOD::MuonSegment* recoSeg : segmentContainer) {
108 const MuonR4::Segment* segment = detailedSegment(*recoSeg);
109 assert(segment != nullptr);
110 std::vector<ObjectMatching>::iterator assoc_itr = allAssociations.end();
111 const xAOD::MuonSegment* truthSeg = getMatchedTruthSegment(*recoSeg);
112 if (truthSeg) {
113 assoc_itr = std::ranges::find_if(allAssociations, [truthSeg](const ObjectMatching& obj){
114 return obj.truthSegment == truthSeg;
115 });
116 }
117 if (assoc_itr == allAssociations.end()) {
118 ObjectMatching & newObj = allAssociations.emplace_back();
119 newObj.chamber = m_detMgr->getSectorEnvelope(recoSeg->chamberIndex(),
120 recoSeg->sector(),
121 recoSeg->etaIndex());
122 newObj.truthSegment = truthSeg;
123 assoc_itr = allAssociations.end() -1;
124 }
125 ObjectMatching& assocObj{*assoc_itr};
126 assocObj.matchedSegments.push_back(segment);
127 if (!truthSeg) {
128 assocObj.matchedSeeds.push_back(segment->parent());
129 }
130 assocObj.matchedSeedFoundSegment.push_back(1);
131 usedSeeds.insert(segment->parent());
132 }
133
134 if (truthSegments) {
135 for (ObjectMatching& assocObj : allAssociations) {
136 if (!assocObj.truthSegment) {
137 continue;
138 }
139 std::ranges::sort(assocObj.matchedSegments,
140 [&](const Segment* a, const Segment* b){
141 return countOnSameSide(gctx,*assocObj.truthSegment, *a) >
142 countOnSameSide(gctx,*assocObj.truthSegment, *b);
143 });
144 std::ranges::transform(assocObj.matchedSegments, std::back_inserter(assocObj.matchedSeeds),
146 }
147
148 }
150 for (const SegmentSeed* seed : seedContainer) {
152 if (usedSeeds.count(seed)) {
153 continue;
154 }
156 std::vector<std::pair<const xAOD::MuonSegment*, std::size_t>> segCounts{};
157 std::unordered_set<const xAOD::MuonSimHit* > matchedHits = getMatchingSimHits(*seed);
158 for (const xAOD::MuonSimHit* hit : matchedHits) {
159 const xAOD::MuonSegment* truthSeg = getMatchedTruthSegment(*hit);
160 if (!truthSeg) {
161 continue;
162 }
163 auto count_itr = std::ranges::find_if(segCounts, [truthSeg](const auto& segCounter){
164 return segCounter.first == truthSeg;
165 });
166 if (count_itr != segCounts.end()) {
167 ++(count_itr->second);
168 } else {
169 segCounts.emplace_back(std::make_pair(truthSeg, 1ul));
170 }
171 }
173 std::ranges::sort(segCounts, [](const auto& a, const auto& b){
174 return a.second > b.second;
175 });
176 // Add a criterion on the number of counts?
177 const xAOD::MuonSegment* truthSeg = segCounts.size()
178 ? segCounts.front().first : nullptr;
179 if (truthSeg) {
180 auto assoc_itr = std::ranges::find_if(allAssociations,
181 [truthSeg](const ObjectMatching& obj){
182 return obj.truthSegment == truthSeg;
183 });
184 if (assoc_itr == allAssociations.end()) {
185 ObjectMatching & newObj = allAssociations.emplace_back();
186 newObj.chamber = seed->msSector();
187 newObj.truthSegment = truthSeg;
188 assoc_itr = allAssociations.end() -1;
189 }
190 assoc_itr->matchedSeeds.push_back(seed);
191 assoc_itr->matchedSeedFoundSegment.push_back(0);
192 } else {
193 ObjectMatching & newObj = allAssociations.emplace_back();
194 newObj.chamber = seed->msSector();
195 newObj.matchedSeeds.push_back(seed);
196 newObj.matchedSeedFoundSegment.push_back(0);
197 }
198 }
200 if (truthSegments) {
201 for (const xAOD::MuonSegment* truthSeg : *truthSegments) {
203 if (std::ranges::any_of(allAssociations, [truthSeg](const auto& assocObj){
204 return assocObj.truthSegment == truthSeg;
205 })) {
206 continue;
207 }
208 ObjectMatching & newObj = allAssociations.emplace_back();
209 newObj.chamber = m_detMgr->getSectorEnvelope(truthSeg->chamberIndex(),
210 truthSeg->sector(),
211 truthSeg->etaIndex());
212 newObj.truthSegment = truthSeg;
213 }
214 }
215 return allAssociations;
216 }
static Double_t a
const SegmentSeed * parent() const
Returns the seed out of which the segment was built.
const xAOD::MuonSegment * getMatchedTruthSegment(const xAOD::MuonSegment &segment)
Returns the truth-matched segment.
const Segment * detailedSegment(const xAOD::MuonSegment &seg)
Helper function to navigate from the xAOD::MuonSegment to the MuonR4::Segment.
MuonSegment_v1 MuonSegment
Reference the current persistent version:

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

578{
579 size_t pos = 0;
580 while((pos = str.find(oldStr, pos)) != std::string::npos)
581 {
582 str.replace(pos, oldStr.length(), newStr);
583 pos += newStr.length();
584 }
585}

◆ 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
93 ATH_CHECK ( histSvc().retrieve() );
94
95 // Configure the underlying AthHistogramming helper
100
101 // Print some setup information into the log file
102 ATH_MSG_DEBUG ("Initializing " << name() << "...");
103 ATH_MSG_DEBUG (" using THistService = " << m_histSvc );
104 ATH_MSG_DEBUG (" using RootStreamName = " << m_prefix );
105 ATH_MSG_DEBUG (" using RootDirName = " << m_rootDir );
106 ATH_MSG_DEBUG (" using HistNamePrefix = " << m_histNamePrefix );
107 ATH_MSG_DEBUG (" using HistNamePostfix = " << m_histNamePostfix );
108 ATH_MSG_DEBUG (" using HistTitlePrefix = " << m_histTitlePrefix );
109 ATH_MSG_DEBUG (" using HistTitlePostfix = " << m_histTitlePostfix );
110
111
112 // re-direct to base class...
114}
virtual StatusCode sysInitialize() override
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.

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

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

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

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

Definition at line 289 of file MuonHoughTransformTester.h.

289{nullptr};

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

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

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

Tracking geometry context.

Definition at line 113 of file MuonHoughTransformTester.h.

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

◆ m_graphMap

GraphMap_t AthHistogramming::m_graphMap
privateinherited

The map of TGraph names to their pointers.

Definition at line 224 of file AthHistogramming.h.

◆ m_histMap

HistMap_t AthHistogramming::m_histMap
privateinherited

The map of histogram names to their pointers.

Definition at line 203 of file AthHistogramming.h.

◆ m_histNamePostfix

std::string AthHistogramAlgorithm::m_histNamePostfix
privateinherited

The postfix for the histogram THx name.

Definition at line 97 of file AthHistogramAlgorithm.h.

◆ m_histNamePrefix

std::string AthHistogramAlgorithm::m_histNamePrefix
privateinherited

The prefix for the histogram THx name.

Definition at line 94 of file AthHistogramAlgorithm.h.

◆ m_histSvc

ServiceHandle<ITHistSvc> AthHistogramAlgorithm::m_histSvc
privateinherited

Default constructor: AthHistogramAlgorithm();.

a handle on the Hist/TTree registration service

Definition at line 83 of file AthHistogramAlgorithm.h.

◆ m_histTitlePostfix

std::string AthHistogramAlgorithm::m_histTitlePostfix
privateinherited

The postfix for the histogram THx title.

Definition at line 103 of file AthHistogramAlgorithm.h.

◆ m_histTitlePrefix

std::string AthHistogramAlgorithm::m_histTitlePrefix
privateinherited

The prefix for the histogram THx title.

Definition at line 100 of file AthHistogramAlgorithm.h.

◆ m_idHelperSvc

ServiceHandle<Muon::IMuonIdHelperSvc> MuonValR4::MuonHoughTransformTester::m_idHelperSvc {this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}
private

Definition at line 115 of file MuonHoughTransformTester.h.

115{this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"};

◆ m_isMC

Gaudi::Property<bool> MuonValR4::MuonHoughTransformTester::m_isMC {this, "isMC", false, "Toggle whether the job is ran on MC or not"}
private

Definition at line 118 of file MuonHoughTransformTester.h.

118{this, "isMC", false, "Toggle whether the job is ran on MC or not"};

◆ m_msg

MsgStream AthHistogramming::m_msg
privateinherited

Cached Message Stream.

Definition at line 250 of file AthHistogramming.h.

◆ m_name

std::string AthHistogramming::m_name
privateinherited

Instance name.

Definition at line 247 of file AthHistogramming.h.

◆ m_out_bucketEnd

MuonVal::ScalarBranch<float>& MuonValR4::MuonHoughTransformTester::m_out_bucketEnd {m_tree.newScalar<float>("bucketEnd", -1)}
private

Definition at line 133 of file MuonHoughTransformTester.h.

133{m_tree.newScalar<float>("bucketEnd", -1)};

◆ m_out_bucketHitGap

MuonVal::ScalarBranch<float>& MuonValR4::MuonHoughTransformTester::m_out_bucketHitGap {m_tree.newScalar<float>("bucketHiGap", 0.)}
private

Definition at line 134 of file MuonHoughTransformTester.h.

134{m_tree.newScalar<float>("bucketHiGap", 0.)};

◆ m_out_bucketStart

MuonVal::ScalarBranch<float>& MuonValR4::MuonHoughTransformTester::m_out_bucketStart {m_tree.newScalar<float>("bucketStart", 1)}
private

Definition at line 132 of file MuonHoughTransformTester.h.

132{m_tree.newScalar<float>("bucketStart", 1)};

◆ m_out_chamberIndex

MuonVal::ScalarBranch<int>& MuonValR4::MuonHoughTransformTester::m_out_chamberIndex {m_tree.newScalar<int>("chamberIndex")}
private

====== Common block: Filled for all entries ===========

chamber index field

Definition at line 126 of file MuonHoughTransformTester.h.

126{m_tree.newScalar<int>("chamberIndex")};

◆ m_out_gen_Eta

MuonVal::ScalarBranch<float>& MuonValR4::MuonHoughTransformTester::m_out_gen_Eta {m_tree.newScalar<float>("genEta",-10.)}
private

global particle properties

Definition at line 163 of file MuonHoughTransformTester.h.

163{m_tree.newScalar<float>("genEta",-10.)};

◆ m_out_gen_maxYhit

MuonVal::ScalarBranch<float>& MuonValR4::MuonHoughTransformTester::m_out_gen_maxYhit {m_tree.newScalar<float>("genMaxYhit", -1.0)}
private

Definition at line 188 of file MuonHoughTransformTester.h.

188{m_tree.newScalar<float>("genMaxYhit", -1.0)};

◆ m_out_gen_minYhit

MuonVal::ScalarBranch<float>& MuonValR4::MuonHoughTransformTester::m_out_gen_minYhit {m_tree.newScalar<float>("genMinYhit", 1.0)}
private

Definition at line 187 of file MuonHoughTransformTester.h.

187{m_tree.newScalar<float>("genMinYhit", 1.0)};

◆ m_out_gen_nHits

MuonVal::ScalarBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_gen_nHits {m_tree.newScalar<unsigned short>("genNHits",0)}
private

Truth - hit count summary.

Definition at line 178 of file MuonHoughTransformTester.h.

178{m_tree.newScalar<unsigned short>("genNHits",0)};

◆ m_out_gen_nMDTHits

MuonVal::ScalarBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_gen_nMDTHits {m_tree.newScalar<unsigned short>("genNMdtHits",0)}
private

Definition at line 180 of file MuonHoughTransformTester.h.

180{m_tree.newScalar<unsigned short>("genNMdtHits",0)};

◆ m_out_gen_nMmHits

MuonVal::ScalarBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_gen_nMmHits {m_tree.newScalar<unsigned short>("genNMmHits",0)}
private

Definition at line 183 of file MuonHoughTransformTester.h.

183{m_tree.newScalar<unsigned short>("genNMmHits",0)};

◆ m_out_gen_nNswHits

MuonVal::ScalarBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_gen_nNswHits {m_tree.newScalar<unsigned short>("genNNswHits",0)}
private

Definition at line 182 of file MuonHoughTransformTester.h.

182{m_tree.newScalar<unsigned short>("genNNswHits",0)};

◆ m_out_gen_nRPCHits

MuonVal::ScalarBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_gen_nRPCHits {m_tree.newScalar<unsigned short>("genNRpcHits",0)}
private

Definition at line 179 of file MuonHoughTransformTester.h.

179{m_tree.newScalar<unsigned short>("genNRpcHits",0)};

◆ m_out_gen_nSTGCHits

MuonVal::ScalarBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_gen_nSTGCHits {m_tree.newScalar<unsigned short>("genNsTgcHits",0)}
private

Definition at line 184 of file MuonHoughTransformTester.h.

184{m_tree.newScalar<unsigned short>("genNsTgcHits",0)};

◆ m_out_gen_nTGCHits

MuonVal::ScalarBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_gen_nTGCHits {m_tree.newScalar<unsigned short>("genNTgcHits",0)}
private

Definition at line 181 of file MuonHoughTransformTester.h.

181{m_tree.newScalar<unsigned short>("genNTgcHits",0)};

◆ m_out_gen_Phi

MuonVal::ScalarBranch<float>& MuonValR4::MuonHoughTransformTester::m_out_gen_Phi {m_tree.newScalar<float>("genPhi",-10.)}
private

Definition at line 164 of file MuonHoughTransformTester.h.

164{m_tree.newScalar<float>("genPhi",-10.)};

◆ m_out_gen_Pt

MuonVal::ScalarBranch<float>& MuonValR4::MuonHoughTransformTester::m_out_gen_Pt {m_tree.newScalar<float>("genPt",-10.)}
private

Definition at line 165 of file MuonHoughTransformTester.h.

165{m_tree.newScalar<float>("genPt",-10.)};

◆ m_out_gen_Q

MuonVal::ScalarBranch<short>& MuonValR4::MuonHoughTransformTester::m_out_gen_Q {m_tree.newScalar<short>("genQ", 0)}
private

Definition at line 166 of file MuonHoughTransformTester.h.

166{m_tree.newScalar<short>("genQ", 0)};

◆ m_out_gen_tanphi

MuonVal::ScalarBranch<float>& MuonValR4::MuonHoughTransformTester::m_out_gen_tanphi {m_tree.newScalar<float>("genTanPhi", 0.0)}
private

Definition at line 170 of file MuonHoughTransformTester.h.

170{m_tree.newScalar<float>("genTanPhi", 0.0)};

◆ m_out_gen_tantheta

MuonVal::ScalarBranch<float>& MuonValR4::MuonHoughTransformTester::m_out_gen_tantheta {m_tree.newScalar<float>("genTanTheta", 0.0)}
private

Definition at line 169 of file MuonHoughTransformTester.h.

169{m_tree.newScalar<float>("genTanTheta", 0.0)};

◆ m_out_gen_time

MuonVal::ScalarBranch<float>& MuonValR4::MuonHoughTransformTester::m_out_gen_time {m_tree.newScalar<float>("genTime", 0.0)}
private

Definition at line 172 of file MuonHoughTransformTester.h.

172{m_tree.newScalar<float>("genTime", 0.0)};

◆ m_out_gen_truthOrigin

MuonVal::ScalarBranch<int>& MuonValR4::MuonHoughTransformTester::m_out_gen_truthOrigin {m_tree.newScalar<int>("genTruthOrigin", -1)}
private

Definition at line 174 of file MuonHoughTransformTester.h.

174{m_tree.newScalar<int>("genTruthOrigin", -1)};

◆ m_out_gen_truthType

MuonVal::ScalarBranch<int>& MuonValR4::MuonHoughTransformTester::m_out_gen_truthType {m_tree.newScalar<int>("genTruthType", -1)}
private

Definition at line 175 of file MuonHoughTransformTester.h.

175{m_tree.newScalar<int>("genTruthType", -1)};

◆ m_out_gen_x0

MuonVal::ScalarBranch<float>& MuonValR4::MuonHoughTransformTester::m_out_gen_x0 {m_tree.newScalar<float>("genX0", 0.0)}
private

Definition at line 171 of file MuonHoughTransformTester.h.

171{m_tree.newScalar<float>("genX0", 0.0)};

◆ m_out_gen_y0

MuonVal::ScalarBranch<float>& MuonValR4::MuonHoughTransformTester::m_out_gen_y0 {m_tree.newScalar<float>("genY0", 0.0)}
private

Truth - segment parameters.

Definition at line 168 of file MuonHoughTransformTester.h.

168{m_tree.newScalar<float>("genY0", 0.0)};

◆ m_out_hasTruth

MuonVal::ScalarBranch<bool>& MuonValR4::MuonHoughTransformTester::m_out_hasTruth {m_tree.newScalar<bool>("hasTruth",false)}
private

======= Truth block: Filled if we have a truth match. ============

existence of a truth match

Definition at line 160 of file MuonHoughTransformTester.h.

160{m_tree.newScalar<bool>("hasTruth",false)};

◆ m_out_nPhiSpacePoints

MuonVal::ScalarBranch<unsigned char>& MuonValR4::MuonHoughTransformTester::m_out_nPhiSpacePoints {m_tree.newScalar<unsigned char>("bucketNPhiMeass",0)}
private

Number of phi hits in the bucket.

Definition at line 148 of file MuonHoughTransformTester.h.

148{m_tree.newScalar<unsigned char>("bucketNPhiMeass",0)};

◆ m_out_nPrecSpacePoints

MuonVal::ScalarBranch<unsigned char>& MuonValR4::MuonHoughTransformTester::m_out_nPrecSpacePoints {m_tree.newScalar<unsigned char>("bucketNPrecMeas",0)}
private

Number of precision hits in the bucket.

Definition at line 146 of file MuonHoughTransformTester.h.

146{m_tree.newScalar<unsigned char>("bucketNPrecMeas",0)};

◆ m_out_nSpacePoints

MuonVal::ScalarBranch<unsigned char>& MuonValR4::MuonHoughTransformTester::m_out_nSpacePoints {m_tree.newScalar<unsigned char>("bucketNHits",0)}
private

Number of all space points in the bucket.

Definition at line 144 of file MuonHoughTransformTester.h.

144{m_tree.newScalar<unsigned char>("bucketNHits",0)};

◆ m_out_nTruePhiSpacePoints

MuonVal::ScalarBranch<unsigned char>& MuonValR4::MuonHoughTransformTester::m_out_nTruePhiSpacePoints {m_tree.newScalar<unsigned char>("bucketNTruePhiMeass",0)}
private

Number of phi hits in the bucket.

Definition at line 155 of file MuonHoughTransformTester.h.

155{m_tree.newScalar<unsigned char>("bucketNTruePhiMeass",0)};

◆ m_out_nTruePrecSpacePoints

MuonVal::ScalarBranch<unsigned char>& MuonValR4::MuonHoughTransformTester::m_out_nTruePrecSpacePoints {m_tree.newScalar<unsigned char>("bucketNTruePrecMeas",0)}
private

Number of precision hits in the bucket.

Definition at line 153 of file MuonHoughTransformTester.h.

153{m_tree.newScalar<unsigned char>("bucketNTruePrecMeas",0)};

◆ m_out_nTrueSpacePoints

MuonVal::ScalarBranch<unsigned char>& MuonValR4::MuonHoughTransformTester::m_out_nTrueSpacePoints {m_tree.newScalar<unsigned char>("bucketNTrueMeas",0)}
private

Number of all space points in the bucket.

Definition at line 151 of file MuonHoughTransformTester.h.

151{m_tree.newScalar<unsigned char>("bucketNTrueMeas",0)};

◆ m_out_seed_hasPhiExtension

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_seed_hasPhiExtension {m_tree.newVector<unsigned short>("seedHasPhiExtension", false)}
private

Definition at line 196 of file MuonHoughTransformTester.h.

196{m_tree.newVector<unsigned short>("seedHasPhiExtension", false)};

◆ m_out_seed_ledToSegment

MuonVal::VectorBranch<unsigned char>& MuonValR4::MuonHoughTransformTester::m_out_seed_ledToSegment {m_tree.newVector<unsigned char>("seedConvertedToSegment",false)}
private

Definition at line 232 of file MuonHoughTransformTester.h.

232{m_tree.newVector<unsigned char>("seedConvertedToSegment",false)};

◆ m_out_seed_maxYhit

MuonVal::VectorBranch<float>& MuonValR4::MuonHoughTransformTester::m_out_seed_maxYhit {m_tree.newVector<float>("seedMaxYhit", -1.0)}
private

Definition at line 208 of file MuonHoughTransformTester.h.

208{m_tree.newVector<float>("seedMaxYhit", -1.0)};

◆ m_out_seed_minYhit

MuonVal::VectorBranch<float>& MuonValR4::MuonHoughTransformTester::m_out_seed_minYhit {m_tree.newVector<float>("seedMinYhit", 1.0)}
private

Definition at line 207 of file MuonHoughTransformTester.h.

207{m_tree.newVector<float>("seedMinYhit", 1.0)};

◆ m_out_seed_n

MuonVal::ScalarBranch<unsigned>& MuonValR4::MuonHoughTransformTester::m_out_seed_n {m_tree.newScalar<unsigned>("nSeeds", 0)}
private

========== Seed block: Filled when we have one or multiple seeds ============= seed count

Definition at line 192 of file MuonHoughTransformTester.h.

192{m_tree.newScalar<unsigned>("nSeeds", 0)};

◆ m_out_seed_nEtaHits

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_seed_nEtaHits {m_tree.newVector<unsigned short>("seedNEtaHits", 0)}
private

Definition at line 212 of file MuonHoughTransformTester.h.

212{m_tree.newVector<unsigned short>("seedNEtaHits", 0)};

◆ m_out_seed_nHits

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_seed_nHits {m_tree.newVector<unsigned short>("seedNHits", 0)}
private

Definition at line 210 of file MuonHoughTransformTester.h.

210{m_tree.newVector<unsigned short>("seedNHits", 0)};

◆ m_out_seed_nMatchedHits

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_seed_nMatchedHits {m_tree.newVector<unsigned short>("seedTruthMatchedHits", false)}
private

Definition at line 198 of file MuonHoughTransformTester.h.

198{m_tree.newVector<unsigned short>("seedTruthMatchedHits", false)};

◆ m_out_seed_nMdt

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_seed_nMdt {m_tree.newVector<unsigned short>("seedNMdtHits", 0)}
private

Definition at line 215 of file MuonHoughTransformTester.h.

215{m_tree.newVector<unsigned short>("seedNMdtHits", 0)};

◆ m_out_seed_nMmEta

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_seed_nMmEta {m_tree.newVector<unsigned short>("seedNMmHitsEta", 0)}
private

Definition at line 222 of file MuonHoughTransformTester.h.

222{m_tree.newVector<unsigned short>("seedNMmHitsEta", 0)};

◆ m_out_seed_nMmStereo

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_seed_nMmStereo {m_tree.newVector<unsigned short>("seedNMmHitsStereo", 0)}
private

Definition at line 223 of file MuonHoughTransformTester.h.

223{m_tree.newVector<unsigned short>("seedNMmHitsStereo", 0)};

◆ m_out_seed_nPhiHits

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_seed_nPhiHits {m_tree.newVector<unsigned short>("seedNPhiHits", 0)}
private

Definition at line 213 of file MuonHoughTransformTester.h.

213{m_tree.newVector<unsigned short>("seedNPhiHits", 0)};

◆ m_out_seed_nPrecHits

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_seed_nPrecHits {m_tree.newVector<unsigned short>("seedNPrecHits", 0)}
private

Definition at line 211 of file MuonHoughTransformTester.h.

211{m_tree.newVector<unsigned short>("seedNPrecHits", 0)};

◆ m_out_seed_nRpc

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_seed_nRpc {m_tree.newVector<unsigned short>("seedNRpcHits", 0)}
private

Definition at line 216 of file MuonHoughTransformTester.h.

216{m_tree.newVector<unsigned short>("seedNRpcHits", 0)};

◆ m_out_seed_nsTgcPad

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_seed_nsTgcPad {m_tree.newVector<unsigned short>("seedNsTgcPadHits", 0)}
private

Definition at line 220 of file MuonHoughTransformTester.h.

220{m_tree.newVector<unsigned short>("seedNsTgcPadHits", 0)};

◆ m_out_seed_nsTgcStrip

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_seed_nsTgcStrip {m_tree.newVector<unsigned short>("seedNsTgcStripHits", 0)}
private

Definition at line 218 of file MuonHoughTransformTester.h.

218{m_tree.newVector<unsigned short>("seedNsTgcStripHits", 0)};

◆ m_out_seed_nsTgcWire

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_seed_nsTgcWire {m_tree.newVector<unsigned short>("seedNsTgcWireHits", 0)}
private

Definition at line 219 of file MuonHoughTransformTester.h.

219{m_tree.newVector<unsigned short>("seedNsTgcWireHits", 0)};

◆ m_out_seed_nTgc

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_seed_nTgc {m_tree.newVector<unsigned short>("seedNTgcHits", 0)}
private

Definition at line 217 of file MuonHoughTransformTester.h.

217{m_tree.newVector<unsigned short>("seedNTgcHits", 0)};

◆ m_out_seed_nTrueEtaHits

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_seed_nTrueEtaHits {m_tree.newVector<unsigned short>("seedNTrueEtaHits", 0)}
private

Definition at line 230 of file MuonHoughTransformTester.h.

230{m_tree.newVector<unsigned short>("seedNTrueEtaHits", 0)};

◆ m_out_seed_nTrueHits

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_seed_nTrueHits {m_tree.newVector<unsigned short>("seedNTrueHits", 0)}
private

Labelled hits from the pattern visualization tool.

Definition at line 227 of file MuonHoughTransformTester.h.

227{m_tree.newVector<unsigned short>("seedNTrueHits", 0)};

◆ m_out_seed_nTruePhiHits

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_seed_nTruePhiHits {m_tree.newVector<unsigned short>("seedNTruePhiHits", 0)}
private

Definition at line 229 of file MuonHoughTransformTester.h.

229{m_tree.newVector<unsigned short>("seedNTruePhiHits", 0)};

◆ m_out_seed_nTruePrecHits

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_seed_nTruePrecHits {m_tree.newVector<unsigned short>("seedNTruePrecHits", 0)}
private

Definition at line 228 of file MuonHoughTransformTester.h.

228{m_tree.newVector<unsigned short>("seedNTruePrecHits", 0)};

◆ m_out_seed_tanphi

MuonVal::VectorBranch<float>& MuonValR4::MuonHoughTransformTester::m_out_seed_tanphi {m_tree.newVector<float>("seedTanPhi", 0.0)}
private

Definition at line 204 of file MuonHoughTransformTester.h.

204{m_tree.newVector<float>("seedTanPhi", 0.0)};

◆ m_out_seed_tantheta

MuonVal::VectorBranch<float>& MuonValR4::MuonHoughTransformTester::m_out_seed_tantheta {m_tree.newVector<float>("seedTanTheta", 0.0)}
private

Definition at line 203 of file MuonHoughTransformTester.h.

203{m_tree.newVector<float>("seedTanTheta", 0.0)};

◆ m_out_seed_x0

MuonVal::VectorBranch<float>& MuonValR4::MuonHoughTransformTester::m_out_seed_x0 {m_tree.newVector<float>("seedX0", 0.0)}
private

Definition at line 202 of file MuonHoughTransformTester.h.

202{m_tree.newVector<float>("seedX0", 0.0)};

◆ m_out_seed_y0

MuonVal::VectorBranch<float>& MuonValR4::MuonHoughTransformTester::m_out_seed_y0 {m_tree.newVector<float>("seedY0", 0.0)}
private

Definition at line 201 of file MuonHoughTransformTester.h.

201{m_tree.newVector<float>("seedY0", 0.0)};

◆ m_out_segment_chi2

MuonVal::VectorBranch<float>& MuonValR4::MuonHoughTransformTester::m_out_segment_chi2 {m_tree.newVector<float>("segmentChi2", -1.)}
private

Definition at line 241 of file MuonHoughTransformTester.h.

241{m_tree.newVector<float>("segmentChi2", -1.)};

◆ m_out_segment_err_tanphi

MuonVal::VectorBranch<float>& MuonValR4::MuonHoughTransformTester::m_out_segment_err_tanphi {m_tree.newVector<float>("segmentErrTanPhi", -1.0)}
private

Definition at line 260 of file MuonHoughTransformTester.h.

260{m_tree.newVector<float>("segmentErrTanPhi", -1.0)};

◆ m_out_segment_err_tantheta

MuonVal::VectorBranch<float>& MuonValR4::MuonHoughTransformTester::m_out_segment_err_tantheta {m_tree.newVector<float>("segmentErrTanTheta", -1.0)}
private

Definition at line 259 of file MuonHoughTransformTester.h.

259{m_tree.newVector<float>("segmentErrTanTheta", -1.0)};

◆ m_out_segment_err_time

MuonVal::VectorBranch<float>& MuonValR4::MuonHoughTransformTester::m_out_segment_err_time {m_tree.newVector<float>("segmentErrTime", -1.0)}
private

Definition at line 261 of file MuonHoughTransformTester.h.

261{m_tree.newVector<float>("segmentErrTime", -1.0)};

◆ m_out_segment_err_x0

MuonVal::VectorBranch<float>& MuonValR4::MuonHoughTransformTester::m_out_segment_err_x0 {m_tree.newVector<float>("segmentErrX0", -1.0)}
private

Definition at line 258 of file MuonHoughTransformTester.h.

258{m_tree.newVector<float>("segmentErrX0", -1.0)};

◆ m_out_segment_err_y0

MuonVal::VectorBranch<float>& MuonValR4::MuonHoughTransformTester::m_out_segment_err_y0 {m_tree.newVector<float>("segmentErrY0", -1.0)}
private

Definition at line 257 of file MuonHoughTransformTester.h.

257{m_tree.newVector<float>("segmentErrY0", -1.0)};

◆ m_out_segment_fitIter

MuonVal::VectorBranch<uint16_t>& MuonValR4::MuonHoughTransformTester::m_out_segment_fitIter {m_tree.newVector<uint16_t>("segmentFitIterations", 0)}
private

Definition at line 244 of file MuonHoughTransformTester.h.

244{m_tree.newVector<uint16_t>("segmentFitIterations", 0)};
setWord1 uint16_t

◆ m_out_segment_hasPhi

MuonVal::VectorBranch<bool>& MuonValR4::MuonHoughTransformTester::m_out_segment_hasPhi {m_tree.newVector<bool>("segmentHasPhiHits", false)}
private

Definition at line 247 of file MuonHoughTransformTester.h.

247{m_tree.newVector<bool>("segmentHasPhiHits", false)};

◆ m_out_segment_hasTimeFit

MuonVal::VectorBranch<bool>& MuonValR4::MuonHoughTransformTester::m_out_segment_hasTimeFit {m_tree.newVector<bool>("segmentHasTimeFit", false)}
private

Definition at line 243 of file MuonHoughTransformTester.h.

243{m_tree.newVector<bool>("segmentHasTimeFit", false)};

◆ m_out_segment_maxYhit

MuonVal::VectorBranch<float>& MuonValR4::MuonHoughTransformTester::m_out_segment_maxYhit {m_tree.newVector<float>("segmentMaxYhit", -1.0)}
private

Definition at line 284 of file MuonHoughTransformTester.h.

284{m_tree.newVector<float>("segmentMaxYhit", -1.0)};

◆ m_out_segment_minYhit

MuonVal::VectorBranch<float>& MuonValR4::MuonHoughTransformTester::m_out_segment_minYhit {m_tree.newVector<float>("segmentMinYhit", 1.0)}
private

Definition at line 283 of file MuonHoughTransformTester.h.

283{m_tree.newVector<float>("segmentMinYhit", 1.0)};

◆ m_out_segment_n

MuonVal::ScalarBranch<unsigned>& MuonValR4::MuonHoughTransformTester::m_out_segment_n {m_tree.newScalar<unsigned>("nSegments", 0)}
private

========== Segment block: Filled when we have one or multiple segments =============

Definition at line 237 of file MuonHoughTransformTester.h.

237{m_tree.newScalar<unsigned>("nSegments", 0)};

◆ m_out_segment_nDoF

MuonVal::VectorBranch<uint16_t>& MuonValR4::MuonHoughTransformTester::m_out_segment_nDoF {m_tree.newVector<uint16_t>("segmentNdoF", 0)}
private

Definition at line 242 of file MuonHoughTransformTester.h.

242{m_tree.newVector<uint16_t>("segmentNdoF", 0)};

◆ m_out_segment_nMdtHits

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_segment_nMdtHits {m_tree.newVector<unsigned short>("segmentNMdtHits", 0)}
private

Definition at line 265 of file MuonHoughTransformTester.h.

265{m_tree.newVector<unsigned short>("segmentNMdtHits", 0)};

◆ m_out_segment_nMmEtaHits

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_segment_nMmEtaHits {m_tree.newVector<unsigned short>("segmentNMmEtaHits", 0)}
private

Definition at line 271 of file MuonHoughTransformTester.h.

271{m_tree.newVector<unsigned short>("segmentNMmEtaHits", 0)};

◆ m_out_segment_nMmStereoHits

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_segment_nMmStereoHits {m_tree.newVector<unsigned short>("segmentNMmStereoHits", 0)}
private

Definition at line 272 of file MuonHoughTransformTester.h.

272{m_tree.newVector<unsigned short>("segmentNMmStereoHits", 0)};

◆ m_out_segment_nRpcEtaHits

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_segment_nRpcEtaHits {m_tree.newVector<unsigned short>("segmentNRpcEtaHits", 0)}
private

Definition at line 266 of file MuonHoughTransformTester.h.

266{m_tree.newVector<unsigned short>("segmentNRpcEtaHits", 0)};

◆ m_out_segment_nRpcPhiHits

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_segment_nRpcPhiHits {m_tree.newVector<unsigned short>("segmentNRpcPhiHits", 0)}
private

Definition at line 267 of file MuonHoughTransformTester.h.

267{m_tree.newVector<unsigned short>("segmentNRpcPhiHits", 0)};

◆ m_out_segment_nsTgcPadpHits

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_segment_nsTgcPadpHits {m_tree.newVector<unsigned short>("segmentNsTgcPadHits", 0)}
private

Definition at line 275 of file MuonHoughTransformTester.h.

275{m_tree.newVector<unsigned short>("segmentNsTgcPadHits", 0)};

◆ m_out_segment_nsTgcStripHits

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_segment_nsTgcStripHits {m_tree.newVector<unsigned short>("segmentNsTgcStripHits", 0)}
private

Definition at line 273 of file MuonHoughTransformTester.h.

273{m_tree.newVector<unsigned short>("segmentNsTgcStripHits", 0)};

◆ m_out_segment_nsTgcWireHits

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_segment_nsTgcWireHits {m_tree.newVector<unsigned short>("segmentNsTgcWireHits", 0)}
private

Definition at line 274 of file MuonHoughTransformTester.h.

274{m_tree.newVector<unsigned short>("segmentNsTgcWireHits", 0)};

◆ m_out_segment_nTgcEtaHits

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_segment_nTgcEtaHits {m_tree.newVector<unsigned short>("segmentNTgcEtaHits", 0)}
private

Definition at line 268 of file MuonHoughTransformTester.h.

268{m_tree.newVector<unsigned short>("segmentNTgcEtaHits", 0)};

◆ m_out_segment_nTgcPhiHits

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_segment_nTgcPhiHits {m_tree.newVector<unsigned short>("segmentNTgcPhiHits", 0)}
private

Definition at line 269 of file MuonHoughTransformTester.h.

269{m_tree.newVector<unsigned short>("segmentNTgcPhiHits", 0)};

◆ m_out_segment_nTrueEtaHits

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_segment_nTrueEtaHits {m_tree.newVector<unsigned short>("segmentNTrueEtaHits", 0)}
private

Definition at line 280 of file MuonHoughTransformTester.h.

280{m_tree.newVector<unsigned short>("segmentNTrueEtaHits", 0)};

◆ m_out_segment_nTrueHits

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_segment_nTrueHits {m_tree.newVector<unsigned short>("segmentNTrueHits", 0)}
private

Labelled hits from the pattern visualization tool.

Definition at line 277 of file MuonHoughTransformTester.h.

277{m_tree.newVector<unsigned short>("segmentNTrueHits", 0)};

◆ m_out_segment_nTruePhiHits

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_segment_nTruePhiHits {m_tree.newVector<unsigned short>("segmentNTruePhiHits", 0)}
private

Definition at line 279 of file MuonHoughTransformTester.h.

279{m_tree.newVector<unsigned short>("segmentNTruePhiHits", 0)};

◆ m_out_segment_nTruePrecHits

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_segment_nTruePrecHits {m_tree.newVector<unsigned short>("segmentNTruePrecHits", 0)}
private

Definition at line 278 of file MuonHoughTransformTester.h.

278{m_tree.newVector<unsigned short>("segmentNTruePrecHits", 0)};

◆ m_out_segment_tanphi

MuonVal::VectorBranch<float>& MuonValR4::MuonHoughTransformTester::m_out_segment_tanphi {m_tree.newVector<float>("segmentTanPhi", 0.0)}
private

Definition at line 253 of file MuonHoughTransformTester.h.

253{m_tree.newVector<float>("segmentTanPhi", 0.0)};

◆ m_out_segment_tantheta

MuonVal::VectorBranch<float>& MuonValR4::MuonHoughTransformTester::m_out_segment_tantheta {m_tree.newVector<float>("segmentTanTheta", 0.0)}
private

Definition at line 252 of file MuonHoughTransformTester.h.

252{m_tree.newVector<float>("segmentTanTheta", 0.0)};

◆ m_out_segment_time

MuonVal::VectorBranch<float>& MuonValR4::MuonHoughTransformTester::m_out_segment_time {m_tree.newVector<float>("segmentTime", 0.)}
private

Definition at line 254 of file MuonHoughTransformTester.h.

254{m_tree.newVector<float>("segmentTime", 0.)};

◆ m_out_segment_truthMatchedHits

MuonVal::VectorBranch<unsigned short>& MuonValR4::MuonHoughTransformTester::m_out_segment_truthMatchedHits {m_tree.newVector<unsigned short>("segmentTruthMatchedHits", 0)}
private

Definition at line 264 of file MuonHoughTransformTester.h.

264{m_tree.newVector<unsigned short>("segmentTruthMatchedHits", 0)};

◆ m_out_segment_x0

MuonVal::VectorBranch<float>& MuonValR4::MuonHoughTransformTester::m_out_segment_x0 {m_tree.newVector<float>("segmentX0", 0.0)}
private

Definition at line 251 of file MuonHoughTransformTester.h.

251{m_tree.newVector<float>("segmentX0", 0.0)};

◆ m_out_segment_y0

MuonVal::VectorBranch<float>& MuonValR4::MuonHoughTransformTester::m_out_segment_y0 {m_tree.newVector<float>("segmentY0", 0.0)}
private

Definition at line 250 of file MuonHoughTransformTester.h.

250{m_tree.newVector<float>("segmentY0", 0.0)};

◆ m_out_stationPhi

MuonVal::ScalarBranch<int>& MuonValR4::MuonHoughTransformTester::m_out_stationPhi {m_tree.newScalar<int>("stationPhi")}
private

phi index of the station

Definition at line 130 of file MuonHoughTransformTester.h.

130{m_tree.newScalar<int>("stationPhi")};

◆ m_out_stationSide

MuonVal::ScalarBranch<short>& MuonValR4::MuonHoughTransformTester::m_out_stationSide {m_tree.newScalar<short>("stationSide")}
private

+1 for A-, -1 of C-side

Definition at line 128 of file MuonHoughTransformTester.h.

128{m_tree.newScalar<short>("stationSide")};

◆ m_patternSeedKeys

SG::ReadHandleKeyArray<MuonR4::SegmentSeedContainer> MuonValR4::MuonHoughTransformTester::m_patternSeedKeys {this, "SegmentSeedKeys", {"MuonHoughStationSegmentSeeds"}}
private

List of the two segment seed containers from which the segments are buiit (Complets the pattern finding step).

Definition at line 109 of file MuonHoughTransformTester.h.

109{this, "SegmentSeedKeys", {"MuonHoughStationSegmentSeeds"}};

◆ m_prefix

std::string AthHistogramAlgorithm::m_prefix
privateinherited

Name of the ROOT output stream (file).

Definition at line 88 of file AthHistogramAlgorithm.h.

◆ m_recoSegKey

SG::ReadHandleKey<xAOD::MuonSegmentContainer> MuonValR4::MuonHoughTransformTester::m_recoSegKey {this, "SegmentKey", "MuonSegmentsFromR4"}
private

Key to the xAOD::MuonSegment container.

Definition at line 105 of file MuonHoughTransformTester.h.

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

◆ m_recoSegLinks

Gaudi::Property<std::vector<std::string> > MuonValR4::MuonHoughTransformTester::m_recoSegLinks {this, "RecoSegLinks", {"truthSegmentLink", "truthParticleLink"}}
private

name of the truth link decorations for the reco segment container

Definition at line 107 of file MuonHoughTransformTester.h.

107{this, "RecoSegLinks", {"truthSegmentLink", "truthParticleLink"}};

◆ m_rootDir

std::string AthHistogramAlgorithm::m_rootDir
privateinherited

Name of the ROOT directory.

Definition at line 91 of file AthHistogramAlgorithm.h.

◆ m_spKeys

SG::ReadHandleKeyArray<MuonR4::SpacePointContainer> MuonValR4::MuonHoughTransformTester::m_spKeys {this, "SpacePointKeys", {"MuonSpacePoints"}}
private

List of the space point containers in the event legacy + NSW containers.

Definition at line 111 of file MuonHoughTransformTester.h.

111{this, "SpacePointKeys", {"MuonSpacePoints"}};

◆ m_spMatchedToPattern

MuonVal::MatrixBranch<unsigned char>& MuonValR4::MuonHoughTransformTester::m_spMatchedToPattern {m_tree.newMatrix<unsigned char>("seedMatchedSp")}
private

Branch indicating which space points in the tree are associated to the i-th pattern.

Definition at line 139 of file MuonHoughTransformTester.h.

139{m_tree.newMatrix<unsigned char>("seedMatchedSp")};

◆ m_spMatchedToSegment

MuonVal::MatrixBranch<unsigned char>& MuonValR4::MuonHoughTransformTester::m_spMatchedToSegment {m_tree.newMatrix<unsigned char>("segmentMatchedSp")}
private

Branch indicating which space points in the tree are associated to the i-th segment.

Definition at line 141 of file MuonHoughTransformTester.h.

141{m_tree.newMatrix<unsigned char>("segmentMatchedSp")};

◆ m_spTester

std::shared_ptr<SpacePointTesterModule> MuonValR4::MuonHoughTransformTester::m_spTester {}
private

Branch dumping all the space points from the difference buckets.

Definition at line 137 of file MuonHoughTransformTester.h.

137{};

◆ m_streamName

std::string AthHistogramming::m_streamName
privateinherited

Name of the ROOT output stream (file).

Definition at line 228 of file AthHistogramming.h.

◆ m_tree

MuonVal::MuonTesterTree MuonValR4::MuonHoughTransformTester::m_tree {"MuonEtaHoughTest","MuonEtaHoughTransformTest"}
private

Definition at line 96 of file MuonHoughTransformTester.h.

96{"MuonEtaHoughTest","MuonEtaHoughTransformTest"};

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

Gaudi::Property<std::vector<std::string> > MuonValR4::MuonHoughTransformTester::m_truthLinks {this, "TruthSegLinks", {"simHitLinks", "truthParticleLink"}}
private

Name of the decorations for the truth segment.

Definition at line 103 of file MuonHoughTransformTester.h.

103{this, "TruthSegLinks", {"simHitLinks", "truthParticleLink"}};

◆ m_truthSegLinkKeys

SG::ReadDecorHandleKeyArray<xAOD::MuonSegmentContainer> MuonValR4::MuonHoughTransformTester::m_truthSegLinkKeys {this, "TruthSegLinkKeys", {}}
private

Declare the dependencies on the decorations.

Definition at line 101 of file MuonHoughTransformTester.h.

101{this, "TruthSegLinkKeys", {}};

◆ m_truthSegmentKey

SG::ReadHandleKey<xAOD::MuonSegmentContainer> MuonValR4::MuonHoughTransformTester::m_truthSegmentKey {this, "TruthSegmentKey","MuonTruthSegments", "truth segment container"}
private

Key to the truth segment.

Definition at line 99 of file MuonHoughTransformTester.h.

99{this, "TruthSegmentKey","MuonTruthSegments", "truth segment container"};

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

◆ m_visionTool

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

Pattern visualization tool.

Definition at line 287 of file MuonHoughTransformTester.h.

287{this, "VisualizationTool", ""};

◆ m_writeSpacePoints

Gaudi::Property<bool> MuonValR4::MuonHoughTransformTester::m_writeSpacePoints
private
Initial value:
{this, "writeSpacePoints", false,
"Toggle whether the particular space poitns shall be written"}

Definition at line 120 of file MuonHoughTransformTester.h.

120 {this, "writeSpacePoints", false,
121 "Toggle whether the particular space poitns shall be written"};

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