|Message: Re: Energy Conservation Problem||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)|
On Tue, 22 Mar 2011 11:53:05 GMT, michel maire wrote:
> On Mon, 21 Mar 2011 18:50:07 GMT, Gianfranco Gargano wrote: > > > Hy, I'm trying to simulate a 150 MeV proton in a water cube. There's a > > strange effect: the sum of all energies deposed in the water cube is not > > 150 MeV. I have considered the energy which escape from the water cube > > (neutrons or gammas), but anyway seems that in the last proton > > interaction near 15 MeV disappear. I'm using an energy cut value of > > 0.01mm and no step max limitation. I've used the physics list of the > > ProtonTherapy example. > > > > I invite you to run TestEm7 (in extended/electromagnetic) with attached macro. > You will get following output and attached plots : > > The run consists of 10000 proton of 150 MeV through 20 cm of G4_WATER > > Projected Range= 15.8575 cm rms= 1.74944 mm > Mean number of primary steps = 66.7172 > Total energy deposit= 150 MeV > niel energy deposit = 0 eV > > cheers, michel > > Attachment: > http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2011/03/22/04.48-47106-gargano.txt > http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2011/03/22/04.49-33789-edep.gif > http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2011/03/22/04.49-69819-range.gif >
It is true that using this example, 150 MeV are deposited, but this is because there is no inelastic hadronic interaction simulated. By adding such processes in the simulation (/testem/phys/addPhysics binary), then there is some missing energy.
The run consists of 10000 proton of 150 MeV through 40 cm of G4_WATER (density: 1 g/cm3 )
Projected Range= 14.5021 cm rms= 3.54837 cm Mean number of primary steps = 25.0488 Total energy deposit= 144.758 MeV niel energy deposit = 413.432 eV
|Inline Depth:||Outline Depth:||Add message:|