From def9e9168dc3f3e3c70ff4bf61d176ed28ca6184 Mon Sep 17 00:00:00 2001 From: tlatorre Date: Tue, 26 Mar 2019 16:15:35 -0500 Subject: forgot to add scattering from acrylic to get_expected_charge() --- src/likelihood.c | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) (limited to 'src') diff --git a/src/likelihood.c b/src/likelihood.c index c53e1fc..6f79325 100644 --- a/src/likelihood.c +++ b/src/likelihood.c @@ -354,7 +354,7 @@ static void get_expected_charge(double beta, double theta0, double *pos, double * Technically we should compute this conditionally on the probability that * a photon is not absorbed, but since the probability of scattering is * pretty low, this is expected to be a very small effect. */ - prob_sct = 1.0 - get_fsct_d2o(l_d2o)*get_fsct_h2o(l_h2o); + prob_sct = 1.0 - get_fsct_d2o(l_d2o)*get_fsct_h2o(l_h2o)*get_fsct_acrylic(AV_THICKNESS); charge = omega*FINE_STRUCTURE_CONSTANT*z*z*(1-(1/(beta*beta*n*n)))*get_probability(beta, cos_theta, theta0); -- cgit