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

#include <HGTDTruthTrackDecorationAlg.h>

Inheritance diagram for ActsTrk::HGTDTruthTrackDecorationAlg:
Collaboration diagram for ActsTrk::HGTDTruthTrackDecorationAlg:

Classes

struct  TruthTrackExtensionData
 Data structure to hold truth information about the HGTD track extension. More...
struct  ClusterTruthInfo

Public Member Functions

virtual ~HGTDTruthTrackDecorationAlg ()=default
virtual StatusCode initialize () override
virtual StatusCode execute (const EventContext &) const override
virtual StatusCode sysInitialize () override
 Override sysInitialize.
virtual bool isClonable () const override
 Specify if the algorithm is clonable.
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
 Get filter decision:
virtual void setFilterPassed (bool state, const EventContext &ctx) const
 Set filter decision:
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

TruthTrackExtensionData createTruthDecoration (const xAOD::TruthParticle *truthParticle, const typename ActsTrk::TrackContainer::ConstTrackProxy trackProxy, const ActsTrk::MeasurementToTruthParticleAssociation *association_map) const
 Generate TruthTrackExtensionData from extension.
StatusCode isPrimaryExpected (const xAOD::TruthParticle *truthParticle, const xAOD::UncalibratedMeasurementContainer measurementContainer, const ActsTrk::MeasurementToTruthParticleAssociation *association_map, std::vector< char > &isPrimaryExistsVec) const
 Checks if truth particle produced hits at each one of the HGTD layers.
std::size_t getHGTDLayerIndex (const Acts::GeometryIdentifier &geoID) const
 returns the index of HGTD layer where surfaces lies.
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::TrackParticleContainerm_trackParticleContainerName {this, "TrackParticleContainerName", "", "Name of the TrackParticle container"}
SG::ReadHandleKey< ActsTrk::TrackContainerm_trackContainerKey {this, "ACTSTracksLocation", "HgtdTracks", "Output track collection (ActsTrk variant)"}
SG::ReadHandleKey< xAOD::UncalibratedMeasurementContainerm_uncalibratedMeasurementContainerKey_HGTD {this, "UncalibratedMeasurementContainerKey_HGTD", "HGTD_Clusters", "input cluster collections for HGTD"}
SG::ReadHandleKey< MeasurementToTruthParticleAssociationm_hgtdClustersToTruth { this, "HgtdClustersToTruthAssociationMap", "HgtdClustersToTruthParticles", "Association map from HGTD measurements to generator particles." }
SG::ReadDecorHandleKey< xAOD::TrackParticleContainerm_hgtdTrackLinkKey {this, "hgtdTrackLink", m_trackParticleContainerName, "hgtdTrackLink", "Link to Hgtd track"}
SG::ReadDecorHandleKey< xAOD::TrackParticleContainerm_truthParticleLinkKey {this, "truthParticleLink", m_trackParticleContainerName, "truthParticleLink", "Link to Truth particle"}
SG::WriteDecorHandleKey< xAOD::TrackParticleContainerm_layerClusterTruthClassKey { this, "HGTD_cluster_truth_class", m_trackParticleContainerName, "HGTD_cluster_truth_class", "Decoration for cluster truth classification" }
SG::WriteDecorHandleKey< xAOD::TrackParticleContainerm_layerClusterShadowedKey { this, "HGTD_cluster_shadowed", m_trackParticleContainerName, "HGTD_cluster_shadowed", "Decoration for shadowed cluster" }
SG::WriteDecorHandleKey< xAOD::TrackParticleContainerm_layerClusterMergedKey { this, "HGTD_cluster_merged", m_trackParticleContainerName, "HGTD_cluster_merged", "Decoration for merged cluster" }
SG::WriteDecorHandleKey< xAOD::TrackParticleContainerm_layerPrimaryExpectedKey { this, "HGTD_primary_expected", m_trackParticleContainerName, "HGTD_primary_expected", "Decoration for primary expected cluster" }
ActsTrk::detail::xAODUncalibMeasSurfAcc m_surfAcc {}
PublicToolHandle< ActsTrk::ITrackingGeometryToolm_trackingGeometryTool {this, "TrackingGeometryTool", ""}
std::unique_ptr< SG::AuxElement::Accessor< int > > m_acc_nHgtdHits
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 42 of file HGTDTruthTrackDecorationAlg.h.

