[Gate-users] Partial PET
Ane Etxebeste
ane.etxebeste at creatis.insa-lyon.fr
Wed Feb 16 13:59:23 CET 2022
Hi Zahra,
Regarding coincidences with a system made of two heads, this command
could solve your problem
/gate/digitizer/Coincidences/minSectorDifference 1
For the configuration you can try these options
/gate/head/repeaters/insert ring
/gate/head/ring/setFirstAngle
/gate/head/ring/setAngularSpan
/gate/head/ring/setRepeatNumber
You can find more details in
https://opengate.readthedocs.io/en/latest/defining_a_geometry.html?highlight=angularspan#ring-repeater.
I hope it helps,
Ane
On 16/02/2022 11:38, zahra ahmadi wrote:
> Dear all
> I used PET_Cylindryal program and changed the dimensions. In the full
> PET ring I have 12 head. I only need two head facing each other.
>
> /gate/head/ring/setRepeatNumber 2
>
> For full ring I have hits, singles and Coincidences. The problem for
> this partial ring is I haven't any Coincidences in output. I have
> Hits and Singles but the Coincidences file is empty.
> Another question is, If I want to make these two heads closer to each
> other, Is there any way for that.
>
> Please let me know if you have any idea to solve that. The program
> attached with this email.
> Best wishes
> Zahra
>
>
>
> =====================================================
> # VISUALISATION
> #=====================================================
>
> ##/vis/disable
>
> /vis/open OGLSQt
> /vis/viewer/set/viewpointThetaPhi 60 60
> /vis/viewer/zoom 3.0
> /vis/drawVolume
> /tracking/storeTrajectory 1
> /vis/scene/add/trajectories
> /vis/scene/add/axes
> /vis/scene/endOfEventAction accumulate
>
> #/control/execute Verbose.mac
>
>
> #=====================================================
> # GEOMETRY
> #=====================================================
>
> /gate/geometry/setMaterialDatabase /home/gate/Desktop/GateMaterials.db
>
>
>
> #
> # W O R L D
> #
> /gate/world/geometry/setXLength 400. cm
> /gate/world/geometry/setYLength 400. cm
> /gate/world/geometry/setZLength 400. cm
>
>
> #-------------------oooooOOOOO00000OOOOOooooo---------------------#
> # #
> # D E F I N I T I O N A N D D E S C R I T I O N #
> # O F TWO LSO PANELS #
> # #
> #-------------------oooooOOOOO00000OOOOOooooo---------------------#
>
> #CYLINDRICAL
> /gate/world/daughters/name cylindricalPET
> /gate/world/daughters/insert cylinder
> /gate/cylindricalPET/placement/setTranslation 0 0.0 0.0 cm
> /gate/cylindricalPET/geometry/setRmax 444.97 mm
> /gate/cylindricalPET/geometry/setRmin 424.97 mm
> /gate/cylindricalPET/geometry/setHeight 222.4 mm
> /gate/cylindricalPET/setMaterial Air
> /gate/cylindricalPET/vis/forceWireframe
> /gate/cylindricalPET/vis/setColor white
> #/gate/cylindricalPET/vis/setVisible True
>
> #HEAD
> /gate/cylindricalPET/daughters/name head
> /gate/cylindricalPET/daughters/insert box
> /gate/head/placement/setTranslation 434.97 0.0 0.0 mm
> /gate/head/geometry/setXLength 2.0 cm
> /gate/head/geometry/setYLength 222.4 mm
> /gate/head/geometry/setZLength 222.4 mm
> /gate/head/setMaterial Air
> /gate/head/vis/setColor white
> /gate/head/vis/setVisible 0
>
>
> #BLOCK
> /gate/head/daughters/name block
> /gate/head/daughters/insert box
> /gate/block/placement/setTranslation 0.0 0.0 0.0 cm
> /gate/block/geometry/setXLength 20.0 mm
> /gate/block/geometry/setYLength 55.6 mm
> /gate/block/geometry/setZLength 55.6 mm
> /gate/block/setMaterial Air
> /gate/block/vis/setVisible 0
>
> #C R Y S T A L
> /gate/block/daughters/name crystal
> /gate/block/daughters/insert box
> /gate/crystal/placement/setTranslation 0.0 0.0 0.0 cm
> /gate/crystal/geometry/setXLength 20 mm
> /gate/crystal/geometry/setYLength 4.0 mm
> /gate/crystal/geometry/setZLength 4.0 mm
> /gate/crystal/setMaterial Air
> /gate/crystal/vis/setVisible 0
>
>
> #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 20.0 mm
> /gate/LSO/geometry/setYLength 4.0 mm
> /gate/LSO/geometry/setZLength 4.0 mm
> /gate/LSO/setMaterial LSO
> /gate/LSO/vis/setColor yellow
> /gate/LSO/vis/forceSolid
>
>
> #R E P E A T C R Y S T A L
> /gate/crystal/repeaters/insert cubicArray
> /gate/crystal/cubicArray/setRepeatNumberX 1
> /gate/crystal/cubicArray/setRepeatNumberY 13
> /gate/crystal/cubicArray/setRepeatNumberZ 13
> /gate/crystal/cubicArray/setRepeatVector 0.0 4.0 4.0 mm
>
>
> #R E P E A T BLOCK
> /gate/block/repeaters/insert cubicArray
> /gate/block/cubicArray/setRepeatNumberX 1
> /gate/block/cubicArray/setRepeatNumberY 4
> /gate/block/cubicArray/setRepeatNumberZ 4
> /gate/block/cubicArray/setRepeatVector 0.0 55.6 55.6 mm
>
>
>
>
> #R E P E A T HEAD
> /gate/head/repeaters/insert ring
> /gate/head/ring/setRepeatNumber 2
>
>
>
> #A T T A C H S Y S T E M
> /gate/systems/cylindricalPET/rsector/attach head
> /gate/systems/cylindricalPET/module/attach block
> /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
>
> /gate/systems/cylindricalPET/describe
>
>
>
> #=====================================================
> # 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
>
> #=====================================================
> # CUTS
> #=====================================================
>
> /gate/physics/Gamma/SetCutInRegion LSO 1.0 cm
> /gate/physics/Electron/SetCutInRegion LSO 1.0 cm
> /gate/physics/Positron/SetCutInRegion LSO 1.0 cm
>
> #/gate/physics/Gamma/SetCutInRegion BGO 1.0 cm
> #/gate/physics/Electron/SetCutInRegion BGO 1.0 cm
> #/gate/physics/Positron/SetCutInRegion BGO 1.0 cm
>
>
> #/gate/physics/Gamma/SetCutInRegion phantom 0.1 mm
> #/gate/physics/Electron/SetCutInRegion phantom 0.1 mm
> #/gate/physics/Positron/SetCutInRegion phantom 0.1 mm
>
> #/gate/physics/SetMaxStepSizeInRegion phantom 0.01 mm
>
> #=====================================================
> # INITIALISATION
> #=====================================================
>
> /gate/run/initialize
>
>
> #=====================================================
> # Show how the geometry moves with time
> #=====================================================
>
> #/control/execute MoveVisu.mac
>
> #=====================================================
> # D I G I T I Z E R: DETECTOR ELECTRONIC RESPONSE
> #=====================================================
>
> /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.26
> /gate/digitizer/Singles/blurring/setEnergyOfReference 511. keV
>
> /gate/digitizer/Singles/insert thresholder
> /gate/digitizer/Singles/thresholder/setThreshold 350. keV
> /gate/digitizer/Singles/insert upholder
> /gate/digitizer/Singles/upholder/setUphold 650. keV
>
>
>
> #=====================================================
> #C O I N C I D E N C E S O R T E R
> #=====================================================
>
> /gate/digitizer/Coincidences/setWindow 10. ns
>
> /gate/digitizer/name delay
> /gate/digitizer/insert coincidenceSorter
> /gate/digitizer/delay/setWindow 10. ns
> /gate/digitizer/delay/setOffset 500. ns
>
> #=====================================================
> # P A R T I C L E S O U R C E
> #=====================================================
>
> /gate/source/addSource posiC11
>
> #/gate/source/posiC11/gps/particle e+
> #/gate/source/posiC11/gps/energytype Carbon11
> /gate/source/posiC11/setType backtoback
>
> # The particles emitted by the source are gammas
> /gate/source/posiC11/gps/particle gamma
>
> # The gammas have an energy of 511 keV
> /gate/source/posiC11/gps/energytype Mono
> /gate/source/posiC11/gps/monoenergy 0.511 MeV
>
> /gate/source/posiC11/setActivity 10000 becquerel
>
> /gate/source/posiC11/setForcedUnstableFlag true
> /gate/source/posiC11/setForcedHalfLife 1223 s
> /gate/source/posiC11/gps/angtype iso
>
> /gate/source/posiC11/visualize red
>
> /gate/source/list
>
> #=====================================================
> # D A T A O U T P U T
> #=====================================================
>
> /gate/output/root/enable
> /gate/output/root/setFileName YourFile
> /gate/output/root/setRootHitFlag 1
> /gate/output/root/setRootSinglesFlag 1
> /gate/output/root/setRootCoincidencesFlag 1
> /gate/output/root/setRootdelayFlag 1
>
> #====================================================
> # R A N D O M E N G I N E A N D S E E D
> #====================================================
> # JamesRandom Ranlux64 MersenneTwister
> /gate/random/setEngineName JamesRandom
> #/gate/random/setEngineSeed default
> #/gate/random/setEngineSeed auto
> #/gate/random/setEngineSeed 123456789
> /gate/random/setEngineSeed default
> #/gate/random/resetEngineFrom fileName
> /gate/random/verbose 1
>
> #=====================================================
> # M E A S U R E M E N T S E T T I N G S
> #=====================================================
>
> /gate/application/setTimeSlice 1 s
> /gate/application/setTimeStart 0 s
> /gate/application/setTimeStop 6 s
>
> /gate/application/startDAQ
>
>
>
>
>
>
>
>
> _______________________________________________
> Gate-users mailing list
> Gate-users at lists.opengatecollaboration.org
> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20220216/0bd4a7fa/attachment-0001.html>
More information about the Gate-users
mailing list