A DICOM-RT-based toolbox for the evaluation and verification of radiotherapy plans

PENELOPE simulations and experiment for 6 MV clinac iX accelerator for standard and small static fields

The goal of this work was to produce accurate data for use as a 'gold standard' and a valid tool for measurements in reference dosimetry for standard/small static field sizes from 0.5 × 0.5 to 10 × 10 cm². It is based on the accuracy of the phase space files (PSFs) as a key quantity. Because the IAEA general public database provides few PSFs for the Varian iX, we simulated the head through Monte Carlo (MC) simulations and calculated validated PSFs for 12 square field sizes including seven for small static fields. The resulting dosimetric calculations allowed us to reach a good level of agreement in comparison to our relative and absolute dose measurements performed on a Varian iX in water phantom. Measured and MC calculated output factors were investigated for different detectors. Based on the TRS 483 formalism and MC (PENELOPE/penEasy), we calculated output correction factors for the unshielded Diode-E (T60017) and the PinPoint-3D (T31016) micro-chamber according to manufacturers' blueprints. Our MC results were in agreement with the recommended data; they compete with recent measurements and MC simulations and in particular the TRS 483 MC data obtained from similar simulations. Moreover, our MC results provide supplemental data in comparison to TRS 483 data in particular for the PinPoint-3D (T31016). We suggest our MC output correction factors as new datasets for future TRS compilations. The work was substantial, used different robust MC strategies depending on the scoring regions, and led in most cases to uncertainties of less than 1%.

Feasibility of a GATE Monte Carlo platform in a clinical pretreatment QA system for VMAT treatment plans using TrueBeam with an HD120 multileaf collimator

Purpose

To evaluate the quality of patient‐specific complicated treatment plans, including commercialized treatment planning systems (TPS) and commissioned beam data, we developed a process of quality assurance (QA) using a Monte Carlo (MC) platform. Specifically, we constructed an interface system that automatically converts treatment plan and dose matrix data in digital imaging and communications in medicine to an MC dose‐calculation engine. The clinical feasibility of the system was evaluated.

Materials and Methods

A dose‐calculation engine based on GATE v8.1 was embedded in our QA system and in a parallel computing system to significantly reduce the computation time. The QA system automatically converts parameters in volumetric‐modulated arc therapy (VMAT) plans to files for dose calculation using GATE. The system then calculates dose maps. Energies of 6 MV, 10 MV, 6 MV flattening filter free (FFF), and 10 MV FFF from a TrueBeam with HD120 were modeled and commissioned. To evaluate the beam models, percentage depth dose (PDD) values, MC calculation profiles, and measured beam data were compared at various depths (D_max, 5 cm, 10 cm, and 20 cm), field sizes, and energies. To evaluate the feasibility of the QA system for clinical use, doses measured for clinical VMAT plans using films were compared to dose maps calculated using our MC‐based QA system.

Results

A LINAC QA system was analyzed by PDD and profile according to the secondary collimator and multileaf collimator (MLC). Values for MC calculations and TPS beam data obtained using CC13 ion chamber (IBA Dosimetry, Germany) were consistent within 1.0%. Clinical validation using a gamma index was performed for VMAT treatment plans using a solid water phantom and arbitrary patient data. The gamma evaluation results (with criteria of 3%/3 mm) were 98.1%, 99.1%, 99.2%, and 97.1% for energies of 6 MV, 10 MV, 6 MV FFF, and 10 MV FFF, respectively.

Conclusions

We constructed an MC‐based QA system for evaluating patient treatment plans and evaluated its feasibility in clinical practice. We observed robust agreement between dose calculations from our QA system and measurements for VMAT plans. Our QA system could be useful in other clinical settings, such as small‐field SRS procedures or analyses of secondary cancer risk, for which dose calculations using TPS are difficult to verify.

Investigating the effect of dental implant materials with different densities on radiotherapy dose distribution using Monte-Carlo simulation and pencil beam convolution algorithm

OBJECTIVES

The aim of this study was to investigate the effect of dental implant materials with different physical densities on dose distribution for head and neck cancer radiotherapy planning.

METHODS

Titanium (Ti), Titanium alloy (Ti-6Al-4V), Zirconia (Y-TZP), Zirconium oxide (ZrO₂), Alumina (Al₂O₃) and polyetheretherketone (PEEK) dental implant materials were used for determination of implant material effect on dose distribution. Dental implant effect was investigated by using pencil beam convolution (PBC) algorithm of Eclipse treatment planning systems (TPS) and Monte Carlo (MC) simulation technique. 6 MV photon beam of the Varian 2300 C/D linear accelerator was simulated by EGSnrc-based BEAMnrc MC code system.

RESULTS

Reasonable consistency was determined for percentage depth dose (PDD) curves between MC simulation and water phantom measurements at 6.4 MeV initial electron energy. The consistency between modelled linear accelerator PDD curve calculations and water-phantom PDD measurements were compatible within 1 % range. The dose increase in front of the dental implant calculated by MC simulation is in the range of 0.4-20.2%. We found by MC and PBC calculations that the differences in dose increase in front of the dental implant materials is in the range of 0.1-17.2% and is dependent on the physical density of the dental implant.

CONCLUSIONS

Dose increase for Zirconia was noted to be maximum while PEEK implant dose increase was minimum among the whole dental implant materials studied. This study revealed that the Eclipse TPS PBC algorithm could not accurately estimate the backscatter radiation from dental implant materials.

A patch-based pseudo-CT approach for MRI-only radiotherapy in the pelvis

PURPOSE

In radiotherapy based only on magnetic resonance imaging (MRI), knowledge about tissue electron densities must be derived from the MRI. This can be achieved by converting the MRI scan to the so-called pseudo-computed tomography (pCT). An obstacle is that the voxel intensities in conventional MRI scans are not uniquely related to electron density. The authors previously demonstrated that a patch-based method could produce accurate pCTs of the brain using conventional T1-weighted MRI scans. The method was driven mainly by local patch similarities and relied on simple affine registrations between an atlas database of the co-registered MRI/CT scan pairs and the MRI scan to be converted. In this study, the authors investigate the applicability of the patch-based approach in the pelvis. This region is challenging for a method based on local similarities due to the greater inter-patient variation. The authors benchmark the method against a baseline pCT strategy where all voxels inside the body contour are assigned a water-equivalent bulk density. Furthermore, the authors implement a parallelized approximate patch search strategy to speed up the pCT generation time to a more clinically relevant level.

