Automated testing platform for radiotherapy treatment planning scripts.

Realizing the potential of user-developed automation software interacting with a treatment planning system (TPS) requires rigorous testing to ensure patient safety and data integrity. We developed an automated test platform to allow comparison of the treatment planning database before and after the execution of a write-enabled script interacting with a commercial TPS (Eclipse, Varian Medical Systems, Palo Alto, CA) using the vendor-provided Eclipse Scripting Application Programming Interface (ESAPI). The C#-application known as Write-Enable Script Testing Engine (WESTE) serializes the treatment planning objects (Patient, Structure Set, PlanSetup) accessible through ESAPI, and then compares the serialization acquired before and after the execution of the script being tested, documenting identified differences to highlight the changes made to the treatment planning data. The first two uses of WESTE demonstrated that the testing platform could acquire and analyze the data quickly (<4 s per test case) and facilitate the clinical implementation of write-enabled scripts.

A fast GPU-accelerated Monte Carlo engine for calculation of MLC-collimated electron fields.

BACKGROUND: Although intensity-modulated radiation therapy and volumetric arc therapy have revolutionized photon external beam therapies, the technological advances associated with electron beam therapy have fallen behind. Modern linear accelerators contain technologies that would allow for more advanced forms of electron treatments, such as beam collimation, using the conventional photon multi-leaf collimator (MLC); however, no commercial solutions exist that calculate dose from such beam delivery modes. Additionally, for clinical adoption to occur, dose calculation times would need to be on par with that of modern dose calculation algorithms. PURPOSE: This work developed a graphics processing unit (GPU)-accelerated Monte Carlo (MC) engine incorporating the Varian TrueBeam linac head geometry for a rapid calculation of electron beams collimated using the conventional photon MLC. METHODS: A compute unified device architecture framework was created for the following: (1) transport of electrons and photons through the linac head geometry, considering multiple scattering, Bremsstrahlung, Møller, Compton, and pair production interactions; (2) electron and photon propagation through the CT geometry, considering all interactions plus the photoelectric effect; and (3) secondary particle cascades through the linac head and within the CT geometry. The linac head collimating geometry was modeled according to the specifications provided by the vendor, who also provided phase-space files. The MC was benchmarked against EGSnrc/DOSXYZnrc/GEANT by simulating individual interactions with simple geometries, pencil, and square beam dose calculations in various phantoms. MC-calculated dose distributions for MLC and jaw-collimated electron fields were compared to measurements in a water phantom and with radiochromic film. RESULTS: Pencil and square beam dose distributions are in good agreement with DOSXYZnrc. Angular and spatial distributions for multiple scattering and secondary particle production in thin slab geometries are in good agreement with EGSnrc and GEANT. Dose profiles for MLC and jaw-collimated 6-20-MeV electron beams showed an average absolute difference of 1.1 and 1.9 mm for the FWHM and 80%-20% penumbra from measured profiles. Percent depth doses showed differences of <5% for as compared to measurement. The computation time on an NVIDIA Tesla V100 card was 2.5 min to achieve a dose uncertainty of <1%, which is ∼300 times faster than published results in a similar geometry using a single-CPU core. CONCLUSIONS: The GPU-based MC can quickly calculate dose for electron fields collimated using the conventional photon MLC. The fast calculation times will allow for a rapid calculation of electron fields for mixed photon and electron particle therapy.

Phantom Verification of AAA and Acuros Dose Calculations for Lung Cancer: Do Tumor Size and Regression Matter?

