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

#include <StripClusterTruthDecoratorAlg.h>

Inheritance diagram for ActsTrk::StripClusterTruthDecoratorAlg:
Collaboration diagram for ActsTrk::StripClusterTruthDecoratorAlg:

Public Member Functions

 StripClusterTruthDecoratorAlg (const std::string &name, ISvcLocator *pSvcLocator)
virtual ~StripClusterTruthDecoratorAlg ()=default
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 labelMeasurementToKeep (const EventContext &ctx, const xAOD::StripClusterContainer &clusters, std::vector< bool > &labels) const
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::StripClusterContainerm_clustercontainer_key {this,"ClusterContainer", "", "Input Strip Cluster container"}
SG::ReadHandleKey< MeasurementToTruthParticleAssociationm_associationMap_key {this,"AssociationMapOut", "", "Association map between measurements and truth particles"}
SG::ReadCondHandleKey< InDetDD::SiDetectorElementCollectionm_stripDetEleCollKey {this, "StripDetectorElements", "ITkStripDetectorElementCollection"}
SG::ReadHandleKeyArray< xAOD::TrackParticleContainerm_trackParticlesKey {this, "TrackParticles", {}, "Input xAOD::TrackParticles"}
SG::WriteHandleKey< xAOD::TrackMeasurementValidationContainerm_write_xaod_key {this,"MeasurementContainer","", "Output Strip Validation Clusters"}
SG::WriteDecorHandleKey< xAOD::StripClusterContainerm_trackMeasurement_link {this, "MeasurementLink", m_clustercontainer_key, "validationMeasurementLink"}
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainerm_measurement_truth_indices {this, "MeasurementTruthIndices", m_write_xaod_key, "truth_index"}
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainerm_measurement_truth_barcodes {this, "MeasurementTruthBarcode", m_write_xaod_key, "truth_barcode"}
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainerm_measurement_detectorElementID {this, "MeasurementDetectorElementID", m_write_xaod_key, "detectorElementID"}
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainerm_measurement_waferID {this, "MeasurementWaferID", m_write_xaod_key, "waferID"}
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainerm_measurement_bec {this, "MeasurementBEC", m_write_xaod_key, "bec"}
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainerm_measurement_layer {this, "MeasurementLayer", m_write_xaod_key, "layer"}
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainerm_measurement_sizePhi {this, "MeasurementSizePhi", m_write_xaod_key, "sizePhi"}
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainerm_measurement_sizeZ {this, "MeasurementSizeZ", m_write_xaod_key, "sizeZ"}
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainerm_measurement_SiWidth {this, "MeasurementSiWidth", m_write_xaod_key, "SiWidth"}
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainerm_measurement_eta_module {this, "MeasurementEtaModule", m_write_xaod_key, "eta_module"}
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainerm_measurement_phi_module {this, "MeasurementPhiModule", m_write_xaod_key, "phi_module"}
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainerm_measurement_omegax {this, "MeasurementOmegaX", m_write_xaod_key, "omegax"}
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainerm_measurement_omegay {this, "MeasurementOmegaY", m_write_xaod_key, "omegay"}
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainerm_measurement_LorentzShift {this, "MeasurementLorentzShift", m_write_xaod_key, "LorentzShift"}
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainerm_measurement_centroid_xphi {this, "MeasurementCentroidXphi", m_write_xaod_key, "centroid_xphi"}
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainerm_measurement_centroid_xeta {this, "MeasurementCentroidXeta", m_write_xaod_key, "centroid_xeta"}
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainerm_measurement_side {this, "MeasurementSide", m_write_xaod_key, "side"}
Gaudi::Property< bool > m_useTruthInfo {this, "UseTruthInfo", true}
Gaudi::Property< bool > m_keepOnlyOnTrackMeasurements {this, "KeepOnlyOnTrackMeasurements", false, "Keep on on-track measurements instead of the full collection"}
const SCT_IDm_stripHelper {nullptr}
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 StripClusterTruthDecoratorAlg.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

◆ StripClusterTruthDecoratorAlg()

ActsTrk::StripClusterTruthDecoratorAlg::StripClusterTruthDecoratorAlg ( const std::string & name,
ISvcLocator * pSvcLocator )

Definition at line 14 of file StripClusterTruthDecoratorAlg.cxx.

