diff options
| -rw-r--r-- | src/likelihood.c | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/src/likelihood.c b/src/likelihood.c index 34f379a..245c20c 100644 --- a/src/likelihood.c +++ b/src/likelihood.c @@ -176,7 +176,7 @@ particle *particle_init(int id, double T0, size_t n) p->x[i] = p->range*i/(n-1); dEdx = muon_get_dEdx(p->T[i-1], HEAVY_WATER_DENSITY); p->T[i] = p->T[i-1] - dEdx*(p->x[i]-p->x[i-1]); - dEdx = muon_get_dEdx_rad(p->T[i-1], WATER_DENSITY); + dEdx = muon_get_dEdx_rad(p->T[i-1], HEAVY_WATER_DENSITY); rad += dEdx*(p->x[i]-p->x[i-1]); } /* Make sure that the energy is zero at the last step. This is so that |