Dosimetry in magnetic fields with dedicated MR-compatible ionization chambers

MR-integrated radiotherapy requires suitable dosimetry detectors to be used in magnetic fields. This study investigates the feasibility of using dedicated MR-compatible ionization chambers at MR-integrated radiotherapy devices. MR-compatible ionization chambers (Exradin A19MR, A1SLMR, A26MR, A28MR) were precisely modeled and their relative response in a 6MV treatment beam in the presence of a magnetic field was simulated using EGSnrc. Monte Carlo simulations were carried out with the magnetic field in three orientations: the magnetic field aligned perpendicular to the chamber and beam axis (transverse orientation), the magnetic field parallel to the chamber as well as parallel to the beam axis. Monte Carlo simulation results were validated with measurements using an electromagnet with magnetic field strength upto 1.1 T with the chambers in transverse orientation. The measurements and simulation results were in good agreement, except for the A26MR ionization chamber in transverse orientation. The maximum increase in response of the ionization chambers observed was 8.6% for the transverse orientation. No appreciable change in chamber response due to the magnetic field was observed for the magnetic field parallel to the ionization chamber and parallel to the photon beam.Polarity and recombination correction factor were experimentally investigated in the transverse orientation. The polarity effect and recombination effect were not altered by a magnetic field.This study further investigates the response of the ionization chambers as a function of the chambers’ rotation around their longitudinal axis. A variation in response was observed when the chamber was not rotationally symmetric, which was independent of the magnetic field.     

High throughput film dosimetry in homogeneous and heterogeneous media for a small animal irradiator

PurposeWe have established a high-throughput Gafchromic film dosimetry protocol for narrow kilovoltage beams in homogeneous and heterogeneous media for small-animal radiotherapy applications. The kV beam characterization is based on extensive Gafchromic film dosimetry data acquired in homogeneous and heterogeneous media. An empirical model is used for parameterization of depth and off-axis dependence of measured data.MethodsWe have modified previously published methods of film dosimetry to suit the specific tasks of the study. Unlike film protocols used in previous studies, our protocol employs simultaneous multi-channel scanning and analysis of up to nine Gafchromic films per scan. A scanner and background correction were implemented to improve accuracy of the measurements. Measurements were taken in homogeneous and inhomogeneous phantoms at 220 kVp and a field size of 5 × 5 mm². The results were compared against Monte Carlo simulations.ResultsDose differences caused by variations in background signal were effectively removed by the corrections applied. Measurements in homogeneous phantoms were used to empirically characterize beam data in homogeneous and heterogeneous media. Film measurements in inhomogeneous phantoms and their empirical parameterization differed by about 2%–3%. The model differed from MC by about 1% (water, lung) to 7% (bone). Good agreement was found for measured and modelled off-axis ratios.ConclusionsEBT2 films are a valuable tool for characterization of narrow kV beams, though care must be taken to eliminate disturbances caused by varying background signals. The usefulness of the empirical beam model in interpretation and parameterization of film data was demonstrated.     

On the use of EBT3 film for relative dosimetry of kilovoltage X ray beams

EBT3 films were evaluated for relative dosimetry in water, in the energy range of therapeutic kV X ray beams. A film batch was calibrated in air for all nine beam qualities of a clinical unit (XStrahl 200). Monte Carlo (MC) simulations using MCNP v.6 facilitated the calculation of the film absorbed dose (f), and beam quality (k_bq) energy dependences in air. Results were found in agreement with corresponding data in the literature. Film samples from the same batch were irradiated in water along the central beam axis for each beam quality. Experimental percentage depth dose (PDD) results obtained using calibration data in air showed quality and depth dependent differences from corresponding MC simulations. These differences increased beyond film dosimetry uncertainty (<3.3%), reaching up to 8% at increased depth. The observed differences reduced only slightly when spectral variation as a function of measurement point was accounted for, using photon effective energy. PDD measurements and corresponding MC results facilitated the determination of f and k_bq in water. Results showed that the origin of the observed differences between experimental and MC PDD results is the difference between film response in air and water, as a result of radiation field perturbation from the film oriented along the central beam axis. This implies a directional dependence of film response which necessitates that the angular distribution of photons impinging on the film is the same in the calibration and measurement geometries.     

