Validation of a precision radiochromic film dosimetry system for quantitative two-dimensional imaging of acute exposure dose distributions

We present an evaluation of the precision and accuracy of image-based radiochromic film (RCF) dosimetry performed using a commercial RCF product (Gafchromic MD-55-2, Nuclear Associates, Inc.) and a commercial high-spatial resolution (100 microm pixel size) He-Ne scanning-laser film-digitizer (Personal Densitometer, Molecular Dynamics, Inc.) as an optical density (OD) imaging system. The precision and accuracy of this dosimetry system are evaluated by performing RCF imaging dosimetry in well characterized conformal external beam and brachytherapy high dose-rate (HDR) radiation fields. Benchmarking of image-based RCF dosimetry is necessary due to many potential errors inherent to RCF dosimetry including: a temperature-dependent time evolution of RCF dose response; nonuniform response of RCF; and optical-polarization artifacts. In addition, laser-densitometer imaging artifacts can produce systematic OD measurement errors as large as 35% in the presence of high OD gradients. We present a RCF exposure and readout protocol that was developed for the accurate dosimetry of high dose rate (HDR) radiation sources. This protocol follows and expands upon the guidelines set forth by the American Association of Physicists in Medicine (AAPM) Task Group 55 report. Particular attention is focused on the OD imaging system, a scanning-laser film digitizer, modified to eliminate OD artifacts that were not addressed in the AAPM Task Group 55 report. RCF precision using this technique was evaluated with films given uniform 6 MV x-ray doses between 1 and 200 Gy. RCF absolute dose accuracy using this technique was evaluated by comparing RCF measurements to small volume ionization chamber measurements for conformal external-beam sources and an experimentally validated Monte Carlo photon-transport simulation code for a 192Ir brachytherapy source. Pixel-to-pixel standard deviations of uniformly irradiated films were less than 1% for doses between 10 and 150 Gy; between 1% and 5% for lower doses down to 1 Gy and 1% and 1.5% for higher doses up to 200 Gy. Pixel averaging to form 200-800 microm pixels reduces these standard deviations by a factor of 2 to 5. Comparisons of absolute dose show agreement within 1.5%-4% of dose benchmarks, consistent with a highly accurate dosimeter limited by its observed precision and the precision of the dose standards to which it is compared. These results provide a comprehensive benchmarking of RCF, enabling its use in the commissioning of novel HDR therapy sources.     

Variation of sensitometric curves of radiographic films in high energy photon beams.

Film dosimetry is an important tool for the verification of irradiation techniques. The shape of the sensitometric curve depends on the type of film as well as on the irradiation and processing conditions. Existing data concerning the influence of irradiation geometry on the sensitometric curve are conflicting. In particular the variation of optical density, OD, with field size and depth in a phantom shows large differences in magnitude between various authors. This variation, as well as the effect of beam energy and film plane orientation on OD, was therefore investigated for two types of film, Kodak X-Omat V and Agfa Structurix D2. Films were positioned in a solid phantom, either perpendicular or (almost) parallel to the beam axis, and irradiated to different dose levels using various photon beams (Co-60, 6 MV, 15 MV, 18 MV, 45 MV). It was found that the sensitometric curves of the Kodak film derived at different depths are almost identical for the four x-ray beams. For the Kodak film the differences in OD with depth are less than 2%, except for the Co-60 beam, where the difference is about 4% at 10 cm depth for a 15 cm x 15 cm field. The slope of the sensitometric curve of the Agfa film is somewhat more dependent on photon beam energy, depth and field size. The sensitometric curves of both types of film are almost independent of the film plane orientation, except for shallow depths. For Co-60 and for the same dose, the Kodak and Agfa films gave at dose maximum an OD lower by 4% and 6%, respectively, for the parallel compared to the perpendicular geometry. Good dosimetric results can be obtained if films from the same batch are irradiated with small to moderate field sizes (up to about 15 cm x 15 cm), at moderate depths (up to about 15 cm), using a single calibration curve, e.g., for a 10 cm x 10 cm field.     

Optically stimulated luminescence (OSL) of carbon-doped aluminum oxide (Al2O3:C) for film dosimetry in radiotherapy

