Monte Carlo simulation using PRIMO code as a tool for checking the credibility of commissioning and quality assurance of 6 MV TrueBeam STx varian LINAC

AimTo validate and implement Monte Carlo simulation using PRIMO code as a tool for checking the credibility of measurements in LINAC initial commissioning and routine Quality Assurance (QA). Relative and absolute doses of 6 MV photon beam from TrueBeam STx Varian Linear Accelerator (LINAC) were simulated and validated with experimental measurement, Analytical Anisotropic Algorithm (AAA) calculation, and golden beam.Methods and MaterialsVarian phase-space files were imported to the PRIMO code and four blocks of jaws were simulated to determine the field size of the photon beam. Water phantom was modeled in the PRIMO code with water equivalent density. Golden beam data, experimental measurement, and AAA calculation results were imported to PRIMO code for gamma comparison.ResultsPRIMO simulations of Percentage Depth Dose (PDD) and in-plane beam profiles had good agreement with experimental measurements, AAA calculations and golden beam. However, PRIMO simulations of cross-plane beam profiles have a better agreement with AAA calculation and golden beam than the experimental measurement. Furthermore, PRIMO simulations of absolute dose agreed well with experimental results with ±0.8% uncertainty.ConclusionThe PRIMO code has good accuracy and is appropriate for use as a tool to check the credibility of beam scanning and output measurement in initial commissioning and routine QA.     

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.     

Experimental assessment of proton dose calculation accuracy in inhomogeneous media

PurposeProton therapy with Pencil Beam Scanning (PBS) has the potential to improve radiotherapy treatments. Unfortunately, its promises are jeopardized by the sensitivity of the dose distributions to uncertainties, including dose calculation accuracy in inhomogeneous media. Monte Carlo dose engines (MC) are expected to handle heterogeneities better than analytical algorithms like the pencil-beam convolution algorithm (PBA). In this study, an experimental phantom has been devised to maximize the effect of heterogeneities and to quantify the capability of several dose engines (MC and PBA) to handle these.MethodsAn inhomogeneous phantom made of water surrounding a long insert of bone tissue substitute (1 × 10 × 10 cm³) was irradiated with a mono-energetic PBS field (10 × 10 cm²). A 2D ion chamber array (MatriXX, IBA Dosimetry GmbH) lied right behind the bone. The beam energy was such that the expected range of the protons exceeded the detector position in water and did not attain it in bone. The measurement was compared to the following engines: Geant4.9.5, PENH, MCsquare, as well as the MC and PBA algorithms of RayStation (RaySearch Laboratories AB).ResultsFor a γ-index criteria of 2%/2 mm, the passing rates are 93.8% for Geant4.9.5, 97.4% for PENH, 93.4% for MCsquare, 95.9% for RayStation MC, and 44.7% for PBA. The differences in γ-index passing rates between MC and RayStation PBA calculations can exceed 50%.ConclusionThe performance of dose calculation algorithms in highly inhomogeneous media was evaluated in a dedicated experiment. MC dose engines performed overall satisfactorily while large deviations were observed with PBA as expected.     

Monte carlo model and output factors of elekta infinity™ 6 and 10 MV photon beam

AimThis study aimed to commission the Elekta Infinity™ working in 6 and 10 MV photon beam installed in Concord International Hospital, Singapore, and compare the OFs between MC simulation and measurement using PTW semiflex and microDiamond detector for small field sizes.Material and MethodsThere are two main steps in this study: modelling of Linac 6 and 10 MV photon beam and analysis of the output factors for field size 2 × 2–10 × 10 cm². The EGSnrc/BEAMnrc-DOSXYZnrc code was used to model and characterize the Linac and to calculate the dose distributions in a water phantom. The dose distribution and OFs were compared to the measurement data in the same condition.ResultsThe commissioning process was only conducted for a 10 × 10 cm² field size. The PDD obtained from MC simulation showed a good agreement with the measurement. The local dose difference of PDDs was less than 2% for 6 and 10 MV. The initial electron energy was 5.2 and 9.4 MeV for 6 and 10 MV photon beam, respectively. This Linac model can be used for dose calculation in other situations and different field sizes because this Linac has been commissioned and validated using Monte Carlo simulation. The 10 MV Linac produces higher electron contamination than that of 6 MV.ConclusionsThe Linac model in this study was acceptable. The most important result in this work comes from OFs resulted from MC calculation. This value was more significant than the OFs from measurement using semiflex and microDiamond for all beam energy and field sizes because of the CPE phenomenon.     

