add field getter for doubles

Author: shahor02

Common/Field/include/Field/MagFieldFast.h

@@ -53,6 +53,7 @@ class MagFieldFast
   bool Field(const double xyz[3], double bxyz[3]) const;
   bool Field(const float xyz[3], float bxyz[3]) const;
   bool Field(const math_utils::Point3D<float> xyz, float bxyz[3]) const;
+  bool Field(const math_utils::Point3D<double> xyz, double bxyz[3]) const;
   bool GetBcomp(EDim comp, const double xyz[3], double& b) const;
   bool GetBcomp(EDim comp, const float xyz[3], float& b) const;
   bool GetBcomp(EDim comp, const math_utils::Point3D<float> xyz, double& b) const;

Common/Field/src/MagFieldFast.cxx

@@ -217,6 +217,22 @@ bool MagFieldFast::Field(const math_utils::Point3D<float> xyz, float bxyz[3]) co
   return true;
 }
 
+//_______________________________________________________________________
+bool MagFieldFast::Field(const math_utils::Point3D<double> xyz, double bxyz[3]) const
+{
+  // get field
+  int zSeg, rSeg, quadrant;
+  if (!GetSegment(xyz.X(), xyz.Y(), xyz.Z(), zSeg, rSeg, quadrant)) {
+    return false;
+  }
+  const SolParam* par = &mSolPar[rSeg][zSeg][quadrant];
+  bxyz[kX] = CalcPol(par->parBxyz[kX], xyz.X(), xyz.Y(), xyz.Z()) * mFactorSol;
+  bxyz[kY] = CalcPol(par->parBxyz[kY], xyz.X(), xyz.Y(), xyz.Z()) * mFactorSol;
+  bxyz[kZ] = CalcPol(par->parBxyz[kZ], xyz.X(), xyz.Y(), xyz.Z()) * mFactorSol;
+  //
+  return true;
+}
+
 //_______________________________________________________________________
 bool MagFieldFast::GetSegment(float x, float y, float z, int& zSeg, int& rSeg, int& quadrant) const
 {