ATLAS Offline Software
Loading...
Searching...
No Matches
MuonR4::sTgcFastDigiTool Class Referencefinal

#include <sTgcFastDigiTool.h>

Inheritance diagram for MuonR4::sTgcFastDigiTool:
Collaboration diagram for MuonR4::sTgcFastDigiTool:

Public Member Functions

StatusCode initialize () override final
StatusCode finalize () override final
StatusCode processBunchXing (int bunchXing, SubEventIterator bSubEvents, SubEventIterator eSubEvents) override final
StatusCode mergeEvent (const EventContext &ctx) override final
StatusCode prepareEvent (const EventContext &ctx, const unsigned int) override final
 When being run from PileUpToolsAlgs, this method is called at the start of the subevts loop.
StatusCode processAllSubEvents (const EventContext &ctx) override final
 alternative interface which uses the PileUpMergeSvc to obtain all the required SubEvents.
StatusCode processAllSubEvents (const EventContext &ctx) const
 Reentrant version of the digitization tool.
 PileUpToolBase (const std::string &type, const std::string &name, const IInterface *parent)

Protected Types

using TimedHit = TimedHitPtr<xAOD::MuonSimHit>
using TimedHits = std::vector<TimedHitPtr<xAOD::MuonSimHit>>
template<class DetType>
using OutDigitCache_t = std::vector<std::unique_ptr<DetType>>
 DigitContainers are sorted by DigitCollections which are the ensemble of all hits in a given MuonChamber.
using DeadTimeMap = std::unordered_map<Identifier, double>

Protected Member Functions

StatusCode digitize (const EventContext &ctx, const TimedHits &hitsToDigit, xAOD::MuonSimHitContainer *sdoContainer) const override final
 Digitize the time ordered hits and write them to the digit format specific for the detector technology.
CLHEP::HepRandomEngine * getRandomEngine (const EventContext &ctx) const
xAOD::MuonSimHitaddSDO (const TimedHit &hit, xAOD::MuonSimHitContainer *sdoContainer) const
 Adds the timed simHit to the output SDO container.
const ActsTrk::GeometryContextgetGeoCtx (const EventContext &ctx) const
 Returns the reference to the ActsTrk::GeometryContext needed to fetch global positions from the Readout geometry.
template<class DigitColl>
DigitColl * fetchCollection (const Identifier &hitId, OutDigitCache_t< DigitColl > &digitCache) const
 Helper function that provides fetches the proper DigitCollection from the DigitCache for a given hit identifier If the Collection is fetched for the first time, it's inserted into the cache first.
template<class DigitCont, class DigitColl>
StatusCode writeDigitContainer (const EventContext &ctx, const SG::WriteHandleKey< DigitCont > &key, OutDigitCache_t< DigitColl > &&digitCache, unsigned int hashMax) const
 Helper function to move the collected digits into the final DigitContainer.

Static Protected Member Functions

static double hitTime (const TimedHit &hit)
 Returns the global time of the hit which is the sum of eventTime & individual hit time.
static bool passDeadTime (const Identifier &channelId, const double hitTime, const double deadTimeWindow, DeadTimeMap &deadTimeMap)
 Returns whether the new digit is within the dead time window.

Protected Attributes

const MuonGMR4::MuonDetectorManagerm_detMgr {nullptr}
ServiceHandle< Muon::IMuonIdHelperSvcm_idHelperSvc

Private Types

using DigiCache = OutDigitCache_t<sTgcDigitCollection>
using PileUpHits = PileUpMergeSvc::TimedList<xAOD::MuonSimHitContainer>::type
using SimHitLocalCopy
 Create a local copy of the sim hits to ensure overlayed hits across the events remain valid.

Private Member Functions

bool digitizeStrip (const EventContext &ctx, const TimedHit &timedHit, const NswErrorCalibData *errorCalibDB, const Muon::DigitEffiData *effiData, CLHEP::HepRandomEngine *rndEngine, sTgcDigitCollection &outCollection) const
bool digitizeWire (const EventContext &ctx, const TimedHit &timedHit, const Muon::DigitEffiData *effiData, CLHEP::HepRandomEngine *rndEngine, sTgcDigitCollection &outCollection) const
bool digitizePad (const EventContext &ctx, const TimedHit &timedHit, const Muon::DigitEffiData *effiData, CLHEP::HepRandomEngine *rndEngine, sTgcDigitCollection &outCollection) const
int associateBCIdTag (const EventContext &ctx, const TimedHit &timedHit) const
 : Associates the global bcIdTag to the digit
StatusCode fillTimedHits (PileUpHits &&hitColl, TimedHits &timedHits) const
 Translates the PileUpHits into the timed hits format.

Private Attributes

SG::WriteHandleKey< sTgcDigitContainerm_writeKey {this, "OutputObjectName", "sTGC_DIGITS"}
SG::ReadCondHandleKey< Muon::DigitEffiDatam_effiDataKey
SG::ReadCondHandleKey< NswErrorCalibDatam_uncertCalibKey
std::array< std::atomic< unsigned >, 3 > m_allHits ATLAS_THREAD_SAFE {}
std::array< std::atomic< unsigned >, 3 > m_acceptedHits ATLAS_THREAD_SAFE {}
Gaudi::Property< bool > m_digitizeStrip {this, "doStrips", true, "Digitize strip hits"}
Gaudi::Property< bool > m_digitizeWire {this, "doWires", true, "Digitize wire hits"}
Gaudi::Property< bool > m_digitizePads {this, "doPads", true, "Digitize pad hits"}
Gaudi::Property< double > m_deadTime {this, "deadTime", 300. * Gaudi::Units::ns}
Gaudi::Property< bool > m_digitizeMuonOnly
Gaudi::Property< double > m_energyDepositThreshold {this,"energyDepositThreshold",300.0*CLHEP::eV}
Gaudi::Property< double > m_limitElectronKineticEnergy {this,"limitElectronKineticEnergy",5.0*CLHEP::MeV}
SG::ReadHandleKey< xAOD::MuonSimHitContainerm_simHitKey {this, "SimHitKey", ""}
SG::ReadHandleKey< ActsTrk::GeometryContextm_geoCtxKey
ServiceHandle< PileUpMergeSvcm_mergeSvc {this, "PileUpMergeSvc", "PileUpMergeSvc", ""}
Gaudi::Property< std::string > m_streamName {this, "StreamName", ""}
ServiceHandle< IAthRNGSvcm_rndmSvc {this, "RndmSvc", "AthRNGSvc", ""}
SG::WriteHandleKey< xAOD::MuonSimHitContainerm_sdoKey {this, "OutputSDOName", ""}
Gaudi::Property< bool > m_onlyUseContainerName
Gaudi::Property< bool > m_includePileUpTruth {this, "IncludePileUpTruth", true, "Include pile-up truth info"}
std::string m_inputObjectName {""}
TimedHits m_timedHits {}
std::vector< SimHitLocalCopym_simHits {}

structors and AlgTool implementation

virtual bool toProcess (int bunchXing) const override
 the method this base class helps implementing
virtual bool filterPassed () const override
 dummy implementation of passing filter
virtual void resetFilter () override
 dummy implementation of filter reset
Gaudi::Property< int > m_firstXing
Gaudi::Property< int > m_lastXing
Gaudi::Property< int > m_vetoPileUpTruthLinks
bool m_filterPassed {true}

Detailed Description

Definition at line 14 of file sTgcFastDigiTool.h.

