[Gate-users] [GATE-users] Multi-system digitizer code
김연경
dusrud026 at gmail.com
Tue Jul 9 07:28:26 CEST 2019
Dear GATE users,
My multi-system consists of two cylindrical PETs.
I try to coincidence sorting by combining single events from two
cylindrical PETs.
Therefore, I would like to get the same results when using the same PET as
multi-system and general cylindrical PET.
However, the digitizer part still is a problem.
Please let me know if you have any sample code or documentation that can
resolve this problem.
Here is my code:
/gate/geometry/setMaterialDatabase
/GATE/GATE_v8.0/Gate_Work/GateMaterials.db
/gate/world/geometry/setXLength 150. cm
/gate/world/geometry/setYLength 150. cm
/gate/world/geometry/setZLength 150. cm
/gate/world/daughters/name cylindricalPET
/gate/world/daughters/systemType PETscanner
/gate/world/daughters/insert cylinder
/gate/cylindricalPET/placement/setTranslation 0.0 0.0 51.6 mm
/gate/cylindricalPET/geometry/setRmax 160 mm
/gate/cylindricalPET/geometry/setRmin 138.5 mm
/gate/cylindricalPET/geometry/setHeight 103.2 mm
/gate/cylindricalPET/setMaterial Air
/gate/cylindricalPET/vis/setVisible 0
/gate/cylindricalPET/daughters/name rsector
/gate/cylindricalPET/daughters/insert box
/gate/rsector/placement/setTranslation 146 0.0 0. mm
/gate/rsector/geometry/setXLength 15.0 mm
/gate/rsector/geometry/setYLength 25.8 mm
/gate/rsector/geometry/setZLength 103.2 mm
/gate/rsector/setMaterial Air
/gate/rsector/vis/setColor red
/gate/rsector/daughters/name module
/gate/rsector/daughters/insert box
/gate/module/placement/setTranslation 0. 0.0 0.0 cm
/gate/module/geometry/setXLength 15.0 mm
/gate/module/geometry/setYLength 25.8 mm
/gate/module/geometry/setZLength 25.8 mm
/gate/module/setMaterial Air
/gate/module/vis/setVisible 0
/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 15.0 mm
/gate/crystal/geometry/setYLength 3.11 mm
/gate/crystal/geometry/setZLength 3.11 mm
/gate/crystal/setMaterial Air
/gate/crystal/daughters/name LYSO
/gate/crystal/daughters/insert box
/gate/LYSO/placement/setTranslation 0.0 0.0 0.0 cm
/gate/LYSO/geometry/setXLength 15. mm
/gate/LYSO/geometry/setYLength 3.11 mm
/gate/LYSO/geometry/setZLength 3.11 mm
/gate/LYSO/setMaterial LYSO
/gate/LYSO/vis/setColor green
/gate/crystal/repeaters/insert cubicArray
/gate/crystal/cubicArray/setRepeatNumberX 1
/gate/crystal/cubicArray/setRepeatNumberY 8
/gate/crystal/cubicArray/setRepeatNumberZ 8
/gate/crystal/cubicArray/setRepeatVector 0.0 3.19 3.19 mm
/gate/module/repeaters/insert cubicArray
/gate/module/cubicArray/setRepeatNumberX 1
/gate/module/cubicArray/setRepeatNumberY 1
/gate/module/cubicArray/setRepeatNumberZ 4
/gate/module/cubicArray/setRepeatVector 0.00 25.8 25.8 mm
/gate/rsector/repeaters/insert ring
/gate/rsector/ring/setRepeatNumber 34
/gate/world/daughters/name Scanner_1
/gate/world/daughters/systemType PETscanner
/gate/world/daughters/insert cylinder
/gate/Scanner_1/placement/setTranslation 0.0 0.0 -51.6 mm
/gate/Scanner_1/geometry/setRmax 160 mm
/gate/Scanner_1/geometry/setRmin 138.5 mm
/gate/Scanner_1/geometry/setHeight 103.2 mm
/gate/Scanner_1/setMaterial Air
/gate/Scanner_1/vis/setVisible 0
/gate/Scanner_1/daughters/name rsector1
/gate/Scanner_1/daughters/insert box
/gate/rsector1/placement/setTranslation 146 0.0 0. mm
/gate/rsector1/geometry/setXLength 15.0 mm
/gate/rsector1/geometry/setYLength 25.8 mm
/gate/rsector1/geometry/setZLength 103.2 mm
/gate/rsector1/setMaterial Air
/gate/rsector1/vis/setColor red
/gate/rsector1/vis/setVisible 0
/gate/rsector1/daughters/name module1
/gate/rsector1/daughters/insert box
/gate/module1/placement/setTranslation 0. 0.0 0.0 cm
/gate/module1/geometry/setXLength 15.0 mm
/gate/module1/geometry/setYLength 25.8 mm
/gate/module1/geometry/setZLength 25.8 mm
/gate/module1/setMaterial Air
/gate/module1/vis/setVisible 0
/gate/module1/daughters/name crystal1
/gate/module1/daughters/insert box
/gate/crystal1/placement/setTranslation 0.0 0.0 0.0 cm
/gate/crystal1/geometry/setXLength 15.0 mm
/gate/crystal1/geometry/setYLength 3.11 mm
/gate/crystal1/geometry/setZLength 3.11 mm
/gate/crystal1/setMaterial Air
/gate/crystal1/daughters/name LYSO1
/gate/crystal1/daughters/insert box
/gate/LYSO1/placement/setTranslation 0.0 0.0 0.0 cm
/gate/LYSO1/geometry/setXLength 15. mm
/gate/LYSO1/geometry/setYLength 3.11 mm
/gate/LYSO1/geometry/setZLength 3.11 mm
/gate/LYSO1/setMaterial LYSO
/gate/LYSO1/vis/setColor green
/gate/crystal1/repeaters/insert cubicArray
/gate/crystal1/cubicArray/setRepeatNumberX 1
/gate/crystal1/cubicArray/setRepeatNumberY 8
/gate/crystal1/cubicArray/setRepeatNumberZ 8
/gate/crystal1/cubicArray/setRepeatVector 0.0 3.19 3.19 mm
/gate/module1/repeaters/insert cubicArray
/gate/module1/cubicArray/setRepeatNumberX 1
/gate/module1/cubicArray/setRepeatNumberY 1
/gate/module1/cubicArray/setRepeatNumberZ 4
/gate/module1/cubicArray/setRepeatVector 0.00 25.8 25.8 mm
/gate/rsector1/repeaters/insert ring
/gate/rsector1/ring/setRepeatNumber 34
/gate/systems/cylindricalPET/level1/attach rsector
/gate/systems/cylindricalPET/level2/attach module
/gate/systems/cylindricalPET/level3/attach crystal
/gate/systems/cylindricalPET/level1/attach LYSO
/gate/systems/Scanner_1/level1/attach rsector1
/gate/systems/Scanner_1/level2/attach module1
/gate/systems/Scanner_1/level3/attach crystal1
/gate/systems/Scanner_1/level1/attach LYSO1
/gate/LYSO/attachCrystalSD
/gate/LYSO1/attachCrystalSD
/gate/world/daughters/name NEMACylinder
/gate/world/daughters/insert sphere
/gate/NEMACylinder/placement/setTranslation 0.0 0.0 -55 mm
/gate/NEMACylinder/setMaterial Water
/gate/NEMACylinder/geometry/setRmax 110. mm
/gate/NEMACylinder/geometry/setRmin 0. mm
/gate/NEMACylinder/geometry/setThetaStart 0. deg
/gate/NEMACylinder/geometry/setDeltaTheta 90. deg
/gate/NEMACylinder/vis/forceWireframe
/gate/NEMACylinder/vis/setColor blue
/gate/NEMACylinder/attachPhantomSD
/gate/physics/addProcess PhotoElectric
/gate/physics/processes/PhotoElectric/setModel StandardModel
/gate/physics/addProcess Compton
/gate/physics/processes/Compton/setModel StandardModel
/gate/physics/addProcess ElectronIonisation
/gate/physics/processes/ElectronIonisation/setModel StandardModel e-
/gate/physics/processes/ElectronIonisation/setModel StandardModel e+
/gate/physics/addProcess Bremsstrahlung
/gate/physics/processes/Bremsstrahlung/setModel StandardModel e-
/gate/physics/processes/Bremsstrahlung/setModel StandardModel e+
/gate/physics/addProcess PositronAnnihilation
/gate/physics/addProcess MultipleScattering e+
/gate/physics/addProcess MultipleScattering e-
/gate/physics/addAtomDeexcitation
/gate/physics/processList Enabled
/gate/physics/processList Initialized
/gate/physics/Gamma/SetCutInRegion NEMACylinder 1.0 cm
/gate/physics/Electron/SetCutInRegion NEMACylinder 1.0 cm
/gate/physics/Positron/SetCutInRegion NEMACylinder 1.0 cm
/gate/physics/Gamma/SetCutInRegion LYSO 1.0 cm
/gate/physics/Electron/SetCutInRegion LYSO 1.0 cm
/gate/physics/Positron/SetCutInRegion LYSO 1.0 cm
/gate/physics/Gamma/SetCutInRegion LYSO1 1.0 cm
/gate/physics/Electron/SetCutInRegion LYSO1 1.0 cm
/gate/physics/Positron/SetCutInRegion LYSO1 1.0 cm
/gate/run/initialize
/gate/digitizer/Singles/insert systemFilter
/gate/digitizer/Singles/systemFilter/selectSystem cylindricalPET
/gate/digitizer/Singles/insert adder
/gate/digitizer/Singles/insert readout
/gate/digitizer/Singles/readout/setDepth 4
/gate/digitizer/Singles/insert blurring
/gate/digitizer/Singles/blurring/setLaw linear
/gate/digitizer/Singles/blurring/linear/setResolution 0.13
/gate/digitizer/Singles/blurring/linear/setEnergyOfReference 511. keV
/gate/digitizer/Singles/insert thresholder
/gate/digitizer/Singles/thresholder/setThreshold 400. keV
/gate/digitizer/Singles/insert upholder
/gate/digitizer/Singles/upholder/setUphold 650. keV
/gate/digitizer/Singles/insert timeResolution
/gate/digitizer/Singles/timeResolution/setTimeResolution 220. ps
/gate/digitizer/Singles/insert deadtime
/gate/digitizer/Singles/deadtime/setDeadTime 150. ns
/gate/digitizer/Singles/deadtime/setMode paralysable
/gate/digitizer/Singles/deadtime/chooseDTVolume crystal
#Scanner_1 digitizer
/gate/digitizer/name Singles_S1
/gate/digitizer/insert singleChain
/gate/digitizer/Singles_S1/insert systemFilter
/gate/digitizer/Singles_S1/systemFilter/selectSystem Scanner_1
/gate/digitizer/Singles_S1/insert adder
/gate/digitizer/Singles_S1/insert readout
/gate/digitizer/Singles_S1/readout/setDepth 4
/gate/digitizer/Singles_S1/insert blurring
/gate/digitizer/Singles_S1/blurring/setLaw linear
/gate/digitizer/Singles_S1/blurring/linear/setResolution 0.13
/gate/digitizer/Singles_S1/blurring/linear/setEnergyOfReference 511. keV
/gate/digitizer/Singles_S1/insert thresholder
/gate/digitizer/Singles_S1/thresholder/setThreshold 400. keV
/gate/digitizer/Singles_S1/insert upholder
/gate/digitizer/Singles_S1/upholder/setUphold 650. keV
/gate/digitizer/Singles_S1/insert timeResolution
/gate/digitizer/Singles_S1/timeResolution/setTimeResolution 220. ps
/gate/digitizer/Singles_S1/insert deadtime
/gate/digitizer/Singles_S1/deadtime/setDeadTime 150. ns
/gate/digitizer/Singles_S1/deadtime/setMode paralysable
/gate/digitizer/Singles_S1/deadtime/chooseDTVolume crystal
/gate/digitizer/Coincidences/setInputName Singles
/gate/digitizer/Coincidences/setOffset 0. ns
/gate/digitizer/Coincidences/setWindow 440. ps
/gate/digitizer/Coincidences/minSectorDifference 2
/gate/digitizer/Coincidences/setInputName Singles_S1
/gate/digitizer/Coincidences/setOffset 0. ns
/gate/digitizer/Coincidences/setWindow 440. ps
/gate/digitizer/Coincidences/minSectorDifference 2
/gate/digitizer/name finalCoinc
/gate/digitizer/insert coincidenceChain
/gate/digitizer/finalCoinc/addInputName Singles
/gate/digitizer/finalCoinc/addInputName Singles_S1
/gate/digitizer/finalCoinc/usePriority true
/gate/digitizer/name delay
/gate/digitizer/insert coincidenceSorter
/gate/digitizer/delay/minSectorDifference 2
/gate/digitizer/delay/setOffset 100. ns
/gate/digitizer/delay/setWindow 8. ns
/gate/digitizer/delay/MultiplesPolicy takeAllGoods
#####################################################
# Source
#
#####################################################
/gate/source/addSource F18LineSource
/gate/source/F18LineSource/setActivity 1000000. becquerel
/gate/source/F18LineSource/gps/particle e+
/gate/source/F18LineSource/setForcedUnstableFlag true
/gate/source/F18LineSource/setForcedHalfLife 6586.2 s
/gate/source/F18LineSource/gps/energytype Fluor18
/gate/source/F18LineSource/gps/type Volume
/gate/source/F18LineSource/gps/shape Cylinder
/gate/source/F18LineSource/gps/radius 200 mm
/gate/source/F18LineSource/gps/halfz 70.0 cm
/gate/source/F18LineSource/gps/angtype iso
/gate/source/F18LineSource/gps/centre 0. 0. -55. mm
/gate/source/F18LineSource/gps/confine NEMACylinder
/gate/source/list
# O U T P U T
/gate/output/root/enable
/gate/output/root/setFileName 190422_nonTOF_257_t20_30s
/gate/output/root/setRootHitFlag 0
/gate/output/root/setRootSinglesFlag 0
/gate/output/root/setRootCoincidencesFlag 1
/gate/output/root/setRootdelayFlag 0
# R A N D O M
# JamesRandom Ranlux64 MersenneTwister
/gate/random/setEngineName MersenneTwister
#/gate/random/setEngineSeed default
#/gate/random/setEngineSeed auto
/gate/random/setEngineSeed 123456789
#/gate/random/resetEngineFrom fileName
/gate/random/verbose 1
Thank you in advance for your help.
Best regards
-
YEONKYEONG KIM
Email : dusrud026 at gmail.com
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20190709/7dca3d79/attachment-0001.html>
More information about the Gate-users
mailing list