15 :
16 AthReentrantAlgorithm(name, pSvcLocator)
17 {}

◆ ~StripClusterTruthDecoratorAlg()

virtual ActsTrk::StripClusterTruthDecoratorAlg::~StripClusterTruthDecoratorAlg ( )
virtualdefault

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}

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

Definition at line 61 of file StripClusterTruthDecoratorAlg.cxx.

61 {
62
63 SG::ReadHandle<xAOD::StripClusterContainer> StripClusterContainer = SG::makeHandle(m_clustercontainer_key,ctx);
65 const xAOD::StripClusterContainer *stripClusters = StripClusterContainer.cptr();
66
67 SG::ReadCondHandle<InDetDD::SiDetectorElementCollection> stripDetEleHandle = SG::makeHandle( m_stripDetEleCollKey, ctx );
68 ATH_CHECK(stripDetEleHandle.isValid());
69 const InDetDD::SiDetectorElementCollection* stripElements = stripDetEleHandle.cptr();
70
71 const ActsTrk::MeasurementToTruthParticleAssociation* measToTruth(nullptr);
72 if (m_useTruthInfo) {
73 SG::ReadHandle<ActsTrk::MeasurementToTruthParticleAssociation> measToTruthHandle = SG::makeHandle(m_associationMap_key,ctx);
74 ATH_CHECK( measToTruthHandle.isValid());
75 measToTruth = measToTruthHandle.ptr();
76 }
77
78 // Setup outputs
79 // Create the xAOD container and its auxiliary store:
80 SG::WriteHandle<xAOD::TrackMeasurementValidationContainer> xaod = SG::makeHandle( m_write_xaod_key,ctx );
81 ATH_CHECK(xaod.record(std::make_unique<xAOD::TrackMeasurementValidationContainer>(),
82 std::make_unique<xAOD::TrackMeasurementValidationAuxContainer>()));
83
84 SG::WriteDecorHandle<xAOD::StripClusterContainer,
85 ElementLink< xAOD::TrackMeasurementValidationContainer > > decorator_measurement_link( m_trackMeasurement_link, ctx );
86
87 SG::WriteDecorHandle<xAOD::TrackMeasurementValidationContainer, std::uint64_t> decor_detectorElementID ( m_measurement_detectorElementID, ctx );
88 SG::WriteDecorHandle<xAOD::TrackMeasurementValidationContainer, int> decor_waferID ( m_measurement_waferID, ctx );
89 SG::WriteDecorHandle<xAOD::TrackMeasurementValidationContainer, int> decor_bec ( m_measurement_bec, ctx );
90 SG::WriteDecorHandle<xAOD::TrackMeasurementValidationContainer, int> decor_layer ( m_measurement_layer, ctx );
91 SG::WriteDecorHandle<xAOD::TrackMeasurementValidationContainer, int> decor_sizePhi ( m_measurement_sizePhi, ctx );
92 SG::WriteDecorHandle<xAOD::TrackMeasurementValidationContainer, int> decor_sizeZ ( m_measurement_sizeZ, ctx );
93 SG::WriteDecorHandle<xAOD::TrackMeasurementValidationContainer, int> decor_SiWidth( m_measurement_SiWidth, ctx );
94 SG::WriteDecorHandle<xAOD::TrackMeasurementValidationContainer, int> decor_eta_module ( m_measurement_eta_module, ctx );
95 SG::WriteDecorHandle<xAOD::TrackMeasurementValidationContainer, int> decor_phi_module ( m_measurement_phi_module, ctx );
96 SG::WriteDecorHandle<xAOD::TrackMeasurementValidationContainer, float> decor_omegax ( m_measurement_omegax, ctx );
97 SG::WriteDecorHandle<xAOD::TrackMeasurementValidationContainer, float> decor_omegay ( m_measurement_omegay, ctx );
98 SG::WriteDecorHandle<xAOD::TrackMeasurementValidationContainer, float> decor_LorentzShift ( m_measurement_LorentzShift, ctx );
99 SG::WriteDecorHandle<xAOD::TrackMeasurementValidationContainer, float> decor_centroid_xphi ( m_measurement_centroid_xphi, ctx );
100 SG::WriteDecorHandle<xAOD::TrackMeasurementValidationContainer, float> decor_centroid_xeta ( m_measurement_centroid_xeta, ctx );
101 SG::WriteDecorHandle<xAOD::TrackMeasurementValidationContainer, int> decor_side ( m_measurement_side, ctx );
102
103 // collection
104 std::vector<bool> keepClusterCollection {};
105 ATH_CHECK( labelMeasurementToKeep(ctx, *stripClusters, keepClusterCollection) );
106 xAOD::TrackMeasurementValidationContainer *measurements = xaod.ptr();
107
108 // loop over collection and convert to xAOD::TrackMeasurementValidation
109 for (std::size_t i(0); i<stripClusters->size(); ++i) {
110 if (not keepClusterCollection[i]) continue;
111 const xAOD::StripCluster* cluster = stripClusters->at(i);
112
113 measurements->push_back( new xAOD::TrackMeasurementValidation() );
114 xAOD::TrackMeasurementValidation* measurement = measurements->back();
115 ElementLink< xAOD::TrackMeasurementValidationContainer > mlink( measurements,
116 measurements->back()->index() );
117 ATH_CHECK( mlink.isValid() );
118 decorator_measurement_link(*cluster) = std::move(mlink);
119
120 xAOD::DetectorIdentType clusterId = cluster->identifier();
121 xAOD::DetectorIDHashType hashId = cluster->identifierHash();
122
123 const InDetDD::SiDetectorElement *element = stripElements->getDetectorElement(hashId);
124 if ( not element ) {
125 ATH_MSG_FATAL( "Invalid strip detector element for hash " << hashId );
126 return StatusCode::FAILURE;
127 }
128
129 const std::vector<Identifier> rdoList = cluster->rdoList();
130 std::vector< std::uint64_t > rdoIdentifierList;
131 rdoIdentifierList.reserve(rdoList.size());
132 for( const Identifier& hitIdentifier : rdoList ){
133 rdoIdentifierList.push_back( hitIdentifier.get_compact() );
134 }
135
136 //Set Identifier
137 measurement->setIdentifier( clusterId );
138 measurement->setRdoIdentifierList(std::move(rdoIdentifierList));
139
140 //Set Global Position
141 auto gpos = cluster->globalPosition();
142 measurement->setGlobalPosition(gpos.x(), gpos.y(), gpos.z());
143
144 // Set local position and error matrix
145 auto locpos = cluster->localPosition<1>();
146 measurement->setLocalPosition(locpos[0], locpos[1]);
147
148 auto localCov = cluster->localCovariance<1>();
149 measurement->setLocalPositionError( localCov(0,0), 0., 0. );
150
151 const Identifier waferId = m_stripHelper->wafer_id(hashId);
152 decor_detectorElementID(*measurement) = hashId;
153 decor_waferID(*measurement) = waferId.get_compact();
154 decor_bec(*measurement) = m_stripHelper->barrel_ec(waferId);
155 decor_layer(*measurement) = m_stripHelper->layer_disk(waferId);
156 decor_sizePhi(*measurement) = cluster->channelsInPhi();
157 decor_sizeZ(*measurement) = 0;
158 decor_SiWidth(*measurement) = cluster->channelsInPhi();
159 decor_eta_module(*measurement) = m_stripHelper->eta_module(waferId);
160 decor_phi_module(*measurement) = m_stripHelper->phi_module(waferId);
161 decor_omegax(*measurement) = 0;
162 decor_omegay(*measurement) = 0;
163 decor_LorentzShift(*measurement) = 0;
164 decor_centroid_xphi(*measurement) = 0;
165 decor_centroid_xeta(*measurement) = 0;
166 decor_side(*measurement) = m_stripHelper->side(waferId);
167 }
168
169
170 // Get a list of all true particle contributing to the cluster
171 if (m_useTruthInfo) {
172 SG::WriteDecorHandle<xAOD::TrackMeasurementValidationContainer, std::vector<unsigned int>> decor_truth_indices( m_measurement_truth_indices, ctx );
173 SG::WriteDecorHandle<xAOD::TrackMeasurementValidationContainer, std::vector<unsigned int>> decor_truth_barcode( m_measurement_truth_barcodes, ctx );
174
175 std::size_t measurementIndex = 0;
176 for (std::size_t i(0); i<stripClusters->size(); ++i) {
177 if (not keepClusterCollection[i]) continue;
178 const xAOD::StripCluster* cluster = stripClusters->at(i);
179 xAOD::TrackMeasurementValidation* measurement = measurements->at(measurementIndex);
180
181 // Check the two match
182 if (cluster->identifier() != measurement->identifier()) {
183 ATH_MSG_ERROR("Cluster and Measurement are not matching!");
184 return StatusCode::FAILURE;
185 }
186
187 if (cluster->index() >= measToTruth->size()) {
188 ATH_MSG_ERROR("PRD index "<< cluster->index() << " not present in the measurement to truth vector with size " << measToTruth->size());
189 return StatusCode::FAILURE;
190 }
191
192 auto tps = measToTruth->at(cluster->index());
193
194 std::vector<unsigned int> tp_indices;
195 std::vector<unsigned int> tp_barcodes;
196 for (auto tp : tps) {
197 tp_indices.push_back(tp->index());
198 tp_barcodes.push_back(HepMC::uniqueID(tp));
199 }
200
201 decor_truth_indices(*measurement) = std::move(tp_indices);
202 decor_truth_barcode(*measurement) = std::move(tp_barcodes);
203 ++measurementIndex;
204 } // loop on clusters
205 }
206
207 ATH_MSG_DEBUG( " recorded StripPrepData objects: size " << measurements->size() );
208 return StatusCode::SUCCESS;
209 }
#define ATH_CHECK
Evaluate an expression and check for errors.
#define ATH_MSG_ERROR(x)
#define ATH_MSG_FATAL(x)
#define ATH_MSG_DEBUG(x)
SG::ReadCondHandleKey< InDetDD::SiDetectorElementCollection > m_stripDetEleCollKey
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainer > m_measurement_phi_module
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainer > m_measurement_SiWidth
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainer > m_measurement_layer
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainer > m_measurement_eta_module
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainer > m_measurement_centroid_xeta
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainer > m_measurement_detectorElementID
SG::ReadHandleKey< xAOD::StripClusterContainer > m_clustercontainer_key
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainer > m_measurement_side
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainer > m_measurement_waferID
SG::WriteHandleKey< xAOD::TrackMeasurementValidationContainer > m_write_xaod_key
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainer > m_measurement_sizePhi
SG::WriteDecorHandleKey< xAOD::StripClusterContainer > m_trackMeasurement_link
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainer > m_measurement_LorentzShift
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainer > m_measurement_sizeZ
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainer > m_measurement_omegay
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainer > m_measurement_centroid_xphi
StatusCode labelMeasurementToKeep(const EventContext &ctx, const xAOD::StripClusterContainer &clusters, std::vector< bool > &labels) const
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainer > m_measurement_omegax
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainer > m_measurement_bec
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainer > m_measurement_truth_barcodes
SG::ReadHandleKey< MeasurementToTruthParticleAssociation > m_associationMap_key
SG::WriteDecorHandleKey< xAOD::TrackMeasurementValidationContainer > m_measurement_truth_indices
const T * back() const
Access the last element in the collection as an rvalue.
const T * at(size_type n) const
Access an element, as an rvalue.
value_type push_back(value_type pElem)
Add an element to the end of the collection.
size_type size() const noexcept
Returns the number of elements in the collection.
value_type get_compact() const
Get the compact id.
const SiDetectorElement * getDetectorElement(const IdentifierHash &hash) const
size_t index() const
Return the index of this element within its container.
const_pointer_type cptr()
const_pointer_type ptr()
Dereference the pointer.
virtual bool isValid() override final
Can the handle be successfully dereferenced?
StatusCode record(std::unique_ptr< T > data)
Record a const object to the store.
pointer_type ptr()
Dereference the pointer.
ConstVectorMap< 3 > globalPosition() const
Returns the global position of the strip cluster.
int channelsInPhi() const
Returns the dimensions of the cluster in numbers of channels in phi (x), respectively.
const std::vector< Identifier > rdoList() const
Returns the list of identifiers of the channels building the cluster.
void setRdoIdentifierList(const std::vector< uint64_t > &rdoIdentifierList)
Sets the list of RDO identifiers.
void setLocalPositionError(float localXError, float localYError, float localXYCorrelation)
Sets the local position error.
void setLocalPosition(float localX, float localY)
Sets the local position.
void setIdentifier(uint64_t identifier)
Sets the identifier.
void setGlobalPosition(float globalX, float globalY, float globalZ)
Sets the global position.
ConstMatrixMap< N > localCovariance() const
Returns the local covariance of the measurement.
DetectorIDHashType identifierHash() const
Returns the IdentifierHash of the measurement (corresponds to the detector element IdentifierHash)
DetectorIdentType identifier() const
Returns the full Identifier of the measurement.
ConstVectorMap< N > localPosition() const
Returns the local position of the measurement.
int uniqueID(const T &p)
SG::ReadCondHandle< T > makeHandle(const SG::ReadCondHandleKey< T > &key, const EventContext &ctx=Gaudi::Hive::currentContext())
TrackMeasurementValidation_v1 TrackMeasurementValidation
Reference the current persistent version:
StripCluster_v1 StripCluster
Define the version of the strip cluster class.
long unsigned int DetectorIdentType
StripClusterContainer_v1 StripClusterContainer
Define the version of the strip cluster container.
TrackMeasurementValidationContainer_v1 TrackMeasurementValidationContainer
Definition of the current "TrackMeasurementValidation container version".
unsigned int DetectorIDHashType
@ detector ID element hash

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

