|Message: Problems with ParticleHP for proton inelastic||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 am using Geant4 10.02 with the vGate interface. I downloaded the TENDL library to study the radioisotope production and neutron emission from the O18+proton reaction. I am using the QGSP_BIC_AllHP physics list.
To understand how Geant4 was managing the proton interactions I looked into the source code. I was surprised by what I found: if i am not mistaken, the inelastic proton interactions are not treated as outgoing channels with multiple particles, but directly sorting possible outgoing particles independantly. To make it clearer with an example: in place of sorting randomly (according to cross sections) for example between the outgoing n+F18 and n+alpha+N15 channels, the code seems to sort randomly if a n will be produced, then if an alpha will be produced, and so on! Maybe this philosophy make sense in high energy physics, but this seems strange at low energy (< 50 MeV).
I think there also is a problem with the energy of outgoing particles in the ParticleHP: the Qvalue is not taken into account ! When testing by forcing the reaction O18(p,n+alpha)N15, I got neutrons up to the kinetic energy of the protons, which should be impossible because the (p,n+alpha) reaction is endothermic (QV=-6.9 MeV). This problem could come from the fact that ParticleHP was derived from NeutronHP for which only positive Qvalues are expected.
Hopefully someone will be able to prove me wrong (browsing through the Geant4 code is no easy task) or maybe corrections have been implemented in more recent releases (I am using the latest release available with vGate).