diff options
Diffstat (limited to 'src/muon.c')
| -rw-r--r-- | src/muon.c | 14 |
1 files changed, 2 insertions, 12 deletions
@@ -216,9 +216,9 @@ double get_dEdx(double T, double rho) return gsl_spline_eval(spline_dEdx, T, acc_dEdx)/rho; } -double get_expected_charge(double x, double T, double T0, double *pos, double *dir, double *pmt_pos, double *pmt_normal, double r) +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, theta0, E, p, beta, z, rho, R, E0, p0, beta0; + double pmt_dir[3], cos_theta, n, wavelength0, omega, E, p, beta, z, R; z = 1.0; @@ -244,23 +244,13 @@ double get_expected_charge(double x, double T, double T0, double *pos, double *d wavelength0 = 400.0; if (R <= AV_RADIUS) { - rho = HEAVY_WATER_DENSITY; n = get_index_snoman_d2o(wavelength0); } else { - rho = WATER_DENSITY; n = get_index_snoman_h2o(wavelength0); } if (beta < 1/n) return 0; - /* Calculate total energy */ - E0 = T0 + MUON_MASS; - p0 = sqrt(E0*E0 - MUON_MASS*MUON_MASS); - beta0 = p0/E0; - - /* FIXME: is this formula valid for muons? */ - theta0 = get_scattering_rms(x,p0,beta0,z,rho); - /* FIXME: add angular response and scattering/absorption. */ return omega*FINE_STRUCTURE_CONSTANT*z*z*(1-(1/(beta*beta*n*n)))*get_probability(beta, cos_theta, theta0); } |