Radiotherapy treatment verification using radiological thickness measured with an amorphous silicon electronic portal imaging device: Monte Carlo simulation and experiment

This work validates the use of an amorphous-silicon, flat-panel electronic portal imaging device (a-Si EPID) for use as a gauge of patient or phantom radiological thickness, as an alternative to dosimetry. The response of the a-Si EPID is calibrated by adapting a technique previously applied to scanning liquid ion chamber EPIDs, and the stability, accuracy and reliability of this calibration are explored in detail. We find that the stability of this calibration, between different linacs at the same centre, is sufficient to justify calibrating only one of the EPIDs every month and using the calibration data thus obtained to perform measurements on all of the other linacs. Radiological thickness is shown to provide a reliable means of relating experimental measurements to the results of BEAMnrc Monte Carlo simulations of the linac-phantom-EPID system. For these reasons we suggest that radiological thickness can be used to verify radiotherapy treatment delivery and identify changes in the treatment field, patient position and target location, as well as patient physical thickness.     

Electron beam treatment verification using measured and Monte Carlo predicted portal images

Electron beam treatments may benefit from techniques to verify patient positioning and dose delivery. This is particularly so for complex techniques such as mixed photon and electron beam radiotherapy and electron beam modulated therapy. This study demonstrates that it is possible to use the bremsstrahlung photons in an electron beam from a dual scattering foil linear accelerator to obtain portal images of electron beam treatments. The possibility of using Monte Carlo (MC) simulations to predict the electron beam treatment portal images was explored. The MC code EGSnrc was used to model a Varian CL21EX linear accelerator (linac) and to characterize the bremsstrahlung photon production in the linac head. It was found that the main sources of photons in the electron beam are the scattering foils, the applicator and the beam-shaping cut-out. Images were acquired using the Varian CL21EX linac and the Varian aS500 electronic portal imager (EPI); four electron energies (6, 9, 12, 16 MeV), and different applicator and cut-out sizes were used. It was possible to acquire images with as little as 10.7 MU per image. The contrast, the contrast-to-noise ratio (CNR), the signal-to-noise ratio (SNR), the resolution and an estimate of the modulated transfer function (MTF) of the electron beam portal images were computed using a quality assurance (QA) phantom and were found to be comparable to those of a 6 MV photon beam. Images were also acquired using a Rando anthropomorphic phantom. MC simulations were used to model the aS500 EPID and to obtain predicted portal images of the QA and Rando phantom. The contrast in simulated and measured portal images agrees within +/-5% for both the QA and the Rando phantom. The measured and simulated images allow for a verification of the phantom positioning by making sure that the structure edges are well aligned. This study suggests that the Varian aS500 portal imager can be used to obtain patient portal images of electron beams in the scattering foil linacs.     

Off-axis dose response characteristics of an amorphous silicon electronic portal imaging device

