diff --git a/catima.cpp b/catima.cpp
index 2b9809a..9c4b0da 100644
--- a/catima.cpp
+++ b/catima.cpp
@@ -211,7 +211,7 @@ double energy_out(double T, double thickness, const Interpolator &range_spline){
     e = T - (thickness*dedx);
     while(1){
         r = range - range_spline(e) - thickness;
-        if(fabs(r)<Eout_epsilon)return e;
+        if(fabs(r)<Eout_th_epsilon)return e;
         double step = -r*dedx;
         e = e-step;
         if(e<Ezero)return 0.0;
diff --git a/constants.h b/constants.h
index 391d27a..9c3a77c 100644
--- a/constants.h
+++ b/constants.h
@@ -11,7 +11,7 @@ constexpr double logEmax = 7.0;  // log of max energy
 constexpr int max_datapoints = 600; // how many datapoints between logEmin and logEmax
 constexpr int max_storage_data = 60; // number of datapoints which can be stored in cache
 constexpr double numeric_epsilon = 10*std::numeric_limits<double>::epsilon();
-constexpr double Eout_epsilon = 1e-5;  //
+constexpr double Eout_th_epsilon = 1e-5;  //
 
 constexpr double thin_target_limit = 1 - 1e-3;