Modeling of couch transmission in the RayStation treatment planning system

PurposeTo present our methods and results regarding the modeling of a carbon fiber couch (Varian Exact IGRT) in the RayStation treatment planning system (TPS).MethodsThree geometrical-models (GMs) were implemented in the TPS to represent the three different regions of the couch (thick, medium and thin). The materials and densities of each GM component were tuned to maximize the agreement between measured and calculated attenuations. Moreover, a couch computed-tomography (CT) scan was acquired and dosimetrically compared with the GMs. For validation, plan-specific quality assurance (QA) of VMAT plans (TG-119 cases, 5 prostate and 5 H&N clinical cases) was performed by comparing measured dose distributions with doses computed with and without including the GMs in the TPS.ResultsCouch attenuations up to 4.3% were measured (energy: 6MV). Compared to couch CT, GMs could be modified to optimize the agreement with measurements and reduce dependence on the dose grid resolution. For both couch CT and GM, absolute deviations between measured and calculated attenuations were within 1.0%. When including the GMs in plan-specific QA, global 2%/2 mm γ-pass rates showed an average improvement of 4.8% (p-value < 0.001, max +18.6%). The couch reduced the mean dose to targets by up to 2.4% of the prescribed dose for prostate cases and up to 1.4% for H&N cases.ConclusionsRayStation accurately considers the implemented couch GMs replicating measured attenuations within an uncertainty of 1.0%. Materials and densities are proposed for the Varian Exact IGRT couch. The results obtained justify introducing couch GMs in clinical routine.     

Monte Carlo determination of a nanoDot OSLD response using quality index for diagnostic kilovoltage X-ray beams

PurposeThis study aims to investigate the energy response of an optically stimulated luminescent dosimeter known as nanoDot for diagnostic kilovoltage X-ray beams via Monte Carlo calculations.MethodsThe nanoDot response is calculated as a function of X-ray beam quality in free air and on a water phantom surface using Monte Carlo simulations. The X-ray fluence spectra are classified using the quality index (QI), which is defined as the ratio of the effective energy to the maximum energy of the photons. The response is calculated for X-ray fluence spectra with QIs of 0.4, 0.5, and 0.6 with tube voltages of 50–137.6 kVp and monoenergetic photon beams. The surface dose estimated using the calculated response is verified by comparing it with that measured using an ionization chamber.ResultsThe nanoDot response in free air for monoenergetic photon beams (QI = 1.0) varies significantly at photon energies below 100 keV and reaches a factor of 3.6 at 25–30 keV. The response differs by up to approximately 6% between QIs of 0.4 and 0.6 for the same half-value layer (HVL). The response at the phantom surface decreases slightly owing to the backscatter effect, and it is almost independent of the field size. The agreement between the surface dose estimated using the nanoDot and that measured using the ionization chamber for assessing X-ray beam qualities is less than 2%.ConclusionsThe nanoDot response is indicated as a function of HVL for the specified QIs, and it enables the direct surface dose measurement.     

Application of Monte Carlo calculations for validation of a treatment planning system in high dose rate brachytherapy

