Spatially fractionated radiotherapy (SFRT) targeting the hypoxic tumor segment for the intentional induction of non-targeted effects: An in silico study to exploit a new treatment paradigm.

INTRODUCTION: The possibility of intentionally triggering non targeted effects (NTEs) using spatially fractionated radiotherapy (SFRT) alone or combined with immunotherapy is an intriguing and fascinating area of research. Among different techniques for SFRT, stereotactic body radiotherapy targeting exclusively the central hypoxic segment of bulky tumors, (SBRT-PATHY) might trigger immunogenic cell death more efficiently. This in silico study aims to identify the best possible dosimetric trade-off for prescribing SFRT with volumetric modulated arc (VMAT) based stereotactic radiotherapy (SRT). MATERIAL AND METHODS: Eight spherical volumes defined "Gross Tumor Volumes" (GTVs) were generated with diameters of 3-10 cm (with incremental steps of 1 cm), simulating tumor lesions. The inner third part of each GTV (GTVcentral) was selected to simulate the central hypoxic area and a ring structure was derived around it to simulate the tumor periphery (GTVperipheral). Volumetric modulated arc radiation treatment (VMAT) plans were calculated to deliver a single fraction of 10 Gy to each GTVcentral with different dose prescription methods: target mean and isodose driven (40, 50, 60, 70, 80 and 90%).The volume of GTVperipheral receiving less than 2 Gy was recorded as dosimetric performance indicator. RESULTS: 56 possible dosimetric scenarios were analyzed. The largest percentage of GTVperipheral spared from the dose of 2 Gy was achieved with dose prescription methods to the 70% isodose line for lesions smaller than 6 cm (range 42.9-48.4%) and to the target mean for larger ones (range 52.9-64.5%). CONCLUSIONS: Optimizing the dose prescription method may reduce the dose to tumor periphery in VMAT-based SFRT, thus potentially sparing tumor infiltrating immune cells. The optimal method may vary according to the size of the lesion. This should be taken into account when designing prospective trials using SFRT.

A validation study of a dedicated software for an automated in vivo dosimetry control in radiotherapy.

In vivo dosimetry (IVD) is the last step of a radiotherapy quality control program aimed to ensure that the dose delivered is in agreement with that prescribed. IVD procedures based on single detectors are time-consuming and impossible to use for the modern radiotherapy techniques, based on static or kinetic beams (modulated in intensity fluence); this means that more efficient and practical methods are highly recommended. The practical method SOFTDISO, based on the use of electronic portal image device (EPID), provides two tests (i) the R ratio between the reconstructed and the planned isocenter doses to verify an agreement within 5% and (ii) the γ-analysis of the EPID images, to verify γ% ≥ 90% and γmean ≤ 0.4. This paper reports the results of 11,357 IVD tests carried out for 823 patients treated by three-dimensional conformal radiation therapy and volumetric modulated arc therapy techniques. In particular, the dose disagreements are reported distinguishing two kinds of causes, those of (i) class 1 that includes the errors due to inadequate quality controls and (ii) the class 2, due to patient morphological changes. About the tests out of tolerance, 6% were by VMAT and 21% by 3DCRT, but taking into account the only class 1 of errors, i.e., removing the causes of class 2, only 7% of patients examined presented at least one of the three mean indexes out of tolerance. The workload for IVD on 9 patients/day per linac is about 52 min/day but recently, a new automated SOFTDISO version has been implemented to reduce the time to about 34 min/day.

Step-and-Shoot IMRT by Siemens Beams: An EPID Dosimetry Verification During Treatment.

