1 files changed, 15 insertions, 0 deletions
diff --git a/src/proton.c b/src/proton.c
index 664d930..4714d76 100644
--- a/src/proton.c
+++ b/src/proton.c
@@ -165,6 +165,21 @@ err:
     return -1;
 }
 
+/* Returns the maximum kinetic energy for a proton in the range tables.
+ *
+ * If you call proton_get_range() or proton_get_dEdx() with a kinetic
+ * energy higher you will get a GSL interpolation error. */
+double proton_get_max_energy(void)
+{
+    if (!initialized) {
+        if (init()) {
+            exit(1);
+        }
+    }
+
+    return x[size-1];
+}
+
 double proton_get_range(double T, double rho)
 {
     /* Returns the approximate range a proton with kinetic energy `T` will travel