Member Typedef Documentation

◆ DeadTimeMap

using MuonR4::MuonDigitizationTool::DeadTimeMap = std::unordered_map<Identifier, double>
protectedinherited

Definition at line 100 of file MuonDigitizationTool.h.

◆ DigiCache

◆ OutDigitCache_t

template<class DetType>
using MuonR4::MuonDigitizationTool::OutDigitCache_t = std::vector<std::unique_ptr<DetType>>
protectedinherited

DigitContainers are sorted by DigitCollections which are the ensemble of all hits in a given MuonChamber.

To fill the final DigitContainer thread-safely, the DigitCollections shall be cached pre cached in a OutDigitCache_t vector which is later moved to the final DigitContainer

Definition at line 80 of file MuonDigitizationTool.h.

◆ PileUpHits

◆ SimHitLocalCopy

Initial value:
std::pair<std::unique_ptr<xAOD::MuonSimHitContainer>,
std::unique_ptr<xAOD::MuonSimHitAuxContainer>>

Create a local copy of the sim hits to ensure overlayed hits across the events remain valid.

Definition at line 143 of file MuonDigitizationTool.h.

◆ TimedHit

Definition at line 57 of file MuonDigitizationTool.h.

◆ TimedHits

using MuonR4::MuonDigitizationTool::TimedHits = std::vector<TimedHitPtr<xAOD::MuonSimHit>>
protectedinherited

Definition at line 58 of file MuonDigitizationTool.h.

Member Function Documentation

◆ addSDO()

xAOD::MuonSimHit * MuonR4::MuonDigitizationTool::addSDO ( const TimedHit & hit,
xAOD::MuonSimHitContainer * sdoContainer ) const
protectedinherited

Adds the timed simHit to the output SDO container.

The hit may be rejected if it's originating from pile-up and the pile-up truth skimming strategy is applied

Definition at line 142 of file MuonDigitizationTool.cxx.

143 {
144 if(!sdoContainer) {
145 ATH_MSG_VERBOSE("No SDO container setup of writing");
146 return nullptr;
147 }
148 if (!m_includePileUpTruth && HepMC::ignoreTruthLink(hit->genParticleLink(), m_vetoPileUpTruthLinks)) {
149 ATH_MSG_VERBOSE("Hit "<<m_idHelperSvc->toString(hit->identify())<<" is a pile-up truth link");
150 return nullptr;
151 }
152 ATH_MSG_VERBOSE(m_idHelperSvc->toString(hit->identify())<<", pdgID: "<<hit->pdgId()<< " genParticleLink :"<<hit->genParticleLink());
153 ATH_MSG_VERBOSE("Genparticle: "<<hit->genParticleLink());
154 xAOD::MuonSimHit* sdoHit = sdoContainer->push_back(std::make_unique<xAOD::MuonSimHit>());
155 (*sdoHit) = (*hit);
156 static const SG::Accessor<float> acc_eventTime{"MuSim_evtTime"};
157 static const SG::Accessor<unsigned short> acc_eventID{"MuSim_evtID"};
158 static const SG::Accessor<unsigned short> acc_puType{"MuSim_puType"};
159 acc_eventTime(*sdoHit) = hit.eventTime();
160 acc_eventID(*sdoHit) = hit.eventId();
161 acc_puType(*sdoHit) = hit.pileupType();
162 return sdoHit;
163 }
#define ATH_MSG_VERBOSE(x)
value_type push_back(value_type pElem)
Add an element to the end of the collection.
ServiceHandle< Muon::IMuonIdHelperSvc > m_idHelperSvc
Gaudi::Property< bool > m_includePileUpTruth
Gaudi::Property< int > m_vetoPileUpTruthLinks
bool ignoreTruthLink(const T &p, bool vetoPileUp)
Helper function for SDO creation in PileUpTools.
MuonSimHit_v1 MuonSimHit
Defined the version of the MuonSimHit.
Definition MuonSimHit.h:12

◆ associateBCIdTag()

int MuonR4::sTgcFastDigiTool::associateBCIdTag ( const EventContext & ctx,
const TimedHit & timedHit ) const
private

: Associates the global bcIdTag to the digit

Parameters
ctxEventContext
timedHitHit to calculate the time from

To be implemented

Definition at line 294 of file sTgcFastDigiTool.cxx.

295 {
297 return 0;
298 }

◆ digitize()

StatusCode MuonR4::sTgcFastDigiTool::digitize ( const EventContext & ctx,
const TimedHits & hitsToDigit,
xAOD::MuonSimHitContainer * sdoContainer ) const
finaloverrideprotectedvirtual

Digitize the time ordered hits and write them to the digit format specific for the detector technology.

A new MuonSimHitContainer pointer is parsed to also create the MuonSDO. If a new SDO should be added to the container plese use the addSDO() method as defined below.

Fetch the conditions for efficiency calculations

ignore radiation for now

Write everything at the end into the final digit container

Implements MuonR4::MuonDigitizationTool.

Definition at line 40 of file sTgcFastDigiTool.cxx.