Amorphous silicon (a-Si) electronic portal imaging devices (EPIDs) have typically been calibrated to dose at central axis (CAX). Division of acquired images by the flood-field (FF) image that corrects for pixel sensitivity variation as well as open field energy-dependent off-axis response variation should result in a flat EPID response over the entire matrix for the same field size. While the beam profile can be reintroduced to the image by an additional correction matrix, the CAX EPID response to dose calibration factor is assumed to apply to all pixels in the detector. The aim of this work was to investigate the dose response of the Varian aS500 amorphous silicon detector across the entire detector area. First it was established that the EPID response across the panel became stable (within approximately 0.2%) for MU settings greater than approximately 200 MU. The EPID was then FF calibrated with a high MU setting of approximately 400 for all subsequent experiments. Whole detector images with varying MU settings from 2-500 were then acquired for two dose rates (300 and 600 MU/min) for 6 MV photons for two EPIDs. The FF corrected EPID response was approximately flat or uniform across the detector for greater than 100 MU delivered (within 0.5%). However, the off-axis EPID response was greater than the CAX response for small MU irradiations, giving a raised EPID profile. Up to 5% increase in response at 20 cm off-axis compared to CAX was found for very small MU settings for one EPID, while it was within 2% for the second (newer) EPID. Off-axis response nonuniformities attributed to detector damage were also found for the older EPID. Similar results were obtained with the EPID at 18 MV energy and operating in asynchronous mode (acquisition not synchronized with beam pulses), however the profiles were flatter and more irregular for the small MU irradiations. By moving the detector laterally and repeating the experiments, the increase in response off-axis was found to depend on the pixel position relative to the beam CAX. When the beam was heavily filtered by a phantom the off-axis response variation was reduced markedly to within 0.5% for all MU settings. Independent measurements of off-axis point doses with ion chamber did not show any change in off-axis factor with MUs. Measurements of beam quality (TMR20-10) for MU settings of 2, 5, and 100 at central axis and at 15 cm off-axis could not explain the effect. The response change is unlikely to be significant for clinical IMRT verification with this imaging/acclerator system where MUs are of the order of 100-300, provided the detector does not exhibit radiation damage artifacts.     

Dose-response and ghosting effects of an amorphous silicon electronic portal imaging device

The purpose of this study was to investigate the dose-response characteristics, including ghosting effects, of an amorphous silicon-based electronic portal imaging device (a-Si EPID) under clinical conditions. EPID measurements were performed using one prototype and two commercial a-Si detectors on two linear accelerators: one with 4 and 6 MV and the other with 8 and 18 MV x-ray beams. First, the EPID signal and ionization chamber measurements in a mini-phantom were compared to determine the amount of buildup required for EPID dosimetry. Subsequently, EPID signal characteristics were studied as a function of dose per pulse, pulse repetition frequency (PRF) and total dose, as well as the effects of ghosting. There was an over-response of the EPID signal compared to the ionization chamber of up to 18%, with no additional buildup layer over an air gap range of 10 to 60 cm. The addition of a 2.5 mm thick copper plate sufficiently reduced this over-response to within 1% at clinically relevant patient-detector air gaps (> 40 cm). The response of the EPIDs varied by up to 8% over a large range of dose per pulse values, PRF values and number of monitor units. The EPID response showed an under-response at shorter beam times due to ghosting effects, which depended on the number of exposure frames for a fixed frame acquisition rate. With an appropriate build-up layer and corrections for dose per pulse, PRF and ghosting, the variation in the a-Si EPID response can be reduced to well within +/- 1%.     

Improving IMRT quality control efficiency using an amorphous silicon electronic portal imager

An amorphous silicon electronic portal imaging device (EPID) has been investigated to determine its usefulness and efficiency for performing linear accelerator quality control checks specific to step and shoot intensity modulated radiation therapy (IMRT). Several dosimetric parameters were measured using the EPID: dose linearity and segment to segment reproducibility of low dose segments, and delivery accuracy of fractions of monitor units. Results were compared to ion chamber measurements. Low dose beam flatness and symmetry were tested by overlaying low dose beam profiles onto the profile from a stable high-dose exposure and visually checking for differences. Beam flatness and symmetry were also calculated and plotted against dose. Start-up reproducibility was tested by overlaying profiles from twenty successive two monitor unit segments. A method for checking the MLC leaf calibration was also tested, designed to be used on a daily or weekly basis, which consisted of summing the images from a series of matched fields. Daily images were coregistered with, then subtracted from, a reference image. A threshold image showing dose differences corresponding to > 0.5 mm positional errors was generated and the number of pixels with such dose differences used as numerical parameter to which a tolerance can be applied. The EPID was found to be a sensitive relative dosemeter, able to resolve dose differences of 0.01 cGy. However, at low absolute doses a reproducible dosimetric nonlinearity of up to 7% due to image lag/ghosting effects was measured. It was concluded that although the EPID is suitable to measure segment to segment reproducibility and fractional monitor unit delivery accuracy, it is still less useful than an ion chamber as a tool for dosimetric checks. The symmetry/flatness test proved to be an efficient method of checking low dose profiles, much faster than any of the alternative methods. The MLC test was found to be extremely sensitive to sudden changes in MLC calibration but works best with a composite reference image consisting of an average of five successive days' images. When used in this way it proved an effective and efficient daily check of MLC calibration. Overall, the amorphous silicon EPID was found to be a suitable device for IMRT QC although it is not recommended for dosimetric tests. Automatic procedures for low monitor unit profile analysis and MLC leaf positioning yield considerable time-savings over traditional film techniques.     

