58 {
59
63
66 const InDetDD::SiDetectorElementCollection* stripElements = stripDetEleHandle.
cptr();
67
68 const ActsTrk::MeasurementToTruthParticleAssociation* measToTruth(nullptr);
72 measToTruth = measToTruthHandle.
ptr();
73 }
74
75
76
78 ATH_CHECK(xaod.
record(std::make_unique<xAOD::TrackMeasurementValidationContainer>(),
79 std::make_unique<xAOD::TrackMeasurementValidationAuxContainer>()));
80
83
85 SG::WriteDecorHandle<xAOD::TrackMeasurementValidationContainer, int> decor_waferID (
m_measurement_waferID, ctx );
86 SG::WriteDecorHandle<xAOD::TrackMeasurementValidationContainer, int> decor_bec (
m_measurement_bec, ctx );
87 SG::WriteDecorHandle<xAOD::TrackMeasurementValidationContainer, int> decor_layer (
m_measurement_layer, ctx );
88 SG::WriteDecorHandle<xAOD::TrackMeasurementValidationContainer, int> decor_sizePhi (
m_measurement_sizePhi, ctx );
89 SG::WriteDecorHandle<xAOD::TrackMeasurementValidationContainer, int> decor_sizeZ (
m_measurement_sizeZ, ctx );
90 SG::WriteDecorHandle<xAOD::TrackMeasurementValidationContainer, int> decor_SiWidth(
m_measurement_SiWidth, ctx );
91 SG::WriteDecorHandle<xAOD::TrackMeasurementValidationContainer, int> decor_eta_module (
m_measurement_eta_module, ctx );
92 SG::WriteDecorHandle<xAOD::TrackMeasurementValidationContainer, int> decor_phi_module (
m_measurement_phi_module, ctx );
93 SG::WriteDecorHandle<xAOD::TrackMeasurementValidationContainer, float> decor_omegax (
m_measurement_omegax, ctx );
94 SG::WriteDecorHandle<xAOD::TrackMeasurementValidationContainer, float> decor_omegay (
m_measurement_omegay, ctx );
98 SG::WriteDecorHandle<xAOD::TrackMeasurementValidationContainer, int> decor_side (
m_measurement_side, ctx );
99
100
101 std::vector<bool> keepClusterCollection {};
104
105
106 for (std::size_t
i(0);
i<stripClusters->
size(); ++
i) {
107 if (not keepClusterCollection[i]) continue;
109
113 measurements->
back()->index() );
115 decorator_measurement_link(*cluster) = std::move(mlink);
116
119
120 const InDetDD::SiDetectorElement *element = stripElements->
getDetectorElement(hashId);
121 if ( not element ) {
122 ATH_MSG_FATAL(
"Invalid strip detector element for hash " << hashId );
123 return StatusCode::FAILURE;
124 }
125
128 std::vector< std::uint64_t > rdoIdentifierList;
129 rdoIdentifierList.reserve(rdoList.size());
130 for( Identifier::value_type hitIdentifierValue : rdoList ){
131 rdoIdentifierList.push_back( hitIdentifierValue );
132 }
133
134
137
138
141
142
145
148
150 decor_detectorElementID(*measurement) = hashId;
151 decor_waferID(*measurement) = waferId.
get_compact();
153 decor_layer(*measurement) =
m_stripHelper->layer_disk(waferId);
155 decor_sizeZ(*measurement) = 0;
157 decor_eta_module(*measurement) =
m_stripHelper->eta_module(waferId);
158 decor_phi_module(*measurement) =
m_stripHelper->phi_module(waferId);
159 decor_omegax(*measurement) = 0;
160 decor_omegay(*measurement) = 0;
161 decor_LorentzShift(*measurement) = 0;
162 decor_centroid_xphi(*measurement) = 0;
163 decor_centroid_xeta(*measurement) = 0;
165 }
166
167
168
170 SG::WriteDecorHandle<xAOD::TrackMeasurementValidationContainer, std::vector<unsigned int>> decor_truth_indices(
m_measurement_truth_indices, ctx );
171 SG::WriteDecorHandle<xAOD::TrackMeasurementValidationContainer, std::vector<unsigned int>> decor_truth_barcode(
m_measurement_truth_barcodes, ctx );
172
173 std::size_t measurementIndex = 0;
174 for (std::size_t
i(0);
i<stripClusters->
size(); ++
i) {
175 if (not keepClusterCollection[i]) continue;
178
179
182 return StatusCode::FAILURE;
183 }
184
185 if (cluster->index() >= measToTruth->size()) {
186 ATH_MSG_ERROR(
"PRD index "<< cluster->index() <<
" not present in the measurement to truth vector with size " << measToTruth->size());
187 return StatusCode::FAILURE;
188 }
189
190 auto tps = measToTruth->at(cluster->index());
191
192 std::vector<unsigned int> tp_indices;
193 std::vector<unsigned int> tp_barcodes;
194 for (auto tp : tps) {
195 tp_indices.push_back(tp->index());
197 }
198
199 decor_truth_indices(*measurement) = std::move(tp_indices);
200 decor_truth_barcode(*measurement) = std::move(tp_barcodes);
201 ++measurementIndex;
202 }
203 }
204
205 ATH_MSG_DEBUG(
" recorded StripPrepData objects: size " << measurements->
size() );
206 return StatusCode::SUCCESS;
207 }
#define ATH_CHECK
Evaluate an expression and check for errors.
ElementLink()
Default constructor.
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
Gaudi::Property< bool > m_useTruthInfo
const SCT_ID * m_stripHelper
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
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.
SG::ConstAccessor< SG::JaggedVecElt< Identifier::value_type > >::element_type rdoList() const
Returns the list of identifiers of the channels building the cluster.
int channelsInPhi() const
Returns the dimensions of the cluster in numbers of channels in phi (x), respectively.
void setRdoIdentifierList(const std::vector< uint64_t > &rdoIdentifierList)
Sets the list of RDO identifiers.
uint64_t identifier() const
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.
SG::ConstAccessor< T, ALLOC > ConstAccessor
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