METHODS

The data consisted of CT and T1-weighted MRI scans of 10 prostate patients. pCTs were generated using an approximate patch search algorithm in a leave-one-out fashion and compared with the CT using frequently described metrics such as the voxel-wise mean absolute error (MAEvox) and the deviation in water-equivalent path lengths. Furthermore, the dosimetric accuracy was tested for a volumetric modulated arc therapy plan using dose-volume histogram (DVH) point deviations and γ-index analysis.

RESULTS

The patch-based approach had an average MAEvox of 54 HU; median deviations of less than 0.4% in relevant DVH points and a γ-index pass rate of 0.97 using a 1%/1 mm criterion. The patch-based approach showed a significantly better performance than the baseline water pCT in almost all metrics. The approximate patch search strategy was 70x faster than a brute-force search, with an average prediction time of 20.8 min.

CONCLUSIONS

The authors showed that a patch-based method based on affine registrations and T1-weighted MRI could generate accurate pCTs of the pelvis. The main source of differences between pCT and CT was positional changes of air pockets and body outline.

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 4MeV 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 4MeV to 16MeV 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 2mm. 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 2mm with the first material set. Larger differences were observed for a large 4MeV beam entering curved patient surface with extended SSD and also in regions of large dose gradients. Still the DTA values were within 3mm. 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.

A design of a DICOM-RT-based tool box for nonrigid 4D dose calculation

The study was aimed to introduce a design of a DICOM-RT-based tool box to facilitate 4D dose calculation based on deformable voxel-dose registration. The computational structure and the calculation algorithm of the tool box were explicitly discussed in the study. The tool box was written in MATLAB in conjunction with CERR. It consists of five main functions which allow a) importation of DICOM-RT-based 3D dose plan, b) deformable image registration, c) tracking voxel doses along breathing cycle, d) presentation of temporal dose distribution at different time phase, and e) derivation of 4D dose. The efficacy of using the tool box for clinical application had been verified with nine clinical cases on retrospective-study basis. The logistic and the robustness of the tool box were tested with 27 applications and the results were shown successful with no computational errors encountered. In the study, the accumulated dose coverage as a function of planning CT taken at end-inhale, end-exhale, and mean tumor position were assessed. The results indicated that the majority of the cases (67%) achieved maximum target coverage, while the planning CT was taken at the temporal mean tumor position and 56% at the end-exhale position. The comparable results to the literature imply that the studied tool box can be reliable for 4D dose calculation. The authors suggest that, with proper application, 4D dose calculation using deformable registration can provide better dose evaluation for treatment with moving target.

Quantification of dose differences between two versions of Acuros XB algorithm compared to Monte Carlo simulations--the effect on clinical patient treatment planning

A commercialized implementation of linear Boltzmann transport equation solver, the Acuros XB algorithm (AXB), represents a class of most advanced type 'c' photon radiotherapy dose calculation algorithms. The purpose of the study was to quantify the effects of the modifications implemented in the more recent version 11 of the AXB (AXB11) compared to the first commercial implementation, version 10 of the AXB (AXB10), in various anatomical regions in clinical treatment planning. Both versions of the AXB were part of Varian's Eclipse clinical treatment planning system and treatment plans for 10 patients were created using intensity-modulated radiotherapy (IMRT) and volumetric-modulated arc radiotherapy (VMAT). The plans were first created with the AXB10 and then recalculated with the AXB11 and full Monte Carlo (MC) simulations. Considering the full MC simulations as reference, a DVH analysis for gross tumor and planning target volumes (GTV and PTV) and organs at risk was performed, and also 3D gamma agreement index (GAI) values within a 15% isodose region and for the PTV were determined. Although differences up to 12% in DVH analysis were seen between the MC simulations and the AXB, based on the results of this study no general conclusion can be drawn that the modifications made in the AXB11 compared to the AXB10 would imply that the dose calculation accuracy of the AXB10 would be inferior to the AXB11 in the clinical patient treatment planning. The only clear improvement with the AXB11 over the AXB10 is the dose calculation accuracy in air cavities. In general, no large deviations are present in the DVH analysis results between the two versions of the algorithm, and the results of 3D gamma analysis do not favor one or the other. Thus it may be concluded that the results of the comprehensive studies assessing the accuracy of the AXB10 may be extended to the AXB11.

The accuracy of Acuros XB algorithm for radiation beams traversing a metallic hip implant - comparison with measurements and Monte Carlo calculations

In this study, the clinical benefit of the improved accuracy of the Acuros XB (AXB) algorithm, implemented in a commercial radiotherapy treatment planning system (TPS), Varian Eclipse, was demonstrated with beams traversing a high‐Z material. This is also the first study assessing the accuracy of the AXB algorithm applying volumetric modulated arc therapy (VMAT) technique compared to full Monte Carlo (MC) simulations. In the first phase the AXB algorithm was benchmarked against point dosimetry, film dosimetry, and full MC calculation in a water‐filled anthropometric phantom with a unilateral hip implant. Also the validity of the full MC calculation used as reference method was demonstrated. The dose calculations were performed both in original computed tomography (CT) dataset, which included artifacts, and in corrected CT dataset, where constant Hounsfield unit (HU) value assignment for all the materials was made. In the second phase, a clinical treatment plan was prepared for a prostate cancer patient with a unilateral hip implant. The plan applied a hybrid VMAT technique that included partial arcs that avoided passing through the implant and static beams traversing the implant. Ultimately, the AXB‐calculated dose distribution was compared to the recalculation by the full MC simulation to assess the accuracy of the AXB algorithm in clinical setting. A recalculation with the anisotropic analytical algorithm (AAA) was also performed to quantify the benefit of the improved dose calculation accuracy of type ‘c’ algorithm (AXB) over type ‘b’ algorithm (AAA). The agreement between the AXB algorithm and the full MC model was very good inside and in the vicinity of the implant and elsewhere, which verifies the accuracy of the AXB algorithm for patient plans with beams traversing through high‐Z material, whereas the AAA produced larger discrepancies.

PACS numbers: 87.55.‐x, 87.55.D‐, 87.55.K‐, 87.55.kd, 87.55.Qr

Three-dimensional gamma analysis of dose distributions in individual structures for IMRT dose verification

