ATLAS Offline Software
Loading...
Searching...
No Matches
Muon::MuonTruthSegmentCreationAlg Class Reference

#include <MuonTruthSegmentCreationAlg.h>

Inheritance diagram for Muon::MuonTruthSegmentCreationAlg:
Collaboration diagram for Muon::MuonTruthSegmentCreationAlg:

Public Types

using ChamberIdMap = std::map<Muon::MuonStationIndex::ChIndex, std::vector<Identifier>>
 This map contains all the hits corresponding to truth muons classified by chamber layer that recorded them.

Public Member Functions

virtual StatusCode initialize () override
virtual StatusCode execute (const EventContext &ctx) const override
virtual StatusCode sysInitialize () override
 Override sysInitialize.
virtual bool isClonable () const override
 Specify if the algorithm is clonable.
virtual unsigned int cardinality () const override
 Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.
virtual StatusCode sysExecute (const EventContext &ctx) override
 Execute an algorithm.
virtual const DataObjIDColl & extraOutputDeps () const override
 Return the list of extra output dependencies.
virtual bool filterPassed (const EventContext &ctx) const
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.
const ServiceHandle< StoreGateSvc > & detStore () const
 The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
virtual StatusCode sysStart () override
 Handle START transition.
virtual std::vector< Gaudi::DataHandle * > inputHandles () const override
 Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * > outputHandles () const override
 Return this algorithm's output handles.
Gaudi::Details::PropertyBase & declareProperty (Gaudi::Property< T, V, H > &t)
void updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream & msg () const
bool msgLvl (const MSG::Level lvl) const

Protected Member Functions

void renounceArray (SG::VarHandleKeyArray &handlesArray)
 remove all handles from I/O resolution
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > renounce (T &h)
void extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
 Add StoreName to extra input/output deps as needed.

Private Types

typedef ServiceHandle< StoreGateSvcStoreGateSvc_t

Private Member Functions

StatusCode fillChamberIdMap (const EventContext &ctx, const xAOD::TruthParticle &truthParticle, ChamberIdMap &ids) const
 This function uses the 6 vectors, contained in.
StatusCode createSegments (const EventContext &ctx, const ElementLink< xAOD::TruthParticleContainer > &truthLink, const ChamberIdMap &ids, xAOD::MuonSegmentContainer &segmentContainer) const
 This function performs, for each truth muon, the actual segment creation and stores segments into a new container.
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes

SG::ReadHandleKey< xAOD::TruthParticleContainerm_muonTruth {this, "muonTruth", "MuonTruthParticles"}
 Key for the truth muon container and muon origin decoration.
SG::ReadDecorHandleKey< xAOD::TruthParticleContainerm_truthOriginKey {this, "truthOriginKey", m_muonTruth, "truthOrigin"}
SG::ReadDecorHandleKeyArray< xAOD::TruthParticleContainer, std::vector< unsigned long long > > m_truthHitsKeyArray {this, "truthHitsKeyArray", { }}
 Keys of the truth muon decorations that we need to read to (re-)fill the chamberIdMap.
SG::WriteHandleKey< xAOD::MuonSegmentContainerm_muonTruthSegmentContainerName {this, "MuonTruthSegmentName", "MuonTruthSegments"}
 Key for segment container that will be populated with segments.
SG::ReadHandleKeyArray< MuonSimDataCollectionm_SDO_TruthNames
 Keys for all çontainers of muon hit simulation data, classified by detector technology.
SG::ReadHandleKey< CscSimDataCollectionm_CSC_SDO_TruthNames {this, "CSCSDOs", "CSC_SDO"}
ServiceHandle< Muon::IMuonIdHelperSvcm_idHelperSvc {this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}
 Handle for the muonIdHelper service.
SG::ReadCondHandleKey< MuonGM::MuonDetectorManagerm_detMgrKey
 MuonDetectorManager from the conditions store.
DataObjIDColl m_extendedExtraObjects
 Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default)
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default)
std::vector< SG::VarHandleKeyArray * > m_vhka
bool m_varHandleArraysDeclared

Detailed Description

Definition at line 38 of file MuonTruthSegmentCreationAlg.h.