Member Typedef Documentation

◆ StoreGateSvc_t

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

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ ~HGTDTruthTrackDecorationAlg()

virtual ActsTrk::HGTDTruthTrackDecorationAlg::~HGTDTruthTrackDecorationAlg ( )
virtualdefault

Member Function Documentation

◆ createTruthDecoration()

HGTDTruthTrackDecorationAlg::TruthTrackExtensionData ActsTrk::HGTDTruthTrackDecorationAlg::createTruthDecoration ( const xAOD::TruthParticle * truthParticle,
const typename ActsTrk::TrackContainer::ConstTrackProxy trackProxy,
const ActsTrk::MeasurementToTruthParticleAssociation * association_map ) const
private

Generate TruthTrackExtensionData from extension.

Parameters
truthParticle- truth particle associated with track
trackProxy- measurement container with HGTD clusters
association_map- hgtd cluster to truth particles map

Definition at line 172 of file HGTDTruthTrackDecorationAlg.cxx.

175 {
176
178
179 std::vector<int> truthClassPerLayer = {-1, -1, -1, -1};
180 std::vector<char> isShadowedPerLayer = {false, false, false, false};
181 std::vector<char> isMergedPerLayer = {false, false, false, false};
182 std::vector<char> isPrimaryExistsVec = {false, false, false, false};
183
184 for (auto state : trackProxy.trackStatesReversed()) {
185 auto flags = state.typeFlags();
186 if (flags.isMeasurement()) {
187 // Check if this is an HGTD hit
188 const auto& surface = state.referenceSurface();
189 Acts::GeometryIdentifier geoID = surface.geometryId();
190
191 std::size_t layerIndex = getHGTDLayerIndex(geoID);
192 ClusterTruthInfo cluster_truth_info;
193
194 assert( state.hasUncalibratedSourceLink() );
195 auto uncalibMeas = detail::xAODUncalibMeasCalibrator::unpack(state.getUncalibratedSourceLink());
196
197 if (association_map->at(uncalibMeas->index()).empty()) {
198 cluster_truth_info.origin = ActsTrk::ClusterTruthOrigin::SECONDARY;
199 ATH_MSG_DEBUG(" \\__Layer "<< layerIndex << " SECONDARY: hit doesnt have truth particle associated with it");
200 }
201 else if(association_map->at(uncalibMeas->index()).size() == 1) {
202 const xAOD::TruthParticle *truth_particle = association_map->at(uncalibMeas->index()).at(0);
203 if(truthParticle->index() == truth_particle->index()){
204 cluster_truth_info.origin = ActsTrk::ClusterTruthOrigin::TRUTH_PARTICLE;
205 ATH_MSG_DEBUG(" \\__Layer "<< layerIndex<< " TRUTH_PARTICLE: hit matches ITk track");
206 }
207 else{
208 cluster_truth_info.origin = ActsTrk::ClusterTruthOrigin::UNRELATED_PARTICLE;
209 ATH_MSG_DEBUG(" \\__Layer "<< layerIndex<< " UNRELATED_PARTICLE: hit doesnt match ITk track");
210 }
211 }
212 else {
213 cluster_truth_info.is_merged = true;
214 cluster_truth_info.origin = ActsTrk::ClusterTruthOrigin::UNRELATED_PARTICLE;
215 ATH_MSG_DEBUG(" \\__Layer "<< layerIndex<< " MERGED: hit with more than one contributing particle");
216 for (const xAOD::TruthParticle *truth_particle : association_map->at(uncalibMeas->index()) ) {
217 if(truthParticle->index() == truth_particle->index()){
218 cluster_truth_info.origin = ActsTrk::ClusterTruthOrigin::TRUTH_PARTICLE;
219 ATH_MSG_DEBUG(" \\__TRUTH_PARTICLE: hit matches ITk track");
220 }
221 }
222 }
223
224
225 truthClassPerLayer[layerIndex] = (int) cluster_truth_info.origin;
226 isShadowedPerLayer[layerIndex] = cluster_truth_info.is_shadowed;
227 isMergedPerLayer[layerIndex] = cluster_truth_info.is_merged;
228 }
229 }
230
231 // Fill the data structure with results
232 data.truthClassVec = std::move(truthClassPerLayer);
233 data.isShadowedVec = std::move(isShadowedPerLayer);
234 data.isMergedVec = std::move(isMergedPerLayer);
235 data.primaryExistsVec = std::move(isPrimaryExistsVec);
236 return data;
237 }
#define ATH_MSG_DEBUG(x)
std::size_t getHGTDLayerIndex(const Acts::GeometryIdentifier &geoID) const
returns the index of HGTD layer where surfaces lies.
static const xAOD::UncalibratedMeasurement * unpack(const Acts::SourceLink &sl)
Helper method to unpack an Acts source link to an uncalibrated measurement.
TruthParticle_v1 TruthParticle
Typedef to implementation.
Data structure to hold truth information about the HGTD track extension.

