[Gate-users] range for diferent elements in water

[Andrea González]

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! 		 	   		  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20150310/053f771b/attachment.html>