Analytic representation of the backscatter correction factor at the exit of high energy photon beams

A dual two dimensional electronic portal imaging device transit dosimetry model based on an empirical quadratic formalism

OBJECTIVE

This study describes a two dimensional electronic portal imaging device (EPID) transit dosimetry model that can predict either: (1) in-phantom exit dose, or (2) EPID transit dose, for treatment verification.

METHODS

The model was based on a quadratic equation that relates the reduction in intensity to the equivalent path length (EPL) of the attenuator. In this study, two sets of quadratic equation coefficients were derived from calibration dose planes measured with EPID and ionization chamber in water under reference conditions. With two sets of coefficients, EPL can be calculated from either EPID or treatment planning system (TPS) dose planes. Consequently, either the in-phantom exit dose or the EPID transit dose can be predicted from the EPL. The model was tested with two open, five wedge and seven sliding window prostate and head and neck intensity-modulated radiation therapy (IMRT) fields on phantoms. Results were analysed using absolute gamma analysis (3%/3 mm).

RESULTS

The open fields gamma pass rates were >96.8% for all comparisons. For wedge and IMRT fields, comparisons between predicted and TPS-computed in-phantom exit dose resulted in mean gamma pass rate of 97.4% (range, 92.3-100%). As for the comparisons between predicted and measured EPID transit dose, the mean gamma pass rate was 97.5% (range, 92.6-100%).

CONCLUSION

An EPID transit dosimetry model that can predict in-phantom exit dose and EPID transit dose was described and proven to be valid.

ADVANCES IN KNOWLEDGE

The described model is practical, generic and flexible to encourage widespread implementation of EPID dosimetry for the improvement of patients' safety in radiotherapy.

Evaluation of six TPS algorithms in computing entrance and exit doses

Entrance and exit doses are commonly measured in in vivo dosimetry for comparison with expected values, usually generated by the treatment planning system (TPS), to verify accuracy of treatment delivery. This report aims to evaluate the accuracy of six TPS algorithms in computing entrance and exit doses for a 6 MV beam. The algorithms tested were: pencil beam convolution (Eclipse PBC), analytical anisotropic algorithm (Eclipse AAA), AcurosXB (Eclipse AXB), FFT convolution (XiO Convolution), multigrid superposition (XiO Superposition), and Monte Carlo photon (Monaco MC). Measurements with ionization chamber (IC) and diode detector in water phantoms were used as a reference. Comparisons were done in terms of central axis point dose, 1D relative profiles, and 2D absolute gamma analysis. Entrance doses computed by all TPS algorithms agreed to within 2% of the measured values. Exit doses computed by XiO Convolution, XiO Superposition, Eclipse AXB, and Monaco MC agreed with the IC measured doses to within 2%‐3%. Meanwhile, Eclipse PBC and Eclipse AAA computed exit doses were higher than the IC measured doses by up to 5.3% and 4.8%, respectively. Both algorithms assume that full backscatter exists even at the exit level, leading to an overestimation of exit doses. Despite good agreements at the central axis for Eclipse AXB and Monaco MC, 1D relative comparisons showed profiles mismatched at depths beyond 11.5 cm. Overall, the 2D absolute gamma (3%/3 mm) pass rates were better for Monaco MC, while Eclipse AXB failed mostly at the outer 20% of the field area. The findings of this study serve as a useful baseline for the implementation of entrance and exit in vivo dosimetry in clinical departments utilizing any of these six common TPS algorithms for reference comparison.

PACS numbers: 87.55.‐x, 87.55.D‐, 87.55.N‐, 87.53.Bn

Backscatter correction factor for megavoltage photon beam

PURPOSE

For routine clinical dosimetry of photon beams, it is often necessary to know the minimum thickness of backscatter phantom material to ensure that full backscatter condition exists.

METHODS

In case of insufficient backscatter thickness, one can determine the backscatter correction factor, BCF(s,d,t), defined as the ratio of absorbed dose measured on the central-axis of a phantom with backscatter thickness of t to that with full backscatter for square field size s and forward depth d. Measurements were performed in SAD geometry for 6 and 15 MV photon beams using a 0.125 cc thimble chamber for field sizes between 10 × 10 and 30 × 30 cm at depths between d(max) (1.5 cm for 6 MV and 3 cm for 15 MV) and 20 cm.

