ATLAS Offline Software
ClusterConversionUtilities.cxx
Go to the documentation of this file.
1 /*
2  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
3 */
4 
6 
9 
11 
15 #include "GeoModelKernel/throwExcept.h"
16 
17 constexpr static double one_over_twelve = 1. / 12.;
18 
19 namespace TrackingUtilities {
20 
21  std::pair<xAOD::MeasVector<3>, xAOD::MeasMatrix<3>> convertHGTD_LocalPosCov(const HGTD_Cluster &cluster) {
22  auto localPos = cluster.localPosition();
23  auto localCov = cluster.localCovariance();
24 
25  const float time = cluster.time();
26  const float timeResolution = cluster.timeResolution();
27 
28  Eigen::Matrix<float,3,1> localPosition = Eigen::Matrix<float,3,1>::Zero();
29  localPosition(0, 0) = localPos.x();
30  localPosition(1, 0) = localPos.y();
31  localPosition(2, 0) = time;
32 
33  Eigen::Matrix<float,3,3> localCovariance = Eigen::Matrix<float,3,3>::Zero();
34  localCovariance(0, 0) = localCov(0, 0);
35  localCovariance(1, 1) = localCov(1, 1);
36  localCovariance(2, 2) = timeResolution * timeResolution;
37 
38  return {localPosition, localCovariance};
39  }
40 
42  const InDetDD::HGTD_DetectorElement& element,
43  xAOD::HGTDCluster& xaodCluster)
44  {
45  IdentifierHash idHash = element.identifyHash();
46 
47  const auto [localPosition, localCovariance] = convertHGTD_LocalPosCov(indetCluster);
48 
49  const auto& RDOs = indetCluster.rdoList();
50  const auto& ToTs = indetCluster.totList();
51 
52  xaodCluster.setMeasurement<3>(idHash, localPosition, localCovariance);
53  xaodCluster.setIdentifier( indetCluster.identify().get_compact() );
54  xaodCluster.setRDOlist(RDOs);
55  xaodCluster.setToTlist(ToTs);
56 
57  return StatusCode::SUCCESS;
58  }
59 
60  std::pair<xAOD::MeasVector<2>, xAOD::MeasMatrix<2>> convertPix_LocalPosCov(const InDet::PixelCluster &cluster) {
61  auto localPos = cluster.localPosition();
62  auto localCov = cluster.localCovariance();
63 
64  Eigen::Matrix<float,2,1> localPosition(localPos.x(), localPos.y());
65 
66  Eigen::Matrix<float,2,2> localCovariance;
67  localCovariance.setZero();
68  localCovariance(0, 0) = localCov(0, 0);
69  localCovariance(1, 1) = localCov(1, 1);
70 
71  return {localPosition, localCovariance};
72  }
73 
75  const InDetDD::SiDetectorElement& element,
76  xAOD::PixelCluster& xaodCluster)
77  {
78  IdentifierHash idHash = element.identifyHash();
79 
80  const auto [localPosition, localCovariance] = convertPix_LocalPosCov(indetCluster);
81 
82  auto globalPos = indetCluster.globalPosition();
83  Eigen::Matrix<float, 3, 1> globalPosition(globalPos.x(), globalPos.y(), globalPos.z());
84 
85  const auto& RDOs = indetCluster.rdoList();
86  const auto& ToTs = indetCluster.totList();
87  const auto& charges = indetCluster.chargeList();
88  const auto& width = indetCluster.width();
89  auto isSplit = indetCluster.isSplit();
90  auto splitProbability1 = indetCluster.splitProbability1();
91  auto splitProbability2 = indetCluster.splitProbability2();
92 
93  xaodCluster.setMeasurement<2>(idHash, localPosition, localCovariance);
94  xaodCluster.setIdentifier( indetCluster.identify().get_compact() );
95  xaodCluster.setRDOlist(RDOs);
96  xaodCluster.globalPosition() = globalPosition;
97  xaodCluster.setToTlist(ToTs);
99  xaodCluster.setChargelist(charges);
101  xaodCluster.setLVL1A(indetCluster.LVL1A());
102  xaodCluster.setChannelsInPhiEta(width.colRow()[0], width.colRow()[1]);
103  xaodCluster.setWidthInEta(static_cast<float>(width.widthPhiRZ()[1]));
104  xaodCluster.setIsSplit(isSplit);
105  xaodCluster.setSplitProbabilities(splitProbability1, splitProbability2);
106 
107  return StatusCode::SUCCESS;
108  }
109 
110  std::pair<xAOD::MeasVector<1>, xAOD::MeasMatrix<1>> convertSCT_LocalPosCov(const InDet::SCT_Cluster &cluster) {
111  const InDetDD::SiDetectorElement& element{*cluster.detectorElement()};
112  auto localPos = cluster.localPosition();
113 
114  float localPosition = 0.f, localCovariance = 0.f;
115  if (element.isBarrel()) {
116  localPosition = localPos.x();
117  localCovariance = element.phiPitch() * element.phiPitch() * one_over_twelve;
118  } else {
119  InDetDD::SiCellId cellId = element.cellIdOfPosition(localPos);
120  const auto* design = dynamic_cast<const InDetDD::StripStereoAnnulusDesign *>(&element.design());
121  if ( design == nullptr ) {
122  THROW_EXCEPTION("Invalid bounds from "<<cluster);
123  }
124  InDetDD::SiLocalPosition localInPolar = design->localPositionOfCellPC(cellId);
125  localPosition = localInPolar.xPhi();
126  localCovariance = design->phiPitchPhi() * design->phiPitchPhi() * one_over_twelve;
127  }
128 
129  return std::make_pair(xAOD::MeasVector<1>{localPosition},
130  xAOD::MeasMatrix<1>{localCovariance});
131  }
132 
134  const InDetDD::SiDetectorElement& element,
135  xAOD::StripCluster& xaodCluster)
136  {
137  IdentifierHash idHash = element.identifyHash();
138 
139  const auto [localPosition, localCovariance] = convertSCT_LocalPosCov(indetCluster);
140 
141  auto globalPos = indetCluster.globalPosition();
142  Eigen::Matrix<float, 3, 1> globalPosition(globalPos.x(), globalPos.y(), globalPos.z());
143 
144  const auto& RDOs = indetCluster.rdoList();
145  const auto& width = indetCluster.width();
146 
147  xaodCluster.setMeasurement<1>(idHash, localPosition, localCovariance);
148  xaodCluster.setIdentifier( indetCluster.identify().get_compact() );
149  xaodCluster.setRDOlist(RDOs);
150  xaodCluster.globalPosition() = globalPosition;
151  xaodCluster.setChannelsInPhi(width.colRow()[0]);
152 
153  return StatusCode::SUCCESS;
154  }
155 
157  const InDetDD::SiDetectorElement& element,
158  const PixelID& pixelID,
159  InDet::PixelCluster*& indetCluster)
160  {
161  const InDetDD::PixelModuleDesign* design(dynamic_cast<const InDetDD::PixelModuleDesign*>(&element.design()));
162  if (design == nullptr) {
163  return StatusCode::FAILURE;
164  }
165 
166  const auto& locPos = xaodCluster.localPosition<2>();
167  Amg::Vector2D localPosition(locPos(0,0), locPos(1,0));
168 
169  InDetDD::SiLocalPosition centroid(localPosition);
170  const Identifier id = element.identifierOfPosition(centroid);
171 
172  const auto& globalPos = xaodCluster.globalPosition();
173  Amg::Vector3D globalPosition(globalPos(0, 0), globalPos(1, 0), globalPos(2, 0));
174 
175  auto errorMatrix = Amg::MatrixX(2,2);
176  errorMatrix.setIdentity();
177  errorMatrix.fillSymmetric(0, 0, xaodCluster.localCovariance<2>()(0, 0));
178  errorMatrix.fillSymmetric(1, 1, xaodCluster.localCovariance<2>()(1, 1));
179 
180  int colmax = std::numeric_limits<int>::min();
181  int rowmax = std::numeric_limits<int>::min();
182  int colmin = std::numeric_limits<int>::max();
183  int rowmin = std::numeric_limits<int>::max();
184 
185  float qRowMin = 0.f;
186  float qRowMax = 0.f;
187  float qColMin = 0.f;
188  float qColMax = 0.f;
189 
190  const std::vector<Identifier>& rod_list_cluster = xaodCluster.rdoList();
191  const std::vector<float>& charge_list_cluster = xaodCluster.chargeList();
192  if (rod_list_cluster.size() != charge_list_cluster.size()) {
193  return StatusCode::FAILURE;
194  }
195 
196  for (std::size_t i(0); i<rod_list_cluster.size(); ++i) {
197  const Identifier& this_rdo = rod_list_cluster.at(i);
198  const float this_charge = charge_list_cluster.at(i);
199 
200  const int row = pixelID.phi_index(this_rdo);
201  if (row > rowmax) {
202  rowmax = row;
203  qRowMax = this_charge;
204  } else if (row == rowmax) {
205  qRowMax += this_charge;
206  }
207 
208  if (row < rowmin) {
209  rowmin = row;
210  qRowMin = this_charge;
211  } else if (row == rowmin) {
212  qRowMin += this_charge;
213  }
214 
215  const int col = pixelID.eta_index(this_rdo);
216  if (col > colmax) {
217  colmax = col;
218  qColMax = this_charge;
219  } else if (col == colmax) {
220  qColMax += this_charge;
221  }
222 
223  if (col < colmin) {
224  colmin = col;
225  qColMin = this_charge;
226  } else if (col == colmin) {
227  qColMin += this_charge;
228  }
229 
230  }
231 
232  // Compute omega for charge interpolation correction (if required)
233  // Two pixels may have charge=0 (very rarely, hopefully)
234  float omegax = -1.f;
235  float omegay = -1.f;
236  if(qRowMin + qRowMax > 0) omegax = qRowMax/(qRowMin + qRowMax);
237  if(qColMin + qColMax > 0) omegay = qColMax/(qColMin + qColMax);
238 
239  double etaWidth = design->widthFromColumnRange(colmin, colmax);
240  double phiWidth = design->widthFromRowRange(rowmin, rowmax);
241  InDet::SiWidth width( Amg::Vector2D(xaodCluster.channelsInPhi(), xaodCluster.channelsInEta()),
242  Amg::Vector2D(phiWidth,etaWidth) );
243 
244  indetCluster = new InDet::PixelCluster(id,
245  localPosition,
246  globalPosition,
247  std::vector<Identifier>(xaodCluster.rdoList()),
248  xaodCluster.lvl1a(),
249  std::vector<int>(xaodCluster.totList()),
250  std::vector<float>(xaodCluster.chargeList()),
251  width,
252  &element,
253  std::move(errorMatrix),
254  omegax,
255  omegay,
256  xaodCluster.isSplit(),
257  xaodCluster.splitProbability1(),
258  xaodCluster.splitProbability2());
259 
260  return StatusCode::SUCCESS;
261  }
262 
264  const InDetDD::SiDetectorElement& element,
265  const SCT_ID& stripID,
266  InDet::SCT_Cluster*& indetCluster,
267  double shift)
268  {
269  bool isBarrel = element.isBarrel();
270  const InDetDD::SCT_ModuleSideDesign* design = nullptr;
271  if (not isBarrel) {
272  design = dynamic_cast<const InDetDD::StripStereoAnnulusDesign*>(&element.design());
273  } else {
274  design = dynamic_cast<const InDetDD::SCT_ModuleSideDesign*>(&element.design());
275  }
276 
277  if (design == nullptr) {
278  return StatusCode::FAILURE;
279  }
280 
281  const auto designShape = design->shape();
282 
283 
284  const auto& rdoList = xaodCluster.rdoList();
285  Identifier id = rdoList.front();
286 
287  const auto& localPos = xaodCluster.localPosition<1>();
288 
289  double pos_x = localPos(0, 0);
290  double pos_y = 0;
291  if (not isBarrel) {
292  const Identifier firstStripId = rdoList.front();
293  int firstStrip = stripID.strip(firstStripId);
294  int stripRow = stripID.row(firstStripId);
295  int clusterSizeInStrips = xaodCluster.channelsInPhi();
296  auto clusterPosition = design->localPositionOfCluster(design->strip1Dim(firstStrip, stripRow), clusterSizeInStrips);
297  pos_x = clusterPosition.xPhi() + shift;
298  pos_y = clusterPosition.xEta();
299  }
300 
301  Amg::Vector2D locpos = Amg::Vector2D( pos_x, pos_y );
302 
303  // Most of the following is taken from what is done in ClusterMakerTool
304  // Need to make this computation instead of using the local pos
305  // with local pos instead some differences w.r.t. reference are observed
306  const auto& firstStrip = stripID.strip(rdoList.front());
307  const auto& lastStrip = stripID.strip(rdoList.back());
308  const auto& row = stripID.row(rdoList.front());
309  const int firstStrip1D = design->strip1Dim (firstStrip, row );
310  const int lastStrip1D = design->strip1Dim( lastStrip, row );
311  const InDetDD::SiCellId cell1(firstStrip1D);
312  const InDetDD::SiCellId cell2(lastStrip1D);
313  const InDetDD::SiLocalPosition firstStripPos( element.rawLocalPositionOfCell(cell1 ));
314  const InDetDD::SiLocalPosition lastStripPos( element.rawLocalPositionOfCell(cell2) );
315  const InDetDD::SiLocalPosition centre( (firstStripPos+lastStripPos) * 0.5 );
316  const double clusterWidth = design->stripPitch() * ( lastStrip - firstStrip + 1 );
317 
318  const std::pair<InDetDD::SiLocalPosition, InDetDD::SiLocalPosition> ends( design->endsOfStrip(centre) );
319  const double stripLength( std::abs(ends.first.xEta() - ends.second.xEta()) );
320 
321  InDet::SiWidth width( Amg::Vector2D(xaodCluster.channelsInPhi(), 1),
322  Amg::Vector2D(clusterWidth, stripLength) );
323 
324  const double col_x = width.colRow().x();
325  const double col_y = width.colRow().y();
326 
327  double scale_factor = 1.;
328  if ( col_x == 1 )
329  scale_factor = 1.05;
330  else if ( col_x == 2 )
331  scale_factor = 0.27;
332 
333  auto errorMatrix = Amg::MatrixX(2,2);
334  errorMatrix.setIdentity();
335  errorMatrix.fillSymmetric(0, 0, scale_factor * scale_factor * width.phiR() * width.phiR() * one_over_twelve);
336  errorMatrix.fillSymmetric(1, 1, width.z() * width.z() / col_y / col_y * one_over_twelve);
337 
338  if( designShape == InDetDD::Trapezoid or
339  designShape == InDetDD::Annulus) {
340  // rotation for endcap SCT
341 
342  // The following is being computed with the local position,
343  // without considering the lorentz shift
344  // So we remove it from the local position
345  Amg::Vector2D local(pos_x - shift, pos_y);
346  double sn = element.sinStereoLocal(local);
347  double sn2 = sn * sn;
348  double cs2 = 1. - sn2;
349  double w = element.phiPitch(local) / element.phiPitch();
350  double v0 = errorMatrix(0,0) * w * w;
351  double v1 = errorMatrix(1,1);
352  errorMatrix.fillSymmetric( 0, 0, cs2 * v0 + sn2 * v1 );
353  errorMatrix.fillSymmetric( 0, 1, sn * std::sqrt(cs2) * (v0 - v1) );
354  errorMatrix.fillSymmetric( 1, 1, sn2 * v0 + cs2 * v1 );
355  }
356 
357  indetCluster = new InDet::SCT_Cluster(id,
358  locpos,
359  std::vector<Identifier>(rdoList),
360  width,
361  &element,
362  std::move(errorMatrix));
363 
364  return StatusCode::SUCCESS;
365  }
366 
368  const InDetDD::HGTD_DetectorElement& element,
369  ::HGTD_Cluster*& indetCluster) {
370 
371  const auto& locPos = xaodCluster.localPosition<3>();
372  Amg::Vector2D localPosition(locPos(0,0), locPos(1,0));
373  float time = locPos(2,0);
374 
375  InDetDD::SiLocalPosition centroid(localPosition);
376  const Identifier id = element.identifierOfPosition(centroid);
377 
378  auto errorMatrix = Amg::MatrixX(2,2);
379  errorMatrix.setIdentity();
380  errorMatrix.fillSymmetric(0, 0, xaodCluster.localCovariance<3>()(0, 0));
381  errorMatrix.fillSymmetric(1, 1, xaodCluster.localCovariance<3>()(1, 1));
382  float time_resolution = std::sqrt(xaodCluster.localCovariance<3>()(2, 2));
383 
384  double etaWidth = 1.3;
385  double phiWidth = 1.3;
386  int channelsPhi = 1;
387  int channelsEta = 1;
388  InDet::SiWidth width( Amg::Vector2D(channelsPhi, channelsEta), Amg::Vector2D(phiWidth, etaWidth) );
389 
390  indetCluster = new ::HGTD_Cluster(id,
391  localPosition,
392  std::vector<Identifier>(xaodCluster.rdoList()),
393  width,
394  &element,
395  std::move(errorMatrix),
396  time,
397  time_resolution,
398  std::vector<int>(xaodCluster.totList()));
399 
400  return StatusCode::SUCCESS;
401  }
402 
403 } // Namespace
404 
405 
InDetDD::SolidStateDetectorElementBase::identifierOfPosition
Identifier identifierOfPosition(const Amg::Vector2D &localPos) const
Full identifier of the cell for a given position: assumes a raw local position (no Lorentz shift)
Definition: SolidStateDetectorElementBase.cxx:217
xAOD::PixelCluster_v1::chargeList
const std::vector< float > & chargeList() const
Returns the list of charges of the channels building the cluster.
InDetDD::SCT_ModuleSideDesign::stripPitch
virtual double stripPitch(const SiLocalPosition &chargePos) const =0
give the strip pitch (dependence on position needed for forward)
TrackingUtilities::convertInDetToXaodCluster
StatusCode convertInDetToXaodCluster(const InDet::PixelCluster &indetCluster, const InDetDD::SiDetectorElement &element, xAOD::PixelCluster &xaodCluster)
Definition: ClusterConversionUtilities.cxx:74
xAOD::PixelCluster_v1::setTotalCharge
void setTotalCharge(float totalCharge)
Sets the total charge.
PixelID::phi_index
int phi_index(const Identifier &id) const
Definition: PixelID.h:654
xAOD::StripCluster_v1::setChannelsInPhi
void setChannelsInPhi(int channelsInPhi)
Sets the dimensions of the cluster in numbers of channels in phi (x)
Amg::MatrixX
Eigen::Matrix< double, Eigen::Dynamic, Eigen::Dynamic > MatrixX
Dynamic Matrix - dynamic allocation.
Definition: EventPrimitives.h:27
InDetDD::PixelModuleDesign::widthFromRowRange
double widthFromRowRange(const int rowMin, const int rowMax) const
Method to calculate phi width from a row range.
Definition: PixelModuleDesign.cxx:179
InDetDD::PixelModuleDesign
Definition: PixelModuleDesign.h:45
xAOD::StripCluster_v1::channelsInPhi
int channelsInPhi() const
Returns the dimensions of the cluster in numbers of channels in phi (x), respectively.
xAOD::PixelCluster_v1::rdoList
const std::vector< Identifier > rdoList() const
Returns the list of identifiers of the channels building the cluster.
Definition: PixelCluster_v1.cxx:42
Amg::Vector2D
Eigen::Matrix< double, 2, 1 > Vector2D
Definition: GeoPrimitives.h:48
InDetDD::SCT_ModuleSideDesign
Definition: SCT_ModuleSideDesign.h:40
max
constexpr double max()
Definition: ap_fixedTest.cxx:33
InDetDD::DetectorDesign::shape
virtual DetectorShape shape() const
Shape of element.
Definition: DetectorDesign.cxx:96
InDet::PixelCluster::isSplit
bool isSplit() const
access mehtods for splitting
Definition: InnerDetector/InDetRecEvent/InDetPrepRawData/InDetPrepRawData/PixelCluster.h:249
min
constexpr double min()
Definition: ap_fixedTest.cxx:26
Trk::PrepRawData::localCovariance
const Amg::MatrixX & localCovariance() const
return const ref to the error matrix
xAOD::HGTDCluster_v1::setRDOlist
void setRDOlist(const std::vector< Identifier > &rdolist)
Sets the list of identifiers of the channels building the cluster.
Definition: HGTDCluster_v1.cxx:15
xAOD::xAODInDetMeasurement::Utilities::computeTotalToT
int computeTotalToT(const SG::AuxElement &cluster)
Definition: Event/xAOD/xAODInDetMeasurement/Root/Utilities.cxx:22
InDetDD::HGTD_DetectorElement
Definition: HGTD_DetectorElement.h:40
InDet::PixelCluster::splitProbability1
double splitProbability1() const
Definition: InnerDetector/InDetRecEvent/InDetPrepRawData/InDetPrepRawData/PixelCluster.h:253
Trk::PrepRawData::rdoList
const std::vector< Identifier > & rdoList() const
return the List of rdo identifiers (pointers)
Identifier::get_compact
value_type get_compact() const
Get the compact id.
xAOD::PixelCluster_v1::setRDOlist
void setRDOlist(const std::vector< Identifier > &rdolist)
Sets the list of identifiers of the channels building the cluster.
Definition: PixelCluster_v1.cxx:32
keylayer_zslicemap.row
row
Definition: keylayer_zslicemap.py:155
xAOD::PixelCluster_v1::setChargelist
void setChargelist(const std::vector< float > &charges)
Sets the list of charges of the channels building the cluster.
TrackingUtilities::convertXaodToInDetCluster
StatusCode convertXaodToInDetCluster(const xAOD::PixelCluster &xaodCluster, const InDetDD::SiDetectorElement &element, const PixelID &pixelID, InDet::PixelCluster *&indetCluster)
Definition: ClusterConversionUtilities.cxx:156
xAOD::HGTDCluster_v1
Definition: HGTDCluster_v1.h:23
xAOD::PixelCluster_v1::setIsSplit
void setIsSplit(bool isSplit)
Sets if the cluster is split or not.
InDetSecVtxTruthMatchUtils::isSplit
bool isSplit(int matchInfo)
Definition: InDetSecVtxTruthMatchTool.h:83
xAOD::PixelCluster_v1::setWidthInEta
void setWidthInEta(float widthInEta)
Sets the width of the cluster in eta (y) direction.
TrackingUtilities::convertHGTD_LocalPosCov
std::pair< xAOD::MeasVector< 3 >, xAOD::MeasMatrix< 3 > > convertHGTD_LocalPosCov(const HGTD_Cluster &cluster)
Definition: ClusterConversionUtilities.cxx:21
InDet::PixelCluster::totList
const std::vector< int > & totList() const
Definition: InnerDetector/InDetRecEvent/InDetPrepRawData/InDetPrepRawData/PixelCluster.h:202
xAOD::PixelCluster_v1::setTotalToT
void setTotalToT(int totalToT)
Sets the total ToT.
InDetDD::SolidStateDetectorElementBase::identifyHash
virtual IdentifierHash identifyHash() const override final
identifier hash (inline)
xAOD::StripCluster_v1::rdoList
const std::vector< Identifier > rdoList() const
Returns the list of identifiers of the channels building the cluster.
Definition: StripCluster_v1.cxx:37
xAOD::PixelCluster_v1::splitProbability2
float splitProbability2() const
HGTD_Cluster.h
Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration.
HGTD_Cluster
Definition: HGTD_Cluster.h:35
xAOD::PixelCluster_v1::setSplitProbabilities
void setSplitProbabilities(float prob1, float prob2)
Sets the splitting probabilities for the cluster.
InDetDD::SiLocalPosition
Definition: SiLocalPosition.h:31
InDetDD::SiLocalPosition::xPhi
double xPhi() const
position along phi direction:
Definition: SiLocalPosition.h:123
InDetDD::SiDetectorElement::phiPitch
double phiPitch() const
Pitch (inline methods)
InDet::PixelCluster::splitProbability2
double splitProbability2() const
Definition: InnerDetector/InDetRecEvent/InDetPrepRawData/InDetPrepRawData/PixelCluster.h:257
xAOD::UncalibratedMeasurement_v1::setIdentifier
void setIdentifier(const DetectorIdentType measId)
Sets the full Identifier of the measurement.
InDetDD::PixelModuleDesign::widthFromColumnRange
double widthFromColumnRange(const int colMin, const int colMax) const
Method to calculate eta width from a column range.
Definition: PixelModuleDesign.cxx:167
parseMapping.v0
def v0
Definition: parseMapping.py:149
SCT_CalibAlgs::firstStrip
@ firstStrip
Definition: SCT_CalibNumbers.h:10
lumiFormat.i
int i
Definition: lumiFormat.py:85
EL::StatusCode
::StatusCode StatusCode
StatusCode definition for legacy code.
Definition: PhysicsAnalysis/D3PDTools/EventLoop/EventLoop/StatusCode.h:22
InDetDD::Annulus
@ Annulus
Definition: DetectorDesign.h:42
SCT_ID::row
int row(const Identifier &id) const
Definition: SCT_ID.h:758
InDetDD::StripStereoAnnulusDesign
Definition: StripStereoAnnulusDesign.h:50
HGTDClusterContainer.h
xAOD::PixelCluster_v1::setToTlist
void setToTlist(const std::vector< int > &tots)
Sets the list of ToT of the channels building the cluster.
HGTDClusterAuxContainer.h
Trk::RIO_OnTrackType::HGTD_Cluster
@ HGTD_Cluster
Definition: RIO_OnTrack.h:63
PixelID::eta_index
int eta_index(const Identifier &id) const
Definition: PixelID.h:660
TrackingUtilities
Definition: ClusterConversionUtilities.h:26
InDet::PixelCluster::LVL1A
int LVL1A() const
Definition: InnerDetector/InDetRecEvent/InDetPrepRawData/InDetPrepRawData/PixelCluster.h:270
xAOD::StripCluster_v1
Definition: StripCluster_v1.h:17
xAOD::HGTDCluster_v1::totList
const std::vector< int > & totList() const
Returns the list of Time Over Threshold of the channels building the cluster.
xAOD::HGTDCluster_v1::setToTlist
void setToTlist(const std::vector< int > &tots)
Sets the list of ToT of the channels building the cluster.
xAOD::UncalibratedMeasurement_v1::localCovariance
ConstMatrixMap< N > localCovariance() const
Returns the local covariance of the measurement.
HGTD_Cluster::timeResolution
virtual float timeResolution() const
Definition: HGTD_Cluster.h:116
InDet::SCT_Cluster
Definition: InnerDetector/InDetRecEvent/InDetPrepRawData/InDetPrepRawData/SCT_Cluster.h:34
ClusterConversionUtilities.h
xAOD::UncalibratedMeasurement_v1::localPosition
ConstVectorMap< N > localPosition() const
Returns the local position of the measurement.
InDet::SiCluster::detectorElement
virtual const InDetDD::SiDetectorElement * detectorElement() const override final
return the detector element corresponding to this PRD The pointer will be zero if the det el is not d...
xAOD::PixelCluster_v1::setLVL1A
void setLVL1A(int lvl1a)
Sets the LVL1 accept.
xAOD::MeasVector
Eigen::Matrix< float, N, 1 > MeasVector
Abrivation of the Matrix & Covariance definitions.
Definition: MeasurementDefs.h:53
InDetDD::local
@ local
Definition: InDetDD_Defs.h:16
xAOD::xAODInDetMeasurement::Utilities::computeTotalCharge
float computeTotalCharge(const SG::AuxElement &cluster)
Definition: Event/xAOD/xAODInDetMeasurement/Root/Utilities.cxx:9
Trk::PrepRawData::identify
Identifier identify() const
return the identifier
Utilities.h
xAOD::PixelCluster_v1::globalPosition
ConstVectorMap< 3 > globalPosition() const
Returns the global position of the pixel cluster.
Definition: PixelCluster_v1.cxx:15
Trk::PrepRawData::localPosition
const Amg::Vector2D & localPosition() const
return the local position reference
InDetDD::SiDetectorElement
Definition: SiDetectorElement.h:109
InDetDD::SiDetectorElement::isBarrel
bool isBarrel() const
xAOD::PixelCluster_v1::totList
const std::vector< int > & totList() const
Returns the list of ToT of the channels building the cluster.
Amg::Vector3D
Eigen::Matrix< double, 3, 1 > Vector3D
Definition: GeoPrimitives.h:47
THROW_EXCEPTION
#define THROW_EXCEPTION(MESSAGE)
Definition: throwExcept.h:10
InDetDD::SCT_ModuleSideDesign::localPositionOfCluster
virtual SiLocalPosition localPositionOfCluster(const SiCellId &cellId, int cluserSize) const =0
InDetDD::SiCellId
Definition: SiCellId.h:29
InDet::PixelCluster
Definition: InnerDetector/InDetRecEvent/InDetPrepRawData/InDetPrepRawData/PixelCluster.h:49
StripStereoAnnulusDesign.h
InDet::SiCluster::globalPosition
const Amg::Vector3D & globalPosition() const
return global position reference
InDet::SiCluster::width
const InDet::SiWidth & width() const
return width class reference
xAOD::StripCluster_v1::setRDOlist
void setRDOlist(const std::vector< Identifier > &rdolist)
Sets the list of identifiers of the channels building the cluster.
Definition: StripCluster_v1.cxx:29
InDetDD::SolidStateDetectorElementBase::rawLocalPositionOfCell
Amg::Vector2D rawLocalPositionOfCell(const SiCellId &cellId) const
Returns position (center) of cell.
Definition: SolidStateDetectorElementBase.cxx:230
xAOD::PixelCluster_v1
Definition: PixelCluster_v1.h:17
xAOD::UncalibratedMeasurement_v1::setMeasurement
void setMeasurement(const DetectorIDHashType idHash, MeasVector< N > locPos, MeasMatrix< N > locCov)
Sets IdentifierHash, local position and local covariance of the measurement.
Base_Fragment.width
width
Definition: Sherpa_i/share/common/Base_Fragment.py:59
CaloSwCorrections.time
def time(flags, cells_name, *args, **kw)
Definition: CaloSwCorrections.py:242
SCT_ID
Definition: SCT_ID.h:68
InDetDD::SCT_ModuleSideDesign::strip1Dim
virtual int strip1Dim(int strip, int row) const override
only relevant for SCT.
Definition: SCT_ModuleSideDesign.h:279
TrackingUtilities::convertSCT_LocalPosCov
std::pair< xAOD::MeasVector< 1 >, xAOD::MeasMatrix< 1 > > convertSCT_LocalPosCov(const InDet::SCT_Cluster &cluster)
Definition: ClusterConversionUtilities.cxx:110
SCT_ID::strip
int strip(const Identifier &id) const
Definition: SCT_ID.h:764
PixelModuleDesign.h
xAOD::PixelCluster_v1::setChannelsInPhiEta
void setChannelsInPhiEta(int channelsInPhi, int channelsInEta)
Sets the dimensions of the cluster in numbers of channels in phi (x) and eta (y) directions.
Definition: PixelCluster_v1.cxx:55
InDetDD::SiDetectorElement::sinStereoLocal
double sinStereoLocal(const Amg::Vector2D &localPos) const
Angle of strip in local frame with respect to the etaAxis.
Definition: SiDetectorElement.cxx:288
SCT_CalibAlgs::lastStrip
@ lastStrip
Definition: SCT_CalibNumbers.h:10
xAOD::phiWidth
phiWidth
Definition: RingSetConf_v1.cxx:612
HGTD_Cluster::totList
virtual const std::vector< int > & totList() const
Definition: HGTD_Cluster.h:118
InDet::SiWidth
Definition: SiWidth.h:25
xAOD::StripCluster_v1::globalPosition
ConstVectorMap< 3 > globalPosition() const
Returns the global position of the strip cluster.
Definition: StripCluster_v1.cxx:15
python.LArCondContChannels.isBarrel
isBarrel
Definition: LArCondContChannels.py:659
TrackingUtilities::convertPix_LocalPosCov
std::pair< xAOD::MeasVector< 2 >, xAOD::MeasMatrix< 2 > > convertPix_LocalPosCov(const InDet::PixelCluster &cluster)
Definition: ClusterConversionUtilities.cxx:60
xAOD::MeasMatrix
Eigen::Matrix< float, N, N > MeasMatrix
Definition: MeasurementDefs.h:55
IdentifierHash
This is a "hash" representation of an Identifier. This encodes a 32 bit index which can be used to lo...
Definition: IdentifierHash.h:25
python.IoTestsLib.w
def w
Definition: IoTestsLib.py:198
xAOD::PixelCluster_v1::channelsInEta
int channelsInEta() const
InDetDD::SiDetectorElement::design
virtual const SiDetectorDesign & design() const override final
access to the local description (inline):
PixelID
Definition: PixelID.h:67
xAOD::HGTDCluster_v1::rdoList
const std::vector< Identifier > rdoList() const
Returns the list of identifiers of the channels building the cluster.
Definition: HGTDCluster_v1.cxx:23
InDetDD::SCT_ModuleSideDesign::endsOfStrip
virtual std::pair< SiLocalPosition, SiLocalPosition > endsOfStrip(const SiLocalPosition &position) const override=0
give the ends of strips
HGTD_Cluster::time
virtual float time() const
Definition: HGTD_Cluster.h:114
NSWL1::centroid
Vertex centroid(const Polygon &p)
Definition: GeoUtils.cxx:59
xAOD::PixelCluster_v1::channelsInPhi
int channelsInPhi() const
Returns the dimensions of the cluster in numbers of channels in phi (x) and eta (y) directions,...
xAOD::PixelCluster_v1::splitProbability1
float splitProbability1() const
Returns the splitting probabilities for the cluster.
InDetDD::Trapezoid
@ Trapezoid
Definition: DetectorDesign.h:42
xAOD::PixelCluster_v1::lvl1a
int lvl1a() const
Return the LVL1 accept.
generate::Zero
void Zero(TH1D *hin)
Definition: generate.cxx:32
InDet::PixelCluster::chargeList
const std::vector< float > & chargeList() const
Definition: InnerDetector/InDetRecEvent/InDetPrepRawData/InDetPrepRawData/PixelCluster.h:210
xAOD::PixelCluster_v1::isSplit
bool isSplit() const
Returns if the cluster is split or not.
Identifier
Definition: IdentifierFieldParser.cxx:14