◆ 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 ActsTrk::StripClusterTruthDecoratorAlg::initialize ( )
overridevirtual

Definition at line 20 of file StripClusterTruthDecoratorAlg.cxx.

20 {
21 ATH_MSG_DEBUG("Initialize " << name() << " ...");
22
23 // Read keys
26 ATH_CHECK(m_stripDetEleCollKey.initialize());
27
28 // Tracks only needed if we want on track clusters only
30
31 // Write keys
32 ATH_CHECK(m_write_xaod_key.initialize());
33
34 // Decorator
36
39
41 ATH_CHECK(m_measurement_waferID.initialize());
42 ATH_CHECK(m_measurement_bec.initialize());
43 ATH_CHECK(m_measurement_layer.initialize());
44 ATH_CHECK(m_measurement_sizePhi.initialize());
45 ATH_CHECK(m_measurement_sizeZ.initialize());
46 ATH_CHECK(m_measurement_SiWidth.initialize());
49 ATH_CHECK(m_measurement_omegax.initialize());
50 ATH_CHECK(m_measurement_omegay.initialize());
54 ATH_CHECK(m_measurement_side.initialize());
55
57
58 return StatusCode::SUCCESS;
59 }
SG::ReadHandleKeyArray< xAOD::TrackParticleContainer > m_trackParticlesKey
const ServiceHandle< StoreGateSvc > & detStore() const
retrieve(aClass, aKey=None)
Definition PyKernel.py:110

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

