aboutsummaryrefslogtreecommitdiff
path: root/src/likelihood.c
diff options
context:
space:
mode:
authortlatorre <tlatorre@uchicago.edu>2019-03-26 16:15:35 -0500
committertlatorre <tlatorre@uchicago.edu>2019-03-26 16:15:35 -0500
commitdef9e9168dc3f3e3c70ff4bf61d176ed28ca6184 (patch)
tree60c8cacd7cac61b609da2e4b6be361cf584af9c1 /src/likelihood.c
parentc92657e33b81e37eb5851217fcdea8767c183e62 (diff)
downloadsddm-def9e9168dc3f3e3c70ff4bf61d176ed28ca6184.tar.gz
sddm-def9e9168dc3f3e3c70ff4bf61d176ed28ca6184.tar.bz2
sddm-def9e9168dc3f3e3c70ff4bf61d176ed28ca6184.zip
forgot to add scattering from acrylic to get_expected_charge()
Diffstat (limited to 'src/likelihood.c')
-rw-r--r--src/likelihood.c2
1 files changed, 1 insertions, 1 deletions
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);