PURPOSE: This work reports the extension of a semiempirical method based on the correlation ratios to convert electronic portal imaging devices transit signals into in vivo doses for the step-and-shoot intensity-modulated radiotherapy Siemens beams. The dose reconstructed at the isocenter point Diso, compared to the planned dose, Diso,TPS, and a γ-analysis between 2-dimensional electronic portal imaging device images obtained day to day, seems to supply a practical method to verify the beam delivery reproducibility. METHOD: The electronic portal imaging device images were obtained by the superposition of many segment fields, and the algorithm for the Diso reconstruction for intensity-modulated radiotherapy step and shoot was formulated using a set of simulated intensity-modulated radiotherapy beams. Moreover, the in vivo dose-dedicated software was integrated with the record and verify system of the centers. RESULTS: Three radiotherapy centers applied the in vivo dose procedure at 30 clinical intensity-modulated radiotherapy treatments, each one obtained with 5 or 7 beams, and planned for patients undergoing radiotherapy for prostatic tumors. Each treatment beam was checked 5 times, obtaining 900 tests of the ratios R = Diso/Diso,TPS. The average R value was equal to 1.002 ± 0.056 (2 standard deviation), while the mean R value for each patient was well within 5%, once the causes of errors were removed. The γ-analysis of the electronic portal imaging device images, with 3% 3 mm acceptance criteria, showed 90% of the tests with Pγ < 1 ≥ 95% and γmean ≤ 0.5. The off-tolerance tests were found due to incorrect setup or presence of morphological changes. This preliminary experience shows the great utility of obtaining the in vivo dose results in quasi real time and close to the linac, where the radiotherapy staff may immediately spot possible causes of errors. The in vivo dose procedure presented here is one of the objectives of a project, for the development of practical in vivo dose procedures, financially supported by the Istituto Nazionale di Fisica Nucleare.

Quasi real time in vivo dosimetry for VMAT.

PURPOSE: Results about the feasibility of a method for quasi real time in vivo dosimetry (IVD) at the isocenter point for volumetric modulated arc therapy (VMAT) are here reported. The method is based on correlations between the EPID signal and the dose on the beam central axis. Moreover, the γ-analysis of EPID images was adopted to verify off-axis reproducibility of fractionated plan delivery. METHODS: An algorithm to reconstruct in vivo the isocenter dose, D(iso), for RapidArc treatments has been developed. 20 VMAT plans, optimized with two opposite arcs, for prostate, pancreas, and head treatments have been delivered by a Varian linac both to a conic PMMA phantom with elliptical section and to patients. The ratios R between reconstructed D(iso) and the planned doses were determined for phantom and patient irradiations adopting an acceptance criterion of ±5%. In total, 40 phantom checks and 400 patient checks were analyzed. Moreover, 3% and 3 mm criteria were adopted for portal image γ-analysis to assess patient irradiation reproducibility. RESULTS: The average ratio R, between reconstructed and planned doses for the PMMA phantom irradiations was equal to 1.007 ± 0.024. When the IVD method was applied to the 20 patients, the average R ratio was equal to 1.003 ± 0.017 and 96% of the tests were within the acceptance criteria. The portal image γ-analysis supplied 88% of the tests within the pass rates γ(mean) ≤ 0.4 and P(γ<1) ≥ 98%. All the warnings were understood comparing the CT and the cone beam CT images and in one case a patient's setup error was detected and corrected for the successive fractions. CONCLUSIONS: This preliminary experience suggests that the method is able to detect dosimetric errors in quasi real time at the end of the therapy session. The authors intend to extend this procedure to other pathologies with the integration of in-room imaging verification by cone beam CT.

An in-vivo dosimetry procedure for Elekta step and shoot IMRT.

PURPOSE: The aim of this work was to extend an in-vivo dosimetry (IVD) method, previously developed by the authors for 3D-conformal radiotherapy, to step and shoot IMRT treatments for pelvic tumors delivered by Elekta linacs. MATERIALS AND METHODS: The algorithm is based on correlation functions to convert EPID transit signals into in-vivo dose values at the isocenter point, Diso. The EPID images were obtained by the so-called "IMRT Dosimetric Weighting" mode as a superposition of many segment fields. This way each integral dosimetric image could be acquired in about 10 s after the end of beam delivery and could be processed while delivering the successive IMRT beams. A specific algorithm for Diso reconstruction especially featured for step and shoot IMRT was implemented using a fluence inhomogeneity index, FI, introduced to describe the degree of beam modulation with respect to open beams. A γ-analysis of 2D-EPID images obtained day to day, resulted rapid enough to verify the plan delivery reproducibility. RESULTS: Fifty clinical IMRT beams, planned for patients undergoing radiotherapy of pelvic tumors, were used to irradiate a homogeneous phantom. For each beam the agreement between the reconstructed dose, Diso, and the TPS computed dose, Diso,TPS, was well within 5%, while the mean ratio R = Diso/Diso,TPS resulted for 250 tests equal to 1.006 ± 0.036. The same beams were checked in vivo, i.e. during patient treatment delivery, obtaining 500 tests whose average R ratio resulted equal to 1.011 ± 0.042. The γ-analysis of the EPID images with 5% 3 mm criteria supplied 85% of the tests with pass rates γ(mean) ≤ 0.5 and P(γ<1) ≥ 90%.

