<div dir="ltr"><div class="gmail_default" style="font-family:comic sans ms,sans-serif;font-size:small"> 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 </div> -- <div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><pre>LORENA VÁZQUEZ CANELAS -- LRI - Datación (Despacho Técnicos) Calle Espejo, s/n, 37008 - Salamanca (Spain) TEL: +34 923 294930 ext:4932 email: <a href="mailto:lorena.vazquez@usal.es" target="_blank">lorena.vazquez@usal.es</a> web: <a href="http://campus.usal.es/%7Elri/" target="_blank">http://campus.usal.es/~lri/</a></pre></div></div></div></div></div></div></div></div></div></div> </div>