scispace - formally typeset
Search or ask a question

Showing papers by "Angelo Schiavi published in 2022"


Journal ArticleDOI
TL;DR: In this article , the authors presented a data-driven model for carbon ion fragmentation and compared it with the FLUKA MC software, and a good agreement was found with respect to the results.
Abstract: The advent of Graphics Processing Units (GPU) has prompted the development of Monte Carlo (MC) algorithms that can significantly reduce the simulation time with respect to standard MC algorithms based on Central Processing Unit (CPU) hardware. The possibility to evaluate a complete treatment plan within minutes, instead of hours, paves the way for many clinical applications where the time-factor is important. FRED (Fast paRticle thErapy Dose evaluator) is a software that exploits the GPU power to recalculate and optimise ion beam treatment plans. The main goal when developing the FRED physics model was to balance accuracy, calculation time and GPU execution guidelines. Nowadays, FRED is already used as a quality assurance tool in Maastricht and Krakow proton clinical centers and as a research tool in several clinical and research centers across Europe. Lately the core software has been updated including a model of carbon ions interactions with matter. The implementation is phenomenological and based on carbon fragmentation data currently available. The model has been tested against the MC FLUKA software, commonly used in particle therapy, and a good agreement was found. In this paper, the new FRED data-driven model for carbon ion fragmentation will be presented together with the validation tests against the FLUKA MC software. The results will be discussed in the context of FRED clinical applications to 12C ions treatment planning.

3 citations


Journal ArticleDOI
TL;DR: In this article , the authors implemented a continuous scoring approach for the production of positron emitting isotopes within FRED version 5.59 and GATE version 9.0 to reduce computational resources and time required for simulation of patient activation during proton therapy using the GPU accelerated Monte Carlo code FRED.
Abstract: Objective. Verification of delivered proton therapy treatments is essential for reaping the many benefits of the modality, with the most widely proposed in vivo verification technique being the imaging of positron emitting isotopes generated in the patient during treatment using positron emission tomography (PET). The purpose of this work is to reduce the computational resources and time required for simulation of patient activation during proton therapy using the GPU accelerated Monte Carlo code FRED, and to validate the predicted activity against the widely used Monte Carlo code GATE. Approach. We implement a continuous scoring approach for the production of positron emitting isotopes within FRED version 5.59.9. We simulate treatment plans delivered to 95 head and neck patients at Centrum Cyklotronowe Bronowice using this GPU implementation, and verify the accuracy using the Monte Carlo toolkit GATE version 9.0. Main results. We report an average reduction in computational time by a factor of 50 when using a local system with 2 GPUs as opposed to a large compute cluster utilising between 200 to 700 CPU threads, enabling simulation of patient activity within an average of 2.9 min as opposed to 146 min. All simulated plans are in good agreement across the two Monte Carlo codes. The two codes agree within a maximum of 0.95σ on a voxel-by-voxel basis for the prediction of 7 different isotopes across 472 simulated fields delivered to 95 patients, with the average deviation over all fields being 6.4 × 10−3 σ. Significance. The implementation of activation calculations in the GPU accelerated Monte Carlo code FRED provides fast and reliable simulation of patient activation following proton therapy, allowing for research and development of clinical applications of range verification for this treatment modality using PET to proceed at a rapid pace.

3 citations



Journal ArticleDOI
TL;DR: In this paper , the authors presented a method to quantify the biological range extension of the proton beam, which results from the application of variable relative biological effectiveness (vRBE) approach in RBE-weighted dose calculation.
Abstract: Variable relative biological effectiveness (vRBE) in proton therapy might significantly modify the prediction of RBE-weighted dose delivered to a patient during proton therapy. In this study we will present a method to quantify the biological range extension of the proton beam, which results from the application of vRBE approach in RBE-weighted dose calculation.The treatment plans of 95 patients (brain and skull base patients) were used for RBE-weighted dose calculation with constant and the McNamara RBE model. For this purpose the Monte Carlo tool FRED was used. The RBE-weighted dose distributions were analysed using indices from dose-volume histograms. We used the volumes receiving at least 95% of the prescribed dose (V95) to estimate the biological range extension resulting from vRBE approach.The vRBE model shows higher median value of relative deposited dose and D95 in the planning target volume by around 1% for brain patients and 4% for skull base patients. The maximum doses in organs at risk calculated with vRBE was up to 14 Gy above dose limit. The mean biological range extension was greater than 0.4 cm.Our method of estimation of biological range extension is insensitive for dose inhomogeneities and can be easily used for different proton plans with intensity-modulated proton therapy (IMPT) optimization. Using volumes instead of dose profiles, which is the common method, is more universal. However it was tested only for IMPT plans on fields arranged around the tumor area.Adopting a vRBE model results in an increase in dose and an extension of the beam range, which is especially disadvantageous in cancers close to organs at risk. Our results support the need to re-optimization of proton treatment plans when considering vRBE.