A Monte Carlo based three-dimensional dose reconstruction method derived from portal dose images

The verification of intensity-modulated radiation therapy (IMRT) is necessary for adequate quality control of the treatment. Pretreatment verification may trace the possible differences between the planned dose and the actual dose delivered to the patient. To estimate the impact of differences between planned and delivered photon beams, a three-dimensional (3-D) dose verification method has been developed that reconstructs the dose inside a phantom. The pretreatment procedure is based on portal dose images measured with an electronic portal imaging device (EPID) of the separate beams, without the phantom in the beam and a 3-D dose calculation engine based on the Monte Carlo calculation. Measured gray scale portal images are converted into portal dose images. From these images the lateral scattered dose in the EPID is subtracted and the image is converted into energy fluence. Subsequently, a phase-space distribution is sampled from the energy fluence and a 3-D dose calculation in a phantom is started based on a Monte Carlo dose engine. The reconstruction model is compared to film and ionization chamber measurements for various field sizes. The reconstruction algorithm is also tested for an IMRT plan using 10 MV photons delivered to a phantom and measured using films at several depths in the phantom. Depth dose curves for both 6 and 10 MV photons are reconstructed with a maximum error generally smaller than 1% at depths larger than the buildup region, and smaller than 2% for the off-axis profiles, excluding the penumbra region. The absolute dose values are reconstructed to within 1.5% for square field sizes ranging from 5 to 20 cm width. For the IMRT plan, the dose was reconstructed and compared to the dose distribution with film using the gamma evaluation, with a 3% and 3 mm criterion. 99% of the pixels inside the irradiated field had a gamma value smaller than one. The absolute dose at the isocenter agreed to within 1% with the dose measured with an ionization chamber. It can be concluded that our new dose reconstruction algorithm is able to reconstruct the 3-D dose distribution in phantoms with a high accuracy. This result is obtained by combining portal dose images measured prior to treatment with an accurate dose calculation engine.     

Dosimetric investigation and portal dose image prediction using an amorphous silicon electronic portal imaging device.

A two step algorithm to predict portal dose images in arbitrary detector systems has been developed recently. The current work provides a validation of this algorithm on a clinically available, amorphous silicon flat panel imager. The high-atomic number, indirect amorphous silicon detector incorporates a gadolinium oxysulfide phosphor scintillating screen to convert deposited radiation energy to optical photons which form the portal image. A water equivalent solid slab phantom and an anthropomorphic phantom were examined at beam energies of 6 and 18 MV and over a range of air gaps (approximately 20-50 cm). In the many examples presented here, portal dose images in the phosphor were predicted to within 5% in low-dose gradient regions, and to within 5 mm (isodose line shift) in high-dose gradient regions. Other basic dosimetric characteristics of the amorphous silicon detector were investigated, such as linearity with dose rate (+/- 0.5%), repeatability (+/- 2%), and response with variations in gantry rotation and source to detector distance. The latter investigation revealed a significant contribution to the image from optical photon spread in the phosphor layer of the detector. This phenomenon is generally known as "glare," and has been characterized and modeled here as a radially symmetric blurring kernel. This kernel is applied to the calculated dose images as a convolution, and is successfully demonstrated to account for the optical photon spread. This work demonstrates the flexibility and accuracy of the two step algorithm for a high-atomic number detector. The algorithm may be applied to improve performance of dosimetric treatment verification applications, such as direct image comparison, backprojected patient dose calculation, and scatter correction in megavoltage computed tomography. The algorithm allows for dosimetric applications of the new, flat panel portal imager technology in the indirect configuration, taking advantage of a greater than tenfold increase in detector sensitivity over a direct configuration.     