Real-time dose reconstruction for wedged photon beams: a generalized procedure.

A practical and accurate generalized procedure to reconstruct the isocenter dose D(iso) for 3D conformal radiotherapy (3DCRT) has been developed for X-ray open beams supplied by linacs of different manufacturers and equipped with aSi electronic portal imaging devices (aSi EPIDs). This paper reports an extension of the method, to be applied at the wedged X-ray beams characterized by the wedge attenuation factor W(AF). Using water-equivalent solid phantoms (SPs) of different thicknesses, w, and photon square fields of sizes, L, the generalized midplane doses D(0)(W(AF), w/2,L) and generalized transit signals s(t)(0)(W(AF),w,L) by 38 beams of six different linacs were determined. The generalized data were fitted by surface equations and used together with the information of the 'record & verify' network of the centers. In this manner, for every beam, the D(iso) reconstruction was obtained in about 25 seconds after the treatment. To test the in vivo dosimetric procedure, six pelvic treatments that used conformed wedged beams were carried out with three linacs of different manufacturers. For every beam, the comparison between the reconstructed D(iso) and the D(iso,TPS) computed by the TPS, resulted in an acceptable tolerance level of ±5%, estimated for this kind of treatment. Generally the in vivo dosimetry methods that use EPIDs require: (i) a special effort for the dosimetric commissioning with SPs of different thicknesses, and (ii) extra time for the analysis of the EPID signals. The proposed procedure simplifies the commissioning step and supplies for Varian, Elekta, and Siemens linacs equipped with the aSi EPIDs a quasi-real time in vivo dosimetry for open and wedged 3DCRT fields.

Calibration of Elekta aSi EPIDs used as transit dosimeter.

The transit in vivo dosimetry performed by the Electronic Portal Imaging Device (EPID), avoids the problem of solid-state detector positioning on the patient. Moreover, the dosimetric characterization of the recent Elekta aSi EPIDs in terms of signal stability and linearity enables these detectors adaptable for the transit in vivo dosimetry with 6, 10 and 15 MV photon beams. However, the implementation of the EPID transit dosimetry requires several measurements. Recently, the present authors have developed an in vivo dosimetry method for the 3D CRT based on correlation functions defined by the ratios between the transit signal, s(t) (w,L), by the EPID and the phantom mid-plane dose, D(m)(w,L), at the Source to Axis Distance (SAD) as a function of the phantom thickness, w, and the square field dimensions, L. When the phantom mid-plane was positioned at distance d from the SAD, the ratios st(w,L)/s't(d,w,L), were used to take into account the variation of the scattered photon contributions on the EPID as a function of, d and L. The aim of this paper was the implementation of a procedure that uses generalized correlation functions obtained by nine Elekta Precise linac beams. The procedure can be used by other Elekta Precise linacs equipped with the same aSi EPIDs assuring the stabilities of the beam output factors and the EPID signals. The calibration procedure of the aSi EPID here reported avoids measurements in solid water equivalent phantoms needed to implement the in vivo dosimetry method in the radiotherapy center. A tolerance level ranging between ±5% and ±6% (depending on the type of tumor) was estimated for the comparison between the reconstructed isocenter dose, D(iso) and the computed dose D(iso,TPS) by the treatment planning system (TPS).

Portal dose measurements by a 2D array.

A 2D array (PTW, type 10024), equipped with 729 vented plane parallel ion-chambers, has been calibrated as a detector for the in vivo comparison between measured and predicted portal doses for head-neck tumors. The comparison of absolute portal doses measured to ones predicted by a commercial treatment planning system within the field of view of the CT scanner, can help the delivered dose verification during different treatment fractions, in particular when the patient's present weight loss. This paper reports the preliminary results of the comparison of the portal doses measured by a PTW 2D array during several radiotherapy fractions and the predicted portal doses for seven patients undergoing head-neck tumor radiotherapy. The gamma index analysis supplied an agreement of more than 95% of the dose-point P(gamma)>95% within acceptance criteria, in terms of dose difference, DeltaD(max), and distance-agreement, Deltad(max), equal to 5% and 4mm, respectively. After the third week, one patient showed a decrease of P(gamma) values due to the markedly reduced patient's thickness. Even if the spatial resolution of the 2D array was 1cm, there were two advantages in the use of this 2D array as a portal dose device for IMRT quality control. The first one was the use of a stable and efficient absolute dosimeter for in vivo verification, although its construction and behavior for other gantry angles need to be tested, and the second one was the time efficiency in verifying the correct dose delivery in several fractions of the therapy. This study presents acceptance criteria for the comparison of TPS-predicted portal dose images with in vivo 2D ion-chamber measurements for IMRT. In particular, portal dose measurements offer clues for additional studies as to which indicators can signal the need for replanning during treatment.

