I am currently creating an Sb-124 BeO photo neutron source. This is built upon the Neutron source example in geant4 version 10.06.
I am using QGSP_BIC_HP for inelastic hadron physics and for the GammaNuclearPhysics I use:
G4ProcessManager* pManager = G4Gamma::Gamma()->GetProcessManager();
G4PhotoNuclearProcess* process = new G4PhotoNuclearProcess();
G4LENDInelastic* lend = new G4LENDInelastic(G4Gamma::Gamma());
G4LENDCombinedCrossSection lendXS = new G4LENDCombinedCrossSection(G4Gamma::Gamma());
While I expect to get the emerging neutrons to have a peak energy around 24 keV, I am getting a peak around 0.1 MeV. Am I missing something when defining the physics processes?