hardcode the density when computing dE/dx

Since we only have the range and dE/dx tables for light water for electrons and
protons it's not correct to use the heavy water density. Also, even though we
have both tables for muons, currently we only load the heavy water table, so we
hardcode the density to that of heavy water.

In the future, it would be nice to load both tables and use the correct one
depending on if we are fitting in the heavy or light water.

1 files changed, 3 insertions, 3 deletions

diff --git a/src/test.c b/src/test.c
index 20a66ec..c1da272 100644
--- a/src/test.c
+++ b/src/test.c
@@ -124,7 +124,7 @@ int test_proton_get_energy(char *err)
 
     /* Assume initial kinetic energy is 1 GeV. */
     T0 = 1000.0;
-    p = particle_init(IDP_PROTON, T0,1.0,10000);
+    p = particle_init(IDP_PROTON, T0, 10000);
     T = particle_get_energy(1e-9,p);
 
     /* At the beginning of the track we should have roughly the same energy. */
@@ -160,7 +160,7 @@ int test_electron_get_energy(char *err)
 
     /* Assume initial kinetic energy is 1 GeV. */
     T0 = 1000.0;
-    p = particle_init(IDP_E_MINUS, T0,1.0,10000);
+    p = particle_init(IDP_E_MINUS, T0, 10000);
     T = particle_get_energy(1e-9,p);
 
     /* At the beginning of the track we should have roughly the same energy. */
@@ -196,7 +196,7 @@ int test_muon_get_energy(char *err)
 
     /* Assume initial kinetic energy is 1 GeV. */
     T0 = 1000.0;
-    p = particle_init(IDP_MU_MINUS, T0,1.0,10000);
+    p = particle_init(IDP_MU_MINUS, T0, 10000);
     T = particle_get_energy(1e-9,p);
 
     /* At the beginning of the track we should have roughly the same energy. */