ATLAS Offline Software
Classes | Public Member Functions | Protected Member Functions | Private Types | Private Member Functions | Private Attributes | List of all members
MuonR4::TruthSegmentMaker Class Reference

#include <TruthSegmentMaker.h>

Inheritance diagram for MuonR4::TruthSegmentMaker:
Collaboration diagram for MuonR4::TruthSegmentMaker:

Classes

struct  WriteDecorHolder
 Helper struct to ship the write DecorHandles and the reference to the output segment container through the class methods. More...
 

Public Member Functions

 ~TruthSegmentMaker ()=default
 
StatusCode initialize () override final
 
StatusCode execute (const EventContext &ctx) const override
 
virtual StatusCode sysInitialize () override
 Override sysInitialize. More...
 
virtual bool isClonable () const override
 Specify if the algorithm is clonable. More...
 
virtual unsigned int cardinality () const override
 Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant. More...
 
virtual StatusCode sysExecute (const EventContext &ctx) override
 Execute an algorithm. More...
 
virtual const DataObjIDColl & extraOutputDeps () const override
 Return the list of extra output dependencies. More...
 
virtual bool filterPassed (const EventContext &ctx) const
 
virtual void setFilterPassed (bool state, const EventContext &ctx) const
 
ServiceHandle< StoreGateSvc > & evtStore ()
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
 
const ServiceHandle< StoreGateSvc > & evtStore () const
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
 
const ServiceHandle< StoreGateSvc > & detStore () const
 The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc. More...
 
virtual StatusCode sysStart () override
 Handle START transition. More...
 
virtual std::vector< Gaudi::DataHandle * > inputHandles () const override
 Return this algorithm's input handles. More...
 
virtual std::vector< Gaudi::DataHandle * > outputHandles () const override
 Return this algorithm's output handles. More...
 
Gaudi::Details::PropertyBase & declareProperty (Gaudi::Property< T, V, H > &t)
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &)
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, T &property, const std::string &doc="none")
 Declare a new Gaudi property. More...
 
void updateVHKA (Gaudi::Details::PropertyBase &)
 
MsgStream & msg () const
 
MsgStream & msg (const MSG::Level lvl) const
 
bool msgLvl (const MSG::Level lvl) const
 

Protected Member Functions

void renounceArray (SG::VarHandleKeyArray &handlesArray)
 remove all handles from I/O resolution More...
 
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. More...
 

Private Types

using HitPosTuple_t = std::tuple< const xAOD::MuonSimHit *, Amg::Vector3D, Amg::Vector3D >
 Tuple consisting out of pointer to the sim hit and the position & direction expressed in the chamber's frame. More...
 
using SimHitVec_t = std::vector< HitPosTuple_t >
 
using EleLink_t = ElementLink< xAOD::MuonSimHitContainer >
 
using HitLinkVec_t = std::vector< EleLink_t >
 
using LinkDecor_t = SG::WriteDecorHandle< xAOD::MuonSegmentContainer, HitLinkVec_t >
 
using FloatDecor_t = SG::WriteDecorHandle< xAOD::MuonSegmentContainer, float >
 
using SegPars_t = xAOD::MeasVector< SegmentFit::toInt(SegmentFit::ParamDefs::nPars)>
 
using SegParDecor_t = SG::WriteDecorHandle< xAOD::MuonSegmentContainer, SegPars_t >
 
typedef ServiceHandle< StoreGateSvcStoreGateSvc_t
 

Private Member Functions

Amg::Transform3D toChamber (const ActsGeometryContext &gctx, const Identifier &chanId) const
 Returns the transform from the local simHit frame -> chamber frame. More...
 
xAOD::MuonSegmentconstructSegmentFromHits (const EventContext &ctx, const Amg::Transform3D &locToGlob, const SimHitVec_t &hits, WriteDecorHolder &writeShip) const
 Takes a list of related sim hits and transforms them into a truth segment. More...
 
void buildSegmentsFromBkg (const EventContext &ctx, const Amg::Transform3D &locToGlob, const SimHitVec_t &simHits, WriteDecorHolder &writeShip) const
 Attempts to assemble truth segments from a list of loose sim hits, i.e., the hits are stemming from a muon but due to the pile-up truth record, no HepMCParticles are recorded. More...
 
float muonPt (const xAOD::MuonSimHit &hit, const Amg::Vector3D &globDir) const
 Returns the muon pt from the sim hit. More...
 
float hitUncertainty (const EventContext &ctx, const xAOD::MuonSimHit &hit) const
 Evaluates the hit uncertainty. More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyArrayType &)
 specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleType &)
 specialization for handling Gaudi::Property<SG::VarHandleBase> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &t, const SG::NotHandleType &)
 specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...
 

Private Attributes

ServiceHandle< Muon::IMuonIdHelperSvcm_idHelperSvc {this, "IdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}
 IdHelperSvc to decode the Identifiers. More...
 
SG::ReadHandleKeyArray< xAOD::MuonSimHitContainerm_readKeys {this, "SimHitKeys", {}}
 List of sim hit containers from which the truth segments shall be retrieved. More...
 
SG::ReadHandleKey< ActsGeometryContextm_geoCtxKey {this, "AlignmentKey", "ActsAlignment", "cond handle key"}
 Key to the geometry context. More...
 
SG::WriteHandleKey< xAOD::MuonSegmentContainerm_segmentKey {this, "WriteKey", "TruthSegmentsR4"}
 Key under which the segment Container will be recorded in StoreGate. More...
 
SG::WriteDecorHandleKey< xAOD::MuonSegmentContainerm_eleLinkKey {this, "SimHitLink", m_segmentKey, "simHitLinks"}
 Decoration key of the associated sim hit links. More...
 
SG::WriteDecorHandleKey< xAOD::MuonSegmentContainerm_ptKey {this, "PtKey", m_segmentKey, "pt"}
 Decoration key of the associated particle pt. More...
 
SG::WriteDecorHandleKey< xAOD::MuonSegmentContainerm_locParKey {this, "LocParKey", m_segmentKey,"localSegPars"}
 Decoration key of the local parameters. More...
 
SG::WriteDecorHandleKey< xAOD::MuonSegmentContainerm_qKey {this, "qKey", m_segmentKey, "charge"}
 Decoration key of the muon charge
More...
 
Gaudi::Property< bool > m_useOnlyMuonHits {this, "useOnlyMuonHits", true}
 Build segments from muon hits only. More...
 
Gaudi::Property< bool > m_includePileUpHits {this, "includePileUpHits", true}
 Construct segments from pile-up hits without GenParticleLink. More...
 
Gaudi::Property< float > m_pileUpHitMinE {this, "energyThresholdPileUp", 1.*Gaudi::Units::GeV}
 Minimum energy threshold for pile up hits to be converted
More...
 
Gaudi::Property< float > m_pileUpHitELoss {this, "pileUpHitELoss", 5.*Gaudi::Units::MeV}
 Maximum energy loss between two pile-up hits. More...
 
Gaudi::Property< float > m_pileUpHitAngleCone {this, "pileUpHitAngleCone", 1.*Gaudi::Units::deg}
 Maximum scattering angle between two pile-up hits. More...
 
Gaudi::Property< float > m_pileUpHitDistance {this, "pileUpHitDistance", 2.*Gaudi::Units::m}
 Maximum separation between two pile-up hits. More...
 
Gaudi::Property< float > m_idCylinderR {this, "IdCylinderR", 1.1*Gaudi::Units::m}
 ID / ITk cylinder radius. More...
 
Gaudi::Property< float > m_idCylinderHalfZ {this, "IdCylinderHalfZ", 3.2*Gaudi::Units::m}
 ID / Itk cylinder half length. More...
 
const MuonGMR4::MuonDetectorManagerm_detMgr {nullptr}
 Pointer to the muon readout geometry. More...
 
SG::ReadCondHandleKey< MuonCalib::MdtCalibDataContainerm_mdtCalibKey
 Data dependency on the Mdt calibration container to calculate the uncertainty. More...
 
SG::ReadCondHandleKey< NswErrorCalibDatam_nswUncertKey
 Data dependency on the Nsw calibration container to estimate the uncertaintys. More...
 
DataObjIDColl m_extendedExtraObjects
 Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks. More...
 
StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default) More...
 
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default) More...
 
