From bf60d08d517e7887417f0aa4068b726a8c749e58 Mon Sep 17 00:00:00 2001 From: tlatorre Date: Mon, 25 Mar 2019 19:29:25 -0500 Subject: update rayleigh scattering calculation This commit updates the optics code to calculate the rayleigh scattering length using the Einstein-Smoluchowski formula instead of using the effective rayleigh scattering lengths from the RSPR bank. --- src/likelihood.c | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) (limited to 'src/likelihood.c') diff --git a/src/likelihood.c b/src/likelihood.c index b0f6b94..442a8b3 100644 --- a/src/likelihood.c +++ b/src/likelihood.c @@ -235,7 +235,7 @@ static void get_expected_charge_delta_ray(particle *p, double *pos, double *dir, get_path_length(pos,pmts[pmt].pos,AV_RADIUS,&l_d2o,&l_h2o); prob_abs = 1.0 - get_fabs_d2o(l_d2o)*get_fabs_h2o(l_h2o)*get_fabs_acrylic(AV_THICKNESS); - 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); constant = get_weighted_quantum_efficiency()*omega/(2*M_PI); @@ -281,7 +281,7 @@ static void get_expected_charge_shower(particle *p, double *pos, double *dir, in get_path_length(pos,pmts[pmt].pos,AV_RADIUS,&l_d2o,&l_h2o); prob_abs = 1.0 - get_fabs_d2o(l_d2o)*get_fabs_h2o(l_h2o)*get_fabs_acrylic(AV_THICKNESS); - 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); constant = get_weighted_quantum_efficiency()*omega/(2*M_PI); @@ -730,7 +730,7 @@ static double get_total_charge_approx(double beta0, double *pos, double *dir, pa get_path_length(tmp,pmts[i].pos,AV_RADIUS,&l_d2o,&l_h2o); prob_abs = 1.0 - get_fabs_d2o(l_d2o)*get_fabs_h2o(l_h2o)*get_fabs_acrylic(AV_THICKNESS); - 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); /* Assume the particle is travelling at the speed of light. */ *t = s/SPEED_OF_LIGHT + l_d2o*avg_index_d2o/SPEED_OF_LIGHT + l_h2o*avg_index_h2o/SPEED_OF_LIGHT; -- cgit