[Gate-users] Offset layer geometry

임원석 malcomlim at naver.com
Thu Sep 23 07:03:37 CEST 2021


Hi
 
I am trying to make dual layer offset structure detector.
I tried many ways to make the geometry and get the simulation done. 
I want to make the top-layer by 13x13 and the bottom-layer by 14x14.
Can anyone let me know how to make it?
 
Here is a code that I made which didn't work
 
#=====================================================
# VISUALISATION
#=====================================================
#/vis/disable
/vis/open OGLSQt
/vis/viewer/set/viewpointThetaPhi 0 90
#/vis/viewer/zoom 6.0
/vis/drawVolume
#/tracking/storeTrajectory 1
#/vis/scene/add/trajectories
#/vis/scene/endOfEventAction accumulate

#=====================================================
# VERBOSE
#=====================================================
/gate/verbose Physic 0
/gate/verbose Cuts 0
/gate/verbose SD 0
/gate/verbose Actions 0
/gate/verbose Actor 0
/gate/verbose Step 0
/gate/verbose Error 0
/gate/verbose Warning 0
/gate/verbose Output 0
/gate/verbose Beam 0
/gate/verbose Volume 0
/gate/verbose Image 0
/gate/verbose Geometry 0
/gate/verbose Core 0
/run/verbose 0
/event/verbose 0
/tracking/verbose 0

/gate/geometry/setMaterialDatabase /GATE/GATE_v8.0/Gate_Work/GateMaterials.db

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

# 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 210 mm
/gate/cylindricalPET/geometry/setRmin 165 mm
/gate/cylindricalPET/geometry/setHeight 26 mm
/gate/cylindricalPET/setMaterial Air
/gate/cylindricalPET/vis/forceWireframe
/gate/cylindricalPET/vis/setColor white
/gate/cylindricalPET/vis/setVisible 0


# HEAD
/gate/cylindricalPET/daughters/name rsector
/gate/cylindricalPET/daughters/insert box
/gate/rsector/placement/setTranslation 175 0.0 0 mm
/gate/rsector/geometry/setXLength  21 mm
/gate/rsector/geometry/setYLength 25.8 mm
/gate/rsector/geometry/setZLength 25.8 mm
/gate/rsector/setMaterial Air
#/gate/rsector/vis/setColor gray
/gate/rsector/vis/forceWireframe
#/gate/rsector/vis/setVisible 0



# MODULE
/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 21 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/vis/forceWireframe

# BLOCK
/gate/module/daughters/name submodule
/gate/module/daughters/insert box
/gate/submodule/placement/setTranslation 0.0 0.0 0.0 cm
/gate/submodule/geometry/setXLength 21 mm
/gate/submodule/geometry/setYLength 25.8 mm
/gate/submodule/geometry/setZLength 25.8 mm
/gate/submodule/setMaterial Air
/gate/submodule/vis/setVisible 1

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



# LSO layer0
/gate/crystal/daughters/name LYSO
/gate/crystal/daughters/insert box
/gate/LYSO/placement/setTranslation -0.6 0.0 0. cm
/gate/LYSO/geometry/setXLength 8.0 mm
/gate/LYSO/geometry/setYLength 1.75 mm
/gate/LYSO/geometry/setZLength 1.75 mm
/gate/LYSO/setMaterial LYSO
/gate/LYSO/vis/setColor green
#/gate/LYSO/vis/setVisible 0



# LSO layer1
/gate/crystal/daughters/name LYSO1
/gate/crystal/daughters/insert box
/gate/LYSO1/placement/setTranslation 0.4 0.0 0.0 cm
/gate/LYSO1/geometry/setXLength 12.0 mm
/gate/LYSO1/geometry/setYLength 1.75 mm
/gate/LYSO1/geometry/setZLength 1.75 mm
/gate/LYSO1/setMaterial LYSO
/gate/LYSO1/vis/setColor red
#/gate/LYSO1/vis/setVisible 0