Member Typedef Documentation

◆ ChamberIdMap

This map contains all the hits corresponding to truth muons classified by chamber layer that recorded them.

ChIndex is the chamber layer index (e.g. BMS, BOL,...), while Identifier is the hit ID, which represents the channel where the hit has been recorded.

Definition at line 52 of file MuonTruthSegmentCreationAlg.h.

◆ StoreGateSvc_t

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

Definition at line 388 of file AthCommonDataStore.h.

Member Function Documentation

◆ cardinality()

unsigned int AthCommonReentrantAlgorithm< Gaudi::Algorithm >::cardinality ( ) const
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.

Definition at line 75 of file AthCommonReentrantAlgorithm.cxx.

64{
65 return 0;
66}

◆ createSegments()

StatusCode Muon::MuonTruthSegmentCreationAlg::createSegments ( const EventContext & ctx,
const ElementLink< xAOD::TruthParticleContainer > & truthLink,
const ChamberIdMap & ids,
xAOD::MuonSegmentContainer & segmentContainer ) const
private

This function performs, for each truth muon, the actual segment creation and stores segments into a new container.

Parameters
segmentContainer.Firstly, we retrieve simulation data from m_SDO_TruthNames and copy it into a new vector, keeping only the data corresponding to the correct detector technologies. Then, we loop over chamber layers and over the corresponding hit IDs, both provided by
ids,andwe store: *- gas-gap for hits with a phi measurement -> phiLayers *- gas-gap for hits recorded by MM or CSC -> precLayers *- 10*multi-layer + tube_layer for hits recorded by MDT -> precLayers *- gas-gap for hits recorded by trigger chambers without phi measurement -> etaLayers Then, in the same loop over hit IDs, we retrieve simulation data for that hit, specifically the global position, needed to retrieve for each chamber layer the pair of hits with the largest distance. If a pair of hits is found and at least three precision measurements are available, we proceed with the truth segment creation. We create a segment obj and we fill it with number and type (see above) of hits available, the
truthLinklink to the corresponding truth muon, and the position and direction of the truth segment, calculated with hit pair global positions. Finally, we append the sector, chamber layer, station eta, and technology of the chamber where the segment was recorded.
ctxtells us in which event we are
truthLinklink to the corresponding truth muons
idschamberIdMap storing truth hit IDs classified by chamber layer
segmentContainersegment container

Definition at line 119 of file MuonTruthSegmentCreationAlg.cxx.