INTRODUCTION AND PURPOSE: Conventional x-ray films and radiochromic films have inherent challenges for high precision radiotherapy dosimetry. Here we have investigated basic characteristics of optically stimulated luminescence (OSL) of irradiated films containing carbon-doped aluminum oxide (Al2O3:C) for dosimetry in therapeutic photon and electron beams. MATERIALS AND METHODS: The OSL films consist of a polystyrene sheet, with a top layer of a mixture of single crystals of Al2O3:C, ground into a powder, and a polyester base. The total thickness of the films is 0.3 mm. Measurements have been performed in a water equivalent phantom, using 4, 6, 10, and 18 MV photon beams, and 6-22 MeV electron beams. The studies include assessment of the film response (acquired OSL signal/delivered dose) on delivered dose (linearity), dose rate (1-6 Gy/min), beam quality, field size and depth (6 MV, ranges 4 x 4-30 x 30 cm2, dmax-35 cm). Doses have been derived from ionization chamber measurements. OSL films have also been compared with conventional x-ray and GafChromic films for dosimetry outside the high dose area, with a high proportion of low dose scattered photons. In total, 787 OSL films have been irradiated. RESULTS: Overall, the OSL response for electron beams was 3.6% lower than for photon beams. Differences between the various electron beam energies were not significant. The 6 and 18 MV photon beams differed in response by 4%. No response dependencies on dose rate were observed. For the 6 MV beam, the field size and depth dependencies of the OSL response were within +/-2.5%. The observed inter-film response variation for films irradiated with the same dose varied from 1% to 3.2% (1 SD), depending on the measurement day. At a depth of 20 cm, 5 cm outside the 20 x 20 cm2 6 and 18 MV beams, an over response of 17% was observed. In contrast to GafChromic and conventional x-ray films, the response of the Al2O3:C films is linear in the clinically relevant dose range 0-200 cGy. CONCLUSIONS: Measurement of the OSL signal of irradiated films containing Al2O3:C is a promising technique for film dosimetry in radiotherapy with no or small response variations with dose rate, beam quality, field size and depth, and a linear response from 0 to 200 cGy.     

The impact of MLC transmitted radiation on EPID dosimetry for dynamic MLC beams

The purpose of this study was to experimentally quantify the change in response of an amorphous silicon (a-Si) electronic portal imaging device (EPID) to dynamic multileaf collimator (dMLC) beams with varying MLC-transmitted dose components and incorporate the response into a commercial treatment planning system (TPS) EPID prediction model. A combination of uniform intensity dMLC beams and static beams were designed to quantify the effect of MLC transmission on EPID response at the central axis of 10 x 10 cm2 beams, at off-axis positions using wide dMLC beam profiles, and at different field sizes. The EPID response to MLC transmitted radiation was 0.79 +/- 0.02 of the response to open beam radiation at the central axis of a 10 x 10 cm2 field. The EPID response to MLC transmitted radiation was further reduced relative to the open beam response with off-axis distance. The EPID response was more sensitive to field size changes for MLC transmitted radiation compared to open beam radiation by a factor of up to 1.17 at large field sizes. The results were used to create EPID response correction factors as a function of the fraction of MLC transmitted radiation, off-axis distance, and field size. Software was developed to apply the correction factors to each pixel in the TPS predicted EPID image. The corrected images agreed more closely with the measured EPID images in areas of intensity modulated fields with a large fraction of MLC transmission and, as a result the accuracy of portal dosimetry with a-Si EPIDs can be improved. Further investigation into the detector response function and the radiation source model are required to achieve improvements in accuracy for the general case.     

Evaluation of Kodak EDR2 film for dose verification of intensity modulated radiation therapy delivered by a static multileaf collimator

