Validation of XiO Electron Monte Carlo-based calculations by measurements in a homogeneous phantom and by EGSnrc calculations in a heterogeneous phantom

The purpose of the present study is to perform a clinical validation of a new commercial Monte Carlo (MC) based treatment planning system (TPS) for electron beams, i.e. the XiO 4.60 electron MC (XiO eMC). Firstly, MC models for electron beams (4, 8, 12 and 18 MeV) have been simulated using BEAMnrc user code and validated by measurements in a homogeneous water phantom. Secondly, these BEAMnrc models have been set as the reference tool to evaluate the ability of XiO eMC to reproduce dose perturbations in the heterogeneous phantom. In the homogeneous phantom calculations, differences between MC computations (BEAMnrc, XiO eMC) and measurements are less than 2% in the homogeneous dose regions and less than 1 mm shifting in the high dose gradient regions. As for the heterogeneous phantom, the accuracy of XiO eMC has been benchmarked with predicted BEAMnrc models. In the lung tissue, the overall agreement between the two schemes lies under 2.5% for the most tested dose distributions at 8, 12 and 18 MeV and is better than the 4 MeV one. In the non-lung tissue, a good agreement has been found between BEAMnrc simulation and XiO eMC computation for 8, 12 and 18 MeV. Results are worse in the case of 4 MeV calculations (discrepancies ≈ 4%). XiO eMC can predict dose perturbation induced by high-density heterogeneities for 8, 12 and 18 MeV. However, significant deviations found in the case of 4 MeV demonstrate that caution is necessary in using XiO eMC at lower electron energies.     

Full Monte Carlo and measurement-based overall performance assessment of improved clinical implementation of eMC algorithm with emphasis on lower energy range

New version 13.6.23 of the electron Monte Carlo (eMC) algorithm in Varian Eclipse™ treatment planning system has a model for 4 MeV electron beam and some general improvements for dose calculation. This study provides the first overall accuracy assessment of this algorithm against full Monte Carlo (MC) simulations for electron beams from 4 MeV to 16 MeV with most emphasis on the lower energy range. Beams in a homogeneous water phantom and clinical treatment plans were investigated including measurements in the water phantom. Two different material sets were used with full MC: (1) the one applied in the eMC algorithm and (2) the one included in the Eclipse™ for other algorithms. The results of clinical treatment plans were also compared to those of the older eMC version 11.0.31. In the water phantom the dose differences against the full MC were mostly less than 3% with distance-to-agreement (DTA) values within 2 mm. Larger discrepancies were obtained in build-up regions, at depths near the maximum electron ranges and with small apertures. For the clinical treatment plans the overall dose differences were mostly within 3% or 2 mm with the first material set. Larger differences were observed for a large 4 MeV beam entering curved patient surface with extended SSD and also in regions of large dose gradients. Still the DTA values were within 3 mm. The discrepancies between the eMC and the full MC were generally larger for the second material set. The version 11.0.31 performed always inferiorly, when compared to the 13.6.23.     

Validation of a commercial TPS based on the VMC++ Monte Carlo code for electron beams: Commissioning and dosimetric comparison with EGSnrc in homogeneous and heterogeneous phantoms

The aim of the present work was the validation of the VMC⁺⁺ Monte Carlo (MC) engine implemented in the Oncentra Masterplan (OMTPS) and used to calculate the dose distribution produced by the electron beams (energy 5-12 MeV) generated by the linear accelerator (linac) Primus (Siemens), shaped by a digital variable applicator (DEVA). The BEAMnrc/DOSXYZnrc (EGSnrc package) MC model of the linac head was used as a benchmark.Commissioning results for both MC codes were evaluated by means of 1D Gamma Analysis (2%, 2 mm), calculated with a home-made Matlab (The MathWorks) program, comparing the calculations with the measured profiles. The results of the commissioning of OMTPS were good [average gamma index (γ) > 97%]; some mismatches were found with large beams (size ≥ 15 cm). The optimization of the BEAMnrc model required to increase the beam exit window to match the calculated and measured profiles (final average γ > 98%).Then OMTPS dose distribution maps were compared with DOSXYZnrc with a 2D Gamma Analysis (3%, 3 mm), in 3 virtual water phantoms: (a) with an air step, (b) with an air insert, and (c) with a bone insert.The OMTPD and EGSnrc dose distributions with the air-water step phantom were in very high agreement (γ ∼ 99%), while for heterogeneous phantoms there were differences of about 9% in the air insert and of about 10–15% in the bone region. This is due to the Masterplan implementation of VMC⁺⁺ which reports the dose as “dose to water”, instead of “dose to medium”.     