A preliminary dosimetric characterization of chemical vapor deposition diamond detector prototypes in photon and electron radiotherapy beams.

Three radiation detectors based on polycrystalline diamond films with different thickness and resistivity, obtained by microwave chemical vapor deposition, were tested to assess their suitability for relative dosimetry of photon and electron beams supplied by clinical linear accelerators. All samples showed a linear response as a function of the absorbed dose. The sensitivity per unit of detector sensitive volume spanned between 7 and 43 nC Gy(-1) mm(-3) with an applied electric field of 40 kV/cm. The dose rate dependence was evaluated following the Fowler theory and delta coefficient values between 0.95 and 1.00 were found for the three samples when polarized at 40 kV/cm. Percentage depth dose curves, output factors, and normalized dose profiles were determined for 6 and 10 MV photon beams and for 6 and 15 MeV electron beams. The results obtained with the diamond detectors were in good agreement with those obtained by reference detector measurements [all the data were within the experimental uncertainty of 1% (1sigma)].

Experimental dosimetry of a 32P catheter-based endovascular brachytherapy source.

The experimental dosimetry in a water phantom of a 32P linear source, 20 mm in length, used for the brachytherapy of coronary vessels is reported. The source content activity, A, was determined by means of a calibrated well ion-chamber and the value was compared with the contained activity reported in the manufacturer's certification. In this field of brachytherapy dosimetry, radiochromic film supplies a high enough spatial resolution. A highly sensitive radiochromic film, that presents only one active layer, was used in this work for the source dosimetry in a water phantom. The radiochromic film was characterized by electron beams produced by a clinical linac. A Monte Carlo calculation of beta spectra in water at different distances along the source transverse bisector axis allowed to take into account the low dependence of film response from the electron beam energy. The adopted experimental set-up, with the source in its catheter positioned on the film plane inside the water phantom, supplies accurate dosimetric information. The measured dose rate to water per unit of source activity at reference distance, D(r0, theta0)/A, in units of cGy s(-1) GBq(-1), was in agreement with the value reported in the manufacturer's certification within the experimental uncertainty. The radial dose function, g(r), is in good agreement with the literature data. The anisotropy function F(r, theta) is also reported. The analysis of the dose profile obtained at 2 mm from the source longitudinal axis shows that the uniformity is within 10% along 75% of the 20 mm treatment length. The adopted experimental set-up seems to be adequate for the quality control procedure of the dose homogeneity distribution in the water medium.

The wall correction factor for a spherical ionization chamber used in brachytherapy source calibration.

The effect of wall chamber attenuation and scattering is one of the most important corrections that must be determined when the linear interpolation method between two calibration factors of an ionization chamber is used. For spherical ionization chambers the corresponding correction factors A(w) have to be determined by a non-linear trend of the response as a function of the wall thickness. The Monte Carlo and experimental data here reported show that the A(w) factors obtained for an Exradin A4 chamber, used in the brachytherapy source calibration, in terms of reference air kerma rate, are up to 1.2% greater than the values obtained by the linear extrapolation method for the studied beam qualities. Using the Aw factors derived from Monte Carlo calculations, the accuracy of the calibration factor N(K,Ir) for the Exradin A4, obtained by the interpolation between two calibration factors, improves about 0.6%. The discrepancy between the new calculated factor and that obtained using the complete calibration curve of the ion-chamber and the 192Ir spectrum is only 0.1%.

A standard dosimetry procedure for 192Ir sources used for endovascular brachytherapy.