A new type of radiographic film, Kodak EDR2 film, was evaluated for dose verification of intensity modulated radiation therapy (IMRT) delivered by a static multileaf collimator (SMLC). A sensitometric curve of EDR2 film irradiated by a 6 MV x-ray beam was compared with that of Kodak X-OMAT V (XV) film. The effects of field size, depth and dose rate on the sensitometric curve were also studied. It is found that EDR2 film is much less sensitive than XV film. In high-energy x-ray beams, the double hit process is the dominant mechanism that renders the grains on EDR2 films developable. As a result, in the dose range that is commonly used for film dosimetry for IMRT and conventional external beam therapy, the sensitometric curves of EDR2 films cannot be approximated as a linear function, OD = c * D. Within experimental uncertainty, the film sensitivity does not depend on the dose rate (50 vs 300 MU/min) or dose per pulse (from 1.0 x 10(-4) to 4.21 x 10(-4) Gy/pulse). Field sizes and depths (up to field size of 10 x 10 cm2 and depth = 10 cm) have little effect on the sensitometric curves. Percent depth doses (PDDs) for both 6 and 23 MV x rays were measured with both EDR2 and XV films and compared with ion chamber data. Film data are within 2.5% of the ion chamber results. Dose profiles measured with EDR2 film are consistent with those measured with an ion chamber. Examples of measured IMRT isodose distributions versus calculated isodoses are presented. We have used EDR2 films for verification of all IMRT patients treated by SMLC in our clinic. In most cases, with EDR2 film, actual clinical daily fraction doses can be used for verification of composite isodose distributions of SMLC-based IMRT.     

Dosimetry limitations and a dose correction methodology for step-and-shoot IMRT

For the step-and-shoot intensity-modulated radiation therapy (IMRT) technique, the combination of high dose rate, multiple beam segments and low dose per segment can lead to significant differences between the planned dose and the dose delivered to the patient. In this technique, a dose delivery inaccuracy known as the 'overshoot' effect is caused by the dose servo control system. This typically occurs in the first and last beam segments and causes an over- and underdose, respectively. Some dose positional error in the segment sequence is also possible there. Commercial ionization chambers (RK-type) and radiographic Kodak films were used for the measurements. The reported results were obtained using the Pinnacle(3)-V6.2 treatment planning system and a Varian Clinac 21 EX linear accelerator equipped with a 120-leaf Millennium MLC. The dose inaccuracy measurements were based on the comparison of the dose and profiles for reference fields and fields irradiated with the step-and-shoot technique. For our linear accelerators, an 'overshoot' effect ranging from 0.1 to 0.6 MU was found, corresponding to a dose rate from 100 to 600 MU min(-1), respectively. For segments with off-axis distances from 0 to 5.5 cm with >3.5 MU per segment and all dose rates, a MLC leaf-position error of <1 mm was measured. For segments with an off-axis distance of 9.5 cm, a positional error >2 mm was measured for 600 MU min(-1) and 1 MU per segment. The purpose of this study was to find a correction method for segmental dose errors caused by the 'overshoot' effect when small monitor unit and high dose rate are used. To better represent the fluctuation of the segment doses in the beam, a dose ratio between reference and step-and-shoot irradiated fields was defined. A method for the correction of segment dose inaccuracies and a quality assurance programme for the 'overshoot' effect were developed. The ordering of the biggest segment shape in the segment sequence was studied for ten randomly selected prostate patients planned for IMRT. The results of this work can be used to improve the agreement between the planned and delivered doses for IMRT.     

Analysis of various beamlet sizes for IMRT with 6 MV photons

Application of intensity modulated radiation therapy (IMRT) using multileaf collimation often requires the use of small beamlets to optimize the delivered radiation distribution. Small-beam dose distribution measurements were compared to dose distributions calculated using a commercial treatment planning system that models its data acquired using measurements from relatively large fields. We wanted to evaluate only the penumbra, percent depth-dose (PDD) and output model, so we avoided dose distribution features caused by rounded leaf ends and interleaf leakage by making measurements using the secondary collimators. We used a validated radiochromic film dosimetry system to measure high-resolution dose distributions of 6 MV photon beams. A commercial treatment planning system using the finite size pencil beam (FSPB) dose calculation algorithm was commissioned using measured central axis outputs from 4.0x4.0 to 40.0x40.0 cm2 beams and radiographic-film profile measurements of a 4.0x4.0 cm2 beam at twice the depth of maximum dose (dmax). Calculated dose distributions for square fields of 0.5x0.5 cm2, and 1.0x1.0 cm2, to 6.0x6.0 cm2, in 1.0x1.0 cm2, increments were compared against radiochromic film measurements taken with the film oriented parallel to the beam central axis in a water equivalent phantom. The PDD of the smaller field sizes exhibited behavior typical of small fields, namely a decrease in dmax with decreasing field size. The FSPB accurately modeled the depth-dose and central axis output for depths deeper than the nominal dmax of 1.5 cm plus 0.5 cm. The dose distribution in the build-up and penumbra regions was not accurately modeled for depths less than 2 cm, especially for the fields of 2.0x2.0 cm2 and smaller. Using the gamma function with 2 mm and 2% criteria, the dose model was shown to accurately predict the penumbra. While for single small beams the compared dose distributions passed the gamma function criteria, the clinical appropriateness of these criteria is not clear for a composite IMRT plan. Further investigation of the cumulative impact of the observed dose discrepancies is warranted. We speculate that the observed differences in the penumbra regions arise from some energy dependent artifact in the radiographic-film profiles used for commissioning. In the future, radiochromic film based commissioning might provide a more accurate data set for dose modeling.     