◆ labelMeasurementToKeep()

StatusCode ActsTrk::StripClusterTruthDecoratorAlg::labelMeasurementToKeep ( const EventContext & ctx,
const xAOD::StripClusterContainer & clusters,
std::vector< bool > & labels ) const
private

Definition at line 211 of file StripClusterTruthDecoratorAlg.cxx.

214{
215 labels.clear();
217 labels.resize(clusters.size(), true);
218 return StatusCode::SUCCESS;
219 }
220 labels.resize(clusters.size(), false);
221
222 static const SG::ConstAccessor< ElementLink<ActsTrk::TrackContainer> > decorator_trackLink("actsTrack");
223
224 // get the tracks
225 for (const SG::ReadHandleKey<xAOD::TrackParticleContainer>& trackParticlesKey : m_trackParticlesKey) {
226 SG::ReadHandle<xAOD::TrackParticleContainer> trackParticleHandle = SG::makeHandle( trackParticlesKey, ctx );
227 ATH_CHECK(trackParticleHandle.isValid());
228 const xAOD::TrackParticleContainer* trackParticles = trackParticleHandle.cptr();
229
230 for (const xAOD::TrackParticle* trackParticle : *trackParticles) {
231 // Get the ACTS track object
232 ATH_CHECK( decorator_trackLink.isAvailable(*trackParticle) );
233 ElementLink<ActsTrk::TrackContainer> trackLink = decorator_trackLink(*trackParticle);
234 ATH_CHECK(trackLink.isValid());
235
236 std::optional<ActsTrk::TrackContainer::ConstTrackProxy> optional_track = *trackLink;
237 if ( not optional_track.has_value() ) {
238 ATH_MSG_ERROR("Invalid track link for particle " << trackParticle->index());
239 return StatusCode::FAILURE;
240 }
241 ActsTrk::TrackContainer::ConstTrackProxy track = optional_track.value();
242
243 // loop on track states
244 track.container().trackStateContainer()
245 .visitBackwards(track.tipIndex(),
246 [&labels]
247 (const typename ActsTrk::TrackContainer::ConstTrackStateProxy& state)
248 {
249 auto flags = state.typeFlags();
250 if (not flags.test(Acts::TrackStateFlag::MeasurementFlag) and
251 not flags.test(Acts::TrackStateFlag::OutlierFlag)) return;
252
253 auto sl = state.getUncalibratedSourceLink().template get<ATLASUncalibSourceLink>();
254 if (sl == nullptr) return;
255
256 const xAOD::UncalibratedMeasurement &cluster = getUncalibratedMeasurement(sl);
257 if (cluster.type() != xAOD::UncalibMeasType::StripClusterType) return;
258 labels.at(cluster.index()) = true;
259 });
260 } // loop on tracks
261 } // loop on read handle keys
262
263 return StatusCode::SUCCESS;
264}
const_pointer_type cptr()
Dereference the pointer.
TrackParticle_v1 TrackParticle
Reference the current persistent version:
TrackParticleContainer_v1 TrackParticleContainer
Definition of the current "TrackParticle container version".

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

