Problem 2236

When users try to use G4QMD model based on JQMD(JAERI QMD) for proton, neutron and pions, there are two types of error messages.

1. When G4QMDReaction is registered in the PhysicsList, other nuclei are good, but not for proton and neutron. If users register the QMD model for proton or neutron, GEANT4 gives us CRITICAL ERROR MESSAGE about Energy non-conservation. I looked up the source code of G4QMD model, I found that the mass number is not conserved. See the code "G4QMDGroundStateNucleus.cc", line 105 and 110. When G4QMDGroundStateNucleus has been initialized, G4QMDGroundStateNucleus sets participants at the beginning of the loop(line 74). After the loop, in the case of proton and neutron, it sets participants one more. I think it should be fixed.

2. For the pions, we also cannot use the QMD model. See the code "G4QMDReaction.cc" line 62. It says when users use the QMD model for pions, GEANT4 gives an error message "There is no GetIsoCrossSection", which G4PiNuclearCrossSection doesn't have. I think it also should be changed to other PionCrossSection classes. I recommend that it would be better to use G4BGGPionElasticXS and G4BGGPionInelasticXS. G4QMDReaction only uses total cross-section to calculate "bmax_0"(line 157).

I was trying to use the G4QMD model for neutron and pions and I found these. I am currently using the fixed GEANT, and hope that these errors should be fixed ASAP.