The experimental dosimetry of a high dose rate (HDR) 192Ir source used for the brachytherapy of peripheral vessels is reported. The direct determination of the reference air kerma rate Kr agrees, within the experimental uncertainty, with the results obtained by a well ionization chamber calibrated at the NIST and the manufacturer's certification. A highly sensitive (HS) radiochromic film (RCF), that presents only one active layer, was used for the source dosimetry in a water phantom. The adopted experimental set-up, with the source in its catheter positioned on the RCF plane, seems to have given better accuracy of the RCF optical density measurements. The agreement between the measurement of the dose rate constant DKr (10 mm, pi/2) and the literature data confirmed the coherence of the HS RCF calibration obtained by the kerma in air measurements. The RCF measurements supplied dosimetric information about the dose to water per reference air kerma rate D(r, theta)/Kr along the source transverse bisector axis, the radial dose function g(r) and the anisotropy function F(r, theta). The value D(2 mm, pi/2)/Kr = 22.4 +/- 1.2 cGy h(-1)/(microGy h(-1)) is supplied with a dose uncertainty that is essentially due to the indeterminacy of the source position in the catheter. The data of the radial and anisotropy functions have been compared with Monte Carlo determinations reported in the literature.

A correction method for diamond detector signal dependence with proton energy.

Small dosimeters as solid state detectors can be useful for the dosimetric characterization and periodic quality control of radiotherapy proton beams. The calibration of solid state detectors for proton beams is not a solved problem especially for ophthalmologic proton beams, where these detectors present a LET-dependent signal. In this work a PTW diamond detector has been selected because of its good signal reproducibility (0.3%) and stable response with accumulated dose. A method that takes into account the LET dependence of the diamond detector signal, at 62 MeV proton beam, is here proposed. In particular an empirical correction factor, kDD(Eo) (Rres), has been determined as a function of the residual range quality index, to correct the diamond detector signal for a proton beam of incident effective energy E0= 62 MeV. A dedicated software allows us to use the diamond detector as an on-line reference dosimeter, where an ionization chamber may be difficult to use, or for periodic quality control procedures. The article also reports a comparison between the signal dependence on proton energy of silicon, diamond, and radiochromic film detectors.

Characteristics of silicon and diamond detectors in a 60 MeV proton beam.

The calibration factor variation for a PTW natural diamond detector and a Scanditronix p-type stereotactic silicon diode (designed for use in photon beams) was studied in the 10-59 MeV range. Irradiations were performed in a water phantom with the 60 MeV ocular therapy beam at the CCO (UK). The diamond detector showed a sensitivity increase with energy, underestimating the dose by about 18% at the Bragg peak, by 7% at the centre and by 17% at the distal end of the SOBP region. The silicon diode did not show any significant sensitivity change with energy. However, a decrease in response of 24% was observed for an accumulated dose of 300 Gy.

PTW-diamond detector: dose rate and particle type dependence.

In this paper the suitability of a PTW natural diamond detector (DD) for relative and reference dosimetry of photon and electron beams, with dose per pulse between 0.068 mGy and 0.472 mGy, was studied and the results were compared with those obtained by a stereotactic silicon detector (SFD). The results show that, in the range of the examined dose per pulse the DD sensitivity changes up to 1.8% while the SFD sensitivity changes up to 4.5%. The fitting parameter, delta, used to correct the dose per pulse dependence of solid state detectors, was delta = 0.993 +/- 0.002 and delta = 1.025 +/- 0.002 for the diamond detector and for the silicon diode, respectively. The delta values were found to be independent of particle type of two conventional beams (a 10 MV x-ray beam and a 21 MeV electron beam). So if delta is determined for a radiotherapy beam, it can be used to correct relative dosimetry for other conventional radiotherapy beams. Moreover the diamond detector shows a calibration factor which is independent of beam quality and particle type, so an empirical dosimetric formalism is proposed here to obtain the reference dosimetry. This formalism is based on a dose-to-water calibration factor and on an empirical coefficient, that takes into account the reading dependence on the dose per pulse.

Dosimetric characterization of silicon and diamond detectors in low-energy proton beams.

The dosimetric behaviour of a Scanditronix p-type silicon diode and a PTW natural diamond detector was studied in low-energy proton beams in the 8.3-21.5 MeV range. The properties investigated were linearity, reproducibility, dose rate dependence, energy and linear energy transfer (LET) dependence. The influence of detector thickness on the results of depth dose measurements was also demonstrated. A Markus parallel plate ionization chamber was used for reference dosimetry. Silicon diode and diamond detectors showed linearity at therapeutic dose level, reproducibility better than 1% (1sigma) and sensitivity variation with dose rate and proton energy.