# R E P E A T    C R Y S T A L

/gate/LYSO/repeaters/insert cubicArray
/gate/LYSO/cubicArray/setRepeatNumberX 1
/gate/LYSO/cubicArray/setRepeatNumberY 13
/gate/LYSO/cubicArray/setRepeatNumberZ 13
/gate/LYSO/cubicArray/setRepeatVector 0.0 1.83 1.83 mm

/gate/LYSO1/repeaters/insert cubicArray
/gate/LYSO1/cubicArray/setRepeatNumberX 1
/gate/LYSO1/cubicArray/setRepeatNumberY 14
/gate/LYSO1/cubicArray/setRepeatNumberZ 14
/gate/LYSO1/cubicArray/setRepeatVector 0.0 1.83 1.83 mm

/gate/crystal/repeaters/insert cubicArray
/gate/crystal/cubicArray/setRepeatNumberX 1
/gate/crystal/cubicArray/setRepeatNumberY 1
/gate/crystal/cubicArray/setRepeatNumberZ 1
/gate/crystal/cubicArray/setRepeatVector 0.0 25.8 25.8 mm



###################################################################################


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


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



#/gate/rsector/repeaters/insert linear
#/gate/rsector/linear/setRepeatNumber 2
#/gate/rsector/linear/setRepeatVector 0. 0. 25.8 mm

# 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/submodule/attach submodule
/gate/systems/cylindricalPET/crystal/attach crystal
/gate/systems/cylindricalPET/layer0/attach LYSO
/gate/systems/cylindricalPET/layer1/attach LYSO1


# A T T A C H    C R Y S T A L  SD
/gate/LYSO/attachCrystalSD
/gate/LYSO1/attachCrystalSD
#/gate/LYSO3/attachCrystalSD
#/gate/LYSO4/attachCrystalSD
/gate/systems/cylindricalPET/describe


#/gate/world/daughters/name phantom
#/gate/world/daughters/insert cylinder
#/gate/phantom/vis/forceWireframe
#/gate/phantom/placement/setTranslation 0. 0. 0. mm
#/gate/phantom/setMaterial Polyethylene
#/gate/phantom/geometry/setRmax 10. cm
#/gate/phantom/geometry/setHeight 5. cm
#/gate/phantom/vis/forceSolid
#/gate/phantom/vis/setColor white

#/gate/phantom/attachPhantomSD


#=====================================================
# PHYSICS
#=====================================================
#/gate/physics/addProcess PhotoElectric
#/gate/physics/processes/PhotoElectric/setModel StandardModel
#/gate/physics/addProcess Compton
#/gate/physics/processes/Compton/setModel StandardModel

#/gate/physics/addProcess RayleighScattering
#/gate/physics/processes/RayleighScattering/setModel PenelopeModel

#/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/processList Enabled
#/gate/physics/processList Initialized


/gate/physics/addPhysicsList emstandard_opt3

#=====================================================
# CUTS
#=====================================================
/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/physics/Gamma/SetCutInRegion LYSO3 1.0 cm
#/gate/physics/Electron/SetCutInRegion LYSO3 1.0 cm
#/gate/physics/Positron/SetCutInRegion LYSO3 1.0 cm
#/gate/physics/Gamma/SetCutInRegion LYSO4 1.0 cm
#/gate/physics/Electron/SetCutInRegion LYSO4 1.0 cm
#/gate/physics/Positron/SetCutInRegion LYSO4 1.0 cm
#/gate/physics/Gamma/SetCutInRegion phantom 1.0 cm
#/gate/physics/Electron/SetCutInRegion phantom 1.0 cm
#/gate/physics/Positron/SetCutInRegion phantom 1.0 cm
/gate/physics/displayCuts

/gate/run/initialize



# D I G I T I Z E R 
/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.15
/gate/digitizer/Singles/blurring/linear/setEnergyOfReference 511. keV

