ATLAS Offline Software
Typedefs | Functions
MuonR4::SegmentFitHelpers Namespace Reference

Typedefs

using State = CalibratedSpacePoint::State
 

Functions

double chiSqTerm (const Amg::Vector3D &posInChamber, const Amg::Vector3D &dirInChamber, const SpacePoint &measurement, MsgStream &msg)
 Calculates the chi2 contribuation to a linear segment line from an uncalibrated measurement. More...
 
double chiSqTermMdt (const Amg::Vector3D &posInChamber, const Amg::Vector3D &dirInChamber, const SpacePoint &measurement, MsgStream &msg)
 Calculates the chi2 contribuation to a linear segment line from an uncalibrated Mdt measurement. More...
 
double chiSqTermStrip (const Amg::Vector3D &posInChamber, const Amg::Vector3D &dirInChamber, const SpacePoint &measurement, MsgStream &msg)
 Calculates the chi2 contribuation to a linear segment line from an uncalibrated strip measurement. More...
 
double chiSqTerm (const Amg::Vector3D &posInChamber, const Amg::Vector3D &dirInChamber, const double timeShift, std::optional< double > arrivalTime, const CalibratedSpacePoint &measurement, MsgStream &msg)
 Calculates the chi2 contribution from the given measurement. More...
 
double chiSqTermMdt (const Amg::Vector3D &posInChamber, const Amg::Vector3D &dirInChamber, const CalibratedSpacePoint &mdtSpacePoint, MsgStream &msg)
 Calculates the chi2 contribution from a mdt space point to the segment line. More...
 
double chiSqTermStrip (const Amg::Vector3D &posInChamber, const Amg::Vector3D &dirInChamber, const double timeShift, std::optional< double > arrivalTime, const CalibratedSpacePoint &strip, MsgStream &msg)
 Calculates the chi2 contribution from a strip measurement to the segment line. More...
 
double chiSqTermBeamspot (const Amg::Vector3D &posInChamber, const Amg::Vector3D &dirInChamber, const CalibratedSpacePoint &beamSpotMeas, MsgStream &msg)
 Calculates the chi2 contribution from an external beam spot constraint. More...
 
std::vector< int > driftSigns (const Amg::Vector3D &posInChamber, const Amg::Vector3D &dirInChamber, const std::vector< const SpacePoint * > &uncalibHits, MsgStream &msg)
 Calculates whether a segment line travereses the tube measurements on the left (-1) or right (1) side of the tube wire. More...
 
std::vector< int > driftSigns (const Amg::Vector3D &posInChamber, const Amg::Vector3D &dirInChamber, const std::vector< std::unique_ptr< CalibratedSpacePoint >> &calibHits, MsgStream &msg)
 Calculates whether a segment line travereses the tube measurements on the left (-1) or right (1) side of the tube wire. More...
 
int driftSign (const Amg::Vector3D &posInChamber, const Amg::Vector3D &dirInChamber, const SpacePoint &uncalibHit, MsgStream &msg)
 Calculates whether a segement line travereses the tube measurement on the left (-1) or right (1) side of the tube wire. More...
 
int driftSign (const Amg::Vector3D &posInChamber, const Amg::Vector3D &dirInChamber, const CalibratedSpacePoint &calibHit, MsgStream &msg)
 Calculates whether a segement line travereses the tube measurement on the left (-1) or right (1) side of the tube wire. More...
 

Typedef Documentation

◆ State

using MuonR4::SegmentFitHelpers::State = typedef CalibratedSpacePoint::State

Definition at line 25 of file SegmentFitHelperFunctions.cxx.

Function Documentation

◆ chiSqTerm() [1/2]

double MuonR4::SegmentFitHelpers::chiSqTerm ( const Amg::Vector3D &  posInChamber,
const Amg::Vector3D &  dirInChamber,
const double  timeShift,
std::optional< double >  arrivalTime,
const CalibratedSpacePoint &  measurement,
MsgStream &  msg 
)

Calculates the chi2 contribution from the given measurement.

Currently, MdtDriftCircles, Rpc & Tgc as well as the Beamspot are supported