1 citations


Journal ArticleDOI
TL;DR: In this article , the authors presented a method to quantify the biological range extension of the proton beam, which results from the application of variable relative biological effectiveness (vRBE) approach in RBE-weighted dose calculation.
Abstract: Variable relative biological effectiveness (vRBE) in proton therapy might significantly modify the prediction of RBE-weighted dose delivered to a patient during proton therapy. In this study we will present a method to quantify the biological range extension of the proton beam, which results from the application of vRBE approach in RBE-weighted dose calculation.The treatment plans of 95 patients (brain and skull base patients) were used for RBE-weighted dose calculation with constant and the McNamara RBE model. For this purpose the Monte Carlo tool FRED was used. The RBE-weighted dose distributions were analysed using indices from dose-volume histograms. We used the volumes receiving at least 95% of the prescribed dose (V95) to estimate the biological range extension resulting from vRBE approach.The vRBE model shows higher median value of relative deposited dose and D95 in the planning target volume by around 1% for brain patients and 4% for skull base patients. The maximum doses in organs at risk calculated with vRBE was up to 14 Gy above dose limit. The mean biological range extension was greater than 0.4 cm.Our method of estimation of biological range extension is insensitive for dose inhomogeneities and can be easily used for different proton plans with intensity-modulated proton therapy (IMPT) optimization. Using volumes instead of dose profiles, which is the common method, is more universal. However it was tested only for IMPT plans on fields arranged around the tumor area.Adopting a vRBE model results in an increase in dose and an extension of the beam range, which is especially disadvantageous in cancers close to organs at risk. Our results support the need to re-optimization of proton treatment plans when considering vRBE.

1 citations