42 {
43 const sTgcIdHelper& idHelper{m_idHelperSvc->stgcIdHelper()};
44 // Prepare the temporary cache
45 DigiCache digitCache{};
47 const Muon::DigitEffiData* efficiencyMap{nullptr};
48 ATH_CHECK(SG::get(efficiencyMap, m_effiDataKey, ctx));
49 const NswErrorCalibData* nswUncertDB{nullptr};
50 ATH_CHECK(SG::get(nswUncertDB, m_uncertCalibKey, ctx));
51
52 CLHEP::HepRandomEngine* rndEngine = getRandomEngine(ctx);
53 xAOD::ChamberViewer viewer{hitsToDigit, m_idHelperSvc.get()};
54 do {
55 for (const TimedHit& simHit : viewer) {
57 if (m_digitizeMuonOnly && !MC::isMuon(simHit)){
58 continue;
59 }
60
61 if (simHit->energyDeposit() < m_energyDepositThreshold){
62 ATH_MSG_VERBOSE("Hit with Energy Deposit of " << simHit->energyDeposit()
63 << " less than " << m_energyDepositThreshold << ". Skip this hit." );
64 continue;
65 }
66
67 const double hitKineticEnergy = simHit->kineticEnergy();
68 if (hitKineticEnergy < m_limitElectronKineticEnergy && MC::isElectron(simHit)) {
69 ATH_MSG_DEBUG("Skip electron hit with kinetic energy " << hitKineticEnergy
70 << ", which is less than the lower limit of " << m_limitElectronKineticEnergy);
71 continue;
72 }
73 sTgcDigitCollection* digiColl = fetchCollection(simHit->identify(), digitCache);
74 const bool digitizedPad = digitizePad(ctx, simHit, efficiencyMap, rndEngine, *digiColl);
75 const std::int16_t padChannel = digitizedPad ? idHelper.channel(digiColl->back()->identify()) : -1;
76 const bool digitizedWire = digitizeWire(ctx, simHit, efficiencyMap, rndEngine, *digiColl);
77 const std::int16_t wireChannel = digitizedWire ? idHelper.channel(digiColl->back()->identify()) : -1;
78 const bool digitizedStrip = digitizeStrip(ctx, simHit, nswUncertDB, efficiencyMap, rndEngine, *digiColl);
79 const std::int16_t stripCh = digitizedStrip ? idHelper.channel(digiColl->back()->identify()) : -1;
80
81 if (digitizedStrip || digitizedPad || digitizedWire) {
82 xAOD::MuonSimHit* sdo = addSDO(simHit, sdoContainer);
83 ChVec_t& stripChV{dec_stripCh(*sdo)};
84 ChVec_t& wireChV{dec_wireCh(*sdo)};
85 ChVec_t& padChV{dec_padCh(*sdo)};
86 if (stripCh > 0) { stripChV.push_back(stripCh); }
87 if (wireChannel > 0) { wireChV.push_back(wireChannel); }
88 if (padChannel > 0) { padChV.push_back(padChannel); }
89 sdo->setIdentifier(digiColl->back()->identify());
90 }
91 }
92 } while (viewer.next());
94 ATH_CHECK(writeDigitContainer(ctx, m_writeKey, std::move(digitCache),
95 idHelper.module_hash_max()));
96 return StatusCode::SUCCESS;
97 }
#define ATH_CHECK
Evaluate an expression and check for errors.
#define ATH_MSG_DEBUG(x)
const T * back() const
Access the last element in the collection as an rvalue.
Identifier identify() const
Definition MuonDigit.h:30
size_type module_hash_max() const
the maximum hash value
xAOD::MuonSimHit * addSDO(const TimedHit &hit, xAOD::MuonSimHitContainer *sdoContainer) const
Adds the timed simHit to the output SDO container.
CLHEP::HepRandomEngine * getRandomEngine(const EventContext &ctx) const
TimedHitPtr< xAOD::MuonSimHit > TimedHit
DigitColl * fetchCollection(const Identifier &hitId, OutDigitCache_t< DigitColl > &digitCache) const
Helper function that provides fetches the proper DigitCollection from the DigitCache for a given hit ...
StatusCode writeDigitContainer(const EventContext &ctx, const SG::WriteHandleKey< DigitCont > &key, OutDigitCache_t< DigitColl > &&digitCache, unsigned int hashMax) const
Helper function to move the collected digits into the final DigitContainer.
SG::ReadCondHandleKey< Muon::DigitEffiData > m_effiDataKey
Gaudi::Property< double > m_limitElectronKineticEnergy
bool digitizeStrip(const EventContext &ctx, const TimedHit &timedHit, const NswErrorCalibData *errorCalibDB, const Muon::DigitEffiData *effiData, CLHEP::HepRandomEngine *rndEngine, sTgcDigitCollection &outCollection) const
Gaudi::Property< double > m_energyDepositThreshold
Gaudi::Property< bool > m_digitizeMuonOnly
SG::WriteHandleKey< sTgcDigitContainer > m_writeKey
SG::ReadCondHandleKey< NswErrorCalibData > m_uncertCalibKey
bool digitizePad(const EventContext &ctx, const TimedHit &timedHit, const Muon::DigitEffiData *effiData, CLHEP::HepRandomEngine *rndEngine, sTgcDigitCollection &outCollection) const
OutDigitCache_t< sTgcDigitCollection > DigiCache
bool digitizeWire(const EventContext &ctx, const TimedHit &timedHit, const Muon::DigitEffiData *effiData, CLHEP::HepRandomEngine *rndEngine, sTgcDigitCollection &outCollection) const
int channel(const Identifier &id) const override
void setIdentifier(const Identifier &id)
Sets the global ATLAS identifier.
bool isElectron(const T &p)
bool isMuon(const T &p)
const T * get(const ReadCondHandleKey< T > &key, const EventContext &ctx)
Convenience function to retrieve an object given a ReadCondHandleKey.

◆ digitizePad()

bool MuonR4::sTgcFastDigiTool::digitizePad ( const EventContext & ctx,
const TimedHit & timedHit,
const Muon::DigitEffiData * effiData,
CLHEP::HepRandomEngine * rndEngine,
sTgcDigitCollection & outCollection ) const
private

Check efficiencies

Definition at line 300 of file sTgcFastDigiTool.cxx.

304 {
305
306 if (!m_digitizePads) {
307 return false;
308 }
309
310 ++m_allHits[channelType::Pad];
311
312 const Identifier hitId{timedHit->identify()};
313 const sTgcIdHelper& idHelper{m_idHelperSvc->stgcIdHelper()};
314 const MuonGMR4::sTgcReadoutElement* readOutEle = m_detMgr->getsTgcReadoutElement(hitId);
315 const int gasGap = idHelper.gasGap(hitId);
316
317 const IdentifierHash hitHash = readOutEle->createHash(gasGap, sTgcIdHelper::sTgcChannelTypes::Pad, 1);
318
319
320 const Amg::Vector2D padPos{readOutEle->stripLayer(hitHash).to2D(xAOD::toEigen(timedHit->localPosition()), true)};
322 const MuonGMR4::PadDesign& design{readOutEle->padDesign(hitHash)};
323
324 const auto [padEta, padPhi] = design.channelNumber(padPos);
325 if (padEta <= 0 || padPhi <= 0) {
326 ATH_MSG_VERBOSE("The pad "<<Amg::toString(padPos)<<" in "<<m_idHelperSvc->toStringGasGap(hitId)
327 <<" is outside of the acceptance of "<<std::endl<<design);
328 return false;
329 }
330 bool isValid{false};
331 const Identifier padId = idHelper.padID(hitId, readOutEle->multilayer(),
332 gasGap, channelType::Pad, padEta, padPhi, isValid);
333
334
335 if (!isValid) {
336 ATH_MSG_WARNING("Failed to deduce a valid pad Identifier from "<<Amg::toString(padPos)
337 <<" in "<<m_idHelperSvc->toStringGasGap(hitId));
338 return false;
339 }
340
342 bool isInnerQ1 = readOutEle->isEtaZero(readOutEle->measurementHash(hitId), xAOD::toEigen(timedHit->localPosition()).block<2,1>(0,0));
343 if (efficiencyMap && efficiencyMap->getEfficiency(hitId, isInnerQ1) < CLHEP::RandFlat::shoot(rndEngine,0.,1.)){
344 ATH_MSG_VERBOSE("Simulated pad hit "<<xAOD::toEigen(timedHit->localPosition())
345 << m_idHelperSvc->toString(hitId) <<" is rejected because of efficency modelling");
346 return false;
347 }
348
349 outCollection.push_back(std::make_unique<sTgcDigit>(padId,
350 associateBCIdTag(ctx, timedHit),
351 hitTime(timedHit), 666, false, false));
352
353 ++m_acceptedHits[channelType::Pad];
354
355 return true;
356 }
#define ATH_MSG_WARNING(x)
bool isValid(const T &p)
Av: we implement here an ATLAS-sepcific convention: all particles which are 99xxxxx are fine.
Definition AtlasPID.h:878
std::pair< int, int > channelNumber(const Amg::Vector2D &hitPos) const
Function to retrieve the pad eta and phi given a local position coordinate.
Amg::Vector2D to2D(const Amg::Vector3D &vec, const bool phiView) const
Transforms a 3D vector from the strip design into a 2D vector.
IdentifierHash measurementHash(const Identifier &measId) const override final
Constructs the identifier hash from the full measurement Identifier.
const PadDesign & padDesign(const IdentifierHash &measHash) const
Retrieves the readoutElement Layer given the Identifier/Hash.
int multilayer() const
Returns the multilayer of the sTgcReadoutElement.
bool isEtaZero(const IdentifierHash &measurementHash, const Amg::Vector2D &localPosition) const
static IdentifierHash createHash(const unsigned gasGap, const unsigned channelType, const unsigned channel, const unsigned wireInGrp=0)
Create a measurement hash from the Identifier fields.
const StripLayer & stripLayer(const IdentifierHash &measId) const
const MuonGMR4::MuonDetectorManager * m_detMgr
static double hitTime(const TimedHit &hit)
Returns the global time of the hit which is the sum of eventTime & individual hit time.
Gaudi::Property< bool > m_digitizePads
int associateBCIdTag(const EventContext &ctx, const TimedHit &timedHit) const
: Associates the global bcIdTag to the digit
Identifier padID(int stationName, int stationEta, int stationPhi, int multilayer, int gasGap, int channelType, int padEta, int padPhi) const
int gasGap(const Identifier &id) const override
get the hashes
std::string toString(const Translation3D &translation, int precision=4)
GeoPrimitvesToStringConverter.
Eigen::Matrix< double, 2, 1 > Vector2D