Parameters
posInChamberPosition of the segment in the sector frame
dirInChamberDirection of flight of the segment expressed in the sector frame
timeShiftShift from the nominal arrival time, that's calculated as R / c.
arrivalTimeArrival time of the particle at the segment's refernce plane. If the arrival time is passed, the timing information of the measurement is explicitly taken into account.
measurementSpace point to which the chi2 term is calculated
msgReference to the callers msgStream. If the level is VERBOSE, then all relevant parameters are printed

Definition at line 99 of file SegmentFitHelperFunctions.cxx.

102  {
103  double chi2{0.};
104  switch (hit.type()) {
105  case xAOD::UncalibMeasType::MdtDriftCircleType:
106  chi2 = chiSqTermMdt(segPos, segDir, hit, msg);
107  break;
108  case xAOD::UncalibMeasType::RpcStripType:
109  case xAOD::UncalibMeasType::TgcStripType:
110  case xAOD::UncalibMeasType::MMClusterType:
111  case xAOD::UncalibMeasType::sTgcStripType:
112  chi2 = chiSqTermStrip(segPos, segDir, timeOffset, arrivalTime, hit, msg);
113  break;
114  case xAOD::UncalibMeasType::Other:
115  chi2 = chiSqTermBeamspot(segPos,segDir, hit, msg);
116  break;
117  default:
118  if (msg.level() <= MSG::WARNING) {
119  msg<<MSG::WARNING<<__FILE__<<":"<<__LINE__<<" - Unsupported measurement: "
120  <<hit.spacePoint()->msSector()->idHelperSvc()->toString(hit.spacePoint()->identify())<<endmsg;
121  }
122  }
123  return chi2;
124  }

◆ chiSqTerm() [2/2]

double MuonR4::SegmentFitHelpers::chiSqTerm ( const Amg::Vector3D &  posInChamber,
const Amg::Vector3D &  dirInChamber,
const SpacePoint &  measurement,
MsgStream &  msg 
)

Calculates the chi2 contribuation to a linear segment line from an uncalibrated measurement.

Currently only Mdt, Tgc & Rpc are supported

Parameters
posInChamberPosition of the chamber crossing expressed at z=0
dirInChamberSegment direction inside the chamber
measurementMdt measurement
msgReference to the callers msgStream. If the level is VERBOSE, then all relevant parameters are printed

Definition at line 28 of file SegmentFitHelperFunctions.cxx.

31  {
32  double chi2{0.};
33  switch (measurement.type()) {
34  case xAOD::UncalibMeasType::MdtDriftCircleType:
35  chi2 = chiSqTermMdt(posInChamber, dirInChamber,
36  measurement, msg);
37  break;
38  case xAOD::UncalibMeasType::RpcStripType:
39  case xAOD::UncalibMeasType::TgcStripType:
40  case xAOD::UncalibMeasType::MMClusterType:
41  case xAOD::UncalibMeasType::sTgcStripType:
42  chi2 = chiSqTermStrip(posInChamber, dirInChamber, measurement, msg);
43  break;
44  default:
45  if (msg.level() <= MSG::WARNING) {
46  msg<<MSG::WARNING<<__FILE__<<":"<<__LINE__<<" - Unsupported measurement: "
47  <<measurement.msSector()->idHelperSvc()->toString(measurement.identify())<<endmsg;
48  }
49  }
50  return chi2;
51  }

◆ chiSqTermBeamspot()

double MuonR4::SegmentFitHelpers::chiSqTermBeamspot ( const Amg::Vector3D &  posInChamber,
const Amg::Vector3D &  dirInChamber,
const CalibratedSpacePoint &  beamSpotMeas,
MsgStream &  msg 
)

Calculates the chi2 contribution from an external beam spot constraint.

Parameters
posInChamberPosition of the segment in the sector frame
dirInChamberDirection of flight of the segment expressed in the sector frame
beamSpotMeasStrip measurement to consider.
msgReference to the callers msgStream. If the level is VERBOSE, then all relevant parameters are printed

Invalid spacepoints do not contribute to the chi2

Invalid spacepoints do not contribute to the chi2

Definition at line 191 of file SegmentFitHelperFunctions.cxx.

