add a separate `dx_shower` parameter for the spacing of the shower track integral

1 files changed, 2 insertions, 2 deletions

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;