◆ 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 ActsTrk::HGTDTruthTrackDecorationAlg::execute ( const EventContext & ctx) const
overridevirtual

Definition at line 60 of file HGTDTruthTrackDecorationAlg.cxx.

61 {
62 ATH_MSG_DEBUG("Executing " << name() << "...");
63
64 const xAOD::TrackParticleContainer* trackParticles {nullptr};
65 ATH_CHECK(SG::get(trackParticles, m_trackParticleContainerName, ctx));
66
67 ATH_MSG_DEBUG("Size of trackParticles collection " << trackParticles->size());
68 // ================================================== //
69 // ============ RETRIEVE MEASUREMENTS =============== //
70 // ================================================== //
71
72 ATH_MSG_DEBUG("Reading input collection with key " << m_uncalibratedMeasurementContainerKey_HGTD.key());
73
74 const xAOD::UncalibratedMeasurementContainer* uncalibratedMeasurementContainer{nullptr};
75 ATH_CHECK(SG::get(uncalibratedMeasurementContainer ,m_uncalibratedMeasurementContainerKey_HGTD, ctx));
76
77 ATH_MSG_DEBUG("Retrieved " << uncalibratedMeasurementContainer->size()
78 << " input elements from key " << m_uncalibratedMeasurementContainerKey_HGTD.key());
79
80
81 // ================================================== //
82 // ============ TRUTH MATCHING BLOCK ================ //
83 // ================================================== //
84
85 SG::ReadDecorHandle<xAOD::TrackParticleContainer, ElementLink<ActsTrk::TrackContainer>> hgtdTrackLink( m_hgtdTrackLinkKey, ctx );
86
87 SG::ReadDecorHandle<xAOD::TrackParticleContainer, ElementLink<xAOD::TruthParticleContainer>> truthParticleLink( m_truthParticleLinkKey, ctx );
88
89 std::array<const ActsTrk::MeasurementToTruthParticleAssociation *,
90 Acts::toUnderlying(xAOD::UncalibMeasType::nTypes)>
91 measurement_to_truth_association_maps{};
92
93 if (!m_hgtdClustersToTruth.key().empty()) {
94 const ActsTrk::MeasurementToTruthParticleAssociation* hgtdClustersToTruthAssociation{};
95 ATH_CHECK(SG::get(hgtdClustersToTruthAssociation, m_hgtdClustersToTruth, ctx));
96 measurement_to_truth_association_maps[Acts::toUnderlying(xAOD::UncalibMeasType::HGTDClusterType)]=hgtdClustersToTruthAssociation;
97 }
98 auto assocSize = [&measurement_to_truth_association_maps](xAOD::UncalibMeasType type) {
99 const ActsTrk::MeasurementToTruthParticleAssociation *assoc = measurement_to_truth_association_maps[Acts::toUnderlying(type)];
100 return assoc ? assoc->size() : 0ul;
101 };
102
103 ATH_MSG_DEBUG("Measurement association entries: " << assocSize(xAOD::UncalibMeasType::PixelClusterType)
104 << " + " << assocSize(xAOD::UncalibMeasType::StripClusterType)
105 << " + " << assocSize(xAOD::UncalibMeasType::HGTDClusterType)
106 );
107
108
109 SG::ReadHandle<ActsTrk::TrackContainer> tracksContainer = SG::makeHandle( m_trackContainerKey, ctx);
110 if (!tracksContainer.isValid()) {
111 ATH_MSG_ERROR("No tracks for key " << m_trackContainerKey.key() );
112 return StatusCode::FAILURE;
113 }
114 // Create WriteDecorHandles for all decorations
115 SG::WriteDecorHandle<xAOD::TrackParticleContainer, std::vector<int>> layerClusterTruthClassHandle(m_layerClusterTruthClassKey, ctx);
116 SG::WriteDecorHandle<xAOD::TrackParticleContainer, std::vector<char>> layerClusterShadowedHandle(m_layerClusterShadowedKey, ctx);
117 SG::WriteDecorHandle<xAOD::TrackParticleContainer, std::vector<char>> layerClusterMergedHandle(m_layerClusterMergedKey, ctx);
118 SG::WriteDecorHandle<xAOD::TrackParticleContainer, std::vector<char>> layerPrimaryExpectedHandle(m_layerPrimaryExpectedKey, ctx);
119
121
122 for (const xAOD::TrackParticle* trackParticle : *trackParticles) {
123
124 // get the truth particle
125 ElementLink<xAOD::TruthParticleContainer> link_to_truth = truthParticleLink(*trackParticle);
126 const xAOD::TruthParticle* truthParticle = nullptr;
127 if (!link_to_truth.isValid()) {
128 layerClusterTruthClassHandle(*trackParticle) = {-1, -1, -1, -1};
129 layerClusterShadowedHandle(*trackParticle) = {false, false, false, false};
130 layerClusterMergedHandle(*trackParticle) = {false, false, false, false};
131 layerPrimaryExpectedHandle(*trackParticle) = {false, false, false, false};
132 ATH_MSG_WARNING("TrackParticle " << trackParticle->index() << ": invalid truth link");
133 continue;
134 }
135 else {
136 truthParticle = *link_to_truth;
137 }
138 // Check if the TrackParticle has a link to an hgtd extension track
139 if(hgtdTrackLink.isAvailable()){
140 ElementLink<ActsTrk::TrackContainer> link_to_track = hgtdTrackLink(*trackParticle);
141 if (!link_to_track.isValid()) {
142 ATH_MSG_DEBUG("TrackParticle " << trackParticle->index() << ": invalid extension link");
143 data.truthClassVec = {-1, -1, -1, -1};
144 data.isShadowedVec = {false, false, false, false};
145 data.isMergedVec = {false, false, false, false};
146 }
147 else{
148 std::optional<ActsTrk::TrackContainer::ConstTrackProxy> optional_track = *link_to_track;
149 const ActsTrk::TrackContainer::ConstTrackProxy& track = optional_track.value();
150 ATH_MSG_DEBUG("TrackParticle " << trackParticle->index() << ": extension with chi2 " << track.chi2() << " and nHits " << track.nMeasurements());
151 data = createTruthDecoration(truthParticle, track, measurement_to_truth_association_maps[Acts::toUnderlying(xAOD::UncalibMeasType::HGTDClusterType)]);
152 }
153 }
154
155
156 ATH_CHECK(isPrimaryExpected(truthParticle,
157 *uncalibratedMeasurementContainer,
158 measurement_to_truth_association_maps[Acts::toUnderlying(xAOD::UncalibMeasType::HGTDClusterType)],
159 data.primaryExistsVec));
160
161 layerClusterTruthClassHandle(*trackParticle) = data.truthClassVec;
162 layerClusterShadowedHandle(*trackParticle) = data.isShadowedVec;
163 layerClusterMergedHandle(*trackParticle) = data.isMergedVec;
164 layerPrimaryExpectedHandle(*trackParticle) = data.primaryExistsVec;
165
166 }
167 ATH_MSG_DEBUG("Finshed truth matching");
168
169 return StatusCode::SUCCESS;
170 }
#define ATH_CHECK
Evaluate an expression and check for errors.
#define ATH_MSG_ERROR(x)
#define ATH_MSG_WARNING(x)
SG::ReadHandleKey< xAOD::TrackParticleContainer > m_trackParticleContainerName
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer > m_layerClusterShadowedKey
SG::ReadHandleKey< MeasurementToTruthParticleAssociation > m_hgtdClustersToTruth
SG::ReadHandleKey< xAOD::UncalibratedMeasurementContainer > m_uncalibratedMeasurementContainerKey_HGTD
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer > m_layerClusterMergedKey
StatusCode isPrimaryExpected(const xAOD::TruthParticle *truthParticle, const xAOD::UncalibratedMeasurementContainer measurementContainer, const ActsTrk::MeasurementToTruthParticleAssociation *association_map, std::vector< char > &isPrimaryExistsVec) const
Checks if truth particle produced hits at each one of the HGTD layers.
SG::ReadHandleKey< ActsTrk::TrackContainer > m_trackContainerKey
SG::ReadDecorHandleKey< xAOD::TrackParticleContainer > m_truthParticleLinkKey
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer > m_layerPrimaryExpectedKey
TruthTrackExtensionData createTruthDecoration(const xAOD::TruthParticle *truthParticle, const typename ActsTrk::TrackContainer::ConstTrackProxy trackProxy, const ActsTrk::MeasurementToTruthParticleAssociation *association_map) const
Generate TruthTrackExtensionData from extension.
SG::ReadDecorHandleKey< xAOD::TrackParticleContainer > m_hgtdTrackLinkKey
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer > m_layerClusterTruthClassKey
size_type size() const noexcept
Returns the number of elements in the collection.
virtual bool isValid() override final
Can the handle be successfully dereferenced?
const T * get(const ReadCondHandleKey< T > &key, const EventContext &ctx)
Convenience function to retrieve an object given a ReadCondHandleKey.
SG::ReadCondHandle< T > makeHandle(const SG::ReadCondHandleKey< T > &key, const EventContext &ctx=Gaudi::Hive::currentContext())
unsigned long ul
TrackParticle_v1 TrackParticle
Reference the current persistent version:
UncalibMeasType
Define the type of the uncalibrated measurement.
TrackParticleContainer_v1 TrackParticleContainer
Definition of the current "TrackParticle container version".
UncalibratedMeasurementContainer_v1 UncalibratedMeasurementContainer
Define the version of the uncalibrated measurement 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 & AthCommonAlgorithm< 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 89 of file AthCommonAlgorithm.cxx.

