diff options
| author | tlatorre <tlatorre@uchicago.edu> | 2018-08-31 14:14:40 -0500 |
|---|---|---|
| committer | tlatorre <tlatorre@uchicago.edu> | 2018-08-31 14:14:40 -0500 |
| commit | 4cb1e319dce25a93719994587cbbe483426536c4 (patch) | |
| tree | fbe00e8422bdf3178eb41c78c14e91865d33bba0 /src/likelihood.c | |
| parent | 7478c6eab0b618817213893794b9ada1d836ac1c (diff) | |
| download | sddm-4cb1e319dce25a93719994587cbbe483426536c4.tar.gz sddm-4cb1e319dce25a93719994587cbbe483426536c4.tar.bz2 sddm-4cb1e319dce25a93719994587cbbe483426536c4.zip | |
add epsrel argument to likelihood function
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]; |