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

Chen, Yu Yu.Chen at umassmed.edu
Tue Dec 11 05:01:28 CET 2007


How did you get those two plots?
Using Singles or Coincidences?
 
It seems that photons you recorded are even not isotropical (peak in y=0) in transverse plane.
I do not understand what are those scattered points around central points in your left plot.
 
You better put your PET geometry and coincidence module macro too.
 
And please plot your sourcePosX,Y,Z in Singles.
   
 
Yu Chen, Ph.D.
University of Massachusetts Medical School
Division of Nuclear Medicine
55 Lake Avenue North
Worcester, MA 01655-0243
Phone:  (508) 856-6123
Fax:    (508) 856-4572

________________________________

From: gate-users-bounces at lists.healthgrid.org on behalf of Joon Young Suk
Sent: Mon 12/10/2007 6:55 PM
To: gate-users at lists.healthgrid.org
Subject: [Gate-users] isotropic emission from a rod source ?



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

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



-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20071210/94b1ded2/attachment.htm>


More information about the Gate-users mailing list