194  {
195 
197  if (beamSpotMeas.fitState() != State::Valid) return 0.;
198  const Amg::Vector3D planeCrossing{segPos + Amg::intersect<3>(segPos, segDir, Amg::Vector3D::UnitZ(),
199  beamSpotMeas.positionInChamber().z()).value_or(0) * segDir};
200 
201  const Amg::Vector2D residual = (planeCrossing - beamSpotMeas.positionInChamber()).block<2,1>(0,0);
202  const double chi2 = contract(inverse(beamSpotMeas.covariance()), residual);
203  if (msg.level() <= printLvl) {
204  msg<<printLvl<<"Print beamspot constraint "<<Amg::toString(beamSpotMeas.positionInChamber())<<endmsg;
205  msg<<printLvl<<"Segment: "<<Amg::toString(segPos) <<" + x "<<Amg::toString(segDir)<<endmsg;
206  msg<<printLvl<<"Plane crossing: "<<Amg::toString(planeCrossing)<<", residual: "<<Amg::toString(residual)
207  <<", covariance: "<<std::endl<<toString(beamSpotMeas.covariance())<<std::endl
208  <<", inverse: "<<std::endl<<toString(inverse(beamSpotMeas.covariance()))<<std::endl
209  <<"chi2 term: "<<chi2<<endmsg;
210  }
211  return chi2;
212  }

◆ chiSqTermMdt() [1/2]

double MuonR4::SegmentFitHelpers::chiSqTermMdt ( const Amg::Vector3D &  posInChamber,
const Amg::Vector3D &  dirInChamber,
const CalibratedSpacePoint &  mdtSpacePoint,
MsgStream &  msg 
)

Calculates the chi2 contribution from a mdt space point to the segment line.

Parameters
posInChamberPosition of the segment in the sector frame
dirInChamberDirection of flight of the segment expressed in the sector frame
mdtSpacePointSpace point to which the chi2 term is calculated
msgReference to the callers msgStream. If the level is VERBOSE, then all relevant parameters are printed

Invalid spacepoints do not contribute to the chi2

Invalid spacepoints do not contribute to the chi2

Definition at line 126 of file SegmentFitHelperFunctions.cxx.

129  {
131  if (hit.fitState() != State::Valid) return 0.;
132  const Amg::Vector3D& hitPos{hit.positionInChamber()};
133  const Amg::Vector3D& hitDir{hit.directionInChamber()};
134 
135  // Calculate the closest point along the segment to the wire
136  const Amg::Vector3D closePointSeg = segPos + Amg::intersect<3>(hitPos, hitDir, segPos, segDir).value_or(0)* segDir;
137  // Closest point along the wire to the segment
138  const Amg::Vector3D closePointWire = hitPos + hitDir.dot(closePointSeg - hitPos) * hitDir;
139  Amg::Vector2D residual{Amg::Vector2D::Zero()};
140  residual[Amg::x] = (hitPos - closePointSeg).x();
141  residual[Amg::y] = (closePointWire - closePointSeg).mag() - hit.driftRadius();
142  const double chi2{contract(inverse(hit.covariance()), residual)};
143  if (msg.level() <= printLvl) {
144  const SpacePoint* sp = hit.spacePoint();
145  msg<<printLvl<<"Calibrated measurement "<<sp->msSector()->idHelperSvc()->toString(sp->identify());
146  msg<<printLvl<<", position: "<<Amg::toString(hit.positionInChamber())<<", driftRadius: "<<hit.driftRadius()
147  <<", dimension: "<<sp->dimension()<<endmsg;
148  msg<<printLvl<<", segment: "<<Amg::toString(segPos) <<" + x "<<Amg::toString(segDir)<<endmsg;
149  msg<<printLvl<<", residual: "<<Amg::toString(residual)<<"/ "<<
150  (Amg::lineDistance<3>(hitPos, hitDir, segPos, segDir) -hit.driftRadius())<<std::endl <<", covariance: "<<std::endl
151  <<toString(hit.covariance())<<", inverse: "
152  <<std::endl<<toString(inverse(hit.covariance()))<<", "<<endmsg;
153  msg<<printLvl<<"Chi2 term: "<<chi2<<endmsg;
154  }
155  assert(chi2 >=0.);
156  return chi2;
157  }

◆ chiSqTermMdt() [2/2]

double MuonR4::SegmentFitHelpers::chiSqTermMdt ( const Amg::Vector3D &  posInChamber,
const Amg::Vector3D &  dirInChamber,
const SpacePoint &  measurement,
MsgStream &  msg 
)