PURPOSE: This study aimed to evaluate dose calculation accuracy for the Eclipse Analytical Anisotropic Algorithm (AAA) and Acuros XB algorithm for various lung tumor sizes and to investigate dosimetric changes associated with treatment of regressing tumors. METHODS AND MATERIALS: A water phantom with cylindrical cork inserts (lung surrogates) was fabricated. Large (202 cm3), medium (54 cm3), and small (3 cm3) solid water tumors were implanted within cork inserts. A plain cork insert was used to simulate a lung without a tumor. The cork inserts and tumors were cut along the long axis, and Gafchromic film was placed between the sections to measure dose distributions. Three-dimensional conformal plans were created using 6 MV and 10 MV beams, and volumetric modulated arc therapy plans were created using 6 MV beams for each tumor size. Doses were calculated using Eclipse AAA and Acuros XB. The measured and calculated dose distributions were compared for each tumor size and treatment algorithm. To simulate a regressing tumor, the original plans created for the large tumor were separately delivered to the phantom that contained a small, medium, or no tumor. The dosimetric effects were evaluated using gamma passing rates with a 2%/2 mm criterion and dose profile comparisons. RESULTS: Agreement between the measurements and AAA calculations decreased as tumor size decreased, but Acuros XB showed better agreement for all tumor sizes. The largest difference was observed for a 6 MV volumetric modulated arc therapy plan created to treat the smallest tumor. The gamma passing rate was 89.7% but that of Acuros was 99.5%. For the tumor regression evaluation, the gamma passing rates ranged from 53% to 99% for AAA. For Acuros XB, the gamma passing rates were >98% for all scenarios. CONCLUSION: Both AAA and Acuros XB calculated the dose accurately for the largest lung tumor. For the smallest and regressing tumors, Acuros XB more accurately modelled the dose distribution compared with AAA.

Treatment planning for metals using an extended CT number scale.

Metal implants which saturate the CT number scale may require dosimetrist and physicist involvement to manually contour and assign an appropriate value to the metal for accurate dose calculation. This study investigated dose calculation based directly on extended CT scale images for different metals and geometries. The aim was to evaluate extended CT accuracy as a suitable alternative to standard CT methods in the presence of high-Z materials and artifacts, despite the reduced HU resolution of extended CT. Gafchromic film measurements were made for comparison to calculated doses. The method of direct dose calculation on extended CT scale was compared to our institution's standard method of manually contouring and assigning metal values on saturated CT images for each of the metal samples. Clinical patient plans with metal implants were investigated and DVHs were compared between standard CT and extended CT dose calculations. Dose calculations showed agreement within 2% between the two methods of metal characterization and the film measurement in the case of the strongest metal attenuator, cobalt-chromium. In the clinical treatment plans, the greatest dose discrepancy between the two methods was 1.2%. This study suggests that direct dose calculation on an extended scale CT image in the presence of metal implants can produce accurate clinically viable treatment plans, thereby improving efficiency of clinical workflow and eliminating a potential source of human error by manual CT number assignment.

Evaluation of the Eclipse eMC algorithm for bolus electron conformal therapy using a standard verification dataset.