SG::ReadHandleKey<MeasurementToTruthParticleAssociation> ActsTrk::StripClusterTruthDecoratorAlg::m_associationMap_key {this,"AssociationMapOut", "", "Association map between measurements and truth particles"}
private

Definition at line 56 of file StripClusterTruthDecoratorAlg.h.

56{this,"AssociationMapOut", "", "Association map between measurements and truth particles"};

◆ m_clustercontainer_key

SG::ReadHandleKey<xAOD::StripClusterContainer> ActsTrk::StripClusterTruthDecoratorAlg::m_clustercontainer_key {this,"ClusterContainer", "", "Input Strip Cluster container"}
private

Definition at line 55 of file StripClusterTruthDecoratorAlg.h.

55{this,"ClusterContainer", "", "Input Strip Cluster container"};

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

Gaudi::Property<bool> ActsTrk::StripClusterTruthDecoratorAlg::m_keepOnlyOnTrackMeasurements {this, "KeepOnlyOnTrackMeasurements", false, "Keep on on-track measurements instead of the full collection"}
private

Definition at line 85 of file StripClusterTruthDecoratorAlg.h.

85{this, "KeepOnlyOnTrackMeasurements", false, "Keep on on-track measurements instead of the full collection"};

◆ m_measurement_bec

SG::WriteDecorHandleKey<xAOD::TrackMeasurementValidationContainer> ActsTrk::StripClusterTruthDecoratorAlg::m_measurement_bec {this, "MeasurementBEC", m_write_xaod_key, "bec"}
private