Monte Carlo dose computation for IMRT optimization

A method which combines the accuracy of Monte Carlo dose calculation with a finite size pencil-beam based intensity modulation optimization is presented. The pencil-beam algorithm is employed to compute the fluence element updates for a converging sequence of Monte Carlo dose distributions. The combination is shown to improve results over the pencil-beam based optimization in a lung tumour case and a head and neck case. Inhomogeneity effects like a broader penumbra and dose build-up regions can be compensated for by intensity modulation.     

An investigation of a new amorphous silicon electronic portal imaging device for transit dosimetry

The relationship between the pixel value and exit dose was investigated for a new commercially available amorphous silicon electronic portal imaging device. The pixel to dose mapping function was established to be linear for detector distances between 116.5 cm to 150 cm from the source, radiation field sizes from 5 x 5 cm2 to 20 x 20 cm2 and beam energies of 6 to 18 MV. Coefficients in the mapping function were found to be dependent on beam energy and field size. Open and wedged field profiles measured with the device showed agreement to a maximum of 5% and 8%, respectively, as compared to film. A comparison of relative transmission measurements between the EPID and ion chamber indicate a maximum deviation of 6% and 2% at 6 and 18 MV, respectively, for an attenuator thickness of 21 cm and SDD > or = 130 cm. It was found that accuracies of better than 1% could be obtained if detector position and field size specific fitting parameters were used to generate unique mapping functions for each configuration.     

Patient-specific dosimetry of conventional and intensity modulated radiation therapy using a novel full Monte Carlo phase space reconstruction method from electronic portal images

Electronic portal imagers have promising dosimetric applications in external beam radiation therapy. In this study a patient dose computation algorithm based on Monte Carlo (MC) simulations and on portal images is developed and validated. The patient exit fluence from primary photons is obtained from the portal image after correction for scattered radiation. The scattered radiation at the portal imager and the spectral energy distribution of the primary photons are estimated from MC simulations at the treatment planning stage. The patient exit fluence and the spectral energy distribution of the primary photons are then used to ray-trace the photons from the portal image towards the source through the CT geometry of the patient. Photon weights which reflect the probability of a photon being transmitted are computed during this step. A dedicated MC code is used to transport back these photons from the source through the patient CT geometry to obtain patient dose. Only Compton interactions are considered. This code also produces a reconstructed portal image which is used as a verification tool to ensure that the dose reconstruction is reliable. The dose reconstruction algorithm is compared against MC dose calculation (MCDC) predictions and against measurements in phantom. The reconstructed absolute absorbed doses and the MCDC predictions in homogeneous and heterogeneous phantoms agree within 3% for simple open fields. Comparison with film-measured relative dose distributions for IMRT fields yields agreement within 3 mm, 5%. This novel dose reconstruction algorithm allows for daily patient-specific dosimetry and verification of patient movement.     

Monte Carlo simulation and dosimetric verification of radiotherapy beam modifiers.

