Please fill out the following information to help in answering your question, and also see tips for posting code snippets. If you don’t provide this information it will take more time to help with your problem!
Geant4 Version: geant4-v11.3.0
Operating System: ubuntu
Compiler/Version:
CMake Version:
I am writing to discuss an issue I am encountering in my simulation study of the 7Li(p,n)7Be^7\mathrm{Li}(p,n)^7\mathrm{Be}7Li(p,n)7Be reaction. I have performed simulations using the physics list G4HadronPhysicsQGSP_BIC_AllHP along with G4IonPhysicsPHP for a proton beam energy of 2.13 MeV.
Subsequently, I converted the neutron data obtained from direct kinematics into inverse kinematics, ensuring the same center-of-mass (COM) energy and adopting the same N_t N_\phi normalization as implemented in the ORSAY experiment. However, I observe that the neutron counts in the angular distribution are lower than those reported in the ORSAY experimental results, which were obtained using inverse kinematics for the p(7Li,7Be)np(^7\mathrm{Li},^7\mathrm{Be})np(7Li,7Be)n reaction.
In the ORSAY experiment, a 15 MeV 7Li^7\mathrm{Li}7Li beam was incident on a polypropylene (proton) target. Since the cross-section in the center-of-mass frame should be invariant between direct and inverse kinematics, I would expect the converted angular distributions to show reasonable agreement. However, a clear discrepancy persists, particularly in the count rates.
I would be grateful if you could suggest possible reasons for this mismatch. Any guidance or insight on this matter would be extremely helpful in improving the accuracy of my analysis.