Ionization chamber dosimetry of small photon fields: a Monte Carlo study on stopping-power ratios for radiosurgery and IMRT beams

Absolute dosimetry with ionization chambers of the narrow photon fields used in stereotactic techniques and IMRT beamlets is constrained by lack of electron equilibrium in the radiation field. It is questionable that stopping-power ratio in dosimetry protocols, obtained for broad photon beams and quasi-electron equilibrium conditions, can be used in the dosimetry of narrow fields while keeping the uncertainty at the same level as for the broad beams used in accelerator calibrations. Monte Carlo simulations have been performed for two 6 MV clinical accelerators (Elekta SL-18 and Siemens Mevatron Primus), equipped with radiosurgery applicators and MLC. Narrow circular and Z-shaped on-axis and off-axis fields, as well as broad IMRT configured beams, have been simulated together with reference 10 x 10 cm2 beams. Phase-space data have been used to generate 3D dose distributions which have been compared satisfactorily with experimental profiles (ion chamber, diodes and film). Photon and electron spectra at various depths in water have been calculated, followed by Spencer-Attix (delta = 10 keV) stopping-power ratio calculations which have been compared to those used in the IAEA TRS-398 code of practice. For water/air and PMMA/air stopping-power ratios, agreements within 0.1% have been obtained for the 10 x 10 cm2 fields. For radiosurgery applicators and narrow MLC beams, the calculated s(w,air) values agree with the reference within +/-0.3%, well within the estimated standard uncertainty of the reference stopping-power ratios (0.5%). Ionization chamber dosimetry of narrow beams at the photon qualities used in this work (6 MV) can therefore be based on stopping-power ratios data in dosimetry protocols. For a modulated 6 MV broad beam used in clinical IMRT, s(w,air) agrees within 0.1% with the value for 10 x 10 cm2, confirming that at low energies IMRT absolute dosimetry can also be based on data for open reference fields. At higher energies (24 MV) the difference in s(w,air) was up to 1.1%, indicating that the use of protocol data for narrow beams in such cases is less accurate than at low energies, and detailed calculations of the dosimetry parameters involved should be performed if similar accuracy to that of 6 MV is sought.     

Experimental investigation of the response of an amorphous silicon EPID to intensity modulated radiotherapy beams

The aim of this work was to experimentally determine the difference in response of an amorphous silicon (a-Si) electronic portal imaging device (EPID) to the open and multileaf collimator (MLC) transmitted beam components of intensity modulated radiation therapy (IMRT) beams. EPID dose response curves were measured for open and MLC transmitted (MLCtr) 10 x 10 cm2 beams at central axis and with off axis distance using a shifting field technique. The EPID signal was obtained by replacing the flood-field correction with a pixel sensitivity variation matrix correction. This signal, which includes energy-dependent response, was then compared to ion-chamber measurements. An EPID calibration method to remove the effect of beam energy variations on EPID response was developed for IMRT beams. This method uses the component of open and MLCtr fluence to an EPID pixel calculated from the MLC delivery file and applies separate radially dependent calibration factors for each component. The calibration procedure does not correct for scatter differences between ion chamber in water measurements and EPID response; these must be accounted for separately with a kernel-based approach or similar method. The EPID response at central axis for the open beam was found to be 1.28 +/- 0.03 of the response for the MLCtr beam, with the ratio increasing to 1.39 at 12.5 cm off axis. The EPID response to MLCtr radiation did not change with off-axis distance. Filtering the beam with copper plates to reduce the beam energy difference between open and MLCtr beams was investigated; however, these were not effective at reducing EPID response differences. The change in EPID response for uniform sliding window IMRT beams with MLCtr dose components from 0.3% to 69% was predicted to within 2.3% using the separate EPID response calibration factors for each dose component. A clinical IMRT image calibrated with this method differed by nearly 30% in high MLCtr regions from an image calibrated with an open beam calibration factor only. Accounting for the difference in EPID response to open and MLCtr radiation should improve IMRT dosimetry with a-Si EPIDs.     

