diff options
Diffstat (limited to 'src/likelihood.c')
| -rw-r--r-- | src/likelihood.c | 10 |
1 files changed, 1 insertions, 9 deletions
diff --git a/src/likelihood.c b/src/likelihood.c index ca0b2b6..ca73902 100644 --- a/src/likelihood.c +++ b/src/likelihood.c @@ -165,7 +165,7 @@ void particle_free(particle *p) static double get_expected_charge_shower(double *pos, double *dir, double *pmt_pos, double *pmt_normal, double r, double *reflected, double n_d2o, double n_h2o, double l_d2o, double l_h2o, double alpha, double beta, double rad) { - double pmt_dir[3], cos_theta, n, omega, R, f, f_reflec, cos_theta_pmt, absorption_length_h2o, absorption_length_d2o, absorption_length_acrylic, l_acrylic, theta_pmt, charge; + double pmt_dir[3], cos_theta, omega, f, f_reflec, cos_theta_pmt, absorption_length_h2o, absorption_length_d2o, absorption_length_acrylic, l_acrylic, theta_pmt, charge; SUB(pmt_dir,pmt_pos,pos); @@ -182,14 +182,6 @@ static double get_expected_charge_shower(double *pos, double *dir, double *pmt_p omega = get_solid_angle_fast(pos,pmt_pos,pmt_normal,r); - R = NORM(pos); - - if (R <= AV_RADIUS) { - n = n_d2o; - } else { - n = n_h2o; - } - theta_pmt = acos(-cos_theta_pmt); f_reflec = get_weighted_pmt_reflectivity(theta_pmt); f = get_weighted_pmt_response(theta_pmt); |