Hello,
I have written a simulation to detect neutron spallation by colliding a 26GeV proton beam with a cylindrical Lead target. To my surprise, I get around 300 neutrons, each of around 1GeV, per event. Adding up all the energies together, the resulting energy is way higher than the energy of the incident particles. Is this because the simulation takes the proton interacting with the lead nucleus into account --ie. the binding energy?
I am using the FTFP_BERT physics list, with G4DecayPhysics and G4RadioactiveDecayPhysics.
Geant4 Version: 11.2.2
Operating System: Kali Linux (WSL2)
Compiler/Version: G++ 14.2.0
CMake Version: 3.30.5