RESULTS

A convolution method was used to calculate BCF using Monte-Carlo simulated point-spread kernels generated for clinical photon beams for energies between Co-60 and 24 MV. The convolution calculation agrees with the experimental measurements to within 0.8% with the same physical trend. The value of BCF deviates more from 1 for lower energies and larger field sizes. According to our convolution calculation, the minimum BCF occurs at forward depth d(max) and 40 × 40 cm field size, 0.970 for 6 MV and 0.983 for 15 MV.

CONCLUSIONS

The authors concluded that backscatter thickness is 6.0 cm for 6 MV and 4.0 cm for 15 MV for field size up to 10 × 10 cm when BCF = 0.998. If 4 cm backscatter thickness is used, BCF is 0.997 and 0.983 for field size of 10 × 10 and 40 × 40 cm for 6 MV, and is 0.998 and 0.990 for 10 × 10 and 40 × 40 cm for 15 MV, respectively.

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.

Transmission measurements in air using the ESTRO mini-phantom

The aim of this work is to study the possibility of using the ESTRO mini-phantom for transmission measurements of primary kerma in water at a point free in air. We discuss in-air measurements in general, with special attention given to in-air equivalent measurements using a water equivalent mini-phantom. The study includes four different photon energies (4, 6, 10 and 18 MV), where scoring of dose and primary kerma inside a mini-phantom in narrow beam geometry is performed with the Monte Carlo code EGS4. The results reveal that relative measurements (i.e. with and without a water absorber present) at 10 cm depth in a mini-phantom do not represent the variation of primary kerma in water at a point free in air (deviations as large as 7% at 4 MV are observed). Minimum deviations are obtained at depths somewhat larger than the depth of dose maximum. The observed deviations are due to a considerable beam hardening in the water absorber, which changes the amount of attenuation and scatter inside the mini-phantom.

Electron contamination and build-up doses in conformal radiotherapy fields

The dose in the build-up region depends upon the primary photon beam, backscattered radiation from the patient and contamination radiation from outside the patient. In this paper, a model based on measured data is proposed which allows the build-up dose for arbitrarily shaped treatment fields to be determined. The dose in the build-up region is assumed to comprise a primary photon component and a contamination component that is a function of the field size and shape. This contamination component, for modelling purposes, is subdivided into contributions that correspond to elements of 1 cm by 1 cm cross-sectional area at the plane of the isocentre. The magnitude of these components has been obtained by fitting measured data to an exponential function. The exponent was found to vary linearly with depth for energies between 4 MV and 20 MV. The coefficient decreased linearly with depth at 4, 6 and 8 MV, but exhibited a broad build-up region at 20 MV. The primary component, in the build-up region, could be approximated by a 100 - (100 - PSD) e(-mu d) function, where PSD is the primary surface dose. The values obtained during the fitting procedure were used to calculate dose in the build-up region for arbitrarily shaped fields. Good agreement was found in each case.

Backscatter correction algorithm for TBI treatment conditions

The accuracy requirements in target dose delivery is, according to ICRU, +/- 5%. This is so not only in standard radiotherapy but also in total body irradiation (TBI). Physical dosimetry plays an important role in achieving this recommended level. The semi-infinite phantoms, customarily used for dosimetry purposes, give scatter conditions different to those of the finite thickness of the patient. So dose calculated in patient's points close to beam exit surface may be overestimated. It is then necessary to quantify the backscatter factor in order to decrease the uncertainty in this dose calculation. The backward scatter has been well studied at standard distances. The present work intends to evaluate the backscatter phenomenon under our particular TBI treatment conditions. As a consequence of this study, a semi-empirical expression has been derived to calculate (within 0.3% uncertainty) the backscatter factor. This factor depends lineally on the depth and exponentially on the underlying tissue. Differences found in the qualitative behavior with respect to standard distances are due to scatter in the bunker wall close to the measurement point.