[Gate-users] issues with the ion source

[Marta.Lai at ca.infn.it]

I'm trying to run a PET simulation with an ion source of Sodium22 with
Gate v6.0.
My source is defined with the command lines:



/gate/source/addSource                                         posiNa22

/gate/source/posiNa22/gps/particle                              ion
/gate/source/posiNa22/gps/ion					11 22 0 0
/gate/source/posiNa22/gps/monoenergy				0. keV
/gate/source/posiNa22/setActivity                               635000.
becquerel
/gate/source/posiNa22/gps/type					Volume
/gate/source/posiNa22/gps/shape					Sphere
/gate/source/posiNa22/gps/radius				0.5 mm
/gate/source/posiNa22/gps/centre				0.0 0.0 0.0 mm
/gate/source/posiNa22/setForcedUnstableFlag                     true
/gate/source/posiNa22/setForcedHalfLife                         984680064. s
/gate/source/posiNa22/gps/angtype                               iso
/gate/source/posiNa22/gps/monoenergy				0. keV

/gate/source/list

when I run the simulation I don't obtain any output (there is the output
.root but it's empty), that is the detector doesn't see any gamma ray from
the decay of the source.
I've tried a twogamma source and it works and it also works with a
positron source. Maybe there's some troubles with the physics list:

/gate/physics/addProcess PhotoElectric
/gate/physics/addProcess Compton
/gate/physics/addProcess GammaConversion
/gate/physics/addProcess LowEnergyRayleighScattering
/gate/physics/addProcess ElectronIonisation
/gate/physics/addProcess Bremsstrahlung
/gate/physics/addProcess PositronAnnihilationStd
/gate/physics/addProcess MultipleScattering e+
/gate/physics/addProcess MultipleScattering e-
/gate/physics/processList Enabled
/gate/physics/processList Initialized

or the cuts:

gate/physics/Gamma/SetCutInRegion      COO 0.1 keV
/gate/physics/Electron/SetCutInRegion   COO 1.0 m
/gate/physics/Gamma/SetCutInRegion      CEE 0.1 keV
/gate/physics/Electron/SetCutInRegion   CEE 1.0 m
/gate/physics/Gamma/SetCutInRegion      COOO 0.1 keV
/gate/physics/Electron/SetCutInRegion   COOO 1.0 m
/gate/physics/Gamma/SetCutInRegion      CEEE 0.1 keV
/gate/physics/Electron/SetCutInRegion   CEEE 1.0 m
/gate/physics/Gamma/SetCutInRegion      phantom 0.1 keV
/gate/physics/Electron/SetCutInRegion   phantom 1.0 m

Can anybody help me?
thanks, Marta