54{
55 // If we didn't find any symlinks to add, just return the collection
56 // from the base class. Otherwise, return the extended collection.
57 if (!m_extendedExtraObjects.empty()) {
59 }
61}
Common base class for algorithms.

◆ filterPassed()

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

Get filter decision:

Definition at line 93 of file AthCommonAlgorithm.h.

93 {
94 return execState( ctx ).filterPassed();
95 }
virtual bool filterPassed(const EventContext &ctx) const
Get filter decision:

◆ getHGTDLayerIndex()

std::size_t ActsTrk::HGTDTruthTrackDecorationAlg::getHGTDLayerIndex ( const Acts::GeometryIdentifier & geoID) const
private

returns the index of HGTD layer where surfaces lies.

This index is used at to locate where in the vectors of TrackExtensionData the hit information should be written Returns 99 if surface is outiside of HGTD.

Parameters
geoID- surface geometry id

Definition at line 266 of file HGTDTruthTrackDecorationAlg.cxx.

266 {
267 // Get volume and layer ID
268 std::uint32_t volume = geoID.volume();
269 std::uint32_t layer = geoID.layer();
270
271 // Check if we're in the positive or negative endcap
272 bool isPositiveEndcap = (volume == 25);
273 bool isNegativeEndcap = (volume == 2);
274
275 // Different mapping for different sides to maintain consistent physical ordering
276 if (isPositiveEndcap) {
277 // Mapping for positive endcap
278 switch(layer) {
279 case 2: return 0; // First HGTD layer (closest to IP)
280 case 4: return 1; // Second HGTD layer
281 case 6: return 2; // Third HGTD layer
282 case 8: return 3; // Fourth HGTD layer (farthest from IP)
283 default: return 99; // Invalid layer
284 }
285 } else if (isNegativeEndcap) {
286 // Mapping for negative endcap - potentially different ordering
287 switch(layer) {
288 case 2: return 3;
289 case 4: return 2;
290 case 6: return 1;
291 case 8: return 0;
292 default: return 99; // Invalid layer
293 }
294 } else {
295 return 99; // Not an HGTD volume
296 }
297 }
@ layer
Definition HitInfo.h:79