◆ digitizeStrip()

bool MuonR4::sTgcFastDigiTool::digitizeStrip ( const EventContext & ctx,
const TimedHit & timedHit,
const NswErrorCalibData * errorCalibDB,
const Muon::DigitEffiData * effiData,
CLHEP::HepRandomEngine * rndEngine,
sTgcDigitCollection & outCollection ) const
private

Check efficiencies

We're using the NSW uncertainty DB to smear the truth-hit positions using the best known uncertainties. In the process of digit -> RDO -> PRD, the smeared hit position is gonna be lost However, recall that the simplest way to generate a sTgc prd strip prd is the clustering by means of center of gravity. On average, a muon passage makes 3-4 strips to fire. Assuming 3 strips, the smeared position can be written as

        x = w_{1}*(c - p) + w_{2}*c + w_{3}*(c + p),

where c is the strip centre, p the pitch and w_{i} are the induced charge fractions on each strip

      1 =  w_{1} + w_{2} + w_{3}
      x = (w_{1} + w_{2} + w_{3}) * c + (w_{3} - w_{1})*p
         --> x = c + (w_{3} - w_{1}) *p
         --> (x-c) / p = (w_{3} - w_{1})  
         --> (x-c) / p + w_{1} = w_{3}

            1 =  2 * w_{1} + w_{2} + (x-c) / p
            1 - (x-c) /p - 2*w_{1} = w_{2}

            0< w_{1} < 0.5 - (x-c) / 2p  

Definition at line 98 of file sTgcFastDigiTool.cxx.

103 {
104
105 if (!m_digitizeStrip) {
106 return false;
107 }
108 ++m_allHits[channelType::Strip];
109 const Identifier hitId{timedHit->identify()};
110 const MuonGMR4::sTgcReadoutElement* readOutEle{m_detMgr->getsTgcReadoutElement(hitId)};
111
112 const sTgcIdHelper& idHelper{m_idHelperSvc->stgcIdHelper()};
113
114 const IdentifierHash meashHash = readOutEle->measurementHash(hitId);
115 const MuonGMR4::StripDesign& design{readOutEle->stripDesign(meashHash)};
116
117 const Amg::Vector2D stripPos = readOutEle->stripLayer(meashHash)
118 .to2D(xAOD::toEigen(timedHit->localPosition()), false);
119
120 const int stripNum = design.stripNumber(stripPos);
121 if (stripNum < 0) {
122 ATH_MSG_VERBOSE("Strip hit "<<Amg::toString(stripPos)<<" "<<m_idHelperSvc->toStringGasGap(hitId)
123 <<" is out of range "<<std::endl<<design);
124 return false;
125 }
126
127
128 bool isValid{false};
129 const int gasGap = idHelper.gasGap(hitId);
130 const Identifier stripId = idHelper.channelID(hitId, readOutEle->multilayer(),
131 gasGap, channelType::Strip, stripNum, isValid);
132
133 if (!isValid) {
134 ATH_MSG_WARNING("Failed to deduce a valid identifier from "
135 <<m_idHelperSvc->toStringGasGap(hitId)<<" strip: "<<stripNum);
136 return false;
137 }
138
140 bool isInnerQ1 = readOutEle->isEtaZero(readOutEle->measurementHash(hitId), stripPos);
141 if (efficiencyMap && efficiencyMap->getEfficiency(hitId, isInnerQ1) < CLHEP::RandFlat::shoot(rndEngine,0.,1.)){
142 ATH_MSG_VERBOSE("Simulated strip hit "<<xAOD::toEigen(timedHit->localPosition())
143 << m_idHelperSvc->toString(hitId) <<" is rejected because of efficency modelling");
144 return false;
145 }
146
147 NswErrorCalibData::Input errorCalibInput{};
148 errorCalibInput.stripId= stripId;
149 errorCalibInput.locTheta = M_PI - timedHit->localDirection().theta();
150 errorCalibInput.clusterAuthor = 3; // centroid
151
152 const double uncert = errorCalibDB->clusterUncertainty(errorCalibInput);
153 const double smearedX = CLHEP::RandGaussZiggurat::shoot(rndEngine, stripPos.x(), uncert);
154
155 const Amg::Vector2D digitPos{smearedX * Amg::Vector2D::UnitX()};
156
157 const int digitStrip = design.stripNumber(digitPos);
158 if (digitStrip < 0) {
159 ATH_MSG_VERBOSE("Smeared strip hit "<<Amg::toString(digitPos)<<" "<<m_idHelperSvc->toStringGasGap(hitId)
160 <<" is out of range "<<std::endl<<design);
161 return false;
162 }
163 const Identifier digitId = idHelper.channelID(hitId, readOutEle->multilayer(),
164 gasGap, channelType::Strip, digitStrip, isValid);
165
166 if (!isValid) {
167 ATH_MSG_WARNING("Failed to deduce a valid identifier from "
168 <<m_idHelperSvc->toStringGasGap(hitId)<<" digit: "<<digitStrip);
169 return false;
170 }
171 constexpr double dummyCharge = 66666;
192 const double pull = (smearedX - (*design.center(digitStrip)).x()) / design.stripPitch();
193 const double w1 = CLHEP::RandFlat::shoot(rndEngine, 0., 0.5 *(1. - pull));
194 const double w2 = 1. - pull -2.*w1;
195 const double w3 = pull + w1;
196
197 if (digitStrip> 1) {
198 const Identifier stripIdB = idHelper.channelID(hitId, readOutEle->multilayer(),
199 gasGap, channelType::Strip, digitStrip -1);
200 outCollection.push_back(std::make_unique<sTgcDigit>(stripIdB,
201 associateBCIdTag(ctx, timedHit),
202 hitTime(timedHit), dummyCharge * w1, false, false));
203 }
204 outCollection.push_back(std::make_unique<sTgcDigit>(digitId,
205 associateBCIdTag(ctx, timedHit),
206 hitTime(timedHit), dummyCharge * w2, false, false));
207
208 if (digitStrip + 1 <= design.numStrips()) {
209 const Identifier stripIdA = idHelper.channelID(hitId, readOutEle->multilayer(),
210 gasGap, channelType::Strip, digitStrip + 1);
211
212 outCollection.push_back(std::make_unique<sTgcDigit>(stripIdA,
213 associateBCIdTag(ctx, timedHit),
214 hitTime(timedHit), dummyCharge * w3, false, false));
215 }
216 ++m_acceptedHits[channelType::Strip];
217 return true;
218 }
#define M_PI
#define x
CheckVector2D center(int stripNumb) const
Returns the bisector of the strip (Global numbering scheme)
double stripPitch() const
Distance between two adjacent strips.
virtual int stripNumber(const Amg::Vector2D &pos) const
Calculates the number of the strip whose center is closest to the given point.
virtual int numStrips() const
Number of strips on the panel.
const StripDesign & stripDesign(const IdentifierHash &measHash) const
Retrieves the readoutElement Layer given the Identifier/Hash.
Gaudi::Property< bool > m_digitizeStrip
double clusterUncertainty(const Input &clustInfo) const
Identifier channelID(int stationName, int stationEta, int stationPhi, int multilayer, int gasGap, int channelType, int channel) const
double locTheta
Direction of the muon in the local coordinate frame.
Identifier stripId
Identifier of the strip.
uint8_t clusterAuthor
Author of the cluster.