AimThe accuracy of treatment planning systems is of vital importance in treatment outcomes in brachytherapy. In the current study the accuracy of dose calculations of a high dose rate (HDR) brachytherapy treatment planning system (TPS) was validated using the Monte Carlo method.Materials and methodsThree ⁶⁰Co sources of the GZP6 afterloading brachytherapy system were modelled using MCNP4C Monte Carlo (MC) code. The dose distribution around all the sources was calculated by MC and a dedicated treatment planning system. The results of both methods were compared.ResultsThere was good agreement (<2%) between TPS and MC calculated dose distributions except at a point near the sources (<1 cm) and beyond the tip of the sources.ConclusionsOur study confirmed the accuracy of TPS calculated dose distributions for clinical use in HDR brachytherapy.     

FRoG dose computation meets Monte Carlo accuracy for proton therapy dose calculation in lung

PurposeTo benchmark and evaluate the clinical viability of novel analytical GPU-accelerated and CPU-based Monte Carlo (MC) dose-engines for spot-scanning intensity-modulated-proton-therapy (IMPT) towards the improvement of lung cancer treatment.MethodsNine patient cases were collected from the CNAO clinical experience and The Cancer Imaging Archive-4D-Lung-Database for in-silico study. All plans were optimized with 2 orthogonal beams in RayStation (RS) v.8. Forward calculations were performed with FRoG, an independent dose calculation system using a fast robust approach to the pencil beam algorithm (PBA), RS-MC (CPU for v.8) and general-purpose MC (gp-MC). Dosimetric benchmarks were acquired via irradiation of a lung-like phantom and ionization chambers for both a single-field-uniform-dose (SFUD) and IMPT plans. Dose-volume-histograms, dose-difference and γ-analyses were conducted.ResultsWith respect to reference gp-MC, the average dose to the GTV was 1.8% and 2.3% larger for FRoG and the RS-MC treatment planning system (TPS). FRoG and RS-MC showed a local γ-passing rate of ~96% and ~93%. Phantom measurements confirmed FRoG’s high accuracy with a deviation < 0.1%.ConclusionsDose calculation performance using the GPU-accelerated analytical PBA, MC-TPS and gp-MC code were well within clinical tolerances. FRoG predictions were in good agreement with both the full gp-MC and experimental data for proton beams optimized for thoracic dose calculations. GPU-accelerated dose-engines like FRoG may alleviate current issues related to deficiencies in current commercial analytical proton beam models. The novel approach to the PBA implemented in FRoG is suitable for either clinical TPS or as an auxiliary dose-engine to support clinical activity for lung patients.     

Determination of initial electron parameters by means of Monte Carlo simulations for the Siemens Artiste Linac 6 MV photon beam

AimIn this study, we investigated initial electron parameters of Siemens Artiste Linac with 6 MV photon beam using the Monte Carlo method.BackgroundIt is essential to define all the characteristics of initial electrons hitting the target, i.e. mean energy and full width of half maximum (FWHM) of the spatial distribution intensity, which is needed to run Monte Carlo simulations. The Monte Carlo is the most accurate method for simulation of radiotherapy treatments.Materials and methodsLinac head geometry was modeled using the BEAMnrc code. The phase space files were used as input file to DOSXYZnrc simulation to determine the dose distribution in a water phantom. We obtained percent depth dose curves and the lateral dose profile. All the results were obtained at 100 cm of SSD and for a 10 × 10 cm² field.ResultsWe concluded that there existed a good conformity between Monte Carlo simulation and measurement data when we used electron mean energy of 6.3 MeV and 0.30 cm FWHM value as initial parameters. We observed that FWHM values had very little effect on PDD and we found that the electron mean energy and FWHM values affected the lateral dose profile. However, these effects are between tolerance values.ConclusionsThe initial parameters especially depend on components of a linac head. The phase space file which was obtained from Monte Carlo Simulation for a linac can be used as calculation of scattering, MLC leakage, to compare dose distribution on patients and in various studies.     

A patient immobilization device for prone breast radiotherapy: Dosimetric effects and inclusion in the treatment planning system