Our purpose in this study was to implement three-dimensional (3D) gamma analysis for structures of interest such as the planning target volume (PTV) or clinical target volume (CTV), and organs at risk (OARs) for intensity-modulated radiation therapy (IMRT) dose verification. IMRT dose distributions for prostate and head and neck (HN) cancer patients were calculated with an analytical anisotropic algorithm in an Eclipse (Varian Medical Systems) treatment planning system (TPS) and by Monte Carlo (MC) simulation. The MC dose distributions were calculated with EGSnrc/BEAMnrc and DOSXYZnrc user codes under conditions identical to those for the TPS. The prescribed doses were 76 Gy/38 fractions with five-field IMRT for the prostate and 33 Gy/17 fractions with seven-field IMRT for the HN. TPS dose distributions were verified by the gamma passing rates for the whole calculated volume, PTV or CTV, and OARs by use of 3D gamma analysis with reference to MC dose distributions. The acceptance criteria for the 3D gamma analysis were 3/3 and 2 %/2 mm for a dose difference and a distance to agreement. The gamma passing rates in PTV and OARs for the prostate IMRT plan were close to 100 %. For the HN IMRT plan, the passing rates of 2 %/2 mm in CTV and OARs were substantially lower because inhomogeneous tissues such as bone and air in the HN are included in the calculation area. 3D gamma analysis for individual structures is useful for IMRT dose verification.

Performance of dose calculation algorithms from three generations in lung SBRT: comparison with full Monte Carlo-based dose distributions

The accuracy of dose calculation is a key challenge in stereotactic body radiotherapy (SBRT) of the lung. We have benchmarked three photon beam dose calculation algorithms — pencil beam convolution (PBC), anisotropic analytical algorithm (AAA), and Acuros XB (AXB) — implemented in a commercial treatment planning system (TPS), Varian Eclipse. Dose distributions from full Monte Carlo (MC) simulations were regarded as a reference. In the first stage, for four patients with central lung tumors, treatment plans using 3D conformal radiotherapy (CRT) technique applying 6 MV photon beams were made using the AXB algorithm, with planning criteria according to the Nordic SBRT study group. The plans were recalculated (with same number of monitor units (MUs) and identical field settings) using BEAMnrc and DOSXYZnrc MC codes. The MC‐calculated dose distributions were compared to corresponding AXB‐calculated dose distributions to assess the accuracy of the AXB algorithm, to which then other TPS algorithms were compared. In the second stage, treatment plans were made for ten patients with 3D CRT technique using both the PBC algorithm and the AAA. The plans were recalculated (with same number of MUs and identical field settings) with the AXB algorithm, then compared to original plans. Throughout the study, the comparisons were made as a function of the size of the planning target volume (PTV), using various dose‐volume histogram (DVH) and other parameters to quantitatively assess the plan quality. In the first stage also, 3D gamma analyses with threshold criteria 3%/3 mm and 2%/2 mm were applied. The AXB‐calculated dose distributions showed relatively high level of agreement in the light of 3D gamma analysis and DVH comparison against the full MC simulation, especially with large PTVs, but, with smaller PTVs, larger discrepancies were found. Gamma agreement index (GAI) values between 95.5% and 99.6% for all the plans with the threshold criteria 3%/3 mm were achieved, but 2%/2 mm threshold criteria showed larger discrepancies. The TPS algorithm comparison results showed large dose discrepancies in the PTV mean dose (D50%), nearly 60%, for the PBC algorithm, and differences of nearly 20% for the AAA, occurring also in the small PTV size range. This work suggests the application of independent plan verification, when the AAA or the AXB algorithm are utilized in lung SBRT having PTVs smaller than 20‐25 cc. The calculated data from this study can be used in converting the SBRT protocols based on type ‘a’ and/or type ‘b’ algorithms for the most recent generation type ‘c’ algorithms, such as the AXB algorithm.

PACS numbers: 87.55.‐x, 87.55.D‐, 87.55.K‐, 87.55.kd, 87.55.Qr

Software for quantitative analysis of radiotherapy: overview, requirement analysis and design solutions

Radiotherapy is a fast-developing discipline which plays a major role in cancer care. Quantitative analysis of radiotherapy data can improve the success of the treatment and support the prediction of outcome. In this paper, we first identify functional, conceptional and general requirements on a software system for quantitative analysis of radiotherapy. Further we present an overview of existing radiotherapy analysis software tools and check them against the stated requirements. As none of them could meet all of the demands presented herein, we analyzed possible conceptional problems and present software design solutions and recommendations to meet the stated requirements (e.g. algorithmic decoupling via dose iterator pattern; analysis database design). As a proof of concept we developed a software library "RTToolbox" following the presented design principles. The RTToolbox is available as open source library and has already been tested in a larger-scale software system for different use cases. These examples demonstrate the benefit of the presented design principles.

Commissioning kilovoltage cone-beam CT beams in a radiation therapy treatment planning system

The feasibility of accounting of the dose from kilovoltage cone‐beam CT in treatment planning has been discussed previously for a single cone‐beam CT (CBCT) beam from one manufacturer. Modeling the beams and computing the dose from the full set of beams produced by a kilovoltage cone‐beam CT system requires extensive beam data collection and verification, and is the purpose of this work. The beams generated by Elekta X‐ray volume imaging (XVI) kilovoltage CBCT (kV CBCT) system for various cassettes and filters have been modeled in the Philips Pinnacle treatment planning system (TPS) and used to compute dose to stack and anthropomorphic phantoms. The results were then compared to measurements made using thermoluminescent dosimeters (TLDs) and Monte Carlo (MC) simulations. The agreement between modeled and measured depth‐dose and cross profiles is within 2% at depths beyond 1 cm for depth‐dose curves, and for regions within the beam (excluding penumbra) for cross profiles. The agreements between TPS‐calculated doses, TLD measurements, and Monte Carlo simulations are generally within 5% in the stack phantom and 10% in the anthropomorphic phantom, with larger variations observed for some of the measurement/calculation points. Dose computation using modeled beams is reasonably accurate, except for regions that include bony anatomy. Inclusion of this dose in treatment plans can lead to more accurate dose prediction, especially when the doses to organs at risk are of importance.

PACS numbers: 87.55.D, 87.55.K, 87.56.bd

CTC-ask: a new algorithm for conversion of CT numbers to tissue parameters for Monte Carlo dose calculations applying DICOM RS knowledge

