The use of radiographic film for linear accelerator stereotactic radiosurgical dosimetry

Similar documents
Limits of Precision and Accuracy of Radiation Delivery Systems

DOSE DISTRIBUTION ANALYZE OF THE BODY STI USED MONTE CARLO METHOD

Small field diode dosimetry

A new geometric and mechanical verification device for medical LINACs

IMRT point dose measurements with a diamond detector

Intensity Modulated Radiation Therapy: Dosimetric Aspects & Commissioning Strategies

Verification of treatment planning system parameters in tomotherapy using EBT Radiochromic Film

Unrivaled, End-to-End

Topics covered 7/21/2014. Radiation Dosimetry for Proton Therapy

Verification of Relative Output Factor (ROF) Measurement for Radiosurgery Small Photon Beams

Multilayer Gafchromic film detectors for breast skin dose determination in vivo

Assessment of Dosimetric Functions of An Equinox 100 Telecobalt Machine

An anthropomorphic head phantom with a BANG polymer gel insert for dosimetric evaluation of IMRT treatment delivery

RADIATION PROTECTION IN DIAGNOSTIC AND INTERVENTIONAL RADIOLOGY. L19: Optimization of Protection in Mammography

Neutron Measurements for Intensity Modulated Radiation Therapy

Assessment of variation of wedge factor with depth, field size and SSD for Neptun 10PC Linac in Mashhad Imam Reza Hospital

DOSIMETRIC COMPARISION FOR RADIATION QUALITY IN HIGH ENERGY PHOTON BEAMS

Rotating and static sources for gamma knife radiosurgery systems: Monte Carlo studies

4 Essentials of CK Physics 8/2/2012. SRS using the CyberKnife. Disclaimer/Conflict of Interest

Introduction. Measurement of Secondary Radiation for Electron and Proton Accelerators. Introduction - Photons. Introduction - Neutrons.

Dosimetric characterization with 62 MeV protons of a silicon segmented detector for 2D dose verifications in radiotherapy

IGRT1 technologies. Paweł Kukołowicz Warsaw, Poland

Application of MCNP4C Monte Carlo code in radiation dosimetry in heterogeneous phantom

Disclosure. Outline. Machine Overview. I have received honoraria from Accuray in the past. I have had travel expenses paid by Accuray in the past.

Absolute Dosimetry. Versatile solid and water phantoms. Introduction Introduction Introduction Introtro Intro Intro

Radiochromic film dosimetry in water phantoms

Dosimetric study of 2D ion chamber array matrix for the modern radiotherapy treatment verification

JOURNAL OF APPLIED CLINICAL MEDICAL PHYSICS, VOLUME 4, NUMBER 1, WINTER 2003

ph fax

EXPERIMENTAL RESULTS IN PERCENTAGE DEPTH DOSE (PDD) DETERMINATION AT THE EXTENDED DISTANCES *

Ionization chamber dosimetry of small photon fields: a Monte Carlo study on stopping-power

Assessment of a three-dimensional (3D) water scanning system for beam commissioning and measurements on a helical tomotherapy unit

Production and dosimetry of simultaneous therapeutic photons and electrons beam by linear accelerator: a monte carlo study

STEREOTACTIC DOSE VERIFICATION PHANTOM VERSATILE STEREOTACTIC QA PHANTOMS

Technical Study. Institution University of Texas Health San Antonio. Location San Antonio, Texas. Medical Staff. Daniel Saenz. Niko Papanikolaou.

An Active Matrix Flat Panel Dosimeter AMFPD for in-phantom dosimetric measurements

Digital film dosimetry in radiotherapy and the development of analytical applications software

PROCEDURE FOR ABSORBED DOSE TO WATER DETERMINATION IN HIGH ENERGY PHOTON AND ELECTRON BEAMS BY FERROUS SULPHATE DOSIMETER AT INMRI-ENEA

Can we hit the target? Can we put the dose where we want it? Quality Assurance in Stereotactic Radiosurgery and Fractionated Stereotactic Radiotherapy

TLD as a tool for remote verification of output for radiotherapy beams: 25 years of experience

Independent corroboration of monitor unit calculations performed by a 3D computerized planning system

Diode calibration for dose determination in total body irradiation

Nuclear Associates

A preliminary clinic dosimetry study for synchrotron radiation therapy at SSRF

Additional Questions for Review 2D & 3D

Normal tissue doses from MV image-guided radiation therapy (IGRT) using orthogonal MV and MV-CBCT

AAPM Task Group 180 Image Guidance Doses Delivered During Radiotherapy: Quantification, Management, and Reduction

The influence of field size on stopping-power ratios in- and out-of-field: quantitative data for the BrainLAB m3 micro-multileaf collimator

MAX-HD SRS PHANTOM THE COMPREHENSIVE END-TO-END SRS PHANTOM SCAN PLAN LOCALIZE TREAT. distributed by:

3D printing of air slab caps for correction free small field dosimetry. Benjamin Perrett, Paul Charles, Tim Markwell, Tanya Kairn, Scott Crowe

