diff options
| author | tlatorre <tlatorre@uchicago.edu> | 2019-07-16 11:17:48 -0500 |
|---|---|---|
| committer | tlatorre <tlatorre@uchicago.edu> | 2019-07-16 11:17:48 -0500 |
| commit | 9f63d7eae7d8be9b599955644aa4c7a92ccb640f (patch) | |
| tree | 1e27e2df25c13dff0664993945c9e062d2d640c7 /src | |
| parent | 6351b9de6c5b96da26d33bb1c303f092194987d7 (diff) | |
| download | sddm-9f63d7eae7d8be9b599955644aa4c7a92ccb640f.tar.gz sddm-9f63d7eae7d8be9b599955644aa4c7a92ccb640f.tar.bz2 sddm-9f63d7eae7d8be9b599955644aa4c7a92ccb640f.zip | |
multiply rayleigh scattering by 2 to account for the 50% PMT coverage
Diffstat (limited to 'src')
| -rw-r--r-- | src/likelihood.c | 15 |
1 files changed, 12 insertions, 3 deletions
diff --git a/src/likelihood.c b/src/likelihood.c index 90fdd82..2681475 100644 --- a/src/likelihood.c +++ b/src/likelihood.c @@ -519,7 +519,10 @@ static void get_expected_charge_shower(particle *p, double *pos, double *dir, in charge = get_weighted_quantum_efficiency()*omega/(2*M_PI); - *reflected = (1.0-prob_abs)*(1.0-prob_sct)*f_reflec*charge + prob_sct*charge; + /* Note: We multiply the amount of Rayleigh scattered light here by 2 since + * we are only calculating the fraction of Rayleigh scattered light that + * would hit the PMT and the coverage in SNO+ is about 50%. */ + *reflected = (1.0-prob_abs)*(1.0-prob_sct)*f_reflec*charge + 2*prob_sct*charge; *t = (l_d2o*get_avg_index_d2o() + l_h2o*get_avg_index_h2o())/SPEED_OF_LIGHT; @@ -605,7 +608,10 @@ static void get_expected_charge(double beta, double theta0, double *pos, double charge = omega*(1-cos_theta_cerenkov*cos_theta_cerenkov)*get_probability(beta, cos_theta, sin_theta, theta0); - *reflected = (1.0-prob_abs)*(1.0-prob_sct)*f_reflec*charge + prob_sct*charge; + /* Note: We multiply the amount of Rayleigh scattered light here by 2 since + * we are only calculating the fraction of Rayleigh scattered light that + * would hit the PMT and the coverage in SNO+ is about 50%. */ + *reflected = (1.0-prob_abs)*(1.0-prob_sct)*f_reflec*charge + 2*prob_sct*charge; *q = (1.0-prob_abs)*(1.0-prob_sct)*f*charge; } @@ -924,7 +930,10 @@ static double get_total_charge_approx(double beta0, double *pos, double *dir, pa if (p->delta_ray_photons > 0) charge += get_weighted_quantum_efficiency()*p->delta_ray_photons*electron_get_angular_pdf_delta_ray(cos_theta,p->delta_ray_a,p->delta_ray_b,1.0/get_avg_index_d2o())*omega/(2*M_PI); - *mu_reflected = (1.0-prob_abs)*(1.0-prob_sct)*f_reflec*charge + prob_sct*charge; + /* Note: We multiply the amount of Rayleigh scattered light here by 2 since + * we are only calculating the fraction of Rayleigh scattered light that + * would hit the PMT and the coverage in SNO+ is about 50%. */ + *mu_reflected = (1.0-prob_abs)*(1.0-prob_sct)*f_reflec*charge + 2*prob_sct*charge; return (1.0-prob_abs)*(1.0-prob_sct)*f*charge; } |