Journal ArticleDOI
Marco Toppi, Alessio Sarti, A. Alexandrov, Behcet Alpat, G. Ambrosi, Stefano Argiro, Raul Arteche Diaz, Mattia Barbanera, Nazar Bartosik, G. Battistoni, Nicola Belcari, Silvia Biondi, M. G. Bisogni, M. Bon, Graziano Bruni, P. Carra, F. Cavanna, P. Cerello, Esther Ciarrocchi, Alberto Clozza, S. Colombi, G. De Lellis, A. De Gregorio, A. Del Guerra, M. De Simoni, A. Di Crescenzo, B. Di Ruzza, M. Donetti, Y. Dong, Maria Durante, Veronica Ferrero, E. Fiandrini, Ch. Finck, Elisa Fiorina, M. Fischetti, Matteo Francesconi, Matteo Franchini, Gaia Franciosini, G. Galati, Luca Galli, Giacomo Giraudo, R. Hetzel, E. Iarocci, Maria Ionica, A. Iuliano, Keida Kanxheri, A. C. Kraan, C. La Tessa, Martina Laurenza, Adele Lauria, E. López Torres, Michela Marafini, Maurizio Massa, Cristian Massimi, Ilaria Mattei, A Meneghetti, Alberto Mengarelli, R. Mirabelli, A. Moggi, Maria Cristina Montesi, M. C. Morone, Matteo Morrocchi, Silvia Muraro, F. Murtas, A. Muscato, A R Pastore, Nadia Pastrone, Vincenzo Patera, Francesco Pennazio, Francesco Peverini, Pisana Placidi, M. Pullia, Luciano Ramello, C.-A. Reidel, R. Ridolfi, Valeria Rosso, Claudio Sanelli, G. Sartorelli, Osamu Sato, Simone Savazzi, L. Scavarda, Angelo Schiavi, Christoph Schuy, Emanuele Scifoni, A. Sciubba, A. Secher, M. Selvi, L. Servoli, G. Silvestre, Mario Sitta, Roberto Spighi, E. Spiriti, Giancarlo Sportelli, Andre Stahl, Sandro Tomassini, Francesco Tommasino, V. Tioukov, Giacomo Traini, Antonio Trigilio, S. M. Valle, M. Vanstalle, Ulrich Weber, R. Zarrella, Antonio Zoccoli, Maria Pia Villa 
TL;DR: In this paper , the analysis of the data collected in the interactions of an oxygen ion beam of 400 MeV/u with a graphite target using a partial FOOT setup, at the GSI Helmholtz Center for Heavy Ion Research facility in Darmstadt, is presented.
Abstract: The study of nuclear fragmentation plays a central role in many important applications: from the study of Particle Therapy (PT) up to radiation protection for space (RPS) missions and the design of shielding for nuclear reactors. The FragmentatiOn Of Target (FOOT) collaboration aims to study the nuclear reactions that describe the interactions with matter of different light ions (like H 1 , H e 4 , C 12 , O 16 ) of interest for such applications, performing double differential fragmentation cross section measurements in the energy range of interest for PT and RPS. In this manuscript, we present the analysis of the data collected in the interactions of an oxygen ion beam of 400 MeV/u with a graphite target using a partial FOOT setup, at the GSI Helmholtz Center for Heavy Ion Research facility in Darmstadt. During the data taking the magnets, the silicon trackers and the calorimeter foreseen in the final FOOT setup were not yet available, and hence precise measurements of the fragments kinetic energy, momentum and mass were not possible. However, using the FOOT scintillator detectors for the time of flight (TOF) and energy loss (ΔE) measurements together with a drift chamber, used as beam monitor, it was possible to measure the elemental fragmentation cross sections. The reduced detector set-up and the limited available statistics allowed anyway to obtain relevant results, providing statistically significant measurements of cross sections eagerly needed for PT and RPS applications. Whenever possible the obtained results have been compared with existing measurements helping in discriminating between conflicting results in the literature and demonstrating at the same time the proper functioning of the FOOT ΔE-TOF system. Finally, the obtained fragmentation cross sections are compared to the Monte Carlo predictions obtained with the FLUKA software.

1 citations



Journal ArticleDOI
TL;DR: In this paper , a GPU-accelerated Monte Carlo software for dose valuation in electron and photon based radiotherapy, developed as an update of the FRED (Fast paRticle thErapy Dose evaluator) software, is presented.
Abstract: Objective. The Monte Carlo simulation software is a valuable tool in radiation therapy, in particular to achieve the needed accuracy in the dose evaluation for the treatment plans optimisation. The current challenge in this field is the time reduction to open the way to many clinical applications for which the computational time is an issue. In this manuscript we present an innovative GPU-accelerated Monte Carlo software for dose valuation in electron and photon based radiotherapy, developed as an update of the FRED (Fast paRticle thErapy Dose evaluator) software. Approach. The code transports particles through a 3D voxel grid, while scoring their energy deposition along their trajectory. The models of electromagnetic interactions in the energy region between 1 MeV—1 GeV available in literature have been implemented to efficiently run on GPUs, allowing to combine a fast tracking while keeping high accuracy in dose assessment. The FRED software has been bench-marked against state-of-art full MC (FLUKA, GEANT4) in the realm of two different radiotherapy applications: Intra-Operative Radio Therapy and Very High Electron Energy radiotherapy applications. Results. The single pencil beam dose-depth profiles in water as well as the dose map computed on non-homogeneous phantom agree with full-MCs at 2% level, observing a gain in processing time from 200 to 5000. Significance. Such performance allows for computing a plan with electron beams in few minutes with an accuracy of ∼%, demonstrating the FRED potential to be adopted for fast plan re-calculation in photon or electron radiotherapy applications.



