A Monte Carlo model of noise components in 3D PET

In this work, a model to estimate noise components (scatter and random coincidences) in 3D PET was developed by Monte Carlo methods. The model allows the amount and the spatial distribution of true, scattered and random coincidences to be independently estimated for any radioactive source (both phantoms and real patients), properly accounting for system dead time. The model was applied to a 3D NaI(Tl) current generation PET scanner for which no methods to estimate scatter and random components in wholebody studies are currently available. The quantitative accuracy of the developed noise model was proved by comparison between simulated and measured PET scanner physical performance. Scatter and random components were assessed for brain and thorax studies giving evidence to the high scatter contribution and to the relatively low random component In whole-body 3D PET studies. The clinical response of the PET system, in terms of signal-to-noise ratio, was assessed and optimized, confirming the suitability of the default energy window, although suggesting a possible improvement by setting a lower energy threshold higher than the current default one. The proposed noise model applies to any current generation 3D PET scanner and was included in the Monte Carlo software package PET-EGS, devoted to 3D PET and freely available from authors