PETSIM: Monte Carlo simulation of all sensitivity and resolution parameters of cylindrical positron imaging systems

Abstract

Monte Carlo simulation techniques are applied to track the annihilation photons from positron decay, and store the photon histories. Reasonably realistic models of the isotope distribution in the brain and heart during typical PET studies, as well as the traditional phantoms used for measuring PET scanner performance can be built out of up to 10 hollow or solid cylinders. Separate programs model the source distribution and its attenuation characteristics, the collimators and the detectors. These modules are connected by compact gamma history files which are stored on disc or tape. Over 50 million gamma ray histories can be saved on a 1 Gbyte disc, representing the decay of several billion atoms. This allows for good precision even for single thin slices in scanners with wide axial acceptance. The simulation results include spectrum analysis, sensitivity to true coincident events, scattered coincident and single rays, and the effects on these parameters of detector dead time. The storage of intermediate results on tape reduces simulation time, since most common source geometries need be generated only once. The sensitivities in multi-slice systems are presented as matrices of coincident crystal planes.