Radiochromic film dosimetry of a low energy proton beam.

In this work some dosimetric characteristics of MD-55-2 GafChromic films were studied in a low energy proton beam (21.5 MeV) directly in a water phantom. The nonlinearity of the optical density was quantified by a factor P(lin). A correction factor P(en), that accounts for optical density dependence on the energy, was empirically determined. The effects of detector thickness in depth dose measurements and of the film orientation with respect to beam direction were investigated. The results show that the MD-55-2 films provide dose measurements with the films positioned perpendicularly to the proton beam. A dosimetric formalizm is proposed to determine the dose to water at depth d, with films oriented perpendicularly to the beam axis. This formalism uses a calibration factor of the radiochromic film determined directly on the proton beam at a reference depth in water, and the P(lin) factor, that takes into account the nonlinearity of the calibration curve and the P(en) factor that, in turn takes into account the change of proton beam energy in water. The MD-55-2 films with their high spatial resolution and the quasiwater equivalent material are attractive, positioned perpendicularly along the beam axis, for the absolute dose determination of very small beam sizes and modulated proton beams.

The saturation loss for plane parallel ionization chambers at high dose per pulse values.

The use of plane parallel ionization chambers with electron beams with high dose per pulse entails dose uncertainties due to the overestimation of the ion recombination factor, k, up to 20% if conventional dosimetric protocols are used. In this work MD-55-2 radiochromic films have been used as reference dosimeters to obtain dose to water per pulse DGAF(w) values for three Novac7 (Hitesys) electron beams of E0 = 5.8 MeV. However, the beam calibration by MD-55-2 films is time consuming and the use of plane parallel chambers is fundamental for a periodic quality control procedure. Three plane parallel chambers have been used and the general formula for the k determination has been tested using the calibration doses, DGAF(w). In particular, consistent ion recombination factors ksat(V0) (with the ion chamber polarized at V0), that follow the Boag theory, have been estimated at different dose per pulse values for the three plane parallel ionization chambers. This means that at present any ion chamber needs a specific ksat (V0) determination by using a reference dosimeter for which the response is independent of the dose rate. An accurate determination of ksat(V0), using a reference quality beam, can be used to determine the dose to water per pulse for electron beams of different quality and geometrical configuration.

Alanine/EPR dosimetry in brachytherapy.

The paper reports the experimental procedure adopted to determine the absorbed dose rate in water per reference air kerma rate, D(Kr), (d, theta), along the transverse bisector axis of a 137Cs brachytherapy source. The dose rate measurements have been carried out at difference distances, d, from the source using alanine dosemeters in a water phantom. The reference air kerma rate, Kr, was determined adopting a 'direct procedure' that uses a spherical ionization chamber in air. The dose rate constant of the source examined was D(Kr) (1, pi/2) = 0.99 +/- 0.03 cGy h(-1) (microGy h(-1))(-1). The values of the radial dose function along the transverse axis, g(d), determined with an uncertainty of 3.4% (1sigma), were found to be in good agreement with the results reported in the literature. The uncertainty in dose rate value has been estimated as 2.8% (1sigma) for distances from the source up to 7 cm. Kr has been determined with 1.2% (1sigma) uncertainty. So D(r) (d, pi/2) values were determined with 3% (1sigma) uncertainty.

Two quality control procedures on radiotherapy beam calibration and treatment planning system implementation.

New challenging dosimetric approaches, such as narrow beams and 3D algorithms, are being used in radiotherapy. In this paper two quality control (QC) procedures are reported. The first one concerns the QC of the dosimetry of small x-ray beams, generally carried out by using silicon detectors. The comparison of dose values obtained by a silicon diode, a diamond detector, and radiochromic films shows that for x-ray beams of high energy, the silicon diode can give an overestimation of the output factors in phantom, up to 4%. This is due to the higher than unit density silicon diode and the surrounding envelope that restore the lateral electron equilibrium. About the 3D algorithms for breast treatment planning, a quality control test has been adopted to verify the accuracy of the computed dosimetry when "loss of scatter" occurs. The results show a sensible agreement (within 1.5%) between computed and experimental data.