One of the building blocks in Monte Carlo (MC) treatment planning is to convert patient CT data to MC compatible phantoms, consisting of density and media matrices. The resulting dose distribution is highly influenced by the accuracy of the conversion. Two major contributing factors are precise conversion of CT number to density and proper differentiation between air and lung. Existing tools do not address this issue specifically. Moreover, their density conversion may depend on the number of media used. Differentiation between air and lung is an important task in MC treatment planning and misassignment may lead to local dose errors on the order of 10%. A novel algorithm, CTC-ask, is presented in this study. It enables locally confined constraints for the media assignment and is independent of the number of media used for the conversion of CT number to density. MC compatible phantoms were generated for two clinical cases using a CT-conversion scheme implemented in both CTC-ask and the DICOM-RT toolbox. Full MC dose calculation was subsequently conducted and the resulting dose distributions were compared. The DICOM-RT toolbox inaccurately assigned lung in 9.9% and 12.2% of the voxels located outside of the lungs for the two cases studied, respectively. This was completely avoided by CTC-ask. CTC-ask is able to reduce anatomically irrational media assignment. The CTC-ask source code can be made available upon request to the authors.

Patient-specific three-dimensional concomitant dose from cone beam computed tomography exposure in image-guided radiotherapy

PURPOSE

The purpose of the present study was to quantify the concomitant dose received by patients undergoing cone beam computed tomography (CBCT) scanning in different clinical scenarios as a part of image-guided radiotherapy (IGRT) procedures.

METHODS AND MATERIALS

We calculated the three-dimensional concomitant dose received as a result of CBCT scans in 6 patients representing different clinical scenarios: two pelvis, two head and neck, and two chest. We assessed the effect that a daily on-line IGRT strategy would have on the patient dose distribution, assuming 40 CBCT scans throughout the treatment course. The additional dose to the planning target volume margin region was also estimated.

RESULTS

In the pelvis, a single CBCT scan delivered a mean dose to the femoral heads of 2-6 cGy and the rectum of 1-2 cGy. An additional dose to the planning target volume was within 1-3 cGy. In the chest, the mean dose to the planning target volume varied from 2.5 to 5 cGy. The lung and spinal cord planning organ at risk volume received ≤4 cGy and ≤5 cGy, respectively. In the head and neck, a single CBCT scan delivered a mean dose of 0.3 cGy, with bony structures receiving 0.5-0.8 cGy. The femoral heads received an additional dose of 1.5-2.5 Gy. A reduction of 20-30% in the mean dose to the organs at risk was achieved using bowtie filtration. In the head and neck, the dose to the eyes and brainstem was eliminated by decreasing the craniocaudal field size.

CONCLUSIONS

The additional dose from on-line IGRT procedures can be clinically relevant. The organ dose can be significantly reduced with the use of appropriate patient-specific settings. The concomitant dose from CBCT should be accounted for and the acquisition settings optimized for optimal IGRT strategies on a patient basis.

Position-probability-sampled Monte Carlo calculation of VMAT, 3DCRT, step-shoot IMRT, and helical tomotherapy dose distributions using BEAMnrc/DOSXYZnrc

PURPOSE

The commercial release of volumetric modulated arc therapy techniques using a conventional linear accelerator and the growing number of helical tomotherapy users have triggered renewed interest in dose verification methods, and also in tools for exploring the impact of machine tolerance and patient motion on dose distributions without the need to approximate time-varying parameters such as gantry position, MLC leaf motion, or patient motion. To this end we have developed a Monte Carlo-based calculation method capable of simulating a wide variety of treatment techniques without the need to resort to discretization approximations.

METHODS

The ability to perform complete position-probability-sampled Monte Carlo dose calculations was implemented in the BEAMnrc/DOSXZYnrc user codes of EGSnrc. The method includes full accelerator head simulations of our tomotherapy and Elekta linacs, and a realistic representation of continous motion via the sampling of a time variable. The functionality of this algorithm was tested via comparisons with both measurements and treatment planning dose distributions for four types of treatment techniques: 3D conformal, step-shoot intensity modulated radiation therapy, helical tomotherapy, and volumetric modulated are therapy.

RESULTS

For static fields, the absolute dose agreement between the EGSnrc Monte Carlo calculations and measurements is within 2%/1 mm. Absolute dose agreement between Monte Carlo calculations and treatment planning system for the four different treatment techniques is within 3%/3 mm. Discrepancies with the tomotherapy TPS on the order of 10%/5 mm were observed for the extreme example of a small target located 15 cm off-axis and planned with a low modulation factor. The increase in simulation time associated with using position-probability sampling, as opposed to the discretization approach, was less than 2% in most cases.

CONCLUSIONS

A single Monte Carlo simulation method can be used to calculate patient dose distribution for various types of treatment techniques delivered with either tomotherapy or a conventional linac. The method simplifies the simulation process, improves dose calculation accuracy, and involves an acceptably small change in computation time.

A rapid communication from the AAPM Task Group 201: recommendations for the QA of external beam radiotherapy data transfer. AAPM TG 201: quality assurance of external beam radiotherapy data transfer

The transfer of radiation therapy data among the various subsystems required for external beam treatments is subject to error. Hence, the establishment and management of a data transfer quality assurance program is strongly recommended. It should cover the QA of data transfers of patient specific treatments, imaging data, manually handled data and historical treatment records. QA of the database state (logical consistency and information integrity) is also addressed to ensure that accurate data are transferred.

PACS numbers: 87.56.bd, 87.55.D, 87.55.Gh, 87.55.km, 87.55.Qr, 87.55.T

Extension of the NCAT phantom for the investigation of intra-fraction respiratory motion in IMRT using 4D Monte Carlo

