update theta0 calculation to use the radiation length in light water

Previously I was using the radiation length in light water but scaling it by
the density of heavy water, which isn't correct. Since the radiation length in
heavy and light water is almost identical, we just use the radiation length in
light water.

1 files changed, 7 insertions, 2 deletions

diff --git a/src/pdg.h b/src/pdg.h
index 98046f5..df73380 100644
--- a/src/pdg.h
+++ b/src/pdg.h
@@ -2,13 +2,18 @@
 #define PDG_H
 
 #define SPEED_OF_LIGHT 29.9792458 /* cm/ns */
-/* From http://pdg.lbl.gov/2017/AtomicNuclearProperties/HTML/water_liquid.html */
+/* From http://pdg.lbl.gov/2017/AtomicNuclearProperties/HTML/water_liquid.html.
+ *
+ * Technically there is a different radiation length for D2O and H2O, and each
+ * should be divided by the density of the water. However, the results for
+ * water and heavy water are very close so instead of having to keep track of
+ * which medium we are in we just use the radiation length in light water. */
 #define RADIATION_LENGTH 36.08 /* g/cm^2 */
 #define ELECTRON_MASS 0.5109989461 /* MeV */
 #define MUON_MASS 105.6583745 /* MeV */
 #define PROTON_MASS 938.272081 /* MeV */
 #define FINE_STRUCTURE_CONSTANT 7.297352566417e-3
 
-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);
 
 #endif