set the starting energy to MAX_ENERGY if it's greater

Also increase the maximum kinetic energy to 10^4 GeV which is approximately the
maximum expected energy for cosmic muons at SNO.

1 files changed, 6 insertions, 2 deletions

diff --git a/src/fit.c b/src/fit.c
index e4e054f..cb8409b 100644
--- a/src/fit.c
+++ b/src/fit.c
@@ -49,8 +49,11 @@
 #define MAX_PARS 100
 /* Maximum number of peaks to search for in Hough transform. */
 #define MAX_NPEAKS 10
-/* Maximum kinetic energy for any particle. */
-#define MAX_ENERGY 10000
+/* Maximum kinetic energy for any particle. Currently set to 10^4 GeV which is
+ * approximately the energy of the most energetic muons in SNO (see the paper
+ * "Measurement of the Cosmic Ray and Neutrino-Induced Muon Flux at the Sudbury
+ * Neutrino Observatory".*/
+#define MAX_ENERGY 1e7
 
 char *GitSHA1(void);
 char *GitDirty(void);
@@ -5388,6 +5391,7 @@ int fit_event2(event *ev, double *xopt, double *fmin, int *id, size_t n, double
              * threshold. */
             double Tmin2 = sqrt(mass*mass/(1-pow(0.9,2)))-mass;
             if (T0 < Tmin2) T0 = Tmin2;
+            if (T0 > MAX_ENERGY) T0 = MAX_ENERGY;
 
             x[4+3*j] = T0;
             ss[4+3*j] = x[4+3*j]*0.1;