Abstract: In order to exploit the advantages of ion-beam therapy in a clinical setting, delivery verification techniques are necessary to detect deviations from the planned treatment. Efforts are currently oriented towards the development of devices for real-time range monitoring. Among the different detector concepts proposed, Compton cameras are employed to detect prompt gammas and represent a valid candidate for real-time range verification. We present the first on-beam test of MACACO, a Compton telescope (multi-layer Compton camera) based on lanthanum bromide crystals and silicon photo-multipliers. The Compton telescope was first characterized through measurements and Monte Carlo simulations. The detector linearity was measured employing Na-22 and Am-Be sources, obtaining about 10% deviation from linearity at 3.44 MeV. A spectral image reconstruction algorithm was tested on synthetic data. Point-like sources emitting gamma rays with energy between 2 and 7 MeV were reconstructed with 3-5 mm resolution. The two-layer Compton telescope was employed to measure radiation emitted from a beam of 150 MeV protons impinging on a cylindrical PMMA target. Bragg-peak shifts were achieved via adjustment of the PMMA target location and the resulting measurements used during image reconstruction. Reconstructed Bragg peak profiles proved sufficient to observe peak-location differences within 10 mm demonstrating the potential of the MACACO Compton Telescope as a monitoring device for ion-beam therapy.

Abstract: Given the strong variations in the sensitivity of Compton cameras for the detection of events originating from different points in the field of view (FoV), sensitivity correction is often necessary in Compton image reconstruction. Several approaches for the calculation of the sensitivity matrix have been proposed in the literature. While most of these models are easily implemented and can be useful in many cases, they usually assume high angular coverage over the scattered photon, which is not the case for our prototype. In this work, we have derived an analytical model that allows us to calculate a detailed sensitivity matrix, which has been compared to other sensitivity models in the literature. Specifically, the proposed model describes the probability of measuring a useful event in a two-plane Compton camera, including the most relevant physical processes involved. The model has been used to obtain an expression for the system and sensitivity matrices for iterative image reconstruction. These matrices have been validated taking Monte Carlo simulations as a reference. In order to study the impact of the sensitivity, images reconstructed with our sensitivity model and with other models have been compared. Images have been reconstructed from several simulated sources, including point-like sources and extended distributions of activity, and also from experimental data measured with Na-22 sources. Results show that our sensitivity model is the best suited for our prototype. Although other models in the literature perform successfully in many scenarios, they are not applicable in all the geometrical configurations of interest for our system. In general, our model allows to effectively recover the intensity of point-like sources at different positions in the FoV and to reconstruct regions of homogeneous activity with minimal variance. Moreover, it can be employed for all Compton camera configurations, including those with low angular coverage over the scatterer.

Abstract: Compton cameras are gamma-ray imaging systems which have been proposed for a wide variety of applications such as medical imaging, nuclear decommissioning or homeland security. In the design and optimization of such a system Monte Carlo simulations play an essential role. In this work, we propose a generic module to perform Monte Carlo simulations and analyses of Compton Camera imaging which is included in the open-source GATE/Geant4 platform. Several digitization stages have been implemented within the module to mimic the performance of the most commonly employed detectors (e.g. monolithic blocks, pixelated scintillator crystals, strip detectors...). Time coincidence sorter and sequence coincidence reconstruction are also available in order to aim at providing modules to facilitate the comparison and reproduction of the data taken with different prototypes. All processing steps may be performed during the simulation (on-the-fly mode) or as a post-process of the output files (offline mode). The predictions of the module have been compared with experimental data in terms of energy spectra, angular resolution, efficiency and back-projection image reconstruction. Consistent results within a 3-sigma interval were obtained for the energy spectra except for low energies where small differences arise. The angular resolution measure for incident photons of 1275 keV was also in good agreement between both data sets with a value close to 13 degrees. Moreover, with the aim of demonstrating the versatility of such a tool the performance of two different Compton camera designs was evaluated and compared.