◆ initialize()

StatusCode ActsTrk::HGTDTruthTrackDecorationAlg::initialize ( )
overridevirtual

Definition at line 37 of file HGTDTruthTrackDecorationAlg.cxx.

38 {
39 ATH_MSG_DEBUG("Initializing " << name() << "...");
40
42
45 ATH_CHECK( m_layerClusterMergedKey.initialize() );
47 ATH_CHECK( m_trackContainerKey.initialize() );
48 ATH_CHECK( m_hgtdTrackLinkKey.initialize() );
49 ATH_CHECK( m_truthParticleLinkKey.initialize() );
52
53 // Initialize surface accessor
54 m_surfAcc = ActsTrk::detail::xAODUncalibMeasSurfAcc{m_trackingGeometryTool.get()};
56
57 return StatusCode::SUCCESS;
58 }
PublicToolHandle< ActsTrk::ITrackingGeometryTool > m_trackingGeometryTool
ActsTrk::detail::xAODUncalibMeasSurfAcc m_surfAcc

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

virtual bool AthCommonAlgorithm< Gaudi::Algorithm >::isClonable ( ) const
inlineoverridevirtualinherited

Specify if the algorithm is clonable.

Only relevant for non-reentrant algorithms. Actual number of clones needs to be set via the "Cardinality" property.

