13 #include <tracking/trackFindingVXD/filterMap/filterFramework/SelectionVariable.h>
14 #include <tracking/trackFindingVXD/filterMap/threeHitVariables/CircleCenterXY.h>
15 #include <tracking/trackFindingVXD/filterMap/twoHitVariables/CosDirectionXY.h>
16 #include <framework/geometry/B2Vector3.h>
19 #define DELTASOVERZ_NAME DeltaSoverZ
29 template <
typename Po
intType >
38 static double value(
const PointType& outerHit,
const PointType& centerHit,
const PointType& innerHit)
41 B2Vector3D circleCenter = CircleCenterXY<PointType>::value(outerHit, centerHit, innerHit);
42 B2Vector3D points2outerHit((outerHit.X() - circleCenter.
X()),
43 (outerHit.Y() - circleCenter.
Y()),
44 (outerHit.Z() - circleCenter.
Z()));
45 B2Vector3D points2centerHit((centerHit.X() - circleCenter.
X()),
46 (centerHit.Y() - circleCenter.
Y()),
47 (centerHit.Z() - circleCenter.
Z()));
48 B2Vector3D points2innerHit((innerHit.X() - circleCenter.
X()),
49 (innerHit.Y() - circleCenter.
Y()),
50 (innerHit.Z() - circleCenter.
Z()));
52 double alfaOC = acos(CosDirectionXY<B2Vector3D>::value(points2outerHit, points2centerHit));
53 double alfaCI = acos(CosDirectionXY<B2Vector3D>::value(points2centerHit, points2innerHit));
56 return (alfaOC *
double(centerHit.Z() - innerHit.Z())) - (alfaCI * double(outerHit.Z() - centerHit.Z()));