|Message: Re: Non energy conservation in a hadron elastic scattering between neutrons and protons||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)|
To show you that there is something strange I used the QGSP_BERT_HP reference physics list instead of what I used before. I saw exactly the same results.
I join you two plots made with the QGSP_BERT_HP reference physics list : PR_10MeV.pdf corresponds to the number of protons recoils vs their vertex kinetic energy for 10MeV neutrons created in a liquid scintillator and PR_100MeV.pdf is for 100MeV neutrons created in the same liquid scintillator.
You can see that if everything is working well for the 10MeV neutrons (10000 neutrons created), there is a huge peak at 100MeV for the 100MeV neutrons (10000 neutrons created). It corresponds to what I put in my first mail : the protons are taking in only one step all the neutrons energies and even more, which is physically impossible :
- neutron 2112 E: 0 keV; E0: 100000 keV; Vertex Position: (0,3230,0) TrackID: 1; ParentID: 0; Step: 1; HadronElastic - proton 2212 E: 100054 keV; E0: 100125 keV; Vertex Position: (0,3058.76,0) TrackID: 2; ParentID: 1; Step: 1; hIoni - proton 2212 E: 99797 keV; E0: 100125 keV; Vertex Position: (0,3058.76,0) TrackID: 2; ParentID: 1; Step: 2; hIoni - proton 2212 E: 99612.8 keV; E0: 100125 keV; Vertex Position: (0,3058.76,0) TrackID: 2; ParentID: 1; Step: 3; hIoni
Do you have any idea about what is happening ?
Attachment: http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2012/05/14/00.03-32114-PR_10MeV.pdf http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2012/05/14/00.03-7671-PR_100MeV.pdf