PET digitization chain for Monte Carlo simulation in GATE

Objective. We introduce a versatile methodology for the accurate modelling of PET imaging systems via Monte Carlo simulations, using the Geant4 application for tomographic emission (GATE) platform. Accurate Monte Carlo modelling involves the incorporation of a complete analytical signal processing chain, called the digitizer in GATE, to emulate the different count rates encountered in actual positron emission tomography (PET) systems. Approach. The proposed approach consists of two steps: (1) modelling the digitizer to replicate the detection chain of real systems, covering all available parameters, whether publicly accessible or supplied by manufacturers; (2) estimating the remaining parameters, i.e. background noise level, detection efficiency, and pile-up, using optimisation techniques based on experimental single and prompt event rates. We show that this two-step optimisation reproduces the other experimental count rates (true, scatter, and random), without the need for additional adjustments. This method has been applied and validated with experimental data derived from the NEMA count losses test for three state-of-the-art SiPM-based time-of-flight (TOF)-PET systems: Philips Vereos, Siemens Biograph Vision 600 and GE Discovery MI 4-ring. Main results. The results show good agreement between experiments and simulations for the three PET systems, with absolute relative discrepancies below 3%, 6%, 6%, 7% and 12% for prompt, random, true, scatter and noise equivalent count rates, respectively, within the 0-10 kBq<middle dot>ml-1 activity concentration range typically observed in whole-body 18F scans. Significance. Overall, the proposed digitizer optimisation method was shown to be effective in reproducing count rates and NECR for three of the latest generation SiPM-based TOF-PET imaging systems. The proposed methodology could be applied to other PET scanners.