The purpose of this work was to create a computational platform for studying motion in intensity modulated radiotherapy (IMRT). Specifically, the non-uniform rational B-spline (NURB) cardiac and torso (NCAT) phantom was modified for use in a four-dimensional Monte Carlo (4D-MC) simulation system to investigate the effect of respiratory-induced intra-fraction organ motion on IMRT dose distributions as a function of diaphragm motion, lesion size and lung density. Treatment plans for four clinical scenarios were designed: diaphragm peak-to-peak amplitude of 1 cm and 3 cm, and two lesion sizes--2 cm and 4 cm diameter placed in the lower lobe of the right lung. Lung density was changed for each phase using a conservation of mass calculation. Further, a new heterogeneous lung model was implemented and tested. Each lesion had an internal target volume (ITV) subsequently expanded by 15 mm isotropically to give the planning target volume (PTV). The PTV was prescribed to receive 72 Gy in 40 fractions. The MLC leaf sequence defined by the planning system for each patient was exported and used as input into the MC system. MC simulations using the dose planning method (DPM) code together with deformable image registration based on the NCAT deformation field were used to find a composite dose distribution for each phantom. These composite distributions were subsequently analyzed using information from the dose volume histograms (DVH). Lesion motion amplitude has the largest effect on the dose distribution. Tumor size was found to have a smaller effect and can be mitigated by ensuring the planning constraints are optimized for the tumor size. The use of a dynamic or heterogeneous lung density model over a respiratory cycle does not appear to be an important factor with a <or=0.6% change in the mean dose received by the ITV, PTV and right lung. The heterogeneous model increases the realism of the NCAT phantom and may provide more accurate simulations in radiation therapy investigations that use the phantom. This work further evaluates the NCAT phantom for use as a tool in radiation therapy research in addition to its extensive use in diagnostic imaging and nuclear medicine research. Our results indicate that the NCAT phantom, combined with 4D-MC simulations, is a useful tool in radiation therapy investigations and may allow the study of relative effects in many clinically relevant situations.

Monte Carlo simulation and patient dosimetry for a kilovoltage cone-beam CT unit

PURPOSE

The purpose of this work is to characterize the x-ray volume imager (XVI), the cone-beam computed tomography (CBCT) unit mounted on the Elekta Synergy linac, with F1 bowtie filter and to calculate the three-dimensional dose delivered to patients using volumetric acquisition.

METHODS

The XVI is modeled in detail using a new Monte Carlo (MC) code, BEAMPP, under development at the National Research Council Canada. In this investigation, a new component module is developed to accurately model the unit's bowtie filter used in conjunction with the available beam collimators at the clinical energy of 120 kV. The modeling is compared against percentage depth dose (PDD) and profile measurements. Kilovoltage radiation beams' phase space files are also analyzed. The authors also describe a method for the absolute dose calibration of the MC model of the CBCT unit when used in a clinical volumetric acquisition mode. Finally, they calculate three-dimensional patient dose from CBCT image acquisition in three clinical cases of interest: Pelvis, lung, and head and neck.

RESULTS

The agreement between measurement and MC is shown to be very good: Within +/- 2% for the PDD and within +/- 3.5% inside the radiation field for all the collimators with the F1 bowtie filter. A full account of the absolute calibration method is given and dose calculation is validated against ion chamber measurements in different locations of a plastic phantom. Calculations and experiments agree within +/- 2% or better in both at the center and the periphery of the phantom, with worst agreement of 4.5% at the surface of the phantom and for one specific combination of collimator and filter. Patient dose from CBCT scan reveals that dose to tissue is between 2 and 2.5 cGy for a pelvis or a lung full acquisition. For H&N dose to tissue is 5 cGy, with the unit presets used in this work. Dose to bony structures can be two to three times higher than dose to tissue.

CONCLUSIONS

The XVI CBCT unit has been fully modeled including the F1 bowtie filter. Absolute dose distribution from the unit has been successfully validated. Full MC patient dose calculation has shown that the three-dimensional dose distribution from CBCT is complex. Patient dose from CBCT exposure cannot be completely accounted for by using a numerical factor as an estimate of the dose at the center of the body. Furthermore, additional dose to bone should be taken into account when adopting any IGRT strategy and weighed vs the unquestionable benefits of the technique in order to optimize treatment. Full three-dimensional dose calculation is recommended if patient dose from CBCT is to be integrated in any adaptive planning strategy.

A preliminary study of in-house Monte Carlo simulations: an integrated Monte Carlo verification system

PURPOSE

To develop an infrastructure for the integrated Monte Carlo verification system (MCVS) to verify the accuracy of conventional dose calculations, which often fail to accurately predict dose distributions, mainly due to inhomogeneities in the patient's anatomy, for example, in lung and bone.

METHODS AND MATERIALS

The MCVS consists of the graphical user interface (GUI) based on a computational environment for radiotherapy research (CERR) with MATLAB language. The MCVS GUI acts as an interface between the MCVS and a commercial treatment planning system to import the treatment plan, create MC input files, and analyze MC output dose files. The MCVS consists of the EGSnrc MC codes, which include EGSnrc/BEAMnrc to simulate the treatment head and EGSnrc/DOSXYZnrc to calculate the dose distributions in the patient/phantom. In order to improve computation time without approximations, an in-house cluster system was constructed.

RESULTS

The phase-space data of a 6-MV photon beam from a Varian Clinac unit was developed and used to establish several benchmarks under homogeneous conditions. The MC results agreed with the ionization chamber measurements to within 1%. The MCVS GUI could import and display the radiotherapy treatment plan created by the MC method and various treatment planning systems, such as RTOG and DICOM-RT formats. Dose distributions could be analyzed by using dose profiles and dose volume histograms and compared on the same platform. With the cluster system, calculation time was improved in line with the increase in the number of central processing units (CPUs) at a computation efficiency of more than 98%.

CONCLUSIONS

Development of the MCVS was successful for performing MC simulations and analyzing dose distributions.

Analysis of outcomes in radiation oncology: an integrated computational platform

Radiotherapy research and outcome analyses are essential for evaluating new methods of radiation delivery and for assessing the benefits of a given technology on locoregional control and overall survival. In this article, a computational platform is presented to facilitate radiotherapy research and outcome studies in radiation oncology. This computational platform consists of (1) an infrastructural database that stores patient diagnosis, IMRT treatment details, and follow-up information, (2) an interface tool that is used to import and export IMRT plans in DICOM RT and AAPM/RTOG formats from a wide range of planning systems to facilitate reproducible research, (3) a graphical data analysis and programming tool that visualizes all aspects of an IMRT plan including dose, contour, and image data to aid the analysis of treatment plans, and (4) a software package that calculates radiobiological models to evaluate IMRT treatment plans. Given the limited number of general-purpose computational environments for radiotherapy research and outcome studies, this computational platform represents a powerful and convenient tool that is well suited for analyzing dose distributions biologically and correlating them with the delivered radiation dose distributions and other patient-related clinical factors. In addition the database is web-based and accessible by multiple users, facilitating its convenient application and use.

A beam-matching concept for medical linear accelerators