An application of GafChromic MD-55 film for 67.5 MeV clinical proton beam dosimetry

The purpose of this study is to explore the use of GafChromic MD-55 (RC) film for 67.5 MeV clinical proton beam dosimetry at the Crocker Nuclear Laboratory, University of California, Davis. Several strips of RC film 6 cm x 6 cm in dimension were irradiated at a depth of 18.2 mm corresponding to the middle of a 24 mm spread-out Bragg peak (SOBP). The films were irradiated to a proton dose in the range of 0.5 Gy to 100 Gy. The beam profiles were also measured at the middle of the 24 mm SOBP. The Bragg peak was measured by using a wedge shaped phantom made of Lucite. The Bragg peak measured with RC film was compared with diode and ionization chamber measurements. After background subtraction, the calibration of the dose response of RC film showed, to a maximum deviation of 10%, a linear increase of optical density (OD) with dose from 0.5 to 100 Gy. The uniformity of OD over a single sheet of film showed a variation of +/-6%. The distal-fall off between 90% and 20% measured with GafChromic film for the Bragg peak was 1.3 mm as compared to 1.1 mm for a diode measurement and 1.4 mm for an ionization chamber measurement. The FWHM of the Bragg peak was 7.5 mm when measured with GafChromic film, 5.3 mm when measured with a diode and 8.1 mm as measured by an ionization chamber. The peak/plateau ratio with GafChromic film was 3.3 as compared to 3.7 with a diode and 3.2 with an ionization chamber. In conclusion, GafChromic MD-55 film may be a useful and convenient detector for dose measurement and quality assurance programmes of proton beams.     

Matchline dosimetry in step and shoot IMRT fields: a film study

The Varian millennium 120 multileaf collimator has curved leaf ends. Transmission through the leaf ends generates a small asymmetric penumbral dose effect. This design can lead to hot spots between neighbouring beam segments during step and shoot IMRT dose delivery. We have observed some matchlines with film for clinical beams optimized using the pinnacle radiotherapy treatment planning system; hence we sought to verify the optimum leaf offset required to minimize the matchline effect. An in-house program was created to control the MLC leaf banks in 2 cm steps with a 2 cm gap. The gap was varied by the following offset values from 0.0 to 0.1 cm. Two types of radiographic films (Kodak EDR and XV films) and a radiochromic film (Gafchromic MD-55-2) were used to measure the optical density maps. The films were positioned in a solid water phantom perpendicular to the beam axis and irradiated at d(max) using a 6 MV photon beam. An ion chamber (IC4) was used to measure point doses for normalization in a beam umbral minima position. The relative mean peak to valley dose ratios measured with no leaf offset were 1.31, 1.30 and 1.31 for the XV, EDR2 and Gafchromic films, respectively. For a 0.07 cm gap per leaf and a performance of end leaf repeatability of 0.01 cm, the central matchline was reduced to about 1.0 for all dosimeters, with two mini-peaks measured as 1.05, 1.05 and 1.08 each side of the matchline, for XV, EDR2 and Gafchromic, respectively. The average relative dose across the umbra for this offset was XO-mat V = 1.01, EDR = 1.01 and radiochromic film = 1.02, respectively. While we expected the beam penumbral tails from segment neighbours to cause overprediction of the dose in the central valley regions due to the energy response of radiographic films, by normalizing all dosimeters to an ion chamber reading in the minimum we could not observe any major shape distortion between the radiographic film and radiochromic film results. In conclusion, relative doses measured by radiographic and radiochromic films agree well with IC4 within +/-2%.     

Underdosage of the upper-airway mucosa for small fields as used in intensity-modulated radiation therapy: a comparison between radiochromic film measurements,Monte Carlo simulations,and collapsed cone convolution calculations

