Abstract: Purpose: In the present work, the authors compare geometrical and Monte Carlo projectors in detail. The geometrical projectors considered were the conventional geometrical Siddon ray-tracer (S-RT) and the orthogonal distance-based ray-tracer (OD-RT), based on computing the orthogonal distance from the center of image voxel to the line-of-response. A comparison of these geometrical projectors was performed using different point spread function (PSF) models. The Monte Carlo-based method under consideration involves an extensive model of the system response matrix based on Monte Carlo simulations and is computed off-line and stored on disk. Methods: Comparisons were performed using simulated and experimental data of the commercial small animal PET scanner rPET. Results: The results demonstrate that the orthogonal distance-based ray-tracer and Siddon ray-tracer using PSF image-space convolutions yield better images in terms of contrast and spatial resolution than those obtained after using the conventional method and the multiray-based S-RT. Furthermore, the Monte Carlo-based method yields slight improvements in terms of contrast and spatial resolution with respect to these geometrical projectors. Conclusions: The orthogonal distance-based ray-tracer and Siddon ray-tracer using PSF image-space convolutions represent satisfactory alternatives to factorizing the system matrix or to the conventional on-the-fly ray-tracing methods for list-mode reconstruction, where an extensive modeling based on Monte Carlo simulations is unfeasible.

Abstract: Random coincidences ('randoms') are one of the main sources of image degradation in PET imaging. In order to correct for this effect, an accurate method to estimate the contribution of random events is necessary. This aspect becomes especially relevant for high-resolution PET scanners where the highest image quality is sought and accurate quantitative analysis is undertaken. One common approach to estimate randoms is the so-called singles rate method (SR) widely used because of its good statistical properties. SR is based on the measurement of the singles rate in each detector element. However, recent studies suggest that SR systematically overestimates the correct random rate. This overestimation can be particularly marked for low energy thresholds, below 250 keV used in some applications and could entail a significant image degradation. In this work, we investigate the performance of SR as a function of the activity, geometry of the source and energy acceptance window used. We also investigate the performance of an alternative method, which we call 'singles trues' (ST) that improves SR by properly modeling the presence of true coincidences in the sample. Nevertheless, in any real data acquisition the knowledge of which singles are members of a true coincidence is lost. Therefore, we propose an iterative method, STi, that provides an estimation based on ST but which only requires the knowledge of measurable quantities: prompts and singles. Due to inter-crystal scatter, for wide energy windows ST only partially corrects SR overestimations. While SR deviations are in the range 86-300% (depending on the source geometry), the ST deviations are systematically smaller and contained in the range 4-60%. STi fails to reproduce the ST results, although for not too high activities the deviation with respect to ST is only a few percent. For conventional energy windows, i.e. those without inter-crystal scatter, the ST method corrects the SR overestimations, and deviations from the true random rate are of the order of 1% or less. In addition, in the case of conventional energy window STi results reproduce ST results and therefore the former can be used to obtain the true random rate.

Abstract: Charged particle therapy, or so-called hadrontherapy, is developing very rapidly. There is large pressure on the scientific community to deliver dedicated accelerators, providing the best possible treatment modalities at the lowest cost. In this context, the Italian research Foundation TERA is developing fast-cycling accelerators, dubbed 'cyclinacs'. These are a combination of a cyclotron (accelerating ions to a fixed initial energy) followed by a high gradient linac boosting the ions energy up to the maximum needed for medical therapy. The linac is powered by many independently controlled klystrons to vary the beam energy from one pulse to the next. This accelerator is best suited to treat moving organs with a 4D multipainting spot scanning technique. A dual proton/carbon ion cyclinac is here presented. It consists of an Electron Beam Ion Source, a superconducting isochronous cyclotron and a high-gradient linac. All these machines are pulsed at high repetition rate (100-400 Hz). The source should deliver both C6+ and H-2(+) ions in short pulses (1.5 μs flat-top) and with sufficient intensity (at least 10(8) fully stripped carbon ions per pulse at 300 Hz). The cyclotron accelerates the ions to 120 MeV/u. It features a compact design (with superconducting coils) and a low power consumption. The linac has a novel C-band high-gradient structure and accelerates the ions to variable energies up to 400 MeV/u. High RF frequencies lead to power consumptions which are much lower than the ones of synchrotrons for the same ion extraction energy. This work is part of a collaboration with the CLIC group, which is working at CERN on high-gradient electron-positron colliders.