[Gate-users] range for diferent elements in water

[J.J.M. Kouwenberg]

Dear Andrea,

I would suggest adding a complete set of physical processes:

gate/physics/addProcess ElectronIonisation
/gate/physics/processes/ElectronIonisation/setModel StandardModel e-
/gate/physics/processes/ElectronIonisation/setModel StandardModel e+

/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 eMultipleScattering e- 
/gate/physics/addProcess eMultipleScattering e+
/gate/physics/processes/eMultipleScattering/setModel Urban95Model

Also make sure that the order of your code is correct:
 - Geometry and phantoms
 - Physics
 - Initialization
 - Digitizers
 - Source
 - Output 
 - Run

With kind regards,

Jasper J.M. Kouwenberg
Medical Physics Graduate student
Nuclear Medicine

E-mail j.kouwenberg at erasmusmc.nl | Phone +316 2532 9809

________________________________________
From: gate-users-bounces at lists.opengatecollaboration.org [gate-users-bounces at lists.opengatecollaboration.org] on behalf of Andrea González [andreita_5_92 at hotmail.com]
Sent: 10 March 2015 11:25
To: gate-users at lists.opengatecollaboration.org
Subject: Re: [Gate-users] range for diferent elements in water

Hi all, I want to simulate positron annihilation in water, and compare the range for diferent elements.
My problem is that I always get the same range regardless of the element I use.
Please, can you help me?
My code:


#=====================================================
# VISUALISATION and VERBOSE
#=====================================================


/control/execute                    ./carpeta/Visualisation.mac
/control/execute                    ./carpeta/Verbose.mac

#=====================================================
# GEOMETRY
#=====================================================

/gate/geometry/setMaterialDatabase    ./GateMaterials.db


#=====================================================
#     W O R L D
#=====================================================

/gate/world/geometry/setXLength       30 cm
/gate/world/geometry/setYLength       30 cm
/gate/world/geometry/setZLength       30 cm
/gate/world/setMaterial               Air

#=======================================================
#                   G E O M E T R Y
#=======================================================

/control/execute                    ./carpeta/geometria.mac

#=====================================================
#                  P H Y S I C S
#=====================================================

/gate/physics/addProcess RadioactiveDecay
/gate/physics/addProcess PositronAnnihilation e+
#/gate/physics/addProcess G4PositronAnnihilation e+
#/gate/physics/processes/G4PositronAnnihilation/setModel StandardModel

/gate/physics/processList Enabled
/gate/physics/processList Initialized


#=====================================================
#         INITIALISATION
#=====================================================

/gate/run/initialize

#=====================================================
#  D A T A   O U T P U T
#=====================================================


/gate/output/root/enable
/gate/output/root/setFileName                     PuntualO
/gate/output/root/setRootHitFlag          1
/gate/output/root/setRootSinglesFlag      1
/gate/output/root/setRootOpticalFlag      1



#/control/execute                        ./archivo/Projection.mac

#=====================================================
#   MEASUREMENT SETTINGS
#=====================================================

# A C Q U I S I T I O N

/control/execute                        ./carpeta/StartDAQ.mac


#=======================================================
#                    Geometry
#=======================================================


/gate/world/daughters/name                      scanner
/gate/world/daughters/insert                    cylinder

#Fijamos las dimensiones del sistema.
#Seleccionamos el material.

/gate/scanner/geometry/setRmax        12. mm
/gate/scanner/geometry/setHeight      6. mm
/gate/scanner/setMaterial             Water

/gate/scanner/attachCrystalSD

#=======================================================
#                    Source
#=======================================================

/control/execute ./carpeta/fuenteO.mac

/gate/source/addSource Oxigen
/gate/source/Oxigen/setActivity 11000 becquerel

/gate/source/Oxigen/gps/particle e+
/gate/source/Oxigen/setForcedUnstableFlag true
/gate/source/Oxigen/setForcedHalfLife 122.24 s
/gate/source/Oxigen/gps/energytype Oxygen15


/gate/source/Oxigen/gps/type Plane
/gate/source/Oxigen/gps/shape Circle
/gate/source/Oxigen/gps/radius 0.5 mm
/gate/source/Oxigen/gps/angtype iso
/gate/source/Oxigen/gps/centre 0. 0. 0. cm



/gate/source/list




Thank you!