Message: Energy deposition via ion particle with RadioactiveDecay process Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Energy deposition via ion particle with RadioactiveDecay process 

Forum: Hits, Digitization and Pileup
Date: 23 Jan, 2012
From: Clare Huang <Clare Huang>

Hi all,

I am working on recording energy deposited in a water phantom.  I have defined sensitive detectors to process hits in all the water cells.  I got the program to compile and run.  However, I am trying to understand the results from the tracking versbose.  Here is the macro file that i used to generate the result below.

/grdm/allVolumes # use all volumes
/grdm/fBeta 0
/grdm/verbose 0
/grdm/analogueMC 1 # (analog monte carlo)
/gps/pos/type Volume
/gps/pos/shape Para
/gps/pos/centre 0 0 0
/gps/ang/type iso
/gps/pos/halfx 1.0 mm 
/gps/pos/halfy 1.0 mm
/gps/pos/halfz 1.0 mm
/gps/source/intensity 1.0   
/gps/particle ion
/gps/ion 53 131 
/grdm/nucleusLimits 131 131 53 53
/tracking/verbose 1
/run/beamOn 1 # number of particles

See below for the tracking verbose for a run with 1 I-131 ion with radioactivedecay process...

G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("World", copy 0) still exists and is being used.
  WARNING: This does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
Starting Run: 0
Starting Event: 0

* G4Track Information:   Particle = I131[0.0],   Track ID = 1,   Parent ID = 0

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0   -0.428   -0.267    0.886         1        0        0         0 PhantomSens initStep
    1   -0.428   -0.268    0.885         0        1  0.00147   0.00147 PhantomSens ionIoni
ix,iy,iz = 49,49,100, indx = 989900, Edep = 1
    2   -0.428   -0.268    0.885         0        0        0   0.00147 PhantomSens RadioactiveDecay

this is the tracking verbose from the very first primary particle.  Why is the kinetic energy (KinE) set at 1 MeV and the energy deposited (dE) is 1 MeV for the following step? It looks there 1 MeV is the default value set somewhere in a function.  Does anyone know what this means? In this regards, do I not count this track as part of the energy deposited in this run?

I attached the codes where the physicslist are defined.  The physics definition is similar to that used in examples/extended/runAndEvent/RE02.  But I will attach for clarify.  let me know if any of you have any thought on this question.  I can include more files if needed.

Thanks a lot in advance for all your time and help.

best regards,


 Add Message Add Message
to: "Energy deposition via ion particle with RadioactiveDecay process"

 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 ]