From 4372a0bb54b374061e1b00bca9da45aa0325e7c5 Mon Sep 17 00:00:00 2001 From: tlatorre Date: Sun, 25 Nov 2018 11:56:26 -0600 Subject: add a separate `dx_shower` parameter for the spacing of the shower track integral --- src/likelihood.c | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) (limited to 'src/likelihood.c') diff --git a/src/likelihood.c b/src/likelihood.c index d0ab20f..dcf681a 100644 --- a/src/likelihood.c +++ b/src/likelihood.c @@ -705,7 +705,7 @@ static double getKineticEnergy(double x, void *p) return particle_get_energy(x, (particle *) p); } -double nll_muon(event *ev, int id, double T0, double *pos, double *dir, double t0, double *z1, double *z2, size_t n, double dx, int fast) +double nll_muon(event *ev, int id, double T0, double *pos, double *dir, double t0, double *z1, double *z2, size_t n, double dx, double dx_shower, int fast) { /* Returns the negative log likelihood for event `ev` given a particle with * id `id`, initial kinetic energy `T0`, position `pos`, direction `dir` and @@ -790,7 +790,7 @@ double nll_muon(event *ev, int id, double T0, double *pos, double *dir, double t /* Number of points to sample along the longitudinal direction for * the electromagnetic shower. */ - npoints_shower = (size_t) ((xhi-xlo)/dx + 0.5); + npoints_shower = (size_t) ((xhi-xlo)/dx_shower + 0.5); if (npoints_shower < MIN_NPOINTS) npoints_shower = MIN_NPOINTS; -- cgit