The flexibility in radiotherapy can be improved if a patient can be moved between any one of the department's medical linear accelerators without the need to change anything in the patient's treatment plan. For this to be possible, the dosimetric characteristics of the various accelerators must be the same, or nearly the same i.e. the accelerators must be beam-matched. During a period of nine months, eight Varian iX accelerators with 6 and 15 MV photon beams and 6-18 MeV electron beams (only four of the eight) were installed at our clinic. All accelerators fulfilled the vendor-defined "fine beam-match" criteria, and a more extensive set of measurements was carried out during commissioning. The measured absorbed dose data for each accelerator were compared with the first accelerator, chosen as reference, and the TPS calculations. Two of the eight accelerators showed a larger discrepancy for the 15 MV beam not revealed by the vendor-defined acceptance criteria, whereas the other six accelerators were satisfactorily matched. The beam-matching acceptance criteria defined by the vendor are not strict enough to guarantee optimal beam-match. Deviations related to dose calculations and to beam-matched accelerators may add up. The safest and most practical way to ensure that all accelerators are within clinical acceptable accuracy is to include TPS calculations in the evaluation. Further, comparisons between measurements and calculations should be done in absolute dose terms.

Experimental measurements and Monte Carlo simulations for dosimetric evaluations of intrafraction motion for gated and ungated intensity modulated arc therapy deliveries

Respiratory gated radiation therapy allows for a smaller margin expansion for the planning target volume (PTV) to account for respiratory induced motion and is emerging as a common method to treat lung and liver tumors. We investigated the dosimetric effect of free motion and gated delivery for intensity modulated arc therapy (IMAT) with experimental measurements and Monte Carlo simulations. The impact of PTV margin and duty cycle for gated delivery is studied with Monte Carlo simulations. A motion phantom is used for this study. Two sets of contours were drawn on the mid-inspiration CT scan of this motion phantom. For each set of contours, an IMAT plan to be delivered with constant dose rate was created. The plans were generated on a CT scan of the phantom in the static condition with 3 mm PTV margin and applied to the motion phantom under four conditions: static, full superior-inferior (SI) motion (A = 1 cm, T = 4 s) and gating conditions (25% and 50% duty cycles) with full SI motion. A 6 by 15 cm piece of radiographic film was placed in the sagittal plane of the phantom and then irradiated under all measurement conditions. Film calibration was performed with a step-wedge method to convert optical density to dose. Gated IMAT delivery was first validated in 2D by comparing static film with that from gating and full motion. A previously verified simulation tool for IMRT that takes the log files from the multileaf collimator (MLC) controller and the gating system were adapted to simulate the delivered IMAT treatment for full 3D dosimetric analysis. The IMAT simulations were validated against the 2D film measurements. The resultant IMAT simulations were evaluated with dose criteria, dose-volume histograms and 3D gamma analysis. We validated gated IMAT deliveries when we compared the static film with the one from gating using 25% duty cycle using 2D gamma analysis. Within experimental and setup uncertainties, film measurements agreed with their corresponding simulated plans using 2D gamma analysis. Finally, when planning with margins designed for gating with 25% duty cycle and applying 50% or no gating during treatment, the dose differences in D(min,) D(99%) and D(95%) of the clinical target volume can be up to 27 cGy, 20 cGy and 18 cGy, respectively, for a plan with 200 cGy prescription dose. We have experimentally delivered gated IMAT with constant dose rate to a motion phantom and assessed their accuracies with film dosimetry and Monte Carlo simulations. Film dosimetry demonstrated that 25% gating and static plans are within 2%, 2 mm. The Monte Carlo simulation method was employed to generate dose delivered in 3D to a motion phantom, and the dosimetric results were reported. Since our film measurements agreed well with Monte Carlo simulations, we can reliably use this simulation tool to further study the dosimetric effects of target motion and effectiveness of gating for IMAT deliveries.

Dosimetric and radiobiological evaluation of dose distribution perturbation due to head heterogeneities for Linac and Gamma Knife stereotactic radiotherapy

INTRODUCTION

In SRT/SRS, dedicated treatment planning systems are used for the calculation of the dose distribution. The majority of these systems utilize the standard TMR/OAR formalism for dose calculation as well as they usually neglect any perturbation due to head heterogeneities. The aim of this study is to examine the errors due to head heterogeneities for both absolute and relative dose distributions in stereotactic radiotherapy.

MATERIALS AND METHODS

Dosimetric measurements in phantoms have been made for linac stereotactic irradiation. CT-based phantoms have been used for Monte Carlo simulations for both linac-based stereotactic system and Gamma Knife unit. Absolute and relative dose distributions have been compared between homogeneous and heterogeneous media. DVH and TCP results are presented for all cases.

RESULTS

The maximum absolute dose difference at the isocenter was 2.2% and 6.9% for the linac and Gamma Knife respectively. The impact of heterogeneity in the target DVH was minor for the linac technique whereas considerable difference was observed for the Gamma Knife treatment. This was reflected also to the radiobiological evaluation, where the maximum TCP difference for the linac system was 2.7% and for the Gamma Knife was 4%.

DISCUSSION AND CONCLUSIONS

The errors rising from the existence of head heterogeneities are not negligible especially for the Gamma Knife which uses lower energy beams. The errors of the absolute dose calculation could be easily eliminated by implementing a simple heterogeneity correction algorithm at the TPS. Nevertheless, the errors for not taking into account the lateral electron transport would require a more sophisticated approach and even direct Monte Carlo calculation.

Monte Carlo simulation of RapidArc radiotherapy delivery

RapidArc radiotherapy technology from Varian Medical Systems is one of the most complex delivery systems currently available, and achieves an entire intensity-modulated radiation therapy (IMRT) treatment in a single gantry rotation about the patient. Three dynamic parameters can be continuously varied to create IMRT dose distributions-the speed of rotation, beam shaping aperture and delivery dose rate. Modeling of RapidArc technology was incorporated within the existing Vancouver Island Monte Carlo (VIMC) system (Zavgorodni et al 2007 Radiother. Oncol. 84 S49, 2008 Proc. 16th Int. Conf. on Medical Physics). This process was named VIMC-Arc and has become an efficient framework for the verification of RapidArc treatment plans. VIMC-Arc is a fully automated system that constructs the Monte Carlo (MC) beam and patient models from a standard RapidArc DICOM dataset, simulates radiation transport, collects the resulting dose and converts the dose into DICOM format for import back into the treatment planning system (TPS). VIMC-Arc accommodates multiple arc IMRT deliveries and models gantry rotation as a series of segments with dynamic MLC motion within each segment. Several verification RapidArc plans were generated by the Eclipse TPS on a water-equivalent cylindrical phantom and re-calculated using VIMC-Arc. This includes one 'typical' RapidArc plan, one plan for dual arc treatment and one plan with 'avoidance' sectors. One RapidArc plan was also calculated on a DICOM patient CT dataset. Statistical uncertainty of MC simulations was kept within 1%. VIMC-Arc produced dose distributions that matched very closely to those calculated by the anisotropic analytical algorithm (AAA) that is used in Eclipse. All plans also demonstrated better than 1% agreement of the dose at the isocenter. This demonstrates the capabilities of our new MC system to model all dosimetric features required for RapidArc dose calculations.

