Protect agains FPE in getXatLabR, redefine Almost1

Author: shahor02

Common/Constants/include/CommonConstants/MathConstants.h

@@ -23,7 +23,7 @@ namespace constants
 namespace math
 {
 constexpr float Almost0 = 1.17549e-38;
-constexpr float Almost1 = 1.f - Almost0;
+constexpr float Almost1 = 1.f - 1.0e-6;
 constexpr float VeryBig = 1.f / Almost0;
 
 constexpr float PI = 3.14159274101257324e+00f;

DataFormats/Reconstruction/src/TrackParametrization.cxx

@@ -585,6 +585,9 @@ GPUd() bool TrackParametrization<value_T>::getXatLabR(value_t r, value_t& x, val
   const auto fy = mP[0], sn = mP[2];
   const value_t kEps = 1.e-6;
   //
+  if (gpu::CAMath::Abs(getSnp()) > constants::math::Almost1) {
+    return false;
+  }
   auto crv = getCurvature(bz);
   while (gpu::CAMath::Abs(crv) > constants::math::Almost0) { // helix ?
     // get center of the track circle