Head-and-neck tumors are often situated at an air-tissue interface what may result in an underdosage of part of the tumor in radiotherapy treatments using megavoltage photons, especially for small fields. In addition to effects of transient electronic disequilibrium, for these small fields, an increased lateral electron range in air will result in an important extra reduction of the central axis dose beyond the cavity. Therefore dose calculation algorithms need to model electron transport accurately. We simulated the trachea by a 2 cm diameter cylindrical air cavity with the rim situated 2 cm beneath the phantom surface. A 6 MV photon beam from an Elekta SLiplus linear accelerator, equipped with the standard multileaf collimator (MLC), was assessed. A 10 x 2 cm2 and a 10 x 1 cm2 field, both widthwise collimated by the MLC, were applied with their long side parallel to the cylinder axis. Central axis dose rebuild-up was studied. Radiochromic film measurements were performed in an in-house manufactured polystyrene phantom with the films oriented either along or perpendicular to the beam axis. Monte Carlo simulations were performed with BEAM and EGSnrc. Calculations were also performed using the pencil beam (PB) algorithm and the collapsed cone convolution (CCC) algorithm of Helax-TMS (MDS Nordion, Kanata, Cahada) version 6.0.2 and using the CCC algorithm of Pinnacle (ADAC Laboratories, Milpitas, CA, USA) version 4.2. A very good agreement between the film measurements and the Monte Carlo simulations was found. The CCC algorithms were not able to predict the interface dose accurately when lateral electronic disequilibrium occurs, but were shown to be a considerable improvement compared to the PB algorithm. The CCC algorithms overestimate the dose in the rebuild-up region. The interface dose was overestimated by a maximum of 31% or 54%, depending on the implementation of the CCC algorithm. At a depth of 1 mm, the maximum dose overestimation was 14% or 24%.     

Normalized sensitometric curves for the verification of hybrid IMRT treatment plans with multiple energies

With the clinical implementation of time-variable dose patterns and intensity modulated radiotherapy (IMRT) film dosimetry has regained popularity. Films are currently the most frequently used dosimetric means for patient specific quality assurance in IMRT. A common method is to verify a so-called hybrid IMRT plan, which is the patient specific treatment plan with unmodified fluence patterns recalculated in a dedicated phantom. For such applications the sensitometric curve, i.e., the relation between optical density (OD) and absorbed dose, should not depend critically on beam energy, field size and depth, or film orientation. In order to minimize the influence of all these variables a normalization of sensitometric curves is performed at various photon beam energies (6 MV, 10 MV, 25 MV). By doing so one unique sensitometric curve can be used for these three beam qualities. This holds for both film types investigated: Kodak X-Omat V films and EDR-2 films. Additionally, the influence of field size, depth and film orientation on a normalized sensitometric curve is determined for both film types. For doses smaller than 0.8 Gy for X-Omat V and doses smaller than 3 Gy for EDR-2 films the field size variation of normalized sensitometric curves is much smaller than 3% for fields up to 20 x 20 cm2. For X-Omat V films all differences between sensitometric curves determined at depths of 5, 10, and 15 cm are smaller than 3%. For EDR-2 films deviations larger than 3% are only observed at low net OD smaller than 0.25. The dependence of film orientation (parallel versus perpendicular) on a normalized sensitometric curve is found to be not critical. However, processing conditions have the largest influence and can result in differences up to 20% for sensitometric curves derived from films of the same batch but using different film processors. When normalizing sensitometric curves to the dose value necessary to obtain a net OD=1 for that respective geometry and energy the large energy dependence of sensitometric curves can be almost eliminated. This becomes especially important for the verification of hybrid IMRT plans with multiple energies. Additionally, such a normalization minimizes other influences such as field size, depth, and film orientation. This method is generally applicable to both Kodak X-Omat V and EDR-2 films. In order to achieve the highest accuracy level an upper dose limit of 0.8 Gy for X-Omat V films and 3 Gy for EDR-2 films should be taken into account. However, these dose limits may vary with film reading instrument and film processor.     

Calibration of EBT2 film by the PDD method with scanner non-uniformity correction

