1 of 31

A Geant4 Extended Example

for the calculation of relevant radiobiological quantities

Giuliana Navarra - 05/11/2021

2 of 31

Overview

  • Hadrontherapy and radiobiological quantities
  • The advanced example “Hadrontherapy”
  • A new example: the “Extended Example”
  • Comparison between Hadrontherapy and Extended Example
  • Validation: Extended Example vs. Experimental results

3 of 31

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]

4 of 31

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]

5 of 31

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:

  • high energy physics
  • astrophysics
  • medical physics

  • space science
  • radiation protection

references: [5]

6 of 31

All aspects of the simulation process have been included in the toolkit:

  • the geometry of the system,
  • the materials involved,
  • the fundamental particles of interest,
  • the generation of primary events,
  • the tracking of particles through materials and electromagnetic fields,
  • the physics processes governing particle interactions,
  • the response of sensitive detector components,
  • the generation of event data,
  • the storage of events and tracks,
  • the visualization of the detector and particle trajectories,
  • the capture and analysis of simulation data at different levels of detail and refinement.

references: [5]

7 of 31

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.

8 of 31

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]

9 of 31

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]

10 of 31

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.

11 of 31

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]

12 of 31

RBE

Input --> tables that contain data about cellular lines and different incident ions:

  • α and β values in function of kinetic energies of particles
  • values Dcut , αx and βx that depends on the cell type

The code interpolates for each particles and find αi and βi

Than calculates and (mean values for each slice)

13 of 31

RBE and Survival Fraction can be calculated:

if D Dcut

if D > Dcut

14 of 31

Comparison between Hadrontherapy and Extended Example

  • Source: 60 MeV protons
  • Phantom: 4x4x4 cm3
  • Voxel 0.1 mm (x direction)
  • n. of particles 1000000
  • Physics: Hadrontherapy_1
  • Material: G4_WATER

macro Hadrontherapy

macro Extended Example

15 of 31

Dose, LETD and LETT

16 of 31

Comparison between Hadrontherapy and Extended Example

  • Source: 150 MeV protons
  • Phantom: 20x20x20 cm3
  • Voxel 0.1 mm (x direction)
  • n. of particles 1000000
  • Physics: Hadrontherapy_1
  • Material: G4_WATER

macro Hadrontherapy

macro Extended Example

17 of 31

Dose, LETD and LETT

18 of 31

Comparison between Hadrontherapy and Extended Example

  • Source: 62 MeV protons
  • Phantom: 4x4x4 cm3
  • Voxel 0.1 mm (x direction)
  • n. of particles 1000000
  • Physics: Hadrontherapy_1
  • Material: G4_WATER
  • Table used for RBE:

Lem1_ARPE19

macro Hadrontherapy

macro Extended Example

19 of 31

RBE

20 of 31

Comparison between Hadrontherapy and Extended Example

  • Source: 62 MeV protons
  • Phantom: 4x4x4 cm3
  • Voxel 0.1 mm (x direction)
  • n. of particles 1000000
  • Physics: Hadrontherapy_1
  • Material: G4_WATER
  • Table used for RBE:

Lem2_ARPE19

macro Hadrontherapy

macro Extended Example

21 of 31

RBE

22 of 31

Comparison between Extended Example and experimental data

  • Source: 58.8 土 0.3 MeV protons
  • Phantom: 4x4x4 cm3
  • Voxel 0.1 mm (x direction)
  • n. of particles 10000
  • Physics: Hadrontherapy_1
  • Material: G4_WATER

Results for cuts: 0.1 mm, 0.01 mm, 0.001 mm.

DOSE

23 of 31

Comparison between Extended Example and experimental data

  • Source: 58.8 土 0.3 MeV protons
  • Phantom: 4x4x4 cm3
  • Voxel 0.1 mm (x direction)
  • n. of particles 10000
  • Physics: Hadrontherapy_1
  • Material: G4_WATER
  • Cut: 0.01 mm

LETD

24 of 31

Comparison between Extended Example and experimental data

  • Source: 58.8 土 0.3 MeV protons
  • Phantom: 4x4x4 cm3
  • Voxel 0.1 mm (x direction)
  • n. of particles 10000
  • Physics: Hadrontherapy_1
  • Material: G4_WATER
  • Cut: 0.01 mm

LETT

25 of 31

Comparison between Extended Example and experimental data

  • Source: 58.8 土 0.3 MeV protons
  • Phantom: 4x4x4 cm3
  • Voxel 0.1 mm (x direction)
  • n. of particles 10000
  • Physics: Hadrontherapy_1
  • Material: G4_WATER
  • Cut: 0.001 mm

LETD

26 of 31

Comparison between Extended Example and experimental data

  • Source: 58.8 土 0.3 MeV protons
  • Phantom: 4x4x4 cm3
  • Voxel 0.1 mm (x direction)
  • n. of particles 10000
  • Physics: Hadrontherapy_1
  • Material: G4_WATER
  • Cut: 0.001 mm

LETT

27 of 31

Comparison between Extended Example and experimental data

  • Source: 58.8 土 0.3 MeV protons
  • Phantom: 4x4x4 cm3
  • Voxel 0.1 mm (x direction)
  • n. of particles 100000
  • Physics: Hadrontherapy_1
  • Material: G4_WATER
  • Cut: 0.01 mm

DOSE

28 of 31

Comparison between Extended Example and experimental data

  • Source: 58.8 土 0.3 MeV protons
  • Phantom: 4x4x4 cm3
  • Voxel 0.1 mm (x direction)
  • n. of particles 100000
  • Physics: Hadrontherapy_1
  • Material: G4_WATER
  • Cut: 0.01 mm

LETD

29 of 31

Comparison between Extended Example and experimental data

  • Source: 58.8 土 0.3 MeV protons
  • Phantom: 4x4x4 cm3
  • Voxel 0.1 mm (x direction)
  • n. of particles 100000
  • Physics: Hadrontherapy_1
  • Material: G4_WATER
  • Cut: 0.01 mm

LETT

30 of 31

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)

31 of 31

THANK YOU