122 {
123
124 const MuonGM::MuonDetectorManager* detMgr{nullptr};
125 ATH_CHECK(SG::get(detMgr, m_detMgrKey, ctx));
126
127 constexpr unsigned techMax = toInt(TechnologyIndex::TechnologyIndexMax);
128 std::array<const MuonSimDataCollection*, techMax> sdoCollections{};
129 bool useSDO = !m_CSC_SDO_TruthNames.empty();
130 for (const SG::ReadHandleKey<MuonSimDataCollection>& k : m_SDO_TruthNames) {
131 const MuonSimDataCollection* coll{nullptr};
132 ATH_CHECK(SG::get(coll, k, ctx));
133 if (coll->empty()) {
134 continue;
135 }
136 Identifier id = coll->begin()->first;
137 sdoCollections[toInt(m_idHelperSvc->technologyIndex(id))] = coll;
138 useSDO = true;
139 }
140 // loop over chamber layers
141 for (const auto& [chIdx, assocIds] : ids) {
142 // skip empty layers
143 Amg::Vector3D firstPos{Amg::Vector3D::Zero()}, secondPos{Amg::Vector3D::Zero()};
144 bool firstPosSet{false}, secondPosSet{false};
145 Identifier chId{};
146 int index = -1;
147 uint8_t nprecLayers{0}, nphiLayers{0}, ntrigEtaLayers{0};
148 std::unordered_set<int> phiLayers{}, etaLayers{}, precLayers{};
149 ATH_MSG_DEBUG(" new chamber layer " << Muon::MuonStationIndex::chName(chIdx) << " hits " << assocIds.size());
150 // loop over hits
151 for (const auto& id : assocIds) {
152 ATH_MSG_VERBOSE(" hit " << m_idHelperSvc->toString(id));
153 bool measPhi = m_idHelperSvc->measuresPhi(id);
154 bool isCsc = m_idHelperSvc->isCsc(id);
155 bool isMM = m_idHelperSvc->isMM(id);
156 bool isTrig = m_idHelperSvc->isTrigger(id);
157 bool isEndcap = m_idHelperSvc->isEndcap(id);
158 if (measPhi) {
159 phiLayers.insert(m_idHelperSvc->gasGap(id));
160 } else {
161 if (!isTrig) {
162 if (!chId.is_valid()) chId = id; // use first precision hit in list
163 if (isCsc || isMM) {
164 precLayers.insert(m_idHelperSvc->gasGap(id));
165 } else {
166 int iid = 10 * m_idHelperSvc->mdtIdHelper().multilayer(id) + m_idHelperSvc->mdtIdHelper().tubeLayer(id);
167 precLayers.insert(iid);
168 // ATH_MSG_VERBOSE("iid " << iid << " precLayers size " << precLayers.size() );
169 }
170 } else {
171 etaLayers.insert(m_idHelperSvc->gasGap(id));
172 }
173 }
174 // use SDO to look-up truth position of the hit
175 if (!useSDO) {
176 continue;
177 }
178 Amg::Vector3D gpos{Amg::Vector3D::Zero()};
179 if (!isCsc) {
180 bool ok = false;
181 TechnologyIndex techIdx = m_idHelperSvc->technologyIndex(id);
182 if (sdoCollections[toInt(techIdx)]) {
183 auto pos = sdoCollections[toInt(techIdx)]->find(id);
184 if (pos != sdoCollections[toInt(techIdx)]->end()) {
185 gpos = pos->second.globalPosition();
186 if (gpos.perp() > 0.1) ok = true; // sanity check
187 }
188 }
189 // look up successful, calculate
190 if (!ok) continue;
191
192 // small comparison function
193 auto isSmaller = [isEndcap](const Amg::Vector3D& p1, const Amg::Vector3D& p2) {
194 if (isEndcap)
195 return std::abs(p1.z()) < std::abs(p2.z());
196 else
197 return p1.perp() < p2.perp();
198 };
199 if (!firstPosSet) {
200 firstPos = gpos;
201 firstPosSet = true;
202 } else if (!secondPosSet) {
203 secondPos = gpos;
204 secondPosSet = true;
205 if (isSmaller(secondPos, firstPos)) std::swap(firstPos, secondPos);
206 } else {
207 // update position if we can increase the distance between the two positions
208 if (isSmaller(gpos, firstPos))
209 firstPos = gpos;
210 else if (isSmaller(secondPos, gpos))
211 secondPos = gpos;
212 }
213 } else {
214 SG::ReadHandle cscCollection(m_CSC_SDO_TruthNames, ctx);
215 ATH_CHECK(cscCollection.isPresent());
216 auto pos = cscCollection->find(id);
217 if (pos == cscCollection->end()) {
218 continue;
219 }
220 const MuonGM::CscReadoutElement* descriptor = detMgr->getCscReadoutElement(id);
221 ATH_MSG_DEBUG("found csc sdo with " << pos->second.getdeposits().size() << " deposits");
222 Amg::Vector3D locpos(0, pos->second.getdeposits()[0].second.ypos(), pos->second.getdeposits()[0].second.zpos());
223 gpos = descriptor->localToGlobalCoords(locpos, m_idHelperSvc->cscIdHelper().elementID(id));
224 ATH_MSG_DEBUG("got CSC global position " << gpos);
225 if (!firstPosSet) {
226 firstPos = gpos;
227 firstPosSet = true;
228 } else if (!secondPosSet) {
229 secondPos = gpos;
230 secondPosSet = true;
231 if (secondPos.perp() < firstPos.perp()) std::swap(firstPos, secondPos);
232 } else {
233 if (gpos.perp() < firstPos.perp())
234 firstPos = gpos;
235 else if (secondPos.perp() < gpos.perp())
236 secondPos = gpos;
237 }
238 }
239 }
240 if (precLayers.size() > 2) {
241 if (!phiLayers.empty()) nphiLayers = phiLayers.size();
242 ntrigEtaLayers = etaLayers.size();
243 nprecLayers = precLayers.size();
244 ATH_MSG_DEBUG(" total counts: precision " << static_cast<int>(nprecLayers) << " phi layers " << static_cast<int>(nphiLayers)
245 << " eta trig layers " << static_cast<int>(ntrigEtaLayers)
246 << " associated reco muon " << index << " unique ID " << HepMC::uniqueID(*truthLink)
247 << " truthLink " << truthLink);
248 xAOD::MuonSegment* segment = segmentContainer.push_back(std::make_unique<xAOD::MuonSegment>());
249
250 segment->setNHits(nprecLayers, nphiLayers, ntrigEtaLayers);
251 static const SG::Accessor<ElementLink<xAOD::TruthParticleContainer> >
252 truthParticleLinkAcc("truthParticleLink");
253 truthParticleLinkAcc(*segment) = truthLink;
254 if (chId.is_valid()) { //we should always enter here if we have precision measurements
255 int eta = m_idHelperSvc->stationEta(chId);
256 int sector = m_idHelperSvc->sector(chId);
257 MuonStationIndex::TechnologyIndex technology = m_idHelperSvc->technologyIndex(chId);
258 MuonStationIndex::ChIndex chIndex = m_idHelperSvc->chamberIndex(chId);
259 segment->setIdentifier(sector, chIndex, eta, technology);
260 }
261 if (firstPosSet && secondPosSet) {
262 Amg::Vector3D gpos = (firstPos + secondPos) / 2.;
263 Amg::Vector3D gdir = (firstPos - secondPos).unit();
264 ATH_MSG_DEBUG(" got position : r " << gpos.perp() << " z " << gpos.z() << " and direction: theta " << gdir.theta()
265 << " phi " << gdir.phi());
266 segment->setPosition(gpos.x(), gpos.y(), gpos.z());
267 segment->setDirection(gdir.x(), gdir.y(), gdir.z());
268 }
269 }
270 }
271 return StatusCode::SUCCESS;
272 }
Scalar eta() const
pseudorapidity method
const PlainObject unit() const
This is a plugin that makes Eigen look like CLHEP & defines some convenience methods.
#define ATH_CHECK
Evaluate an expression and check for errors.
#define ATH_MSG_VERBOSE(x)
#define ATH_MSG_DEBUG(x)
value_type push_back(value_type pElem)
Add an element to the end of the collection.
bool is_valid() const
Check if id is in a valid state.
Amg::Vector3D localToGlobalCoords(const Amg::Vector3D &x, const Identifier &id) const
localToGlobalCoords and Transf connect the Gas Gap Frame (defined as a Sensitive Detector) to the Glo...
const CscReadoutElement * getCscReadoutElement(const Identifier &id) const
access via extended identifier (requires unpacking)
ServiceHandle< Muon::IMuonIdHelperSvc > m_idHelperSvc
Handle for the muonIdHelper service.
SG::ReadCondHandleKey< MuonGM::MuonDetectorManager > m_detMgrKey
MuonDetectorManager from the conditions store.
SG::ReadHandleKey< CscSimDataCollection > m_CSC_SDO_TruthNames
SG::ReadHandleKeyArray< MuonSimDataCollection > m_SDO_TruthNames
Keys for all çontainers of muon hit simulation data, classified by detector technology.
void setDirection(float px, float py, float pz)
Sets the direction.
void setNHits(int nPrecisionHits, int nPhiLayers, int nTrigEtaLayers)
Set the number of hits/layers.
void setIdentifier(int sector, ::Muon::MuonStationIndex::ChIndex chamberIndex, int etaIndex, ::Muon::MuonStationIndex::TechnologyIndex technology)
Set the identifier.
void setPosition(float x, float y, float z)
Sets the global position.
Eigen::Matrix< double, 3, 1 > Vector3D
str index
Definition DeMoScan.py:362
int uniqueID(const T &p)
ChIndex chIndex(const std::string &index)
convert ChIndex name string to enum
TechnologyIndex
enum to classify the different layers in the muon spectrometer
constexpr int toInt(const EnumType enumVal)
const std::string & chName(ChIndex index)
convert ChIndex into a string
ChIndex
enum to classify the different chamber layers in the muon spectrometer
const T * get(const ReadCondHandleKey< T > &key, const EventContext &ctx)
Convenience function to retrieve an object given a ReadCondHandleKey.
void swap(ElementLinkVector< DOBJ > &lhs, ElementLinkVector< DOBJ > &rhs)
MuonSegment_v1 MuonSegment
Reference the current persistent version:

◆ declareGaudiProperty()

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

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

Definition at line 156 of file AthCommonDataStore.h.

158 {
160 hndl.value(),
161 hndl.documentation());
162
163 }
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)

◆ declareProperty()

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 }
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>

◆ detStore()

const ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< 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.

◆ evtStore()

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

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

Definition at line 85 of file AthCommonDataStore.h.

◆ execute()

StatusCode Muon::MuonTruthSegmentCreationAlg::execute ( const EventContext & ctx) const
overridevirtual

Definition at line 59 of file MuonTruthSegmentCreationAlg.cxx.

59 {
60
61 const xAOD::TruthParticleContainer* muonTruthContainer{nullptr};
62 ATH_CHECK(SG::get(muonTruthContainer, m_muonTruth, ctx));
63
64 SG::ReadDecorHandle<xAOD::TruthParticleContainer, int> truthOrigin(m_truthOriginKey, ctx);
65 ATH_CHECK(truthOrigin.isPresent());
66
67 // create output container
68 SG::WriteHandle segmentContainer(m_muonTruthSegmentContainerName, ctx);
69 ATH_CHECK(segmentContainer.record(std::make_unique<xAOD::MuonSegmentContainer>(),
70 std::make_unique<xAOD::MuonSegmentAuxContainer>()));
71 ATH_MSG_DEBUG("Recorded MuonSegmentContainer with key: " << segmentContainer.name());
72
73 size_t itr = 0;
74 for (const xAOD::TruthParticle* truthParticle : *muonTruthContainer) {
75
76 const int iOrigin = truthOrigin(*truthParticle);
77 bool goodMuon = bad_origins.find(iOrigin) == bad_origins.end();
78
79 // create segments
80 if (goodMuon) {
81 ElementLink<xAOD::TruthParticleContainer> truthLink(*muonTruthContainer, itr);
82 truthLink.toPersistent();
83
85 ATH_CHECK(fillChamberIdMap(ctx, *truthParticle, ids));
86 ATH_CHECK(createSegments(ctx, truthLink, ids, *segmentContainer));
87 }
88 itr++;
89 }
90
91 ATH_MSG_DEBUG("Registered " << segmentContainer->size() << " truth muon segments ");
92
93 return StatusCode::SUCCESS;
94 }
SG::ReadDecorHandleKey< xAOD::TruthParticleContainer > m_truthOriginKey
StatusCode fillChamberIdMap(const EventContext &ctx, const xAOD::TruthParticle &truthParticle, ChamberIdMap &ids) const
This function uses the 6 vectors, contained in.
std::map< Muon::MuonStationIndex::ChIndex, std::vector< Identifier > > ChamberIdMap
This map contains all the hits corresponding to truth muons classified by chamber layer that recorded...
SG::WriteHandleKey< xAOD::MuonSegmentContainer > m_muonTruthSegmentContainerName
Key for segment container that will be populated with segments.
StatusCode createSegments(const EventContext &ctx, const ElementLink< xAOD::TruthParticleContainer > &truthLink, const ChamberIdMap &ids, xAOD::MuonSegmentContainer &segmentContainer) const
This function performs, for each truth muon, the actual segment creation and stores segments into a n...
SG::ReadHandleKey< xAOD::TruthParticleContainer > m_muonTruth
Key for the truth muon container and muon origin decoration.
int truthOrigin(const U &p)
TruthParticle_v1 TruthParticle
Typedef to implementation.
TruthParticleContainer_v1 TruthParticleContainer
Declare the latest version of the truth particle container.