A dosimetry method in the transverse plane of HDR Ir-192 brachytherapy source using gafchromic EBT2 film

Radiochromic film dosimetry is increasingly used in brachytherapy applications for its higher resolution ability as compared to other experimental methods. The present study was aimed to assess the accuracy and suitability of use of the improved radiochromic film model, Gafchromic EBT2, to evaluate the dose distribution in the transverse plane of microselectron HDR ¹⁹²Ir source.A specially designed and locally fabricated Polymethyl methacrylate (PMMA) phantom was used in this work for the experimental measurement of dose distribution around the source in its transverse plane. The AAPM TG-43U1 recommended radial dose function, g (r), and dose rate constant, Λ, for the source were measured using Gafchromic EBT2 film and thermoluminescent dosimeters (TLD). The EBT2 film measured dosimetric quantities were validated against their values obtained from the TLD measurements and previously published values for the same source available in literature.The dose rate constant and radial dose function for microselectron HDR ¹⁹²Ir source obtained from Gafchromic EBT2 film measurements are in agreement with their TLD measured results within 3.9% and 2.8% respectively. They also agree within the accepted range of uncertainty with their experimental and Monte Carlo calculated results reported in literature.This work demonstrates the suitability of using Gafchromic EBT2 film dosimetry in characterization of dose distribution in the transverse plane of HDR Ir-192 source. This is a more efficient method than TLD dosimetry at discrete and distant positions. Relative to TLD dosimetry, it is found to be better reproducible, easy to use and a less expensive method of dosimetry.     

Correction of measured Gamma-Knife output factors for angular dependence of diode detectors and PinPoint ionization chamber

Dosimetry for Gamma-Knife requires detectors with high spatial resolution and minimal angular dependence of response. Angular dependence and end effect time for p-type silicon detectors (PTW Diode P and Diode E) and PTW PinPoint ionization chamber were measured with Gamma-Knife beams. Weighted angular dependence correction factors were calculated for each detector. The Gamma-Knife output factors were corrected for angular dependence and end effect time. For Gamma-Knife beams angle range of 84°–54°. Diode P shows considerable angular dependence of 9% and 8% for the 18 mm and 14, 8, 4 mm collimator, respectively. For Diode E this dependence is about 4% for all collimators. PinPoint ionization chamber shows angular dependence of less than 3% for 18, 14 and 8 mm helmet and 10% for 4 mm collimator due to volumetric averaging effect in a small photon beam. Corrected output factors for 14 mm helmet are in very good agreement (within ±0.3%) with published data and values recommended by vendor (Elekta AB, Stockholm, Sweden). For the 8 mm collimator diodes are still in good agreement with recommended values (within ±0.6%), while PinPoint gives 3% less value. For the 4 mm helmet Diodes P and E show over-response of 2.8% and 1.8%, respectively. For PinPoint chamber output factor of 4 mm collimator is 25% lower than Elekta value which is generally not consequence of angular dependence, but of volumetric averaging effect and lack of lateral electronic equilibrium. Diodes P and E represent good choice for Gamma-Knife dosimetry.     

Application of Gafchromic EBT2 film for intraoperative radiation therapy quality assurance