Use of Bubble Detectors to Characterize Neutron Dose Distribution in a Radiotherapy Treatment Room used for IMRT treatments

Protons Monte Carlo water-equivalence study of two PRESAGE formulations for proton beam dosimetry J. Phys.: Conf. Ser.


Monte Carlo simulation of 192 Ir radioactive source in a phantom designed for brachytherapy dosimetry and source position evaluation.

Mania Aspradakis John Byrne Hugo Palmans John Conway Jim Warrington Karen Rosser Simon Duane

Absorbed Dose Response in Water of Kilovoltage X-rays Beams of Radiochromic Film and Thermoluminescent for Brachytherapy Dosimetry

Evaluation of Dosimetry Check software for IMRT patient-specific quality assurance

SRS Uncertainty: Linac and CyberKnife Uncertainties

Neutron dose evaluation in radiotherapy

Traceability and absorbed dose standards for small fields, IMRT and helical tomotherapy

On the impact of ICRU report 90 recommendations on k Q factors for high-energy photon beams

ROPES eye plaque dosimetry: commissioning and verification of an ophthalmic brachytherapy treatment planning system

Quality assurance program for prototype stereotactic system developed for Neptun 10 PC linac

Hot Topics in SRS: Small Field Dosimetry & Other Treatment Uncertainties. Sonja Dieterich, PhD University of California Davis

Calibration of dosimeters for small mega voltage photon fields at ARPANSA

INTRODUCTION. Material and Methods

Outline. Chapter 12 Treatment Planning Combination of Beams. Opposing pairs of beams. Combination of beams. Opposing pairs of beams

Implementing New Technologies for Stereotactic Radiosurgery and Stereotactic Body Radiation Therapy

X-Ray Guided Robotic Radiosurgery for Solid Tumors

Evaluation of Dosimetric Characteristics of a Double-focused Dynamic Micro-Multileaf Collimator (DMLC)

A spiral phantom for IMRT and tomotherapy treatment delivery verification

Plan-Specific Correction Factors for Small- Volume Ion Chamber Dosimetry in Modulated Treatments on a Varian Trilogy

Measurement of Dose to Implanted Cardiac Devices in Radiotherapy Patients

Varian Edge Experience. Jinkoo Kim, Ph.D Henry Ford Health System

Inhomogeneity effect in Varian Trilogy Clinac ix 10 MV photon beam using EGSnrc and Geant4 code system

A Dosimetric study of different MLC expansion aperture For the radiotherapy of pancreas cancer

Applications of Monte Carlo simulations to radiation dosimetry

D DAVID PUBLISHING. Uncertainties of in vivo Dosimetry Using Semiconductors. I. Introduction. 2. Methodology

Electron Beam ET - Reloaded Therapy - Reloaded

Diode Dosimetric Characteristics Assessment for In-Vivo Dosimetry in Radiotherapy Treatment

Radiosurgery by Leksell gamma knife. Josef Novotny, Na Homolce Hospital, Prague

IMRT QUESTIONNAIRE. Address: Physicist: Research Associate: Dosimetrist: Responsible Radiation Oncologist(s)

Dose distribution and dosimetry parameters calculation of MED3633 palladium-103 source in water phantom using MCNP

Non-target dose from radiotherapy: Magnitude, Evaluation, and Impact. Stephen F. Kry, Ph.D., D.ABR.

Patient dosimetry for total body irradiation using single-use MOSFET detectors

Dosimetric characterization with 62 MeV protons of a silicon segmented detector for 2D dose verifications in radiotherapy

ABSTRACTS FOR RADIOTHERAPY STANDARDS USERS MEETING. 5 th June 2007

Investigation of the clinical performance of a novel solid-state diagnostic dosimeter

CyberKnife Technology in Ablative Radiation Therapy. Jun Yang PhD Cyberknife Center of Philadelphia Drexel University Jan 2017

ARCCHECK: COMPREHENSIVE EVALUATION OF THE DIODE ARRAY PHANTOM - RULES OF THUMB FOR PHANTOM USE FOR QA. By Vibha Chaswal, Ph.D.

Dosisverifikation mit Delta 4 Discover während der Behandlung

Y FILMS DOSIMETR Nederland België / Belgique

Ion recombination corrections of ionization chambers in flattening filter-free photon radiation

SHIELDING TECHNIQUES FOR CURRENT RADIATION THERAPY MODALITIES

Use of Electronic Portal Image Detectors for Quality Assurance of Advanced Treatments

Quality Assurance of TPS: comparison of dose calculation for stereotactic patients in Eclipse and iplan RT Dose

Semiflex 3D. Always perfectly oriented. 3D Thimble Ionization Chamber for Relative and Absolute Dosimetry

Practical Reference Dosimetry Course April 2015 PRDC Program, at a glance. Version 1.0. Day 1 Day 2 Day 3 Day 4

A comparison of dose distributions measured with two types of radiochromic film dosimeter MD55 and EBT for proton beam of energy 175 MeV

The Journey of Cyberknife Commissioning

