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, 1 insertions, 0 deletions

diff --git a/src/muon.h b/src/muon.h
index a8d2ea1..8e1d39e 100644
--- a/src/muon.h
+++ b/src/muon.h
@@ -30,6 +30,7 @@
  * FIXME: Actually determine what this is. */
 #define MUON_PHOTONS_PER_MEV 7368.0
 
+double muon_get_max_energy(void);
 double muon_get_shower_photons(double T0, double rad);
 void muon_get_position_distribution_parameters(double T0, double *a, double *b);
 double muon_get_angular_distribution_alpha(double T0);