Message: Re: Energy Conservation Problem Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: Energy Conservation Problem 

Forum: Event and Track Management
Re: Question Energy Conservation Problem (Gianfranco Gargano)
Re: None Re: Energy Conservation Problem (michel maire)
Date: 22 Mar, 2011
From: Loc Grevillot <Loc Grevillot>

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:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 None: Re: Energy Conservation Problem   (Gianfranco Gargano - 22 Mar, 2011)
1 None: Re: Energy Conservation Problem   (Gianfranco Gargano - 22 Mar, 2011)
2 None: Re: Energy Conservation Problem   (Loc Grevillot - 22 Mar, 2011)
1 None: Re: Energy Conservation Problem   (Gianfranco Gargano - 22 Mar, 2011)
2 None: Re: Energy Conservation Problem   (michel maire - 22 Mar, 2011)
1 Warning: Re: Energy Conservation Problem   (Vladimir Ivanchenko - 22 Mar, 2011)
 Add Message Add Message
to: "Re: Energy Conservation Problem"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]