Characterization of irregular electron beam for boost dose after whole breast irradiation

Irradiating a tumor bed with boost dose after whole breast irradiation helps reducing the probability of local recurrence. However, the success of electron beam treatment with a small area aiming to cover a superficial lesion is a dual challenge as it requires an adequate dosimetry beside a double check for dose coverage with an estimation of various combined uncertainty of tumor location and losing lateral electron equilibrium within small field dimensions.Aim of workthis work aims to measure the electron beam fluence within different field dimensions and the deviation from measurement performed in standard square electron applicator beam flatness and symmetry, then to calculate the average range of the correction factor required to overcome the loss of lateral electron equilibrium.Material and methodthe electron beam used in this work generated from the linear accelerator model ELEKTA Precise and dosimetry system used were a pair of PTW Pin Point ion chambers for electron beam dosimetry at standard conditions and assessment of beam quality at a reference depth of measurement, with an automatic water phantom, then a Roos ion chamber was used for absolute dose measurement, and PTW 2Darray to investigate the beam fluence of four applicators 6, 10, 14 and 20 cm² and 4 rectangular cutouts 6 × 14, 8 × 14, 6 × 17 and 8 × 17 cm², the second part was clinical application which was performed in a precise treatment planning system and examined boost dose after whole breast irradiation.Resultsrevealed that lower energy (6MeV and 8MeV) showed the loss of lateral electron equilibrium and deviation from measurements of a standard applicator more than the high energy (15 MeV) which indicated that the treatment of superficial dose with 6MeV required higher monitor unit to allow for the loss of lateral electron equilibrium and higher margin as well.     

Towards breast cancer rotational radiotherapy with synchrotron radiation

PurposeWe performed the first investigations, via measurements and Monte Carlo simulations on phantoms, of the feasibility of a new technique for synchrotron radiation rotational radiotherapy for breast cancer (SR³T).MethodsA Monte Carlo (MC) code based on Geant4 toolkit was developed in order to simulate the irradiation with the SR³T technique and to evaluate the skin sparing effect in terms of centre-to-periphery dose ratio at different energies in the range 60–175 keV. Preliminary measurements were performed at the Australian Synchrotron facility. Radial dose profiles in a 14-cm diameter polyethylene phantom were measured with a 100-mm pencil ionization chamber for different beam sizes and compared with the results of MC simulations. Finally, the dose painting feasibility was demonstrated with measurements with EBT3 radiochromic films in a phantom and collimating the SR beam at 1.5 cm in the horizontal direction.ResultsMC simulations showed that the SR³T technique assures a tumour-to-skin absorbed dose ratio from about 7:1 (at 60 keV photon energy) to about 10:1 (at 175 keV), sufficient for skin sparing during radiotherapy. The comparison between the results of MC simulations and measurements showed an agreement within 5%. Two off-centre foci were irradiated shifting the rotation centre in the horizontal direction.ConclusionsThe SR³T technique permits to obtain different dose distributions in the target with multiple rotations and can be guided via synchrotron radiation breast computed tomography imaging, in propagation based phase-contrast conditions. Use of contrast agents like iodinated solutions or gold nanoparticles for dose enhancement (DE-SR³T) is foreseen and will be investigated in future work.     

Dosimetric shield evaluation with tungsten sheet in 4, 6, and 9 MeV electron beams