/gate/digitizer/Singles/insert thresholder
/gate/digitizer/Singles/thresholder/setThreshold 450. keV
/gate/digitizer/Singles/insert upholder
/gate/digitizer/Singles/upholder/setUphold 650. keV

# D E A D  T I M E
/gate/digitizer/Singles/insert deadtime
/gate/digitizer/Singles/deadtime/setDeadTime 150. ns
/gate/digitizer/Singles/deadtime/setMode nonparalysable
/gate/digitizer/Singles/deadtime/chooseDTVolume module



# C O I N C I    S O R T E R 
/gate/digitizer/Coincidences/setWindow 6. ns
/gate/digitizer/Coincidences/MultiplesPolicy  takeWinnerOfGoods

#/gate/digitizer/name delay
#/gate/digitizer/insert coincidenceSorter
#/gate/digitizer/delay/minSectorDifference 2
#/gate/digitizer/delay/setWindow 6. ns
#/gate/digitizer/delay/setOffset 500. ns
#/gate/digitizer/delay/MultiplesPolicy  takeWinnerOfGoods




#################
#               #
#  S O U R C E  #
#               #
#################

/gate/source/addSource SourceConfinement1
/gate/source/SourceConfinement1/gps/type Volume
/gate/source/SourceConfinement1/gps/shape Sphere
/gate/source/SourceConfinement1/gps/radius 0.5 mm
#/gate/source/SourceConfinement1/gps/radius 100. mm
#/gate/source/SourceConfinement1/gps/halfz 25. mm
/gate/source/SourceConfinement1/gps/centre 0. 0. 0. mm
/gate/source/SourceConfinement1/gps/particle gamma
#/gate/source/rod1/gps/energytype Mono #Fluor18
/gate/source/SourceConfinement1/setType backtoback
/gate/source/SourceConfinement1/gps/monoenergy 511. keV
/gate/source/SourceConfinement1/setActivity 1000000. Bq
/gate/source/SourceConfinement1/gps/angtype iso
#/gate/source/SourceConfinement1/visualize 1000 white 1
#/gate/source/SourceConfinement1/gps/confine phantom

/gate/source/addSource SourceConfinement2
/gate/source/SourceConfinement2/gps/type Volume
/gate/source/SourceConfinement2/gps/shape Sphere
/gate/source/SourceConfinement2/gps/radius 0.5 mm
#/gate/source/SourceConfinement2/gps/radius 100. mm
#/gate/source/SourceConfinement2/gps/halfz 25. mm
/gate/source/SourceConfinement2/gps/centre 20. 0. 0. mm
/gate/source/SourceConfinement2/gps/particle gamma
#/gate/source/rod1/gps/energytype Mono #Fluor18
/gate/source/SourceConfinement2/setType backtoback
/gate/source/SourceConfinement2/gps/monoenergy 511. keV
/gate/source/SourceConfinement2/setActivity 1000000. Bq
/gate/source/SourceConfinement2/gps/angtype iso
#/gate/source/SourceConfinement2/visualize 1000 white 1
#/gate/source/SourceConfinement2/gps/confine phantom


/gate/source/addSource SourceConfinement3
/gate/source/SourceConfinement3/gps/type Volume
/gate/source/SourceConfinement3/gps/shape Sphere
/gate/source/SourceConfinement3/gps/radius 0.5 mm
#/gate/source/SourceConfinement3/gps/radius 100. mm
#/gate/source/SourceConfinement3/gps/halfz 25. mm
/gate/source/SourceConfinement3/gps/centre 40. 0. 0. mm
/gate/source/SourceConfinement3/gps/particle gamma
#/gate/source/rod1/gps/energytype Mono #Fluor18
/gate/source/SourceConfinement3/setType backtoback
/gate/source/SourceConfinement3/gps/monoenergy 511. keV
/gate/source/SourceConfinement3/setActivity 1000000. Bq
/gate/source/SourceConfinement3/gps/angtype iso
/gate/source/SourceConfinement3/visualize 100 white 1
#/gate/source/SourceConfinement3/gps/confine phantom

