10 #include <tracking/trackFindingCDC/filters/axialSegmentPair/HitGapAxialSegmentPairVarSet.h>
12 #include <tracking/trackFindingCDC/eventdata/tracks/CDCAxialSegmentPair.h>
13 #include <tracking/trackFindingCDC/eventdata/segments/CDCSegment2D.h>
14 #include <tracking/trackFindingCDC/eventdata/hits/CDCRecoHit2D.h>
16 #include <tracking/trackFindingCDC/numerics/Angle.h>
19 using namespace TrackFindingCDC;
23 if (not ptrAxialSegmentPair)
return false;
45 const Vector2D hitPosGap = toFirstHitPos - fromLastHitPos;
46 const Vector2D longHitPosGap = toLastHitPos - fromFirstHitPos;
48 const double hitDistance = hitPosGap.
norm();
49 const double longHitDistance = longHitPosGap.
norm();
54 finitevar<
named(
"delta_hit_pos_phi")>() = fromLastHitPos.angleWith(toFirstHitPos);
57 double fromLastHitAlpha = fromLastHit.
getAlpha();
58 double toFirstHitAlpha = toFirstHit.
getAlpha();
64 finitevar<
named(
"delta_hit_distance")>() = longHitDistance - hitDistance;
70 const Vector2D fromStretch = fromLastHitPos - fromFirstHitPos;
71 const Vector2D toStretch = toLastHitPos - toFirstHitPos;
73 const double fromLength = fromStretch.
norm();
74 const double toLength = toStretch.
norm();
76 const Vector2D firstPosGap = toFirstHitPos - fromFirstHitPos;
77 const Vector2D lastPosGap = toLastHitPos - fromLastHitPos;
79 const double firstOffset = firstPosGap.
norm();
80 const double lastOffset = lastPosGap.
norm();
83 firstOffset * lastOffset - longHitDistance * hitDistance - fromLength * toLength;