In electron radiotherapy, shielding material is required to attenuate beam and scatter. A newly introduced shielding material, tungsten functional paper (TFP), has been anticipated to become a very useful device that is lead-free, light, flexible, and easily processed, containing very fine tungsten powder at as much as 80% by weight. The purpose of this study was to investigate the dosimetric changes due to TFP shielding for electron beams. TFP (thickness 0–15 mm) was placed on water or a water-equivalent phantom. Percentage depth ionization and transmission were measured for 4, 6, and 9 MeV electron beams. Off-center ratio was also measured using film dosimetry at depth of dose maximum under similar conditions. Then, beam profiles and transmission with two shielding materials, TFP and lead, were evaluated. Reductions of 95% by using TFP at 0.5 cm depth occurred at 4, 9, and 15 mm with 4, 6, and 9 MeV electron beams, respectively. It is found that the dose tend to increase at the field edge shaped with TFP, which might be influenced by the thickness. TFP has several unique features and is very promising as a useful tool for radiation protection for electron beams, among others.     

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.     

Measurement of the influence of titanium hip prosthesis on therapeutic electron beam dose distributions in a novel pelvic phantom

PurposeTo investigate the degree of 18 and 22 MeV electron beam dose perturbations caused by unilateral hip titanium (Ti) prosthesis.MethodsMeasurements were acquired using Gafchromic EBT2 film in a novel pelvic phantom made out of Nylon-12 slices in which a Ti-prosthesis is embedded. Dose perturbations were measured and compared using depth doses for 8 × 8, 10 × 10 and 11 × 11 cm² applicator-defined field sizes at 95 cm source-surface-distance (SSD). Comparisons were also made between film data at 100 cm SSD for a 10 × 10 cm² field and dose calculations made on CMS XiO treatment planning system utilizing the pencil beam algorithm. The extent of dose deviations caused by the Ti prosthesis based on film data was quantified through the dose enhancement factor (DEF), defined as the ratio of the dose influenced by the prosthesis and the unchanged beam.ResultsAt the interface between Nylon-12 and the Ti implant on the prosthesis entrance side, the dose increased to values of 21 ± 1% and 23 ± 1% for 18 and 22 MeV electron beams, respectively. DEFs increased with increasing electron energy and field size, and were found to fall off quickly with distance from the nylon-prosthesis interface. A comparison of film and XiO depth dose data for 18 and 22 MeV gave relative errors of 20% and 25%, respectively.ConclusionThis study outlines the lack of accuracy of the XiO TPS for electron planning in highly heterogeneous media. So a dosimetric error of 20–25% could influence clinical outcome.     

Dosimetric characteristics of a MOSFET dosimeter for clinical electron beams

The fundamental dosimetric characteristics of commercially available metal oxide semiconductor field effect transistor (MOSFET) detectors were studied for clinical electron beam irradiations. MOSFET showed excellent linearity against doses measured using an ion chamber in the dose range of 20–630 cGy. MOSFET reproducibility is better at high doses compared to low doses. The output factors measured with the MOSFET were within ±3% when compared with those measured with a parallel plate chamber. From 4 to 12 MeV, MOSFETs showed a large angular dependence in the tilt directions and less in the axial directions. MOSFETs do not show any dose-rate dependence between 100 and 600 MU/min. However, MOSFETs have shown under-response when the dose per pulse of the beam is decreased. No measurable effect in MOSFET response was observed in the temperature range of 23–40 °C. The energy dependence of a MOSFET dosimeter was within ±3.0% for 6–18 MeV electron beams and 5.5% for 4 MeV ones. This study shows that MOSFET detectors are suitable for dosimetry of electron beams in the energy range of 4–18 MeV.     

Characterization of prompt gamma-ray emission with respect to the Bragg peak for proton beam range verification: A Monte Carlo study