PurposeIntraoperative radiation therapy (IORT) using electron beam is commonly done by mobile dedicated linacs that have a variable range of electron energies. This paper focuses on the evaluation of the EBT2 film response in the green and red colour channels for IORT quality assurance (QA).MethodsThe calibration of the EBT2 films was done in two ranges; 0–8 Gy for machine QA by red channel and 8–24 Gy for patient-specific QA by green channel analysis. Irradiation of calibration films and relative dosimetries were performed in a water phantom. To evaluate the accuracy of the film response in relative dosimetry, gamma analysis was used to compare the results of the Monte Carlo simulation and ionometric dosimetry. Ten patients with early stage breast cancer were selected for in-vivo dosimetry using the green channel of the EBT2 film.ResultsThe calibration curves were obtained by linear fitting of the green channel and a third-order polynomial function in the red channel (R2 = 0.99). The total dose uncertainty was up to 4.2% and 4.7% for the red and green channels, respectively. There was a good agreement between the relative dosimetries of films by the red channel, Monte Carlo simulations and ionometric values. The mean dose difference of the in-vivo dosimetry by green channel of this film and the expected values was about 1.98% ± 0.75.ConclusionThe results of this study showed that EBT2 film can be considered as an appropriate tool for machine and patient-specific QA in IORT.     

Evaluation of Gafchromic EBT2 film for the measurement of anisotropy function for high-dose-rate 192Ir brachytherapy source with respect to thermoluminescent dosimetry

Aim

The aim of this work was to assess the suitability of the use of a Gafchromic EBT2 film for the measurement of anisotropy function for microSelectron HDR ¹⁹²Ir (classic) source with a comparative dosimetry method using a Gafchromic EBT2 film and thermoluminescence dosimeters (TLDs).

Background

Sealed linear radiation sources are commonly used for high dose rate (HDR) brachytherapy treatments. Due to self-absorption and oblique filtration of radiation in the source capsule material, an inherent anisotropy is present in the dose distribution around the source which can be described by a measurable two-dimensional anisotropy function, F(r, θ).

Materials and methods

Measurements were carried out in a specially designed and locally fabricated PMMA phantom with provisions to accommodate miniature LiF TLD rods and EBT2 film dosimeters at identical radial distances with respect to the ¹⁹²Ir source.

Results

The data of anisotropy function generated by the use of the Gafchromic EBT2 film method are in agreement with their TLD measured values within 4%. The produced data are also consistent with their experimental and Monte Carlo calculated results for this source available in the literature.

Conclusion

Gafchromic EBT2 film was found to be a feasible dosimeter in determining anisotropy in the dose distribution of ¹⁹²Ir source. It offers high resolution and is a viable alternative to TLD dosimetry at discrete points. The method described in this paper is useful for comparing the performances of detectors and can be applied for other brachytherapy sources as well.     

Comparative analysis of several detectors for the measurement of radiation transmission and leakage from a multileaf collimator

The multileaf collimator (MLC) is the standard device used to shape radiation beams for 3-d conformal and intensity-modulated radiation therapy (IMRT). Due to the inherent properties of MLC, there is a small amount of radiation transmitted through the leaves, called radiation transmission (RT). Accurate measurements of this radiation are required to commission and validate IMRT-capable treatment planning systems because this radiation may impact the dosimetry of IMRT-calculated dose distributions. This work compares several detectors in the measurement of RT for a micro-multileaf collimation system. The results show that there are statistically significant differences in the measured RT values between detectors from 3.5 to 12.5% for the same MLC model and less than 0.2% relative to the isocentre dose for an open reference field. However, although small in magnitude, these differences may impact the dosimetry of IMRT treatment planning by up to 1.78 Gy to the healthy tissue surrounding the target for a treatment of 60 Gy in 30 fractions. By the later, these differences must be included as a source of uncertainty in IMRT dose delivery. Also, it must be established which detector offers the most reliable results in the measurement of the RT by using Monte Carlo simulation methods.     

Comparison of dosimetric characteristics of physical wedge and enhanced dynamic wedge in inhomogeneous medium using Monte Carlo simulations