◆ 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 ( ) const
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 }
97}
An algorithm that can be simultaneously executed in multiple threads.

◆ fillChamberIdMap()

StatusCode Muon::MuonTruthSegmentCreationAlg::fillChamberIdMap ( const EventContext & ctx,
const xAOD::TruthParticle & truthParticle,
ChamberIdMap & ids ) const
private

This function uses the 6 vectors, contained in.

Parameters
hitsCollection,carryingtruth muon hit IDs categorized by detector technology (MDT, TGC, RPC, sTGC, CSC, MM) to (re-)fill the
idsChamberIdMap, which contains all truth muon hit IDs classified by chamber layer. This function retrieves the chamber layer from the hit ID by means of the muonIdHelper tool.
hitsCollectionhit collection
idsempty chamberIdMap

Definition at line 96 of file MuonTruthSegmentCreationAlg.cxx.

98 {
99
100 for (SG::ReadDecorHandle<xAOD::TruthParticleContainer, std::vector<unsigned long long>>& hitCollection : m_truthHitsKeyArray.makeHandles(ctx)){
101 for (const unsigned long long& hit_compID : hitCollection(truthParticle)){
102 const Identifier id{hit_compID};
103 if (m_idHelperSvc->isTgc(id)) { // TGCS should be added to both EIL and EIS
104 if (m_idHelperSvc->phiIndex(id) == PhiIndex::T4) {
105 ids[ChIndex::EIS].push_back(id);
106 ids[ChIndex::EIL].push_back(id);
107 } else {
108 ids[ChIndex::EMS].push_back(id);
109 ids[ChIndex::EML].push_back(id);
110 }
111 } else {
112 ids[m_idHelperSvc->chamberIndex(id) ].push_back(id);
113 }
114 }
115 }
116 return StatusCode::SUCCESS;
117 }
SG::ReadDecorHandleKeyArray< xAOD::TruthParticleContainer, std::vector< unsigned long long > > m_truthHitsKeyArray
Keys of the truth muon decorations that we need to read to (re-)fill the chamberIdMap.