std::vector< SG::VarHandleKeyArray * > m_vhka
 
bool m_varHandleArraysDeclared
 

Detailed Description

Definition at line 32 of file TruthSegmentMaker.h.

Member Typedef Documentation

◆ EleLink_t

Definition at line 55 of file TruthSegmentMaker.h.

◆ FloatDecor_t

Definition at line 58 of file TruthSegmentMaker.h.

◆ HitLinkVec_t

using MuonR4::TruthSegmentMaker::HitLinkVec_t = std::vector<EleLink_t>
private

Definition at line 56 of file TruthSegmentMaker.h.

◆ HitPosTuple_t

Tuple consisting out of pointer to the sim hit and the position & direction expressed in the chamber's frame.

Definition at line 50 of file TruthSegmentMaker.h.

◆ LinkDecor_t

Definition at line 57 of file TruthSegmentMaker.h.

◆ SegParDecor_t

Definition at line 61 of file TruthSegmentMaker.h.

◆ SegPars_t

Definition at line 60 of file TruthSegmentMaker.h.

◆ SimHitVec_t

Definition at line 51 of file TruthSegmentMaker.h.

◆ StoreGateSvc_t

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

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ ~TruthSegmentMaker()

MuonR4::TruthSegmentMaker::~TruthSegmentMaker ( )
default

Member Function Documentation

◆ buildSegmentsFromBkg()

void MuonR4::TruthSegmentMaker::buildSegmentsFromBkg ( const EventContext &  ctx,
const Amg::Transform3D locToGlob,
const SimHitVec_t simHits,
WriteDecorHolder writeShip 
) const
private

Attempts to assemble truth segments from a list of loose sim hits, i.e., the hits are stemming from a muon but due to the pile-up truth record, no HepMCParticles are recorded.

Hits are paired based on kinetic energy, and line compatibility

Parameters
ctxEventContext to access the uncertainty conditions data during segment building
locToGlobTransform from the sector -> global frame
simHitsCollection of all sim hits in the sector without truth link
writeShipHelper struct carrying the output segment container & the associated decorators

Closest beam spot approach

First compare the energies

Calculate the approach to the cylinder

Definition at line 111 of file TruthSegmentMaker.cxx.

114  {
115 
116  ATH_MSG_VERBOSE("Assemble segments from "<<simHits.size()<<" background hits.");
117  std::vector<char> alreadyUsed(simHits.size(), 0);
118  for (std::size_t h = 0 ;h < simHits.size(); ++h) {
119  if (alreadyUsed[h]) {
120  continue;
121  }
122  const auto& [refHit, refPos, refDir] = simHits[h];
124  ATH_MSG_VERBOSE("Try to find other hits on trajectory given by "<<m_idHelperSvc->toString(refHit->identify())
125  <<", "<<Amg::toString(refPos)<<", "<<Amg::toString(refDir)<<", E: "<<refHit->kineticEnergy() / Gaudi::Units::GeV<<" [GeV] "
126  <<", pt: "<<muonPt(*refHit, locToGlob.linear()* refDir)/ Gaudi::Units::GeV <<" [GeV].");
127  std::vector<std::size_t> indicesOnSeg{h};
128  for (std::size_t h1 = h +1; h1 < simHits.size(); ++h1) {
129  if (alreadyUsed[h1]) {
130  continue;
131  }
132  const auto&[testHit, testPos, testDir] = simHits[h1];
133 
135  if (std::abs(refHit->kineticEnergy() - testHit->kineticEnergy()) > m_pileUpHitELoss) {
136  continue;
137  }
138  if (MC::charge(refHit) != MC::charge(testHit)){
139  continue;
140  }
141  const double angleDir = Amg::angle(testDir, refDir);
142  const double hitSep = std::abs(Amg::signedDistance(refPos, refDir, testPos,testDir));
143  ATH_MSG_VERBOSE("Test hit "<<m_idHelperSvc->toString(testHit->identify())<<", "
144  <<Amg::toString(testPos)<<", testDir: "<<Amg::toString(testDir)
145  <<", E: "<<testHit->kineticEnergy() / (Gaudi::Units::GeV)<<", angle: "
146  <<(angleDir / Gaudi::Units::deg) <<", distance: "<<hitSep<<".");
147 
148  if (std::abs(angleDir) > m_pileUpHitAngleCone || hitSep > m_pileUpHitDistance) {
149  continue;
150  }
151  indicesOnSeg.push_back(h1);
152  }
154  const Amg::Vector3D globPos = locToGlob * refPos;
155  const Amg::Vector3D globDir = locToGlob.linear() * refDir;
156  const Amg::Vector3D perigee = globPos + Amg::intersect<3>(Amg::Vector3D::Zero(), Amg::Vector3D::UnitZ(),
157  globPos, globDir).value_or(0.) * globDir;
158  ATH_MSG_VERBOSE("Found "<<indicesOnSeg.size()<<" matching hits. Closest perigee "
159  <<Amg::toString(perigee)<<", "<<m_idCylinderR<<", "<<m_idCylinderHalfZ);
160  SimHitVec_t hitsOnSeg{};
161  std::ranges::transform(indicesOnSeg,
162  std::back_inserter(hitsOnSeg),
163  [&simHits](const auto idx) {
164  return simHits[idx];
165  });
166  if (perigee.perp() > m_idCylinderR || std::abs(perigee.z()) > m_idCylinderHalfZ ||
167  constructSegmentFromHits(ctx, locToGlob, hitsOnSeg, out)){
168  std::ranges::for_each(indicesOnSeg,
169  [&alreadyUsed](const auto idx) {
170  alreadyUsed[idx] = 1;
171  });
172  }
173  }
174  }

◆ cardinality()

unsigned int AthCommonReentrantAlgorithm< Gaudi::Algorithm >::cardinality
overridevirtualinherited

Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.

Override this to return 0 for reentrant algorithms.

Override this to return 0 for reentrant algorithms.

Definition at line 75 of file AthCommonReentrantAlgorithm.cxx.

64 {
65  return 0;
66 }

◆ constructSegmentFromHits()

xAOD::MuonSegment * MuonR4::TruthSegmentMaker::constructSegmentFromHits ( const EventContext &  ctx,
const Amg::Transform3D locToGlob,
const SimHitVec_t hits,
WriteDecorHolder writeShip 
) const
private

Takes a list of related sim hits and transforms them into a truth segment.

If there're less than 3 precision hits, then no segment is created.

Parameters
ctxEventContext to fetch the hit uncertainties from store gate
locToGlobCoordinate transform from sector frame -> global frame
hitsList of hits to be combined onto a single segment
writeShipHelper struct carrying the output segment container & the associated decorators

Find the hit which is closest to the refernce plane to express the segment at the centre plane

Express the simulated hit in the center of the chamber

Associate the hits to the segment and calculate the expected hits

Definition at line 176 of file TruthSegmentMaker.cxx.