Calculates the chi2 contribuation to a linear segment line from an uncalibrated Mdt measurement.

Parameters
posInChamberPosition of the chamber crossing expressed at z=0
dirInChamberSegment direction inside the chamber
measurementMdt measurement
msgReference to the callers msgStream. If the level is VERBOSE, then all relevant parameters are printed

Definition at line 52 of file SegmentFitHelperFunctions.cxx.

55  {
56 
57  Amg::Vector2D residual{Amg::Vector2D::Zero()};
58 
59  residual[Amg::y] = std::abs(Amg::signedDistance(mdtSP.positionInChamber(), mdtSP.directionInChamber(), segPos, segDir))
60  - mdtSP.driftRadius();
61  const double chi2{residual.dot(mdtSP.covariance().inverse()* residual)};
62  if (msg.level() <= printLvl) {
63  msg<<printLvl<<"Measurement "<<mdtSP.msSector()->idHelperSvc()->toString(mdtSP.identify());
64  msg<<printLvl<<", position: "<<Amg::toString(mdtSP.positionInChamber())
65  <<", driftRadius: "<<mdtSP.driftRadius()
66  <<", status: "<<static_cast<const xAOD::MdtDriftCircle*>(mdtSP.primaryMeasurement())->status()<<endmsg;
67  msg<<printLvl<<"segment: "<<Amg::toString(segPos) <<" + x "<<Amg::toString(segDir)<<","<<endmsg;
68  msg<<printLvl<<"residual: "<<Amg::toString(residual) <<", covariance: "<<std::endl
69  <<Amg::toString(mdtSP.covariance())<<std::endl<<", inverse: "
70  <<std::endl<<Amg::toString(mdtSP.covariance().inverse())<<", "<<endmsg;
71  msg<<printLvl<<"chi2 term: "<<chi2<<endmsg;
72  }
73  assert(chi2 >=0.);
74  return chi2;
75  }

◆ chiSqTermStrip() [1/2]

double MuonR4::SegmentFitHelpers::chiSqTermStrip ( const Amg::Vector3D &  posInChamber,
const Amg::Vector3D &  dirInChamber,
const double  timeShift,
std::optional< double >  arrivalTime,
const CalibratedSpacePoint &  strip,
MsgStream &  msg 
)

Calculates the chi2 contribution from a strip measurement to the segment line.

Parameters
posInChamberPosition of the segment in the sector frame
dirInChamberDirection of flight of the segment expressed in the sector frame
timeShiftShift from the nominal arrival time, that's calculated as R / c.
arrivalTimeArrival time of the particle at the segment's refernce plane. If the arrival time is passed, the timing information of the measurement is explicitly taken into account.
stripStrip measurement to consider.
msgReference to the callers msgStream. If the level is VERBOSE, then all relevant parameters are printed

Invalid spacepoints do not contribute to the chi2

Invalid spacepoints do not contribute to the chi2

Definition at line 158 of file SegmentFitHelperFunctions.cxx.

163  {
164 
166  if (strip.fitState() != State::Valid) return 0.;
167 
168  const SpacePoint* spacePoint = strip.spacePoint();
169  const Amg::Vector3D normal = spacePoint->planeNormal();
170  std::optional<double> travelledDist = Amg::intersect<3>(segPos, segDir, normal, normal.dot(spacePoint->positionInChamber()));
171  const Amg::Vector3D planeCrossing = segPos + travelledDist.value_or(0) * segDir;
172  Amg::Vector3D residual{(planeCrossing - strip.positionInChamber())};
173  if (strip.measuresTime() && arrivalTime) {
174  residual[Amg::z] = strip.time() - (*arrivalTime) - travelledDist.value_or(0) * c_inv - offsetTime;
175  } else {
176  residual[Amg::z] = 0;
177  }
178  const double chi2 = contract(inverse(strip.covariance()), residual);
179  if (msg.level() <= printLvl) {
180  msg<<printLvl<<"Calibrated measurement "<<spacePoint->msSector()->idHelperSvc()->toString(spacePoint->identify())
181  <<", position: "<<Amg::toString(strip.positionInChamber())<<endmsg;
182  msg<<printLvl<<"Segment: "<<Amg::toString(segPos) <<" + x "<<Amg::toString(segDir)<<endmsg;
183  msg<<printLvl<<"Plane crossing: "<<Amg::toString(planeCrossing)<<", residual: "<<Amg::toString(residual)
184  <<", covariance: "<<std::endl<<toString(strip.covariance())<<", inverse: "
185  <<std::endl<<toString(inverse(strip.covariance()))<<endmsg;
186  msg<<printLvl<<" -- final chi2 term: "<<chi2<<endmsg;
187  }
188  assert(chi2 >=0.);
189  return chi2;
190  }

