I am trying to simulate a 40 MeV deuteron beam onto a thick graphite target.
the simulation I have currently seems to be producing much less neutrons than expected.
So I am wondering if potentially my physics are incorrect.
Currently I have:
I think have figured out the issue with the 95+% of the processes called being ionIoni…
But I am still not producing the expected amount of neutrons.
Is there any more ideal physics lists to use for deuterons?
Did you try G4IonINCLXXPhysics ?
Also, if you have downloaded G4TENDL data and set env variable G4PARTICLEHPDATA to point on it, you can try to play with G4IonPhysicsPHP.
I do not know if they are " better " or not than QMD …
I have tried G4IonINCLXXPhysics before with similar results but I may try it again.
I have had the G4TENDL data but i have the individual G4 HPData for each particle assigned.
i.e. I have G4ProtonHPData pointing to …/G4TENDL1.3.2/Proton, with alpha, he3, deuteron, trition all pointing to there respective Tendl Folders.
Would it make a difference to have the G4ParticleHPData instead of the particles separate??
=============================================================================
INCL++ WARNING
INCL++/G4ExcitationHandler could not use its own level-density parameter for fission
=============================================================================
It is just a warning but it repeats this multiple times. is this warning an issue?
Hello Davide, thank you for the confirmation!
The INCL physics seem to be the best I’ve tried so far!
I appreciate all the help given to me in this thread!