diff options
Diffstat (limited to 'src/likelihood.c')
| -rw-r--r-- | src/likelihood.c | 6 |
1 files changed, 3 insertions, 3 deletions
diff --git a/src/likelihood.c b/src/likelihood.c index c90914b..66e72e4 100644 --- a/src/likelihood.c +++ b/src/likelihood.c @@ -122,7 +122,7 @@ double getKineticEnergy(double x, double T0) return get_T(T0, x, HEAVY_WATER_DENSITY); } -double nll_muon(event *ev, double T0, double *pos, double *dir, double t0, double *z1, double *z2, size_t n) +double nll_muon(event *ev, double T0, double *pos, double *dir, double t0, double *z1, double *z2, size_t n, double epsrel) { size_t i, j; intParams params; @@ -202,14 +202,14 @@ double nll_muon(event *ev, double T0, double *pos, double *dir, double t0, doubl x = R*sin(theta_cerenkov-theta)/sin(theta_cerenkov); F.function = &gsl_muon_charge; - gsl_integration_cquad(&F, 0, range, 0, 1e-2, w, &result, &error, &nevals); + gsl_integration_cquad(&F, 0, range, 0, epsrel, w, &result, &error, &nevals); mu_direct[i] = result; total_charge += mu_direct[i]; if (mu_direct[i] > 0.001) { F.function = &gsl_muon_time; - gsl_integration_cquad(&F, 0, range, 0, 1e-2, w, &result, &error, &nevals); + gsl_integration_cquad(&F, 0, range, 0, epsrel, w, &result, &error, &nevals); ts[i] = result; ts[i] /= mu_direct[i]; |