Definition at line 70 of file StripClusterTruthDecoratorAlg.h.

70{this, "MeasurementBEC", m_write_xaod_key, "bec"};

◆ m_measurement_centroid_xeta

SG::WriteDecorHandleKey<xAOD::TrackMeasurementValidationContainer> ActsTrk::StripClusterTruthDecoratorAlg::m_measurement_centroid_xeta {this, "MeasurementCentroidXeta", m_write_xaod_key, "centroid_xeta"}
private

Definition at line 81 of file StripClusterTruthDecoratorAlg.h.

81{this, "MeasurementCentroidXeta", m_write_xaod_key, "centroid_xeta"};

◆ m_measurement_centroid_xphi

SG::WriteDecorHandleKey<xAOD::TrackMeasurementValidationContainer> ActsTrk::StripClusterTruthDecoratorAlg::m_measurement_centroid_xphi {this, "MeasurementCentroidXphi", m_write_xaod_key, "centroid_xphi"}
private

Definition at line 80 of file StripClusterTruthDecoratorAlg.h.

80{this, "MeasurementCentroidXphi", m_write_xaod_key, "centroid_xphi"};

◆ m_measurement_detectorElementID

SG::WriteDecorHandleKey<xAOD::TrackMeasurementValidationContainer> ActsTrk::StripClusterTruthDecoratorAlg::m_measurement_detectorElementID {this, "MeasurementDetectorElementID", m_write_xaod_key, "detectorElementID"}
private