The EBT2 film together with a flatbed scanner is a convenient dosimetry QA tool for verification of clinical radiotherapy treatments. However, it suffers from a relatively high degree of uncertainty and a tedious film calibration process for every new lot of films, including cutting the films into several small pieces, exposing with different doses, restoring them back and selecting the?proper region of interest (ROI) for each piece for curve fitting. In this work, we present a percentage depth dose (PDD) method that can accurately calibrate the EBT2 film together with the scanner non-uniformity correction and provide an easy way to perform film dosimetry. All films were scanned before and after the irradiation in one of the two homemade 2?mm thick acrylic frames (one portrait and the other landscape), which was located at a fixed position on the scan bed of an Epson 10?000XL scanner. After the pre-irradiated scan, the film was placed parallel to the beam central axis and sandwiched between six polystyrene plates (5?cm thick each), followed by irradiation of a 20?×?20?cm(2)?6 MV photon beam. Two different beams on times were used on two different films to deliver a dose to the film ranging from 32 to 320 cGy. After the post-irradiated scan, the net optical densities for a total of 235 points on the beam central axis on the films were auto-extracted and compared with the corresponding depth doses that were calculated through the measurement of a 0.6 cc farmer chamber and the related PDD table to perform the curve fitting. The portrait film location was selected for routine calibration, since the central beam axis on the film is parallel to the scanning direction, where non-uniformity correction is not needed (Ferreira et al 2009 Phys. Med. Biol. 54 1073-85). To perform the scanner non-uniformity calibration, the cross-beam profiles of the film were analysed by referencing the measured profiles from a Profiler?. Finally, to verify our method, the films were exposed to 60° physical wedge fields and the compositive fields, and their relative dose profiles were compared with those from the water phantom measurement. The fitting uncertainty was less than 0.5% due to the many calibration points, and the overall calibration uncertainty was within 3% for doses above 50 cGy, when the average of four films were used for the calibration. According to our study, the non-uniformity calibration factor was found to be independent of the given dose for the EBT2 film and the relative dose differences between the profiles measured by the film and the Profiler were within 1.5% after applying the non-uniformity correction. For the verification tests, the relative dose differences between the measurements by films and in the water phantom, when the average of three films were used, were generally within 3% for the 60° wedge fields and compositive fields, respectively. In conclusion, our method is convenient, time-saving and cost-effective, since no film cutting is needed and only two films with two exposures are needed.     

Film dosimetry of small elongated electron beams for treatment planning

The characteristics of 5, 7, 10, 12, 15, and 18 Mev electron beams for small elongated fields of dimensions L X W (where L = 1, 2, 3, 4, 5, and 10 cm; and W = 1, 2, 3, 4, 5, and 10 cm) have been studied. Film dosimetry and parallel-plate ion chamber measurements have been used to obtain various dose parameters. Selective results of a series of systematic measurements for central axis depth dose data, uniformity index, field flatness, and relative output factors of small elongated electron beams are reported. The square-root method is employed to predict the beam data of small elongated electron fields from corresponding small square electron fields using film dosimetry. The single parameter area/perimeter radio A/P is used to characterize the relative output factors of elongated electron beams. It is our conclusion that for clinical treatment planning square-root method may be applied with caution in determining the beam characteristics of small elongated electron fields from film dosimetry. The calculated and estimated relative output factors from square-root method and A/P ratio are in good agreement and show agreement to within 1% with the measured film values.     

Monte Carlo simulation of a dynamic MLC based on a multiple source model.

Detailed knowledge of the characteristics of the radiation field shaped by a multileaf collimator (MLC) is essential in intensity modulated radiotherapy (IMRT). A previously developed multiple source model (MSM) for a 6 MV beam was extended to a 15 MV beam and supplemented with an accurate model of an 80-leaf dynamic MLC. Using the supplemented MSM and the MC code GEANT, lateral dose distributions were calculated in a water phantom and a portal water phantom. A field which is normally used for the validation of the step and shoot technique and a field from a realistic IMRT treatment plan delivered with dynamic MLC are investigated. To assess possible spectral changes caused by the modulation of beam intensity by an MLC, the energy spectra in five portal planes were calculated for moving slits of different widths. The extension of the MSM to 15 MV was validated by analysing energy fluences, depth doses and dose profiles. In addition, the MC-calculated primary energy spectrum was verified with an energy spectrum which was reconstructed from transmission measurements. MC-calculated dose profiles using the MSM for the step and shoot case and for the dynamic MLC case are in very good agreement with the measured data from film dosimetry. The investigation of a 13 cm wide field shows an increase in mean photon energy of up to 16% for the 0.25 cm slit compared to the open beam for 6 MV and of up to 6% for 15 MV, respectively. In conclusion, the MSM supplemented with the dynamic MLC has proven to be a powerful tool for investigational and benchmarking purposes or even for dose calculations in IMRT.     