Reimplemented in AFP_DigiTop, AlgB, AlgT, BCM_Digitization, CscDigitBuilder, CscDigitToCscRDO, G4AtlasAlg, G4RunAlg, HGTD_Digitization, HiveAlgBase, InDet::GNNSeedingTrackMaker, InDet::SCT_Clusterization, InDet::SiSPGNNTrackMaker, InDet::SiSPSeededTrackFinder, InDet::SiTrackerSpacePointFinder, ISF::SimKernelMT, ITk::StripDigitization, ITkPixelCablingAlg, ITkStripCablingAlg, LArHitEMapMaker, LArTTL1Maker, LUCID_DigiTop, LVL1::L1TopoSimulation, MergeCalibHits, MergeGenericMuonSimHitColl, MergeHijingPars, MergeMcEventCollection, MergeTrackRecordCollection, MergeTruthJets, MergeTruthParticles, MuonDigitizer, PileUpMTAlg, PixelDigitization, RoIBResultToxAOD, SCT_ByteStreamErrorsTestAlg, SCT_CablingCondAlgFromCoraCool, SCT_CablingCondAlgFromText, SCT_ConditionsParameterTestAlg, SCT_ConditionsSummaryTestAlg, SCT_ConfigurationConditionsTestAlg, SCT_Digitization, SCT_FlaggedConditionTestAlg, SCT_LinkMaskingTestAlg, SCT_MajorityConditionsTestAlg, SCT_ModuleVetoTestAlg, SCT_MonitorConditionsTestAlg, SCT_PrepDataToxAOD, SCT_RawDataToxAOD, SCT_ReadCalibChipDataTestAlg, SCT_ReadCalibDataTestAlg, SCT_RODVetoTestAlg, SCT_SensorsTestAlg, SCT_SiliconConditionsTestAlg, SCT_StripVetoTestAlg, SCT_TdaqEnabledTestAlg, SCT_TestCablingAlg, SCTEventFlagWriter, SCTRawDataProvider, SCTSiLorentzAngleTestAlg, SCTSiPropertiesTestAlg, SGInputLoader, Simulation::BeamEffectsAlg, TileHitVecToCnt, TileMuonFitter, TilePulseForTileMuonReceiver, TileRawChannelMaker, TRTDigitization, and ZDC_DigiTop.

Definition at line 68 of file AthCommonAlgorithm.h.

68 {
69 return true;
70 }

◆ isPrimaryExpected()

StatusCode ActsTrk::HGTDTruthTrackDecorationAlg::isPrimaryExpected ( const xAOD::TruthParticle * truthParticle,
const xAOD::UncalibratedMeasurementContainer measurementContainer,
const ActsTrk::MeasurementToTruthParticleAssociation * association_map,
std::vector< char > & isPrimaryExistsVec ) const
private

Checks if truth particle produced hits at each one of the HGTD layers.

Parameters
truthParticle- truth particle associated with track
measurementContainer- measurement container with HGTD clusters
association_map- hgtd cluster to truth particles map
isPrimaryExistsVec- vector to be return with the information about the

Definition at line 239 of file HGTDTruthTrackDecorationAlg.cxx.