Definition at line 68 of file StripClusterTruthDecoratorAlg.h.

68{this, "MeasurementDetectorElementID", m_write_xaod_key, "detectorElementID"};

◆ m_measurement_eta_module

SG::WriteDecorHandleKey<xAOD::TrackMeasurementValidationContainer> ActsTrk::StripClusterTruthDecoratorAlg::m_measurement_eta_module {this, "MeasurementEtaModule", m_write_xaod_key, "eta_module"}
private

Definition at line 75 of file StripClusterTruthDecoratorAlg.h.

75{this, "MeasurementEtaModule", m_write_xaod_key, "eta_module"};

◆ m_measurement_layer

SG::WriteDecorHandleKey<xAOD::TrackMeasurementValidationContainer> ActsTrk::StripClusterTruthDecoratorAlg::m_measurement_layer {this, "MeasurementLayer", m_write_xaod_key, "layer"}
private

Definition at line 71 of file StripClusterTruthDecoratorAlg.h.

71{this, "MeasurementLayer", m_write_xaod_key, "layer"};

◆ m_measurement_LorentzShift

SG::WriteDecorHandleKey<xAOD::TrackMeasurementValidationContainer> ActsTrk::StripClusterTruthDecoratorAlg::m_measurement_LorentzShift {this, "MeasurementLorentzShift", m_write_xaod_key, "LorentzShift"}
private

Definition at line 79 of file StripClusterTruthDecoratorAlg.h.

79{this, "MeasurementLorentzShift", m_write_xaod_key, "LorentzShift"};

◆ m_measurement_omegax

SG::WriteDecorHandleKey<xAOD::TrackMeasurementValidationContainer> ActsTrk::StripClusterTruthDecoratorAlg::m_measurement_omegax {this, "MeasurementOmegaX", m_write_xaod_key, "omegax"}
private

Definition at line 77 of file StripClusterTruthDecoratorAlg.h.

77{this, "MeasurementOmegaX", m_write_xaod_key, "omegax"};

◆ m_measurement_omegay

SG::WriteDecorHandleKey<xAOD::TrackMeasurementValidationContainer> ActsTrk::StripClusterTruthDecoratorAlg::m_measurement_omegay {this, "MeasurementOmegaY", m_write_xaod_key, "omegay"}
private

Definition at line 78 of file StripClusterTruthDecoratorAlg.h.

78{this, "MeasurementOmegaY", m_write_xaod_key, "omegay"};

◆ m_measurement_phi_module

SG::WriteDecorHandleKey<xAOD::TrackMeasurementValidationContainer> ActsTrk::StripClusterTruthDecoratorAlg::m_measurement_phi_module {this, "MeasurementPhiModule", m_write_xaod_key, "phi_module"}
private

Definition at line 76 of file StripClusterTruthDecoratorAlg.h.

76{this, "MeasurementPhiModule", m_write_xaod_key, "phi_module"};

◆ m_measurement_side

SG::WriteDecorHandleKey<xAOD::TrackMeasurementValidationContainer> ActsTrk::StripClusterTruthDecoratorAlg::m_measurement_side {this, "MeasurementSide", m_write_xaod_key, "side"}
private

Definition at line 82 of file StripClusterTruthDecoratorAlg.h.

82{this, "MeasurementSide", m_write_xaod_key, "side"};

◆ m_measurement_SiWidth

SG::WriteDecorHandleKey<xAOD::TrackMeasurementValidationContainer> ActsTrk::StripClusterTruthDecoratorAlg::m_measurement_SiWidth {this, "MeasurementSiWidth", m_write_xaod_key, "SiWidth"}
private

Definition at line 74 of file StripClusterTruthDecoratorAlg.h.

74{this, "MeasurementSiWidth", m_write_xaod_key, "SiWidth"};

◆ m_measurement_sizePhi

SG::WriteDecorHandleKey<xAOD::TrackMeasurementValidationContainer> ActsTrk::StripClusterTruthDecoratorAlg::m_measurement_sizePhi {this, "MeasurementSizePhi", m_write_xaod_key, "sizePhi"}
private

Definition at line 72 of file StripClusterTruthDecoratorAlg.h.

