[Gate-users] about GE Discovery MI sensitive

袁波 gogenbyts at sina.com
Tue Jul 9 10:08:58 CEST 2019 

hi , everyone :
          when  i simulate  GE discovery  MI PET/CT  https://aapm.onlinelibrary.wiley.com/doi/abs/10.1002/mp.13576
 i  simulation as below, but the result is 9.8cps/kbq, but  is 20.8cps/kbq.so may you tell me where i did wrong ?by the way , i got the SF（Scattering fraction） is 39.7%，which is right。thank you for your time。by the way, i got sensitive by this formula ：  number of true event / radioactivity *timein here the time is 1s,radioactivity is 10000000bq


here is gate code
/vis/disable

/gate/geometry/setMaterialDatabase 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

/control/execute camera.mac

/control/execute phantom.mac

/control/execute physics.mac



#    I N I T I A L I Z E 
/gate/run/initialize

/control/execute digitizer.mac

#    S O U R C E
/control/execute sources.mac



#    V E R B O S I T Y
#/gate/verbose Physic    0
#/gate/verbose Cuts      0
#/gate/verbose Actor     0
#/gate/verbose SD        0
#/gate/verbose Actions   0
#/gate/verbose Step      0
#/gate/verbose Error     0
#/gate/verbose Warning   0
#/gate/verbose Output    0
#/gate/verbose Core      0

/run/verbose 0
/event/verbose 0
/tracking/verbose 0

#    O U T P U T 

#root
/gate/output/root/enable
/gate/output/root/setFileName    benchmarkPET
/gate/output/root/setRootHitFlag 0
/gate/output/root/setRootSinglesFlag 0
/gate/output/root/setRootCoincidencesFlag 1

#lmf
#/gate/output/lmf/enable
#/gate/output/lmf/setFileName              10min
#/gate/output/lmf/setDetectorIDBool        1
#/gate/output/lmf/setEnergyBool            0
#/gate/output/lmf/setGantryAxialPosBool    0
#/gate/output/lmf/setGantryAngularPosBool  0
#/gate/output/lmf/setSourcePosBool         0
#/gate/output/lmf/setNeighbourBool         0
#/gate/output/lmf/setNeighbourhoodOrder    0
#/gate/output/lmf/setCoincidenceBool       0
#/gate/output/lmf/setGateDigiBool          1 
#/gate/output/lmf/setComptonBool           0 
#/gate/output/lmf/setSourceIDBool          0
#/gate/output/lmf/setSourceXYZPosBool      1
#/gate/output/lmf/setGlobalXYZPosBool      0
#/gate/output/lmf/setEventIDBool           0
#/gate/output/lmf/setRunIDBool             0

#ascii
#/gate/output/ascii/enable
#/gate/output/ascii/setFileName                  10min
#/gate/output/ascii/setOutFileHitsFlag           0
#/gate/output/ascii/setOutFileSinglesFlag        0
#/gate/output/ascii/setOutFileCoincidencesFlag   1
#/gate/output/ascii/setCoincidenceMask 0 0 0 1 1 1 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 1 1 1 0 0 0 0 1 0 0 0 0 0 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      1. s

/gate/application/startDAQ

# C Y L I N D R I C A L 
/gate/world/daughters/name cylindricalPET 
/gate/world/daughters/insert cylinder 
/gate/cylindricalPET/setMaterial Air 
/gate/cylindricalPET/geometry/setRmax 380 mm 
/gate/cylindricalPET/geometry/setRmin 350 mm 
/gate/cylindricalPET/geometry/setHeight 250 mm 
/gate/cylindricalPET/vis/forceWireframe

# R S E C T O R 
/gate/cylindricalPET/daughters/name rsector 
/gate/cylindricalPET/daughters/insert box 
/gate/rsector/placement/setTranslation 363 0 0 mm 
/gate/rsector/geometry/setXLength 25. mm 
/gate/rsector/geometry/setYLength 66.2 mm 
/gate/rsector/geometry/setZLength 250.0 mm 
/gate/rsector/setMaterial Air 
/gate/rsector/vis/forceWireframe

# M O D U L E 
/gate/rsector/daughters/name module 
/gate/rsector/daughters/insert box 
/gate/module/geometry/setXLength 25. mm 
/gate/module/geometry/setYLength 66.2 mm 
/gate/module/geometry/setZLength 49.3 mm 
/gate/module/setMaterial Air 
/gate/module/vis/forceWireframe 
/gate/module/vis/setColor gray