In this paper we report a Geant4 simulation study to investigate the characteristic prompt gamma (PG) emission in a water phantom for real-time monitoring of the Bragg peak (BP) during proton beam irradiation. The PG production, emission spatial correlation with the BP, and position preference for detection with respect to the BP have been quantified in different PG energy windows as a function of proton pencil-beam energy from 100 to 200 MeV. The PG response to small BP shifts was evaluated using a 2 cm-thick slab with different human body materials embedded in a water phantom. Our results show that the prominent characteristic PG emissions of 4.44, 5.21 and 6.13 MeV exhibit distinctive correlation with the dose deposition curve. The accuracy in BP position identification using these characteristic PG rays is highly consistent as the beam energy increases from 100 to 200 MeV. There exists a position preference for PG detection with respect to the BP position, which has a strong dependence on the proton beam energy and PG energies. It was also observed that a submillimeter shift of the BP position can be realized by using PG signals. These results indicate that the characteristic PG signal is sensitive and reliable for BP tracking. Although the maximization of the PG measurement associated with the BP is difficult, it can be optimized with energy and detection position preferences.     

A depth dose study between AAA and AXB algorithm against Monte Carlo simulation using AIP CT of a 4D dataset from a moving phantom

Aim

To identifying depth dose differences between the two versions of the algorithms using AIP CT of a 4D dataset.

Dosimetric comparison between two MLC systems commonly used for stereotactic radiosurgery and radiotherapy: A Monte Carlo and experimental study

In this work dosimetric parameters of two multi-leaf collimator (MLC) systems, namely the beam modulator (BM), which is the MLC commercial name for Elekta “Synergy S” linear accelerator and Radionics micro-MLC (MMLC), are compared using measurements and Monte Carlo simulations. Dosimetric parameters, such as percentage depth doses (PDDs), in-plane and cross-plane dose profiles, and penumbras for different depths and field sizes of the 6 MV photon beams were measured using ionization chamber and a water tank. The collimator leakages were measured using radiographic films. MMLC and BM were modeled using the EGSnrc-based BEAMnrc Monte Carlo code and above dosimetric parameters were calculated. The energy fluence spectra for the two MLCs were also determined using the BEAMnrc and BEAMDP. Dosimetric parameters of the two MLCs were similar, except for penumbras. Leaf-side and leaf-end 80–20% dose penumbras at 10 cm depth for a 10 × 10 cm² field size were 4.8 and 5.1 mm for MMLC and 5.3 mm and 6.3 mm for BM, respectively. Both Radionics MMLC and Elekta BM can be used effectively based on their dosimetric characteristics for stereotactic radiosurgery and radiotherapy, although the former showed slightly sharper dose penumbra especially in the leaf-end direction.     

The determination of a dose deposited in reference medium due to (p,n) reaction occurring during proton therapy

AimThe aim of the investigation was to determine the undesirable dose coming from neutrons produced in reactions (p,n) in irradiated tissues represented by water.BackgroundProduction of neutrons in the system of beam collimators and in irradiated tissues is the undesirable phenomenon related to the application of protons in radiotherapy. It makes that proton beams are contaminated by neutrons and patients receive the undesirable neutron dose.Materials and methodsThe investigation was based on the Monte Carlo simulations (GEANT4 code). The calculations were performed for five energies of protons: 50 MeV, 55 MeV, 60 MeV, 65 MeV and 75 MeV. The neutron doses were calculated on the basis of the neutron fluence and neutron energy spectra derived from simulations and by means of the neutron fluence–dose conversion coefficients taken from the ICRP dosimetry protocol no. 74 for the antero-posterior irradiation geometry.ResultsThe obtained neutron doses are much less than the proton ones. They do not exceed 0.1%, 0.4%, 0.5%, 0.6% and 0.7% of the total dose at a given depth for the primary protons with energy of 50 MeV, 55 MeV, 60 MeV, 65 MeV and 70 MeV, respectively.ConclusionsThe neutron production takes place mainly along the central axis of the beam. The maximum neutron dose appears at about a half of the depth of the maximum proton dose (Bragg peak), i.e. in the volume of a healthy tissue. The doses of neutrons produced in the irradiated medium (water) are about two orders of magnitude less than the proton doses for the considered range of energy of protons.     