◆ 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 }
virtual bool filterPassed(const EventContext &ctx) const

◆ initialize()

StatusCode Muon::MuonTruthSegmentCreationAlg::initialize ( )
overridevirtual

Definition at line 37 of file MuonTruthSegmentCreationAlg.cxx.

37 {
38 ATH_CHECK(m_muonTruth.initialize());
39
41 ATH_CHECK(m_SDO_TruthNames.initialize());
43
44 ATH_CHECK(m_detMgrKey.initialize());
45 ATH_CHECK(m_idHelperSvc.retrieve());
46
47 ATH_CHECK(m_truthOriginKey.initialize());
48 if(m_idHelperSvc->hasRPC()) {m_truthHitsKeyArray.emplace_back(m_muonTruth, "truthRpcHits");}
49 if(m_idHelperSvc->hasTGC()) {m_truthHitsKeyArray.emplace_back(m_muonTruth, "truthTgcHits");}
50 if(m_idHelperSvc->hasMDT()) {m_truthHitsKeyArray.emplace_back(m_muonTruth, "truthMdtHits");}
51 if(m_idHelperSvc->hasMM()) {m_truthHitsKeyArray.emplace_back(m_muonTruth, "truthMMHits");}
52 if(m_idHelperSvc->hasSTGC()) {m_truthHitsKeyArray.emplace_back(m_muonTruth, "truthStgcHits");}
53 if(m_idHelperSvc->hasCSC()) {m_truthHitsKeyArray.emplace_back(m_muonTruth, "truthCscHits");}
54 ATH_CHECK(m_truthHitsKeyArray.initialize());
55 return StatusCode::SUCCESS;
56 }

◆ 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()

◆ msg()

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

Definition at line 24 of file AthCommonMsg.h.

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

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

◆ 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();
384 }
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > renounce(T &h)

◆ renounceArray()

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

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364 {
366 }

◆ setFilterPassed()

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

Definition at line 100 of file AthCommonReentrantAlgorithm.h.

100 {
102 }
virtual void setFilterPassed(bool state, const EventContext &ctx) const