◆ digitizeWire()

bool MuonR4::sTgcFastDigiTool::digitizeWire ( const EventContext & ctx,
const TimedHit & timedHit,
const Muon::DigitEffiData * effiData,
CLHEP::HepRandomEngine * rndEngine,
sTgcDigitCollection & outCollection ) const
private

Check efficiencies

Check efficiencies

Definition at line 220 of file sTgcFastDigiTool.cxx.

224 {
225
226 if (!m_digitizeWire) {
227 return false;
228 }
229
230 ++m_allHits[channelType::Wire];
231
232 const Identifier hitId{timedHit->identify()};
234 if (efficiencyMap && efficiencyMap->getEfficiency(hitId) < CLHEP::RandFlat::shoot(rndEngine,0.,1.)){
235 ATH_MSG_VERBOSE("Simulated wire hit "<<xAOD::toEigen(timedHit->localPosition())
236 << m_idHelperSvc->toString(hitId) <<" is rejected because of efficency modelling");
237 return false;
238 }
239 const sTgcIdHelper& idHelper{m_idHelperSvc->stgcIdHelper()};
240 const MuonGMR4::sTgcReadoutElement* readOutEle = m_detMgr->getsTgcReadoutElement(hitId);
241 const IdentifierHash hitHash = readOutEle->measurementHash(hitId);
242 const int gasGap = idHelper.gasGap(hitId);
243
244 const Amg::Vector2D wirePos = readOutEle->stripLayer(hitHash).to2D(xAOD::toEigen(timedHit->localPosition()), true);
245 // do not digitise wires that are never read out in reality
246 bool isInnerQ1 = readOutEle->isEtaZero(hitHash, xAOD::toEigen(timedHit->localPosition()).block<2,1>(0,0));
247 if(isInnerQ1) {
248 return false;
249 }
251 if (efficiencyMap && efficiencyMap->getEfficiency(hitId, isInnerQ1) < CLHEP::RandFlat::shoot(rndEngine,0.,1.)){
252 ATH_MSG_VERBOSE("Simulated wire hit "<<xAOD::toEigen(timedHit->localPosition())
253 << m_idHelperSvc->toString(hitId) <<" is rejected because of efficency modelling");
254 return false;
255 }
256
257 const MuonGMR4::WireGroupDesign& design{readOutEle->wireDesign(hitHash)};
258
259 const int wireGrpNum = design.stripNumber(wirePos);
260 if (wireGrpNum < 0) {
261 ATH_MSG_VERBOSE("The wire "<<Amg::toString(wirePos)<<" in "<<m_idHelperSvc->toStringGasGap(hitId)
262 <<" is outside of the acceptance of "<<std::endl<<design);
263 return false;
264 }
265 const double uncert = design.stripPitch() * design.numWiresInGroup(wireGrpNum);
266
267 const double smearedX = CLHEP::RandGaussZiggurat::shoot(rndEngine, wirePos.x(), uncert);
268
269 const Amg::Vector2D digitPos{smearedX * Amg::Vector2D::UnitX()};
270
271 const int digitWire = design.stripNumber(digitPos);
272 if (digitWire < 0) {
273 ATH_MSG_VERBOSE("Strip hit "<<Amg::toString(digitPos)<<" "<<m_idHelperSvc->toStringGasGap(hitId)
274 <<" is out of range "<<std::endl<<design);
275 return false;
276 }
277 bool isValid{false};
278 const Identifier digitId = idHelper.channelID(hitId, readOutEle->multilayer(),
279 gasGap, channelType::Wire, digitWire, isValid);
280
281 if (!isValid) {
282 ATH_MSG_WARNING("Failed to deduce a valid identifier from "
283 <<m_idHelperSvc->toStringGasGap(hitId)<<" digit: "<<digitWire);
284 return false;
285 }
286 outCollection.push_back(std::make_unique<sTgcDigit>(digitId,
287 associateBCIdTag(ctx, timedHit),
288 hitTime(timedHit), 666, false, false));
289
290
291 ++m_acceptedHits[channelType::Wire];
292 return true;
293 }
int stripNumber(const Amg::Vector2D &pos) const override
Calculates the number of the strip whose center is closest to the given point.
unsigned int numWiresInGroup(unsigned int groupNum) const
Returns the number of wires in a given group.
const WireGroupDesign & wireDesign(const IdentifierHash &measHash) const
Retrieves the readoutElement Layer given the Identifier/Hash.
Gaudi::Property< bool > m_digitizeWire

◆ fetchCollection()

template<class DigitColl>
DigitColl * MuonR4::MuonDigitizationTool::fetchCollection ( const Identifier & hitId,
OutDigitCache_t< DigitColl > & digitCache ) const
protectedinherited

Helper function that provides fetches the proper DigitCollection from the DigitCache for a given hit identifier If the Collection is fetched for the first time, it's inserted into the cache first.

◆ fillTimedHits()

StatusCode MuonR4::MuonDigitizationTool::fillTimedHits ( PileUpHits && hitColl,
TimedHits & timedHits ) const
privateinherited

Translates the PileUpHits into the timed hits format.

Definition at line 44 of file MuonDigitizationTool.cxx.

44 {
45 for (const auto& [timeIndex, simHitColl] : hitColl) {
46 timedHits.reserve(timedHits.capacity() + simHitColl->size());
47 for (const xAOD::MuonSimHit* simHit : *simHitColl) {
48 timedHits.emplace_back(timeIndex.time(), timeIndex.index(), simHit, timeIndex.type());
49 }
50 }
51 std::stable_sort(timedHits.begin(), timedHits.end(),
52 [](const TimedHit& a, const TimedHit& b){
53 if (a->identify() != b->identify()){
54 return a->identify() < b->identify();
55 }
56 if (a.eventId() != b.eventId()) {
57 return a.eventId() < b.eventId();
58 }
59 return a.eventTime() < b.eventTime();
60 });
61 return StatusCode::SUCCESS;
62 }
static Double_t a
if(febId1==febId2)
void stable_sort(DataModel_detail::iterator< DVL > beg, DataModel_detail::iterator< DVL > end)
Specialization of stable_sort for DataVector/List.

◆ filterPassed()

virtual bool PileUpToolBase::filterPassed ( ) const
inlineoverridevirtualinherited

dummy implementation of passing filter

Definition at line 49 of file PileUpToolBase.h.

49{ return m_filterPassed; }

◆ finalize()

StatusCode MuonR4::sTgcFastDigiTool::finalize ( )
finaloverride

Definition at line 30 of file sTgcFastDigiTool.cxx.