DOI
G. Silvestre, Francesca Peverini, L. Servoli, Keida Kanxheri, Mattia Barbanera, Pisana Placidi, Maria Ionica, M. Caprai, G. Ambrosi, Lucia Salvi, E. Fiandrini, Behcet Alpat, A. Alexandrov, Stefano Argiro, Diaz Raul Arteche, Nazar Bartosik, G. Battistoni, Nicola Belcari, E.V. Bellinzona, Silvia Biondi, Maria Giuseppina Bisogni, Graziano Bruni, P. Carra, F. Cavanna, P. Cerello, Esther Ciarrocchi, Alberto Clozza, S. Colombi, A. De Gregorio, G. De Lellis, A. Del Guerra, M. De Simoni, A. Di Crescenzo, B. Di Ruzza, M. Donetti, Yu Dong, Marco Durante, R. Faccini, Veronica Ferrero, Ch. Finck, Elisa Fiorina, M. Fischetti, Marco Francesconi, Matteo Franchini, Gaia Franciosini, Giuliana Galati, Luca Galli, V. Gentile, Giuseppe Giraudo, R. Hetzel, E. Iarocci, A. C. Kraan, C. La Tessa, Martina Laurenza, Adele Lauria, E. López Torres, Alice Manna, Michela Marafini, Maurizio Massa, Cristian Massimi, Ilaria Mattei, Alessio Mereghetti, Alberto Mengarelli, A. Moggi, Maria Cristina Montesi, M. C. Morone, Matteo Morrocchi, Silvia Muraro, Alessandro Pastore, Nadia Pastrone, Vincenzo Patera, Francesco Pennazio, M. Pullia, Luciano Ramello, C.-A. Reidel, R. Ridolfi, Valeria Rosso, Claudio Sanelli, Alessio Sarti, G. Sartorelli, Osamu Sato, Simone Savazzi, L. Scavarda, Angelo Schiavi, Christoph Schuy, Emanuele Scifoni, A. Sciubba, A. Secher, M. Selvi, Mario Sitta, Roberto Spighi, E. Spiriti, Giancarlo Sportelli, Achim Stahl, Sandro Tomassini, Francesco Tommasino, Marco Toppi, Giacomo Traini, Antonio Trigilio, Giacomo Ubaldi, V. Tioukov, A. Valetti, S. M. Valle, M. Vanstalle, Mauro Villa, Ulrich Weber, R. Zarrella, Antonio Zoccoli 
TL;DR: In this article , the MSD (Microstrip Silicon Detector) was used for the FOOT (FragmentatiOn Of Target) experiment, which measured the proton double differential fragmentation cross-section on H, C, O targets at beam energies of interest for hadrontherapy (50-250 MeV for protons and 50-400 MeV/u for carbon ions), and also at higher energy, up to 1 GeV/U for radioprotection in space.
Abstract: The goals of the FOOT (FragmentatiOn Of Target) experiment are to measure the proton double differential fragmentation cross-section on H, C, O targets at beam energies of interest for hadrontherapy (50–250 MeV for protons and 50–400 MeV/u for carbon ions), and also at higher energy, up to 1 GeV/u for radioprotection in space. Given the short range of the fragments, an inverse kinematic approach has been chosen, requiring precise tracking capabilities for charged particles. One of the subsystems designed for the experiment will be the MSD (Microstrip Silicon Detector), consisting of three x-y measurement planes, each one made by two single sided silicon microstrip sensors. In this document, we will present a detailed description of the first MSD prototype assembly, developed by INFN Perugia group and the subsequent characterization of the detector performance. The prototype is a wide area (∼ 100 cm2) single sensor, 150 μm thick to reduce material budget and fragmentation probability along the beam path, with 50 μm strip pitch and 2 floating strip readout approach. The pitch adapter to connect strips with the readout channels of the ASIC has been implemented directly on the silicon surface. Beside the interest for the FOOT experiment, the results in terms of cluster signal, signal-to-noise ratio, dynamic range of the readout chips, as well as long-term stability studies in terms of noise, are relevant also for other experiments where the use of thin sensors is crucial.