13 #include <tracking/trackFindingVXD/filterMap/filterFramework/SelectionVariable.h>
14 #include <framework/geometry/B2Vector3.h>
18 #define ANGLE3DFULL_NAME Angle3DFull
31 template <
typename Po
intType >
40 static double value(
const PointType& outerHit,
const PointType& centerHit,
const PointType& innerHit)
43 B2Vector3D outerVector(outerHit.X() - centerHit.X(), outerHit.Y() - centerHit.Y(), outerHit.Z() - centerHit.Z());
44 B2Vector3D innerVector(centerHit.X() - innerHit.X(), centerHit.Y() - innerHit.Y(), centerHit.Z() - innerHit.Z());
46 double result = acos(outerVector.
Dot(innerVector) / (outerVector.
Mag() * innerVector.
Mag()));
47 result = (result * (180. / M_PI));
48 return (std::isnan(result) || std::isinf(result)) ?
double(0) : result;