30 {
31 ATH_MSG_INFO("Tried to convert "<<m_allHits[channelType::Strip]<<"/"
32 <<m_allHits[channelType::Wire]<<"/"
33 <<m_allHits[channelType::Pad]<<" strip/wire/pad hits. In, "
34 <<percentage(m_acceptedHits[channelType::Strip], m_allHits[channelType::Strip]) <<"/"
35 <<percentage(m_acceptedHits[channelType::Wire], m_allHits[channelType::Wire])<<"/"
36 <<percentage(m_acceptedHits[channelType::Pad], m_allHits[channelType::Pad])
37 <<"% of the cases, the conversion was successful");
38 return StatusCode::SUCCESS;
39 }
#define ATH_MSG_INFO(x)

◆ getGeoCtx()

const ActsTrk::GeometryContext & MuonR4::MuonDigitizationTool::getGeoCtx ( const EventContext & ctx) const
protectedinherited

Returns the reference to the ActsTrk::GeometryContext needed to fetch global positions from the Readout geometry.

Definition at line 135 of file MuonDigitizationTool.cxx.

135 {
136 const ActsTrk::GeometryContext* gctx{};
137 if (!SG::get(gctx, m_geoCtxKey, ctx).isSuccess()) {
138 THROW_EXCEPTION("Failed to retrieve the geometry context "<<m_geoCtxKey.fullKey());
139 }
140 return *gctx;
141 }
SG::ReadHandleKey< ActsTrk::GeometryContext > m_geoCtxKey
#define THROW_EXCEPTION(MESSAGE)
Definition throwExcept.h:10

◆ getRandomEngine()

CLHEP::HepRandomEngine * MuonR4::MuonDigitizationTool::getRandomEngine ( const EventContext & ctx) const
protectedinherited

Definition at line 129 of file MuonDigitizationTool.cxx.

129 {
130 ATHRNG::RNGWrapper* rngWrapper = m_rndmSvc->getEngine(this, m_streamName);
131 std::string rngName = m_streamName;
132 rngWrapper->setSeed(rngName, ctx);
133 return rngWrapper->getEngine(ctx);
134 }
void setSeed(const std::string &algName, const EventContext &ctx)
Set the random seed using a string (e.g.
Definition RNGWrapper.h:169
CLHEP::HepRandomEngine * getEngine(const EventContext &ctx) const
Retrieve the random engine corresponding to the provided EventContext.
Definition RNGWrapper.h:134
ServiceHandle< IAthRNGSvc > m_rndmSvc
Gaudi::Property< std::string > m_streamName

◆ hitTime()

double MuonR4::MuonDigitizationTool::hitTime ( const TimedHit & hit)
staticprotectedinherited

Returns the global time of the hit which is the sum of eventTime & individual hit time.

Definition at line 63 of file MuonDigitizationTool.cxx.

63 {
64 return hit.eventTime() + hit->globalTime();
65 }

◆ initialize()

StatusCode MuonR4::sTgcFastDigiTool::initialize ( )
finaloverridevirtual

Reimplemented from MuonR4::MuonDigitizationTool.

Definition at line 23 of file sTgcFastDigiTool.cxx.

23 {
25 ATH_CHECK(m_writeKey.initialize());
26 ATH_CHECK(m_effiDataKey.initialize(!m_effiDataKey.empty()));
27 ATH_CHECK(m_uncertCalibKey.initialize());
28 return StatusCode::SUCCESS;
29 }

◆ mergeEvent()

StatusCode MuonR4::MuonDigitizationTool::mergeEvent ( const EventContext & ctx)
finaloverrideinherited

Definition at line 92 of file MuonDigitizationTool.cxx.

92 {
93 ATH_MSG_DEBUG("mergeEvent()");
94
95 SG::WriteHandle<xAOD::MuonSimHitContainer> sdoContainer{};
96 if (!m_sdoKey.empty()) {
97 sdoContainer = SG::WriteHandle<xAOD::MuonSimHitContainer>{m_sdoKey, ctx};
98 ATH_CHECK(sdoContainer.record(std::make_unique<xAOD::MuonSimHitContainer>(),
99 std::make_unique<xAOD::MuonSimHitAuxContainer>()));
100 }
101 ATH_CHECK(digitize(ctx, m_timedHits, !m_sdoKey.empty() ? sdoContainer.ptr() : nullptr));
102 m_timedHits.clear();
103 m_simHits.clear();
104 return StatusCode::SUCCESS;
105 }
std::vector< SimHitLocalCopy > m_simHits
virtual StatusCode digitize(const EventContext &ctx, const TimedHits &hitsToDigit, xAOD::MuonSimHitContainer *sdoContainer) const =0
Digitize the time ordered hits and write them to the digit format specific for the detector technolog...
SG::WriteHandleKey< xAOD::MuonSimHitContainer > m_sdoKey
StatusCode record(std::unique_ptr< T > data)
Record a const object to the store.
pointer_type ptr()
Dereference the pointer.

◆ passDeadTime()

bool MuonR4::MuonDigitizationTool::passDeadTime ( const Identifier & channelId,
const double hitTime,
const double deadTimeWindow,
DeadTimeMap & deadTimeMap )
staticprotectedinherited

Returns whether the new digit is within the dead time window.

channelId: Identifier of the fired channel

hitTime: Current hit time

deadTimeWindow: Dead time to pass before the next hit may be accepted

deadTimeMap: Map storing the last hit times from each Identifier

Channel not seen before

Update dead time map & accept hit

Definition at line 164 of file MuonDigitizationTool.cxx.

167 {
168 auto insertItr = deadTimeMap.insert(std::make_pair(channelId,hitTime));
170 if (insertItr.second) {
171 return true;
172 }
173 if (hitTime - insertItr.first->second < deadTimeWindow) {
174 return false;
175 }
177 insertItr.first->second = hitTime;
178 return true;
179 }

◆ PileUpToolBase()

PileUpToolBase::PileUpToolBase ( const std::string & type,
const std::string & name,
const IInterface * parent )
inherited

Definition at line 22 of file PileUpToolBase.cxx.

10 : base_class(type, name, parent)
11{
12}

◆ prepareEvent()

StatusCode MuonR4::MuonDigitizationTool::prepareEvent ( const EventContext & ctx,
const unsigned int nInputEvents )
finaloverrideinherited

When being run from PileUpToolsAlgs, this method is called at the start of the subevts loop.

Not able to access SubEvents

Definition at line 35 of file MuonDigitizationTool.cxx.

36 {
37
38 ATH_MSG_DEBUG("prepareEvent() called for " << nInputEvents << " input events");
39 m_timedHits.clear();
40 m_simHits.clear();
41 return StatusCode::SUCCESS;
42 }

◆ processAllSubEvents() [1/2]

StatusCode MuonR4::MuonDigitizationTool::processAllSubEvents ( const EventContext & ctx) const
inherited

Reentrant version of the digitization tool.

In case of single hits container just load the collection using read handles

Definition at line 71 of file MuonDigitizationTool.cxx.

71 {
72 PileUpHits hitCollList{};
73 TimedHits timedHits{};
76 const xAOD::MuonSimHitContainer* hitCollection{nullptr};
77 ATH_CHECK(SG::get(hitCollection, m_simHitKey, ctx));
78 hitCollList.emplace_back(PileUpTimeEventIndex(0), hitCollection);
79 } else {
80 ATH_CHECK(m_mergeSvc->retrieveSubEvtsData(m_inputObjectName, hitCollList));
81 }
82 ATH_CHECK(fillTimedHits(std::move(hitCollList), timedHits));
83 SG::WriteHandle<xAOD::MuonSimHitContainer> sdoContainer{};
84 if (!m_sdoKey.empty()) {
85 sdoContainer = SG::WriteHandle<xAOD::MuonSimHitContainer>{m_sdoKey, ctx};
86 ATH_CHECK(sdoContainer.record(std::make_unique<xAOD::MuonSimHitContainer>(),
87 std::make_unique<xAOD::MuonSimHitAuxContainer>()));
88 }
89 ATH_CHECK(digitize(ctx, timedHits, !m_sdoKey.empty() ? sdoContainer.ptr() : nullptr));
90 return StatusCode::SUCCESS;
91 }
ServiceHandle< PileUpMergeSvc > m_mergeSvc
std::vector< TimedHitPtr< xAOD::MuonSimHit > > TimedHits
Gaudi::Property< bool > m_onlyUseContainerName
SG::ReadHandleKey< xAOD::MuonSimHitContainer > m_simHitKey
StatusCode fillTimedHits(PileUpHits &&hitColl, TimedHits &timedHits) const
Translates the PileUpHits into the timed hits format.
PileUpMergeSvc::TimedList< xAOD::MuonSimHitContainer >::type PileUpHits
MuonSimHitContainer_v1 MuonSimHitContainer
Define the version of the pixel cluster container.

