[Gate-users] LOW_VALUE_SENSITIVITY_PET
LORENA VAZQUEZ CANELAS
lorena.vazquez at usal.es
Mon Mar 9 14:29:14 CET 2015
Hello all!
I am trying to calculate the value of the sensitivity of PET, but the
problem is that it gives me a very low value and I don't know that I can be
doing wrong. I should get a value of 0.95 % = 9.5 but all it gives me is
0.6 %!!
The geometry have it well implemented .
I put a twogamma source , so that the simulation is faster. I also changed
the value of the cut, the threshold and dead-time and gives me the same
thing: S = 0.6 %
Send two photos of the PET manufacturing parameters I want to simulate.
Any ideas???
Thank you
I send you the code if you see anything wrong:
/gate/geometry/setMaterialDatabase
/home/lri/opengate-public/GateMaterials.db
# WORLD:
/gate/world/geometry/setXLength 400. cm
/gate/world/geometry/setYLength 400. cm
/gate/world/geometry/setZLength 400. cm
#-------------------oooooOOOOO00000OOOOOooooo---------------------#
# #
# DEFINITION AND DESCRITION OF YOUR PET DEVICE #
# #
#-------------------oooooOOOOO00000OOOOOooooo---------------------#
# 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 47.3 cm
/gate/cylindricalPET/geometry/setRmin 42.1 cm
/gate/cylindricalPET/geometry/setHeight 21.4 cm
/gate/cylindricalPET/setMaterial Air
#/gate/cylindricalPET/vis/forceWireframe 1
/gate/cylindricalPET/vis/setColor white
# HEAD (rsector)
/gate/cylindricalPET/daughters/name head
/gate/cylindricalPET/daughters/insert box
/gate/head/placement/setTranslation 42.1 0 0 cm
/gate/head/geometry/setXLength 2 cm
/gate/head/geometry/setYLength 5.2 cm
/gate/head/geometry/setZLength 20.8 cm
/gate/head/setMaterial Air
/gate/head/vis/setVisible 1
/gate/head/vis/setColor red
# BLOCK (module)
/gate/head/daughters/name block
/gate/head/daughters/insert box
/gate/block/placement/setTranslation 0 0 0 cm
/gate/block/geometry/setXLength 2 cm
/gate/block/geometry/setYLength 5.2 cm
/gate/block/geometry/setZLength 5.2 cm
/gate/block/setMaterial Air
/gate/block/vis/setVisible 1
/gate/block/vis/setColor blue
# CRYSTAL
/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 2 cm
/gate/crystal/geometry/setYLength 0.4 cm
/gate/crystal/geometry/setZLength 0.4 cm
/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 2 cm
/gate/LSO/geometry/setYLength 0.4 cm
/gate/LSO/geometry/setZLength 0.4 cm
/gate/LSO/setMaterial LSO
/gate/LSO/vis/setVisible 1
/gate/LSO/vis/setColor yellow
# REPEAT CRYSTAL
/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 cm
# REPEAT BLOCK
/gate/block/repeaters/insert cubicArray
/gate/block/cubicArray/setRepeatNumberX 1
/gate/block/cubicArray/setRepeatNumberY 1
/gate/block/cubicArray/setRepeatNumberZ 4
/gate/block/cubicArray/setRepeatVector 0 0 5.2 cm
# REPEAT HEAD
/gate/head/repeaters/insert ring
/gate/head/ring/setRepeatNumber 48
# ATTACH SYSTEM
/gate/systems/cylindricalPET/rsector/attach head
/gate/systems/cylindricalPET/module/attach block
/gate/systems/cylindricalPET/crystal/attach crystal
/gate/systems/cylindricalPET/layer0/attach LSO
# ATTACH CRYSTAL SD
/gate/LSO/attachCrystalSD
/gate/systems/cylindricalPET/describe
#/gate/cylindricalPET/moves/insert orbiting
#/gate/cylindricalPET/orbiting/setSpeed 30 deg/s
#/gate/cylindricalPET/orbiting/setPoint1 0 0 0 cm
#/gate/cylindricalPET/orbiting/setPoint2 0 0 1 cm
#====================================================================
# PHANTOMS
#====================================================================
# CYLINDER 5
/gate/world/daughters/name phantom5
/gate/world/daughters/insert cylinder
/gate/phantom5/setMaterial Air
/gate/phantom5/vis/setColor grey
/gate/phantom5/geometry/setRmax 0.955 cm
/gate/phantom5/geometry/setRmin 0.83 cm
/gate/phantom5/geometry/setHeight 70 cm
/gate/phantom5/placement/setTranslation 0 0 0 mm
#/gate/phantom5/attachPhantomSD
# CYLINDER 4
/gate/phantom5/daughters/name phantom4
/gate/phantom5/daughters/insert cylinder
/gate/phantom4/setMaterial Air
/gate/phantom4/vis/setColor blue
/gate/phantom4/geometry/setRmax 0.795 cm
/gate/phantom4/geometry/setRmin 0.67 cm
/gate/phantom4/geometry/setHeight 70 cm
/gate/phantom4/placement/setTranslation 0 0 0 mm
#/gate/phantom4/attachPhantomSD
# CYLINDER 3
/gate/phantom4/daughters/name phantom3
/gate/phantom4/daughters/insert cylinder
/gate/phantom3/setMaterial Air
/gate/phantom3/vis/setColor yellow
/gate/phantom3/geometry/setRmax 0.635 cm
/gate/phantom3/geometry/setRmin 0.51 cm
/gate/phantom3/geometry/setHeight 70 cm
/gate/phantom3/placement/setTranslation 0 0 0 mm
#/gate/phantom3/attachPhantomSD
# CYLINDER 2
/gate/phantom3/daughters/name phantom2
/gate/phantom3/daughters/insert cylinder
/gate/phantom2/setMaterial Air
/gate/phantom2/vis/setColor green
/gate/phantom2/geometry/setRmax 0.475 cm
/gate/phantom2/geometry/setRmin 0.35 cm
/gate/phantom2/geometry/setHeight 70 cm
/gate/phantom2/placement/setTranslation 0 0 0 mm
#/gate/phantom2/attachPhantomSD
# CYLINDER 1
/gate/phantom2/daughters/name phantom1
/gate/phantom2/daughters/insert cylinder
/gate/phantom1/setMaterial Air
/gate/phantom1/vis/setColor red
/gate/phantom1/geometry/setRmax 0.32 cm
/gate/phantom1/geometry/setRmin 0.195 cm
/gate/phantom1/geometry/setHeight 70 cm
/gate/phantom1/placement/setTranslation 0 0 0 mm
#/gate/phantom1/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
#=====================================================
# CUTS
#=====================================================
/gate/physics/Gamma/SetCutInRegion LSO 1 cm
/gate/physics/Electron/SetCutInRegion LSO 1 cm
/gate/physics/Positron/SetCutInRegion LSO 1 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
#========================================================================
# DIGITALIZADOR: DETECTOR ELECTRONIC RESPONSE
#========================================================================
# ADDER:
/gate/digitizer/convertor/verbose 0
/gate/digitizer/Singles/insert adder
/gate/digitizer/Singles/adder/verbose 0
# READOUT:
gate/digitizer/Singles/insert readout
/gate/digitizer/Singles/readout/setDepth 1
# BLURRING:
/gate/digitizer/Singles/insert blurring
/gate/digitizer/Singles/blurring/setResolution 0.12
/gate/digitizer/Singles/blurring/setCrystalEnergyOfReference 511. keV
# THRESHOLDER/UPHOLDER:
/gate/digitizer/Singles/insert thresholder
/gate/digitizer/Singles/thresholder/setThreshold 435. keV
/gate/digitizer/Singles/insert upholder
/gate/digitizer/Singles/upholder/setUphold 700. keV
#====================================================================================
# COINCIDENCE SORTER
#====================================================================================
/gate/digitizer/Coincidences/setWindow 4.1 ns
/gate/digitizer/name delay
/gate/digitizer/insert coincidenceSorter
/gate/digitizer/delay/setWindow 4.1 ns
/gate/digitizer/delay/setOffset 500. ns
/gate/systems/cylindricalPET/describe
#=====================================================
# PARTICLE SOURCE
#=====================================================
/gate/source/addSource twogamma
/gate/source/twogamma/setActivity 1800 becquerel
/gate/source/twogamma/setType backtoback
/gate/source/twogamma/gps/centre 0. 0. 0. cm
/gate/source/twogamma/gps/particle gamma
/gate/source/twogamma/gps/energytype Mono
/gate/source/twogamma/gps/monoenergy 0.511 MeV
/gate/source/twogamma/gps/shape Cylinder
/gate/source/twogamma/gps/type Volume
/gate/source/twogamma/gps/radius 0.3 cm
/gate/source/twogamma/gps/halfz 35 cm
/gate/source/twogamma/gps/angtype iso
/gate/source/twogamma/gps/mintheta 0. deg
/gate/source/twogamma/gps/maxtheta 180. deg
/gate/source/twogamma/gps/minphi 0. deg
/gate/source/twogamma/gps/maxphi 360. deg
/gate/source/list
#=====================================================
# DATA OUTPUT
#=====================================================
/gate/output/root/enable
/gate/output/root/setFileName SALIDA_ROOT
/gate/output/root/setRootHitFlag 0
/gate/output/root/setRootSinglesFlag 0
/gate/output/root/setRootCoincidencesFlag 1
/gate/output/root/setRootdelayFlag 0
#====================================================
# RANDOM ENGINE AND SEED
#====================================================
/gate/random/setEngineName JamesRandom
/gate/random/setEngineSeed default
/gate/random/verbose 0
#=====================================================
# MEASUREMENT SETTINGS
#=====================================================
/gate/application/setTimeSlice 1 s
/gate/application/setTimeStart 0 s
/gate/application/setTimeStop 1 s
#=====================================================
# VISUALISATION
#=====================================================
#/vis/disable
/vis/open DAWNFILE
/vis/viewer/set/viewpointThetaPhi 60 60
/vis/viewer/zoom 8.0
/vis/drawVolume
/tracking/storeTrajectory 1
/vis/scene/add/trajectories
/vis/scene/endOfEventAction accumulate
/vis/viewer/update
/vis/viewer/refresh
/control/execute Verbose.mac
/gate/application/startDAQ
--
LORENA VÁZQUEZ CANELAS
--
LRI - Datación (Despacho Técnicos)
Calle Espejo, s/n, 37008 - Salamanca (Spain)
TEL: +34 923 294930 ext:4932 email: lorena.vazquez at usal.es
web: http://campus.usal.es/~lri/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20150309/79f442f5/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: PET1.jpg
Type: image/jpeg
Size: 57233 bytes
Desc: not available
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20150309/79f442f5/attachment-0002.jpg>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: PET2.jpg
Type: image/jpeg
Size: 93446 bytes
Desc: not available
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20150309/79f442f5/attachment-0003.jpg>
More information about the Gate-users
mailing list