BackgroundWidely used physical wedges in clinical radiotherapy lead to beam intensity attenuation as well as the beam hardening effect, which must be considered. Dynamic wedges devised to overcome the physical wedges (PWs) problems result in dosimetry complications due to jaw movement while the beam is on. This study was aimed to investigate the usability of physical wedge data instead of enhanced dynamic wedge due to the enhanced dynamic wedge (EDW) dosimetry measurement hardships of Varian 2100CD in inhomogeneous phantom by Monte Carlo code as a reliable method in radiation dosimetry.Materials and methodsA PW and EDW-equipped-linac head was simulated using BEAMnrc code. DOSXYZnrc was used for three-dimensional dosimetry calculation in the CIRS phantom.ResultsBased on the isodose curves, EDW generated a less scattered as well as lower penumbra width compared to the PW. The depth dose variations of PWs and EDWs were more in soft tissue than the lung tissue. Beam profiles of PW and EDW indicated good coincidence in all points, except for the heel area.ConclusionResults demonstrated that it is possible to apply PW data instead of EDW due to the dosimetry and commissioning hardships caused by EDW in inhomogeneous media.     

Towards MR-guided electron therapy: Measurement and simulation of clinical electron beams in magnetic fields

PurposeIn the current era of MRI-linac radiotherapy, dose optimization with arbitrary dose distributions is a reality. For the first time, we present new and targeted experiments and modeling to aid in evaluating the potential dose improvements offered with an electron beam mode during MRI-linac radiotherapy.MethodsSmall collimated (1 cm diameter and 1.5 × 1.5 cm² square) electron beams (6, 12 and 20 MeV) from a clinical linear accelerator (Varian Clinac 2100C) are incident perpendicular and parallel to the strong and localized magnetic fields (0–0.7 T) generated by a permanent magnet device. Gafchromic EBT3 film is placed inside a slab phantom to measure two-dimensional dose distributions. A benchmarked and comprehensive Monte Carlo model (Geant4) is established to directly compare with experiments.ResultsWith perpendicular fields a 5% narrowing of the beam FWHM and a 10 mm reduction in the 15% isodose penetration is seen for the 20 MeV beam. In the inline setup the penumbral width is reduced by up to 20%, and a local central dose enhancement of 100% is observed. Monte Carlo simulations are in agreement with the measured dose distributions (2% or 2 mm).ConclusionA new range of experiments have been performed to offer insight into how an electron beam mode could offer additional choices in MRI-linac radiotherapy. The work extends on historic studies to bring a successful unified experimental and Monte Carlo modeling approach for studying small field electron beam dosimetry inside magnetic fields. The results suggest further work, particularly on the inline magnetic field scenario.     

Calibration of MTT assay in proton beams using radiochromic films

PurposeThis study provides methodology of calibrating as well as controlling the output for an MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide) colorimetric assay irradiated in a low energy proton beam using EBT3-model GAFCHROMIC^TM film, without correcting for quenching effect.MethodsA calibrated Markus ionization chamber was used to measure the depth dose and beam output for 26.5 MeV protons produced by a CS30 cyclotron. A time-controlled aluminum cylinder was added in front of the horizontal beam-exit serving as a radiation shutter. Following the TRS-398 reference dosimetry protocol for proton beams, the output was calibrated in water at a reference depth of 3 mm. EBT3 film was calibrated for doses up to 8 Gy at the same depth. To verify the dose distribution for each 96-well MTT assay plate, EBT3 film was placed at the reference depth during irradiation and cell doses were scaled by measured percent depth dose (PDD) data.ResultsThe radiochromic film dosimetry system in this study provides dose measurements with an uncertainty better than 3.3% for doses higher than 1 Gy. From a single exposure and utilizing the Gaussian shape of the beam, multiple dose points can be obtained within different wells of the same plate ranging from 6.9 Gy (sigma ∼4%) in the central well, and 2 Gy (sigma ∼8%) for wells positioned closer to the periphery.ConclusionsWe described a methodology for radiochromic film-based dose monitoring system, using low-energy protons, which can be used for the MTT assay in any proton beam, except within Bragg peak region.     

