61 {
62
66
69 const InDetDD::SiDetectorElementCollection* stripElements = stripDetEleHandle.
cptr();
70
71 const ActsTrk::MeasurementToTruthParticleAssociation* measToTruth(nullptr);
75 measToTruth = measToTruthHandle.
ptr();
76 }
77
78
79
81 ATH_CHECK(xaod.
record(std::make_unique<xAOD::TrackMeasurementValidationContainer>(),
82 std::make_unique<xAOD::TrackMeasurementValidationAuxContainer>()));
83
85 ElementLink< xAOD::TrackMeasurementValidationContainer > > decorator_measurement_link(
m_trackMeasurement_link, ctx );
86
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 );
101 SG::WriteDecorHandle<xAOD::TrackMeasurementValidationContainer, int> decor_side (
m_measurement_side, ctx );
102
103
104 std::vector<bool> keepClusterCollection {};
107
108
109 for (std::size_t
i(0);
i<stripClusters->
size(); ++
i) {
110 if (not keepClusterCollection[i]) continue;
112
115 ElementLink< xAOD::TrackMeasurementValidationContainer > mlink( measurements,
118 decorator_measurement_link(*cluster) = std::move(mlink);
119
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
139
140
143
144
147
150
152 decor_detectorElementID(*measurement) = hashId;
153 decor_waferID(*measurement) = waferId.
get_compact();
155 decor_layer(*measurement) =
m_stripHelper->layer_disk(waferId);
157 decor_sizeZ(*measurement) = 0;
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;
167 }
168
169
170
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;
180
181
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());
199 }
200
201 decor_truth_indices(*measurement) = std::move(tp_indices);
202 decor_truth_barcode(*measurement) = std::move(tp_barcodes);
203 ++measurementIndex;
204 }
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.
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
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.
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::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