/gate/source/addSource SourceConfinement4
/gate/source/SourceConfinement4/gps/type Volume
/gate/source/SourceConfinement4/gps/shape Sphere
/gate/source/SourceConfinement4/gps/radius 0.5 mm
#/gate/source/SourceConfinement4/gps/radius 100. mm
#/gate/source/SourceConfinement4/gps/halfz 25. mm
/gate/source/SourceConfinement4/gps/centre 60. 0. 0. mm
/gate/source/SourceConfinement4/gps/particle gamma
#/gate/source/rod1/gps/energytype Mono #Fluor18
/gate/source/SourceConfinement4/setType backtoback
/gate/source/SourceConfinement4/gps/monoenergy 511. keV
/gate/source/SourceConfinement4/setActivity 1000000. Bq
/gate/source/SourceConfinement4/gps/angtype iso
#/gate/source/SourceConfinement4/visualize 1000 white 1
#/gate/source/SourceConfinement4/gps/confine phantom

/gate/source/addSource SourceConfinement5
/gate/source/SourceConfinement5/gps/type Volume
/gate/source/SourceConfinement5/gps/shape Sphere
/gate/source/SourceConfinement5/gps/radius 0.5 mm
#/gate/source/SourceConfinement5/gps/radius 100. mm
#/gate/source/SourceConfinement5/gps/halfz 25. mm
/gate/source/SourceConfinement5/gps/centre 80. 0. 0. mm
/gate/source/SourceConfinement5/gps/particle gamma
#/gate/source/rod1/gps/energytype Mono #Fluor18
/gate/source/SourceConfinement5/setType backtoback
/gate/source/SourceConfinement5/gps/monoenergy 511. keV
/gate/source/SourceConfinement5/setActivity 1000000. Bq
/gate/source/SourceConfinement5/gps/angtype iso
/gate/source/SourceConfinement5/visualize 100 white 1
#/gate/source/SourceConfinement5/gps/confine phantom

/gate/source/addSource SourceConfinement6
/gate/source/SourceConfinement6/gps/type Volume
/gate/source/SourceConfinement6/gps/shape Sphere
/gate/source/SourceConfinement6/gps/radius 0.5 mm
#/gate/source/SourceConfinement6/gps/radius 100. mm
#/gate/source/SourceConfinement6/gps/halfz 25. mm
/gate/source/SourceConfinement6/gps/centre 100. 0. 0. mm
/gate/source/SourceConfinement6/gps/particle gamma
#/gate/source/rod1/gps/energytype Mono #Fluor18
/gate/source/SourceConfinement6/setType backtoback
/gate/source/SourceConfinement6/gps/monoenergy 511. keV
/gate/source/SourceConfinement6/setActivity 1000000. Bq
/gate/source/SourceConfinement6/gps/angtype iso
/gate/source/SourceConfinement6/visualize 100 white 1



/gate/source/list

# O U T P U T 
#/gate/output/ascii/disable
#/gate/output/ascii/setFileName 1
#/gate/output/ascii/setOutFileHitsFlag 0
#/gate/output/ascii/setOutFileSinglesFlag 0
#/gate/output/ascii/setOutFileCoincidencesFlag 0
#/gate/output/ascii/setOutFiledelayFlag 0
#/gate/output/ascii/setSingleMask 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0  
#/gate/output/ascii/setOutFileSizeLimit 0

/gate/output/root/enable
/gate/output/root/setFileName 0923_6p_test
/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

# S T A R T 
/gate/application/setTimeSlice     1. s
/gate/application/setTimeStart     0. s
/gate/application/setTimeStop      60. s
#/gate/application/startDAQ


Thank you

Lim Wonsuk 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20210923/e2dbcc5c/attachment-0001.html>


More information about the Gate-users mailing list