[Gate-users] counting rates again

j.hartmann at dkfz-heidelberg.de j.hartmann at dkfz-heidelberg.de
Fri Mar 2 18:17:18 CET 2007


Hello again!
You are right, I should have sent the details straight away, here the macro of
the full simulation:

#
#	V I S U A L I S A T I O N
#
/vis/open OGLSX
/vis/viewer/reset
/vis/drawView 0 90
#/vis/viewer/zoom 4
/vis/viewer/set/style surface
/vis/drawVolume
/tracking/storeTrajectory 1
/vis/scene/endOfEventAction accumulate
/vis/viewer/update
/gate/geometry/enableAutoUpdate

#
#     W O R L D
#
/gate/world/geometry/setXLength 81. cm
/gate/world/geometry/setYLength 81. cm
/gate/world/geometry/setZLength 81. cm


#     D E F I N I T I O N  O F  T H E  S Y S T E M
#
#     C Y L I N D R I C A L 

/gate/world/daughters/name cylindricalPET
/gate/world/daughters/insert cylinder
/gate/cylindricalPET/setMaterial Vacuum
/gate/cylindricalPET/geometry/setRmax 40 cm
/gate/cylindricalPET/geometry/setRmin 33 cm
/gate/cylindricalPET/geometry/setHeight 22 cm
/gate/cylindricalPET/vis/forceWireframe

#
#     H E A D
#
/gate/cylindricalPET/daughters/name rsector
/gate/cylindricalPET/daughters/insert box
/gate/rsector/placement/setTranslation 0 35.5  0 cm
/gate/rsector/geometry/setXLength 21.12 cm
/gate/rsector/geometry/setYLength 1. cm
/gate/rsector/geometry/setZLength 0.33 cm
/gate/rsector/setMaterial TiO2
/gate/rsector/vis/forceWireframe

#
#	C R Y S T A L
#
/gate/rsector/daughters/name crystal
/gate/rsector/daughters/insert box
/gate/crystal/geometry/setXLength 3.1 mm
/gate/crystal/geometry/setYLength 10. mm
/gate/crystal/geometry/setZLength 3.1 mm
/gate/crystal/setMaterial LSO
/gate/crystal/vis/forceWireframe


#	R E P E A T    C R Y S T A L  
/gate/crystal/repeaters/insert cubicArray
/gate/crystal/cubicArray/setRepeatNumberX 64
/gate/crystal/cubicArray/setRepeatNumberY 1
/gate/crystal/cubicArray/setRepeatNumberZ 1
/gate/crystal/cubicArray/setRepeatVector 3.3 0 0 mm

#	R E P E A T   H E A D
/gate/rsector/repeaters/insert ring
/gate/rsector/ring/setRepeatNumber 2


#	A T T A C H    S Y S T E M
/gate/systems/cylindricalPET/rsector/attach rsector
/gate/systems/cylindricalPET/module/attach crystal

#
#	M O V E M E N T S   O F   T H E   C A M E R A
#
/gate/cylindricalPET/moves/insert orbiting
/gate/cylindricalPET/orbiting/setSpeed 0.2 deg/s
/gate/cylindricalPET/orbiting/setPoint1 0 0 0 cm
/gate/cylindricalPET/orbiting/setPoint2 0 0 1 cm


#     D E F I N I T I O N  O F  T H E   P H Y S I C S                       #

#       EM  P R O C E S S
/gate/physics/gamma/selectRayleigh lowenergy
/gate/physics/gamma/selectPhotoelectric standard
/gate/physics/gamma/selectCompton standard

#	I N A C T I V E    S E C O N D A R Y   E L E C T R O N S 
/gate/physics/setElectronCut 1. m

#	I N A C T I V E    X - R A Y S
/gate/physics/setXRayCut 1. GeV
/gate/physics/setDeltaRayCut 1. GeV


#     I N I T I A L I Z A T I O N  

/gate/geometry/enableAutoUpdate
/run/initialize

