[Gate-users] Radioactive Decay problem

[daniel at if.usp.br]

Hello everybody,

I am trying to simulate the dose deposited by a I-131 source in a  
phantom. I am using GPS and RadioactiveDecay proccess. Here is an very  
simple example (without the phantom) that I can use as reference to  
describe my problem:
/vis/disable
/gate/verbose Error     2
/gate/verbose Warning   2
/run/verbose 1
/event/verbose 1
/tracking/verbose 0
/gate/geometry/setMaterialDatabase GateMaterials.db
/gate/world/setMaterial                            Water
/gate/world/geometry/setXLength                    3 cm
/gate/world/geometry/setYLength                    3 cm
/gate/world/geometry/setZLength                    3 cm
/gate/physics/addProcess PhotoElectric
/gate/physics/processes/PhotoElectric/setModel StandardModel
/gate/physics/addProcess Compton
/gate/physics/processes/Compton/setModel StandardModel
/gate/physics/addProcess GammaConversion
/gate/physics/processes/GammaConversion/setModel StandardModel
/gate/physics/addProcess ElectronIonisation
/gate/physics/processes/ElectronIonisation/setModel StandardModel e-
/gate/physics/processes/ElectronIonisation/setModel StandardModel e+
/gate/physics/processes/ElectronIonisation/setStepFunction e+ 0.2 0.1 mm
/gate/physics/processes/ElectronIonisation/setStepFunction e- 0.2 0.1 mm
/gate/physics/addProcess Bremsstrahlung
/gate/physics/processes/Bremsstrahlung/setModel StandardModel e-
/gate/physics/processes/Bremsstrahlung/setModel StandardModel e+
/gate/physics/addProcess PositronAnnihilation
/gate/physics/addProcess MultipleScattering
/gate/physics/processes/MultipleScattering/setGeometricalStepLimiterType e-  
distanceToBoundary
/gate/physics/processes/MultipleScattering/setGeometricalStepLimiterType e+  
distanceToBoundary
/gate/physics/addProcess RadioactiveDecay
/gate/physics/Gamma/SetCutInRegion      world 1 mm
/gate/physics/Electron/SetCutInRegion   world 1 mm
/gate/physics/Positron/SetCutInRegion   world 1 mm
/gate/source/addSource iodine131
/gate/source/iodine131/setActivity 10. becquerel
/gate/source/iodine131/gps/particle ion
/gate/source/iodine131/gps/ion 53 131 0 0
/gate/source/iodine131/setForcedUnstableFlag true
/gate/source/iodine131/setForcedHalfLife 692902.08 s
/gate/source/iodine131/gps/energytype Mono
/gate/source/iodine131/gps/monoenergy 0. MeV
/gate/source/iodine131/gps/angtype iso
/gate/source/iodine131/gps/number 1
/gate/source/iodine131/gps/centre 0. 0. 0. cm
/gate/source/iodine131/gps/type Volume
/gate/source/iodine131/gps/shape Sphere
/gate/source/iodine131/gps/radius 1. cm
/gate/run/initialize
/gate/EnableGeneralParticleSource
/run/beamOn 10
#/gate/application/noGlobalOutput
#/gate/application/setTotalNumberOfPrimaries 10
#/gate/application/start
#/gate/application/describe

And I am getting the following message for each event:
[G4] =====================================
[G4]   G4EventManager::ProcessOneEvent()
[G4] =====================================
[G4] 0 primaries are passed from G4EventTransformer.
[G4] !!!!!!! Now start processing an event !!!!!!!
[G4] NULL returned from G4StackManager.
[G4] Terminate current event processing.

It looks like the Radioactive Decay proccess is not generating the  
particles. Am I missing something in that script?

Best regards,
Daniel

----------------------------------------------------------------
This message was sent using IMP, the Internet Messaging Program.