[Gate-users] isotropic emission from a rod source ?

Joon Young Suk joonyoung.suk at gmail.com
Tue Dec 11 00:55:01 CET 2007


I built two rod sources which set as isotropic sources, but  
simulation result show most data have been acquired from the  
horizontal plane.
Please check my result in the following link :

http://www.joonyoung.net/root_result_2007-12-10.jpg

Could you please find out what is my problem.
Here is a part of my macro ;


##### PHANTOM GEOMETRY ##############
## Cylindrical phantom has two rod sources ####

# H E A D
/gate/world/daughters/name phantomHead
/gate/world/daughters/insert cylinder
/gate/phantomHead/placement/setTranslation {phXposition} 0. 0. mm
/gate/phantomHead/setMaterial Air
/gate/phantomHead/geometry/setRmax 22.5 mm
/gate/phantomHead/geometry/setHeight 50. mm
/gate/phantomHead/vis/forceWireframe
/gate/phantomHead/vis/setColor blue
/gate/geometry/update

# Add Rods for source confinement

# ROD0 : rod at right-side (positive x-axis)
# Diameter : 1.0 mm
/gate/phantomHead/daughters/name rod0
/gate/phantomHead/daughters/insert cylinder
/gate/rod0/geometry/setRmax 0.5 mm
/gate/rod0/geometry/setRmin 0.0 mm
/gate/rod0/geometry/setHeight 50. mm
/gate/rod0/placement/setTranslation {rod0Xposition} 0. 0. mm
/gate/rod0/setMaterial Water
/gate/rod0/vis/setColor red

# ROD1 : rod at left-side(negative x-axis)
/gate/phantomHead/daughters/name rod1
/gate/phantomHead/daughters/insert cylinder
/gate/rod1/geometry/setRmax 0.5 mm
/gate/rod1/geometry/setRmin 0.0 mm
/gate/rod1/geometry/setHeight 50. mm
/gate/rod1/placement/setTranslation {rod1Xposition} 0. 0. mm
/gate/rod1/setMaterial Water
/gate/rod1/vis/setColor red

# PHANTOM defined as SENSITIVE
/gate/phantomHead/attachPhantomSD


###################################################################
# S O U R C E S
#### two rod sources
## rodsource0 -> confines to rod0 in the phantom
## rodsource1 -> confiines  to rod1 in the phantom
###################################################################
/gate/source/addSource rodsource0
/gate/source/rodsource0/setActivity {activity} becquerel
/gate/source/rodsource0/setType backtoback
/gate/source/rodsource0/gps/particle gamma
/gate/source/rodsource0/gps/energytype Mono
/gate/source/rodsource0/gps/monoenergy 0.511 MeV
/gate/source/rodsource0/gps/type Volume
/gate/source/rodsource0/gps/shape Cylinder
/gate/source/rodsource0/gps/centre {r0sourceXposition} 0. 0. mm
/gate/source/rodsource0/gps/radius 0.5 mm
/gate/source/rodsource0/gps/halfz 25. mm
/gate/source/rodsource0/gps/angtype iso
/gate/source/rodsource0/gps/confine rod0_P

/gate/source/addSource rodsource1
/gate/source/rodsource1/setActivity {activity} becquerel
/gate/source/rodsource1/setType backtoback
/gate/source/rodsource1/gps/particle gamma
/gate/source/rodsource1/gps/energytype Mono
/gate/source/rodsource1/gps/monoenergy 0.511 MeV
/gate/source/rodsource1/gps/type Volume
/gate/source/rodsource1/gps/shape Cylinder
/gate/source/rodsource1/gps/centre {r1sourceXposition} 0. 0. mm
/gate/source/rodsource1/gps/radius 0.5 mm
/gate/source/rodsource1/gps/halfz 25. mm
/gate/source/todsource1/gps/angtype iso
/gate/source/rodsource1/gps/confine rod1_P

### ROOT ####
/gate/output/root/setSaveRndmFlag 0
/gate/output/root/setFileName {outFile}
/gate/output/root/setRootHitFlag 0
/gate/output/root/setRootSinglesFlag 1
/gate/output/root/setRootCoincidencesFlag 1
/gate/output/root/setRootNtupleFlag 0

# D A T A  A C Q U I S I T I O N
/gate/application/setTimeSlice {timeSlice} s
/gate/application/setTimeStart 0. s
/gate/application/setTimeStop {timeStop} s

###### START DAQ #######
/gate/application/startDAQ




Thank you




More information about the Gate-users mailing list