diff options
Diffstat (limited to 'src/pdg.c')
| -rw-r--r-- | src/pdg.c | 5 |
1 files changed, 3 insertions, 2 deletions
@@ -1,7 +1,8 @@ #include "pdg.h" #include "math.h" +#include "sno.h" -double get_scattering_rms(double x, double p, double beta, double z, double rho) +double get_scattering_rms(double x, double p, double beta, double z) { /* Returns the RMS width of the scattering angle for a particle deflected * by many small-angle scatters after a distance `x`. `p` is the momentum @@ -22,5 +23,5 @@ double get_scattering_rms(double x, double p, double beta, double z, double rho) * See Equation 33.15 in * http://pdg.lbl.gov/2018/reviews/rpp2018-rev-passage-particles-matter.pdf. */ if (x == 0.0) return 0.0; - return (13.6/(beta*p))*z*sqrt(x*rho/RADIATION_LENGTH)*(1+0.038*log((x*z*z)/(RADIATION_LENGTH*beta*beta/rho))); + return (13.6/(beta*p))*z*sqrt(x*WATER_DENSITY/RADIATION_LENGTH)*(1+0.038*log((x*z*z)/(RADIATION_LENGTH*beta*beta/WATER_DENSITY))); } |