Monte Carlo simulation of beam modifiers such as physical wedges and compensating filters has been performed with a rectilinear voxel geometry module. A modified version of the EGS4/DOSXYZ code has been developed for this purpose. The new implementations have been validated against the BEAM Monte Carlo code using its standard component modules (CMs) in several geometrical conditions. No significant disagreements were found within the statistical errors of 0.5% for photons and 2% for electrons. The clinical applicability and flexibility of the new version of the code has been assessed through an extensive verification versus dosimetric data. Both Varian multi-leaf collimator (MLC) wedges and standard wedges have been simulated and compared against experiments for 6MV photon beams and different field sizes. Good agreement was found between calculated and measured depth doses and lateral dose profiles along both wedged and unwedged directions for different depths and focus-to-surface distances. Furthermore, Monte Carlo-generated output factors for both open and wedged fields agreed with linac commissioning beam data within statistical uncertainties of the calculations (<3% at largest depths). Compensating filters of both low-density and high-density materials have also been successfully simulated. As a demonstration, a wax compensating filter with a complex three-dimensional concave and convex geometry has been modelled through a CT scan import. Calculated depth doses and lateral dose profiles for different field sizes agreed well with experiments. The code was used to investigate the performance of a commercial treatment planning system in designing compensators. Dose distributions in a heterogeneous water phantom emulating the head and neck region were calculated with the convolution-superposition method (pencil beam and collapsed cone implementations) and compared against those from the MC code developed herein. The new technique presented in this work is versatile, DICOM-RT compliant and accurate in the simulation of beam modulators. This paper addresses the need to reduce the sources of error in the modelling of beam modifiers since they remain a viable alternative to the MLC technique in the delivery of IMRT beams.     

The effect of the Varian amorphous silicon electronic portal imaging device on exit skin dose

Measurements have been made of the increase in exit surface dose resulting from backscattered radiation generated by the Varian amorphous silicon electronic portal imaging device (EPID). An increase of < or = 14% was demonstrated at both 6 MV and 10 MV, in a manner which suggests that backscatter from the EPID acts to re-establish electronic equilibrium at the exit surface, normally absent in the build-down region. The magnitude of this effect was influenced by field size, measurement depth and exit surface to EPID distance. Assuming typical constraints of portal imaging frequency and geometry, the results suggest that EPID generated backscatter is unlikely to alter the frequency or severity of exit skin reactions. However, the results do suggest that a limit on the minimum separation between the EPID and the exit surface should be set, and that similar investigations should be made for other EPID models.     

Monte Carlo dosimetric study of the BEBIG Co-60 HDR source

Although not as widespread as Ir-192, Co-60 is also available on afterloading equipment devoted to high dose rate brachytherapy, mainly addressed to the treatment of gynaecological lesions. The purpose of this study is to obtain the dosimetric parameters of the Co-60 source used by the BEBIG MultiSource remote afterloader (BEBIG GmbH, Germany) for which there are no dosimetric data available in the literature. The Monte Carlo code GEANT4 has been used to obtain the TG43 parameters and the 2D dose rate table in Cartesian coordinates of the BEBIG Co-60 HDR source. The dose rate constant, radial dose function and anisotropy function have been calculated and are presented in a tabular form as well as a detailed 2D dose rate table in Cartesian coordinates. These dosimetric datasets can be used as input data and to validate the treatment planning system calculations.     

An integrated Monte Carlo dosimetric verification system for radiotherapy treatment planning

An integrated Monte Carlo (MC) dose calculation system, MCRTV (Monte Carlo for radiotherapy treatment plan verification), has been developed for clinical treatment plan verification, especially for routine quality assurance (QA) of intensity-modulated radiotherapy (IMRT) plans. The MCRTV system consists of the EGS4/PRESTA MC codes originally written for particle transport through the accelerator, the multileaf collimator (MLC), and the patient/phantom, which run on a 28-CPU Linux cluster, and the associated software developed for the clinical implementation. MCRTV has an interface with a commercial treatment planning system (TPS) (Eclipse, Varian Medical Systems, Palo Alto, CA, USA) and reads the information needed for MC computation transferred in DICOM-RT format. The key features of MCRTV have been presented in detail in this paper. The phase-space data of our 15 MV photon beam from a Varian Clinac 2300C/D have been developed and several benchmarks have been performed under homogeneous and several inhomogeneous conditions (including water, aluminium, lung and bone media). The MC results agreed with the ionization chamber measurements to within 1% and 2% for homogeneous and inhomogeneous conditions, respectively. The MC calculation for a clinical prostate IMRT treatment plan validated the implementation of the beams and the patient/phantom configuration in MCRTV.     