Radiochromic films for dental CT dosimetry: A feasibility study

Dental CT dose evaluations are commonly performed using thermoluminescent dosimeters (TLD) inside anthropomorphic phantoms. Radiochromic films with good sensitivity in the X-ray diagnostic field have recently been developed and are commercially available as GAFCHROMIC XR-QA. There are potential advantages in the use of radiochromic films such as a more comprehensive dosimetry thanks to the adjustable size of the film samples. The purpose of this study was to investigate the feasibility of using radiochromic films for dental CT dose evaluations.Film samples were cut with a width of 5 mm and a length of 25 mm (strips), the same size as the Alderson Rando anthropomorphic phantom holes used in this study. Dental CT dose measurements were performed using simultaneously both TLD and radiochromic strips in the same phantom sites. Two equipment types were considered for dental CT examinations: a 16 slice CT and a cone beam CT. Organ equivalent doses were then obtained averaging the measurements from the sites of the same organ and effective doses were calculated using ICRP 103 weighting factors. The entire procedure was repeated four times for each CT in order to compare also the repeatability of the two dosimeter types.A linear correlation was found between the absorbed dose evaluated with radiochromic films and with TLD, with slopes of 0.930 and 0.944 (correlation r > 0.99). The maximum difference between the two dosimeter’s measurements was 25%, whereas the average difference was 7%. The measurement repeatability was comparable for the two dosimeters at cumulative doses above 15 mGy (estimated uncertainty at 1 sigma level of about 5%), whereas below this threshold radiochromic films show a greater dispersion of data, of about 10% at 1 sigma level. We obtained, using respectively Gafchromic and TLD measurements, effective dose values of 107 μSv and 117 μSv (i.e. difference of 8.6%) for the cone beam CT and of 523 μSv and 562 μSv (i.e. difference of 7%) for the multislice CT.This study demonstrates the feasibility of radiochromic films for dental CT dosimetry, pointing out a good agreement with the results obtained using TLD, with potential advantages and the chance of a more extensive dose investigation.     

Commissioning Monte Carlo algorithm for robotic radiosurgery using cylindrical 3D-array with variable density inserts

IntroductionTo commission the Monte Carlo (MC) algorithm based model of CyberKnife robotic stereotactic system (CK) and evaluate the feasibility of patient specific QA using the ArcCHECK cylindrical 3D-array (AC) with Multiplug inserts (MP).ResultsFour configurations were used for simple beam setup and two for patient QA, replacing water equivalent inserts by lung. For twelve collimators (5–60 mm) in simple setup, mean (SD) differences between MC and RayTracing algorithm (RT) of the number of points failing the 3%/1 mm gamma criteria were 1(1), 1(3), 1(2) and 1(2) for the four MP configurations. Tracking fiducials were placed within AC for patient QA. Single lung insert setup resulted in mean gamma-index 2%/2 mm of 90.5% (range [74.3–95.9]) and 82.3% ([66.8–94.5]) for MC and RT respectively, while 93.5% ([86.8–98.2]) and 86.2% ([68.7–95.4]) in presence of largest inhomogeneities, showing significant differences (p < 0.05).DiscussionAfter evaluating the potential effects, 1.12 g/cc PMMA and 0.09 g/cc lung material assignment showed the best results. Overall, MC-based model showed superior results compared to RT for simple and patient specific testing, using a 2%/2 mm criteria. Results are comparable with other reported commissionings for flattening filter free (FFF) delivery. Further improvement of MC calculation might be challenging as Multiplan has limited material library.ConclusionsThe AC with Multiplug allowed for comprehensive commissioning of CyberKnife MC algorithm and is useful for patient specific QA for stereotactic body radiation therapy. MC calculation accuracy might be limited due to Multiplan’s insufficient material library; still results are comparable with other reported commissioning measurements using FFF beams.     

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.     