179  {
182  auto refHit_itr = std::ranges::min_element(simHits,
183  [](const HitPosTuple_t& a, const HitPosTuple_t& b){
184  return std::abs(std::get<1>(a).z()) <std::abs(std::get<1>(b).z());
185  });
186 
187  const auto [simHit, localPos, chamberDir] = (*refHit_itr);
188  ATH_MSG_VERBOSE("Create segement from hit: "<<m_idHelperSvc->toString(simHit->identify())<<
189  " pdgId: "<<simHit->pdgId()<<", energy: "<<simHit->kineticEnergy()
190  <<", genParticle: "<<simHit->genParticleLink().cptr());
191  const Identifier segId{simHit->identify()};
192 
194  const double distance = Amg::intersect<3>(localPos, chamberDir, Amg::Vector3D::UnitZ(), 0.).value_or(0.);
195  const Amg::Vector3D chamberPos = localPos + distance*chamberDir;
196 
197  const Amg::Vector3D globPos = locToGlob * chamberPos;
198  const Amg::Vector3D globDir = locToGlob.linear() * chamberDir;
199  HitLinkVec_t associatedHits{};
200  unsigned nMdt{0}, nRpcEta{0}, nRpcPhi{0}, nTgcEta{0}, nTgcPhi{0};
201  unsigned nMm{0}, nStgcEta{0}, nStgcPhi{0};
203  float chi2{0.f};
204  for (const auto& [assocMe, pos, dir] : simHits) {
205  chi2 += std::pow(Amg::signedDistance(pos,dir, chamberPos, chamberDir) / hitUncertainty(ctx,*assocMe), 2);
206  const MuonGMR4::MuonReadoutElement* assocRE = m_detMgr->getReadoutElement(assocMe->identify());
207  switch (assocRE->detectorType()) {
209  ++nMdt;
210  break;
212  auto castRE{static_cast<const MuonGMR4::RpcReadoutElement*>(assocRE)};
213  if (castRE->nEtaStrips()) ++nRpcEta;
214  if (castRE->nPhiStrips()) ++nRpcPhi;
215  break;
217  auto castRE{static_cast<const MuonGMR4::TgcReadoutElement*>(assocRE)};
218  const IdentifierHash gapHash = assocRE->measurementHash(assocMe->identify());
219  if (castRE->numStrips(gapHash)){
220  ++nTgcPhi;
221  }
222  if (castRE->numWireGangs(gapHash)) {
223  ++nTgcEta;
224  }
225  break;
227  ++nStgcEta;
228  ++nStgcPhi;
229  break;
231  ++nMm;
232  break;
233  } default:
234  ATH_MSG_WARNING("Csc are not defined "<<m_idHelperSvc->toString(simHit->identify()));
235  }
236  ATH_MSG_VERBOSE("Associate hit "<<m_idHelperSvc->toString(assocMe->identify())
237  <<" pdgId: "<<assocMe->pdgId()<<", energy: "<<assocMe->kineticEnergy()
238  <<", genParticle: "<<assocMe->genParticleLink().cptr()<<", beta: "<<simHit->beta()
239  <<" global time: "<<simHit->globalTime()<<", pos: "<<Amg::toString(pos)
240  <<", dir: "<<Amg::toString(dir));
241  EleLink_t link{*static_cast<const xAOD::MuonSimHitContainer*>(assocMe->container()), assocMe->index()};
242  associatedHits.push_back(std::move(link));
243  }
244  int nPrecisionHits = nMdt + nMm + nStgcEta;
245  int nPhiLayers = nTgcPhi + nRpcPhi + nStgcPhi;
246  // if nMdt + nMm + nStgcEta < 3, do not create a segment
247  if (nPrecisionHits < 3) {
248  return nullptr;
249  }
250 
251  xAOD::MuonSegment* truthSegment = out.segments.push_back(std::make_unique<xAOD::MuonSegment>());
252  out.ptDecor(*truthSegment) = muonPt(*simHit, globDir);
253  out.chargeDecor(*truthSegment) = MC::charge(simHit);
254  SegPars_t& locPars{out.paramDecor(*truthSegment)};
255  locPars[toInt(ParamDefs::x0)] = chamberPos.x();
256  locPars[toInt(ParamDefs::y0)] = chamberPos.y();
257  constexpr float betaLowLimit = 1.e-6;
258  locPars[toInt(ParamDefs::time)] = simHit->globalTime() + distance *c_inv / std::max(simHit->beta(), betaLowLimit);
259  locPars[toInt(ParamDefs::theta)] = chamberDir.theta();
260  locPars[toInt(ParamDefs::phi)] = chamberDir.phi();
261 
262  truthSegment->setPosition(globPos.x(), globPos.y(), globPos.z());
263  truthSegment->setDirection(globDir.x(), globDir.y(), globDir.z());
264  truthSegment->setT0Error(locPars[toInt(ParamDefs::time)], 0.);
265 
266  truthSegment->setNHits(nPrecisionHits, nPhiLayers, nTgcEta + nRpcEta);
267  truthSegment->setIdentifier(m_idHelperSvc->sector(segId),
268  m_idHelperSvc->chamberIndex(segId),
269  m_idHelperSvc->stationEta(segId),
270  m_idHelperSvc->technologyIndex(segId));
271  // adding chi2 and ndof (nHits - 5 for 2 position, 2 direction and 1 time)
272  if (nPhiLayers == 0){
273  truthSegment->setFitQuality(chi2, (nPrecisionHits + nTgcEta + nRpcEta - 3));
274  } else {
275  truthSegment->setFitQuality(chi2, (nPrecisionHits + nPhiLayers + nTgcEta + nRpcEta - 5));
276  }
277  out.hitLinkDecor(*truthSegment) = std::move(associatedHits);
278  return truthSegment;
279  }

◆ declareGaudiProperty() [1/4]

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

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

Definition at line 170 of file AthCommonDataStore.h.

172  {
173  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
174  hndl.value(),
175  hndl.documentation());
176 
177  }

◆ declareGaudiProperty() [2/4]

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

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

Definition at line 156 of file AthCommonDataStore.h.

158  {
159  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
160  hndl.value(),
161  hndl.documentation());
162 
163  }

◆ declareGaudiProperty() [3/4]

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

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

Definition at line 184 of file AthCommonDataStore.h.

186  {
187  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
188  hndl.value(),
189  hndl.documentation());
190  }

◆ declareGaudiProperty() [4/4]

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

specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray>

Definition at line 199 of file AthCommonDataStore.h.

200  {
201  return PBASE::declareProperty(t);
202  }