◆ processAllSubEvents() [2/2]

StatusCode MuonR4::MuonDigitizationTool::processAllSubEvents ( const EventContext & ctx)
finaloverridevirtualinherited

alternative interface which uses the PileUpMergeSvc to obtain all the required SubEvents.

Reimplemented from PileUpToolBase.

Definition at line 66 of file MuonDigitizationTool.cxx.

66 {
67 const MuonDigitizationTool* digiTool = this;
68 return digiTool->processAllSubEvents(ctx);
69 }

◆ processBunchXing()

StatusCode MuonR4::MuonDigitizationTool::processBunchXing ( int bunchXing,
SubEventIterator bSubEvents,
SubEventIterator eSubEvents )
finaloverridevirtualinherited

Reimplemented from PileUpToolBase.

Definition at line 108 of file MuonDigitizationTool.cxx.

110 {
111 ATH_MSG_DEBUG("processBunchXing()" << bunchXing);
112 PileUpHits hitList{}, hitListPermanent{};
113 ATH_CHECK(m_mergeSvc->retrieveSubSetEvtData(m_inputObjectName, hitList, bunchXing, bSubEvents, eSubEvents));
114 ATH_MSG_VERBOSE(hitList.size() << " hits in xAODMuonSimHitContainer " << m_inputObjectName << " found");
115 for (auto& [hitPtr, hitContainer] : hitList) {
116 auto copyContainer = std::make_unique<xAOD::MuonSimHitContainer>();
117 auto copyAuxContainer = std::make_unique<xAOD::MuonSimHitAuxContainer>();
118 copyContainer->setStore(copyAuxContainer.get());
119 for (const xAOD::MuonSimHit* copyMe : *hitContainer) {
120 (*copyContainer->push_back(std::make_unique<xAOD::MuonSimHit>())) = (*copyMe);
121 }
122 hitListPermanent.emplace_back(hitPtr, copyContainer.get());
123 m_simHits.emplace_back(std::move(copyContainer), std::move(copyAuxContainer));
124 }
125 ATH_CHECK(fillTimedHits(std::move(hitListPermanent), m_timedHits));
126 return StatusCode::SUCCESS;
127 }

◆ resetFilter()

virtual void PileUpToolBase::resetFilter ( )
inlineoverridevirtualinherited

dummy implementation of filter reset

Reimplemented in MergeTruthJetsTool.

Definition at line 51 of file PileUpToolBase.h.

51{ m_filterPassed=true; }

◆ toProcess()

virtual bool PileUpToolBase::toProcess ( int bunchXing) const
inlineoverridevirtualinherited

the method this base class helps implementing

Reimplemented in MergeHijingParsTool, and MergeTrackRecordCollTool.

Definition at line 32 of file PileUpToolBase.h.

32 {
33 //closed interval [m_firstXing,m_lastXing]
34 return !((m_firstXing > bunchXing) || (bunchXing > m_lastXing));
35 }
Gaudi::Property< int > m_firstXing
Gaudi::Property< int > m_lastXing

◆ writeDigitContainer()

template<class DigitCont, class DigitColl>
StatusCode MuonR4::MuonDigitizationTool::writeDigitContainer ( const EventContext & ctx,
const SG::WriteHandleKey< DigitCont > & key,
OutDigitCache_t< DigitColl > && digitCache,
unsigned int hashMax ) const
protectedinherited

Helper function to move the collected digits into the final DigitContainer.

The function needs the maximal size of the container in advance which is provided by calling the module_hash_max() function of the corresponding MuonIdHelper.

Member Data Documentation

◆ ATLAS_THREAD_SAFE [1/2]

std::array<std::atomic<unsigned>, 3> m_allHits MuonR4::sTgcFastDigiTool::ATLAS_THREAD_SAFE {}
mutableprivate

Definition at line 64 of file sTgcFastDigiTool.h.

64{};

◆ ATLAS_THREAD_SAFE [2/2]

std::array<std::atomic<unsigned>, 3> m_acceptedHits MuonR4::sTgcFastDigiTool::ATLAS_THREAD_SAFE {}
mutableprivate

Definition at line 65 of file sTgcFastDigiTool.h.

65{};

◆ m_deadTime

Gaudi::Property<double> MuonR4::sTgcFastDigiTool::m_deadTime {this, "deadTime", 300. * Gaudi::Units::ns}
private

Definition at line 71 of file sTgcFastDigiTool.h.

71{this, "deadTime", 300. * Gaudi::Units::ns};

◆ m_detMgr

const MuonGMR4::MuonDetectorManager* MuonR4::MuonDigitizationTool::m_detMgr {nullptr}
protectedinherited

Definition at line 111 of file MuonDigitizationTool.h.

111{nullptr};

◆ m_digitizeMuonOnly

Gaudi::Property<bool> MuonR4::sTgcFastDigiTool::m_digitizeMuonOnly
private
Initial value:
{this, "ProcessTrueMuonsOnly", false,
"If set to true hit with pdgId != 13 are skipped"}

Definition at line 73 of file sTgcFastDigiTool.h.

73 {this, "ProcessTrueMuonsOnly", false,
74 "If set to true hit with pdgId != 13 are skipped"};

◆ m_digitizePads

Gaudi::Property<bool> MuonR4::sTgcFastDigiTool::m_digitizePads {this, "doPads", true, "Digitize pad hits"}
private

Definition at line 69 of file sTgcFastDigiTool.h.

69{this, "doPads", true, "Digitize pad hits"};

◆ m_digitizeStrip

Gaudi::Property<bool> MuonR4::sTgcFastDigiTool::m_digitizeStrip {this, "doStrips", true, "Digitize strip hits"}
private

Definition at line 67 of file sTgcFastDigiTool.h.

67{this, "doStrips", true, "Digitize strip hits"};

◆ m_digitizeWire

Gaudi::Property<bool> MuonR4::sTgcFastDigiTool::m_digitizeWire {this, "doWires", true, "Digitize wire hits"}
private

