[Gate-users] How to make dual offset layer

임원석 malcomlim at naver.com
Tue Sep 21 13:32:00 CEST 2021 

Hi
 
I am trying to make dual layer offset structure detector.
I tried to make it work by repeating the layers of LYSO. However there is an error saying
 
-------- EEEE ------- G4Exception-START -------- EEEE -------
*** G4Exception : Cache001
      issued by : G4CacheReference<V>::Destroy
Internal fatal error. Invalid G4Cache size (requested id: 8 but cache has size: 6 Possibly client created G4Cache object in a thread and tried to delete it from another thread!
*** Fatal Exception *** core dump ***
-------- EEEE -------- G4Exception-END --------- EEEE -------
 
[G4-cerr] 
[G4-cerr] *** G4Exception: Aborting execution ***
 
I don't have any idea what it means.
Can any one help?
 
Here is my code for dual offset layer 
 
/vis/disable
 
#/vis/open                           OGLSX
#/vis/viewer/set/viewpointThetaPhi   90 90
#/vis/viewer/zoom                    6.0
#/vis/drawVolume
#/tracking/storeTrajectory           1
#/vis/scene/add/trajectories 
#/vis/scene/endOfEventAction         accumulate
 
#/vis/open OGLSQt
#/vis/viewer/set/viewpointThetaPhi 0 0
#/vis/viewer/zoom 1.5
#/vis/viewer/set/projection perspective
#/vis/viewer/set/lightsMove camera
#/vis/drawVolume
#/tracking/storeTrajectory 1
#/vis/scene/add/trajectories
#/vis/scene/endOfEventAction accumulate -1
#/vis/scene/add/axes 
 
/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 25.8 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  20.0 mm
/gate/rsector/geometry/setYLength 25.8 mm
/gate/rsector/geometry/setZLength 25.8 mm
/gate/rsector/setMaterial Air
#/gate/rsector/vis/setColor blue
#/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 cm
/gate/module/geometry/setXLength 20.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
 
# 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 20 mm
/gate/submodule/geometry/setYLength 25.8 mm
/gate/submodule/geometry/setZLength 25.8 mm
/gate/submodule/setMaterial Air
/gate/submodule/vis/setVisible 0
 
# 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. 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 0
 
 
# LSO layer0
/gate/crystal/daughters/name LYSO
/gate/crystal/daughters/insert box
/gate/LYSO/placement/setTranslation -0.6 0.0 0.0 cm
/gate/LYSO/geometry/setXLength 8. 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 LYSO2
/gate/crystal/daughters/insert box
/gate/LYSO2/placement/setTranslation 0.4 0.0 0.0 cm
/gate/LYSO2/geometry/setXLength 12. mm
/gate/LYSO2/geometry/setYLength 1.75 mm
/gate/LYSO2/geometry/setZLength 1.75 mm
/gate/LYSO2/setMaterial LYSO
/gate/LYSO2/vis/setColor red
#/gate/LYSO2/vis/setVisible 0
 
 
 
# R E P E A T    L A Y E R
 
/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/LYSO2/repeaters/insert cubicArray
/gate/LYSO2/cubicArray/setRepeatNumberX 1
/gate/LYSO2/cubicArray/setRepeatNumberY 14
/gate/LYSO2/cubicArray/setRepeatNumberZ 14
/gate/LYSO2/cubicArray/setRepeatVector 0.0 1.83 1.83 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
 
 
 
/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
 
 
/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/LYSO/attachCrystalSD
/gate/LYSO2/attachCrystalSD
/gate/systems/cylindricalPET/describe
 
 
 
#       P H Y S I C S
#
/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 G4PositronAnnihilation e+
#/gate/physics/processes/G4PositronAnnihilation/setModel StandardModel
 
/gate/physics/addProcess MultipleScattering e+
/gate/physics/addProcess MultipleScattering e-
 
/gate/physics/addProcess RadioactiveDecay
 
/gate/physics/addAtomDeexcitation
#/process/em/fluo  false
#/process/em/auger false
#/process/em/pixe  false
#/gate/physics/addPhysicsList emstandard_opt3
 
/gate/physics/processList Enabled
/gate/physics/processList Initialized
 
 
/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 LYSO2 1.0 cm
/gate/physics/Electron/SetCutInRegion LYSO2 1.0 cm
/gate/physics/Positron/SetCutInRegion LYSO2 1.0 cm
 
#/gate/physics/Gamma/SetCutInRegion GAGG 1.0 cm
#/gate/physics/Electron/SetCutInRegion GAGG 1.0 cm
#/gate/physics/Positron/SetCutInRegion GAGG 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 2000. ps
/gate/digitizer/Coincidences/minSectorDifference 2
/gate/digitizer/Coincidences/MultiplesPolicy  takeWinnerOfGoods
/gate/digitizer/Coincidences/setOffset 0. ns
#/gate/digitizer/Coincidences/describe
 
 
/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  takeWinnerOfGoods
 
# S O U R C E
#################
#               #
/gate/source/addSource source1
/gate/source/source1/setActivity     1000000. Bq
/gate/source/source1/setType         backtoback
/gate/source/source1/gps/particle    gamma
/gate/source/source1/gps/energytype  Mono
/gate/source/source1/gps/energy      511 keV
/gate/source/source1/gps/angtype     iso
/gate/source/source1/gps/type        Volume
/gate/source/source1/gps/shape Sphere
/gate/source/source1/gps/radius 0.5 mm
#/gate/source/source1/gps/halfz 25 mm
/gate/source/source1/gps/centre      0. 0. 0. mm
/gate/source/source1/visualize 100 red 1
 
/gate/source/addSource source2
/gate/source/source1/setActivity     1000000. Bq
/gate/source/source1/setType         backtoback
/gate/source/source1/gps/particle    gamma
/gate/source/source1/gps/energytype  Mono
/gate/source/source1/gps/energy      511 keV
/gate/source/source1/gps/angtype     iso
/gate/source/source1/gps/type        Volume
/gate/source/source1/gps/shape Sphere
/gate/source/source1/gps/radius 0.5 mm
#/gate/source/source1/gps/halfz 25 mm
/gate/source/source1/gps/centre      20. 0. 0. mm
/gate/source/source1/visualize 100 red 1
 
/gate/source/addSource source3
/gate/source/source1/setActivity     1000000. Bq
/gate/source/source1/setType         backtoback
/gate/source/source1/gps/particle    gamma
/gate/source/source1/gps/energytype  Mono
/gate/source/source1/gps/energy      511 keV
/gate/source/source1/gps/angtype     iso
/gate/source/source1/gps/type        Volume
/gate/source/source1/gps/shape Sphere
/gate/source/source1/gps/radius 0.5 mm
#/gate/source/source1/gps/halfz 25 mm
/gate/source/source1/gps/centre      40. 0. 0. mm
/gate/source/source1/visualize 100 red 1
 
/gate/source/addSource source4
/gate/source/source1/setActivity     1000000. Bq
/gate/source/source1/setType         backtoback
/gate/source/source1/gps/particle    gamma
/gate/source/source1/gps/energytype  Mono
/gate/source/source1/gps/energy      511 keV
/gate/source/source1/gps/angtype     iso
/gate/source/source1/gps/type        Volume
/gate/source/source1/gps/shape Sphere
/gate/source/source1/gps/radius 0.5 mm
#/gate/source/source1/gps/halfz 25 mm
/gate/source/source1/gps/centre      60. 0. 0. mm
/gate/source/source1/visualize 100 red 1
 
/gate/source/addSource source5
/gate/source/source1/setActivity     1000000. Bq
/gate/source/source1/setType         backtoback
/gate/source/source1/gps/particle    gamma
/gate/source/source1/gps/energytype  Mono
/gate/source/source1/gps/energy      511 keV
/gate/source/source1/gps/angtype     iso
/gate/source/source1/gps/type        Volume
/gate/source/source1/gps/shape Sphere
/gate/source/source1/gps/radius 0.5 mm
#/gate/source/source1/gps/halfz 25 mm
/gate/source/source1/gps/centre      80. 0. 0. mm
/gate/source/source1/visualize 100 red 1
 
/gate/source/addSource source6
/gate/source/source1/setActivity     1000000. Bq
/gate/source/source1/setType         backtoback
/gate/source/source1/gps/particle    gamma
/gate/source/source1/gps/energytype  Mono
/gate/source/source1/gps/energy      511 keV
/gate/source/source1/gps/angtype     iso
/gate/source/source1/gps/type        Volume
/gate/source/source1/gps/shape Sphere
/gate/source/source1/gps/radius 0.5 mm
#/gate/source/source1/gps/halfz 25 mm
/gate/source/source1/gps/centre     100. 0. 0. mm
/gate/source/source1/visualize 100 red 1
 
/gate/source/list
 
 
 
# O U T P U T 
 
/gate/output/root/enable
/gate/output/root/setFileName center_DOI_2layer_20mm_1.5mm
/gate/output/root/setRootHitFlag 0
/gate/output/root/setRootSinglesFlag 0
/gate/output/root/setRootCoincidencesFlag 1
/gate/output/root/setRootdelayFlag 1
 
 
#   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     180 s
 
/gate/application/startDAQ
 
 
Thank you 
Lim Wonsuk
 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20210921/8190176d/attachment-0001.html>

More information about the Gate-users mailing list