|Message: 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)|
I am facing non energy conservation in my MC and I need your help to understand what is happening.
Here is the thing : I am creating neutrons in a tank filled with a liquid scintillator. I am using a simple gun : /generator/gun neutron 0 3230 0 0 -1 0 100 0 0 0. I am looking for the protons recoils due to the hadron elastic scattering between neutrons and protons. If everything is working well below ~20MeV I meet strange things above :
- 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
I am using the following HadronElastic Models : - NeutronHPElastic : Emin[GeV]=0; Emax[GeV]=0.02; - G4LElastic : Emin[GeV]=0.02; Emax[GeV]=1.79769e+305;
I don't understand how the neutron can be directly at rest after a scatter with a proton and how this proton can have an energy higher than what was available at the beginning of the process. About 4.4% of all my protons recoils are in the situation showed above.
Is it a known problem ? How can I deal with it ?
|Inline Depth:||Outline Depth:||Add message:|