<p>hi , everyone :</p><p><br></p><p>          when  i simulate  GE discovery  MI PET/CT  <a target="_blank" href="https://aapm.onlinelibrary.wiley.com/doi/abs/10.1002/mp.13576">https://aapm.onlinelibrary.wiley.com/doi/abs/10.1002/mp.13576</a></p><p><br></p><p><font size="4"> i  simulation as below, but the result is 9.8cps/kbq, but  is 20.8cps/kbq.</font></p><div><font size="4">so may you tell me where i did wrong ?</font></div><div><font size="4">by the way , i got the SF（Scattering fraction</font><span style="font-size: large;">） is 39.7%，which is right。</span></div><div><span style="font-size: large;">thank you for your time。</span></div><div><span style="font-size: large;"></span>by the way, i got sensitive by this formula ：  number of true event / radioactivity *time</div><div>in here the time is 1s,radioactivity is 10000000bq<span style="font-size: large;"></span></div><div><span style="font-size: large;"><br></span></div><div><span style="font-size: large;"><br></span></div><div><span style="font-size: large;"><br></span></div><div><div><span style="font-size: large;">here is gate code</span></div><div><span style="font-size: large;"><br></span></div>/vis/disable<br><br>/gate/geometry/setMaterialDatabase GateMaterials.db<br><br>#     W O R L D<br>/gate/world/geometry/setXLength 150. cm<br>/gate/world/geometry/setYLength 150. cm<br>/gate/world/geometry/setZLength 150. cm<br><br>/control/execute camera.mac<br><br>/control/execute phantom.mac<br><br>/control/execute physics.mac<br><br><br><br>#    I N I T I A L I Z E <br>/gate/run/initialize<br><br>/control/execute digitizer.mac<br><br>#    S O U R C E<br>/control/execute sources.mac<br><br><br><br>#    V E R B O S I T Y<br>#/gate/verbose Physic    0<br>#/gate/verbose Cuts      0<br>#/gate/verbose Actor     0<br>#/gate/verbose SD        0<br>#/gate/verbose Actions   0<br>#/gate/verbose Step      0<br>#/gate/verbose Error     0<br>#/gate/verbose Warning   0<br>#/gate/verbose Output    0<br>#/gate/verbose Core      0<br><br>/run/verbose 0<br>/event/verbose 0<br>/tracking/verbose 0<br><br>#    O U T P U T <br><br>#root<br>/gate/output/root/enable<br>/gate/output/root/setFileName    benchmarkPET<br>/gate/output/root/setRootHitFlag 0<br>/gate/output/root/setRootSinglesFlag 0<br>/gate/output/root/setRootCoincidencesFlag 1<br><br>#lmf<br>#/gate/output/lmf/enable<br>#/gate/output/lmf/setFileName              10min<br>#/gate/output/lmf/setDetectorIDBool        1<br>#/gate/output/lmf/setEnergyBool            0<br>#/gate/output/lmf/setGantryAxialPosBool    0<br>#/gate/output/lmf/setGantryAngularPosBool  0<br>#/gate/output/lmf/setSourcePosBool         0<br>#/gate/output/lmf/setNeighbourBool         0<br>#/gate/output/lmf/setNeighbourhoodOrder    0<br>#/gate/output/lmf/setCoincidenceBool       0<br>#/gate/output/lmf/setGateDigiBool          1 <br>#/gate/output/lmf/setComptonBool           0 <br>#/gate/output/lmf/setSourceIDBool          0<br>#/gate/output/lmf/setSourceXYZPosBool      1<br>#/gate/output/lmf/setGlobalXYZPosBool      0<br>#/gate/output/lmf/setEventIDBool           0<br>#/gate/output/lmf/setRunIDBool             0<br><br>#ascii<br>#/gate/output/ascii/enable<br>#/gate/output/ascii/setFileName                  10min<br>#/gate/output/ascii/setOutFileHitsFlag           0<br>#/gate/output/ascii/setOutFileSinglesFlag        0<br>#/gate/output/ascii/setOutFileCoincidencesFlag   1<br>#/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 <br><br><br>#   R A N D O M<br># JamesRandom Ranlux64 MersenneTwister<br>/gate/random/setEngineName MersenneTwister<br>#/gate/random/setEngineSeed default<br>#/gate/random/setEngineSeed auto<br>/gate/random/setEngineSeed 123456789<br>#/gate/random/resetEngineFrom fileName<br>/gate/random/verbose 1<br><br>#    S T A R T <br>/gate/application/setTimeSlice     1. s<br>/gate/application/setTimeStart     0. s<br>/gate/application/setTimeStop      1. s<br><br><div>/gate/application/startDAQ</div><div><br></div><div><br></div><div># C Y L I N D R I C A L <br>/gate/world/daughters/name cylindricalPET <br>/gate/world/daughters/insert cylinder <br>/gate/cylindricalPET/setMaterial Air <br>/gate/cylindricalPET/geometry/setRmax 380 mm <br>/gate/cylindricalPET/geometry/setRmin 350 mm <br>/gate/cylindricalPET/geometry/setHeight 250 mm <br>/gate/cylindricalPET/vis/forceWireframe<br><br># R S E C T O R <br>/gate/cylindricalPET/daughters/name rsector <br>/gate/cylindricalPET/daughters/insert box <br>/gate/rsector/placement/setTranslation 363 0 0 mm <br>/gate/rsector/geometry/setXLength 25. mm <br>/gate/rsector/geometry/setYLength 66.2 mm <br>/gate/rsector/geometry/setZLength 250.0 mm <br>/gate/rsector/setMaterial Air <br>/gate/rsector/vis/forceWireframe<br><br># M O D U L E <br>/gate/rsector/daughters/name module <br>/gate/rsector/daughters/insert box <br>/gate/module/geometry/setXLength 25. mm <br>/gate/module/geometry/setYLength 66.2 mm <br>/gate/module/geometry/setZLength 49.3 mm <br>/gate/module/setMaterial Air <br>/gate/module/vis/forceWireframe <br>/gate/module/vis/setColor gray<br><br># C R Y S T A L <br>/gate/module/daughters/name crystal <br>/gate/module/daughters/insert box <br>/gate/crystal/geometry/setXLength 25. mm <br>/gate/crystal/geometry/setYLength 3.95 mm <br>/gate/crystal/geometry/setZLength 5.3 mm <br>/gate/crystal/setMaterial Air <br>/gate/crystal/vis/forceWireframe <br>/gate/crystal/vis/setColor magenta<br><br># L A Y E R <br>/gate/crystal/daughters/name LYSO <br>/gate/crystal/daughters/insert box <br>/gate/LYSO/geometry/setXLength 25. mm <br>/gate/LYSO/geometry/setYLength 3.95 mm <br>/gate/LYSO/geometry/setZLength 5.3 mm <br>/gate/LYSO/placement/setTranslation 0 0 0 mm <br>/gate/LYSO/setMaterial LYSO <br>/gate/LYSO/vis/setColor yellow<br><br># R E P E A T C R Y S T A L <br>/gate/crystal/repeaters/insert cubicArray <br>/gate/crystal/cubicArray/setRepeatNumberX 1 <br>/gate/crystal/cubicArray/setRepeatNumberY 16 <br>/gate/crystal/cubicArray/setRepeatNumberZ 9 <br>/gate/crystal/cubicArray/setRepeatVector 0.0 4.15 5.5 mm<br><br># R E P E A T M O D U L E <br>/gate/module/repeaters/insert cubicArray <br>/gate/module/cubicArray/setRepeatNumberZ 5 <br>/gate/module/cubicArray/setRepeatVector 0. 0. 50.0 mm<br><br># R E P E A T R S E C T O R <br>/gate/rsector/repeaters/insert ring <br>/gate/rsector/ring/setRepeatNumber 34<br><br># A T T A C H S Y S T E M <br>/gate/systems/cylindricalPET/rsector/attach rsector <br>/gate/systems/cylindricalPET/module/attach module <br>/gate/systems/cylindricalPET/crystal/attach crystal <br>/gate/systems/cylindricalPET/layer0/attach LYSO<br><br># A T T A C H L A Y E R SD <br>/gate/LYSO/attachCrystalSD <br><div><br></div><div><br></div><div>### Sum all pulses<br>/gate/digitizer/Singles/insert adder<br>/gate/digitizer/Singles/insert readout<br>/gate/digitizer/Singles/readout/setDepth 2<br><br>### Model LYSO energy resolution<br>/gate/digitizer/Singles/insert blurring<br>/gate/digitizer/Singles/blurring/setResolution 0.12<br>/gate/digitizer/Singles/blurring/setEnergyOfReference 511 keV<br><br>### Model the 425-650 keV energy window<br>/gate/digitizer/Singles/insert thresholder<br>/gate/digitizer/Singles/thresholder/setThreshold 425 keV<br>/gate/digitizer/Singles/insert upholder<br>/gate/digitizer/Singles/upholder/setUphold 650 keV<br><br><br>### Model the 4.9 ns coincidence time window and the minimum required module<br><br>/gate/digitizer/Coincidences/setWindow 4.9 ns<br><div>/gate/digitizer/Coincidences/MultiplesPolicy keepIfAllAreGoods</div><div><br></div><div>#    P H A N T O M<br>/gate/world/daughters/name NEMACylinder<br>/gate/world/daughters/insert cylinder<br>/gate/NEMACylinder/setMaterial Water<br>/gate/NEMACylinder/geometry/setRmax   2 mm<br>/gate/NEMACylinder/geometry/setHeight 700. mm<br><br>/gate/NEMACylinder/vis/forceSolid<br>/gate/NEMACylinder/vis/setColor red<br><br><br><div>/gate/NEMACylinder/attachPhantomSD</div><div><br></div><div><br></div><div>#################<br>#               #<br>#  S O U R C E  #<br>#               #<br>#################<br><br>/gate/source/addSource F18LineSource<br>/gate/source/F18LineSource/setType backtoback<br>/gate/source/F18LineSource/gps/particle gamma<br>/gate/source/F18LineSource/gps/energytype Mono<br>/gate/source/F18LineSource/gps/monoenergy 511 keV<br>/gate/source/F18LineSource/setActivity 10000000. becquerel<br>/gate/source/F18LineSource/setForcedUnstableFlag true<br>/gate/source/F18LineSource/setForcedHalfLife 6586.2 s<br>/gate/source/F18LineSource/gps/angtype iso<br>/gate/source/F18LineSource/gps/type Volume<br>/gate/source/F18LineSource/gps/shape Cylinder<br>/gate/source/F18LineSource/gps/radius 2 mm<br>/gate/source/F18LineSource/gps/halfz 35.0 cm<br>/gate/source/F18LineSource/gps/centre 0.0 0.0 0.0 cm<br><br><br><br>/gate/source/list<br><br></div></div></div></div></div><div id="">--------------------------------<br></div><div id="">
<div style="text-align: left;"><font size="2" face="Times New Roman"><span style="background-color: rgb(255, 255, 255);"><b>From Yuan Bo</b></span></font></div>
<div style="text-align: left;"><span style="text-align: justify; background-color: rgb(242, 242, 242);"><b><font size="2" face="Times New Roman">School of Physical Engineering </font></b></span><font size="2" face="Times New Roman"><span style="background-color: rgb(255, 255, 255);"><b>,Zhengzhou University</b></span></font></div>
<div style="text-align: left;"><font size="2" face="Times New Roman"><span style="background-color: rgb(255, 255, 255);"><b>No.100 Science Avenue,Zhengzhou City,HeNan Province,P.R.China.</b></span></font></div>
<div style="text-align: left;"><font size="2" face="Times New Roman"><span style="background-color: rgb(255, 255, 255);"><b>来自袁波</b></span></font></div>
<div style="text-align: left;"><font size="2" face="Times New Roman"><span style="background-color: rgb(255, 255, 255);"><b>郑州大学物理工程学院C楼4层核医学仪器实验室</b></span></font></div>
<div style="text-align: left;"><font size="2" face="Times New Roman"><span style="background-color: rgb(255, 255, 255);"><b>河南省郑州市科学大道100号</b></span></font></div>
<br></div>