diff options
| author | tlatorre <tlatorre@uchicago.edu> | 2019-03-25 19:29:25 -0500 |
|---|---|---|
| committer | tlatorre <tlatorre@uchicago.edu> | 2019-03-25 19:29:25 -0500 |
| commit | bf60d08d517e7887417f0aa4068b726a8c749e58 (patch) | |
| tree | 30425882b75dc154e406457fa70e55f359085c8d /src/likelihood.c | |
| parent | 40e79187eb1037966fc723291936b96b7847f4fb (diff) | |
| download | sddm-bf60d08d517e7887417f0aa4068b726a8c749e58.tar.gz sddm-bf60d08d517e7887417f0aa4068b726a8c749e58.tar.bz2 sddm-bf60d08d517e7887417f0aa4068b726a8c749e58.zip | |
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.
Diffstat (limited to 'src/likelihood.c')
| -rw-r--r-- | src/likelihood.c | 6 |
1 files changed, 3 insertions, 3 deletions
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; |