[Gate-users] ecatAccel - No Coincidences
Peter Bloomfield
peter.bloomfield at camhpet.ca
Fri Apr 27 22:28:18 CEST 2007
Dear Gaters
I have a macro that I hope simulates Siemens Biograph XVI HiRez, attached. The
problem I get is that no coincidences are recorded in the root output files.
If anyone has the time, can you look at it and correct my mistakes I have
made. If you could run it and check the root file to check for coincidences
that would be helpful to me as well. It is set to run for 10 s and it is a
37000 Bq f-18 point souce.
Thanks in advance for any suggestions,
Peter
-------------- next part --------------
# *********************************************************************
# * DISCLAIMER *
# * *
# * Neither the authors of this software system, nor their employing *
# * institutes, nor the agencies providing financial support for this *
# * work make any representation or warranty, express or implied, *
# * regarding this software system or assume any liability for its *
# * use. *
# * *
# * This code implementation is the intellectual property of the *
# * OpenGATE collaboration. *
# * By copying, distributing or modifying the Program (or any work *
# * based on the Program) you indicate your acceptance of this *
# * statement, and all its terms. *
# *********************************************************************
#
# This macro has been modified by:
# Peter B
#
# To execute:
# CheckIt
#-------------------oooooOOOOO00000OOOOOooooo---------------------#
# #
# P E T S C A N N E R G E O M E T R Y W I T H E C A T #
# A C C E l S Y S T E M D E F I N I T I O N #
# #
#-------------------oooooOOOOO00000OOOOOooooo---------------------#
#-------------------oooooOOOOO00000OOOOOooooo---------------------#
# #
# R U N T H E V I S U A L I S A T I O N A N D #
# C O N S T R U C T Y O U R E X P E R I M E N T A L A R E A #
# #
#-------------------oooooOOOOO00000OOOOOooooo---------------------#
#
# V I S U A L I S A T I O N
#
/vis/open OGLIX
/vis/viewer/reset
/vis/drawView 45 45
/vis/viewer/zoom 3
/vis/viewer/set/style surface
/vis/drawVolume
/tracking/storeTrajectory 1
/vis/scene/endOfEventAction accumulate
/vis/viewer/update
/gate/geometry/enableAutoUpdate
#
# W O R L D
#
/gate/world/geometry/setXLength 200. cm
/gate/world/geometry/setYLength 200. cm
/gate/world/geometry/setZLength 100. 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 Y O U R P E T D E V I C E #
# #
#-------------------oooooOOOOO00000OOOOOooooo---------------------#
#
# E C A T A C C E L
#
/gate/world/daughters/name ecatAccel
/gate/world/daughters/insert cylinder
/gate/ecatAccel/setMaterial Air
/gate/ecatAccel/geometry/setRmax 450.0 mm
/gate/ecatAccel/geometry/setRmin 420.0 mm
/gate/ecatAccel/geometry/setHeight 166.494 mm
/gate/ecatAccel/setTranslation 0.0 0.0 0.0 mm
#
# B L O C K
#
/gate/ecatAccel/daughters/name block
/gate/ecatAccel/daughters/insert box
/gate/block/geometry/setXLength 55.566 mm
/gate/block/geometry/setYLength 20.0 mm
/gate/block/geometry/setZLength 55.566 mm
/gate/block/setMaterial Air
#
# 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 mm
/gate/crystal/geometry/setXLength 3.904 mm
/gate/crystal/geometry/setYLength 20.0 mm
/gate/crystal/geometry/setZLength 3.904 mm
/gate/crystal/setMaterial LSO
#
# R E P E A T C R Y S T A L
#
/gate/crystal/repeaters/insert cubicArray
/gate/crystal/cubicArray/setRepeatNumberX 13
/gate/crystal/cubicArray/setRepeatNumberY 1
/gate/crystal/cubicArray/setRepeatNumberZ 13
/gate/crystal/cubicArray/setRepeatVector 3.969 0.0 3.969 mm
#
# R E P E A T BLOCK
#
/gate/block/repeaters/insert sphere
/gate/block/sphere/setRadius 424.493 mm
/gate/block/sphere/setRepeatNumberWithTheta 3
/gate/block/sphere/setRepeatNumberWithPhi 48
/gate/block/sphere/setThetaAngle 7.5 deg
#
# A T T A C H S Y S T E M
#
/gate/systems/ecatAccel/block/attach block
/gate/systems/ecatAccel/crystal/attach crystal
#
# A T T A C H S Y S T E M
#
/gate/crystal/attachCrystalSD
#-------------------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 Y O U R P H A N T O M #
# #
#-------------------oooooOOOOO00000OOOOOooooo---------------------#
#
# W A T E R F I L L E D C Y L I N D E R
#
/gate/world/daughters/name phantom
/gate/world/daughters/insert cylinder
/gate/phantom/setMaterial Water
/gate/phantom/vis/setColor green
/gate/phantom/geometry/setRmax 100.0 mm
/gate/phantom/geometry/setHeight 200.0 mm
/gate/phantom/attachPhantomSD
#-------------------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 Y O U R P H Y S I C S #
# #
#-------------------oooooOOOOO00000OOOOOooooo---------------------#
#
# 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 1. m
#
# I N A C T I V E X - R A Y S
#
/gate/physics/setXRayCut 1. GeV
/gate/physics/setDeltaRayCut 1. GeV
#-------------------oooooOOOOO00000OOOOOooooo---------------------#
# #
# I N I T I A L I Z A T I O N O F Y O U R #
# S I M U L A T I O N #
# #
#-------------------oooooOOOOO00000OOOOOooooo---------------------#
/gate/geometry/enableAutoUpdate
/run/initialize
#-------------------oooooOOOOO00000OOOOOooooo---------------------#
# #
# D E F I N I T I O N O F Y O U R A C Q U I S I T I O N #
# D I G I T I Z E R & C O I N C I D E N C E S H O R T E R #
# #
#-------------------oooooOOOOO00000OOOOOooooo---------------------#
#
# A D D E R
#
/gate/digitizer/Singles/insert adder
#
# R E A D O U T
#
/gate/digitizer/Singles/insert readout
/gate/digitizer/Singles/readout/setDepth 1
#
# E N E R G Y B L U R R I N G
#
/gate/digitizer/Singles/insert blurring
/gate/digitizer/Singles/blurring/setResolution 0.15
/gate/digitizer/Singles/blurring/setEnergyOfReference 511. keV
#
# E N E R G Y C U T
#
/gate/digitizer/Singles/insert thresholder
/gate/digitizer/Singles/thresholder/setThreshold 350. 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 3000000. ps
#/gate/digitizer/Singles/deadtime/setMode paralysable
#/gate/digitizer/Singles/deadtime/chooseDTVolume block
#
# C O I N C I S O R T E R
#
/gate/digitizer/Coincidences/setWindow 6.0 ns
# /gate/digitizer/Coincidences/setOffset 256.0 ns
/gate/digitizer/Coincidences/describe
# /gate/digitizer/name delay
# /gate/digitizer/insert coincidenceSorter
# /gate/digitizer/delay/setWindow 6.0 ns
# /gate/digitizer/delay/setOffset 500. ns
# /gate/digitizer/delay/describe
# /gate/digitizer/name finalCoinc
# /gate/digitizer/insert coincidenceChain
# /gate/digitizer/finalCoinc/addInputName delay
# /gate/digitizer/finalCoinc/addInputName Coincidences
# /gate/digitizer/finalCoinc/usePriority true
# /gate/digitizer/finalCoinc/describe
#-------------------oooooOOOOO00000OOOOOooooo---------------------#
# #
# D E F I N I T I O N O F #
# Y O U R O U T P U T F I L E #
# #
#-------------------oooooOOOOO00000OOOOOooooo---------------------#
#
# C H A N G E T H E S E E D (1) O R N O T (0)
#
/gate/output/root/setSaveRndmFlag 1
#
# S E T U P - R O O T _ F I L E
#
/gate/output/root/setFileName UniCyl-2
/gate/output/root/setRootSinglesAdderFlag 1
/gate/output/root/setRootSinglesReadoutFlag 1
/gate/output/root/setRootHitFlag 1
/gate/output/root/setRootSinglesFlag 1
/gate/output/root/setRootNtupleFlag 1
/gate/output/root/setRootCoincidencesFlag 1
#
# S E T U P - A S C I I _ F I L E
#
/gate/output/ascii/setOutFileSinglesAdderFlag 0
/gate/output/ascii/setOutFileSinglesReadoutFlag 0
/gate/output/ascii/setOutFileSinglesThresholderFlag 0
/gate/output/ascii/setOutFileSinglesUpholderFlag 0
/gate/output/ascii/setOutFileSinglesFlag 0
/gate/output/ascii/setOutFileHitsFlag 0
/gate/output/ascii/setOutFileCoincidencesFlag 0
###################################################
# E C A T A C C E L O U T P U T
#
# W I T H N E W F U N C T I O N A L I T I E S
#
#
###################################################
/gate/output/sinoAccel/disable
# /gate/output/sinoAccel/setTangCrystalBlurring 1.8 mm
# /gate/output/sinoAccel/setAxialCrystalBlurring 1.8 mm
# /gate/output/sinoAccel/verbose 2
# /gate/output/sinoAccel/RawOutputEnable false
# /gate/output/sinoAccel/setFileName UniCyl
# /gate/output/sinoAccel/TruesOnly true
# /gate/output/sinoAccel/RadialBins 336
# /gate/output/sinoAccel/StoreDelayeds
# /gate/output/sinoAccel/StoreScatters
# /gate/output/sinoAccel/setInputDataName finalCoinc
# /gate/output/sinoAccel/describe
/gate/output/ecat7/disable
# /gate/output/ecat7/verbose 2
# /gate/output/ecat7/setFileName UniCyl
# /gate/output/ecat7/describe
# /gate/output/ecat7/mashing 1
# /gate/output/ecat7/span 11
# /gate/output/ecat7/maxringdiff 27
# /gate/output/ecat7/system 1080
# /gate/output/ecat7/IsotopeCode C-11
# /gate/output/ecat7/IsotopeHalflife 1224.0 second
# /gate/output/ecat7/IsotopeBranchingFraction 1.0
# /gate/output/ecat7/describe
#-------------------oooooOOOOO00000OOOOOooooo---------------------#
# #
# D E F I N I T I O N O F #
# Y O U R V E R B O S I T Y L E V E L #
# #
#-------------------oooooOOOOO00000OOOOOooooo---------------------#
/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
#-------------------oooooOOOOO00000OOOOOooooo---------------------#
# #
# D E F I N I T I O N O F #
# Y O U R S O U R C E S #
# #
#-------------------oooooOOOOO00000OOOOOooooo---------------------#
#
# F A S T P R O D U C T I O N M O D U L E F O R F - 1 8 I S O T O P E S
#
/gate/source/addSource F18CylinderSource
/gate/source/F18CylinderSource/setActivity 37000. becquerel
/gate/source/F18CylinderSource/gps/particle e+
/gate/source/F18CylinderSource/setForcedUnstableFlag true
/gate/source/F18CylinderSource/setForcedHalfLife 6586 s
/gate/source/F18CylinderSource/gps/energytype Fluor18
/gate/source/F18CylinderSource/gps/type Volume
/gate/source/F18CylinderSource/gps/shape Cylinder
/gate/source/F18CylinderSource/gps/radius 0.5 mm
/gate/source/F18CylinderSource/gps/halfz 0.5 mm
/gate/source/F18CylinderSource/gps/angtype iso
/gate/source/F18CylinderSource/gps/centre 0.0 0.0 0.0 mm
/gate/source/F18CylinderSource/gps/confine phantom
/gate/source/list
#-------------------oooooOOOOO00000OOOOOooooo---------------------#
# #
# S T A R T A C Q U I S I T I O N #
# T I M E P A R A M E T E R S #
# #
#-------------------oooooOOOOO00000OOOOOooooo---------------------#
#
# EXPERIMENT
#
/gate/application/setTimeSlice 10.0 s
/gate/application/setTimeStart 0. s
/gate/application/setTimeStop 10.0 s
#
# RUN THE SIMULATION!
#
/gate/application/startDAQ
#
# E X I T
#
exit
More information about the Gate-users
mailing list