The purpose of this study was to evaluate the accuracy and calculation speed of electron dose distributions calculated by the Eclipse electron Monte Carlo (eMC) algorithm for use with bolus electron conformal therapy (ECT). The recent com-mercial availability of bolus ECT technology requires further validation of the eMC dose calculation algorithm. eMC-calculated electron dose distributions for bolus ECT have been compared to previously measured TLD-dose points throughout patient-based cylindrical phantoms (retromolar trigone and nose), whose axial cross sections were based on the mid-PTV (planning treatment volume) CT anatomy. The phantoms consisted of SR4 muscle substitute, SR4 bone substitute, and air. The treatment plans were imported into the Eclipse treatment planning system, and electron dose distributions calculated using 1% and < 0.2% statistical uncertainties. The accuracy of the dose calculations using moderate smoothing and no smooth-ing were evaluated. Dose differences (eMC-calculated less measured dose) were evaluated in terms of absolute dose difference, where 100% equals the given dose, as well as distance to agreement (DTA). Dose calculations were also evaluated for calculation speed. Results from the eMC for the retromolar trigone phantom using 1% statistical uncertainty without smoothing showed calculated dose at 89% (41/46) of the measured TLD-dose points was within 3% dose difference or 3 mm DTA of the measured value. The average dose difference was -0.21%, and the net standard deviation was 2.32%. Differences as large as 3.7% occurred immediately distal to the mandible bone. Results for the nose phantom, using 1% statistical uncertainty without smoothing, showed calculated dose at 93% (53/57) of the measured TLD-dose points within 3% dose difference or 3 mm DTA. The average dose difference was 1.08%, and the net standard deviation was 3.17%. Differences as large as 10% occurred lateral to the nasal air cavities. Including smoothing had insignificant effects on the accuracy of the retromolar trigone phantom calculations, but reduced the accuracy of the nose phantom calculations in the high-gradient dose areas. Dose calculation times with 1% statistical uncertainty for the retromolar trigone and nose treatment plans were 30 s and 24 s, respectively, using 16 processors (Intel Xeon E5-2690, 2.9 GHz) on a framework agent server (FAS). In comparison, the eMC was significantly more accurate than the pencil beam algorithm (PBA). The eMC has comparable accuracy to the pencil beam redefinition algorithm (PBRA) used for bolus ECT planning and has acceptably low dose calculation times. The eMC accuracy decreased when smoothing was used in high-gradient dose regions. The eMC accuracy was consistent with that previously reported for accuracy of the eMC electron dose algorithm and shows that the algorithm is suitable for clinical implementation of bolus ECT.

Cadaveric verification of the Eclipse AAA algorithm for spine SBRT treatments with titanium hardware.

PURPOSE: To assess the accuracy of the Eclipse Analytical Anisotropic Algorithm when calculating dose for spine stereotactic body radiation therapy treatments involving surgically implanted titanium hardware. METHODS AND MATERIALS: A human spine was removed from a cadaver, cut sagittally along the midline, and then separated into thoracic and lumbar sections. The thoracic section was implanted with titanium stabilization hardware; the lumbar section was not implanted. Spine sections were secured in a water phantom and simulated for treatment planning using both standard and extended computed tomography (CT) scales. Target volumes were created on both spine sections. Dose calculations were performed using (1) the standard CT scale with relative electron density (RED) override of image artifacts and hardware, (2) the extended CT scale with RED override of image artifacts only, and (3) the standard CT scale with no RED overrides for hardware or artifacts. Plans were delivered with volumetric modulated arc therapy using a 6-MV beam with and without a flattening filter. A total of 3 measurements for each plan were made with Gafchromic film placed between the spine sections and compared with Eclipse dose calculations using gamma analysis with a 2%/2 mm passing criteria. A single measurement in a homogeneous phantom was made for each plan before actual delivery. RESULTS: Gamma passing rates for measurements in the homogeneous phantom were 99.6% or greater. Passing rates for measurements made in the lumbar spine section without hardware were 99.3% or greater; measurements made in the thoracic spine containing titanium were 98.6 to 99.5%. CONCLUSIONS: Eclipse Analytical Anisotropic Algorithm can adequately model the effects of titanium implants for spine stereotactic body radiation therapy treatments using volumetric modulated arc therapy. Calculations with standard or extended CT scales give similarly accurate results.

Monitor unit calculations for external photon and electron beams: Report of the AAPM Therapy Physics Committee Task Group No. 71.

