38 static double value(
const PointType&
outerHit,
const PointType& centerHit,
const PointType& innerHit)
42 double circleRadius = CircleRadius<PointType>::calcAvgDistanceXY(
outerHit, centerHit, innerHit, cCenter);
46 B2Vector3D vecCenter2cC((centerHit.X() - cCenter.
X()),
47 (centerHit.Y() - cCenter.
Y()),
48 (centerHit.Z() - cCenter.
Z()));
49 B2Vector3D vecInner2cC((innerHit.X() - cCenter.
X()),
50 (innerHit.Y() - cCenter.
Y()),
51 (innerHit.Z() - cCenter.
Z()));
55 double alfaOCr = acos(CosDirectionXY<B2Vector3D>::value(vecOuter2cC, vecCenter2cC)) * circleRadius ;
56 double alfaCIr = acos(CosDirectionXY<B2Vector3D>::value(vecCenter2cC, vecInner2cC)) * circleRadius ;
59 double result = (asin(
double(
outerHit.Z() - centerHit.Z()) / alfaOCr)) - asin(
double(centerHit.Z() - innerHit.Z()) / alfaCIr);
61 return (std::isnan(result) || std::isinf(result)) ? double(0) : result;