|Message: Re: Apparent Non-conservation of Energy||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)|
That did reduce the number of unphysical gammas, however, it also drastically reduced the number of gammas absorbed by my detectors (from around 100k to 7(E0)). This could be due to a mistake on my side and I am working to test that. I am still seeing some gammas that don't make sense. I don't currently have the information from Geant4's verbose tracking option, but I print out a message any time the total gamma energy in a singular nCapture event exceeds the allowed value for that material. Below is my output. As you can see, the events all occur in LiF, making me wonder if there could be something specifically wrong with that. From what I have seen about G4NEUTRONHP_DO_NOT_ADJUST_FINAL_STATE, it has caused some issues when the reaction has other products. The example given here: http://hypernews.slac.stanford.edu/HyperNews/geant4/get/hadronprocess/1249/1/2.html?inline=-1
uses n+35Cl -> S* + p -> S0 + gamma_excited + p. It also references n+6Li-->t+alpha, which occurs in my simulation. As Malachi Schram notes, "if you set the G4NEUTRONHP_DO_NOT_ADJUST_FINAL_STATE variable, you will have the proper excited state gamma (individual gamma) but the wrong total gamma energy for a given event." This is precisely the issue I started with, before I knew about this new environmental variable.
These are the total energies of gammas produced so far. As you can see, there were at least two 6.773 MeV gammas, along with other energies.
Unphysical LiF: 6.773 MeV Unphysical LiF: 13.546 MeV Unphysical LiF: 20.796 MeV