From efa1929196659ca24eb7d40e9b7532c16fb5e20a Mon Sep 17 00:00:00 2001 From: tlatorre Date: Mon, 1 Oct 2018 16:03:32 -0500 Subject: add absorption length for acrylic --- src/muon.c | 7 +++++-- 1 file changed, 5 insertions(+), 2 deletions(-) (limited to 'src/muon.c') diff --git a/src/muon.c b/src/muon.c index 09f1dd1..575158b 100644 --- a/src/muon.c +++ b/src/muon.c @@ -258,7 +258,7 @@ double get_dEdx(double T, double rho) double get_expected_charge(double x, double T, double theta0, double *pos, double *dir, double *pmt_pos, double *pmt_normal, double r) { - double pmt_dir[3], cos_theta, n, wavelength0, omega, E, p, beta, z, R, f, cos_theta_pmt, absorption_length_h2o, absorption_length_d2o, l_h2o, l_d2o; + double pmt_dir[3], cos_theta, n, wavelength0, omega, E, p, beta, z, R, f, cos_theta_pmt, absorption_length_h2o, absorption_length_d2o, absorption_length_acrylic, l_h2o, l_d2o, l_acrylic; z = 1.0; @@ -298,8 +298,11 @@ double get_expected_charge(double x, double T, double theta0, double *pos, doubl absorption_length_d2o = get_weighted_absorption_length_snoman_d2o(); absorption_length_h2o = get_weighted_absorption_length_snoman_h2o(); + absorption_length_acrylic = get_weighted_absorption_length_snoman_acrylic(); get_path_length(pos,pmt_pos,AV_RADIUS,&l_d2o,&l_h2o); - return f*exp(-l_d2o/absorption_length_d2o-l_h2o/absorption_length_h2o)*omega*FINE_STRUCTURE_CONSTANT*z*z*(1-(1/(beta*beta*n*n)))*get_probability(beta, cos_theta, theta0); + l_acrylic = AV_RADIUS_OUTER - AV_RADIUS_INNER; + + return f*exp(-l_d2o/absorption_length_d2o-l_h2o/absorption_length_h2o-l_acrylic/absorption_length_acrylic)*omega*FINE_STRUCTURE_CONSTANT*z*z*(1-(1/(beta*beta*n*n)))*get_probability(beta, cos_theta, theta0); } -- cgit