Dosimetric and geometric evaluation of an open low-field magnetic resonance simulator for radiotherapy treatment planning of brain tumours

BACKGROUND AND PURPOSE

Magnetic resonance (MR) imaging is superior to computed tomography (CT) in radiotherapy of brain tumours. In this study an open low-field MR-simulator is evaluated in order to eliminate the cost of and time spent on additional CT scanning.

MATERIALS AND METHODS

Eleven patients with brain tumours are both CT and MR scanned and the defined tumour volumes are compared. Image distortions and dose calculations based on CT density correction, MR unit density and MR bulk density, bone segmentation are performed. Monte Carlo simulations using 4 and 8 MV beams on homogeneous and bone segmented mediums are performed.

RESULTS

Mean MR and CT tumour volumes of approximately the same size (V MR =55+/-34 cm3 and V CT =51+/-32 cm3) are observed, but for individual patients, small intersection volumes are observed. The MR images show negligible distortion within radial distances below 12 cm (<1.5 mm). On unit density mediums, dose errors above 2% are observed in low dose areas. Monte Carlo simulations with 4 MV photons show large deviations in dose (>2%) just behind the skull if bone is not segmented.

CONCLUSIONS

It is feasible to use an MR-simulator for radiotherapy planning of brain tumours if bone is segmented or a careful choice of beam energy (>4 MV) is selected.

Interface software for DOSXYZnrc Monte Carlo dose evaluation on a commercial radiation treatment planning system

PURPOSE

As the conventional graphical user interface (GUI) associated with DOSXYZnrc or BEAMnrc is unable to define specific structures such as gross tumor volume (GTV) on computed tomography (CT) data, the quantitative analysis of doses in the form of dose-volume histograms (DVHs) is difficult. The purpose of this study was to develop an interface that enables us to analyze the results of DOSXYZnrc output with a commercial radiation treatment planning (RTP) system and to investigate the validity of the system.

MATERIALS AND METHODS

Interface software to visualize three-dimensional radiotherapy Monte Carlo (MC) dose data from DOSXYZnrc on the XiO RTP system was developed. To evaluate the interface, MC doses for a variety of photon energies were calculated using the CT data of a thorax phantom and a uniform phantom as well as data from patients with lung tumors.

RESULTS

The dose files were analyzed on the XiO RTP system in the form of isodose distributions and DVHs. In all cases, the XiO RTP system perfectly displayed the MC doses for quantitative evaluation in the form of differential and integral DVHs.

CONCLUSION

Three-dimensional display of DOSXYZnrc doses on a dedicated RTP system could provide all the existing facilities of the system for quantitative dose analysis.

A virtual-accelerator-based verification of a Monte Carlo dose calculation algorithm for electron beam treatment planning in clinical situations

BACKGROUND AND PURPOSE

The introduction of Monte Carlo (MC) techniques for treatment planning and also for verification purposes will have considerable impact on the radiation therapy planning process. The aim of this work was to use a virtual accelerator to study the performance of a MC-based electron dose calculation algorithm, implemented in a commercial treatment planning system.

METHODS

The performance in phantoms containing air and bone as well as in patient-specific geometries (thorax wall, nose, parotid gland and spinal cord) has been studied.

RESULTS

The agreement between the virtual accelerator and the MC dose calculation algorithm is generally very good. A gamma-evaluation with criteria of 0.03 Gy/3 mm (per Gy at the depth of maximum dose) shows that, even for the worst cases, only a small volume of about 1.5% has gamma>1.0. In the worst case, with the 0.02 Gy/2 mm criteria, about 92% of the volume receiving more than 0.85 Gy per 100 monitor units (MU) has gamma-values <1.0. The corresponding value for the volume receiving more than 0.10 Gy/100 MU is about 98%. For the 18 MeV spinal-cord case, where a 6 x 20 cm2 insert is used, the TPS underestimates the dose outside the primary field due to inadequate modelling of the insert.

CONCLUSION

The possibility of dose calculations in typical patient cases makes the virtual accelerator a powerful tool for validation and evaluation of dose calculation algorithms present in treatment planning systems.

Coordinate transformations for BEAM/EGSnrc Monte Carlo dose calculations of non-coplanar fields received from a DICOM-compliant treatment planning system

The Monte Carlo (MC) method provides the most accurate to-date dose calculations in heterogeneous media and complex geometries, and this spawns increasing interest in incorporating MC calculations in the treatment planning quality assurance process. This process involves MC dose calculations for the treatment plans produced clinically. Commonly used in radiotherapy, MC codes are BEAMnrc and DOSXYZnrc, which transport particles in a coordinate system (c.s.) that has been established historically and does not correspond to the c.s. of treatment planning systems (TPSs). Relative rotations of these c.s. are not straightforward, especially for non-coplanar treatments. Transformation equations are therefore required to re-calculate a treatment plan using BEAM/DOSXYZnrc codes. This paper presents such transformations for beam angles defined in a DICOM-compliant treatment planning coordinate system. Verification of the derived transformations with two three-field plans simulated on a phantom using TPS as well as MC codes has been performed demonstrating exact geometrical agreement of the MC treatment fields' placement.

Dosimetric verification for intensity-modulated radiotherapy of thoracic cancers using experimental and Monte Carlo approaches

PURPOSE

To investigate the dosimetric accuracy of commercial treatment planning systems used in intensity-modulated radiotherapy (IMRT) for thoracic cancer.

METHODS AND MATERIALS

Clinical IMRT plans for lung and esophageal cancers and mesothelioma were used to investigate the accuracy of dose calculations from two commercial treatment planning systems (Pinnacle and Corvus systems). Dose distributions were measured with ion chambers and thermoluminescent dosimeters for individual IMRT fields and composite treatment plans in water phantoms and anthropomorphic phantoms. A Monte Carlo-based system was established to compute three-dimensional dose distributions to compare with the treatment planning system calculations.

