F. Romano (a)
Past authors: G. Cuttone (b), J. Pipek (c), M.G.Sabini (d), G.A.P. Cirrone (b), G.Russo (e), M.Russo (b)
(a) Sezione di Catania - INFN, Catania, Italy
(b) Laboratori Nazionali del Sud - INFN, Catania, Italy
(c) Faculty of Nuclear Sciences and Physical Engineering, Czech Technical University, Czech Republic
(d) UOS Fisica Sanitaria, Ospedale Cannizzaro, Catania, Italy
(e) Fondazione Istituto San Raffaele G.Giglio, Cefalù (Palermo), Italy
*Corresponding author, email to Francesco Romano
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 Co-60 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. This example allows the simulation of 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 deposited energy distribution is computed at the end of the simulation in a voxelized geometry, using a scoring mesh.
The 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.
The GAMMAKNIFE example contains the following subfolders and files:
- /src: where source .cc files are stored
- /include: where header .hh files are stored
- README file
- MachineAngle.in, where the angular positions of the specific unit model are tabulated
- defaultMacro.mac, where visualization is included and a simple run with few events is done
- 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.
The Gamma-Knife unit consists of a cast-iron body containing the Co-60 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. The elements simulated are:
- cobalt cylinder source (20 mm length, 1 mm radius)
- a stainless steel covering of the source
- a first tungsten cylindrical collimator (65 mm length)
- a second truncate cone lead collimator (92.5 mm lenght)
- final truncated cone tungsten collimator (60 mm lenght), with four available final apertures simulating the mobile helmet
- tungsten, iron and stainless steel surrounding the collimation system for shielding
- a spherical water phantom (8 cm radius) simulating the phantom typically used for dosimetrical measurements and routinely quality check.
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.
Important note: the GAMMAKNIFE application is meant as an example of simulation of a device for Stereotactic Radiosurgery. Therefore, all the main geometrical features of the Gamma-Knife device are accurately simulated but the geometrical displacement of 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.
A local PhysicsList is activated by default (GammaKnifePhysicsList.cc), which includes the G4EmStandardPhysics_option3. Anyhow, it is recommended the use of the Reference Physics List QGSP_BIC_EMY, which is the suggested one for medical applications. It is possible to activate that defining the environment variable PHYSLIST in such a way:
The user can visualise the experimental set-up with OpenGL, DAWN and vrml.
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.
Alternatively, the dose deposited can be also stored using the appropriate scoring commands.