<html><head><meta http-equiv=Content-Type content="text/html; charset=utf-8"><META name="Author" content="Novell GroupWise WebAccess"></head><body style='font-family: Tahoma, sans-serif; font-size: 13px; '>Hello again,<br><br>I have determined my problem with the coincidence window command. The /gate/digitizer/Coincidences branch is not available until the proper system has been defined. I was attempting to simulate a two-plate PET scan, but had my system defined with systemType scanner. I just created a cylinder volume and used a ring repeater, but repeated the ring twice. Once I changed my systemType to cylindricalPET and attached the daughter volumes as required, I was able to set the coincidence window.<br><br>I'm still curious about my crosstalk question I posted previously, though. <br><br>Now I have a different issue. When reading the root files, the leaves in the Hits and Singles branches are all filled, but my leaves for my Coincidences and delay branches are all empty. I can't figure out why. When I run the benchPET macro and PET example macro, it's all there.<br><br>I am simulating two plates of LYSO crystals with a point source of backtoback 511keV gamma rays. I made the emission angles of the source to be concentrated over the plates. I know this works because I'm getting plenty of entries in my Hits and Singles output branches. It's just the Coincidence and delay that I'm having trouble with. I made my phantom a sphere of water for now. Here is my code so far:<br><br><br>### databases<br>/gate/geometry/setMaterialDatabase /home/djvanel/GateMaterials.db<br>/gate/physics/addPhysicsList QGSP_BERT_EMV<br><br>### world<br>/gate/world/geometry/setXLength 100. cm<br>/gate/world/geometry/setYLength 100. cm<br>/gate/world/geometry/setZLength 100. cm<br>/gate/world/setMaterial Air<br><br>### visual<br>/vis/open OGLSX<br>/vis/drawVolume<br>/vis/viewer/set/viewpointThetaPhi 15 30<br>/vis/scene/add/axes<br>#/tracking/storeTrajectory 1<br>#/vis/scene/add/trajectories<br>/vis/scene/endOfEventAction accumulate<br>#/vis/viewer/zoom 1.5<br>/gate/world/vis/setVisible false<br><br>### cylinder<br>/gate/world/daughters/name cylindricalPET<br>#/gate/world/daughters/systemType cylindricalPET<br>/gate/world/daughters/insert cylinder<br>/gate/cylindricalPET/setMaterial Air<br>/gate/cylindricalPET/geometry/setRmin 6. cm<br>/gate/cylindricalPET/geometry/setRmax 15. cm<br>/gate/cylindricalPET/geometry/setHeight 25. cm<br>/gate/cylindricalPET/vis/forceWireframe<br>/gate/cylindricalPET/vis/setColor white<br>/gate/geometry/rebuild<br><br>### board<br>/gate/cylindricalPET/daughters/name board<br>/gate/cylindricalPET/daughters/insert box<br>/gate/board/geometry/setXLength 2. cm<br>/gate/board/geometry/setYLength 20. cm<br>/gate/board/geometry/setZLength 20. cm<br>/gate/board/vis/forceWireframe<br>/gate/board/vis/setColor white<br>/gate/board/placement/setTranslation 10. 0. 0. cm<br>/gate/geometry/rebuild<br><br>### box1<br>/gate/board/daughters/name box1<br>/gate/board/daughters/insert box<br>/gate/box1/geometry/setXLength 2. cm<br>/gate/box1/geometry/setYLength 5. cm<br>/gate/box1/geometry/setZLength 5. cm<br>/gate/box1/setMaterial Air<br>/gate/box1/vis/forceWireframe<br>/gate/box1/vis/setColor red<br>/gate/geometry/rebuild<br><br>### pixel<br>/gate/box1/daughters/name pixel<br>/gate/box1/daughters/insert box<br>/gate/pixel/geometry/setXLength 2. cm<br>/gate/pixel/geometry/setYLength 1. cm<br>/gate/pixel/geometry/setZLength 1. cm<br>/gate/pixel/setMaterial Air<br>/gate/pixel/vis/forceWireframe<br>/gate/pixel/vis/setColor red<br>/gate/geometry/rebuild<br><br>### crystal<br>/gate/pixel/daughters/name crystal<br>/gate/pixel/daughters/insert box<br>/gate/crystal/geometry/setXLength 2. cm<br>/gate/crystal/geometry/setYLength 1. cm<br>/gate/crystal/geometry/setZLength 1. cm<br>/gate/crystal/setMaterial LYSO<br>/gate/crystal/vis/forceSolid<br>/gate/crystal/vis/setColor green<br>/gate/geometry/rebuild<br><br>### repeat pixel<br>/gate/pixel/repeaters/insert cubicArray<br>/gate/pixel/cubicArray/setRepeatNumberY 5<br>/gate/pixel/cubicArray/setRepeatNumberZ 5<br>/gate/pixel/cubicArray/setRepeatVector 0. 1. 1. cm<br>/gate/geometry/rebuild<br><br>### repeat box<br>/gate/box1/repeaters/insert cubicArray<br>/gate/box1/cubicArray/setRepeatNumberY 4<br>/gate/box1/cubicArray/setRepeatNumberZ 4<br>/gate/box1/cubicArray/setRepeatVector 0. 5. 5. cm<br>/gate/geometry/rebuild<br><br>### repeat board<br>/gate/board/repeaters/insert ring<br>/gate/board/ring/setRepeatNumber 2<br><br>### phantom<br>/gate/world/daughters/name phantom<br>/gate/world/daughters/insert sphere<br>/gate/phantom/geometry/setRmin 0. cm<br>/gate/phantom/geometry/setRmax 3. cm<br>/gate/phantom/setMaterial Water<br>/gate/phantom/vis/forceSolid<br>/gate/phantom/vis/setColor blue<br>/gate/geometry/rebuild<br><br>### system<br>/gate/systems/cylindricalPET/rsector/attach board<br>/gate/systems/cylindricalPET/module/attach box1<br>/gate/systems/cylindricalPET/crystal/attach pixel<br>/gate/systems/cylindricalPET/layer0/attach crystal<br>/gate/systems/cylindricalPET/describe<br><br>### attach layer SD<br>/gate/phantom/attachPhantomSD<br>/gate/crystal/attachCrystalSD<br><br><br>### physics<br>/gate/physics/addProcess PhotoElectric<br>/gate/physics/processes/PhotoElectric/setModel StandardModel<br><br>/gate/physics/addProcess Compton<br>/gate/physics/processes/Compton/setModel StandardModel<br><br>/gate/physics/addProcess RayleighScattering gamma<br>/gate/physics/processes/RayleighScattering/setModel PenelopeModel<br><br>/gate/physics/addProcess GammaConversion<br>/gate/physics/processes/GammaConversion/setModel StandardModel<br><br>/gate/physics/addProcess ElectronIonisation<br>/gate/physics/processes/ElectronIonisation/setModel StandardModel e-<br>/gate/physics/processes/ElectronIonisation/setModel StandardModel e+<br><br>/gate/physics/addProcess Bremsstrahlung<br>/gate/physics/processes/Bremsstrahlung/setModel StandardModel e-<br>/gate/physics/processes/Bremsstrahlung/setModel StandardModel e+<br><br>/gate/physics/addProcess PositronAnnihilation<br><br>/gate/physics/addProcess RadioactiveDecay<br><br>/gate/physics/addProcess eMultipleScattering<br>/gate/physics/processes/eMultipleScattering/setGeometricalStepLimiterType e- distanceToBoundary<br>/gate/physics/processes/eMultipleScattering/setGeometricalStepLimiterType e+ distanceToBoundary<br><br>/gate/physics/processList Enabled<br>/gate/physics/processList Initialized<br><br>### cut in regions<br>#/gate/physics/Gamma/SetCutInRegion world 10. cm<br>#/gate/physics/Gamma/SetCutInRegion phantom 5. mm<br>#/gate/physics/SetMaxStepSizeInRegion phantom 5. mm<br><br>### initialization<br>/gate/run/initialize<br><br>### digitizer<br>/gate/digitizer/Singles/insert adder<br>/gate/digitizer/Singles/insert readout<br>/gate/digitizer/Singles/readout/setPolicy TakeEnergyWinner<br>/gate/digitizer/Singles/readout/setDepth 1<br>/gate/digitizer/Singles/insert crystalblurring<br>/gate/digitizer/Singles/crystalblurring/setCrystalResolutionMin 0.15<br>/gate/digitizer/Singles/crystalblurring/setCrystalResolutionMax 0.35<br>/gate/digitizer/Singles/crystalblurring/setCrystalQE 0.9<br>/gate/digitizer/Singles/crystalblurring/setCrystalEnergyOfReference 511. keV<br>/gate/digitizer/Singles/insert thresholder<br>/gate/digitizer/Singles/thresholder/setThreshold 0. keV<br>/gate/digitizer/Singles/insert upholder<br>/gate/digitizer/Singles/upholder/setUphold 750. keV<br><br>/gate/digitizer/Coincidences/setWindow 50. ns<br>/gate/digitizer/Coincidences/minSectorDifference 2<br>/gate/digitizer/Coincidences/MultiplesPolicy takeWinnerOfGoods<br><br>/gate/digitizer/name delay<br>/gate/digitizer/insert coincidenceSorter<br>/gate/digitizer/delay/setWindow 50. ns<br>/gate/digitizer/delay/setOffset 500. ns<br>/gate/digitizer/delay/MultiplesPolicy takeWinnerOfGoods<br>/gate/digitizer/delay/setDepth 1<br><br>### source<br>/gate/source/addSource Gamma<br>/gate/source/Gamma/gps/centre 0. 0. 0. cm<br>/gate/source/Gamma/setActivity 100000. becquerel<br>/gate/source/Gamma/gps/particle gamma<br>/gate/source/Gamma/setType backtoback<br>/gate/source/Gamma/gps/energytype Mono<br>/gate/source/Gamma/gps/monoenergy 511. keV<br>/gate/source/Gamma/setForcedUnstableFlag true<br>/gate/source/Gamma/setForcedHalfLife 1000 s<br>/gate/source/Gamma/gps/angtype iso<br>/gate/source/Gamma/gps/mintheta 60. deg<br>/gate/source/Gamma/gps/maxtheta 120. deg<br>/gate/source/Gamma/gps/minphi 120. deg<br>/gate/source/Gamma/gps/maxphi 240. deg<br><br>### output<br>/gate/output/root/enable<br>/gate/output/root/setFileName dualplate<br>/gate/output/root/setRootHitFlag 1<br>/gate/output/root/setRootSinglesFlag 1<br>/gate/output/root/setRootNtupleFlag 1<br>/gate/output/root/setRootCoincidencesFlag 1<br>/gate/output/root/setRootdelayFlag 1<br><br>### random<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 0<br><br>### simulation length<br>/gate/application/setTimeSlice 1 s<br>/gate/application/setTimeStart 0. s<br>/gate/application/setTimeStop 10 s<br><br><br><br>Any help is appreciated! Thank you<br>Devin Van Elburg<br><br>---------------<br>Physics w/ c/ Biomedical Physics<br>Lakehead University<br><br></body></html>