Belle II Software  release-05-01-25
HitGapTrackRelationVarSet.cc
1 /**************************************************************************
2  * BASF2 (Belle Analysis Framework 2) *
3  * Copyright(C) 2016 - Belle II Collaboration *
4  * *
5  * Author: The Belle II Collaboration *
6  * Contributors: Oliver Frost *
7  * *
8  * This software is provided "as is" without any warranty. *
9  **************************************************************************/
10 #include <tracking/trackFindingCDC/filters/trackRelation/HitGapTrackRelationVarSet.h>
11 
12 #include <tracking/trackFindingCDC/eventdata/tracks/CDCTrack.h>
13 
14 #include <tracking/trackFindingCDC/numerics/Angle.h>
15 
16 using namespace Belle2;
17 using namespace TrackFindingCDC;
18 
20 {
21  if (not ptrTrackRelation) return false;
22 
23  const Relation<const CDCTrack>& trackPair = *ptrTrackRelation;
24 
25  const CDCTrack* ptrFromTrack = trackPair.getFrom();
26  const CDCTrack* ptrToTrack = trackPair.getTo();
27 
28  const CDCTrack& fromTrack = *ptrFromTrack;
29  const CDCTrack& toTrack = *ptrToTrack;
30 
31  const CDCRecoHit3D& fromFirstHit = fromTrack.front();
32  const CDCRecoHit3D& fromLastHit = fromTrack.back();
33 
34  const CDCRecoHit3D& toFirstHit = toTrack.front();
35  const CDCRecoHit3D& toLastHit = toTrack.back();
36 
37  const Vector3D fromLastHitPos = fromLastHit.getRecoPos3D();
38  const Vector3D fromFirstHitPos = fromFirstHit.getRecoPos3D();
39 
40  const Vector3D toFirstHitPos = toFirstHit.getRecoPos3D();
41  const Vector3D toLastHitPos = toLastHit.getRecoPos3D();
42 
43  const Vector3D hitPosGap = toFirstHitPos - fromLastHitPos;
44  const Vector3D longHitPosGap = toLastHitPos - fromFirstHitPos;
45 
46  const double hitDistance = hitPosGap.norm();
47  const double longHitDistance = longHitPosGap.norm();
48 
49  const Vector2D fromLastHitMom = fromLastHit.getFlightDirection2D();
50  const Vector2D toFirstHitMom = toFirstHit.getFlightDirection2D();
51 
52  finitevar<named("delta_hit_pos_phi")>() = fromLastHitPos.angleWith(toFirstHitPos);
53  finitevar<named("delta_hit_mom_phi")>() = fromLastHitMom.angleWith(toFirstHitMom);
54 
55  double fromLastHitAlpha = fromLastHit.getAlpha();
56  double toFirstHitAlpha = toFirstHit.getAlpha();
57  finitevar<named("delta_hit_alpha")>() = AngleUtil::normalised(toFirstHitAlpha - fromLastHitAlpha);
58 
59  finitevar<named("delta_hit_z")>() = fromLastHitPos.z() - toFirstHitPos.z();
60 
61  finitevar<named("hit_distance")>() = hitDistance;
62  finitevar<named("hit_long_distance")>() = longHitDistance;
63 
64  finitevar<named("delta_hit_distance")>() = longHitDistance - hitDistance;
65 
66  finitevar<named("from_hit_forward")>() = hitPosGap.xy().dot(fromLastHitMom);
67  finitevar<named("to_hit_forward")>() = hitPosGap.xy().dot(toFirstHitMom);
68  finitevar<named("hit_forward")>() = hitPosGap.xy().dot(Vector2D::average(fromLastHitMom, toFirstHitMom));
69 
70  const Vector3D fromStretch = fromLastHitPos - fromFirstHitPos;
71  const Vector3D toStretch = toLastHitPos - toFirstHitPos;
72 
73  const double fromLength = fromStretch.norm();
74  const double toLength = toStretch.norm();
75 
76  const Vector3D firstPosGap = toFirstHitPos - fromFirstHitPos;
77  const Vector3D lastPosGap = toLastHitPos - fromLastHitPos;
78 
79  const double firstOffset = firstPosGap.norm();
80  const double lastOffset = lastPosGap.norm();
81 
82  finitevar<named("hit_ptolemy")>() =
83  firstOffset * lastOffset - longHitDistance * hitDistance - fromLength * toLength;
84 
85  return true;
86 }
Belle2::TrackFindingCDC::Relation::getFrom
From * getFrom() const
Getter for the pointer to the from side object.
Definition: Relation.h:69
Belle2::TrackFindingCDC::CDCRecoHit3D
Class representing a three dimensional reconstructed hit.
Definition: CDCRecoHit3D.h:62
Belle2::TrackFindingCDC::CDCTrack
Class representing a sequence of three dimensional reconstructed hits.
Definition: CDCTrack.h:51
Belle2::TrackFindingCDC::Relation
Type for two related objects.
Definition: CDCSegment2D.h:37
Belle2::TrackFindingCDC::Vector2D
A two dimensional vector which is equipped with functions for correct handeling of orientation relat...
Definition: Vector2D.h:37
Belle2::TrackFindingCDC::VarSet< HitGapTrackRelationVarNames >::finitevar
AssignFinite< Float_t > finitevar()
Reference getter for the value of the ith variable. Transforms non-finite values to finite value.
Definition: VarSet.h:140
Belle2::TrackFindingCDC::Vector2D::dot
double dot(const Vector2D &rhs) const
Calculates the two dimensional dot product.
Definition: Vector2D.h:172
Belle2::TrackFindingCDC::Vector3D::norm
double norm() const
Calculates the length of the vector.
Definition: Vector3D.h:204
Belle2::TrackFindingCDC::CDCRecoHit3D::getAlpha
double getAlpha() const
Getter for the direction of flight relative to the position.
Definition: CDCRecoHit3D.h:329
Belle2::TrackFindingCDC::CDCRecoHit3D::getFlightDirection2D
Vector2D getFlightDirection2D() const
Getter for the direction of flight.
Definition: CDCRecoHit3D.h:322
Belle2
Abstract base class for different kinds of events.
Definition: MillepedeAlgorithm.h:19
Belle2::TrackFindingCDC::Vector3D
A three dimensional vector.
Definition: Vector3D.h:34
Belle2::TrackFindingCDC::HitGapTrackRelationVarSet::extract
bool extract(const Relation< const CDCTrack > *ptrTrackRelation) final
Generate and assign the contained variables.
Definition: HitGapTrackRelationVarSet.cc:19
Belle2::TrackFindingCDC::Vector3D::xy
const Vector2D & xy() const
Getter for the xy projected vector ( reference ! )
Definition: Vector3D.h:500
Belle2::TrackFindingCDC::VarSet< HitGapTrackRelationVarNames >::named
constexpr static int named(const char *name)
Getter for the index from the name.
Definition: VarSet.h:88
Belle2::TrackFindingCDC::Relation::getTo
To * getTo() const
Getter for the pointer to the to side object.
Definition: Relation.h:75
Belle2::TrackFindingCDC::Vector2D::angleWith
double angleWith(const Vector2D &rhs) const
The angle between this and rhs.
Definition: Vector2D.h:211
Belle2::TrackFindingCDC::Vector3D::z
double z() const
Getter for the z coordinate.
Definition: Vector3D.h:488
Belle2::TrackFindingCDC::CDCRecoHit3D::getRecoPos3D
const Vector3D & getRecoPos3D() const
Getter for the 3d position of the hit.
Definition: CDCRecoHit3D.h:295
Belle2::TrackFindingCDC::AngleUtil::normalised
static double normalised(const double angle)
Normalise an angle to lie in the range from [-pi, pi].
Definition: Angle.h:43
Belle2::TrackFindingCDC::Vector2D::average
static Vector2D average(const Vector2D &one, const Vector2D &two)
Constructs the average of two vectors.
Definition: Vector2D.h:95