Harmonize error treatment of PropagationUtility::GetUpperLim

Fixes #387
tudo-astroparticlephysics · Jan 31, 2024 · 87e74fb · 87e74fb
1 parent c3b01b7
commit 87e74fb
Showing 1 changed file with 5 additions and 1 deletion.
diff --git a/src/PROPOSAL/detail/PROPOSAL/propagation_utility/PropagationUtilityIntegral.cxx b/src/PROPOSAL/detail/PROPOSAL/propagation_utility/PropagationUtilityIntegral.cxx
@@ -26,7 +26,11 @@ double UtilityIntegral::GetUpperLimit(double energy_initial, double rnd)
     auto sum = integral.IntegrateWithRandomRatio(
         energy_initial, lower_lim, FunctionToIntegral, 4, -rnd);
 
-    assert(sum > rnd); // searched energy is below lower_lim
+    if (rnd > sum) {
+        throw std::logic_error("Unable to calculate GetUpperLimit since result"
+                               "is below lower_lim. rnd was " + std::to_string(rnd)
+                               + " with rnd_max " + std::to_string(sum));
+    }
     (void)sum;
 
     return integral.GetUpperLimit();