Monte Carlo simulation is an essential tool in emission tomography that can assist in the design of new medical imaging devices, the optimization of acquisition protocols and the development or assessment of image reconstruction algorithms and correction techniques. GATE, the Geant4 Application for Tomographic Emission, encapsulates the Geant4 libraries to achieve a modular, versatile, scripted simulation toolkit adapted to the field of nuclear medicine. In particular, GATE allows the description of time-dependent phenomena such as source or detector movement, and source decay kinetics. This feature makes it possible to simulate time curves under realistic acquisition conditions and to test dynamic reconstruction algorithms. This paper gives a detailed description of the design and development of GATE by the OpenGATE collaboration, whose continuing objective is to improve, document and validate GATE by simulating commercially available imaging systems for PET and SPECT. Large effort is also invested in the ability and the flexibility to model novel detection systems or systems still under design. A public release of GATE licensed under the GNU Lesser General Public License can be downloaded at http:/www-lphe.epfl.ch/GATE/. Two benchmarks developed for PET and SPECT to test the installation of GATE and to serve as a tutorial for the users are presented. Extensive validation of the GATE simulation platform has been started, comparing simulations and measurements on commercially available acquisition systems. References to those results are listed. The future prospects towards the gridification of GATE and its extension to other domains such as dosimetry are also discussed.

GATE (Geant4 Application for Emission Tomography) is a Monte Carlo simulation platform developed by the OpenGATE collaboration since 2001 and first publicly released in 2004. Dedicated to the modelling of planar scintigraphy, single photon emission computed tomography (SPECT) and positron emission tomography (PET) acquisitions, this platform is widely used to assist PET and SPECT research. A recent extension of this platform, released by the OpenGATE collaboration as GATE V6, now also enables modelling of x-ray computed tomography and radiation therapy experiments. This paper presents an overview of the main additions and improvements implemented in GATE since the publication of the initial GATE paper (Jan et al 2004 Phys. Med. Biol. 49 4543-61). This includes new models available in GATE to simulate optical and hadronic processes, novelties in modelling tracer, organ or detector motion, new options for speeding up GATE simulations, examples illustrating the use of GATE V6 in radiotherapy applications and CT simulations, and preliminary results regarding the validation of GATE V6 for radiation therapy applications. Upon completion of extensive validation studies, GATE is expected to become a valuable tool for simulations involving both radiotherapy and imaging.

In this paper, the authors' review the applicability of the open-source GATE Monte Carlo simulation platform based on the GEANT4 toolkit for radiation therapy and dosimetry applications. The many applications of GATE for state-of-the-art radiotherapy simulations are described including external beam radiotherapy, brachytherapy, intraoperative radiotherapy, hadrontherapy, molecular radiotherapy, and in vivo dose monitoring. Investigations that have been performed using GEANT4 only are also mentioned to illustrate the potential of GATE. The very practical feature of GATE making it easy to model both a treatment and an imaging acquisition within the same framework is emphasized. The computational times associated with several applications are provided to illustrate the practical feasibility of the simulations using current computing facilities.

Silicon photomultipliers (SiPMs) are of great interest to positron emission tomography (PET), as they enable new detector geometries, for e.g., depthof-interaction (DOI) determination, are MR compatible, and offer faster response and higher gain than other solid-state photosensors such as avalanche photodiodes. Here we present a novel detector design with DOI correction, in which a position-sensitive SiPM array is used to read out a monolithic scintillator. Initial characterization of a prototype detector consisting of a 4 × 4 SiPM array coupled to either the front or back surface of a 13.2 mm × 13.2 mm × 10 mm LYSO:Ce 3+ crystal shows that front-side readout results in significantly better performance than conventional back-side readout. Spatial resolutions <1.6 mm full-width-at-half-maximum (FWHM) were measured at the detector centre in response to an ∼0.54 mm FWHM diameter test beam. Hardly any resolution losses were observed at angles of incidence up to 45• , demonstrating excellent DOI correction. About 14% FWHM energy resolution was obtained. The timing resolution, measured in coincidence with a BaF 2 detector, equals 960 ps FWHM.

The timing performance of scintillation detectors is ultimately limited by photon counting statistics. In fact, photon counting statistics form a dominant contribution to the overall timing resolution of many state-of-the-art detectors. A common approach to investigate this contribution is to calculate the variance in the registration times of individual scintillation photons within the photosensor. However, in general the single-photon variance is not equal to the intrinsic limit on the timing resolution, since in principle one can make use of the timing information carried by all photons detected. In this work, the Cramér-Rao lower bound on the timing resolution of a scintillation detector, based on the information contained in the full set of registered photons, is calculated. The results appear to be in good agreement with trends observed in the literature. Furthermore, it is shown that the timestamp obtained from any single scintillation photon never yields the optimum timing resolution for realistic scintillation detectors. Yet, it appears that the intrinsic timing resolution limit can be approached closely by making use of the timestamps from a relatively small number of photons emitted during the initial part of the scintillation pulse.