PurposeTo assess the dosimetric impact of a patient positioning device for prone breast radiotherapy and assess the accuracy of a treatment planning system (TPS) in predicting this impact.MethodsBeam attenuation and build-up dose perturbations, quantified by ionization chamber and radiochromic film dosimetry, were evaluated for 3 components of the patient positioning device: the carbon fiber baseplate, the support cushions and the support wedge for the contralateral breast. Dose calculations were performed using the XVMC dose engine implemented in the Monaco TPS. All components were included during planning CT acquisition.ResultsBeam attenuation amounted to 7.57% (6 MV) and 5.33% (15 MV) for beams obliquely intersecting the couchtop–baseplate combination. Beams traversing large sections of the support wedge were attenuated by 12.28% (6 MV) and 9.37% (15 MV). For the support cushion foam, beam attenuation remained limited to 0.11% (6 MV) and 0.08% (15 MV) per centimeter thickness. A substantial loss of dose build-up was detected when irradiating through any of the investigated components. TPS dose calculations accurately predicted beam attenuation by the baseplate and support wedge. A manual density overwrite was needed to model attenuation by the support cushion foam. TPS dose calculations in build-up regions differed considerably from measurements for both open beams and beams traversing the device components.ConclusionsIrradiating through the components of the positioning device resulted in a considerable degradation of skin sparing. Inclusion of the device components in the treatment planning CT allowed to accurately model the most important attenuation effect, but failed to accurately predict build-up doses.     

Assessment of out-of-field doses in radiotherapy treatments of paediatric patients using Monte Carlo methods and measurements

PurposeTo assess out-of-field doses in radiotherapy treatments of paediatric patients, using Monte Carlo methods to implement a new model of the linear accelerator validated against measurements and developing a voxelized anthropomorphic paediatric phantom.MethodsCT images of a physical anthropomorphic paediatric phantom were acquired and a dosimetric planning using a TPS was obtained. The CT images were used to perform the voxelization of the physical phantom using the ImageJ software and later implemented in MCNP. In order to validate the Monte Carlo model, dose measurements of the 6 MV beam and Linac with 120 MLC were made in a clinical setting, using ionization chambers and a water phantom. Afterwards TLD measurements in the physical anthropomorphic phantom were performed in order to assess the out-of-field doses in the eyes, thyroid, c-spine, heart and lungs.ResultsThe Monte Carlo model was validated for in-field and out-of-field doses with average relative differences below 3%. The average relative differences between TLD measurements and Monte Carlo is 14,3% whilst the average relative differences between TLD and TPS is 55,8%. Moreover, organs up to 22.5 cm from PTV center show TLD and MCNP6 relative differences and TLD and TPS relative differences up to 21.2% and 92.0%, respectively.ConclusionsOur study provides a novel model that could be used in clinical research, namely in dose evaluation outside the treatment fields. This is particularly relevant, especially in pediatric patients, for studying new radiotherapy treatment techniques, since it can be used to estimate the development of secondary tumours.     

Evaluation of PENFAST – A fast Monte Carlo code for dose calculations in photon and electron radiotherapy treatment planning

The aim of the present study is to demonstrate the potential of accelerated dose calculations, using the fast Monte Carlo (MC) code referred to as PENFAST, rather than the conventional MC code PENELOPE, without losing accuracy in the computed dose. For this purpose, experimental measurements of dose distributions in homogeneous and inhomogeneous phantoms were compared with simulated results using both PENELOPE and PENFAST. The simulations and experiments were performed using a Saturne 43 linac operated at 12 MV (photons), and at 18 MeV (electrons). Pre-calculated phase space files (PSFs) were used as input data to both the PENELOPE and PENFAST dose simulations. Since depth–dose and dose profile comparisons between simulations and measurements in water were found to be in good agreement (within ±1% to 1 mm), the PSF calculation is considered to have been validated. In addition, measured dose distributions were compared to simulated results in a set of clinically relevant, inhomogeneous phantoms, consisting of lung and bone heterogeneities in a water tank. In general, the PENFAST results agree to within a 1% to 1 mm difference with those produced by PENELOPE, and to within a 2% to 2 mm difference with measured values. Our study thus provides a pre-clinical validation of the PENFAST code. It also demonstrates that PENFAST provides accurate results for both photon and electron beams, equivalent to those obtained with PENELOPE. CPU time comparisons between both MC codes show that PENFAST is generally about 9–21 times faster than PENELOPE.     

