A Geant4 Extended Example
for the calculation of relevant radiobiological quantities
Giuliana Navarra - 05/11/2021
Overview
Proton and ion beams therapy
The first cancer treatment with a proton beam was carried out about 70 years ago.
During these decades the clinical use of proton and ion beams for cancer treatment has improved a lot as you can see from the spread of the treatment centers and from the big amount of scientific publication about this topic. ��Unlike electrons, protons allow a high dose deposition in a small deep-seated space; that causes a major precision on irradiating the target and less damage of surrounding and healthy tissues.
From the first time more than 160,000 people are been treated with protons worldwide.
The CATANA (Centro di AdroTerapia ed Applicazioni Nucleari Avanzate) facility, at LNS-INFN in Catania, started its activity on 2002 and more than 300 patients have been treated there until now.
references: [1], [2]
Radiobiological quantities
DOSE [Gy]
the energy deposited by ionizing radiation per unit mass.
RBE
The Relative Biological Effectiveness is the ratio between the dose of a standard radiation that produces a certain biological damage and the dose of the test radiation needed to produce the same damage
LET [keV/μm]
The Linear Energy Transfer is a measure of the energy transferred to the matter by an ionising particle.
references: [3], [4]
GEANT4
GEANT4 is a software toolkit for the simulation of the passage of particles through matter. It is written in C++ and exploits advanced software-engineering techniques and object-oriented technology.
It is used in different application fields:
references: [5]
All aspects of the simulation process have been included in the toolkit:
references: [5]
The advanced example: “Hadrontherapy”
references: [1], [4]
Handrontherapy is a free and open source Monte Carlo application. It uses the libraries of the Geant4 simulation toolkit and it is included in the advanced examples.�The geometry is divided in two main blocks: one consists in the geometry of a typical beam line for proton/ion therapy, another is the detection region.�All the transport elements are included.
CATANA
View of the CATANA beam line: 1. treatment chair for patient immobilization; 2. final collimator; 3. positioning laser; 4. light field simulator; 5. monitor chambers; 6. intermediate collimator; 7. box for the location of modulator wheel and range shifter.
The CATANA facility is in Catania at INFN - LNS. It is a passive proton beam line (shown in the figure).
The proton beam go through two scattering foils made of tantalum: the first one is 15 μm thick and it’s placed in vacuum before the beam exits in air; the second one is 25 μm thick and it has a central brass stopper with a diameter of 4 mm.
references: [6]
The scattering system is designed to obtain a homogenous lateral off-axis dose distribution and to minimize the energy loss.
The exit window is made of 50 μm kapton and it is about 3 meters far from isocenter.
Inside a box placed downstream of the second scattering foil, there are a range shifter and a range modulator.
The facility is provided with two transmission monitor ionization chamber to obtain an on-line control of the dose delivered. A light field allows to simulate the radiation field.
The isocenter identification and the patient centering are obtained by using two diode lasers: one is orthogonal to beam line and the other is coaxial to beam line. Before isocenter (83 mm upstream) there is a brass shaped patient collimator.
references: [6]
A new example: “The Extended Example”
The new application enable the user to simulate, in a simple way, the effects of radiation on a detector region (a parallelepiped).
The classes Dose, LET and RBE allow to obtain these quantities as outputs. An interesting characteristic of the example is the presence of the class Radiobiomanager that makes easy adding new classes in the code to implement the calculation of others radiobiological quantities.
The particle source can be customized by using the Primary Generator Action.
The utilization is straightforward for an inexpert user because of the simplicity of geometry and of the code.
LET dose total and LET track total
N = number of a certain type of ion ; n = number of different type of ion
Also, in both examples, it is possible obtaining as output the contribute of the different type of ions.
references: [7]
RBE
Input --> tables that contain data about cellular lines and different incident ions:
The code interpolates for each particles and find αi and βi
Than calculates and (mean values for each slice)
RBE and Survival Fraction can be calculated:
if D ≤ Dcut
if D > Dcut
Comparison between Hadrontherapy and Extended Example
macro Hadrontherapy
macro Extended Example
Dose, LETD and LETT
Comparison between Hadrontherapy and Extended Example
macro Hadrontherapy
macro Extended Example
Dose, LETD and LETT
Comparison between Hadrontherapy and Extended Example
Lem1_ARPE19
macro Hadrontherapy
macro Extended Example
RBE
Comparison between Hadrontherapy and Extended Example
Lem2_ARPE19
macro Hadrontherapy
macro Extended Example
RBE
Comparison between Extended Example and experimental data
Results for cuts: 0.1 mm, 0.01 mm, 0.001 mm.
DOSE
Comparison between Extended Example and experimental data
LETD
Comparison between Extended Example and experimental data
LETT
Comparison between Extended Example and experimental data
LETD
Comparison between Extended Example and experimental data
LETT
Comparison between Extended Example and experimental data
DOSE
Comparison between Extended Example and experimental data
LETD
Comparison between Extended Example and experimental data
LETT
References
[1] Cirrone GAP, Cuttone G, Raffaele L et al. , “Clinical and Research Activities at the CATANA Facility of INFN-LNS: From the Conventional Hadrontherapy to the Laser-Driven Approach”. Front. Oncol. 7:223 (2017)
[2] https://www.protominternational.com/2018/05/history-of-proton-therapy/
[3] G.Petringa, F.Romano, L.Manti et al. , “Radiobiological quantities in proton therapy: Estimation and validation using Geant4-based Monte Carlo simulation”, Physica Medica , Volume 58, February 2019, Pages 72-80
[4] G.A.P. Cirrone, G. Cuttone, S.E. Mazzaglia et al. , “Hadrontherapy: a Geant4-Based Tool for Proton/Ion-Therapy Studies”, Progress in NUCLEAR SCIENCE and TECHNOLOGY, Vol. 2, pp.207-212 (2011)
[5] https://geant4-userdoc.web.cern.ch/UsersGuides/IntroductionToGeant4/html/IntroductionToG4.html
[6] G. Cuttone, G.A.P. Cirrone, G. Di Franco et al. , “CATANA protontherapy facility: The state of art of clinical and dosimetric experience” , Eur. Phys. J. Plus (2011) 126: 65
[7] G.Petringa, L. Pandola, S. Agosteo et al. , “Monte Carlo implementation of new algorithms for the evaluation of averaged-dose and -track linear energy transfers in 62 MeV clinical proton beams”, Phys. Med. Biol. 65 235043 (2020)
THANK YOU