# C R Y S T A L 
/gate/module/daughters/name crystal 
/gate/module/daughters/insert box 
/gate/crystal/geometry/setXLength 25. mm 
/gate/crystal/geometry/setYLength 3.95 mm 
/gate/crystal/geometry/setZLength 5.3 mm 
/gate/crystal/setMaterial Air 
/gate/crystal/vis/forceWireframe 
/gate/crystal/vis/setColor magenta

# L A Y E R 
/gate/crystal/daughters/name LYSO 
/gate/crystal/daughters/insert box 
/gate/LYSO/geometry/setXLength 25. mm 
/gate/LYSO/geometry/setYLength 3.95 mm 
/gate/LYSO/geometry/setZLength 5.3 mm 
/gate/LYSO/placement/setTranslation 0 0 0 mm 
/gate/LYSO/setMaterial LYSO 
/gate/LYSO/vis/setColor yellow

# 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 16 
/gate/crystal/cubicArray/setRepeatNumberZ 9 
/gate/crystal/cubicArray/setRepeatVector 0.0 4.15 5.5 mm

# R E P E A T M O D U L E 
/gate/module/repeaters/insert cubicArray 
/gate/module/cubicArray/setRepeatNumberZ 5 
/gate/module/cubicArray/setRepeatVector 0. 0. 50.0 mm

# R E P E A T R S E C T O R 
/gate/rsector/repeaters/insert ring 
/gate/rsector/ring/setRepeatNumber 34

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

# A T T A C H L A Y E R SD 
/gate/LYSO/attachCrystalSD 


### Sum all pulses
/gate/digitizer/Singles/insert adder
/gate/digitizer/Singles/insert readout
/gate/digitizer/Singles/readout/setDepth 2

### Model LYSO energy resolution
/gate/digitizer/Singles/insert blurring
/gate/digitizer/Singles/blurring/setResolution 0.12
/gate/digitizer/Singles/blurring/setEnergyOfReference 511 keV

### Model the 425-650 keV energy window
/gate/digitizer/Singles/insert thresholder
/gate/digitizer/Singles/thresholder/setThreshold 425 keV
/gate/digitizer/Singles/insert upholder
/gate/digitizer/Singles/upholder/setUphold 650 keV


### Model the 4.9 ns coincidence time window and the minimum required module

/gate/digitizer/Coincidences/setWindow 4.9 ns
/gate/digitizer/Coincidences/MultiplesPolicy keepIfAllAreGoods
#    P H A N T O M
/gate/world/daughters/name NEMACylinder
/gate/world/daughters/insert cylinder
/gate/NEMACylinder/setMaterial Water
/gate/NEMACylinder/geometry/setRmax   2 mm
/gate/NEMACylinder/geometry/setHeight 700. mm

/gate/NEMACylinder/vis/forceSolid
/gate/NEMACylinder/vis/setColor red


/gate/NEMACylinder/attachPhantomSD

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

/gate/source/addSource F18LineSource
/gate/source/F18LineSource/setType backtoback
/gate/source/F18LineSource/gps/particle gamma
/gate/source/F18LineSource/gps/energytype Mono
/gate/source/F18LineSource/gps/monoenergy 511 keV
/gate/source/F18LineSource/setActivity 10000000. becquerel
/gate/source/F18LineSource/setForcedUnstableFlag true
/gate/source/F18LineSource/setForcedHalfLife 6586.2 s
/gate/source/F18LineSource/gps/angtype iso
/gate/source/F18LineSource/gps/type Volume
/gate/source/F18LineSource/gps/shape Cylinder
/gate/source/F18LineSource/gps/radius 2 mm
/gate/source/F18LineSource/gps/halfz 35.0 cm
/gate/source/F18LineSource/gps/centre 0.0 0.0 0.0 cm



/gate/source/list

--------------------------------

From Yuan Bo
School of Physical Engineering ,Zhengzhou University
No.100 Science Avenue,Zhengzhou City,HeNan Province,P.R.China.
来自袁波
郑州大学物理工程学院C楼4层核医学仪器实验室
河南省郑州市科学大道100号

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20190709/2fcf07f5/attachment.html>

More information about the Gate-users mailing list