Multiple-source models for electron beams of a medical linear accelerator using BEAMDP computer code

AimThe aim of this work was to develop multiple-source models for electron beams of the NEPTUN 10PC medical linear accelerator using the BEAMDP computer code.BackgroundOne of the most accurate techniques of radiotherapy dose calculation is the Monte Carlo (MC) simulation of radiation transport, which requires detailed information of the beam in the form of a phase-space file. The computing time required to simulate the beam data and obtain phase-space files from a clinical accelerator is significant. Calculation of dose distributions using multiple-source models is an alternative method to phase-space data as direct input to the dose calculation system.Materials and methodsMonte Carlo simulation of accelerator head was done in which a record was kept of the particle phase-space regarding the details of the particle history. Multiple-source models were built from the phase-space files of Monte Carlo simulations. These simplified beam models were used to generate Monte Carlo dose calculations and to compare those calculations with phase-space data for electron beams.ResultsComparison of the measured and calculated dose distributions using the phase-space files and multiple-source models for three electron beam energies showed that the measured and calculated values match well each other throughout the curves.ConclusionIt was found that dose distributions calculated using both the multiple-source models and the phase-space data agree within 1.3%, demonstrating that the models can be used for dosimetry research purposes and dose calculations in radiotherapy.     

Measuring the dose in bone for spine stereotactic body radiotherapy

PurposeCurrent quality assurance of radiotherapy involving bony regions generally utilises homogeneous phantoms and dose calculations, ignoring the challenges of heterogeneities with dosimetry problems likely occurring around bone. Anthropomorphic phantoms with synthetic bony materials enable realistic end-to-end testing in clinical scenarios. This work reports on measurements and calculated corrections required to directly report dose in bony materials in the context of comprehensive end-to-end dosimetry audit measurements (63 plans, 6 planning systems).Materials and methodsRadiochromic film and microDiamond measurements were performed in an anthropomorphic spine phantom containing bone equivalent materials. Medium dependent correction factors, k_med, were established using 6 MV and 10 MV Linear Accelerator Monte Carlo simulations to account for the detectors being calibrated in water, but measuring in regions of bony material. Both cortical and trabecular bony material were investigated for verification of dose calculations in dose-to-medium (D_m,m) and dose-to-water (D_w,w) scenarios.ResultsFor D_m,m calculations, modelled correction factors for cortical and trabecular bone in film measurements, and for trabecular bone in microDiamond measurements were 0.875(±0.1%), 0.953(±0.3%) and 0.962(±0.4%), respectively. For D_w,w calculations, the corrections were 0.920(±0.1%), 0.982(±0.3%) and 0.993(±0.4%), respectively. In the audit, application of the correction factors improves the mean agreement between treatment plans and measured microDiamond dose from −2.4%(±3.9%) to 0.4%(±3.7%).ConclusionMonte Carlo simulations provide a method for correcting the dose measured in bony materials allowing more accurate comparison with treatment planning system doses. In verification measurements, algorithm specific correction factors should be applied to account for variations in bony material for calculations based on D_m,m and D_w,w.     

Integration of the M6 Cyberknife in the Moderato Monte Carlo platform and prediction of beam parameters using machine learning

