Protect against calculating exponents of large numbers

diff --git a/TTherminator/Therminator/Integrator.cxx b/TTherminator/Therminator/Integrator.cxx
index cdd7c30a45817a6a778f76b3a3b064d790cc020b..da5e84d56bbc6dba42e0ccc5fd7874aa2e35c65c 100644 (file)
--- a/TTherminator/Therminator/Integrator.cxx
+++ b/TTherminator/Therminator/Integrator.cxx
@@ -54,11 +54,13 @@ Integrator::Integrator(int aNpart)
 
 double Integrator::CalcBE(double aX)
 {
+  if (aX>200.0) return 0.0;
   return 1.0/(TMath::Exp(aX)-1);
 }
 
 double Integrator::CalcFD(double aX)
 {
+  if (aX>200.0) return 0.0;
   return 1.0/(TMath::Exp(aX)+1);
 }