Use of an amorphous silicon electronic portal imaging device for multileaf collimator quality control and calibration

Multileaf collimator (MLC) calibration and quality control is a time-consuming procedure typically involving the processing, scanning and analysis of films to measure leaf and collimator positions. Faster and more reliable calibration procedures are required for these tasks, especially with the introduction of intensity modulated radiotherapy which requires more frequent checking and finer positional leaf tolerances than previously. A routine quality control (QC) technique to measure MLC leaf bank gain and offset, as well as minor offsets (individual leaf position relative to a reference leaf), using an amorphous silicon electronic portal imaging device (EPID) has been developed. The technique also tests the calibration of the primary and back-up collimators. A detailed comparison between film and EPID measurements has been performed for six linear accelerators (linacs) equipped with MLC and amorphous silicon EPIDs. Measurements of field size from 4 to 24 cm with the EPID were systematically smaller than film measurements over all field sizes by 0.4 mm for leaves/back-up collimators and by 0.2 mm for conventional collimators. This effect is due to the gain calibration correction applied by the EPID, resulting in a 'flattening' of primary beam profiles. Linac dependent systematic differences of up to 0.5 mm in individual leaf/collimator positions were also found between EPID and film measurements due to the difference between the mechanical and radiation axes of rotation. When corrections for these systematic differences were applied, the residual random differences between EPID and film were 0.23 mm and 0.26 mm (1 standard deviation) for field size and individual leaf/back-up collimator position, respectively. Measured gains (over a distance of 220 mm) always agreed within 0.4 mm with a standard deviation of 0.17 mm. Minor offset measurements gave a mean agreement between EPID and film of 0.01+/-0.10 mm (1 standard deviation) after correction for the tilt of the EPID and small rotational misalignments between leaf banks and the back-up collimators used as a reference straight edge. Reproducibility of EPID measurements was found to be very high, with a standard deviation of <0.05 mm for field size and <0.1 mm for individual leaf/collimator positions for a 10x10 cm2 field. A standard set of QC images (three field sizes defined both by leaves only and collimators only) can be acquired in less than 20 min and analysed in 5 min.     

Potential and role of a prototype amorphous silicon array electronic portal imaging device in breathing synchronized radiotherapy

Current electronic portal imaging devices (EPID) are limited in their ability to provide direct and quick verification and monitoring of patients during both setup and treatment of breathing synchronized radiotherapy (BSRT, including breathing gated, voluntary and forced breath-hold radiotherapy treatment.) These limitations are largely due to their slow image capture rate and poor image quality. An amorphous silicon array flat panel electronic portal imaging device (si-EPID) is emerging to meet the challenge. The purpose of this study is threefold: (1) to characterize the performance of a prototype si-EPID; (2) to compare image quality against that of digitized films; and (3) to evaluate the device in terms of verification of patient setup and monitoring during BSRT. In this study a Varian prototype si-EPID detector array and Clinic accelerator at the University of California Davis Cancer Center were used for imaging. Three quality assurance phantoms: a Lutz PVC phantom, a modified "Las Vegas" phantom, and a RMI model 1151 phantom, were used to characterize the imaging system. A Rando head phantom was used for anthropomorphic imaging tests. Images were obtained with the si-EPID and a Fuji RX film in a Kodak X-Omatic cassette. To investigate the clinical application, two sets of si-EPID images were collected from a lung cancer patient during a 22 s breath-hold and normal breathing. The quality of images obtained with the fast mode was found to be comparable to that obtained with the digitized films. The images with the standard mode were found to be better than the digitized film images. With this prototype si-EPID, it is possible to collect the images at the beginning, middle, and end of each breath-hold for those patients who can hold their breath for longer than 15 s. The si-EPID images can provide a quick verification of the initial patient setup and subsequent treatment position throughout the daily fractionation.     