Transcription:

The use of radiographic film for linear accelerator stereotactic radiosurgical dosimetry J. L. Robar a) and B. G. Clark The University of British Columbia, Department of Physics, BC Cancer Agency, 600 West 10th Avenue, Vancouver, British Columbia, V5Z 4E6, Canada Received 24 February 1999; accepted for publication 30 June 1999 The measurement of stereotactic radiosurgical dose distributions requires an integrating, highresolution dosimeter capable of providing a spatial map of absorbed dose. Although radiographic film is an accessible dosimeter fulfilling these criteria, for larger radiotherapy photon fields the sensitivity of film emulsion exhibits significant dependencies on both depth in phantom and field size. We have examined the variation of film sensitivity over the ranges of depths and field sizes of interest in radiosurgery with a6mvphoton beam. While for large 20 cm 20 cm fields the potential error in dose due to the variation of the film response with depth reaches 15%, the corresponding maximum error for a 2.5 cm diameter radiosurgical beam is 1.5%. This uncertainty was observed to be comparable in magnitude to that produced by variation in processing conditions 1.1% and by varying the orientation of the film plane relative to the beam central axis 1.5%. The dependence of emulsion sensitivity on field size has been observed to be negligible for fields ranging in diameter from 1.0 cm to 4.0 cm. The source of the dependence of film sensitivity has been illustrated by using an EGS4 Monte Carlo simulation for large fields to illustrate significant increases in the photon spectrum below 400 kev with depth in phantom. In contrast, relative increase of this low-energy component is negligible for radiosurgical photon fields. 1999 American Association of Physicists in Medicine. S0094-2405 99 00110-8 Key words: film dosimetry, stereotactic radiosurgery, dose verification, Monte Carlo, sensitometry I. INTRODUCTION Linear accelerator-based stereotactic radiosurgery is a treatment modality employing narrowly-collimated, 6 10 MV photon beams to deliver a high dose to intracranial lesions with a sharp dose gradient at the periphery of the treatment volume. Conventionally, a number of circular convergent beams or beam arcs are used to create an approximately spherical or ellipsoidal distribution of dose. Many lesions are highly irregular in shape, however, and recent advances in radiosurgery have focused upon improving the conformity between the target volume and the delivered dose distribution. While software treatment planning systems provide analytic, three-dimensional calculations of even complex dose distributions, the required spatial and numerical tolerances for dose delivery in radiosurgery are demanding 1 mm and 5%, respectively 1 and dose calculations must be verified through direct measurement. Dosimetric measurements are required both in commissioning the radiosurgical hardware and software and in the course of periodic quality assurance. Dosimetry for conventional circular-beam radiosurgery is challenging 2 and the criteria for an appropriate dosimeter become more stringent for measurement of irregular dose distributions. First, the required spatial accuracy in localizing and treating the lesion necessitates a dosimeter with high sub-millimeter spatial resolution. The sensitive volume of the detector must be sufficiently low in order to record the dose over the sharp fall-off region of the distribution while minimizing inaccuracies caused by volume averaging of the detector size 3 and the presence of lateral electronic disequilibrium. 4 Second, the dosimeter must provide a spatial map of the dose distribution rather than the dose at a single point, particularly for cases where it is necessary to verify conformal distributions. Finally, the dosimeter must integrate over time to reflect the total dose administered by a number of beams or beam arcs. The requirement for a simultaneous measurement of a spatial distribution of dose alone precludes the use of common dosimeters such as ionization chambers, thermoluminescent dosimeters TLDs, and diodes. In order to verify conformal radiosurgical dose distributions, we have evaluated the accuracy of radiographic film. Radiographic film is an accessible, integrating dosimeter capable of recording high spatial-resolution up to 20 line-pairs per mm dosimetric information in two dimensions. Radiographic film is affordable and the optical density recorded on the film is readily digitized using scanning densitometers found commonly in radiotherapy centers. The use of radiographic film for the dosimetry of larger i.e., 10 cm 10 cm photon radiotherapy beams is complicated due to a dependence of film emulsion sensitivity on depth within the phantom. This dependence is caused by a relative increase with depth in the population of lowerenergy, scattered photons in the spectrum and the subsequent photoelectric absorption of these photons by the film emulsion. This effect becomes most pronounced with increases in the photon population in the energy region below approximately 400 kev where the mass attenuation coefficient of 2144 Med. Phys. 26 10, October 1999 0094-2405/99/26 10 /2144/7/$15.00 1999 Am. Assoc. Phys. Med. 2144