◆ chiSqTermStrip() [2/2]

double MuonR4::SegmentFitHelpers::chiSqTermStrip ( const Amg::Vector3D &  posInChamber,
const Amg::Vector3D &  dirInChamber,
const SpacePoint &  measurement,
MsgStream &  msg 
)

Calculates the chi2 contribuation to a linear segment line from an uncalibrated strip measurement.

Currently only Tgc & Rpc are supported

Parameters
posInChamberPosition of the chamber crossing expressed at z=0
dirInChamberSegment direction inside the chamber
measurementMdt measurement
msgReference to the callers msgStream. If the level is VERBOSE, then all relevant parameters are printed

Definition at line 76 of file SegmentFitHelperFunctions.cxx.

79  {
80  const Amg::Vector3D normal = stripSP.planeNormal();
81  std::optional<double> travelledDist = Amg::intersect<3>(segPos, segDir, normal, normal.dot(stripSP.positionInChamber()));
82  const Amg::Vector3D planeCrossing = segPos + travelledDist.value_or(0) * segDir;
83 
84  const Amg::Vector2D residual{(planeCrossing - stripSP.positionInChamber()).block<2,1>(0,0)};
85  const double chi2 = residual.dot(stripSP.covariance().inverse()* residual);
86  if (msg.level() <= printLvl) {
87  msg<<printLvl<<"Measurement "<<stripSP.msSector()->idHelperSvc()->toString(stripSP.identify())
88  <<", position: "<<Amg::toString(stripSP.positionInChamber())<<endmsg;
89  msg<<printLvl<<"Segment: "<<Amg::toString(segPos) <<" + x "<<Amg::toString(segDir)<<endmsg;
90  msg<<printLvl<<"Plane crossing: "<<Amg::toString(planeCrossing)<<", residual: "<<Amg::toString(residual)
91  <<", covariance: " <<std::endl<<Amg::toString(stripSP.covariance())<<", inverse: "
92  <<std::endl<<Amg::toString(stripSP.covariance().inverse())<<","<<endmsg;
93  msg<<printLvl<<"chi2 term: "<<chi2<<endmsg;
94  }
95  assert(chi2 >=0.);
96  return chi2;
97  }

◆ driftSign() [1/2]

int MuonR4::SegmentFitHelpers::driftSign ( const Amg::Vector3D &  posInChamber,
const Amg::Vector3D &  dirInChamber,
const CalibratedSpacePoint &  calibHit,
MsgStream &  msg 
)

Calculates whether a segement line travereses the tube measurement on the left (-1) or right (1) side of the tube wire.

Strip measurements & nullptrs are assigned with 0

Parameters
posInChamberPosition of the segment in the sector frame
dirInChamberDirection of flight of the segment expressed in the sector frame
calibHitCalibrated hit to consider
msgReference to the callers msgStream. If the level is VERBOSE, then all relevant parameters are printed

Definition at line 248 of file SegmentFitHelperFunctions.cxx.

249  {
250  if (calibHit.type() != xAOD::UncalibMeasType::MdtDriftCircleType){
251  return 0;
252  }
253  const double signedDist = Amg::signedDistance(segPos, segDir, calibHit.positionInChamber(), calibHit.directionInChamber());
254  if (msg.level() <= printLvl) {
255  const SpacePoint* sp = calibHit.spacePoint();
256  msg<<printLvl<<"Hit "<<sp->msSector()->idHelperSvc()->toString(sp->identify())<<" drift radius "<<calibHit.driftRadius()
257  <<", signed distance: "<<signedDist<<", unsigned distance: "
258  <<Amg::lineDistance<3>(segPos, segDir, calibHit.positionInChamber(), calibHit.directionInChamber())<<endmsg;
259  }
260  return sign(signedDist);
261  }