Dosimetric properties of an amorphous silicon electronic portal imaging device for verification of dynamic intensity modulated radiation therapy

Dosimetric properties of an amorphous silicon electronic portal imaging device (EPID) for verification of dynamic intensity modulated radiation therapy (IMRT) delivery were investigated. The EPID was utilized with continuous frame-averaging during the beam delivery. Properties studied included effect of buildup, dose linearity, field size response, sampling of rapid multileaf collimator (MLC) leaf speeds, response to dose-rate fluctuations, memory effect, and reproducibility. The dependence of response on EPID calibration and a dead time in image frame acquisition occurring every 64 frames were measured. EPID measurements were also compared to ion chamber and film for open and wedged static fields and IMRT fields. The EPID was linear with dose and dose rate, and response to MLC leaf speeds up to 2.5 cm s(-1) was found to be linear. A field size dependent response of up to 5% relative to dmax ion-chamber measurement was found. Reproducibility was within 0.8% (1 standard deviation) for an IMRT delivery recorded at intervals over a period of one month. The dead time in frame acquisition resulted in errors in the EPID that increased with leaf speed and were over 20% for a 1 cm leaf gap moving at 1.0 cm s(-1). The EPID measurements were also found to depend on the input beam profile utilized for EPID flood-field calibration. The EPID shows promise as a device for verification of IMRT, the major limitation currently being due to dead-time in frame acquisition.     

Amorphous silicon EPID calibration for dosimetric applications: comparison of a method based on Monte Carlo prediction of response with existing techniques

For EPID dosimetry, the calibration should ensure that all pixels have a similar response to a given irradiation. A calibration method (MC), using an analytical fit of a Monte Carlo simulated flood field EPID image to correct for the flood field image pixel intensity shape, was proposed. It was compared with the standard flood field calibration (FF), with the use of a water slab placed in the beam to flatten the flood field (WS) and with a multiple field calibration where the EPID was irradiated with a fixed 10x10 field for 16 different positions (MF). The EPID was used in its normal configuration (clinical setup) and with an additional 3 mm copper slab (modified setup). Beam asymmetry measured with a diode array was taken into account in MC and WS methods. For both setups, the MC method provided pixel sensitivity values within 3% of those obtained with the MF and WS methods (mean difference<1%, standard deviation<2%). The difference of pixel sensitivity between MC and FF methods was up to 12.2% (clinical setup) and 11.8% (modified setup). MC calibration provided images of open fields (5x5 to 20x20 cm2) and IMRT fields to within 3% of that obtained with WS and MF calibrations while differences with images calibrated with the FF method for fields larger than 10x10 cm2 were up to 8%. MC, WS and MF methods all provided a major improvement on the FF method. Advantages and drawbacks of each method were reviewed.     

Photon scatter in portal images: physical characteristics of pencil beam kernels generated using the EGS Monte Carlo code

Pencil beam kernels describing scattered photon fluence behind homogeneous water slabs at various air gap distances were generated using the EGS Monte Carlo code. Photon scatter fluence was scored in separate bins based on the particle's history: singly scattered, multiply scattered, and bremsstrahlung and positron annihilation photons. Simultaneously, the mean energy and mean angle with respect to the incident photon pencil beam were tallied. Kernels were generated for incident photon pencil beams exhibiting monoenergetic spectra of 2.0 and 10.0 MeV, and polyenergetic spectra representative of 6 and 24 MV beams. Reciprocity was used to generate scatter fractions on the central axis for various field sizes, phantom thicknesses, and air gaps. The scatter kernels were further characterized by full width at half-maximum estimates. Modulation transfer functions were calculated, providing theoretical estimates of the limit of performance of portal imaging systems due to the intrinsic scattering of photon radiation through the patient.