[gate-users] Energy Spectrum Woes - Dual Head and 4 Head Designs & Number of Heads

Angela M K Foudray afoudray at stanford.edu
Thu Sep 2 21:58:39 CEST 2004


I forgot to mention why I used the subject I did for the last post - the energy spectrum problems at the 
high energy (700 keV) that I mentioned is most apparent in the dual head design, goes slightly away in 
the 4 head design, and is not apparent at all in a ring design of 20+ "heads".  In other words, the severity 
of this anomalous high energy peak seems to be a function of the number of heads in the system.  (In case 
you're wondering, we have /gate/digitizer/Coincidences/minSectorDifference 1)

I'm attaching spectra (lower counts than the previous spectra I sent) for 4,5,6 heads and a ring design.
I'll also send you the .mac file for the dual head ("dhci_sen10cm.mac").

You'll notice, even for the ring, the spectra are not similar in the low energy region ("sasens_ring_np60_oe.mac")...


Ange
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 4head_nrg_spec.tif
Type: image/tiff
Size: 3243840 bytes
Desc: not available
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20040902/4f445059/attachment.tif>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 5head_nrg_spec.tif
Type: image/tiff
Size: 3243840 bytes
Desc: not available
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20040902/4f445059/attachment-0001.tif>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 6head_nrg_spec.tif
Type: image/tiff
Size: 3243840 bytes
Desc: not available
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20040902/4f445059/attachment-0002.tif>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ring_nrg_spec.tif
Type: image/tiff
Size: 3243840 bytes
Desc: not available
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20040902/4f445059/attachment-0003.tif>
-------------- next part --------------

# VISUALISATION
#/control/execute visu.mac

# OpenGL online visualization

/vis/open OGLSX
/vis/viewer/reset
/vis/viewer/viewpointThetaPhi 60 60
/vis/viewer/zoom 2.5
/vis/viewer/set/style surface
/vis/drawVolume
/tracking/storeTrajectory 1
/vis/scene/endOfEventAction accumulate
/vis/viewer/update
/vis/verbose 0
/gate/geometry/enableAutoUpdate


# World
# Define the world dimensions
/gate/world/geometry/setXLength 40 cm
/gate/world/geometry/setYLength 40 cm
/gate/world/geometry/setZLength 40 cm

#/gate/world/daughters/insert 3axes

#	cylindricalPET
/gate/world/daughters/name cylindricalPET
/gate/world/daughters/insert cylinder
/gate/cylindricalPET/placement/setTranslation 0.0 0.0 0.0 cm
/gate/cylindricalPET/geometry/setRmax 16 cm
/gate/cylindricalPET/geometry/setRmin 0 cm
/gate/cylindricalPET/geometry/setHeight 16 cm
/gate/cylindricalPET/setMaterial Air
/gate/cylindricalPET/setVisible 0
/gate/cylindricalPET/vis/forceWireframe
/gate/cylindricalPET/vis/setColor blue

#	rsector
/gate/cylindricalPET/daughters/name rsector
/gate/cylindricalPET/daughters/insert box
/gate/rsector/placement/setTranslation 6.1 0.0 0.0 cm
/gate/rsector/geometry/setXLength  2.2 cm
/gate/rsector/geometry/setYLength 15.4 cm
/gate/rsector/geometry/setZLength 11.05 cm
/gate/rsector/setMaterial Air
/gate/rsector/vis/forceWireframe
/gate/rsector/vis/setVisible 1
/gate/rsector/vis/setColor red

#       MODULE
/gate/rsector/daughters/name module
/gate/rsector/daughters/insert box
/gate/module/placement/setTranslation 0. 0.0 0.0 cm
/gate/module/geometry/setXLength 9 mm
/gate/module/geometry/setYLength 8 mm
/gate/module/geometry/setZLength 1.3 mm
/gate/module/setMaterial Air
/gate/module/vis/setVisible 0


#	crystal
/gate/module/daughters/name crystal
/gate/module/daughters/insert box
/gate/crystal/placement/setTranslation 0. 0.0 0.0 cm
/gate/crystal/geometry/setXLength 3. mm
/gate/crystal/geometry/setYLength 1. mm
/gate/crystal/geometry/setZLength 1.3 mm
/gate/crystal/setMaterial Air
/gate/crystal/vis/setVisible 0
/gate/crystal/vis/forceWireframe
/gate/crystal/vis/setColor yellow 

#	LSO C R Y S T A L  
/gate/crystal/daughters/name LSO
/gate/crystal/daughters/insert box
/gate/crystal/placement/setTranslation 0.0 0.0 0.0 cm
/gate/LSO/geometry/setXLength 3. mm
/gate/LSO/geometry/setYLength 1. mm
/gate/LSO/geometry/setZLength 1. mm
/gate/LSO/setMaterial LSO
/gate/LSO/vis/setVisible 1
/gate/LSO/vis/setColor green


#	R E P E A T  Crystal   
/gate/crystal/repeaters/insert cubicArray
/gate/crystal/cubicArray/setRepeatNumberX 3
/gate/crystal/cubicArray/setRepeatNumberY 8
/gate/crystal/cubicArray/setRepeatNumberZ 1
/gate/crystal/cubicArray/setRepeatVector 3.0 1. 0.0 mm


#       R E P E A T Module
/gate/module/repeaters/insert cubicArray
/gate/module/cubicArray/setRepeatNumberX 2
/gate/module/cubicArray/setRepeatNumberY 14
/gate/module/cubicArray/setRepeatNumberZ 85
/gate/module/cubicArray/setRepeatVector 11.0 11. 1.3 mm


#	R E P E A T rsector
/gate/rsector/repeaters/insert ring
/gate/rsector/ring/setRepeatNumber 2
#/gate/rsector/ring/setAngularPitch 180. deg

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

#	A T T A C H    C R Y S T A L  SD

/gate/LSO/attachCrystalSD

#	D I G I T I Z E R 
/gate/digitizer/Singles/insert adder

/gate/digitizer/Singles/insert readout
/gate/digitizer/Singles/readout/setDepth 1

/gate/digitizer/Singles/insert blurring
/gate/digitizer/Singles/blurring/setResolution 0.12
/gate/digitizer/Singles/blurring/setEnergyOfReference 511.0 keV

/gate/digitizer/Singles/insert thresholder
/gate/digitizer/Singles/thresholder/setThreshold 150 keV
/gate/digitizer/Singles/insert upholder
/gate/digitizer/Singles/upholder/setUphold 850 keV

#####Temporal resolution#####
/gate/digitizer/Singles/insert timeResolution
/gate/digitizer/Singles/timeResolution/setTimeResolution 2 ns

################
# DEADTIME
/gate/digitizer/Singles/insert deadtime
/gate/digitizer/Singles/deadtime/setDeadTime 5000. ns
/gate/digitizer/Singles/deadtime/setMode paralysable
/gate/digitizer/Singles/deadtime/chooseDTVolume module


#	C O I N C I    S O R T E R 
/gate/digitizer/Coincidences/minSectorDifference 1
/gate/digitizer/Coincidences/setWindow 20. ns

#/gate/systems/cylindricalPET/describe

#######################################################
##   		P H A N T O M 
#######################################################

#/gate/world/daughters/name sphere
#/gate/world/daughters/insert sphere
#/gate/sphere/setMaterial Water
#/gate/sphere/vis/setColor grey
#/gate/sphere/geometry/setRmax   2. mm
#/gate/sphere/geometry/setRmin   0.0 cm
#/gate/sphere/attachPhantomSD
#/gate/sphere/placement/setTranslation 0.0 0.0 0.0 cm


########################################
##        Physics.mac
########################################
#       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 30. cm

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

/gate/physics/verbose 2


#	I N I T I A L I Z E 
/gate/systems/cylindricalPET/verbose 0
/gate/geometry/enableAutoUpdate
/run/initialize


#	V E R B O S I T Y
/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

#/control/execute source_single.mac
####################################################
#Back to back point source 
####################################################
# 1 mCi = 37000000 becquerel

/gate/source/addSource na22
/gate/source/na22/setActivity 0.0001 Ci
/gate/source/na22/setType backtoback
/gate/source/na22/gps/particle gamma
/gate/source/na22/gps/energytype Mono
/gate/source/na22/gps/monoenergy 511. keV
/gate/source/na22/gps/type Point
/gate/source/na22/gps/confine NULL
/gate/source/na22/gps/angtype iso
/gate/source/na22/gps/mintheta 0. deg
/gate/source/na22/gps/maxtheta 180. deg
/gate/source/na22/gps/minphi 0. deg
/gate/source/na22/gps/maxphi 360. deg
/gate/source/na22/gps/centre 0. 0. 0. cm

####################################################
#Single F18 point source
####################################################

#/gate/source/addSource F18PointSource1
#/gate/source/F18PointSource1/setActivity 0.0001 Ci
#/gate/source/F18PointSource1/gps/particle e+
#/gate/source/F18PointSource1/setForcedUnstableFlag true
#/gate/source/F18PointSource1/setForcedHalfLife 6586.2 s
#/gate/source/F18PointSource1/gps/energytype Fluor18
#/gate/source/F18PointSource1/gps/type Point
#/gate/source/F18PointSource1/gps/angtype iso
#/gate/source/F18PointSource1/gps/centre 0. 0. 0. cm
#/gate/source/F18PointSource1/gps/confine NULL

/gate/source/list



#	O U T P U T 

/gate/output/verbose 0

# Setup Ascii file
/gate/output/ascii/setFileName dhci_Sen10cm_bkb_test
/gate/output/ascii/setOutFileSinglesAdderFlag 0
/gate/output/ascii/setOutFileSinglesReadoutFlag 0
/gate/output/ascii/setOutFileSinglesThresholderFlag 0
/gate/output/ascii/setOutFileSinglesUpholderFlag 0
/gate/output/ascii/setOutFileHitsFlag 0
/gate/output/ascii/setOutFileSinglesFlag 0
/gate/output/ascii/setOutFileCoincidencesFlag 1

#Setup root file
#/gate/output/root/setFileName dhci_
#/gate/output/root/setRootHitFlag 0
#/gate/output/root/setRootSinglesFlag 0
#/gate/output/root/setRootCoincidencesFlag 0
#/gate/output/root/setRootNtupleFlag 0
#/gate/output/root/setRootSinglesAdderFlag 0
#/gate/output/root/setRootSinglesReadoutFlag 0
#

# Chnage the seed(1)
/gate/output/root/setSaveRndmFlag 1

#	Experiment 
/gate/application/setTimeSlice     .005 s
/gate/application/setTimeStart     0.  s
/gate/application/setTimeStop      .05 s

# Start acquistion
/gate/application/startDAQ
exit
-------------- next part --------------
/control/verbose 0

#######################################
###               V I Z             ###
#######################################

#       V I S U A L I S A T I O N
/vis/open OGLSX
/vis/viewer/reset
/vis/viewer/viewpointThetaPhi 60 60
/vis/viewer/zoom 2
/vis/viewer/set/style surface
/vis/drawVolume
/tracking/storeTrajectory 1
/vis/scene/endOfEventAction accumulate
/vis/viewer/update
/vis/verbose 2
/gate/geometry/enableAutoUpdate


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

#######################################
###           C A M E R A           ###
#######################################

       CYLINDRICAL
/gate/world/daughters/name cylindricalPET
/gate/world/daughters/insert cylinder
/gate/cylindricalPET/placement/setTranslation 0.0 0.0 0.0 cm
/gate/cylindricalPET/geometry/setRmax 7 cm
/gate/cylindricalPET/geometry/setRmin 3.9 cm
/gate/cylindricalPET/geometry/setHeight 8.1 cm
/gate/cylindricalPET/setMaterial Air
/gate/cylindricalPET/vis/forceWireframe
/gate/cylindricalPET/vis/setColor white

#       rsector
/gate/cylindricalPET/daughters/name rsector
/gate/cylindricalPET/daughters/insert box
/gate/rsector/placement/setTranslation 5.1 0.0 0.0 cm
/gate/rsector/geometry/setXLength 2.2 cm
/gate/rsector/geometry/setYLength 1.1 cm
/gate/rsector/geometry/setZLength 8.1 cm
/gate/rsector/setMaterial Air
/gate/rsector/vis/setColor red

#       M O D U L E
/gate/rsector/daughters/name module
/gate/rsector/daughters/insert box
/gate/module/placement/setTranslation 0.0 0.0 0.0 cm
/gate/module/geometry/setXLength 1.1 cm
/gate/module/geometry/setYLength 1.1 cm
/gate/module/geometry/setZLength 1.3 mm
/gate/module/setMaterial Air
/gate/module/vis/setColor red


#       C R Y S T A L
/gate/module/daughters/name crystal
/gate/module/daughters/insert box
/gate/crystal/placement/setTranslation 0.0 0.0 0.0 cm
/gate/crystal/geometry/setXLength 3.0 mm
/gate/crystal/geometry/setYLength 1.0 mm
/gate/crystal/geometry/setZLength 1.3 mm
/gate/crystal/setMaterial Air
/gate/crystal/vis/forceWireframe
/gate/crystal/vis/setColor yellow


#       LSO layer
/gate/crystal/daughters/name LSO
/gate/crystal/daughters/insert box
/gate/LSO/placement/setTranslation 0.0 0.0 0.0 cm
/gate/LSO/geometry/setXLength 3.0 mm
/gate/LSO/geometry/setYLength 1.0 mm
/gate/LSO/geometry/setZLength 1.0 mm
/gate/LSO/setMaterial LSO
/gate/LSO/vis/setColor green

#       R E P E A T    C R Y S T A L
/gate/crystal/repeaters/insert cubicArray
/gate/crystal/cubicArray/setRepeatNumberX 3
/gate/crystal/cubicArray/setRepeatNumberY 8
/gate/crystal/cubicArray/setRepeatNumberZ 1
/gate/crystal/cubicArray/setRepeatVector 3.0 1.0 0.0 mm

#       R E P E A T    M O D U L E
/gate/module/repeaters/insert cubicArray
/gate/module/cubicArray/setRepeatNumberX 2
/gate/module/cubicArray/setRepeatNumberY 1
/gate/module/cubicArray/setRepeatNumberZ 62
/gate/module/cubicArray/setRepeatVector 11.0 0.0 1.3 mm

#       R E P E A T rsector
/gate/rsector/repeaters/insert ring
/gate/rsector/ring/setRepeatNumber 23

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

#       A T T A C H    C R Y S T A L  SD

/gate/LSO/attachCrystalSD

#       D I G I T I Z E R
/gate/digitizer/Singles/insert adder
/gate/digitizer/Singles/insert readout
/gate/digitizer/Singles/readout/setDepth 1
/gate/digitizer/Singles/insert blurring
/gate/digitizer/Singles/blurring/setResolution 0.12
/gate/digitizer/Singles/blurring/setEnergyOfReference 511. keV
/gate/digitizer/Singles/insert thresholder
/gate/digitizer/Singles/thresholder/setThreshold 10. keV
/gate/digitizer/Singles/insert upholder
/gate/digitizer/Singles/upholder/setUphold 1050. keV

#####Temporal resolution#####
/gate/digitizer/Singles/insert timeResolution
/gate/digitizer/Singles/timeResolution/setTimeResolution 2 ns
/gate/digitizer/Singles/timeResolution/verbose 0

################
# DEADTIME
/gate/digitizer/Singles/insert deadtime
/gate/digitizer/Singles/deadtime/setDeadTime 5000. ns
/gate/digitizer/Singles/deadtime/setMode paralysable
/gate/digitizer/Singles/deadtime/chooseDTVolume module
/gate/digitizer/Singles/deadtime/verbose 0

/gate/digitizer/Singles/verbose 0
/gate/digitizer/verbose 0


#       C O I N C I    S O R T E R
/gate/digitizer/Coincidences/setWindow 20. ns
/gate/digitizer/Coincidences/minSectorDifference 1

/gate/systems/cylindricalPET/describe


#/control/execute phantom.mac
##########################################
###        P H A N T O M               ###
##########################################
#       P H A N T O M
/gate/world/daughters/name Breastrod
/gate/world/daughters/insert sphere
/gate/Breastrod/setMaterial Water
/gate/Breastrod/vis/setColor grey
/gate/Breastrod/geometry/setRmax   0.4 cm
/gate/Breastrod/geometry/setRmin   0.0 cm

/gate/Breastrod/attachPhantomSD
/gate/Breastrod/placement/setTranslation 0.0 0.0 -0.60 cm


#/control/execute physics.mac
##########################################
###        P H Y S I C S               ###
##########################################
#
#       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 30. cm

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

/gate/physics/verbose 2


       I N I T I A L I Z E
/gate/systems/cylindricalPET/verbose 0
/gate/geometry/enableAutoUpdate
/run/initialize

#       V E R B O S I T Y
/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


#/control/execute sources.mac
##########################################
###          S O u R C E S             ###
##########################################
#       S O U R C E
#  1 millicurie (mCi) = 37 megabecquerel (MBq)

# Source ZERO
/gate/source/addSource F18LineSource0
/gate/source/F18LineSource0/setActivity 0.0001 Ci

/gate/source/F18LineSource0/gps/particle e+
/gate/source/F18LineSource0/setForcedUnstableFlag true
/gate/source/F18LineSource0/setForcedHalfLife 6586.2 s
/gate/source/F18LineSource0/gps/energytype Fluor18
/gate/source/F18LineSource0/gps/type Point
/gate/source/F18LineSource0/gps/angtype iso
/gate/source/F18LineSource0/gps/centre 0.0 0.0 -0.60 cm
/gate/source/F18LineSource0/gps/confine NULL

/gate/source/list



#	O U T P U T 

/gate/output/verbose 0

#/gate/output/lmf1/disable

/gate/output/ascii/setFileName sasens_ring_np60
/gate/output/ascii/setOutFileHitsFlag 0
/gate/output/ascii/setOutFileSinglesFlag 0
/gate/output/ascii/setOutFileCoincidencesFlag 1

#/gate/output/root/setFileName bNEC_0275
#/gate/output/root/setRootHitFlag 0
#/gate/output/root/setRootSinglesFlag 0
#/gate/output/root/setRootCoincidencesFlag 1
#/gate/output/root/setRootNtupleFlag 0
/gate/output/root/setSaveRndmFlag 1


#       S T A R T
/gate/application/setTimeSlice     .03 s
/gate/application/setTimeStart     0.  s
/gate/application/setTimeStop      .3 s
/gate/application/startDAQ

exit
-------------- next part --------------
/control/verbose 0

#######################################
###               V I Z             ###
#######################################

#       V I S U A L I S A T I O N
/vis/open OGLSX
/vis/viewer/reset
/vis/viewer/viewpointThetaPhi 60 60
/vis/viewer/zoom 2
/vis/viewer/set/style surface
/vis/drawVolume
/tracking/storeTrajectory 1
/vis/scene/endOfEventAction accumulate
/vis/viewer/update
/vis/verbose 2
/gate/geometry/enableAutoUpdate


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

#######################################
###           C A M E R A           ###
#######################################

       CYLINDRICAL
/gate/world/daughters/name cylindricalPET
/gate/world/daughters/insert cylinder
/gate/cylindricalPET/placement/setTranslation 0.0 0.0 0.0 cm
/gate/cylindricalPET/geometry/setRmax 7 cm
/gate/cylindricalPET/geometry/setRmin 3.9 cm
/gate/cylindricalPET/geometry/setHeight 8.1 cm
/gate/cylindricalPET/setMaterial Air
/gate/cylindricalPET/vis/forceWireframe
/gate/cylindricalPET/vis/setColor white

#       rsector
/gate/cylindricalPET/daughters/name rsector
/gate/cylindricalPET/daughters/insert box
/gate/rsector/placement/setTranslation 5.1 0.0 0.0 cm
/gate/rsector/geometry/setXLength 2.2 cm
/gate/rsector/geometry/setYLength 1.1 cm
/gate/rsector/geometry/setZLength 8.1 cm
/gate/rsector/setMaterial Air
/gate/rsector/vis/setColor red

#       M O D U L E
/gate/rsector/daughters/name module
/gate/rsector/daughters/insert box
/gate/module/placement/setTranslation 0.0 0.0 0.0 cm
/gate/module/geometry/setXLength 1.1 cm
/gate/module/geometry/setYLength 1.1 cm
/gate/module/geometry/setZLength 1.3 mm
/gate/module/setMaterial Air
/gate/module/vis/setColor red


#       C R Y S T A L
/gate/module/daughters/name crystal
/gate/module/daughters/insert box
/gate/crystal/placement/setTranslation 0.0 0.0 0.0 cm
/gate/crystal/geometry/setXLength 3.0 mm
/gate/crystal/geometry/setYLength 1.0 mm
/gate/crystal/geometry/setZLength 1.3 mm
/gate/crystal/setMaterial Air
/gate/crystal/vis/forceWireframe
/gate/crystal/vis/setColor yellow


#       LSO layer
/gate/crystal/daughters/name LSO
/gate/crystal/daughters/insert box
/gate/LSO/placement/setTranslation 0.0 0.0 0.0 cm
/gate/LSO/geometry/setXLength 3.0 mm
/gate/LSO/geometry/setYLength 1.0 mm
/gate/LSO/geometry/setZLength 1.0 mm
/gate/LSO/setMaterial LSO
/gate/LSO/vis/setColor green

#       R E P E A T    C R Y S T A L
/gate/crystal/repeaters/insert cubicArray
/gate/crystal/cubicArray/setRepeatNumberX 3
/gate/crystal/cubicArray/setRepeatNumberY 8
/gate/crystal/cubicArray/setRepeatNumberZ 1
/gate/crystal/cubicArray/setRepeatVector 3.0 1.0 0.0 mm

#       R E P E A T    M O D U L E
/gate/module/repeaters/insert cubicArray
/gate/module/cubicArray/setRepeatNumberX 2
/gate/module/cubicArray/setRepeatNumberY 1
/gate/module/cubicArray/setRepeatNumberZ 62
/gate/module/cubicArray/setRepeatVector 11.0 0.0 1.3 mm

#       R E P E A T rsector
/gate/rsector/repeaters/insert ring
/gate/rsector/ring/setRepeatNumber 23

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

#       A T T A C H    C R Y S T A L  SD

/gate/LSO/attachCrystalSD

#       D I G I T I Z E R
/gate/digitizer/Singles/insert adder
/gate/digitizer/Singles/insert readout
/gate/digitizer/Singles/readout/setDepth 1
/gate/digitizer/Singles/insert blurring
/gate/digitizer/Singles/blurring/setResolution 0.12
/gate/digitizer/Singles/blurring/setEnergyOfReference 511. keV
/gate/digitizer/Singles/insert thresholder
/gate/digitizer/Singles/thresholder/setThreshold 10. keV
/gate/digitizer/Singles/insert upholder
/gate/digitizer/Singles/upholder/setUphold 1050. keV

#####Temporal resolution#####
/gate/digitizer/Singles/insert timeResolution
/gate/digitizer/Singles/timeResolution/setTimeResolution 2 ns
/gate/digitizer/Singles/timeResolution/verbose 0

################
# DEADTIME
/gate/digitizer/Singles/insert deadtime
/gate/digitizer/Singles/deadtime/setDeadTime 5000. ns
/gate/digitizer/Singles/deadtime/setMode paralysable
/gate/digitizer/Singles/deadtime/chooseDTVolume module
/gate/digitizer/Singles/deadtime/verbose 0

/gate/digitizer/Singles/verbose 0
/gate/digitizer/verbose 0


#       C O I N C I    S O R T E R
/gate/digitizer/Coincidences/setWindow 20. ns
/gate/digitizer/Coincidences/minSectorDifference 1

/gate/systems/cylindricalPET/describe


#/control/execute phantom.mac
##########################################
###        P H A N T O M               ###
##########################################
#       P H A N T O M
/gate/world/daughters/name Breastrod
/gate/world/daughters/insert sphere
/gate/Breastrod/setMaterial Water
/gate/Breastrod/vis/setColor grey
/gate/Breastrod/geometry/setRmax   0.4 cm
/gate/Breastrod/geometry/setRmin   0.0 cm

/gate/Breastrod/attachPhantomSD
/gate/Breastrod/placement/setTranslation 0.0 0.0 -0.60 cm


#/control/execute physics.mac
##########################################
###        P H Y S I C S               ###
##########################################
#
#       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 30. cm

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

/gate/physics/verbose 2


       I N I T I A L I Z E
/gate/systems/cylindricalPET/verbose 0
/gate/geometry/enableAutoUpdate
/run/initialize

#       V E R B O S I T Y
/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


#/control/execute sources.mac
##########################################
###          S O u R C E S             ###
##########################################
#       S O U R C E
#  1 millicurie (mCi) = 37 megabecquerel (MBq)

# Source ZERO
/gate/source/addSource F18LineSource0
/gate/source/F18LineSource0/setActivity 0.0001 Ci

/gate/source/F18LineSource0/gps/particle e+
/gate/source/F18LineSource0/setForcedUnstableFlag true
/gate/source/F18LineSource0/setForcedHalfLife 6586.2 s
/gate/source/F18LineSource0/gps/energytype Fluor18
/gate/source/F18LineSource0/gps/type Point
/gate/source/F18LineSource0/gps/angtype iso
/gate/source/F18LineSource0/gps/centre 0.0 0.0 -0.60 cm
/gate/source/F18LineSource0/gps/confine NULL

/gate/source/list



#	O U T P U T 

/gate/output/verbose 0

#/gate/output/lmf1/disable

/gate/output/ascii/setFileName sasens_ring_np60
/gate/output/ascii/setOutFileHitsFlag 0
/gate/output/ascii/setOutFileSinglesFlag 0
/gate/output/ascii/setOutFileCoincidencesFlag 1

#/gate/output/root/setFileName bNEC_0275
#/gate/output/root/setRootHitFlag 0
#/gate/output/root/setRootSinglesFlag 0
#/gate/output/root/setRootCoincidencesFlag 1
#/gate/output/root/setRootNtupleFlag 0
/gate/output/root/setSaveRndmFlag 1


#       S T A R T
/gate/application/setTimeSlice     .03 s
/gate/application/setTimeStart     0.  s
/gate/application/setTimeStop      .3 s
/gate/application/startDAQ

exit


More information about the Gate-users mailing list