From 89b3fd9989c8e19fe1ffd4d18e0c2e17cdee790a Mon Sep 17 00:00:00 2001 From: tlatorre Date: Fri, 14 Jun 2019 12:26:10 -0500 Subject: 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. --- src/fit.c | 8 ++++++-- 1 file changed, 6 insertions(+), 2 deletions(-) (limited to 'src') 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; -- cgit