A protocol is presented for the calculation of monitor units (MU) for photon and electron beams, delivered with and without beam modifiers, for constant source-surface distance (SSD) and source-axis distance (SAD) setups. This protocol was written by Task Group 71 of the Therapy Physics Committee of the American Association of Physicists in Medicine (AAPM) and has been formally approved by the AAPM for clinical use. The protocol defines the nomenclature for the dosimetric quantities used in these calculations, along with instructions for their determination and measurement. Calculations are made using the dose per MU under normalization conditions, D'0, that is determined for each user's photon and electron beams. For electron beams, the depth of normalization is taken to be the depth of maximum dose along the central axis for the same field incident on a water phantom at the same SSD, where D'0 = 1 cGy/MU. For photon beams, this task group recommends that a normalization depth of 10 cm be selected, where an energy-dependent D'0 ≤ 1 cGy/MU is required. This recommendation differs from the more common approach of a normalization depth of dm, with D'0 = 1 cGy/MU, although both systems are acceptable within the current protocol. For photon beams, the formalism includes the use of blocked fields, physical or dynamic wedges, and (static) multileaf collimation. No formalism is provided for intensity modulated radiation therapy calculations, although some general considerations and a review of current calculation techniques are included. For electron beams, the formalism provides for calculations at the standard and extended SSDs using either an effective SSD or an air-gap correction factor. Example tables and problems are included to illustrate the basic concepts within the presented formalism.

Total skin electron therapy in the lying-on-the-floor position using a customized flattening filter to eliminate field junctions.

A total skin electron (TSE) floor technique is presented for treating patients who are unable to safely stand for extended durations. A customized flattening filter is used to eliminate the need for field junctioning, improve field uniformity, and reduce setup time. The flattening filter is constructed from copper and polycarbonate, fits into the linac's accessory slot, and is optimized to extend the useful height and width of the beam such that no field junctions are needed during treatment. A TSE floor with flattening filter (TSE FF) treatment course consisted of six patient positions: three supine and three prone. For all treatment fields, electron beam energy was 6 MeV; collimator settings were an x of 30 cm, y of 40 cm, and θcoll of 0°; and a 0.4 cm thick polycarbonate spoiler was positioned in front of the patient. Percent depth dose (PDD) and photon contamination for the TSE FF technique were compared with our standard technique, which is similar to the Stanford technique. Beam profiles were measured using radiochromic film, and dose uniformity was verified using an anthropomorphic radiological phantom. The TSE FF technique met field uniformity requirements specified by the American Association of Physicists in Medicine Task Group 30. TSE FF R80 ranges from 4 to 4.8 mm. TSE FF photon contamination was ~ 3%. Anthropomorphic radiological phantom verification demonstrated that dose to the entire skin surface was expected to be within about ±15% of the prescription dose, except for the perineum, scalp vertex, top of shoulder, and soles of the feet. The TSE floor technique presented herein eliminates field junctioning, is suitable for patients who cannot safely stand during treatment, and provides comparable quality and uniformity to the Stanford technique.

Mathematical solutions of the TG-43 geometry function for curved line, ring, disk, sphere, dome and annulus sources, and applications for quality assurance.

Analytic solutions for the TG-43 geometry function for curved line, ring, disk, sphere, dome and annulus shapes containing uniform distributions of air-kerma are derived. These geometry functions describe how dose distributions vary strictly due to source geometry and not including attenuation or scatter effects. This work extends the use of geometry functions for individual sources to applicators containing multiple sources. Such geometry functions may be used to verify dose distributions computed using advanced techniques, including QA of model-based dose calculation algorithms. The impact of source curvature on linear and planar implants is considered along with the specific clinical case of brachytherapy eye plaques. For eye plaques, the geometry function for a domed distribution is used with published Monte Carlo dose distributions to determine a radial dose function and anisotropy function which includes all the scatter and attenuation effects due to the phantom, eye plaque and sources. This TG-43 model of brachytherapy eye plaques exactly reproduces azimuthally averaged Monte Carlo calculations, both inside and outside the eye.

Recommendations for clinical electron beam dosimetry: supplement to the recommendations of Task Group 25.

