51 ATH_MSG_DEBUG(
"Executing MeasurementToTrackParticleDecoration...");
88 std::optional<ActsTrk::TrackContainer::ConstTrackProxy> >::
value);
89 std::optional<ActsTrk::TrackContainer::ConstTrackProxy> optional_track = *link_to_track;
91 if ( not optional_track.has_value() ) {
92 ATH_MSG_WARNING(
"Invalid track link for particle " << track_particle->index() <<
". Skipping track..");
96 ATH_MSG_DEBUG(
"Track link found for track particle with index " << track_particle->index());
97 ActsTrk::TrackContainer::ConstTrackProxy
track = optional_track.value();
99 std::vector<int>& regions{measurementRegionHandle(*track_particle)};
100 regions.reserve(
track.nMeasurements());
101 std::vector<int>&
detectors{measurementDetectorHandle(*track_particle)};
103 std::vector<int>&
layers{measurementLayerHandle(*track_particle)};
105 std::vector<int>&
types{measurementTypeHandle(*track_particle)};
107 std::vector<int>& sizesPhi{measurementPhiWidthHandle(*track_particle)};
108 sizesPhi.reserve(
track.nMeasurements());
109 std::vector<int>& sizesEta{measurementEtaWidthHandle(*track_particle)};
110 sizesEta.reserve(
track.nMeasurements());
111 std::vector<float>& residualsLocX{residualLocXhandle(*track_particle)};
112 residualsLocX.reserve(
track.nMeasurements());
113 std::vector<float>& pullsLocX{pullLocXhandle(*track_particle)};
114 pullsLocX.reserve(
track.nMeasurements());
115 std::vector<float>& measurementsLocX{measurementLocXhandle(*track_particle)};
116 measurementsLocX.reserve(
track.nMeasurements());
117 std::vector<float>& trackParametersLocX{trackParameterLocXhandle(*track_particle)};
118 trackParametersLocX.reserve(
track.nMeasurements());
119 std::vector<float>& measurementsLocCovX{measurementLocCovXhandle(*track_particle)};
120 measurementsLocCovX.reserve(
track.nMeasurements());
121 std::vector<float>& trackParametersLocCovX{trackParameterLocCovXhandle(*track_particle)};
122 trackParametersLocCovX.reserve(
track.nMeasurements());
123 std::vector<float>& residualsLocY{residualLocYhandle(*track_particle)};
124 residualsLocY.reserve(
track.nMeasurements());
125 std::vector<float>& pullsLocY{pullLocYhandle(*track_particle)};
126 pullsLocY.reserve(
track.nMeasurements());
127 std::vector<float>& measurementsLocY{measurementLocYhandle(*track_particle)};
128 measurementsLocY.reserve(
track.nMeasurements());
129 std::vector<float>& trackParametersLocY{trackParameterLocYhandle(*track_particle)};
130 trackParametersLocY.reserve(
track.nMeasurements());
131 std::vector<float>& measurementsLocCovY{measurementLocCovYhandle(*track_particle)};
132 measurementsLocCovY.reserve(
track.nMeasurements());
133 std::vector<float>& trackParametersLocCovY{trackParameterLocCovYhandle(*track_particle)};
134 trackParametersLocCovY.reserve(
track.nMeasurements());
136 for (
const auto& state :
track.trackStatesReversed()) {
138 auto flag = state.typeFlags();
140 bool anyHit =
flag.test(Acts::TrackStateFlag::HoleFlag) or
flag.test(Acts::TrackStateFlag::MeasurementFlag);
142 ATH_MSG_DEBUG(
"--- This is not a hit measurement, skipping...");
154 float residualLocX = -1.;
155 float pullLocX = -1.;
156 float measurementLocX = -1.;
157 float trackParameterLocX = -1.;
158 float measurementLocCovX = -1.;
159 float trackParameterLocCovX = -1;
160 float residualLocY = -1.;
161 float pullLocY = -1.;
162 float measurementLocY = -1.;
163 float trackParameterLocY = -1.;
164 float measurementLocCovY = -1.;
165 float trackParameterLocCovY = -1;
166 bool isAnnulusBound =
false;
169 if (
flag.test(Acts::TrackStateFlag::HoleFlag)) {
170 type = MeasurementType::HOLE;
172 }
else if (
flag.test(Acts::TrackStateFlag::OutlierFlag)) {
173 type = MeasurementType::OUTLIER;
176 type = MeasurementType::HIT;
181 if (state.hasReferenceSurface() and state.referenceSurface().associatedDetectorElement()) {
183 if (!detectorElement) {
184 ATH_MSG_WARNING(
"--- TrackState reference surface returned an invalid associated detector element");
188 if (!siliconDetectorElement) {
189 ATH_MSG_WARNING(
"--- TrackState associated detector element is not silicon");
193 const Acts::AnnulusBounds* annulusBounds =
dynamic_cast<const Acts::AnnulusBounds*
>(&(state.referenceSurface().bounds()));
194 isAnnulusBound = annulusBounds ? true :
false;
196 if (siliconDetectorElement) {
198 if (siliconDetectorElement->
isPixel()) {
200 int layer_disk = pixel_id->
layer_disk(detectorIdentifier);
202 if (pixel_id->
barrel_ec(detectorIdentifier) == 0) {
203 if (layer_disk == 0) {
204 detector = Subdetector::INNERMOST_PIXEL;
211 }
else if (siliconDetectorElement->
isSCT()) {
214 region = sct_id->
barrel_ec(detectorIdentifier) == 0 ?
217 }
else ATH_MSG_WARNING(
"--- Unknown detector type - It is not pixel nor strip detecor element!");
219 }
else ATH_MSG_WARNING(
"--- Missing reference surface or associated detector element!");
223 if (
type == MeasurementType::OUTLIER ||
type == MeasurementType::HIT) {
226 if (!state.hasSmoothed() || !state.hasProjector())
230 const auto &calibratedParameters = state.effectiveCalibrated();
231 const auto &calibratedCovariance = state.effectiveCalibratedCovariance();
236 bool evaluateUnbiased = (!
flag.test(Acts::TrackStateFlag::OutlierFlag));
238 if (evaluateUnbiased) {
239 ATH_MSG_DEBUG(
"--- Good for unbiased parameters evaluation!");
240 type = MeasurementType::UNBIASED;
242 if (state.hasUncalibratedSourceLink()) {
243 ATLASUncalibSourceLink sourceLink = state.getUncalibratedSourceLink().template get<ATLASUncalibSourceLink>();
254 ATH_MSG_DEBUG(
"xAOD::UncalibratedMeasurement is neither xAOD::PixelCluster nor xAOD::StripCluster");
259 const auto& [unbiasedParameters, unbiasedCovariance] =
262 measurementLocX = calibratedParameters[Acts::eBoundLoc0];
263 measurementLocCovX = calibratedCovariance(Acts::eBoundLoc0, Acts::eBoundLoc0);
265 if (!isAnnulusBound) {
267 trackParameterLocX = unbiasedParameters[Acts::eBoundLoc0];
268 residualLocX = (measurementLocX - trackParameterLocX) / 1_um;
269 trackParameterLocCovX = unbiasedCovariance(Acts::eBoundLoc0, Acts::eBoundLoc0);
274 float locR = unbiasedParameters[Acts::eBoundLoc0];
275 float covR = unbiasedCovariance(Acts::eBoundLoc0,Acts::eBoundLoc0);
276 float locphi = unbiasedParameters[Acts::eBoundLoc1];
277 float covphi = unbiasedCovariance(Acts::eBoundLoc1,Acts::eBoundLoc1);
278 float covRphi = unbiasedCovariance(Acts::eBoundLoc0,Acts::eBoundLoc1);
280 trackParameterLocX = locphi;
281 residualLocX =
locR * (measurementLocX - trackParameterLocX) / 1_um;
283 trackParameterLocCovX =
locR*
locR*covphi + locphi*locphi*covR + 2*locphi*
locR*covRphi;
286 measurementLocCovX =
locR*
locR * measurementLocCovX;
289 pullLocX =
evaluatePull(residualLocX, measurementLocCovX,
290 trackParameterLocCovX, evaluateUnbiased);
292 if (state.calibratedSize() == 2) {
293 measurementLocY = calibratedParameters[Acts::eBoundLoc1];
294 trackParameterLocY = unbiasedParameters[Acts::eBoundLoc1];
295 residualLocY = (measurementLocY - trackParameterLocY) / 1_um;
296 measurementLocCovY = calibratedCovariance(Acts::eBoundLoc1, Acts::eBoundLoc1);
297 trackParameterLocCovY = unbiasedCovariance(Acts::eBoundLoc1, Acts::eBoundLoc1);
298 pullLocY =
evaluatePull(residualLocY, measurementLocCovY,
299 trackParameterLocCovY, evaluateUnbiased);
306 regions.push_back(region);
310 sizesPhi.push_back(sizePhi);
311 sizesEta.push_back(sizeEta);
312 residualsLocX.push_back(residualLocX);
313 pullsLocX.push_back(pullLocX);
314 measurementsLocX.push_back(measurementLocX);
315 trackParametersLocX.push_back(trackParameterLocX);
316 measurementsLocCovX.push_back(measurementLocCovX);
317 trackParametersLocCovX.push_back(trackParameterLocCovX);
318 residualsLocY.push_back(residualLocY);
319 pullsLocY.push_back(pullLocY);
320 measurementsLocY.push_back(measurementLocY);
321 trackParametersLocY.push_back(trackParameterLocY);
322 measurementsLocCovY.push_back(measurementLocCovY);
323 trackParametersLocCovY.push_back(trackParameterLocCovY);
329 return StatusCode::SUCCESS;
1.8.18