When I print out the material table, or when I use the GetNuclearInterLength() function of the liquidHydrogen and liquidDeuterium materials, I receive the wrong nuclear interaction length. The value is correct within about 1% for other materials, such as iron. According to the particle data group, the nuclear interaction length for liquid hydrogen should be 7.346 m and 424.7 m for liquid deuterium. Below is the output of printing the material table for liquidHydrogen and liquidDeuterium, where the values are incorrect by a large amount. Material: liquidDeuterium density: 169.000 mg/cm3 RadL: 7.454 m Nucl.Int.Length: 2.628 m Imean: 19.200 eV ---> Element: liquidDeuterium (liqD) Z = 1.0 N = 2.0 A = 2.01 g/mole ---> Isotope: liqD1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.99 % ---> Isotope: liqD2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.01 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: liquidHydrogen density: 70.800 mg/cm3 RadL: 8.923 m Nucl.Int.Length: 4.993 m Imean: 19.200 eV ---> Element: Hydrogen (H) Z = 1.0 N = 1.0 A = 1.01 g/mole ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.99 % ---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.01 % ElmMassFraction: 100.00 % ElmAbundance 100.00 %
Same problem as 1563.
*** This problem has been marked as a duplicate of problem 1563 ***