|Message: Physics List for n+diamond reactions||Not Logged In (login)|
Click on the Forum title, e.g. on the "Forums by Category" page, to read a sequence of postings to the Forum and its threads all in one page. If you are only interested in one thread or the thread following a specific posting, click the thread or the posting, which takes you to a smaller page, which contains only the part you are interested in and may be easier to navigate.
Messages are "chained" if there are only replies at the first level, i.e. 1/1.html, 1/1/1.html etc. In case of "chained" messages the message number is replaced by the icon and there is no indentation.
Inline: Display the subject line only or also the text of the posting(s); for the choice "All" the "Outline" choices are switched off.
|1||0||1||no text / full text of posting|
|2||1||All||text for level 1 only / text for All postings|
Outline: Choose the depth of the posting thread, successive toggle controls provide increasing detail.
|1||2||1||2 levels / 1 level (original posting)|
|2||3||2||3 levels / 2 levels|
|3||3||All||3 levels / all levels (all postings)|
I have been running a simulations with neutrons of kinetic energy around 18 MeV hitting a diamond to compare experimental results and simulations, so that a better background extimation could be achieved. Having no idea of which physics list to use, I tried different ones, thought I am not really sure about the differences among them.
I obtain the best results using the file that I renamed PhysicsList1.cc. I created that copying and pasting from online and,as far as I understood, it uses the G4EmLivermorePhysics.hh for the elm part, G4HadronInelasticQBBC() for inelastic hadronic and G4HadronElasticPhysics() for elastic hadronic. The resulting energy spectrum with the (n,alpa), (n,p0),.. peaks is PhysicsList1.pdf.
I tried then using the high precision neutron physics, with the G4NDL (I.e. ENDF-BVII) cross sections. The results is that several peaks disappear, while they should not. The PhysicsList for this case is PhysicsList2.cc and the results PhysicsList2.pdf The PhysicsList2.cc, as far as I understood, still uses Livermore for elm and uses G4HadronPhysicsQGSP_BIC_HP for hadronic physics. Is it using this last both for elastic and inelastic hadronic reactions, pls?
The plot obtained from the PhysicsList1 is the most similar to the experimental results, but does not really matche them. My geometry models already includes materials close to the diamond, so I doubt the problem are the scattered neutrons.
I have been playing arond changing few things in the physics list, but since I do not really understan how to build them, I am not aware of their differences and I do not even know if I am including both elestic and inelastic hadron interactions in all cases or not.
Could anyone tell me which is the best physics list to use for neutrons + diamond reactions between 14 and 25 MeV, pls?
Thanks, regards, Franciuska
Attachment: http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2016/08/04/07.13-24578-PhysicsList1.cc http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2016/08/04/07.19-2191-PhysicsList1.pdf http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2016/08/04/07.25-95760-PhysicsList2.pdf http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2016/08/04/07.25-20179-PhysicsList2.cc
|Inline Depth:||Outline Depth:||Add message:|