#                                                                 #
#  D E F I N I T I O N   O F   T H E   A C Q U I S I T I O N    #
#  D I G I T I Z E R   &   C O I N C I D E N C E   S O R T E R  #
#                                                                 #

#	A D D E R
/gate/digitizer/Singles/insert adder

#       R E A D O U T 
/gate/digitizer/Singles/insert readout
/gate/digitizer/Singles/readout/setDepth 2

#      E N E R G Y   B L U R R I N G
/gate/digitizer/Singles/insert blurring
/gate/digitizer/Singles/blurring/setResolution 0.26
/gate/digitizer/Singles/blurring/setEnergyOfReference 511. keV

#	C O I N C I D E N C E    S O R T E R 
/gate/digitizer/Coincidences/minSectorDifference 1
/gate/digitizer/Coincidences/setWindow 10. ns
/gate/digitizer/Coincidences/setDepth 2


#        D E F I N I T I O N   O F                     #
#                T H E   O U T P U T   F I L E         #
#                                                       #

#       C H A N G E    T H E     S E E D (1)  O R   N O T (0)
/gate/output/root/setSaveRndmFlag 1

#       S E T U P  - R O O T _ F I L E
/gate/output/root/setFileName Your_PET_Simu
/gate/output/root/setRootSinglesAdderFlag 1
/gate/output/root/setRootSinglesReadoutFlag 1
/gate/output/root/setRootHitFlag 1
/gate/output/root/setRootSinglesFlag 1
/gate/output/root/setRootNtupleFlag 1
/gate/output/root/setRootCoincidencesFlag 1

#       S E T U P  - A S C I I _ F I L E
/gate/output/ascii/setOutFileSinglesAdderFlag 1
/gate/output/ascii/setOutFileSinglesReadoutFlag 1
#/gate/output/ascii/setOutFileSinglesThresholderFlag 0
#/gate/output/ascii/setOutFileSinglesUpholderFlag 0
/gate/output/ascii/setOutFileSinglesFlag 0
/gate/output/ascii/setOutFileHitsFlag 0
/gate/output/ascii/setOutFileCoincidencesFlag 1
/gate/output/ascii/setOutFileSizeLimit 200000000
/gate/output/ascii/setCoincidenceMask 1 0 0 0 0 0 0 1 0 0 0 1 0 0 0 0 0 0

#                                                                 #
#      D E F I N I T I O N   O F  V E R B O S I T Y   L E V E L   #
#                                                                 #
/control/verbose 0
/grdm/verbose 0
/run/verbose 0
/event/verbose 0
/tracking/verbose 0
/gate/application/verbose 0
/gate/generator/verbose 0
/gate/stacking/verbose 0
/gate/event/verbose 0
/gate/source/verbose 0

#                                                                 #
#          D E F I N I T I O N   O F  T H E   S O U R C E         #

##### simple particle source twogamma backtoback

/gate/source/addSource 2gamma1
/gate/source/2gamma1/gps/particle gamma
/gate/source/2gamma1/setType backtoback
/gate/source/2gamma1/setActivity 5000. becquerel
/gate/source/2gamma1/gps/energytype Mono
/gate/source/2gamma1/gps/monoenergy 511 keV
/gate/source/2gamma1/gps/angtype iso
/gate/source/2gamma1/gps/mintheta 90 deg
/gate/source/2gamma1/gps/maxtheta 90. deg
/gate/source/2gamma1/gps/minphi 0. deg
/gate/source/2gamma1/gps/maxphi 360. deg
/gate/source/2gamma1/gps/type Plane
/gate/source/2gamma1/gps/shape Circle
/gate/source/2gamma1/gps/radius 1 cm
#### placement of the source
/gate/source/2gamma1/gps/centre 0 0 0 cm
/gate/source/2gamma1/gps/confine NULL

#                                                                      #
#    S T A R T   A C Q U I S I T I O N  T I M E   P A R A M E T E R S  #
#                                                                      #
/gate/application/setTimeSlice     5.  s
/gate/application/setTimeStart     0.   s
/gate/application/setTimeStop      900.  s

/gate/application/startDAQ




More information about the Gate-users mailing list