PurposeThis work describes the integration of the M6 Cyberknife in the Moderato Monte Carlo platform, and introduces a machine learning method to accelerate the modelling of a linac.MethodsThe MLC-equipped M6 Cyberknife was modelled and integrated in Moderato, our in-house platform offering independent verification of radiotherapy dose distributions. The model was validated by comparing TPS dose distributions with Moderato and by film measurements. Using this model, a machine learning algorithm was trained to find electron beam parameters for other M6 devices, by simulating dose curves with varying spot size and energy. The algorithm was optimized using cross-validation and tested with measurements from other institutions equipped with a M6 Cyberknife.ResultsOptimal agreement in the Monte Carlo model was reached for a monoenergetic electron beam of 6.75 MeV with Gaussian spatial distribution of 2.4 mm FWHM. Clinical plan dose distributions from Moderato agreed within 2% with the TPS, and film measurements confirmed the accuracy of the model. Cross-validation of the prediction algorithm produced mean absolute errors of 0.1 MeV and 0.3 mm for beam energy and spot size respectively. Prediction-based simulated dose curves for other centres agreed within 3% with measurements, except for one device where differences up to 6% were detected.ConclusionsThe M6 Cyberknife was integrated in Moderato and validated through dose re-calculations and film measurements. The prediction algorithm was successfully applied to obtain electron beam parameters for other M6 devices. This method would prove useful to speed up modelling of new machines in Monte Carlo systems.     

La simulation Monte Carlo en médecine nucléaire

The Monte Carlo method allows for simulating random processes by using series of pseudo-random numbers. It became an important tool in nuclear medicine to assist in the design of new medical imaging devices, optimise their use and analyse their data. Presently, the sophistication of the simulation tools allows the introduction of Monte Carlo predictions in data correction and image reconstruction processes. The availability to simulate time dependent processes opens up new horizons for Monte Carlo simulation in nuclear medicine. In a near future, these developments will allow to tackle simultaneously imaging and dosimetry issues and soon, case system Monte Carlo simulations may become part of the nuclear medicine diagnostic process. This paper describes some Monte Carlo method basics and the sampling methods that were developed for it. It gives a referenced list of different simulation software used in nuclear medicine and enumerates some of their present and prospective applications.     

Small animal irradiator dose distribution verification using radioluminescence imaging

The main objective of this work was the development of a novel 2D dosimetry approach for small animal external radiotherapy using radioluminescence imaging (RLI) with a commercial complementary metal oxide semiconductor detector. Measurements of RLI were performed on the small animal image‐guided platform SmART, RLI data were corrected for perspective distortion using Matlab. Four irradiation fields were tested and the planar 2D dose distributions and dose profiles were compared against dose calculations performed with a Monte Carlo based treatment planning system and gafchromic film. System linearity and RLI image noise against dose were also measured. The maximum difference between beam size measured with RLI and nominal beam size was less than 8% for all the tested beams. The image correction procedure was able to reduce perspective distortion. A novel RLI approach for quality assurance of a small animal irradiator was presented and tested. Results are in agreement with MC dose calculations and gafchromic film measurements.     

Field output factors for small fields: A large multicentre study

PurposeThe determination of output factors in small field dosimetry is a crucial point, especially when implementing stereotactic radiotherapy (SRT). Herein, a working group of the French medical physicist society (SFPM) was created to collect small field output factors. The objective was to gather and disseminate information on small field output factors based on different detectors for various clinical SRT equipment and measurement configurations.MethodParticipants were surveyed for information about their SRT equipment, including the type of linear particle accelerator (linac), collimator settings, measurement conditions for the output factors and the detectors used. Participants had to report both the ratio of detector readings and the correction factors applied as described in the IAEA TRS-483 code of practice for nominal field sizes smaller or equal to 3 cm. Mean field output factors and their associated standard deviations were calculated when data from at least 3 linacs were available.Results23 centres were enrolled in the project. Standard deviations of the mean field output factors were systematically smaller than 1.5% for field sizes larger or equal to 1 cm and reached 5% for the smallest field size (0.5 cm). Deviations with published data were smaller than 2% except for the 0.5 cm circular fixed aperture collimator of the CyberKnife where it reached 3.5%.ConclusionThese field output factor values obtained via a large multicentre study can be considered as an external cross verification for any radiotherapy centre starting a SRT program and should help minimize systematic errors when determining small field output factors.     

Assessment of ionization chamber correction factors in photon beams using a time saving strategy with PENELOPE code

