| Summary: | Nuclear interaction length for liquid deuterium and hydrogen is incorrect | ||
|---|---|---|---|
| Product: | Geant4 | Reporter: | Bryan Kerns <BryanMKerns> |
| Component: | materials | Assignee: | Marc.Verderi |
| Status: | RESOLVED DUPLICATE | ||
| Severity: | normal | CC: | BryanMKerns |
| Priority: | P5 | ||
| Version: | 9.6 | ||
| Hardware: | All | ||
| OS: | All | ||
Same problem as 1563. *** This problem has been marked as a duplicate of problem 1563 *** |
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 %