Verification of total body photon irradiation dosimetry techniques

A method of verifying the dosimetry of patients undergoing total body irradiation (TBI) with photon beams having energies from cobalt-60 to 25 MV is presented. A simple set of spot checks at the TBI axis has been used to verify data used for TBI dosimetry. Calculations to verify dose delivered to TBI patients are done in the same manner as those irradiated at standard treatment distances. A simple method of effective field size determination for various anatomical locations in a typical adult is presented. Measurements in an Alderson phantom with thermoluminescent dosimeters and an ion chamber at several anatomical locations indicate that this calculational method can predict the dose along the patient axis to within 4% for 60Co and 18-MV photon beams, provided the dosimetry data are appropriate (as determined by the spot checks). Results of intercomparisons of TBI beam calibration, off-axis and depth-dose data at various institutions visited by the Radiological Physics Center are also presented.     

An independent dose calculation algorithm for MLC-based stereotactic radiotherapy

We have developed an algorithm to calculate dose in a homogeneous phantom for radiotherapy fields defined by multi-leaf collimator (MLC) for both static and dynamic MLC delivery. The algorithm was developed to supplement the dose algorithms of the commercial treatment planning systems (TPS). The motivation for this work is to provide an independent dose calculation primarily for quality assurance (QA) and secondarily for the development of static MLC field based inverse planning. The dose calculation utilizes a pencil-beam kernel. However, an explicit analytical integration results in a closed form for rectangular-shaped beamlets, defined by single leaf pairs. This approach reduces spatial integration to summation, and leads to a simple method of determination of model parameters. The total dose for any static or dynamic MLC field is obtained by summing over all individual rectangles from each segment which offers faster speed to calculate two-dimensional dose distributions at any depth in the phantom. Standard beam data used in the commissioning of the TPS was used as input data for the algorithm. The calculated results were compared with the TPS and measurements for static and dynamic MLC. The agreement was very good (<2.5%) for all tested cases except for very small static MLC sizes of 0.6 cm x 0.6 cm (<6%) and some ion chamber measurements in a high gradient region (<4.4%). This finding enables us to use the algorithm for routine QA as well as for research developments.     

Comparison of inhomogeneity correction algorithms in small photon fields

Algorithms such as convolution superposition, Batho, and equivalent pathlength which were originally developed and validated for conventional treatments under conditions of electronic equilibrium using relatively large fields greater than 5 x 5 cm2 are routinely employed for inhomogeneity corrections. Modern day treatments using intensity modulated radiation therapy employ small beamlets characterized by the resolution of the multileaf collimator. These beamlets, in general, do not provide electronic equilibrium even in a homogeneous medium, and these effects are exaggerated in media with inhomogenieties. Monte Carlo simulations are becoming a tool of choice in understanding the dosimetry of small photon fields as they encounter low density media. In this study, depth dose data from the Monte Carlo simulations are compared to the results of the convolution superposition, Batho, and equivalent pathlength algorithms. The central axis dose within the low-density inhomogeneity as calculated by Monte Carlo simulation and convolution superposition decreases for small field sizes whereas it increases using the Batho and equivalent pathlength algorithms. The dose perturbation factor (DPF) is defined as the ratio of dose to a point within the inhomogeneity to the same point in a homogeneous phantom. The dose correction factor is defined as the ratio of dose calculated by an algorithm at a point to the Monte Carlo derived dose at the same point, respectively. DPF is noted to be significant for small fields and low density for all algorithms. Comparisons of the algorithms with Monte Carlo simulations is reflected in the DCF, which is close to 1.0 for the convolution-superposition algorithm. The Batho and equivalent pathlength algorithms differ significantly from Monte Carlo simulation for most field sizes and densities. Convolution superposition shows better agreement with Monte Carlo data versus the Batho or equivalent pathlength corrections. As the field size increases the DCF's for all algorithms converge toward 1.0. The largest differences in DCF are at the interface where changes in electron transport are greatest. For a 6 MV photon beam, electronic equilibrium is restored at field sizes above 3 cm diameter and all of the algorithms predict dose in and beyond the inhomogeneous region equally well. For accurate dosimetry of small fields within and near inhomogeneities, however, simple algorithms such as Batho and equivalent pathlength should be avoided.