Definition at line 68 of file sTgcFastDigiTool.h.

68{this, "doWires", true, "Digitize wire hits"};

◆ m_effiDataKey

SG::ReadCondHandleKey<Muon::DigitEffiData> MuonR4::sTgcFastDigiTool::m_effiDataKey
private
Initial value:
{this, "EffiDataKey", "sTgcDigitEff",
"Efficiency constants of the individual sTgc gasGaps"}

Definition at line 57 of file sTgcFastDigiTool.h.

57 {this, "EffiDataKey", "sTgcDigitEff",
58 "Efficiency constants of the individual sTgc gasGaps"};

◆ m_energyDepositThreshold

Gaudi::Property<double> MuonR4::sTgcFastDigiTool::m_energyDepositThreshold {this,"energyDepositThreshold",300.0*CLHEP::eV}
private

Definition at line 75 of file sTgcFastDigiTool.h.

75{this,"energyDepositThreshold",300.0*CLHEP::eV};

◆ m_filterPassed

bool PileUpToolBase::m_filterPassed {true}
protectedinherited

Definition at line 60 of file PileUpToolBase.h.

60{true};

◆ m_firstXing

Gaudi::Property<int> PileUpToolBase::m_firstXing
protectedinherited
Initial value:
{this, "FirstXing", -999,
"First bunch-crossing in which det is live"}

Definition at line 54 of file PileUpToolBase.h.

54 {this, "FirstXing", -999,
55 "First bunch-crossing in which det is live"};

◆ m_geoCtxKey

SG::ReadHandleKey<ActsTrk::GeometryContext> MuonR4::MuonDigitizationTool::m_geoCtxKey
privateinherited
Initial value:
{this, "AlignmentKey", "ActsAlignment",
"Geometry context"}

Definition at line 123 of file MuonDigitizationTool.h.

123 {this, "AlignmentKey", "ActsAlignment",
124 "Geometry context"};

◆ m_idHelperSvc

ServiceHandle<Muon::IMuonIdHelperSvc> MuonR4::MuonDigitizationTool::m_idHelperSvc
protectedinherited
Initial value:
{this, "MuonIdHelperSvc",
"Muon::MuonIdHelperSvc/MuonIdHelperSvc"}

Definition at line 113 of file MuonDigitizationTool.h.

113 {this, "MuonIdHelperSvc",
114 "Muon::MuonIdHelperSvc/MuonIdHelperSvc"};

◆ m_includePileUpTruth

Gaudi::Property<bool> MuonR4::MuonDigitizationTool::m_includePileUpTruth {this, "IncludePileUpTruth", true, "Include pile-up truth info"}
privateinherited

Definition at line 137 of file MuonDigitizationTool.h.

137{this, "IncludePileUpTruth", true, "Include pile-up truth info"};

◆ m_inputObjectName

std::string MuonR4::MuonDigitizationTool::m_inputObjectName {""}
privateinherited

Definition at line 139 of file MuonDigitizationTool.h.

139{""};

◆ m_lastXing

Gaudi::Property<int> PileUpToolBase::m_lastXing
protectedinherited
Initial value:
{this, "LastXing", 999,
"Last bunch-crossing in which det is live"}

Definition at line 56 of file PileUpToolBase.h.

56 {this, "LastXing", 999,
57 "Last bunch-crossing in which det is live"};

◆ m_limitElectronKineticEnergy

Gaudi::Property<double> MuonR4::sTgcFastDigiTool::m_limitElectronKineticEnergy {this,"limitElectronKineticEnergy",5.0*CLHEP::MeV}
private

Definition at line 76 of file sTgcFastDigiTool.h.

76{this,"limitElectronKineticEnergy",5.0*CLHEP::MeV};

◆ m_mergeSvc

ServiceHandle<PileUpMergeSvc> MuonR4::MuonDigitizationTool::m_mergeSvc {this, "PileUpMergeSvc", "PileUpMergeSvc", ""}
privateinherited

Definition at line 126 of file MuonDigitizationTool.h.

126{this, "PileUpMergeSvc", "PileUpMergeSvc", ""};

◆ m_onlyUseContainerName

Gaudi::Property<bool> MuonR4::MuonDigitizationTool::m_onlyUseContainerName
privateinherited
Initial value:
{this, "OnlyUseContainerName", false,
"Don't use the ReadHandleKey directly. Just extract the container name from it."}

Definition at line 134 of file MuonDigitizationTool.h.

134 {this, "OnlyUseContainerName", false,
135 "Don't use the ReadHandleKey directly. Just extract the container name from it."};

◆ m_rndmSvc

ServiceHandle<IAthRNGSvc> MuonR4::MuonDigitizationTool::m_rndmSvc {this, "RndmSvc", "AthRNGSvc", ""}
privateinherited

Definition at line 130 of file MuonDigitizationTool.h.

130{this, "RndmSvc", "AthRNGSvc", ""}; // Random number service

◆ m_sdoKey

SG::WriteHandleKey<xAOD::MuonSimHitContainer> MuonR4::MuonDigitizationTool::m_sdoKey {this, "OutputSDOName", ""}
privateinherited

Definition at line 132 of file MuonDigitizationTool.h.

132{this, "OutputSDOName", ""};

◆ m_simHitKey

SG::ReadHandleKey<xAOD::MuonSimHitContainer> MuonR4::MuonDigitizationTool::m_simHitKey {this, "SimHitKey", ""}
privateinherited

Definition at line 121 of file MuonDigitizationTool.h.

121{this, "SimHitKey", ""};

◆ m_simHits

std::vector<SimHitLocalCopy> MuonR4::MuonDigitizationTool::m_simHits {}
privateinherited

Definition at line 145 of file MuonDigitizationTool.h.

145{};

◆ m_streamName

Gaudi::Property<std::string> MuonR4::MuonDigitizationTool::m_streamName {this, "StreamName", ""}
privateinherited

Definition at line 128 of file MuonDigitizationTool.h.

128{this, "StreamName", ""};

◆ m_timedHits

TimedHits MuonR4::MuonDigitizationTool::m_timedHits {}
privateinherited

Definition at line 141 of file MuonDigitizationTool.h.

141{};

◆ m_uncertCalibKey

SG::ReadCondHandleKey<NswErrorCalibData> MuonR4::sTgcFastDigiTool::m_uncertCalibKey
private
Initial value:
{this, "ErrorCalibKey", "NswUncertData",
"Key of the parametrized NSW uncertainties"}

Definition at line 61 of file sTgcFastDigiTool.h.

61 {this, "ErrorCalibKey", "NswUncertData",
62 "Key of the parametrized NSW uncertainties"};

◆ m_vetoPileUpTruthLinks

Gaudi::Property<int> PileUpToolBase::m_vetoPileUpTruthLinks
protectedinherited
Initial value:
{this, "VetoPileUpTruthLinks", true,
"Ignore links to suppressed pile-up truth"}

Definition at line 58 of file PileUpToolBase.h.

58 {this, "VetoPileUpTruthLinks", true,
59 "Ignore links to suppressed pile-up truth"};

◆ m_writeKey

SG::WriteHandleKey<sTgcDigitContainer> MuonR4::sTgcFastDigiTool::m_writeKey {this, "OutputObjectName", "sTGC_DIGITS"}
private

Definition at line 55 of file sTgcFastDigiTool.h.

55{this, "OutputObjectName", "sTGC_DIGITS"};

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