◆ driftSign() [2/2]

int MuonR4::SegmentFitHelpers::driftSign ( const Amg::Vector3D &  posInChamber,
const Amg::Vector3D &  dirInChamber,
const SpacePoint &  uncalibHit,
MsgStream &  msg 
)

Calculates whether a segement line travereses the tube measurement on the left (-1) or right (1) side of the tube wire.

Strip measurements & nullptrs are assigned with 0

Parameters
posInChamberPosition of the segment in the sector frame
dirInChamberDirection of flight of the segment expressed in the sector frame
uncalibHitUncalibrated hit to consider
msgReference to the callers msgStream. If the level is VERBOSE, then all relevant parameters are printed

Definition at line 235 of file SegmentFitHelperFunctions.cxx.

236  {
237  if (sp.type() != xAOD::UncalibMeasType::MdtDriftCircleType) {
238  return 0;
239  }
240  const double signedDist = Amg::signedDistance(segPos, segDir, sp.positionInChamber(), sp.directionInChamber());
241  if (msg.level() <= printLvl) {
242  msg<<printLvl<<"Hit "<<sp.msSector()->idHelperSvc()->toString(sp.identify())<<" drift radius "<<sp.driftRadius()
243  <<", signed distance: "<<signedDist<<", unsigned distance: "
244  <<Amg::lineDistance<3>(segPos, segDir, sp.positionInChamber(), sp.directionInChamber())<<endmsg;
245  }
246  return sign(signedDist);
247  }

◆ driftSigns() [1/2]

std::vector< int > MuonR4::SegmentFitHelpers::driftSigns ( const Amg::Vector3D &  posInChamber,
const Amg::Vector3D &  dirInChamber,
const std::vector< const SpacePoint * > &  uncalibHits,
MsgStream &  msg 
)

Calculates whether a segment line travereses the tube measurements on the left (-1) or right (1) side of the tube wire.

Strip measurements & nullptrs are assigned with 0

Parameters
posInChamberPosition of the segment in the sector frame
dirInChamberDirection of flight of the segment expressed in the sector frame
uncalibHitsList of uncalibrated measurements
msgReference to the callers msgStream. If the level is VERBOSE, then all relevant parameters are printed

Definition at line 213 of file SegmentFitHelperFunctions.cxx.

216  {
217 
218  std::vector<int> signs{};
219  signs.reserve(uncalibHits.size());
220  for (const SpacePoint* sp : uncalibHits) {
221  signs.push_back( sp ? driftSign(segPos,segDir, *sp, msg) : 0);
222  }
223  return signs;
224  }

◆ driftSigns() [2/2]

std::vector< int > MuonR4::SegmentFitHelpers::driftSigns ( const Amg::Vector3D &  posInChamber,
const Amg::Vector3D &  dirInChamber,
const std::vector< std::unique_ptr< CalibratedSpacePoint >> &  calibHits,
MsgStream &  msg 
)

Calculates whether a segment line travereses the tube measurements on the left (-1) or right (1) side of the tube wire.

Strip measurements & nullptrs are assigned with 0

Parameters
posInChamberPosition of the segment in the sector frame
dirInChamberDirection of flight of the segment expressed in the sector frame
calibHitsList of calibrated measurements
msgReference to the callers msgStream. If the level is VERBOSE, then all relevant parameters are printed

Definition at line 225 of file SegmentFitHelperFunctions.cxx.