The goal of Task Group 25 (TG-25) of the Radiation Therapy Committee of the American Association of.Physicists in Medicine (AAPM) was to provide a methodology and set of procedures for a medical physicist performing clinical electron beam dosimetry in the nominal energy range of 5-25 MeV. Specifically, the task group recommended procedures for acquiring basic information required for acceptance testing and treatment planning of new accelerators with therapeutic electron beams. Since the publication of the TG-25 report, significant advances have taken place in the field of electron beam dosimetry, the most significant being that primary standards laboratories around the world have shifted from calibration standards based on exposure or air kerma to standards based on absorbed dose to water. The AAPM has published a new calibration protocol, TG-51, for the calibration of high-energy photon and electron beams. The formalism and dosimetry procedures recommended in this protocol are based on the absorbed dose to water calibration coefficient of an ionization chamber at 60Co energy, N60Co(D,w), together with the theoretical beam quality conversion coefficient k(Q) for the determination of absorbed dose to water in high-energy photon and electron beams. Task Group 70 was charged to reassess and update the recommendations in TG-25 to bring them into alignment with report TG-51 and to recommend new methodologies and procedures that would allow the practicing medical physicist to initiate and continue a high quality program in clinical electron beam dosimetry. This TG-70 report is a supplement to the TG-25 report and enhances the TG-25 report by including new topics and topics that were not covered in depth in the TG-25 report. These topics include procedures for obtaining data to commission a treatment planning computer, determining dose in irregularly shaped electron fields, and commissioning of sophisticated special procedures using high-energy electron beams. The use of radiochromic film for electrons is addressed, and radiographic film that is no longer available has been replaced by film that is available. Realistic stopping-power data are incorporated when appropriate along with enhanced tables of electron fluence data. A larger list of clinical applications of electron beams is included in the full TG-70 report available at http://www.aapm.org/pubs/reports. Descriptions of the techniques in the clinical sections are not exhaustive but do describe key elements of the procedures and how to initiate these programs in the clinic. There have been no major changes since the TG-25 report relating to flatness and symmetry, surface dose, use of thermoluminescent dosimeters or diodes, virtual source position designation, air gap corrections, oblique incidence, or corrections for inhomogeneities. Thus these topics are not addressed in the TG-70 report.

Influence of source parameters on large-field electron beam profiles calculated using Monte Carlo methods.

The purpose of this paper was to study the source model for a Monte Carlo simulation of electron beams from a medical linear accelerator. In a prior study, a non-divergent Gaussian source with a full-width at half-maximum (FWHM) of 0.15 cm was successful in predicting relative dose distributions for electron beams with applicators. However, for large fields with the applicator removed, discrepancies were found between measured and calculated profiles, particularly in the shoulder region. In this work, the source was changed to a divergent Gaussian spatial distribution and the FWHM parameter was varied to produce better agreement with measured data. The influence of the FWHM source parameter on profiles was observed at multiple locations in the simulation geometry including in-air fluence profiles at a 95 cm source-to-surface distance (SSD), percent depth dose profiles and off-axis profiles (OARs) in a water phantom for two SSDs, 80 and 100 cm. For a 6 MeV 40 x 40 cm(2) OAR profile, discrepancies in the shoulder region were reduced from 15% to 4% using a FWHM value of 0.45 cm. The optimal FWHM values for the other energies were 0.45 cm for 9 MeV, 0.22 for 12 MeV, 0.25 for 16 MeV and 0.2 cm for 20 MeV. Although this range of values was larger than measured focal spot sizes reported by other researchers, using the increased FWHM values improved the fit at most locations in the simulation geometry, giving confidence that the model could be used with a variety of SSDs and field sizes.

Calculating percent depth dose with the electron pencil-beam redefinition algorithm.

