I am trying to study the effect of water composition on the Bragg peak for heavy ions. I have noticed that implementing the “G4_WATER” gives the accurate result compared to experimental data but I cannot get the same result with a user defined water like below:
G4Material* H2O = new G4Material(“Water”, 1.000g/cm3, nelements=2);
Additionally, I can see the effect of density variation but the effect of “I_mean” variation is negligible. I was wondering if this is the correct way to define the “I_mean” of water.