Generate and assign the contained variables.
21{
22 if (not ptrSegmentPair) return false;
23
24 const CDCSegmentPair& segmentPair = *ptrSegmentPair;
25
27 const CDCSegment2D* ptrToSegment = segmentPair.
getToSegment();
28
29 const CDCSegment2D& fromSegment = *ptrFromSegment;
30 const CDCSegment2D& toSegment = *ptrToSegment;
31
33 const CDCSegment2D& axialSegment = *ptrAxialSegment;
34
36 const CDCSegment2D& stereoSegment = *ptrStereoSegment;
37
38
42
43
46
49
52
53
54 using namespace NPerigeeParameterIndices;
57
60
61 const double deltaCurvVar = fromCurvVar + toCurvVar;
62 const double avgPrecision = 1 / fromCurvVar + 1 / toCurvVar;
63 const double deltaCurvSigma = std::sqrt(deltaCurvVar);
64
65 finitevar<
named(
"abs_avg_curv")>() = std::fabs(toCurv / toCurvVar + fromCurv / fromCurvVar) / avgPrecision;
68 finitevar<
named(
"delta_curv_pull")>() = (toCurv - fromCurv) / deltaCurvSigma;
69
70
71 const CDCRecoHit2D& fromFirstHit = fromSegment.front();
72 const CDCRecoHit2D& fromLastHit = fromSegment.back();
73 const CDCRecoHit2D& toFirstHit = toSegment.front();
74 const CDCRecoHit2D& toLastHit = toSegment.back();
75
78
79
80 const Vector2D fromFitPos = fromFit.
getClosest(fromHitPos);
81 const Vector2D toFitPos = toFit.
getClosest(toHitPos);
84
87
88 const double deltaPosPhi = fromFitPos.
angleWith(toFitPos);
89 const double deltaMomPhi = fromFitMom.
angleWith(toFitMom);
91
94
98
99
100
101 const CDCRecoHit2D& nearAxialHit = toFirstHit.
isAxial() ? toFirstHit : fromLastHit;
102 const CDCRecoHit2D& farStereoHit = not fromFirstHit.
isAxial() ? fromFirstHit : toLastHit;
103 const CDCRecoHit2D& nearStereoHit = not toFirstHit.
isAxial() ? toFirstHit : fromLastHit;
104
105 const CDCWire& farStereoWire = farStereoHit.
getWire();
106 const WireLine& farWireLine = farStereoWire.
getWireLine();
107
108 const CDCWire& nearStereoWire = nearStereoHit.
getWire();
109 const WireLine& nearWireLine = nearStereoWire.
getWireLine();
110
111 const Vector3D nearAxialRecoPos = nearAxialHit.reconstruct3D(axialFit);
113 const double farZ = farStereoRecoPos.
z();
114
116 const double nearZ = nearStereoRecoPos.
z();
117
118 const double stereoArcLength2D =
120 farStereoRecoPos.
xy());
121
122 const double arcLength2DGap =
124 nearStereoRecoPos.
xy());
125
127 finitevar<
named(
"stereo_arc_length")>() = fabs(stereoArcLength2D);
128
131
134
135 finitevar<
named(
"coarse_tanl")>() = (farZ - nearZ) / stereoArcLength2D;
136
137 finitevar<
named(
"stereo_rel_size")>() = fabs(stereoSegment.size() / stereoArcLength2D);
138
142
146
149
153
154 return true;
155}
bool isAxial() const
Indicator if the underlying wire is axial.
Vector3D reconstruct3D(const CDCTrajectory2D &trajectory2D, const double z=0) const
Reconstruct the three dimensional position (especially of stereo hits) by determining the z coordinat...
const CDCWire & getWire() const
Getter for the wire the reconstructed hit associated to.
Vector2D getRecoPos2D() const
Getter for the position in the reference plane.
const CDCSegment2D * getAxialSegment() const
Getter for the axial segment.
const CDCSegment2D * getToSegment() const
Getter for the to segment.
const CDCSegment2D * getStereoSegment() const
Getter for the stereo segment.
const CDCSegment2D * getFromSegment() const
Getter for the from segment.
CDCTrajectory2D & getTrajectory2D() const
Getter for the two dimensional trajectory fitted to the segment.
double getTotalArcLength2D(const AHits &hits) const
Calculates the perpendicular travel distance from the first position of the hits to the last position...
double getChi2() const
Getter for the chi2 value of the circle fit.
double calcArcLength2DBetween(const Vector2D &fromPoint, const Vector2D &toPoint) const
Calculate the travel distance between the two given positions Returns the travel distance on the traj...
double getPValue() const
Getter for p-value.
bool isFitted() const
Checks if the circle is already set to a valid value.
double getLocalVariance(EPerigeeParameter i) const
Getter for an individual diagonal element of the covariance matrix of the local helix parameters.
double getArcLength2DGap(const AFromHits &fromHits, const AToHits &toHits) const
Calculates the perpendicular travel distance from the last position of the fromHits to the first posi...
Vector2D getClosest(const Vector2D &point) const
Calculates the closest approach on the trajectory to the given point.
double getArcLength2DBackOffset(const AFromHits &fromHits, const AToHits &toHits) const
Calculates the perpendicular travel distance from the last position of the fromHits to the last posit...
Vector2D getFlightDirection2D(const Vector2D &point) const
Get the unit direction of flight at the given point, where arcLength2D = 0.
double getArcLength2DFrontOffset(const AFromHits &fromHits, const AToHits &toHits) const
Calculates the perpendicular travel distance from the first position of the fromHits to the first pos...
size_t getNDF() const
Getter for the number of degrees of freedom of the circle fit.
double getCurvature() const
Getter for the curvature as seen from the xy projection.
const WireLine & getWireLine() const
Getter for the wire line representation of the wire.
AssignFinite< Float_t > finitevar()
static constexpr int named(const char *name)
double angleWith(const Vector2D &rhs) const
The angle between this and rhs.
const Vector2D & xy() const
Getter for the xy projected vector ( reference ! )
double z() const
Getter for the z coordinate.
double outOfZBoundsFactor(double z) const
Returns the amount how much the given z position is outside the bounds in units of the wire length.
HepGeom::Vector3D< double > Vector3D
3D Vector
static double normalised(const double angle)
Normalise an angle to lie in the range from [-pi, pi].