◆ declareProperty() [1/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( const std::string &  name,
SG::VarHandleBase hndl,
const std::string &  doc,
const SG::VarHandleType  
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
hndlObject holding the property value.
docDocumentation string for the property.

This is the version for types that derive from SG::VarHandleBase. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 245 of file AthCommonDataStore.h.

249  {
250  this->declare(hndl.vhKey());
251  hndl.vhKey().setOwner(this);
252 
253  return PBASE::declareProperty(name,hndl,doc);
254  }

◆ declareProperty() [2/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( const std::string &  name,
SG::VarHandleKey hndl,
const std::string &  doc,
const SG::VarHandleKeyType  
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
hndlObject holding the property value.
docDocumentation string for the property.

This is the version for types that derive from SG::VarHandleKey. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 221 of file AthCommonDataStore.h.

225  {
226  this->declare(hndl);
227  hndl.setOwner(this);
228 
229  return PBASE::declareProperty(name,hndl,doc);
230  }

◆ declareProperty() [3/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( const std::string &  name,
SG::VarHandleKeyArray hndArr,
const std::string &  doc,
const SG::VarHandleKeyArrayType  
)
inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

263  {
264 
265  // std::ostringstream ost;
266  // ost << Algorithm::name() << " VHKA declareProp: " << name
267  // << " size: " << hndArr.keys().size()
268  // << " mode: " << hndArr.mode()
269  // << " vhka size: " << m_vhka.size()
270  // << "\n";
271  // debug() << ost.str() << endmsg;
272 
273  hndArr.setOwner(this);
274  m_vhka.push_back(&hndArr);
275 
276  Gaudi::Details::PropertyBase* p = PBASE::declareProperty(name, hndArr, doc);
277  if (p != 0) {
278  p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
279  } else {
280  ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray "
281  << name);
282  }
283 
284  return p;
285 
286  }

◆ declareProperty() [4/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( const std::string &  name,
T &  property,
const std::string &  doc,
const SG::NotHandleType  
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
propertyObject holding the property value.
docDocumentation string for the property.

This is the generic version, for types that do not derive from SG::VarHandleKey. It just forwards to the base class version of declareProperty.

Definition at line 333 of file AthCommonDataStore.h.

337  {
338  return PBASE::declareProperty(name, property, doc);
339  }

◆ declareProperty() [5/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( const std::string &  name,
T &  property,
const std::string &  doc = "none" 
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
propertyObject holding the property value.
docDocumentation string for the property.

This dispatches to either the generic declareProperty or the one for VarHandle/Key/KeyArray.

Definition at line 352 of file AthCommonDataStore.h.

355  {
356  typedef typename SG::HandleClassifier<T>::type htype;
357  return declareProperty (name, property, doc, htype());
358  }

◆ declareProperty() [6/6]

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

Definition at line 145 of file AthCommonDataStore.h.

145  {
146  typedef typename SG::HandleClassifier<T>::type htype;
148  }

◆ detStore()

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

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

Definition at line 95 of file AthCommonDataStore.h.

95 { return m_detStore; }

◆ evtStore() [1/2]

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

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

Definition at line 85 of file AthCommonDataStore.h.

85 { return m_evtStore; }

◆ evtStore() [2/2]

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

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

Definition at line 90 of file AthCommonDataStore.h.

90 { return m_evtStore; }

◆ execute()

StatusCode MuonR4::TruthSegmentMaker::execute ( const EventContext &  ctx) const
override

Definition at line 281 of file TruthSegmentMaker.cxx.

281  {
282  const ActsGeometryContext* gctx{nullptr};
283  ATH_CHECK(SG::get(gctx, m_geoCtxKey, ctx));
284 
285  using HitsPerParticle = std::unordered_map<HepMC::ConstGenParticlePtr, SimHitVec_t>;
286  using HitCollector = std::unordered_map<const MuonGMR4::SpectrometerSector*, HitsPerParticle>;
287  HitCollector hitCollector{};
288 
290  const xAOD::MuonSimHitContainer* simHits{nullptr};
291  ATH_CHECK(SG::get(simHits, key, ctx));
292  for (const xAOD::MuonSimHit* simHit : *simHits) {
293  const MuonGMR4::MuonReadoutElement* reElement = m_detMgr->getReadoutElement(simHit->identify());
294  const MuonGMR4::SpectrometerSector* id{reElement->msSector()};
295  auto genLink = simHit->genParticleLink();
296  HepMC::ConstGenParticlePtr genParticle = nullptr;
297  if (genLink.isValid()){
298  genParticle = genLink.cptr();
299  }
300  if ( (!genParticle && (!m_includePileUpHits || simHit->kineticEnergy() < m_pileUpHitMinE
301  || !MC::isMuon(simHit))) || (m_useOnlyMuonHits && !MC::isMuon(simHit))) {
302  ATH_MSG_VERBOSE("Skip hit "<<m_idHelperSvc->toString(simHit->identify())<<
303  " pdgId: "<<simHit->pdgId()<<", energy: "<<simHit->kineticEnergy()
304  <<", genParticle: "<<genParticle);
305  continue;
306  }
307  const Amg::Transform3D toChTrf{toChamber(*gctx, simHit->identify())};
308  hitCollector[id][genParticle].emplace_back(simHit,
309  toChTrf *xAOD::toEigen(simHit->localPosition()),
310  toChTrf.linear()* xAOD::toEigen(simHit->localDirection()));
311  }
312  }
313 
314  SG::WriteHandle writeHandle{m_segmentKey, ctx};
315  ATH_CHECK(writeHandle.record(std::make_unique<xAOD::MuonSegmentContainer>(),
316  std::make_unique<xAOD::MuonSegmentAuxContainer>()));
317 
318  WriteDecorHolder writerHolder{*writeHandle,*this, ctx};
319 
320  for (auto& [chamber, collectedParts] : hitCollector) {
321  const Amg::Transform3D& locToGlob{chamber->localToGlobalTrans(*gctx)};
322  for (auto& [particle, simHits]: collectedParts) {
323  /* Sort hits by local z */
324  std::ranges::stable_sort(simHits,[](const HitPosTuple_t& a, const HitPosTuple_t& b){
325  return std::get<1>(a).z() < std::get<1>(b).z();
326  });
327  if (!particle) {
328  buildSegmentsFromBkg(ctx, locToGlob, simHits, writerHolder);
329  continue;
330  }
331  constructSegmentFromHits(ctx, locToGlob, simHits, writerHolder);
332  }
333  }
334  ATH_MSG_DEBUG("Constructed "<<writeHandle->size()<<" truth segments in total ");
335  return StatusCode::SUCCESS;
336  }

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase &  ExtraDeps)
protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

◆ extraOutputDeps()

const DataObjIDColl & AthCommonReentrantAlgorithm< Gaudi::Algorithm >::extraOutputDeps
overridevirtualinherited

Return the list of extra output dependencies.

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

Definition at line 94 of file AthCommonReentrantAlgorithm.cxx.

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

◆ filterPassed()

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

Definition at line 96 of file AthCommonReentrantAlgorithm.h.

96  {
97  return execState( ctx ).filterPassed();
98  }

◆ hitUncertainty()

float MuonR4::TruthSegmentMaker::hitUncertainty ( const EventContext &  ctx,
const xAOD::MuonSimHit hit 
) const
private

Evaluates the hit uncertainty.

For the precision detectors, the uncertainties are fetched from the COOL DB. For the remaining detectors, the uncertainties are based on the strip pitch

Parameters
ctxEventContext to fetch the constants from store gate
hitReference to the sim hit

Definition at line 52 of file TruthSegmentMaker.cxx.

52  {
53  const Identifier hitId{hit.identify()};
54  switch (const auto techIdx = m_idHelperSvc->technologyIndex(hitId)) {
56  case MDT: {
57  const MuonCalib::MdtCalibDataContainer* calibCont{nullptr};
58  if (!SG::get(calibCont, m_mdtCalibKey, ctx).isSuccess()) {
59  THROW_EXCEPTION("Failed to retrieve Mdt calib constants");
60  }
61  const auto& rtCalib{calibCont->getCalibData(hitId, msgStream())->rtRelation};
62  const double driftTime = rtCalib->tr()->driftTime(hit.localPosition().perp()).value_or(rtCalib->tr()->maxRadius());
63  return rtCalib->rtRes()->resolution(driftTime);
64  } case RPC: {
65  const auto* re = m_detMgr->getRpcReadoutElement(hitId);
66  return re->stripEtaPitch() / std::sqrt(12.);
67  } case TGC: {
68  const auto* re = m_detMgr->getTgcReadoutElement(hitId);
69  const IdentifierHash measHash = re->measurementHash(hitId);
70  const auto& design = re->wireGangLayout(measHash);
71  return design.stripPitch() / std::sqrt(12.);
72  } case STGC:
73  case MM: {
74  const NswErrorCalibData* errorCalibDB{nullptr};
75  if (!SG::get(errorCalibDB, m_nswUncertKey, ctx).isSuccess()) {
76  THROW_EXCEPTION("Failed to retrieve the STGC calibration constants");
77  }
78  NswErrorCalibData::Input errorCalibInput{};
79  errorCalibInput.stripId= hitId;
80  errorCalibInput.locTheta = M_PI - hit.localDirection().theta();
81  if (techIdx == STGC) {
82  errorCalibInput.clusterAuthor = 3; // centroid
83  } else {
84  errorCalibInput.clusterAuthor=66; // cluster time projection method
85  }
86  return errorCalibDB->clusterUncertainty(errorCalibInput);
87  }
88  default:
89  break;
90  }
91  return 0.;
92  }

◆ initialize()

StatusCode MuonR4::TruthSegmentMaker::initialize ( )
finaloverride

Definition at line 31 of file TruthSegmentMaker.cxx.

31  {
32  ATH_CHECK(m_idHelperSvc.retrieve());
33  ATH_CHECK(m_readKeys.initialize());
34  if (m_readKeys.empty()){
35  ATH_MSG_ERROR("No simulated hit containers have been parsed to build the segments from ");
36  return StatusCode::FAILURE;
37  }
40 
44 
47 
50  return StatusCode::SUCCESS;
51  }

◆ inputHandles()

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

Return this algorithm's input handles.

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

◆ isClonable()

bool AthCommonReentrantAlgorithm< Gaudi::Algorithm >::isClonable
overridevirtualinherited

Specify if the algorithm is clonable.

Reentrant algorithms are clonable.

Definition at line 68 of file AthCommonReentrantAlgorithm.cxx.

52 {
53  // Reentrant algorithms are clonable.
54  return true;
55 }

◆ msg() [1/2]

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

Definition at line 24 of file AthCommonMsg.h.

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

◆ msg() [2/2]

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

Definition at line 27 of file AthCommonMsg.h.

27  {
28  return this->msgStream(lvl);
29  }

◆ msgLvl()

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

Definition at line 30 of file AthCommonMsg.h.

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

◆ muonPt()

float MuonR4::TruthSegmentMaker::muonPt ( const xAOD::MuonSimHit hit,
const Amg::Vector3D globDir 
) const
private

Returns the muon pt from the sim hit.

Definition at line 93 of file TruthSegmentMaker.cxx.

93  {
94  const auto& link = hit.genParticleLink();
95  if (link.isValid()) {
96  return link->momentum().perp();
97  }
98  const float e{hit.kineticEnergy()},m{hit.mass()};
99  return std::sqrt(std::max(e*e - m*m, 0.f)) * std::sin(globDir.theta());
100  }

◆ outputHandles()

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

Return this algorithm's output handles.

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

◆ renounce()

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

Definition at line 380 of file AthCommonDataStore.h.

381  {
382  h.renounce();
383  PBASE::renounce (h);
384  }

◆ renounceArray()

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

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364  {
365  handlesArray.renounce();
366  }

◆ setFilterPassed()

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

Definition at line 100 of file AthCommonReentrantAlgorithm.h.

100  {
101  execState( ctx ).setFilterPassed( state );
102  }

◆ sysExecute()

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

Execute an algorithm.

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

Definition at line 85 of file AthCommonReentrantAlgorithm.cxx.

77 {
78  return BaseAlg::sysExecute (ctx);
79 }

◆ sysInitialize()

StatusCode AthCommonReentrantAlgorithm< Gaudi::Algorithm >::sysInitialize
overridevirtualinherited

Override sysInitialize.

Override sysInitialize from the base class.

Loop through all output handles, and if they're WriteCondHandles, automatically register them and this Algorithm with the CondSvc

Scan through all outputHandles, and if they're WriteCondHandles, register them with the CondSvc

Reimplemented from AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >.

Reimplemented in InputMakerBase, and HypoBase.

Definition at line 61 of file AthCommonReentrantAlgorithm.cxx.

107  {
109 
110  if (sc.isFailure()) {
111  return sc;
112  }
113 
114  ServiceHandle<ICondSvc> cs("CondSvc",name());
115  for (auto h : outputHandles()) {
116  if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
117  // do this inside the loop so we don't create the CondSvc until needed
118  if ( cs.retrieve().isFailure() ) {
119  ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
120  return StatusCode::SUCCESS;
121  }
122  if (cs->regHandle(this,*h).isFailure()) {
123  sc = StatusCode::FAILURE;
124  ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
125  << " with CondSvc");
126  }
127  }
128  }
129  return sc;
130 }

◆ sysStart()

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

Handle START transition.

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

◆ toChamber()

Amg::Transform3D MuonR4::TruthSegmentMaker::toChamber ( const ActsGeometryContext gctx,
const Identifier chanId 
) const
private

Returns the transform from the local simHit frame -> chamber frame.

Parameters
gctxGeometry context to align the chambers within ATLAS
chanIdIdentifier of the channel for which the transform shall be fetched

Definition at line 101 of file TruthSegmentMaker.cxx.

102  {
104  const IdentifierHash trfHash{reEle->detectorType() == ActsTrk::DetectorType::Mdt ?
105  reEle->measurementHash(chanId) :
106  reEle->layerHash(chanId)};
107  return reEle->msSector()->globalToLocalTrans(gctx) * reEle->localToGlobalTrans(gctx, trfHash);
108 
109  }

◆ updateVHKA()

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

Definition at line 308 of file AthCommonDataStore.h.

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

Member Data Documentation

◆ m_detMgr

const MuonGMR4::MuonDetectorManager* MuonR4::TruthSegmentMaker::m_detMgr {nullptr}
private

Pointer to the muon readout geometry.

Definition at line 155 of file TruthSegmentMaker.h.

◆ m_detStore

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

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_eleLinkKey

SG::WriteDecorHandleKey<xAOD::MuonSegmentContainer> MuonR4::TruthSegmentMaker::m_eleLinkKey {this, "SimHitLink", m_segmentKey, "simHitLinks"}
private

Decoration key of the associated sim hit links.

Definition at line 131 of file TruthSegmentMaker.h.

◆ m_evtStore

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

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_extendedExtraObjects

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

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

Empty if no symlinks were found.

Definition at line 114 of file AthCommonReentrantAlgorithm.h.

◆ m_geoCtxKey

SG::ReadHandleKey<ActsGeometryContext> MuonR4::TruthSegmentMaker::m_geoCtxKey {this, "AlignmentKey", "ActsAlignment", "cond handle key"}
private

Key to the geometry context.

Needed to align the hits inside ATLAS

Definition at line 127 of file TruthSegmentMaker.h.

◆ m_idCylinderHalfZ

Gaudi::Property<float> MuonR4::TruthSegmentMaker::m_idCylinderHalfZ {this, "IdCylinderHalfZ", 3.2*Gaudi::Units::m}
private

ID / Itk cylinder half length.

Definition at line 153 of file TruthSegmentMaker.h.

◆ m_idCylinderR

Gaudi::Property<float> MuonR4::TruthSegmentMaker::m_idCylinderR {this, "IdCylinderR", 1.1*Gaudi::Units::m}
private

ID / ITk cylinder radius.

Definition at line 151 of file TruthSegmentMaker.h.

◆ m_idHelperSvc

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

IdHelperSvc to decode the Identifiers.

Definition at line 123 of file TruthSegmentMaker.h.

◆ m_includePileUpHits

Gaudi::Property<bool> MuonR4::TruthSegmentMaker::m_includePileUpHits {this, "includePileUpHits", true}
private

Construct segments from pile-up hits without GenParticleLink.

Definition at line 141 of file TruthSegmentMaker.h.

◆ m_locParKey

SG::WriteDecorHandleKey<xAOD::MuonSegmentContainer> MuonR4::TruthSegmentMaker::m_locParKey {this, "LocParKey", m_segmentKey,"localSegPars"}
private

Decoration key of the local parameters.

Definition at line 135 of file TruthSegmentMaker.h.

◆ m_mdtCalibKey

SG::ReadCondHandleKey<MuonCalib::MdtCalibDataContainer> MuonR4::TruthSegmentMaker::m_mdtCalibKey
private
Initial value:
{this, "MdtCalibKey", "MdtCalibConstants",
"Conditions object containing the calibrations"}

Data dependency on the Mdt calibration container to calculate the uncertainty.

Definition at line 157 of file TruthSegmentMaker.h.

◆ m_nswUncertKey

SG::ReadCondHandleKey<NswErrorCalibData> MuonR4::TruthSegmentMaker::m_nswUncertKey
private
Initial value:
{this, "NswErrorKey", "NswUncertData",
"Key of the parametrized NSW uncertainties"}

Data dependency on the Nsw calibration container to estimate the uncertaintys.

Definition at line 160 of file TruthSegmentMaker.h.

◆ m_pileUpHitAngleCone

Gaudi::Property<float> MuonR4::TruthSegmentMaker::m_pileUpHitAngleCone {this, "pileUpHitAngleCone", 1.*Gaudi::Units::deg}
private

Maximum scattering angle between two pile-up hits.

Definition at line 147 of file TruthSegmentMaker.h.

◆ m_pileUpHitDistance

Gaudi::Property<float> MuonR4::TruthSegmentMaker::m_pileUpHitDistance {this, "pileUpHitDistance", 2.*Gaudi::Units::m}
private

Maximum separation between two pile-up hits.

Definition at line 149 of file TruthSegmentMaker.h.

◆ m_pileUpHitELoss

Gaudi::Property<float> MuonR4::TruthSegmentMaker::m_pileUpHitELoss {this, "pileUpHitELoss", 5.*Gaudi::Units::MeV}
private

Maximum energy loss between two pile-up hits.

Definition at line 145 of file TruthSegmentMaker.h.

◆ m_pileUpHitMinE

Gaudi::Property<float> MuonR4::TruthSegmentMaker::m_pileUpHitMinE {this, "energyThresholdPileUp", 1.*Gaudi::Units::GeV}
private

Minimum energy threshold for pile up hits to be converted

Definition at line 143 of file TruthSegmentMaker.h.

◆ m_ptKey

SG::WriteDecorHandleKey<xAOD::MuonSegmentContainer> MuonR4::TruthSegmentMaker::m_ptKey {this, "PtKey", m_segmentKey, "pt"}
private

Decoration key of the associated particle pt.

Definition at line 133 of file TruthSegmentMaker.h.

◆ m_qKey

SG::WriteDecorHandleKey<xAOD::MuonSegmentContainer> MuonR4::TruthSegmentMaker::m_qKey {this, "qKey", m_segmentKey, "charge"}
private

Decoration key of the muon charge

Definition at line 137 of file TruthSegmentMaker.h.

◆ m_readKeys

SG::ReadHandleKeyArray<xAOD::MuonSimHitContainer> MuonR4::TruthSegmentMaker::m_readKeys {this, "SimHitKeys", {}}
private

List of sim hit containers from which the truth segments shall be retrieved.

Definition at line 125 of file TruthSegmentMaker.h.

◆ m_segmentKey

SG::WriteHandleKey<xAOD::MuonSegmentContainer> MuonR4::TruthSegmentMaker::m_segmentKey {this, "WriteKey", "TruthSegmentsR4"}
private

Key under which the segment Container will be recorded in StoreGate.

Definition at line 129 of file TruthSegmentMaker.h.

◆ m_useOnlyMuonHits

Gaudi::Property<bool> MuonR4::TruthSegmentMaker::m_useOnlyMuonHits {this, "useOnlyMuonHits", true}
private

Build segments from muon hits only.

Definition at line 139 of file TruthSegmentMaker.h.

◆ m_varHandleArraysDeclared

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

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

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

Definition at line 398 of file AthCommonDataStore.h.


The documentation for this class was generated from the following files:
python.PyKernel.retrieve
def retrieve(aClass, aKey=None)
Definition: PyKernel.py:110
NswErrorCalibData
Definition: NswErrorCalibData.h:19
AllowedVariables::e
e
Definition: AsgElectronSelectorTool.cxx:37
xAOD::MuonSimHit_v1
Definition: MuonSimHit_v1.h:18
MuonR4::TruthSegmentMaker::m_eleLinkKey
SG::WriteDecorHandleKey< xAOD::MuonSegmentContainer > m_eleLinkKey
Decoration key of the associated sim hit links.
Definition: TruthSegmentMaker.h:131
STGC
@ STGC
Definition: RegSelEnums.h:39
GeV
#define GeV
Definition: PhysicsAnalysis/TauID/TauAnalysisTools/Root/HelperFunctions.cxx:18
Trk::ParticleSwitcher::particle
constexpr ParticleHypothesis particle[PARTICLEHYPOTHESES]
the array of masses
Definition: ParticleHypothesis.h:79
MuonGMR4::MuonReadoutElement::msSector
const SpectrometerSector * msSector() const
Returns the pointer to the envelope volume enclosing all chambers in the sector.
MuonGMR4::SpectrometerSector
A spectrometer sector forms the envelope of all chambers that are placed in the same MS sector & laye...
Definition: SpectrometerSector.h:40
xAOD::MuonSimHit_v1::identify
Identifier identify() const
Returns the global ATLAS identifier of the SimHit.
Definition: xAODMuonSimHit_V1.cxx:42
MuonR4::TruthSegmentMaker::m_geoCtxKey
SG::ReadHandleKey< ActsGeometryContext > m_geoCtxKey
Key to the geometry context.
Definition: TruthSegmentMaker.h:127
Amg::angle
double angle(const Amg::Vector3D &v1, const Amg::Vector3D &v2)
calculates the opening angle between two vectors
Definition: GeoPrimitivesHelpers.h:41
xAOD::MuonSegment_v1::setPosition
void setPosition(float x, float y, float z)
Sets the global position.
Definition: MuonSegment_v1.cxx:26
MuonR4::TruthSegmentMaker::buildSegmentsFromBkg
void buildSegmentsFromBkg(const EventContext &ctx, const Amg::Transform3D &locToGlob, const SimHitVec_t &simHits, WriteDecorHolder &writeShip) const
Attempts to assemble truth segments from a list of loose sim hits, i.e., the hits are stemming from a...
Definition: TruthSegmentMaker.cxx:111
calibdata.chamber
chamber
Definition: calibdata.py:31
max
constexpr double max()
Definition: ap_fixedTest.cxx:33
MuonR4::TruthSegmentMaker::SimHitVec_t
std::vector< HitPosTuple_t > SimHitVec_t
Definition: TruthSegmentMaker.h:51
AthCommonDataStore::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 > renounce(T &h)
Definition: AthCommonDataStore.h:380
Muon::MuonStationIndex::TechnologyIndex
TechnologyIndex
enum to classify the different layers in the muon spectrometer
Definition: MuonStationIndex.h:54
ActsTrk::DetectorType::Tgc
@ Tgc
Resitive Plate Chambers.
MuonGMR4::MuonReadoutElement
The MuonReadoutElement is an abstract class representing the geometry representing the muon detector.
Definition: MuonPhaseII/MuonDetDescr/MuonReadoutGeometryR4/MuonReadoutGeometryR4/MuonReadoutElement.h:38
python.AthDsoLogger.out
out
Definition: AthDsoLogger.py:70
M_PI
#define M_PI
Definition: ActiveFraction.h:11
MuonR4::TruthSegmentMaker::m_idCylinderR
Gaudi::Property< float > m_idCylinderR
ID / ITk cylinder radius.
Definition: TruthSegmentMaker.h:151
AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_evtStore
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
Definition: AthCommonDataStore.h:390
AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_vhka
std::vector< SG::VarHandleKeyArray * > m_vhka
Definition: AthCommonDataStore.h:398
deg
#define deg
Definition: SbPolyhedron.cxx:17
MM
@ MM
Definition: RegSelEnums.h:38
xAOD::MuonSegment_v1
Class describing a MuonSegment.
Definition: MuonSegment_v1.h:33
MuonR4::TruthSegmentMaker::m_detMgr
const MuonGMR4::MuonDetectorManager * m_detMgr
Pointer to the muon readout geometry.
Definition: TruthSegmentMaker.h:155
MuonR4::TruthSegmentMaker::m_pileUpHitMinE
Gaudi::Property< float > m_pileUpHitMinE
Minimum energy threshold for pile up hits to be converted
Definition: TruthSegmentMaker.h:143
ActsTrk::IDetectorElementBase::detectorType
virtual DetectorType detectorType() const =0
Returns the detector element type.
MuonGMR4::MuonReadoutElement::layerHash
virtual IdentifierHash layerHash(const Identifier &measId) const =0
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
ATH_MSG_VERBOSE
#define ATH_MSG_VERBOSE(x)
Definition: AthMsgStreamMacros.h:28
MuonR4::TruthSegmentMaker::m_qKey
SG::WriteDecorHandleKey< xAOD::MuonSegmentContainer > m_qKey
Decoration key of the muon charge
Definition: TruthSegmentMaker.h:137
SG::ReadHandleKey
Property holding a SG store/key/clid from which a ReadHandle is made.
Definition: StoreGate/StoreGate/ReadHandleKey.h:39
ActsTrk::DetectorType::sTgc
@ sTgc
Micromegas (NSW)
read_hist_ntuple.h1
h1
Definition: read_hist_ntuple.py:21
MuonGMR4::RpcReadoutElement
Definition: MuonPhaseII/MuonDetDescr/MuonReadoutGeometryR4/MuonReadoutGeometryR4/RpcReadoutElement.h:17
MuonR4::TruthSegmentMaker::SegPars_t
xAOD::MeasVector< SegmentFit::toInt(SegmentFit::ParamDefs::nPars)> SegPars_t
Definition: TruthSegmentMaker.h:60
xAOD::MuonSegment_v1::setT0Error
void setT0Error(float t0, float t0Error)
Sets the time error.
Definition: MuonSegment_v1.cxx:51
AthenaPoolTestRead.sc
sc
Definition: AthenaPoolTestRead.py:27
MuonR4::SegmentFit::ParamDefs::phi
@ phi
MuonR4::TruthSegmentMaker::m_ptKey
SG::WriteDecorHandleKey< xAOD::MuonSegmentContainer > m_ptKey
Decoration key of the associated particle pt.
Definition: TruthSegmentMaker.h:133
AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::detStore
const ServiceHandle< StoreGateSvc > & detStore() const
The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
Definition: AthCommonDataStore.h:95
SG::VarHandleKeyArray::setOwner
virtual void setOwner(IDataHandleHolder *o)=0
IDTPMcnv.htype
htype
Definition: IDTPMcnv.py:29
MuonR4::TruthSegmentMaker::constructSegmentFromHits
xAOD::MuonSegment * constructSegmentFromHits(const EventContext &ctx, const Amg::Transform3D &locToGlob, const SimHitVec_t &hits, WriteDecorHolder &writeShip) const
Takes a list of related sim hits and transforms them into a truth segment.
Definition: TruthSegmentMaker.cxx:176
xAOD::MuonSegment_v1::setNHits
void setNHits(int nPrecisionHits, int nPhiLayers, int nTrigEtaLayers)
Set the number of hits/layers.
Definition: MuonSegment_v1.cxx:88
MuonR4::TruthSegmentMaker::m_useOnlyMuonHits
Gaudi::Property< bool > m_useOnlyMuonHits
Build segments from muon hits only.
Definition: TruthSegmentMaker.h:139
MuonR4::TruthSegmentMaker::m_locParKey
SG::WriteDecorHandleKey< xAOD::MuonSegmentContainer > m_locParKey
Decoration key of the local parameters.
Definition: TruthSegmentMaker.h:135
AthCommonDataStore::declareGaudiProperty
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition: AthCommonDataStore.h:156
AthCommonReentrantAlgorithm::extraOutputDeps
virtual const DataObjIDColl & extraOutputDeps() const override
Return the list of extra output dependencies.
Definition: AthCommonReentrantAlgorithm.cxx:89
MuonR4::TruthSegmentMaker::m_idCylinderHalfZ
Gaudi::Property< float > m_idCylinderHalfZ
ID / Itk cylinder half length.
Definition: TruthSegmentMaker.h:153
xAOD::MuonSimHit_v1::kineticEnergy
float kineticEnergy() const
Returns the kinetic energy of the traversing particle.
python.utils.AtlRunQueryDQUtils.p
p
Definition: AtlRunQueryDQUtils.py:209
AthCommonDataStore
Definition: AthCommonDataStore.h:52
Amg::toString
std::string toString(const Translation3D &translation, int precision=4)
GeoPrimitvesToStringConverter.
Definition: GeoPrimitivesToStringConverter.h:40
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
TRT::Hit::driftTime
@ driftTime
Definition: HitInfo.h:43
xAOD::MuonSimHit_v1::genParticleLink
const HepMcParticleLink & genParticleLink() const
Returns the link to the HepMC particle producing this hit.
Definition: xAODMuonSimHit_V1.cxx:68
AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::outputHandles
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override
Return this algorithm's output handles.
ActsTrk::DetectorType::Mm
@ Mm
Maybe not needed in the migration.
z
#define z
SG::get
const T * get(const ReadCondHandleKey< T > &key, const EventContext &ctx)
Convenience function to retrieve an object given a ReadCondHandleKey.
Definition: ReadCondHandle.h:287
TGC
Definition: TgcBase.h:6
h
EL::StatusCode
::StatusCode StatusCode
StatusCode definition for legacy code.
Definition: PhysicsAnalysis/D3PDTools/EventLoop/EventLoop/StatusCode.h:22
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
extractSporadic.h
list h
Definition: extractSporadic.py:96
AthCommonDataStore::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)
Definition: AthCommonDataStore.h:145
MuonR4::TruthSegmentMaker::HitLinkVec_t
std::vector< EleLink_t > HitLinkVec_t
Definition: TruthSegmentMaker.h:56
MuonGMR4::SpectrometerSector::globalToLocalTrans
Amg::Transform3D globalToLocalTrans(const ActsGeometryContext &gctx) const
Returns the global -> local transformation from the ATLAS global.
Definition: SpectrometerSector.cxx:78
AthCommonReentrantAlgorithm::sysExecute
virtual StatusCode sysExecute(const EventContext &ctx) override
Execute an algorithm.
Definition: AthCommonReentrantAlgorithm.cxx:76
Amg::Transform3D
Eigen::Affine3d Transform3D
Definition: GeoPrimitives.h:46
MuonR4::TruthSegmentMaker::m_pileUpHitDistance
Gaudi::Property< float > m_pileUpHitDistance
Maximum separation between two pile-up hits.
Definition: TruthSegmentMaker.h:149
Amg::transform
Amg::Vector3D transform(Amg::Vector3D &v, Amg::Transform3D &tr)
Transform a point from a Trasformation3D.
Definition: GeoPrimitivesHelpers.h:156
chi2
double chi2(TH1 *h0, TH1 *h1)
Definition: comparitor.cxx:525
xAOD::MuonSimHit_v1::localDirection
ConstVectorMap< 3 > localDirection() const
Returns the local direction of the traversing particle.
Definition: xAODMuonSimHit_V1.cxx:66
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
hist_file_dump.f
f
Definition: hist_file_dump.py:140
AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_detStore
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
Definition: AthCommonDataStore.h:393
SG::VarHandleKey::initialize
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
Definition: AthToolSupport/AsgDataHandles/Root/VarHandleKey.cxx:103
MuonR4::SegmentFit::ParamDefs::x0
@ x0
DataVector
Derived DataVector<T>.
Definition: DataVector.h:794
ActsGeometryContext
Include the GeoPrimitives which need to be put first.
Definition: ActsGeometryContext.h:27
MuonR4::SegmentFit::ParamDefs::time
@ time
beamspotman.dir
string dir
Definition: beamspotman.py:621
SG::VarHandleKeyArray::renounce
virtual void renounce()=0
SG::HandleClassifier::type
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
Definition: HandleClassifier.h:54
NswErrorCalibData::Input::stripId
Identifier stripId
Identifier of the strip.
Definition: NswErrorCalibData.h:27
MuonR4::SegmentFit::ParamDefs::y0
@ y0
HepMC::ConstGenParticlePtr
const GenParticle * ConstGenParticlePtr
Definition: GenParticle.h:38
MuonR4::TruthSegmentMaker::m_pileUpHitELoss
Gaudi::Property< float > m_pileUpHitELoss
Maximum energy loss between two pile-up hits.
Definition: TruthSegmentMaker.h:145
ActsTrk::DetectorType::Mdt
@ Mdt
MuonSpectrometer.
RPC
@ RPC
Definition: RegSelEnums.h:32
merge_scale_histograms.doc
string doc
Definition: merge_scale_histograms.py:9
id
SG::auxid_t id
Definition: Control/AthContainers/Root/debug.cxx:239
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:240
plotBeamSpotMon.b
b
Definition: plotBeamSpotMon.py:76
MuonR4::SegmentFit::toInt
constexpr int toInt(const ParamDefs p)
Definition: MuonHoughDefs.h:42
charge
double charge(const T &p)
Definition: AtlasPID.h:986
DataVector::push_back
value_type push_back(value_type pElem)
Add an element to the end of the collection.
SG::CondHandleKey::initialize
StatusCode initialize(bool used=true)
Amg::Vector3D
Eigen::Matrix< double, 3, 1 > Vector3D
Definition: GeoPrimitives.h:47
THROW_EXCEPTION
#define THROW_EXCEPTION(MESSAGE)
Definition: throwExcept.h:10
python.LumiBlobConversion.pos
pos
Definition: LumiBlobConversion.py:16
SG::WriteHandle
Definition: StoreGate/StoreGate/WriteHandle.h:73
a
TList * a
Definition: liststreamerinfos.cxx:10
SG::WriteDecorHandleKey::initialize
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
h
MuonR4::TruthSegmentMaker::m_includePileUpHits
Gaudi::Property< bool > m_includePileUpHits
Construct segments from pile-up hits without GenParticleLink.
Definition: TruthSegmentMaker.h:141
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
MuonR4::TruthSegmentMaker::toChamber
Amg::Transform3D toChamber(const ActsGeometryContext &gctx, const Identifier &chanId) const
Returns the transform from the local simHit frame -> chamber frame.
Definition: TruthSegmentMaker.cxx:101
MuonCalib::MdtCalibDataContainer
Definition: MdtCalibDataContainer.h:20
xAOD::MuonSegment_v1::setIdentifier
void setIdentifier(int sector, ::Muon::MuonStationIndex::ChIndex chamberIndex, int etaIndex, ::Muon::MuonStationIndex::TechnologyIndex technology)
Set the identifier.
Definition: MuonSegment_v1.cxx:73
re
const boost::regex re(r_e)
MuonR4::TruthSegmentMaker::m_readKeys
SG::ReadHandleKeyArray< xAOD::MuonSimHitContainer > m_readKeys
List of sim hit containers from which the truth segments shall be retrieved.
Definition: TruthSegmentMaker.h:125
MuonGMR4::MuonReadoutElement::measurementHash
virtual IdentifierHash measurementHash(const Identifier &measId) const =0
Constructs the identifier hash from the full measurement Identifier.
xAOD::MuonSimHit_v1::localPosition
ConstVectorMap< 3 > localPosition() const
Returns the local postion of the traversing particle.
Definition: xAODMuonSimHit_V1.cxx:60
ActsTrk::DetectorType::Rpc
@ Rpc
Monitored Drift Tubes.
LArNewCalib_DelayDump_OFC_Cali.idx
idx
Definition: LArNewCalib_DelayDump_OFC_Cali.py:69
MuonR4::TruthSegmentMaker::m_pileUpHitAngleCone
Gaudi::Property< float > m_pileUpHitAngleCone
Maximum scattering angle between two pile-up hits.
Definition: TruthSegmentMaker.h:147
xAOD::MuonSegment_v1::setDirection
void setDirection(float px, float py, float pz)
Sets the direction.
Definition: MuonSegment_v1.cxx:39
SG::VarHandleBase::vhKey
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Definition: StoreGate/src/VarHandleBase.cxx:629
MuonGMR4::MuonReadoutElement::localToGlobalTrans
const Amg::Transform3D & localToGlobalTrans(const ActsGeometryContext &ctx) const
Returns the local to global transformation into the ATLAS coordinate system.
Definition: MuonPhaseII/MuonDetDescr/MuonReadoutGeometryR4/src/MuonReadoutElement.cxx:81
MuonR4::TruthSegmentMaker::HitPosTuple_t
std::tuple< const xAOD::MuonSimHit *, Amg::Vector3D, Amg::Vector3D > HitPosTuple_t
Tuple consisting out of pointer to the sim hit and the position & direction expressed in the chamber'...
Definition: TruthSegmentMaker.h:50
Amg::signedDistance
double signedDistance(const Amg::Vector3D &posA, const Amg::Vector3D &dirA, const Amg::Vector3D &posB, const Amg::Vector3D &dirB)
Calculates the signed distance between two lines in 3D space.
Definition: GeoPrimitivesHelpers.h:329
AthCommonReentrantAlgorithm< Gaudi::Algorithm >::sysInitialize
virtual StatusCode sysInitialize() override
Override sysInitialize.
Definition: AthCommonReentrantAlgorithm.cxx:107
MuonR4::TruthSegmentMaker::m_segmentKey
SG::WriteHandleKey< xAOD::MuonSegmentContainer > m_segmentKey
Key under which the segment Container will be recorded in StoreGate.
Definition: TruthSegmentMaker.h:129
python.Bindings.keys
keys
Definition: Control/AthenaPython/python/Bindings.py:801
MuonR4::TruthSegmentMaker::m_idHelperSvc
ServiceHandle< Muon::IMuonIdHelperSvc > m_idHelperSvc
IdHelperSvc to decode the Identifiers.
Definition: TruthSegmentMaker.h:123
MuonR4::TruthSegmentMaker::m_nswUncertKey
SG::ReadCondHandleKey< NswErrorCalibData > m_nswUncertKey
Data dependency on the Nsw calibration container to estimate the uncertaintys.
Definition: TruthSegmentMaker.h:160
LHEF::Writer
Pythia8::Writer Writer
Definition: Prophecy4fMerger.cxx:12
NswErrorCalibData::Input
Helper struct to be parsed to the object to derive the specific error of the cluster.
Definition: NswErrorCalibData.h:25
MuonR4::TruthSegmentMaker::m_mdtCalibKey
SG::ReadCondHandleKey< MuonCalib::MdtCalibDataContainer > m_mdtCalibKey
Data dependency on the Mdt calibration container to calculate the uncertainty.
Definition: TruthSegmentMaker.h:157
drawFromPickle.sin
sin
Definition: drawFromPickle.py:36
MuonR4::TruthSegmentMaker::hitUncertainty
float hitUncertainty(const EventContext &ctx, const xAOD::MuonSimHit &hit) const
Evaluates the hit uncertainty.
Definition: TruthSegmentMaker.cxx:52
AthCommonReentrantAlgorithm< Gaudi::Algorithm >::m_extendedExtraObjects
DataObjIDColl m_extendedExtraObjects
Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
Definition: AthCommonReentrantAlgorithm.h:114
pow
constexpr int pow(int base, int exp) noexcept
Definition: ap_fixedTest.cxx:15
Amg::distance
float distance(const Amg::Vector3D &p1, const Amg::Vector3D &p2)
calculates the distance between two point in 3D space
Definition: GeoPrimitivesHelpers.h:54
xAOD::MuonSegment_v1::setFitQuality
void setFitQuality(float chiSquared, float numberDoF)
Set the 'Fit Quality' information.
Definition: MuonSegment_v1.cxx:61
MuonR4::TruthSegmentMaker::muonPt
float muonPt(const xAOD::MuonSimHit &hit, const Amg::Vector3D &globDir) const
Returns the muon pt from the sim hit.
Definition: TruthSegmentMaker.cxx:93
MuonGMR4::TgcReadoutElement
Definition: MuonPhaseII/MuonDetDescr/MuonReadoutGeometryR4/MuonReadoutGeometryR4/TgcReadoutElement.h:19
xAOD::MuonSimHit_v1::mass
float mass() const
Returns the rest-mass of the traversing particle.
MuonR4::SegmentFit::ParamDefs::theta
@ theta
MDT
@ MDT
Definition: RegSelEnums.h:31
fitman.k
k
Definition: fitman.py:528
python.SystemOfUnits.m
float m
Definition: SystemOfUnits.py:106
generate::Zero
void Zero(TH1D *hin)
Definition: generate.cxx:32
MuonR4::TruthSegmentMaker::EleLink_t
ElementLink< xAOD::MuonSimHitContainer > EleLink_t
Definition: TruthSegmentMaker.h:55
ServiceHandle< ICondSvc >
MuonGMR4::MuonDetectorManager::getReadoutElement
const MuonReadoutElement * getReadoutElement(const Identifier &id) const
Returns a generic Muon readout element.
mapkey::key
key
Definition: TElectronEfficiencyCorrectionTool.cxx:37
isMuon
bool isMuon(const T &p)
Definition: AtlasPID.h:202
Identifier
Definition: IdentifierFieldParser.cxx:14