In the reaction alpha (5.4 MeV) on Be9, ParticleHP does not produce any neutron. Below a macro to run Hadr03 and the results for ParticleHPInelastic and BinaryLightIonCascade ---------------------------------------------------------------- /control/verbose 2 /run/verbose 1 # /testhadr/det/setIsotopeMat Be9 4 9 1.85 g/cm3 # /run/initialize # /process/list /process/inactivate hadElastic # /gun/particle alpha /gun/energy 5.4 MeV # /run/printProgress 10000 /run/beamOn 100000 ************************************************************* Geant4 version Name: geant4-10-02-ref-06 (30-June-2016) ************************************************************* Run Summary Number of events processed : 100000 User=176.57s Real=177.08s Sys=0.15s The run is 100000 alpha of 5.4 MeV through 10 m of Be9 (density: 1.85 g/cm3 ) Process calls frequency: alphaInelastic= 100000 MeanFreePath: 14.833 cm +- 14.795 cm massic: 27.441 g/cm2 CrossSection: 0.067417 cm^-1 massic: 3.6442 mm2/g crossSection per atom: 545.35 mbarn Verification: crossSections from G4HadronicProcessStore: alphaInelastic= 3.6669 mm2/g 548.76 mbarn total= 3.6669 mm2/g 548.76 mbarn List of nuclear reactions: alpha + Be9 --> N gamma + C13: 100000 Q = 10.647 MeV Number of gamma: N = 2 --> 3 List of generated particles: C13: 100000 Emean = 1.6632 MeV ( 1.4732 MeV --> 1.6645 MeV) gamma: 200002 Emean = 7.1917 MeV ( 692.57 keV --> 13.189 MeV) ************************************************************* Geant4 version Name: geant4-10-02-ref-06 (30-June-2016) ************************************************************* Run Summary Number of events processed : 100000 User=3.54s Real=3.62s Sys=0.01s The run is 100000 alpha of 5.4 MeV through 10 m of Be9 (density: 1.85 g/cm3 ) Process calls frequency: alphaInelastic= 100000 MeanFreePath: 9.3076 cm +- 9.2834 cm massic: 17.219 g/cm2 CrossSection: 0.10744 cm^-1 massic: 5.8075 mm2/g crossSection per atom: 869.1 mbarn Verification: crossSections from G4HadronicProcessStore: alphaInelastic= 5.7697 mm2/g 863.44 mbarn total= 5.7697 mm2/g 863.44 mbarn List of nuclear reactions: alpha + Be9 --> alpha + Be9: 49586 Q = -4.3914e-07 eV alpha + Be9 --> neutron + 3 alpha: 179 Q = -1.5731 MeV alpha + Be9 --> neutron + C12: 10743 Q = 5.7006 MeV alpha + Be9 --> neutron + N gamma + C12: 39492 Q = 5.7006 MeV Number of gamma: N = 1 --> 1 List of generated particles: Be9: 49586 Emean = 2.2954 MeV ( 197.95 eV --> 4.6003 MeV) C12: 50235 Emean = 1.9938 MeV ( 144.56 keV --> 4.3883 MeV) alpha: 50123 Emean = 3.0811 MeV ( 18.152 keV --> 5.3998 MeV) gamma: 39492 Emean = 4.4391 MeV ( 4.3293 MeV --> 4.5496 MeV) neutron: 50414 Emean = 5.6009 MeV ( 70.146 keV --> 10.956 MeV)
Hi Michel, As you probably know, I've taken over Tatsumi's bug reports since he is no longer at SLAC. That's why it's taking me some time to get them done. For this one, I checked with Pedro Arce, who says things seem to be OK with G4TENDL1.3.2/Alpha: neutron yield data for 5.4 MeV alpha on Be9 exist and seem to be correct. Are you using this version of TENDL, and Geant4 10.4 or later? If so, and if this is still a problem for you, I'll start looking at the code. Dennis
Hello Dennis, In Geant4 10.5 neutrons are there. But with the usual behavior of ParticleHP : reasonable neutrons spectrum but many spurious nuclear channels if we do_not_adjust_final_state, and the opposite if we adjust_final_state.
Hi Michel, I'm going to move this problem to development. It is the same problem we have had for a long time, and there is no easy solution. One thing we might do is split HP into "database mode" in which we get the right spectrum, but spurious channels, and "conservation mode" in which we get the correct channels but the wrong spectrum. Right now, I don't see a solution in which everyone is happy. I'm proposing that we have a discussion at Wednesday's hadronic group meeting. Dennis