243 {
244
245 isPrimaryExistsVec = {false, false, false, false};
246
247 for(auto uncalibMeas: measurementContainer) {
248
249 auto measurementTruthParticles = association_map->at(uncalibMeas->index());
250 const Acts::Surface* surface = m_surfAcc.get(uncalibMeas);
251 Acts::GeometryIdentifier geoID = surface->geometryId();
252 std::size_t layerIndex = getHGTDLayerIndex(geoID);
253 if(measurementTruthParticles.size() > 0){
254 for(auto measTruthParticle : measurementTruthParticles){
255 if ( truthParticle->index() == measTruthParticle->index()){
256 isPrimaryExistsVec[layerIndex] = true;
257 ATH_MSG_DEBUG(" \\__HIT Exepected at " << layerIndex);
258 }
259 }
260 }
261 }
262
263 return StatusCode::SUCCESS;
264 }

◆ 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 AthCommonAlgorithm< Gaudi::Algorithm >::setFilterPassed ( bool state,
const EventContext & ctx ) const
inlinevirtualinherited

Set filter decision:

Reimplemented in AthFilterAlgorithm.

Definition at line 99 of file AthCommonAlgorithm.h.

99 {
101 }
virtual void setFilterPassed(bool state, const EventContext &ctx) const
Set filter decision:

◆ sysExecute()

StatusCode AthCommonAlgorithm< 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.

Reimplemented in AthAnalysisAlgorithm.

Definition at line 80 of file AthCommonAlgorithm.cxx.

41{
42 return BaseAlg::sysExecute (ctx);
43}

◆ sysInitialize()

StatusCode AthCommonAlgorithm< 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 AthAnalysisAlgorithm, AthFilterAlgorithm, AthHistogramAlgorithm, HypoBase, InputMakerBase, and PyAthena::Alg.

Definition at line 60 of file AthCommonAlgorithm.cxx.

71 {
73
74 if (sc.isFailure()) {
75 return sc;
76 }
77
78 ServiceHandle<ICondSvc> cs("CondSvc",name());
79 for (auto h : outputHandles()) {
80 if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
81 // do this inside the loop so we don't create the CondSvc until needed
82 if ( cs.retrieve().isFailure() ) {
83 ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
85 }
86 if (cs->regHandle(this,*h).isFailure()) {
88 ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
89 << " with CondSvc");
90 }
91 }
92 }
93 return sc;
94}
virtual StatusCode sysInitialize() override
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_acc_nHgtdHits

std::unique_ptr<SG::AuxElement::Accessor<int> > ActsTrk::HGTDTruthTrackDecorationAlg::m_acc_nHgtdHits
private

Definition at line 72 of file HGTDTruthTrackDecorationAlg.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_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 AthCommonAlgorithm< Gaudi::Algorithm >::m_extendedExtraObjects
privateinherited

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

Empty if no symlinks were found.

Definition at line 108 of file AthCommonAlgorithm.h.

◆ m_hgtdClustersToTruth

SG::ReadHandleKey<MeasurementToTruthParticleAssociation> ActsTrk::HGTDTruthTrackDecorationAlg::m_hgtdClustersToTruth { this, "HgtdClustersToTruthAssociationMap", "HgtdClustersToTruthParticles", "Association map from HGTD measurements to generator particles." }
private

Definition at line 58 of file HGTDTruthTrackDecorationAlg.h.

58{ this, "HgtdClustersToTruthAssociationMap", "HgtdClustersToTruthParticles", "Association map from HGTD measurements to generator particles." };

◆ m_hgtdTrackLinkKey

SG::ReadDecorHandleKey<xAOD::TrackParticleContainer> ActsTrk::HGTDTruthTrackDecorationAlg::m_hgtdTrackLinkKey {this, "hgtdTrackLink", m_trackParticleContainerName, "hgtdTrackLink", "Link to Hgtd track"}
private

Definition at line 60 of file HGTDTruthTrackDecorationAlg.h.

60{this, "hgtdTrackLink", m_trackParticleContainerName, "hgtdTrackLink", "Link to Hgtd track"};

◆ m_layerClusterMergedKey

SG::WriteDecorHandleKey<xAOD::TrackParticleContainer> ActsTrk::HGTDTruthTrackDecorationAlg::m_layerClusterMergedKey { this, "HGTD_cluster_merged", m_trackParticleContainerName, "HGTD_cluster_merged", "Decoration for merged cluster" }
private

Definition at line 66 of file HGTDTruthTrackDecorationAlg.h.

66{ this, "HGTD_cluster_merged", m_trackParticleContainerName, "HGTD_cluster_merged", "Decoration for merged cluster" };