Development of an x-ray-opaque-marker system for quantitative phantom positioning in patient-specific quality assurance

PurposeWe developed an x-ray-opaque-marker (XOM) system with inserted fiducial markers for patient-specific quality assurance (QA) in CyberKnife (Accuray) and a general-purpose linear accelerator (linac). The XOM system can be easily inserted or removed from the existing patient-specific QA phantom. Our study aimed to assess the utility of the XOM system by evaluating the recognition accuracy of the phantom position error and estimating the dose perturbation around a marker.MethodsThe recognition accuracy of the phantom position error was evaluated by comparing the known error values of the phantom position with the values measured by matching the images with target locating system (TLS; Accuray) and on-board imager (OBI; Varian). The dose perturbation was evaluated for 6 and 10 MV single-photon beams through experimental measurements and Monte Carlo simulations.ResultsThe root mean squares (RMSs) of the residual position errors for the recognition accuracy evaluation in translations were 0.07 mm with TLS and 0.30 mm with OBI, and those in rotations were 0.13° with TLS and 0.15° with OBI. The dose perturbation was observed within 1.5 mm for 6 MV and 2.0 mm for 10 MV from the marker.ConclusionsSufficient recognition accuracy of the phantom position error was achieved using our system. It is unnecessary to consider the dose perturbation in actual patient-specific QA. We concluded that the XOM system can be utilized to ensure quantitative and accurate phantom positioning in patient-specific QA with CyberKnife and a general-purpose linac.     

A practical tool to evaluate dose distributions using radiochromic film in radiation oncology

PurposeTriple channel algorithm and specific procedures make more reliable radiochromic dosimetry for treatment planning verification and quality assurance in radiation therapy. A tool to obtain radiochromic dose distributions and compare them with the ones resulting from a treatment planning system was developed and applied.MethodsThe tool was developed as Microsoft Excel macro; it builds dose calibration curves against net optical density of Gafchromic EBT3 film, produces axial, coronal and sagittal dose maps and allows to evaluate them against dose distributions calculated by the Varian treatment planning system Eclipse using gamma index and gamma angle.ResultsThe net optical density standard errors of estimate of calibration curves at 6 MV Varian DBX600 linac energy were 0.2%, 0.4% and 0.2% for the red, green and blue channels. Tests of these curves by means of three independent eight dose points measurement series, at 15 MV and 6 MV Varian 2100C linac and at 6 MV DBX600 linac energies, showed less than 2% of dose errors for the red channel and less than 3% for the green channel in the range 100–450 cGy. The comparisons between dose distributions from Gafchromic EBT3 triple channel algorithm and the ones from Eclipse analytic anisotropic algorithm (AAA) showed values of gamma index 95th percentile between 0.6 and 1.0.ConclusionThe obtained results encourage the application of this tool in radiation therapy quality assurance.     

Dosimetric and Monte Carlo verification of jaws-only IMRT plans calculated by the Collapsed Cone Convolution algorithm for head and neck cancers