The purpose of this study was to investigate Monte Carlo-based perturbation and beam quality correction factors for ionization chambers in photon beams using a saving time strategy with PENELOPE code. Simulations for calculating absorbed doses to water using full spectra of photon beams impinging the whole water phantom and those using a phase-space file previously stored around the point of interest were performed and compared. The widely used NE2571 ionization chamber was modeled with PENELOPE using data from the literature in order to calculate absorbed doses to the air cavity of the chamber. Absorbed doses to water at reference depth were also calculated for providing the perturbation and beam quality correction factors for that chamber in high energy photon beams. Results obtained in this study show that simulations with phase-space files appropriately stored can be up to ten times shorter than using a full spectrum of photon beams in the input-file. Values of k_Q and its components for the NE2571 ionization chamber showed good agreement with published values in the literature and are provided with typical statistical uncertainties of 0.2%. Comparisons to k_Q values published in current dosimetry protocols such as the AAPM TG-51 and IAEA TRS-398 showed maximum percentage differences of 0.1% and 0.6% respectively. The proposed strategy presented a significant efficiency gain and can be applied for a variety of ionization chambers and clinical photon beams.     

Clinical investigations of a CVD diamond detector for radiotherapy dosimetry

X-ray detectors based on single crystal diamond film made via chemical vapor deposition were investigated to evaluate their performance under clinically relevant conditions for radiotherapy dosimetry. Studies focused on repeatability, dose rate dependence, tissue phantom ratios, output factors and beam profiling. Repeatability experiments revealed a temporary loss in sensitivity due to charge detrapping effects following irradiation, which was modeled to make corrections that improved short-term precision. Dose rate dependence was observed (Fowler fitting parameter Δ = 0.96 ± 0.2) using dose rates up to ～2 Gy min^?1. The detector statistically distinguished (n = 5, P < 0.05) between dose values separated by 7.7 × 10^?3 Gy (1 MU). Depth dose measurements from 1 to 15 cm and output factors using 3 × 3 to 10 × 10 cm^2 field sizes compared well with a Farmer ion chamber (<1.3% difference). Output factor measurements indicate encouraging results for fields sizes <4 × 4 cm^2. Off-axis measurements showed that perturbation of the beam could be reduced when the detector is used in the edge-on orientation due to its thin-film sandwich configuration and ～200 nm thick Ag contacts. This relatively inexpensive detector has potential to be used for routine dosimetry using conventional radiotherapy instrumentation.     

Validation of Monte carlo Geant4 multithreading code for a 6 MV photon beam of varian linac on the grid computing

AimTo evaluate the computation time efficiency of the multithreaded code (G4Linac-MT) in the dosimetry application, using the high performance of the HPC-Marwan grid to determine with high accuracy the initial parameters of the 6 MV photon beam of Varian CLINAC 2100C.BackgroundThe difficulty of Monte Carlo methods is the long computation time, this is one of the disadvantages of the Monte Carlo methods.Materials and methodsCalculations are performed by the multithreaded code G4Linac-MT and Geant4.10.04.p02 using the HPC-Marwan computing grid to evaluate the computing speed for each code. The multithreaded version is tested in several CPUs to evaluate the computing speed according to the number of CPUs used. The results were compared to the measurements using different types of comparisons, TPR_20.10, penumbra, mean dose error and gamma index.ResultsThe results obtained for this work indicate a much higher computing time saving for the G4Linac-MT version compared to the Geant4.10.04 version, the computing time decreases with the number of CPUs used, can reach about 12 times if 64CPUs are used. After optimization of the initial electron beam parameters, the results of the dose simulations obtained for this work are in very good agreement with the experimental measurements with a mean dose error of up to 0.41% on the PDDs and 1.79% on the lateral dose.ConclusionsThe gain in computation time leads us to perform Monte Carlo simulations with a large number of events which gives a high accuracy of the dosimetry results obtained in this work.