|Message: energy conservation fails in TestEm5 using Geant4-DNA physics list.||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)|
Dear Geant4 experts/users
The TestEm5 example comes with a stacking action class which is able to print the energy of all secondaries created. When I process the txt file using a simple Python code, I realize that the sum of the energy of all secondaries is higher than the input energy in the system; nevertheless, the energy deposited in the phantom is exactly the same than the input energy. Please see details of the simulation in the file attached. Besides, the TestEm5 has integrated a kill secondaries command line, when I use this feature, the sum of all the secondaries created is less than the input energy, as expected and the energy deposited in the phantom is exactly the input energy. I also have plotted a graph that compares the energies of the secondaries at creation turning on/off the kill secondaries command line (see attached), it seems that Geant4 is creating a higher number of very low energy electrons that it should so it breaks the energy conservation.
Simulation set-up: 10 MeV protons irradiating a water phantom of 2 mm thickness, i.e all protons are absorbed in the phantom. Geant4-DNA physics list is active.
The output of the Geant4 after processing shows the following ( 10 MeV protons of 10 MeV, i.e 100 MeV total energy input in the system):
no kill active: energy for e- (MeV)=', 120.3801 energy for protons (MeV)=', 633.3470 energy for Hydrogen (MeV)=', 633.9289 energy for gamma (MeV)=', 0.0167, total energy=', 1387.6728
kill secs On: energy for e- (MeV)=', 53.5953 energy for protons (MeV)=', 0 energy for Hydrogen (MeV)=', 32.32403 energy for gamma (MeV)=', 0.012, total energy (MeV)=', 85.9313
Attachment: http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2019/02/12/04.24-14392-killVsNoKill.png http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2019/02/12/04.25-40551-Mydna.mac