The FLUKA Monte Carlo Radiation Transport and Interaction code package is widely used to simulate the interaction of particles with matter in a variety of fields, including high energy physics, space radiation, medical applications, radiation protection and shielding assessments, accelerator studies, astrophysical studies and well logging. This paper gives a brief overview of the FLUKA program and describes recent developments, in particular, improvements in the modelling of particle interactions and transport are described in detail. In addition, an overview of selected applications is given.

During particle therapy irradiation, positron emitters with half-lives ranging from 2 to 20 min are generated from nuclear processes. The half-lives are such that it is possible either to detect the positron signal in the treatment room using an in-beam positron emission tomography (PET) system, right after the irradiation, or to quickly transfer the patient to a close PET/CT scanner. Since the activity distribution is spatially correlated with the dose, it is possible to use PET imaging as an indirect method to assure the quality of the dose delivery. In this work, we present a new dedicated PET system able to operate in-beam. The PET apparatus consists in two 10 cm × 10 cm detector heads. Each detector is composed of four scintillating matrices of 23 × 23 LYSO crystals. The crystal size is 1.9 mm × 1.9 mm × 16 mm. Each scintillation matrix is read out independently with a modularized acquisition system. The distance between the two opposing detector heads was set to 20 cm. The system has very low dead time per detector area and a 3 ns coincidence window, which is capable to sustain high single count rates and to keep the random counts relatively low. This allows a new full-beam monitoring modality that includes data acquisition also while the beam is on. The PET system was tested during the irradiation at the CATANA (INFN, Catania, Italy) cyclotron-based proton therapy facility. Four acquisitions with different doses and dose rates were analysed. In all cases the random to total coincidences ratio was equal or less than 25%. For each measurement we estimated the accuracy and precision of the activity range on a set of voxel lines within an irradiated PMMA phantom. Results show that the inclusion of data acquired during the irradiation, referred to as beam-on data, improves both the precision and accuracy of the range measurement with respect to data acquired only after irradiation. Beam-on data alone are enough to give precisions better than 1 mm when at least 5 Gy are delivered.

Objective . Monolithic scintillator crystals coupled to silicon photomultiplier (SiPM) arrays are promising detectors for PET applications, offering spatial resolution around 1 mm and depth-of-interaction information. However, their timing resolution has always been inferior to that of pixellated crystals, while the best results on spatial resolution have been obtained with algorithms that cannot operate in real-time in a PET detector. In this study, we explore the capabilities of monolithic crystals with respect to spatial and timing resolution, presenting new algorithms that overcome the mentioned problems. Approach. Our algorithms were tested first using a simulation framework, then on experimentally acquired data. We tested an event timestamping algorithm based on neural networks which was then integrated into a second neural network for simultaneous estimation of the event position and timestamp. Both algorithms are implemented in a low-cost field-programmable gate array that can be integrated in the detector and can process more than 1 million events per second in real-time. Results. Testing the neural network for the simultaneous estimation of the event position and timestamp on experimental data we obtain 0.78 2D FWHM on the ( x , y ) plane, 1.2 depth-of-interaction FWHM and 156 coincidence time resolution on a 25 mm × 25 mm × 8 mm × LYSO monolith read-out by 64 3 mm × 3 mm Hamamatsu SiPMs. Significance. Our results show that monolithic crystals combined with artificial intelligence can rival pixellated crystals performance for time-of-flight PET applications, while having better spatial resolution and DOI resolution. Thanks to the use of very light neural networks, event characterization can be done on-line directly in the detector, solving the issues of scalability and computational complexity that up to now were preventing the use of monolithic crystals in clinical PET scanners.

Next generation PET scanners should fulfill very high requirements in terms of spatial, energy and timing resolution. Modern scanner performances are inherently limited by the use of standard photomultiplier tubes. The use of Silicon Photomultipliers (SiPMs) is proposed for the construction of a 4D-PET module of 4.8×4.8 cm² aimed to replace the standard PMT based PET block detector. The module will be based on a LYSO continuous crystal read on two faces by Silicon Photomultipliers. A high granularity detection surface made by SiPM matrices of 1.5 mm pitch will be used for the x–y photon hit position determination with submillimetric accuracy, while a low granularity surface constituted by 16 mm² SiPM pixels will provide the fast timing information (t) that will be used to implement the Time of Flight technique (TOF). The spatial information collected by the two detector layers will be combined in order to measure the Depth of Interaction (DOI) of each event (z). The use of large area multi-pixel Silicon Photomultiplier (SiPM) detectors requires the development of a multichannel Data Acquisition system (DAQ) as well as of a dedicated front-end in order not to degrade the intrinsic detector capabilities and to manage many channels. The paper describes the progress made on the development of the proof of principle module under construction at the University of Pisa.