2145 J. L. Robar and B. G. Clark: The use of radiographic film 2145 emulsion increases rapidly, approximately as 1/(h ) 3. As the low-energy spectral component increases, the tissue equivalence of film emulsion diverges from that of tissueequivalent materials due to an approximate Z 3 dependence of the photoelectric interaction cross-section. In practice, this dependence causes the sensitometric calibration curve of the film to vary with the depth in phantom at which the film is exposed. To overcome this problem, Hale et al. 5 have measured a series of sensitometric curves, one for each of a range of depths, and have achieved 2% accuracy by using the appropriate curve for dosimetric studies. Williamson et al. 6 have parametrized film sensitivity as a function of depth in order to compensate for this dependence, and have obtained 3% accuracy with this technique. Rather than correcting measured optical densities for the effect of phantom scatter, Burch et al. 7 have proposed a method to remove low-energy, lateral scatter by placing lead foils on either side of the film plane. These types of corrections increase the overall complexity of film dosimetry with the necessity of an increased quantity of calibration data, special modifications to the experimental setup or post-processing of acquired data to correct for the depth dependence of emulsion sensitivity. While the data provided by Burch et al. 7 demonstrate the potential inaccuracy of film dosimetry for larger 6 cm 6 cm radiotherapy fields, these results have also indicated that the magnitude of the depth dependence of film sensitivity is diminished as the field size is reduced. The present work draws upon these findings with the reasoning that because stereotactic radiosurgery employs the smallest field sizes produced commonly by clinical linear accelerators, the corrections required in converting the recorded film optical density to dose should be minimized. There are examples in the literature of the use of radiographic film for accurate measurement 8,9 for narrow photon beams, but the depth and field size dependencies of film sensitivity have not been quantitated for radiosurgical dosimetry. The focus of the current work is the measurement of the variation of radiographic emulsion sensitivity over the range of depths up to 20 cm and field sizes 1.0 cm to 4.0 cm diameter of interest for stereotactic radiosurgery. This variation is quantified by comparing optical density-to-dose conversion curves as a function of depth in phantom and field size. Since the dependence of emulsion sensitivity on depth has been attributed to shifts in the photon spectrum, 10 we have interpreted experimental results by using the EGS4 FLURZ Monte Carlo simulation 11,12 to examine the variation of the spectrum with depth in phantom. In particular this simulation has permitted a comparison of spectral increases in the energy region below 400 kev for both large i.e., radiotherapy-sized and small radiosurgical-sized beams. Additional experiments were conducted to examine whether emulsion sensitivity depends on the film orientation relative to the central axis of the beam. This possible dependence becomes important when measuring dose distributions created by multiple beams or beam-arcs, where the orientation of the beams relative to the film plane may be arbitrary. Finally, the reproducibility of film dosimetry was assessed by comparing identically-exposed films. II. METHODS AND MATERIALS A. Film exposure technique In order to examine possible dependencies of emulsion exposed by radiosurgical beams, a series of sensitometric curves was established for a range of depths and field sizes. For each exposure a 10 in. 12 in. sheet of Kodak X-Omat V film from a single batch #194 05 2 was sealed within a light tight Solid Water Gammex RMI, Inc. cassette. The cassette consists of two 2 cm-thick slabs of solid water sealed around three edges by nylon screws and a rubber O-ring. In order to minimize the occurrence of air gaps, the film was removed from its envelope and paper liner under safelight before insertion into the unsealed edge cassette. On the linear accelerator couch, 10 cm of solid water phantom material was added both above and below the cassette to simulate the amount of scattering material that would be adjacent to the film if it were placed at the center an anthropomorphic head phantom. The 6 MV photon beam used for all exposures was identical to that generated routinely for radiosurgical treatment. The linear accelerator Clinac 2100C/D, Varian Associates was used with titanium alloy tertiary collimators BrainLAB, GmbH, which are tapered according to beam divergence and range in size to produce fields from 1.0 cm to 4.0 cm diameter at the isocenter. The linear accelerator secondary collimator jaws were fixed at 5 cm 5 cm for all tertiary collimator sizes. B. Film processing and measurement of optical density To minimize possible variations due to film processing conditions, a Kodak X-Omat RP processor was used for which the throughput is very high and quality assurance is performed daily. Developer temperature fluctuated by less than 0.5 F between processing sessions. Optical density was measured using a scanning densitometer Wellhöfer WP102, Schwarzenbruck, Germany. The densitometer infrared light source beam width is approximately 0.8 mm in diameter and is mounted on a translation stage capable of approximately 0.2 mm precision. Base-plus-fog optical density was subtracted from scanned optical densities for each film. C. Depth-dependence of film sensitivity In order to measure dose in phantom along the depth axis, the gantry was rotated to 90 from the vertical and room lasers were used to align the beam central axis and plane of the film to within approximately 0.05 cm. All films were exposed using the median stereotactic collimator field size of 2.5 cm. To acquire data to establish sensitometric curves in sufficient detail, 14 films were exposed with different monitor unit MU settings resulting in maximum doses to phantom ranging between 5 and 385 cgy. Although the goal of