In the present work, we investigated the accuracy of the electron pencil-beam redefinition algorithm (PBRA) in calculating central-axis percent depth dose in water for rectangular fields. The PBRA energy correction factor C(E) was determined so that PBRA-calculated percent depth dose best matched the percent depth dose measured in water. The hypothesis tested was that a method can be implemented into the PBRA that will enable the algorithm to calculate central-axis percent depth dose in water at a 100-cm source-to-surface distance (SSD) with an accuracy of 2% or 1-mm distance to agreement for rectangular field sizes > or = 2 x 2 cm. Preliminary investigations showed that C(E), determined using a single percent depth dose for a large field (that is, having side-scatter equilibrium), was insufficient for the PBRA to accurately calculate percent depth dose for all square fields > or = 2 x 2 cm. Therefore, two alternative methods for determining C(E) were investigated. In Method 1, C(E), modeled as a polynomial in energy, was determined by fitting the PBRA calculations to individual rectangular-field percent depth doses. In Method 2, C(E) for square fields, described by a polynomial in both energy and side of square W [that is, C = C(E,W)], was determined by fitting the PBRA calculations to measured percent depth dose for a small number of square fields. Using the function C(E,W), C(E) for other square fields was determined, and C(E) for rectangular field sizes was determined using the geometric mean of C(E) for the two measured square fields of the dimension of the rectangle (square root method). Using both methods, PBRA calculations were evaluated by comparison with measured square-field and derived rectangular-field percent depth doses at 100-cm SSD for the Siemens Primus radiotherapy accelerator equipped with a 25 x 25-cm applicator at 10 MeV and 15 MeV. To improve the fit of C(E) and C(E,W) to the electron component of percent depth dose, it was necessary to modify the PBRA's photon depth dose model to include dose buildup. Results showed that, using both methods, the PBRA was able to predict percent depth dose within criteria for all square and rectangular fields. Results showed that second- or third-order polynomials in energy (Methods 1 and 2) and in field size (Method 2) were typically required. Although the time for dose calculation using Method 1 is approximately twice that using Method 2, we recommend that Method 1 be used for clinical implementation of the PBRA because it is more accurate (most measured depth doses predicted within approximately 1%) and simpler to implement.

Evaluation of the accuracy of fetal dose estimates using TG-36 data.

The American Association of Physicists in Medicine Radiation Therapy Committee Task Group 36 report (TG-36) provides guidelines for managing radiation therapy of pregnant patients. Included in the report are data that can be used to estimate the dose to the fetus. The purpose of this study is to evaluate the accuracy of these fetal dose estimates as compared to clinically measured values. TG-36 calculations were performed and compared with measurements of the fetal dose made in vivo or in appropriately-designed phantoms. Calculation and measurement data was collected for eight pregnant patients who underwent radiation therapy at the MD Anderson Cancer Center as well as for several fetal dose studies in the literature. The maximum measured unshielded fetal dose was 47 cGy, which was 1.5% of the prescription dose. For all cases, TG-36 calculations and measured fetal doses differed by up to a factor of 3--the ratio of the calculated to measured dose ranged from 0.34 to 2.93. On average, TG-36 calculations underestimated the measured dose by 31%. No significant trends in the relationship between the calculated and measured fetal doses were found based on the distance from, or the size of, the treatment field.

Application of the electron pencil beam redefinition algorithm to electron arc therapy.

