Version:
[ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]
1
2 ==========================================================================================
3 Geant4 - GAMMAKNIFE example
4 ==========================================================================================
5
6 README file
7 --------------------------------------------
8
9 AUTHORS:
10 F.Romano* (a)
11
12 PAST AUTHORS:
13 J. Pipek (c), A. Varisano (b), G.Russo (e), G.A.P. Cirrone (b), M.Russo (e), G. Cuttone (b), M.G.Sabini (d)
14
15 (a) INFN - Sezione di Catania, Catania, Italy
16 (b) INFN - Laboratori Nazionali del Sud, Catania, Italy
17 (c) Faculty of Nuclear Sciences and Physical Engineering, Czech Technical University, Czech Republic
18 (d) UOS Fisica Sanitaria, Ospedale Cannizzaro, Catania, Italy
19 (e) Fondazione Istituto San Raffaele G.Giglio, Cefalù (Palermo), Italy
20
21
22 *Corresponding author, email to francesco.romano@ct.infn.it
23
24 -------------------------------------------------------------------------------------------------
25
26
27 ----> INTRODUCTION.
28
29 The GAMMAKNIFE example is a Geant4 application specifically developed to simulate an advanced device for Stereotactic Radiosurgery. In particular, the example reproduces a Leksell Gamma-Knife unit model C, currently active at the Cannizzaro Hospital in Catania. Gamma-Knife Stereotactic Radiosurgery allows in a single session accurate external irradiations with high dose and sharp dose gradient to minimize radiation doses given to adjacent critical brain structures. The unit consists of a cast iron body containing 201 Co60 sources in a hemispherical arrangement. Beams of gamma radiation, coming from different directions focus through a collimation system at the isocentre, where the target volume is positioned.
30
31 This example allows to simulate the main components of a Gamma-Knife unit: the source, collimation system and also a spherical phantom used for routinely measurements. All the four available final collimators are reproduced and its dimension can be changed run-time by a macro command. The energy deposited distribution is computed at the end of the simulation in a voxellized geometry, using scoring mesh. GAMMAKNIFE example exploits original solutions for the sources device simulation, providing a very simple configuration of the geometry and primary events: a single source is, indeed, simulated (instead of 201) and the scoring mesh (placed in a spherical water phantom) is rotated at 201 angles correspondent to the real angular displacement of the sources in the device. This is possible thanks to the symmetry of the source positions respect to the isocentre.
32
33 The GAMMAKNIFE example has been recently migrated to multi-thread. Thus, according to the number of cores, you can choose the number of threads using the following command:
34 /run/numberOfThreads <nThread>
35
36 The GAMMAKNIFE example contains the following subfolders and files:
37 - /src: where source .cc files are stored
38 - /include: where header .hh files are stored
39 - README file
40 - MachineAngle.in, where the angular positions of the specific unit model are tabulated
41 - defaultMacro.mac, where visualization is included and a simple run with few events is done
42 - simulation.mac, where a whole simulation is performed, changing all the 201 angular positions, storing the deposited energy for each run and saving the total energy deposted at the end of the simulation.
43
44
45 ----> INSTALLATION
46
47 GAMMAKNIFE source is currently released inside the official distribution of the Geant4 toolkit in the $G4INSTALL/examples/advanced folder.
48
49 To run GAMMAKNIFE you must first install the Geant4 package. Once Geant4 is installed, the example must be first compiled. When compilation has been completed the program can be executed.
50
51 A complete guide for the Geant4 installation in different operating systems can be found inside the official installation Geant4 pages.
52
53 If some troubles are faced with the Geant4 installation please send an e-mail to us.
54
55
56 ----> HOW TO RUN THE EXAMPLE
57
58 The defaultMacro.mac is executed if no macro file is specified.
59
60 For a realistic simulation of the Gamma-Knife device, including the irradiation of all the sources, the macro simulation.mac has to be used, defining the number of particles shooted per each angle (that is, for each "source"). The following commands have to be used to this aim:
61 /gammaknife/loadAngles MachineAngle.in --> rotation agles file is called
62 /gammaknife/beamOn <number-of-events> --> <number-of-events> are shooted for each angle
63
64
65 ----> PRIMARY EVENTS GENERATION
66
67 The primary particle beam parameter are:
68 - energy of gamma: 1.17 and 1.33 MeV with the same probability
69 - X Y Z positions sampled in a cylinder 20 mm lenght and 1 mm of radius
70 - angular biasing within an angle of 3 degree: it has been checked not to influence final results (respect to a 4p isotropic sampling) drastically decreasing computatuin time
71 GPS (General Particle Source) class is used to generate the primary events.The GPS.in macro is called and executed when yoo use bot simulation.mac and defaultMacro.mac macro.
72
73
74 ----> GEOMETRICAL SET-UP.
75
76 The Gamma-Knife unit consists of a cast-iron body containing the 60Co sources and the collimation system. Sources are displaced along five parallel rings at the same distance from the isocentre. Depending on the specific ring, sources are located with a different spacing, respectively with an angle of 7.6°, 8°, 9°, 9° and 10°. These angles represent the rotation angles of the scoring mesh (contained in the MachineAngle.in file). Gamma radiation (1.17 and 1.33 MeV), emitted by Co60 sources arranged in a hemispherical shape, focuses on the isocentre, located 401 mm far from the source. The radiation beam is collimated by means of fixed and mobile collimation systems in order to obtain four possible beam diameters at the isocentre: 4, 8, 14 and 18 mm.
77
78 The elements simulated are:
79
80 - cobalt cylinder source (20 mm length, 1 mm radius)
81 - a stainless steel covering of the source
82 - a first tungsten cylindrical collimator (65 mm length)
83 - a second truncate cone lead collimator (92.5 mm lenght)
84 - final truncated cone tungsten collimator (60 mm lenght), with four available final apertures simulating the mobile helmet
85 - tungsten, iron and stainless steel surrounding the collimation system for shielding
86 - a spherical water phantom (8 mm radius) simulating the phantom typically used for dosimetrical measurements and routinely quality check
87
88 To change the diameter of the final truncated collimator (mobile helmet) the following command has to be used for chaching the diameter:
89 /calorimeter/helmetSize <diameter in mm>
90 Notice that, as in the real device, only four final diameters are allowed: 4, 8, 14 and 18 mm. These are not the physical apertures of the diameter but represent the final diameter of the beam at the isocentre. In the GammaKnifeDetectorConstruction::UpdateHelmet() method more details can be found in this concern.
91
92 Important note: the GAMMAKNIFE application is meant as an example of simulation of a device for Stereotactic Radiosurgery. Thus, all the main geometrical features of the Gamma-Knife device are accurately simulated but the geometrical displecement of the some elements is not exactly correspondent to the real one, because covered by a patent. In case you need detailed geometrical information, please contact the Elekta Company.
93
94
95 ----> PHYSICS
96
97 Local PhysicsList is activated by default (GammaKnifePhysicsList.cc), which includes the G4EmStandardPhysics_option3.
98
99 Anyhow, it is recommended the use of the Reference Physics List QGSP_BIC_EMY, which is the suggested on for medical applications.
100
101 ----> VISUALISATION
102
103 The user can visualise the experimental set-up with OpenGL, DAWN and vrml.
104
105
106 ----> SIMULATION OUTPUT
107
108 At the end of the simulation the energy deposited in each voxel of the scoring mesh inside the spherical phantom is registered in the eDep_scorer.out file. Default voxelization is 45 x 45 x 45 cubic voxels 1 mm of side. User can change dimension and shape of box the mesh and the voxel.
109 Alternatively, the dose deposited can be also stored using the appropriate scoring commands
110
111 ---------------------------------------------------------------------------
112
113
114 For comments, advices and questions please contact:
115 francesco.romano@ct.infn.it
116
117 last modified: F.Romano, 03/11/2020
118
119