2146 J. L. Robar and B. G. Clark: The use of radiographic film 2146 this series of exposures was to establish the sensitometric curves over a dose range of approximately 0 140 cgy, the higher dose settings used were required in order to obtain a sufficient range of dose at greater depths in phantom. At least three films were exposed for each dose setting in order to minimize variations due to processing and scanning. This procedure was then repeated using a 20 cm 20 cm field size for comparison of small-field and large-field results. For each exposure made with the film in this parallel orientation, a two-dimensional profile of optical density is recorded along the depth dimension of the phantom. In scanning the film, it is necessary to ensure that the densitometer travels along the recorded location of the central axis i.e., so that the film is not tilted relative to the direction of scanning. To align the film, the densitometer was first used to record beam profiles perpendicular to the central axis at depths of 5 cm and 10 cm, and the geometric centers of each profile were calculated from the scanned data. A discrepancy between the calculated location of these profile centers indicated that the film was angled relative to the scanning direction. Once aligned, the film was scanned down the central axis location and the optical density was measured along the depth axis from the surface to a depth of 20.0 cm in 0.1 cm increments. Each scan produces a series of optical density values recorded at each of 200 depths. In order to relate these optical densities to absolute doses, the dose to water corresponding to each of these depths and for each MU setting was calculated from D d,a MU PDD d,a S t A, 1 100 where D(d,A) is the dose at depth d for field size A at the phantom surface, MU is the number of monitor units given, PDD(d,A) is the percent depth dose and S t (A) is the total scatter factor. 13 Both S t (A) and PDD(d,A) were measured in a water phantom model WP 600 Wellhöfer, Schwarzenbruck, Germany using a p-type silicon electron diode 14 Scanditronix, Uppsala, Sweden. The PDD was measured in 0.1 cm increments to a maximum depth of 25.0 cm and was linearly interpolated when required to obtain values at arbitrary depths. After scanning all 14 films the resultant set of data consisted of 14 optical density, dose pairs at each of 200 depths in phantom. These data were therefore sufficient to generate sensitometric curves at each of 200 different depths in phantom, with each curve defined by 14 points. D. Field size-dependence of film sensitivity Similar sets of sensitometric curves were established, as described above, for both the minimum 1.0 cm-diameter and maximum 4.0 cm-diameter field sizes. These data facilitate the examination of variation between the sensitometric curves, at any arbitrary depth, for the range of field sizes used routinely in radiosurgery. E. Effect of film orientation In order to quantify a possible dependence of film orientation, exposures were also made with the film plane perpendicular to the central axis at a single depth of 10 cm in phantom. Five separate films were exposed with given doses ranging from 5 to 385 cgy. At a depth of 10 cm this dose range extends into the saturation region of the film, but it was useful to define the sensitometric curve over as wide a range as possible for accurate curve-fitting. For each exposure the optical density was measured at the center of the spot pattern recorded on the film. Using Eq. 1 to calculate absolute doses for each exposure, a sensitometric curve was established for this perpendicular geometry and compared with the sensitometric curve for the parallel geometry. Films for both orientations were developed in the same batch to minimize the effect of processing variability. F. Film reproducibility During each experimental session a minimum of three control films were exposed. The control exposure was made using 2.5 cm-diameter collimator by administering a maximum dose of 87 cgy to a film oriented parallel to the central axis. The optical density curves obtained by scanning the film along the depth dimension at the location of the central axis were compared. Because these control exposures were repeated during a single session and also over the course of approximately six months, realistic measures of both film-to-film and session-to-session variation were obtained. In addition to processing variations, these values also include contributions such as drift and fluctuation of the densitometer signal and nonuniformity of the film emulsion. G. Monte Carlo simulation To examine the variation of the photon spectrum with depth and to compare large- and small-field sizes in terms of this variation, the EGS4 Monte Carlo FLURZ code was used to model a cylindrical water phantom as shown in Fig. 1. The phantom is 40.0 cm in diameter and 40.0 cm deep. Both large 25 cm-diameter and small 2.5 cm-diameter beams were modeled for comparison. A published 6 MV incident spectrum 15 was selected to be representative of that used in the experiments. The full photon spectrum was scored along the depth axis at 1.0 cm increments within a 1.0 cm-diameter region at the center of the phantom. 15 10 6 incident histories were recorded for each beam. The PRESTA algorithm 16 was used to model electron transport, and the electron cut-off ECUT and photon cut-off PCUT parameters were set to 0.521 MeV and 0.010 MeV, respectively. III. RESULTS A. Depth-dependence of film sensitivity Figure 2 a shows the sensitometric curves obtained for the large 20 cm 20 cm photon beam at depths of 1.0 cm,

