13 #include <tracking/trackFindingVXD/filterMap/filterFramework/SelectionVariable.h>
14 #include <framework/geometry/B2Vector3.h>
18 #define ANGLE3DSIMPLE_NAME Angle3DSimple
30 template <
typename Po
intType >
37 static double value(
const PointType& outerHit,
const PointType& centerHit,
const PointType& innerHit)
40 B2Vector3<double> outerVector(outerHit.X() - centerHit.X(), outerHit.Y() - centerHit.Y(), outerHit.Z() - centerHit.Z());
41 B2Vector3<double> innerVector(centerHit.X() - innerHit.X(), centerHit.Y() - innerHit.Y(), centerHit.Z() - innerHit.Z());
44 double result = outerVector.
Dot(innerVector) / (outerVector.
Mag2() * innerVector.
Mag2());
45 return (std::isnan(result) || std::isinf(result)) ?
double(0) : result;