[Gate-users] Re : issues with the ion source

[Samuel BURG]

Of course ... with a 0keV energy your particle will not go far away ...

 

Samuel Burg


CHU Bichat Claude Bernard
75018 Paris



________________________________
De : "Marta.Lai at ca.infn.it" <Marta.Lai at ca.infn.it>
À : gate-users at lists.opengatecollaboration.org
Envoyé le : Jeu 21 octobre 2010, 10h 40min 39s
Objet : [Gate-users] issues with the ion source

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


_______________________________________________
Gate-users mailing list
Gate-users at lists.opengatecollaboration.org
http://lists.opengatecollaboration.org/mailman/listinfo/gate-users