Abstract: Compton Cameras are electronically collimated photon imagers suitable for sub-MeV to few MeV gamma-ray detection. Such features are desirable to enablein vivorange verification in hadron therapy, through the detection of secondary Prompt Gammas. A major concern with this technique is the poor image quality obtained when the incoming gamma-ray energy is unknown. Compton Cameras with more than two detector planes (multi-layer Compton Cameras) have been proposed as a solution, given that these devices incorporate more signal sequences of interactions to the conventional two interaction events. In particular, three interaction events convey more spectral information as they allow inferring directly the incident gamma-ray energy. A three-layer Compton Telescope based on continuous Lanthanum (III) Bromide crystals coupled to Silicon Photomultipliers is being developed at the IRIS group of IFIC-Valencia. In a previous work we proposed a spectral reconstruction algorithm for two interaction events based on an analytical model for the formation of the signal. To fully exploit the capabilities of our prototype, we present here an extension of the model for three interaction events. Analytical expressions of the sensitivity and the System Matrix are derived and validated against Monte Carlo simulations. Implemented in a List Mode Maximum Likelihood Expectation Maximization algorithm, the proposed model allows us to obtain four-dimensional (energy and position) images by using exclusively three interaction events. We are able to recover the correct spectrum and spatial distribution of gamma-ray sources when ideal data are employed. However, the uncertainties associated to experimental measurements result in a degradation when real data from complex structures are employed. Incorrect estimation of the incident gamma-ray interaction positions, and missing deposited energy associated with escaping secondaries, have been identified as the causes of such degradation by means of a detailed Monte Carlo study. As expected, our current experimental resolution and efficiency to three interaction events prevents us from correctly recovering complex structures of radioactive sources. However, given the better spectral information conveyed by three interaction events, we expect an improvement of the image quality of conventional Compton imaging when including such events. In this regard, future development includes the incorporation of the model assessed in this work to the two interaction events model in order to allow using simultaneously two and three interaction events in the image reconstruction.

Abstract: Sensitivity and spatial resolution are the main parameters to maximize in the performance of a PET scanner. For this purpose, detectors consisting of a combination of continuous crystals optically coupled to segmented photodetectors have been employed. With the use of continuous crystals the sensitivity is increased with respect to the pixelated crystals. In addition, spatial resolution is no longer limited to the crystal size. The main drawback is the difficulty in determining the interaction position. In this work, we present the characterization of the performance of a full ring based on cuboid continuous crystals coupled to SiPMs. To this end, we have employed the simulations developed in a previous work for our experimental detector head. Sensitivity could be further enhanced by using tapered crystals. This enhancement is obtained by increasing the solid angle coverage, reducing the wedge-shaped gaps between contiguous detectors. The performance of the scanners based on both crystal geometries was characterized following NEMA NU 4-2008 standardized protocol in order to compare them. An average sensitivity gain over the entire axial field of view of 13.63% has been obtained with tapered geometry while similar performance of the spatial resolution has been proven with both scanners. The activity at which NECR and true peak occur is smaller and the peak value is greater for tapered crystals than for cuboid crystals. Moreover, a higher degree of homogeneity was obtained in the sensitivity map due to the tighter packing of the crystals, which reduces the gaps and results in a better recovery of homogeneous regions than for the cuboid configuration. Some of the results obtained, such as spatial resolution, depend on the interaction position estimation and may vary if other method is employed.

Abstract: The applicability extent of hadron therapy for tumor treatment is currently limited by the lack of reliable online monitoring techniques. An active topic of investigation is the research of monitoring systems based on the detection of secondary radiation produced during treatment. MACACO, a multi-layer Compton camera based on LaBr3 scintillator crystals and SiPMs, is being developed at IFIC-Valencia for this purpose. This work reports the results obtained from measurements of a 150 MeV proton beam impinging on a PMMA target. A neural network trained on Monte Carlo simulations is used for event selection, increasing the signal to background ratio before image reconstruction. Images of the measured prompt gamma distributions are reconstructed by means of a spectral reconstruction code, through which the 4.439 MeV spectral line is resolved. Images of the emission distribution at this energy are reconstructed, allowing calculation of the distal fall-off and identification of target displacements of 3 mm.

Abstract: The characterization of a PET detector head based on continuous LYSO crystals and silicon photomultiplier (SiPM) arrays as photodetectors has been carried out for its use in the development of a small animal PET prototype. The detector heads are composed of a continuous crystal and a SiPM matrix with 64 pixels in a common substrate, fabricated specifically for this project. Three crystals of 12 mm x 12 mm x 5 mm size with different types of painting have been tested: white, black and black on the sides but white on the back of the crystal. The best energy resolution, obtained with the white crystal, is 16% FWHM. The detector response is linear up to 1275 keV. Tests with different position determination algorithms have been carried out with the three crystals. The spatial resolution obtained with the center of gravity algorithm is around 0.9 mm FWHM for the three crystals. As expected, the use of this algorithm results in the displacement of the reconstructed position toward the center of the crystal, more pronounced in the case of the white crystal. A maximum likelihood algorithm has been tested that can reconstruct correctly the interaction position of the photons also in the case of the white crystal.

Abstract: The interest in using continuous monolithic crystals in positron emission tomography (PET) has grown in the last years. Coupled to silicon photomultipliers (SiPMs), the detector can combine high sensitivity and high resolution, the two main factors to be maximized in a positron emission tomograph. In this work, the position determination capability of a detector comprised of a 12 x 12 x 10 mm(3) LYSO crystal coupled to an 8 x 8-pixel array of SiPMs is evaluated. The 3D interaction position of.-rays is estimated using an analytical model of the light distribution including reflections on the facets of the crystal. Monte Carlo simulations have been performed to evaluate different crystal reflectors and geometries. The method has been characterized and applied to different cases. Intrinsic resolution obtained with the position estimation method used in this work, applied to experimental data, achieves sub-millimetre resolution values. Average resolution over the detector surface for 5 mm thick crystal is similar to 0.9 mm FWHM and similar to 1.2 mm FWHM for 10 mm thick crystal. Depth of interaction resolution is close to 2 mm FWHM in both cases, while the FWTM is similar to 5.3 mm for 5 mm thick crystal and similar to 9.6 mm for 10 mm thick crystal.