◆ m_layerClusterShadowedKey

SG::WriteDecorHandleKey<xAOD::TrackParticleContainer> ActsTrk::HGTDTruthTrackDecorationAlg::m_layerClusterShadowedKey { this, "HGTD_cluster_shadowed", m_trackParticleContainerName, "HGTD_cluster_shadowed", "Decoration for shadowed cluster" }
private

Definition at line 65 of file HGTDTruthTrackDecorationAlg.h.

65{ this, "HGTD_cluster_shadowed", m_trackParticleContainerName, "HGTD_cluster_shadowed", "Decoration for shadowed cluster" };

◆ m_layerClusterTruthClassKey

SG::WriteDecorHandleKey<xAOD::TrackParticleContainer> ActsTrk::HGTDTruthTrackDecorationAlg::m_layerClusterTruthClassKey { this, "HGTD_cluster_truth_class", m_trackParticleContainerName, "HGTD_cluster_truth_class", "Decoration for cluster truth classification" }
private

Definition at line 64 of file HGTDTruthTrackDecorationAlg.h.

64{ this, "HGTD_cluster_truth_class", m_trackParticleContainerName, "HGTD_cluster_truth_class", "Decoration for cluster truth classification" };

◆ m_layerPrimaryExpectedKey

SG::WriteDecorHandleKey<xAOD::TrackParticleContainer> ActsTrk::HGTDTruthTrackDecorationAlg::m_layerPrimaryExpectedKey { this, "HGTD_primary_expected", m_trackParticleContainerName, "HGTD_primary_expected", "Decoration for primary expected cluster" }
private

Definition at line 67 of file HGTDTruthTrackDecorationAlg.h.

67{ this, "HGTD_primary_expected", m_trackParticleContainerName, "HGTD_primary_expected", "Decoration for primary expected cluster" };

◆ m_surfAcc

ActsTrk::detail::xAODUncalibMeasSurfAcc ActsTrk::HGTDTruthTrackDecorationAlg::m_surfAcc {}
private

Definition at line 69 of file HGTDTruthTrackDecorationAlg.h.

69{};

◆ m_trackContainerKey

SG::ReadHandleKey<ActsTrk::TrackContainer> ActsTrk::HGTDTruthTrackDecorationAlg::m_trackContainerKey {this, "ACTSTracksLocation", "HgtdTracks", "Output track collection (ActsTrk variant)"}
private

Definition at line 56 of file HGTDTruthTrackDecorationAlg.h.

56{this, "ACTSTracksLocation", "HgtdTracks", "Output track collection (ActsTrk variant)"};

◆ m_trackingGeometryTool

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

Definition at line 70 of file HGTDTruthTrackDecorationAlg.h.

70{this, "TrackingGeometryTool", ""};

◆ m_trackParticleContainerName

SG::ReadHandleKey<xAOD::TrackParticleContainer> ActsTrk::HGTDTruthTrackDecorationAlg::m_trackParticleContainerName {this, "TrackParticleContainerName", "", "Name of the TrackParticle container"}
private

Definition at line 55 of file HGTDTruthTrackDecorationAlg.h.

55{this, "TrackParticleContainerName", "", "Name of the TrackParticle container"};

◆ m_truthParticleLinkKey

SG::ReadDecorHandleKey<xAOD::TrackParticleContainer> ActsTrk::HGTDTruthTrackDecorationAlg::m_truthParticleLinkKey {this, "truthParticleLink", m_trackParticleContainerName, "truthParticleLink", "Link to Truth particle"}
private

Definition at line 61 of file HGTDTruthTrackDecorationAlg.h.

61{this, "truthParticleLink", m_trackParticleContainerName, "truthParticleLink", "Link to Truth particle"};

◆ m_uncalibratedMeasurementContainerKey_HGTD

SG::ReadHandleKey<xAOD::UncalibratedMeasurementContainer> ActsTrk::HGTDTruthTrackDecorationAlg::m_uncalibratedMeasurementContainerKey_HGTD {this, "UncalibratedMeasurementContainerKey_HGTD", "HGTD_Clusters", "input cluster collections for HGTD"}
private

Definition at line 57 of file HGTDTruthTrackDecorationAlg.h.

57{this, "UncalibratedMeasurementContainerKey_HGTD", "HGTD_Clusters", "input cluster collections for HGTD"};

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