◆ 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 HypoBase, and InputMakerBase.

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()) {
124 ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
125 << " with CondSvc");
126 }
127 }
128 }
129 return sc;
130}
#define ATH_MSG_ERROR(x)
#define ATH_MSG_WARNING(x)
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override

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

◆ 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) {
313 for (auto k : keys) {
314 k->setOwner(this);
315 }
316 }
317 }

Member Data Documentation

◆ m_CSC_SDO_TruthNames

SG::ReadHandleKey<CscSimDataCollection> Muon::MuonTruthSegmentCreationAlg::m_CSC_SDO_TruthNames {this, "CSCSDOs", "CSC_SDO"}
private

Definition at line 103 of file MuonTruthSegmentCreationAlg.h.

103{this, "CSCSDOs", "CSC_SDO"};

◆ m_detMgrKey

SG::ReadCondHandleKey<MuonGM::MuonDetectorManager> Muon::MuonTruthSegmentCreationAlg::m_detMgrKey
private
Initial value:
{this, "DetectorManagerKey", "MuonDetectorManager",
"Key of input MuonDetectorManager condition data"}

MuonDetectorManager from the conditions store.

Definition at line 109 of file MuonTruthSegmentCreationAlg.h.

109 {this, "DetectorManagerKey", "MuonDetectorManager",
110 "Key of input MuonDetectorManager condition data"};

◆ 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_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_idHelperSvc

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

Handle for the muonIdHelper service.

Definition at line 106 of file MuonTruthSegmentCreationAlg.h.

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

◆ m_muonTruth

SG::ReadHandleKey<xAOD::TruthParticleContainer> Muon::MuonTruthSegmentCreationAlg::m_muonTruth {this, "muonTruth", "MuonTruthParticles"}
private

Key for the truth muon container and muon origin decoration.

Definition at line 90 of file MuonTruthSegmentCreationAlg.h.

90{this, "muonTruth", "MuonTruthParticles"};

◆ m_muonTruthSegmentContainerName

SG::WriteHandleKey<xAOD::MuonSegmentContainer> Muon::MuonTruthSegmentCreationAlg::m_muonTruthSegmentContainerName {this, "MuonTruthSegmentName", "MuonTruthSegments"}
private

Key for segment container that will be populated with segments.

Definition at line 98 of file MuonTruthSegmentCreationAlg.h.

98{this, "MuonTruthSegmentName", "MuonTruthSegments"};

◆ m_SDO_TruthNames

SG::ReadHandleKeyArray<MuonSimDataCollection> Muon::MuonTruthSegmentCreationAlg::m_SDO_TruthNames
private
Initial value:
{
this, "SDOs", {"RPC_SDO", "TGC_SDO", "MDT_SDO" }, "remove NSW by default for now, can always be changed in the configuration"}

Keys for all çontainers of muon hit simulation data, classified by detector technology.

Definition at line 101 of file MuonTruthSegmentCreationAlg.h.

101 {
102 this, "SDOs", {"RPC_SDO", "TGC_SDO", "MDT_SDO" }, "remove NSW by default for now, can always be changed in the configuration"};

◆ m_truthHitsKeyArray

SG::ReadDecorHandleKeyArray<xAOD::TruthParticleContainer, std::vector<unsigned long long> > Muon::MuonTruthSegmentCreationAlg::m_truthHitsKeyArray {this, "truthHitsKeyArray", { }}
private

Keys of the truth muon decorations that we need to read to (re-)fill the chamberIdMap.

Each key corresponds to the decorator containing hit IDs for each muon for a specific detector technology.

Definition at line 95 of file MuonTruthSegmentCreationAlg.h.

95{this, "truthHitsKeyArray", { }};

◆ m_truthOriginKey

SG::ReadDecorHandleKey<xAOD::TruthParticleContainer> Muon::MuonTruthSegmentCreationAlg::m_truthOriginKey {this, "truthOriginKey", m_muonTruth, "truthOrigin"}
private

Definition at line 91 of file MuonTruthSegmentCreationAlg.h.

91{this, "truthOriginKey", m_muonTruth, "truthOrigin"};

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