From 9641594eb070aaa4af399a311927b159e22e62b0 Mon Sep 17 00:00:00 2001 From: tlatorre Date: Tue, 18 Sep 2018 18:44:03 -0500 Subject: speed likelihood calculation up a bit --- src/muon.c | 3 +-- 1 file changed, 1 insertion(+), 2 deletions(-) (limited to 'src/muon.c') diff --git a/src/muon.c b/src/muon.c index 0ab5cd8..9114dda 100644 --- a/src/muon.c +++ b/src/muon.c @@ -303,6 +303,5 @@ double get_expected_charge(double x, double T, double theta0, double *pos, doubl get_path_length(pos,pmt_pos,AV_RADIUS,&l_d2o,&l_h2o); - /* FIXME: add angular response and scattering/absorption. */ - return f*exp(-l_d2o/absorption_length_d2o)*exp(-l_h2o/absorption_length_h2o)*omega*FINE_STRUCTURE_CONSTANT*z*z*(1-(1/(beta*beta*n*n)))*get_probability(beta, cos_theta, theta0); + return f*exp(-l_d2o/absorption_length_d2o-l_h2o/absorption_length_h2o)*omega*FINE_STRUCTURE_CONSTANT*z*z*(1-(1/(beta*beta*n*n)))*get_probability(beta, cos_theta, theta0); } -- cgit