AimThe aim of this study is to verify the Prowess Panther jaws-only intensity modulated radiation therapy (JO-IMRT) treatment planning (TP) by comparing the TP dose distributions for head-and-neck (H&N) cancer with the ones simulated by Monte Carlo (MC).BackgroundTo date, dose distributions planned using JO-IMRT for H&N patients were found superior to the corresponding three-dimensional conformal radiotherapy (3D-CRT) plans. Dosimetry of the JO-IMRT plans were also experimentally verified using an ionization chamber, MapCHECK 2, and Octavius 4D and good agreements were shown.Materials and methodsDose distributions of 15 JO-IMRT plans of nasopharyngeal patients were recalculated using the EGSnrc Monte Carlo code. The clinical photon beams were simulated using the BEAMnrc. The absorbed dose to patients treated by fixed-field IMRT was computed using the DOSXYZnrc. The simulated dose distributions were then compared with the ones calculated by the Collapsed Cone Convolution (CCC) algorithm on the TPS, using the relative dose error comparison and the gamma index using global methods implemented in PTW-VeriSoft with 3%/3 mm, 2%/2 mm, 1%/1 mm criteria.ResultsThere is a good agreement between the MC and TPS dose. The average gamma passing rates were 93.3 ± 3.1%, 92.8 ± 3.2%, 92.4 ± 3.4% based on the 3%/3 mm, 2%/2 mm, 1%/1 mm criteria, respectively.ConclusionsAccording to the results, it is concluded that the CCC algorithm was adequate for most of the IMRT H&N cases where the target was not immediately adjacent to the critical structures.     

Experimental benchmarking of RayStation proton dose calculation algorithms inside and outside the target region in heterogeneous phantom geometries

PurposeThe aim of the presented study was to complement existing literature on benchmarking proton dose by comparing dose calculations with experimental measurements in heterogeneous phantom. Points of interest inside and outside the target were considered to quantify the magnitude of calculation uncertainties in current and previous proton therapy practice that might especially have an impact on the dose in organs at risk (OARs).MethodsThe RayStation treatment planning system (RaySearch Laboratories), offering two dose calculation algorithms for pencil beam scanning in proton therapy, i.e., Pencil Beam (PB) and Monte Carlo (MC), was utilized. Treatment plans for a target located behind the interface of the heterogeneous tissues were generated. Dose measurements within and behind the target were performed in a water phantom with embedded slabs of various tissue equivalent materials and 24 PinPoint ionization chambers (PTW). In total 12 test configurations encompassing two different target depths, oblique beam incidence of 30 degrees and range shifter, were considered.ResultsPB and MC calculated doses agreed equally well with the measurements for all test geometries within the target, including the range shifter (mean dose differences ± 3%). Outside the target, the maximum dose difference of 9% (19%) was observed for MC (PB) for the oblique beam incidence and inserted range shifter.ConclusionThe accuracy of MC dose algorithm was superior compared to the PB algorithm, especially outside the target volumes. MC based dose calculation should therefore be preferred in treatment scenarios with heterogeneities, especially to reduce clinically relevant uncertainties for OARs.     

Monte Carlo comparison of superficial dose between flattening filter free and flattened beams

This study investigates the superficial dose from FFF beams in comparison with the conventional flattened ones using a Monte Carlo (MC) method. Published phase-space files which incorporated real geometry of a TrueBeam accelerator were used for the dose calculation in phantom and clinical cases. The photon fluence on the central axis is 3 times that of a flattened beam for a 6 MV FFF beam and 5 times for a 10 MV beam. The mean energy across the field in air at the phantom surface is 0.92–0.95 MeV for the 6 MV FFF beam and 1.18–1.30 MeV for the corresponding flattened beam. At 10 MV, the values are 1.52–1.72 and 2.15–2.87 MeV for the FFF and flattened beams, respectively. The phantom dose at the depth of 1 mm in the 6 MV FFF beam is 6% ± 2.5% (of the maximum dose) higher compared to the flattened beam for a 25 × 25 cm² field and 14.6% ± 1.9% for the 2 × 2 cm² field. For the 10 MV beam, the corresponding differences are 3.4% ± 1.5% and 10.7% ± 0.6%. The skin dose difference at selected points on the patient's surface between the plans using FFF and flattened beams in the head-and-neck case was 6.5% ± 2.3% (1SD), and for the breast case it was 6.4% ± 2.3%. The Monte Carlo simulations showed that due to the lower mean energy in the FFF beam, the clinical superficial dose is higher without the flattening filter compared to the flattened beam.     

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.