This project investigated the potential of summing fixed-beam dose distributions calculated using the pencil-beam redefinition algorithm (PBRA) at small angular steps (1 degree) to model an electron arc therapy beam. The PRBA, previously modified to model skin collimation, was modified further by incorporating two correction factors. One correction factor that is energy, SSD (source-to-surface distance), and field-width dependent constrained the calculated dose output to be the same as the measured dose output for fixed-beam geometries within the range of field widths and SSDs encountered in arc therapy. Another correction factor (single field-width correction factor for each energy) compensated for large-angle scattering not being modeled, allowing a more accurate calculation of dose output at mid arc. The PBRA was commissioned to accurately calculate dose in a water phantom for fixed-beam geometries typical of electron arc therapy. Calculated central-axis depth doses agreed with measured doses to within 2% in the low-dose gradient regions and within 1-mm in the high-dose gradient regions. Off-axis doses agreed to within 2 mm in the high-dose gradient regions and within 3% in the low-dose gradient regions. Arced-beam calculations of dose output and depth dose at mid arc were evaluated by comparing to data measured using two cylindrical water phantoms with radii of 12 and 15 cm at 10 and 15 MeV. Dose output was measured for all combinations of phantom radii of curvature, collimator widths (4, 5, and 6 cm), and arc angles (0 degrees, 20 degrees, 40 degrees, 60 degrees, 80 degrees, and 90 degrees) for both beam energies. Results showed the calculated mid-arc dose output to agree within 2% of measurement for all combinations. For a 90 degree arc angle and 5 x 20 cm2 field size, the calculated mid-arc depth dose in the low-dose gradient region agreed to within 2% of measurement for all depths at 10 MeV and for depths greater than depth of dose maximum R100 at 15 MeV. For depths in the buildup region at 15 MeV the calculations overestimated the measured dose by as much as 3.4%. Mid-arc depth dose in the high-dose gradient region agreed to within 2.2 mm of measured dose. Calculated two-dimensional relative dose distributions in the plane of rotation were compared to dose measurements using film in a cylindrical polystyrene phantom for a 90 degree arc angle and field widths of 4, 5, and 6 cm at 10 and 15 MeV. Results showed that off-axis dose at the ends of arc (without skin collimation) agreed to within 2% in the low-dose gradient region and to within 1.2 mm in the high-dose gradient region. This work showed that the accuracy of the PBRA arced-beam dose model met the criteria specified by Van Dyk et al. [Int. J. Radiat. Oncol. Biol. Phys. 26, 261-273 (1993)] with the exception of the buildup region of the 15 MeV beam. Based on the present results, results of a previous study showing acceptable accuracy in the presence of skin collimation, and results of a previous study showing acceptable accuracy in the presence of internal heterogeneities, it is concluded that the PBRA arced-beam dose model should be adequate for planning electron arc therapy.

Comprehensive evaluation of a commercial macro Monte Carlo electron dose calculation implementation using a standard verification data set.

A commercial electron dose calculation software implementation based on the macro Monte Carlo algorithm has recently been introduced. We have evaluated the performance of the system using a standard verification data set comprised of two-dimensional (2D) dose distributions in the transverse plane of a 15 X 15 cm2 field. The standard data set was comprised of measurements performed for combinations of 9-MeV and 20-MeV beam energies and five phantom geometries. The phantom geometries included bone and air heterogeneities, and irregular surface contours. The standard verification data included a subset of the data needed to commission the dose calculation. Additional required data were obtained from a dosimetrically equivalent machine. In addition, we performed 2D dose measurements in a water phantom for the standard field sizes, a 4 cm X 4 cm field, a 3 cm diameter circle, and a 5 cm X 13 cm triangle for the 6-, 9-, 12-, 15-, and 18-MeV energies of a Clinac 21EX. Output factors were also measured. Synthetic CT images and structure contours duplicating the measurement configurations were generated and transferred to the treatment planning system. Calculations for the standard verification data set were performed over the range of each of the algorithm parameters: statistical precision, grid-spacing, and smoothing. Dose difference and distance-to-agreement were computed for the calculation points. We found that the best results were obtained for the highest statistical precision, for the smallest grid spacing, and for smoothed dose distributions. Calculations for the 21EX data were performed using parameters that the evaluation of the standard verification data suggested would produce clinically acceptable results. The dose difference and distance-to-agreement were similar to that observed for the standard verification data set except for the portion of the triangle field narrower than 3 cm for the 6- and 9-MeV electron beams. The output agreed with measurements to within 2%, with the exception of the 3-cm diameter circle and the triangle for 6 MeV, which were within 5%. We conclude that clinically acceptable results may be obtained using a grid spacing that is no larger than approximately one-tenth of the distal falloff distance of the electron depth dose curve (depth from 80% to 20% of the maximum dose) and small relative to the size of heterogeneities. For judicious choices of parameters, dose calculations agree with measurements to better than 3% dose difference and 3-mm distance-to-agreement for fields with dimensions no less than about 3 cm.