227  {
228  std::vector<int> signs{};
229  signs.reserve(calibHits.size());
230  for (const std::unique_ptr<CalibratedSpacePoint>& hit : calibHits) {
231  signs.push_back(driftSign(segPos,segDir, *hit, msg));
232  }
233  return signs;
234  }
ClusterSeg::residual
@ residual
Definition: ClusterNtuple.h:20
Amg::Vector2D
Eigen::Matrix< double, 2, 1 > Vector2D
Definition: GeoPrimitives.h:48
Amg::y
@ y
Definition: GeoPrimitives.h:35
xAOD::UncalibMeasType::MMClusterType
@ MMClusterType
xAOD::UncalibMeasType::sTgcStripType
@ sTgcStripType
x
#define x
keylayer_zslicemap.strip
strip
Definition: keylayer_zslicemap.py:151
SpacePoint
Definition: Trigger/TrigAccel/TrigCudaFitter/src/SpacePoint.h:7
xAOD::UncalibMeasType::TgcStripType
@ TgcStripType
MuonR4::SegmentFitHelpers::driftSign
int driftSign(const Amg::Vector3D &posInChamber, const Amg::Vector3D &dirInChamber, const SpacePoint &uncalibHit, MsgStream &msg)
Calculates whether a segement line travereses the tube measurement on the left (-1) or right (1) side...
Definition: SegmentFitHelperFunctions.cxx:235
Amg::z
@ z
Definition: GeoPrimitives.h:36
MuonR4::SegmentFitHelpers::chiSqTermStrip
double chiSqTermStrip(const Amg::Vector3D &posInChamber, const Amg::Vector3D &dirInChamber, const SpacePoint &measurement, MsgStream &msg)
Calculates the chi2 contribuation to a linear segment line from an uncalibrated strip measurement.
Definition: SegmentFitHelperFunctions.cxx:76
Amg::toString
std::string toString(const Translation3D &translation, int precision=4)
GeoPrimitvesToStringConverter.
Definition: GeoPrimitivesToStringConverter.h:40
MuonR4::SegmentFitHelpers::chiSqTermMdt
double chiSqTermMdt(const Amg::Vector3D &posInChamber, const Amg::Vector3D &dirInChamber, const SpacePoint &measurement, MsgStream &msg)
Calculates the chi2 contribuation to a linear segment line from an uncalibrated Mdt measurement.
Definition: SegmentFitHelperFunctions.cxx:52
endmsg
#define endmsg
Definition: AnalysisConfig_Ntuple.cxx:63
MuonR4::inverse
CalibratedSpacePoint::Covariance_t inverse(const CalibratedSpacePoint::Covariance_t &mat)
Inverts the parsed matrix.
Definition: MuonSpectrometer/MuonPhaseII/Event/MuonSpacePoint/src/UtilFunctions.cxx:65
Amg::x
@ x
Definition: GeoPrimitives.h:34
chi2
double chi2(TH1 *h0, TH1 *h1)
Definition: comparitor.cxx:525
sign
int sign(int a)
Definition: TRT_StrawNeighbourSvc.h:107
xAOD::Other
@ Other
Amg::Vector3D
Eigen::Matrix< double, 3, 1 > Vector3D
Definition: GeoPrimitives.h:47
Amg::signedDistance
double signedDistance(const Amg::Vector3D &posA, const Amg::Vector3D &dirA, const Amg::Vector3D &posB, const Amg::Vector3D &dirB)
Calculates the signed distance between two lines in 3D space.
Definition: GeoPrimitivesHelpers.h:329
merge.status
status
Definition: merge.py:17
xAOD::MdtDriftCircle_v1
https://gitlab.cern.ch/atlas/athena/-/blob/master/MuonSpectrometer/MuonReconstruction/MuonRecEvent/Mu...
Definition: MdtDriftCircle_v1.h:21
xAOD::UncalibMeasType::RpcStripType
@ RpcStripType
python.AutoConfigFlags.msg
msg
Definition: AutoConfigFlags.py:7
xAOD::UncalibMeasType::MdtDriftCircleType
@ MdtDriftCircleType
mag
Scalar mag() const
mag method
Definition: AmgMatrixBasePlugin.h:26
generate::Zero
void Zero(TH1D *hin)
Definition: generate.cxx:32
MuonR4::SegmentFitHelpers::chiSqTermBeamspot
double chiSqTermBeamspot(const Amg::Vector3D &posInChamber, const Amg::Vector3D &dirInChamber, const CalibratedSpacePoint &beamSpotMeas, MsgStream &msg)
Calculates the chi2 contribution from an external beam spot constraint.
Definition: SegmentFitHelperFunctions.cxx:191
MuonR4::contract
double contract(const CalibratedSpacePoint::Covariance_t &mat, const Amg::Vector3D &a, const Amg::Vector3D &b)
Definition: MuonSpectrometer/MuonPhaseII/Event/MuonSpacePoint/src/UtilFunctions.cxx:13
Generated on Sat Apr 19 2025 22:04:49 for ATLAS Offline Software by   doxygen 1.8.18