2147 J. L. Robar and B. G. Clark: The use of radiographic film 2147 FIG. 1. Illustration of the phantom simulated using the EGS4 FLURZ user code. Two 6 MV photon beams, with diameters of 25 cm and 2.5 cm, were simulated separately for a comparison of the results. Total photon spectra were scored in a1cm-diameter region along the central axis at 1 cm depth increments. 10.0 cm and 20.0 cm. For each depth a curve was fitted to the measured sensitometric data using the single-target/ single-hit theory equation, 7 OD OD sat 1 10 D, where OD and D are the measured optical density and given dose, respectively. The saturation density of the film, OD sat, was estimated by delivering a large dose 500 cgy to a film in phantom, and was held constant in the fitting algorithm. The fitting parameter represents emulsion sensitivity and was allowed to vary. For this large field, the curves diverge markedly, indicating that using a sensitometric curve that is not depth-specific would introduce significant error in converting optical density to dose. In order to minimize this error to below approximately 10%, for example, it is necessary to confine the dose range to less than 80 cgy. In contrast, Fig. 2 b shows sensitometric curves corresponding to the same depths for the radiosurgical 2.5 cm-diameter field. For this small field the curves agree to within the reproducibility of film development and scanning discussed below. The fitted curves using the single-hit/single-target equation correspond to the data closely, in all cases resulting in correlation coefficient values greater than 0.9957. By determining the value of for each curve, film sensitivity was obtained as a function of depth as illustrated for the 20 cm 20 cm and 2.5 cm-diameter fields Fig. 3. As expected from the disparity in the sensitometric curves in Fig. 2 a, the 2 FIG. 2. The sensitometric curves established experimentally for a the 20 cm 20 cm radiotherapy field and b the 2.5 cm-diameter radiosurgical field, shown for depths in solid water phantom of 1.0 cm, 10.0 cm and 20.0 cm. film sensitivity for the large 20 cm 20 cm field increases systematically with depth. While small fluctuation of the values of is apparent for the radiosurgical field, no systematic variation of sensitivity with depth is apparent. FIG. 3. The curve-fitting parameter, representing emulsion sensitivity, as a function of depth for the 20 cm 20 cm radiotherapy field and the 2.5 cm-diameter radiosurgical field.