72{this, "MeasurementSizePhi", m_write_xaod_key, "sizePhi"};

◆ m_measurement_sizeZ

SG::WriteDecorHandleKey<xAOD::TrackMeasurementValidationContainer> ActsTrk::StripClusterTruthDecoratorAlg::m_measurement_sizeZ {this, "MeasurementSizeZ", m_write_xaod_key, "sizeZ"}
private

Definition at line 73 of file StripClusterTruthDecoratorAlg.h.

73{this, "MeasurementSizeZ", m_write_xaod_key, "sizeZ"};

◆ m_measurement_truth_barcodes

SG::WriteDecorHandleKey<xAOD::TrackMeasurementValidationContainer> ActsTrk::StripClusterTruthDecoratorAlg::m_measurement_truth_barcodes {this, "MeasurementTruthBarcode", m_write_xaod_key, "truth_barcode"}
private

Definition at line 66 of file StripClusterTruthDecoratorAlg.h.

66{this, "MeasurementTruthBarcode", m_write_xaod_key, "truth_barcode"};

◆ m_measurement_truth_indices

SG::WriteDecorHandleKey<xAOD::TrackMeasurementValidationContainer> ActsTrk::StripClusterTruthDecoratorAlg::m_measurement_truth_indices {this, "MeasurementTruthIndices", m_write_xaod_key, "truth_index"}
private

Definition at line 65 of file StripClusterTruthDecoratorAlg.h.

65{this, "MeasurementTruthIndices", m_write_xaod_key, "truth_index"};

◆ m_measurement_waferID

SG::WriteDecorHandleKey<xAOD::TrackMeasurementValidationContainer> ActsTrk::StripClusterTruthDecoratorAlg::m_measurement_waferID {this, "MeasurementWaferID", m_write_xaod_key, "waferID"}
private

Definition at line 69 of file StripClusterTruthDecoratorAlg.h.

69{this, "MeasurementWaferID", m_write_xaod_key, "waferID"};

◆ m_stripDetEleCollKey

SG::ReadCondHandleKey<InDetDD::SiDetectorElementCollection> ActsTrk::StripClusterTruthDecoratorAlg::m_stripDetEleCollKey {this, "StripDetectorElements", "ITkStripDetectorElementCollection"}
private

Definition at line 57 of file StripClusterTruthDecoratorAlg.h.

57{this, "StripDetectorElements", "ITkStripDetectorElementCollection"};

◆ m_stripHelper

const SCT_ID* ActsTrk::StripClusterTruthDecoratorAlg::m_stripHelper {nullptr}
private

Definition at line 87 of file StripClusterTruthDecoratorAlg.h.

87{nullptr};

◆ m_trackMeasurement_link

SG::WriteDecorHandleKey< xAOD::StripClusterContainer > ActsTrk::StripClusterTruthDecoratorAlg::m_trackMeasurement_link {this, "MeasurementLink", m_clustercontainer_key, "validationMeasurementLink"}
private

Definition at line 63 of file StripClusterTruthDecoratorAlg.h.

63{this, "MeasurementLink", m_clustercontainer_key, "validationMeasurementLink"};

◆ m_trackParticlesKey

SG::ReadHandleKeyArray< xAOD::TrackParticleContainer > ActsTrk::StripClusterTruthDecoratorAlg::m_trackParticlesKey {this, "TrackParticles", {}, "Input xAOD::TrackParticles"}
private

Definition at line 59 of file StripClusterTruthDecoratorAlg.h.

59{this, "TrackParticles", {}, "Input xAOD::TrackParticles"};

◆ m_useTruthInfo

Gaudi::Property<bool> ActsTrk::StripClusterTruthDecoratorAlg::m_useTruthInfo {this, "UseTruthInfo", true}
private

Definition at line 84 of file StripClusterTruthDecoratorAlg.h.

84{this, "UseTruthInfo", true};

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

◆ m_write_xaod_key

SG::WriteHandleKey<xAOD::TrackMeasurementValidationContainer> ActsTrk::StripClusterTruthDecoratorAlg::m_write_xaod_key {this,"MeasurementContainer","", "Output Strip Validation Clusters"}
private

Definition at line 61 of file StripClusterTruthDecoratorAlg.h.

61{this,"MeasurementContainer","", "Output Strip Validation Clusters"};

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