set the maximum kinetic energy in the fit dynamically based on particle ID

The range and energy loss tables have different maximum values for electrons,
muons, and protons so we have to dynamically set the maximum energy of the fit
in order to avoid a GSL interpolation error.

This commit adds {electron,muon,proton}_get_max_energy() functions to return
the maximum energy in the tables and that is then used to set the maximum value
in the fit.

1 files changed, 15 insertions, 0 deletions

diff --git a/src/electron.c b/src/electron.c
index 5b4133a..3f5d4d8 100644
--- a/src/electron.c
+++ b/src/electron.c
@@ -365,6 +365,21 @@ err:
     return -1;
 }
 
+/* Returns the maximum kinetic energy for an electron in the range tables.
+ *
+ * If you call electron_get_range() or electron_get_dEdx() with a kinetic
+ * energy higher you will get a GSL interpolation error. */
+double electron_get_max_energy(void)
+{
+    if (!initialized) {
+        if (init()) {
+            exit(1);
+        }
+    }
+
+    return x[size-1];
+}
+
 double electron_get_range(double T, double rho)
 {
     /* Returns the approximate range a electron with kinetic energy `T` will travel