2148 J. L. Robar and B. G. Clark: The use of radiographic film 2148 FIG. 4. Sensitometric curves shown for the 1.0 cm, 25.0 cm and 40.0 cmdiameter radiosurgical beams. All curves correspond to a single depth of 10.0 cm, chosen as typical for radiosurgery treatment. FIG. 6. The effect of changing the orientation of the film from parallel circles to perpendicular squares to the central axis. The maximum disparity is 1.5%. B. Field-size dependence of film sensitivity Figure 4 illustrates the measured sensitometric curves for the range of field sizes used typically for linear accelerator stereotactic radiosurgery. These data are shown for a single depth of 10.0 cm. Curves fitted using Eq. 2 are congruent, indicating that there is no significant variation in sensitivity with field size over this range. This congruence was observed for all depths in phantom ranging from 1.0 cm to 20.0 cm. C. Film reproducibility A subset of the film reproducibility data is shown in Fig. 5 which compares the optical density recorded on six films exposed identically and processed in three batches of two films, where each batch was processed on a separate day over the course of approximately six months. At a typical depth for radiosurgery of 10.0 cm, the intra-session maximum variation in optical density between films processed in one batch is 0.85% where variation has been quantified by calulating OD max OD min /OD mean. The maximum inter-session variation in optical density between films processed on separate occasions is 2.3%. Using the fitted sensitometric curve shown in Fig. 4, this translates into intrasession and inter-session variations in the calculated dose of 1.1% and 2.9%, respectively. To determine the component of this fluctuation resulting from the film digitization process, a single film was scanned repeatedly with the densitometer. This source of variation was found to be negligible compared to the total variation in optical density between control films. For lower OD values, however 0.5, the presence of noise became apparent. D. Effect of film orientation Figure 6 compares the sensitometric curves obtained with the film positioned parallel and perpendicular to the central axis at a depth of 10.0 cm in phantom. Over this wide dose range, the maximum discrepancy between these curves is 1.5%, slightly greater than the inter-session variation due to processing. FIG. 5. A subset of the reproducibility data, showing the optical density measured as a function of depth along the location of the central axis for six films processed and measured in three separate sessions. The inter-session and intra-session uncertainties in dose caused by this variation are 2.9% and 1.1%, respectively. E. Monte Carlo simulation The photon spectra computed using EGS4 FLURZ are shown in Fig. 7 and serve to explain the source of the increase in emulsion sensitivity with depth for larger fields, as well as the comparatively low depth-dependence for radiosurgical fields. Figure 7 a shows the total primary plus scattered photon spectra for a large 25 cm-diameter field on central axis at depths of 1.0 cm, 5.0 cm, and 20.0 cm. Each spectrum has been normalized to its own maximum to illustrate relative spectral changes. These data depict a consistent increase in the lower-energy ( 400 kev component of the spectrum with depth. Since the mass attenuation coefficient for film emulsion increases rapidly in this region, this would result in an increase with depth in the proportion of absorbed photons and therefore an enhanced sensitivity. Changing only the field size in the simulation to a radiosurgical-sized diameter of 2.5 cm produces the spectra

2149 J. L. Robar and B. G. Clark: The use of radiographic film 2149 FIG. 8. The 6 MV incident energy spectrum solid line shown on the same energy scale as the mass attenuation coefficients for Kodak emulsion, polyester film base and water. The shaded region in the spectrum indicates the population of photons below 400 kev. FIG. 7. The 6 MV energy spectra shown at depths in phantom of 1.0 cm, 5.0 cm and 20.0 cm for a the 25.0 cm-diameter field and b the 2.5 cmdiameter field. The spectra were in finer energy bins below 1 MeV in order to visualize the variation of this low-energy region with depth. illustrated in Fig. 7 b. In contrast to the large-field case, the spectra for the radiosurgical field show a negligible increase in the low-energy population of photons. IV. DISCUSSION Although film dosimetry is most problematic for lowerenergy photon beams e.g., cobalt-60 and 4 MV due to their lower primary-to-scatter ratio, Fig. 2 a emphasizes the importance of corrections in converting scanned optical density to dose even for the 6 MV beam. For this 20 cm 20 cm field, using a depth-unspecific curve may introduce an error in dose of up to 15% for the dose range examined in this study. This error is confined to below 10% only for doses below 80 cgy. For the 2.5 cm-diameter radiosurgical field, using a sensitometric curve established for a depth of 1.0 cm to calculate a dose at 20.0 cm depth would result in a maximum error of 1.5%, which is comparable in magnitude to the intra-session variability. This invariance of the calibration curve with depth greatly simplifies the task of film dosimetry for 6 MV radiosurgical beams in two respects. First, corrections to compensate for a depth-dependence of the film emulsion sensitivity are not required. Second, since a single calibration curve can be used over a range of depths, the quantity of data required in calibrating the film is greatly reduced. As depicted by the data in Fig. 7, the photon spectra for both the large 25.0 cm-diameter and small 2.5 cmdiameter fields exhibit variation with depth in phantom. These results agree with the simulation findings by Heydarian et al., 2 which indicate that the mean photon energy of 6 MV radiosurgical beams increases with depth due to beam hardening, but for large beams 20 cm-diameter the mean photon energy decreases with depth up to 15 cm due to an increase of phantom scatter. In terms of the potential overresponse of film emulsion, however, only spectral changes below approximately 400 kev are crucial. Figure 8 shows the 6 MV spectrum on the same energy scale as the mass attenuation coefficients for water, film base and film emulsion. The shaded region of this incident spectrum indicates the small population of photons ( 2% of the total spectrum with energies below 400 kev. As indicated by the Monte Carlo results in Fig. 7 b for radiosurgical fields, this lowenergy tail does not increase significantly with depth in phantom. Over the range of depths of interest in radiosurgery to 20 cm the vast majority of photons in the spectrum exist in the energy range between 400 kev to 6 MeV. It should be noted that this simulation was performed only for the median radiosurgical collimator size of 2.5 cm, and potential increased softening of the spectrum in phantom for the larger radiosurgical collimators was not examined. However, the recent Monte Carlo simulation of a linac radiosurgery unit by Verhaegen et al. 17 suggests that the mean photon energy at depth of dose maximum decreases by only 0.45 MeV when the aperture diameter is increased from 0.5 cm to 5.0 cm. Figures 2 and 4 show that sensitometric curve is approximately linear for doses below 40 cgy. As indicated by Evans and Schreiner, 18 restricting the dose to this region offers the advantage of direct conversion of optical density to dose. As illustrated by the data in Fig. 2 a for the larger radiotherapy field, this constraint would also minimize the error caused by the variation of film sensitivity with depth, and thus would alleviate the need to acquire a family of calibration curves over a range of depths. In measuring radiosurgical dose distributions, however, in many cases it may be difficult to con-

2150 J. L. Robar and B. G. Clark: The use of radiographic film 2150 fine the total dose to the linear region of the sensitometric curve. For example, for multiple convergent-arc techniques, a minimum achievable monitor unit-per-degree rate is imposed by the linear accelerator, and thus the total dose may necessarily exceed 40 cgy. Similarly, treatments involving multiple, static conformal fields would necessitate a very low dose per field, and thus the uncertainty inherent in the dose per beam would become appreciable. Therefore, to facilitate dosimetry over a larger dose range, we generally acquire sufficient calibration data to extend the curve-fitting over the nonlinear region. It should be noted, however, that the nonlinear shape of the sensitometric curve is specific to Kodak XOmat-V film. Cheng and Das 19 and Cadman 20 have reported on a relatively new type of silver-halide radiographic film CEA TVS which is characterized by a linear sensitometric curve up to 90 cgy, followed by abrupt saturation. While Cheng and Das have shown that the relationship between optical density and dose does not vary for 4 18 MV photon spectra for measurements made at the depth of maximum dose in phantom, possible variation of sensitivity with depth in phantom has not been examined. The low variation of emulsion sensitivity with depth for small ( 4.0 cm-diameter 6 MV beams will offer an advantage in terms of dosimetric accuracy only if the variation due to processing conditions does not prevail. In this study, although the developer temperature was held constant to within 0.5 F and the replenishment rate was kept as continuous as possible, the data in Fig. 5 suggest that in order to achieve a variation of 2%, it is necessary to process recalibration films in the same batch as the dosimetric films. V. CONCLUSION This work was motivated by a need for an accessible, high-resolution dosimeter for stereotactic radiosurgery, and in particular, for small conformal fields. The results indicate that radiographic film, while providing the required spatial resolution, is particularly applicable to small 6 MV photon beams. Dependencies of emulsion sensitivity on depth and field size which limit film dosimetry for larger photon fields are sufficiently minor for radiosurgical fields to warrant the use of a single sensitometric curve. Monte Carlo results indicate that even at depth in phantom, the majority of photons in the spectrum have energies above 400 kev, and therefore film does not greatly compromise tissue equivalence. Other potential sources of uncertainty, such as the variability of processing conditions and dependence of film sensitivity with orientation do not preclude the use of radiographic film for the measurement of radiosurgical dose distributions produced by a6mvbeam. ACKNOWLEDGMENTS The authors would like to thank Kurt Luchka for valuable input concerning film dosimetry and Dr. Corey Zankowski for his assistance in Monte Carlo calculations. This work forms part of a doctoral dissertation by J. L. R. and it has been funded by the Natural Sciences and Engineering Council of Canada Scholarship No. PGSB-207772-1998, The BC Cancer Agency and The University of British Columbia. a Electronic mail: jrobar@bccancer.bc.ca 1 Task Group 42 Radiation Therapy Committee, American Association of Physicists in Medicine, AAPM Report No. 54: Stereotactic Radiosurgery American Institute of Physics, Woodbury, NY, 1995. 2 M. Heydarian, P. W. Hoban, and A. H. Beddoe, A comparison of techniques in stereotactic radiosurgery, Phys. Med. Biol. 41, 93 110 1996. 3 P. D. Higgins, C. H. Sibata, L. Siskind, and J. W. Sohn, Deconvolution of detector size effect for small field measurement, Med. Phys. 22, 1663 1666 1995. 4 R. K. Rice, J. L. Hansen, G. K. Svensson, and R. L. Siddon, Measurements of dose distributions in small beams of 6 MV x-rays, Phys. Med. Biol. 32, 1087 1099 1987. 5 J. I. Hale, A. T. Kerr, and P. C. Shragge, Calibration of film for accurate megavoltage photon dosimetry, Med. Dosimetry 19, 43 46 1994. 6 J. F. Williamson, F. M. Khan, and S. C. Sharma, Film dosimetry of megavoltage photon beams: A practical method of isodensity-to-isodose curve conversion, Med. Phys. 8, 94 99 1981. 7 S. E. Burch, K. J. Kearfott, J. H. Trueblood, W. C. Sheils, J. Inhwan Yeo, and C. K. C. Wang, A new approach to film dosimetry for high energy photon beams: Lateral scatter filtering, Med. Phys. 24, 775 783 1997. 8 B. E. Bjärngard, J.-S. Tsai, and R. K. Rice, Doses on the central axes of narrow 6-MV x-ray beams, Med. Phys. 17, 794 799 1990. 9 D. A. Mazal, J. C. Rosenwald, G. Gaboriaud, and M. C. Plazas, Charactéristiques dosimétriques des petits faisceaux de photons he haute énergie: measure et modélisation, Proceedings of the XXVII Congres Société Francaise des Physiciens d Hôpital, 1988. 10 J. Greening, The photographic action of x-rays, Proc. Phys. Soc. London, Sect. B 64, 977 992 1951. 11 W. R. Nelson, H. Hirayama, and D. W. O. Rogers, The EGS4 Code System, Stanford Linear Accelerator Center, Stanford, CA, 1985, Report SLAC-265. 12 C. Malamut, D. W. O. Rogers, and A. F. Bielajew, Calculation of water/air stopping-power ratios using EGS4 with explicit treatment of electron positron differences, Med. Phys. 18, 1222 1228 1991. 13 F. M. Khan, The Physics of Radiation Therapy, 2nd ed. Williams and Wilkins, Baltimore, MD, 1994. 14 A. S. Beddar, D. J. Mason, and P. F. O Brien, Absorbed dose perturbation caused by diodes for small field dosimetry, Med. Phys. 21, 1075 1079 1994. 15 R. Mohan, C. Chui, and L. Lidofsky, Energy and angular distributions of photons from medical linear accelerators, Med. Phys. 12, 592 597 1985. 16 A. Bielajew and D. W. O. Rogers, PRESTA: the parameter reduced electron-step transport algorithm for electron Monte Carlo transport, Nucl. Instrum. Methods Phys. Res. B 18, 165 181 1987. 17 F. Verhaegen, I. J. Das, and H. Palmans, Monte Carlo dosimetry study of a 6 MV stereotactic radiosurgery unit, Phys. Med. Biol. 43, 2755 2768 1998. 18 M. D. C. Evans and L. J. Schreiner, A simple technique for film dosimetry, Radiother. Oncol. 23, 265 267 1992. 19 C.-W. Cheng and I. J. Das, Dosimetry of high energy photon and electron beams with CEA films, Med. Phys. 23, 1225 1232 1996. 20 P. Cadman, Use of CEA TVS film for measuring high energy dose distribution photon beams, Med. Phys. 25, 1435 1437 1998.