RESULTS

Dose calculations from the Pinnacle system were acceptable within 5% of the local dose or a 5-mm distance-to-agreement for 80% of the points measured with ion chambers, 74% of the points measured with thermoluminescent dosimeters, and 96% of the points compared with the Monte Carlo calculations. For the Corvus system, 89% of the points agreed with the measured dose and 98% agreed with the Monte Carlo calculations. Underestimation of the dose from the treatment planning system was found in the low-dose regions (<50% of the prescribed dose), possibly caused by inadequate modeling of the multileaf collimators.

CONCLUSION

The Pinnacle and Corvus dose calculations were acceptable for thoracic IMRT in high-dose regions. Beam modeling is likely the most critical factor for the accuracy of IMRT dose calculations.

Comparison of dose calculation algorithms for treatment planning in external photon beam therapy for clinical situations

A study of the performance of five commercial radiotherapy treatment planning systems (TPSs) for common treatment sites regarding their ability to model heterogeneities and scattered photons has been performed. The comparison was based on CT information for prostate, head and neck, breast and lung cancer cases. The TPSs were installed locally at different institutions and commissioned for clinical use based on local procedures. For the evaluation, beam qualities as identical as possible were used: low energy (6 MV) and high energy (15 or 18 MV) x-rays. All relevant anatomical structures were outlined and simple treatment plans were set up. Images, structures and plans were exported, anonymized and distributed to the participating institutions using the DICOM protocol. The plans were then re-calculated locally and exported back for evaluation. The TPSs cover dose calculation techniques from correction-based equivalent path length algorithms to model-based algorithms. These were divided into two groups based on how changes in electron transport are accounted for ((a) not considered and (b) considered). Increasing the complexity from the relatively homogeneous pelvic region to the very inhomogeneous lung region resulted in less accurate dose distributions. Improvements in the calculated dose have been shown when models consider volume scatter and changes in electron transport, especially when the extension of the irradiated volume was limited and when low densities were present in or adjacent to the fields. A Monte Carlo calculated algorithm input data set and a benchmark set for a virtual linear accelerator have been produced which have facilitated the analysis and interpretation of the results. The more sophisticated models in the type b group exhibit changes in both absorbed dose and its distribution which are congruent with the simulations performed by Monte Carlo-based virtual accelerator.

Gamma histograms for radiotherapy plan evaluation

BACKGROUND AND PURPOSE

The technique known as the 'gamma evaluation method' incorporates pass-fail criteria for both distance-to-agreement and dose difference analysis of 3D dose distributions and provides a numerical index (gamma) as a measure of the agreement between two datasets. As the gamma evaluation index is being adopted in more centres as part of treatment plan verification procedures for 2D and 3D dose maps, the development of methods capable of encapsulating the information provided by this technique is recommended.

PATIENTS AND METHODS

In this work the concept of gamma index was extended to create gamma histograms (GH) in order to provide a measure of the agreement between two datasets in two or three dimensions. Gamma area histogram (GAH) and gamma volume histogram (GVH) graphs were produced using one or more 2D gamma maps generated for each slice of the irradiated volume. GHs were calculated for IMRT plans, evaluating the 3D dose distribution from a commercial treatment planning system (TPS) compared to a Monte Carlo (MC) calculation used as reference dataset.

RESULTS

The extent of local anatomical inhomogenities in the plans under consideration was strongly correlated with the level of difference between reference and evaluated calculations. GHs provided an immediate visual representation of the proportion of the treated volume that fulfilled the gamma criterion and offered a concise method for comparative numerical evaluation of dose distributions.

CONCLUSIONS

We have introduced the concept of GHs and investigated its applications to the evaluation and verification of IMRT plans. The gamma histogram concept set out in this paper can provide a valuable technique for quantitative comparison of dose distributions and could be applied as a tool for the quality assurance of treatment planning systems.

Monte Carlo dosimetric evaluation of high energy vs low energy photon beams in low density tissues

BACKGROUND AND PURPOSE

Low megavoltage photon beams are often the treatment choice in radiotherapy when low density heterogeneities are involved, because higher energies show some undesirable dosimetric effects. This work is aimed at investigating the effects of different energy selection for low density tissues.

PATIENTS AND METHODS

BEAMnrc was used to simulate simple treatment set-ups in a simple and a CT reconstructed lung phantom and an air-channel phantom. The dose distribution of 6, 15 and 20 MV photon beams was studied using single, AP/PA and three-field arrangements.

RESULTS

Our results showed no significant changes in the penumbra width in lung when a pair of opposed fields were used. The underdosage at the anterior/posterior tumor edge caused by the dose build-up at the lung-tumor interface reached 7% for a 5 x 5 cm AP/PA set-up. Shrinkage of the 90% isodose volume was noticed for the same set-up, which could be rectified by adding a lateral field. For the CT reconstructed phantom, the AP/PA set-up offered better tumor coverage when lower energies were used but for the three field set-up, higher energies resulted to better sparing of the lung tissue. For the air-channel set-up, adding an opposed field reduced the penumbra width. Using higher energies resulted in a 7% cold spot around the air-tissue interface for a 5 x 5 cm field.

CONCLUSIONS

The choice of energy for treatment in the low density areas is not a straightforward decision but depends on a number of parameters such as the beam set-up and the dosimetric criteria. Updated calculation algorithms should be used in order to be confident for the choice of energy of treatment.

A virtual-accelerator-based verification of a Monte Carlo dose calculation algorithm for electron beam treatment planning in homogeneous phantoms

By introducing Monte Carlo (MC) techniques to the verification procedure of dose calculation algorithms in treatment planning systems (TPSs), problems associated with conventional measurements can be avoided and properties that are considered unmeasurable can be studied. The aim of the study is to implement a virtual accelerator, based on MC simulations, to evaluate the performance of a dose calculation algorithm for electron beams in a commercial TPS. The TPS algorithm is MC based and the virtual accelerator is used to study the accuracy of the algorithm in water phantoms. The basic test of the implementation of the virtual accelerator is successful for 6 and 12 MeV (gamma < 1.0, 0.02 Gy/2 mm). For 18 MeV, there are problems in the profile data for some of the applicators, where the TPS underestimates the dose. For fields equipped with patient-specific inserts, the agreement is generally good. The exception is 6 MeV where there are slightly larger deviations. The concept of the virtual accelerator is shown to be feasible and has the potential to be a powerful tool for vendors and users.