Realistic PET Monte Carlo Simulation With Pixelated Block Detectors, Light Sharing, Random Coincidences and Dead-Time Modeling

We have developed and validated a realistic simulation of random coincidences, pixelated block detectors, light sharing among crystal elements and dead-time in 2D and 3D positron emission tomography (PET) imaging based on the SimSET Monte Carlo simulation software. Our simulation was validated by comparison to a Monte Carlo transport code widely used for PET modeling, GATE, and to measurements made on a PET scanner. Methods: We have modified the SimSET software to allow independent tracking of single photons in the object and septa while taking advantage of existing voxel based attenuation and activity distributions and validated importance sampling techniques implemented in SimSET. For each single photon interacting in the detector, the energy-weighted average of interaction points was computed, a blurring model applied to account for light sharing and the associated crystal identified. Detector dead-time was modeled in every block as a function of the local single rate using a variance reduction technique. Electronic dead-time was modeled for the whole scanner as a function of the prompt coincidences rate. Energy spectra predicted by our simulation were compared to GATE. NEMA NU-2 2001 performance tests were simulated with the new simulation as well as with SimSET and compared to measurements made on a Discovery ST (DST) camera. Results: Errors in simulated spatial resolution (full width at half maximum, FWHM) were 5.5% (6.1%) in 2D (3D) with the new simulation, compared with 42.5% (38.2%) with SimSET. Simulated (measured) scatter fractions were 17.8% (21.3%) in 2D and 45.8% (45.2%) in 3D. Simulated and measured sensitivities agreed within 2.3% in 2D and 3D for all planes and simulated and acquired count rate curves (including NEC) were within 12.7% in 2D in the [0: 80 kBq/cc] range and in 3D in the [0: 35 kBq/cc] range. The new simulation yielded significantly more realistic singles´ and coincidences´ spectra, spatial resolution, global sensitivity and lesion c- - ontrasts than the SimSET software.