<div dir="ltr"><br><div><br></div><div>is there anyway to model an ideal PET system where the (x,y,z) coordinates in the detector crystals have no variance? that is, each coincidence event contains two points where the connecting line passes through the source up to round-off error</div>
<div><br></div><div>i tried to do this with the following:</div>
<div><div># setup simple physics - turn off scattering - this is horribly inefficient since it now requires 2 simultaneous PE at 511 keV to generate coincidence, but i get about 5000 in 1M primaries</div><div>/gate/physics/addProcess PhotoElectric</div>
<div>/gate/physics/processes/PhotoElectric/setModel StandardModel</div><div>#/gate/physics/addProcess Compton<br></div>
<div>#/gate/physics/processes/Compton/setModel StandardModel</div><div>#/gate/physics/addProcess RayleighScattering<br></div><div>#/gate/physics/processes/RayleighScattering/setModel PenelopeModel</div><div>/gate/physics/addProcess ElectronIonisation<br>
</div><div>/gate/physics/processes/ElectronIonisation/setModel StandardModel e-</div><div>/gate/physics/processes/ElectronIonisation/setModel StandardModel e+</div></div>
<div><div>/gate/physics/processList Enabled</div><div>/gate/physics/processList Initialized</div></div><div><br></div><div style>#disable adder and blurring</div><div style><div>#/gate/digitizer/Singles/insert adder</div>
<div>/gate/digitizer/Singles/insert readout</div><div>/gate/digitizer/Singles/readout/setDepth 1</div><div>#/gate/digitizer/Singles/insert blurring<br></div><div>#/gate/digitizer/Singles/blurring/setResolution 0.01</div><div>
#/gate/digitizer/Singles/blurring/setEnergyOfReference 511. keV</div><div>/gate/digitizer/Singles/insert thresholder<br></div><div>/gate/digitizer/Singles/thresholder/setThreshold 350. keV</div><div>/gate/digitizer/Singles/insert upholder</div>
<div>/gate/digitizer/Singles/upholder/setUphold 650. keV</div><div>/gate/digitizer/Coincidences/setWindow 10. ns<br></div><div>/gate/digitizer/Coincidences/minSectorDifference 1</div><div>/gate/digitizer/name delay<br></div>
<div>/gate/digitizer/insert coincidenceSorter</div><div>/gate/digitizer/delay/setWindow 10. ns</div><div>/gate/digitizer/delay/setOffset 500. ns</div><div><br></div><div style>source is backtoback gamma 0.511 MeV at 0 0 0 </div>
<div style><br></div><div style>If I look at the coincidence events (globalPosX1, etc.), I was hoping the LORs would all go through the origin.</div><div style>the mean of all the intercepts is close to 0, but there is a residual 0.2 mm PSF.</div>
<div style>Can anyone tell me where this is coming from? something in the digitizer?</div><div style><br></div><div style>mlewis</div><div style><br></div></div></div>