Surface dose measurements with GafChromic EBT film for 6 and 18 MV photon beams

The aim of this study was to determine the surface doses using GafChromic EBT films and compare them with plane-parallel ionization chamber measurements for 6 and 18 MV high energy photon beams. The measurements were made in a water equivalent solid phantom in the build-up region of the 6 and 18 MV photon beams at 100 cm SSD for various field sizes. Markus type plane-parallel ion chamber with fixed-separation between collecting electrodes was used to measure the percent depth doses. GafChromic EBT film measurements were performed both on the phantom surface and maximum dose depth at the same geometry with ion chamber measurements. The surface doses found using GafChromic EBT film were 15%, 20%, 29%and 39% ± 2% (1SD) for 6 MV photons, 6%, 11%, 23% and 32% ± 2% (1SD) for 18 MV photons at 5, 10, 20 and 30 cm² field sizes, respectively. GafChromic EBT film provides precise measurements for surface dose in the high energy photons. Agreement between film and plane-parallel chamber measurements was found to be within ±3% for 18 MV photon beams. There was 5% overestimate on the surface doses when compared with the plane-parallel chamber measurements for all field sizes in the 6 MV photon beams.     

Determination of inflection points of CyberKnife dose profiles within acceptability criteria of deviations in measurements

AimThe aim of this study was to determine the Inflection Points (IPs) of flattening filter free (FFF) CyberKnife dose profiles for cone-based streotactic radiotherapy. In addition, dosimetric field sizes were determined.BackgroundThe increased need for treatment in the early stages of cancer necessitated the treatment of smaller tumors. However, efforts in that direction required the modeling accuracy of the beam. Removal of the flattening filter (FF) from the path of x-ray beam has provided the solution to those efforts, but required a different normalization approach for the beam to ensure the delivery of the dose accurately. As a solution, researchers proposed a normalization factor based on IPs.Materials and methodsMeasurements using microDiamond (PTW 60019), Diode SRS (PTW 60018) and Monte Carlo (MC) calculations of dose profiles were completed at SAD 80 cm and 5 cm depth for 15–60 mm cones. Performance analysis of detectors with respect to MC calculation was carried out. Gamma evaluation method was used to determine achievable acceptability criteria for FFF CyberKnife beams.ResultsAcceptability within (3%–0.5 mm) was found to be anachievable criterion for all dose profile measurements of the cone beams used in this study. To determine the IP, the first and second derivatives of the dose profile were determined via the cubic spline interpolation technique.ConclusionDerivatives of the interpolated profiles showed that locations of IPs and 50% isodose points coincide.     

Verification in the water phantom of the irradiation time calculation done by the algorithm used in intraoperative radiotherapy

AimThe investigation of the irradiation time calculation accuracy of the GGPB algorithm used for IORT.BackgroundConventionally, breast conserving therapy consists of breast conserving surgery followed by postoperative whole breast irradiation and boost. The use of intraoperative radiotherapy (IORT) enables the boost to be delivered already during the surgery. In this case, the treatment dose for IORT can be calculated by use of General Gaussian Pencil Beam (GGPB) algorithm, which is implemented in TPS Eclipse.Materials and methodsPDDs and OFs for electron beams from Mobetron and all available applicators were measured in order to configure the GGPB algorithm. Afterwards, the irradiation times for the prescribed dose of 3 Gy were calculated by means of it. The results of calculations were verified in the water phantom using the Marcus ionization chamber.ResultsThe results differed between energies. For 6 MeV the irradiation times calculated by the GGPB algorithm were correct, for the energy of 9 MeV they were too small and for the energy of 4 MeV they were too large for applicators with smaller diameters, while acceptable for the remaining ones.ConclusionThe GGPB algorithm can be used in intraoperative radiotherapy for energy and applicator sets for which no significant difference between the measured and the prescribed dose was obtained. For the rest of energy-applicator sets the configuration should be verified and possibly repeated.