Displacement-based binning of time-dependent computed tomography image data sets.

Respiration can cause tumors in the thorax or abdomen to move by as much as 3 cm; this movement can adversely affect the planning and delivery of radiation treatment. Several techniques have been used to compensate for respiratory motion, but all have shortcomings. Manufacturers of computed tomography (CT) equipment have recently used a technique developed for cardiac CT imaging to track respiratory-induced anatomical motion and to sort images according to the phase of the respiratory cycle they represent. Here we propose a method of generating CT images that accounts for respiratory-induced anatomical motion on the basis of displacement, i.e., displacement-binned CT image sets. This technique has shown great promise, however, it is not fully supported by currently used CT image reconstruction software. As an interim solution, we have developed a method for extracting displacement-binned CT image data sets from data sets assembled on the basis of a prospectively determined breathing phase acquired on a multislice helical CT scanner. First, the projection data set acquired from the CT scanner was binned at small phase intervals before reconstruction. The manufacturer's software then generated image sets identified as belonging to particular phases of the respiratory cycle. All images were then individually correlated to the displacement of an external fiducial marker. Next, CT image data sets were resorted on the basis of the displacement and assigned an appropriate phase. Finally, displacement-binned image data sets were transferred to a treatment-planning system for analysis. Although the technique is currently limited by the phase intervals allowed by the CT software, some improvement in image reconstruction was seen, indicating that this technique is useful at least as an interim measure.

Modeling skin collimation using the electron pencil beam redefinition algorithm.

Skin collimation is an important tool for electron beam therapy that is used to minimize the penumbra when treating near critical structures, at extended treatment distances, with bolus, or using arc therapy. It is usually made of lead or lead alloy material that conforms to and is placed on patient surface. Presently, commercially available treatment-planning systems lack the ability to model skin collimation and to accurately calculate dose in its presence. The purpose of the present work was to evaluate the use of the pencil beam redefinition algorithm (PBRA) in calculating dose in the presence of skin collimation. Skin collimation was incorporated into the PBRA by terminating the transport of electrons once they enter the skin collimator. Both fixed- and arced-beam dose calculations for arced-beam geometries were evaluated by comparing them with measured dose distributions for 10- and 15-MeV beams. Fixed-beam dose distributions were measured in water at 88-cm source-to-surface distance with an air gap of 32 cm. The 6 x 20-cm2 field (dimensions projected to isocenter) had a 10-mm thick lead collimator placed on the surface of the water with its edge 5 cm inside the field's edge located at +10 cm. Arced-beam dose distributions were measured in a 13.5-cm radius polystyrene circular phantom. The beam was arced 90 degrees (-45 degrees to +45 degrees), and 10-mm thick lead collimation was placed at +/- 30 degrees. For the fixed beam at 10 MeV, the PBRA- calculated dose agreed with measured dose to within 2.0-mm distance to agreement (DTA) in the regions of high-dose gradient and 2.0% in regions of low dose gradient. At 15 MeV, the PBRA agreed to within a 2.0-mm DTA in the regions of high-dose gradient; however, the PBRA underestimated the dose by as much as 5.3% over small regions at depths less than 2 cm because it did not model electrons scattered from the edge of the skin collimation. For arced beams at 10 MeV, the agreement was 1-mm DTA in the high-dose gradient regions, and 2% in the low-dose gradient regions. For arced beams at 15 MeV, the agreement was 1 mm in the high-dose gradient regions, and in the low-dose gradient region at depth less than 2 cm, as much as 5% dose difference was observed. This study demonstrated the ease with which skin collimation can be incorporated into the PBRA. The good agreement of PBRA calculated with measured dose shows that the PBRA is likely sufficiently accurate for clinical use in the presence of skin collimation for electron arc therapy. To further improve the accuracy of the PBRA in regions having significant electrons scattered from the edge of the skin collimation would require transporting the electrons through the lead skin collimation near its edges.