# Evaluation of Dose Calculation Algorithms of Isogray Treatment Planning System Using Measurement in Heterogeneous Phantom

*
*
AUTHORS

Asiyeh Golestani
^{
1
}
,
Mohammad Houshyari
^{
2
}
,
Ahmad Mostaar
^{
1
, *
}
,
Ali Jabbari Arfaie
^{
3
}

1 Department of Medical Physics and Biomedical Engineering, School of Medicine, Shahid Beheshti University of Medical Sciences, Tehran, IR Iran

2 Department of Radiotherapy, Shohada-e-Tajrish Hospital, School of Medicine, Shahid Beheshti University of Medical Sciences, Tehran, IR Iran

3 Shohada-e-Tajrish Hospital, Shahid Beheshti University of Medical Sciences, Tehran, IR Iran

**How to Cite:**
Golestani
A, Houshyari
M, Mostaar
A, Jabbari Arfaie
A. Evaluation of Dose Calculation Algorithms of Isogray Treatment Planning System Using Measurement in Heterogeneous Phantom,
Rep Radiother Oncol.
2015
; 2(3):e5320.
doi: 10.17795/rro-5320.

ARTICLE INFORMATION

**Reports of Radiotherapy and Oncology:**2 (3); e5320

**Published Online:**September 26, 2015

**Article Type:**Research Article

**Received:**July 11, 2015

**Accepted:**September 12, 2015

**DOI**: 10.17795/rro-5320

### Abstract

**Background:**
Radiotherapy is a high-energy ionizing radiation treatment for some kinds of cancers. Accuracy and the quality of radiotherapy treatment planning systems depend on the type of dose calculation algorithms that are utilized. There are uncertainties in the calculation of dose distribution, especially in a heterogeneous situation.

**Objectives:**
This study calculated and compared dose with model-based algorithms (superposition, collapsed cone convolution (CCC) and fast Fourier transform (FFT) and a measure-based algorithm (Clarkson in two modes: heterogeneity active and inactive).

**Materials and Methods:**
A heterogeneous phantom was used based on a semi-anthropomorphic phantom CIRS thorax 002 LFC. All of the tests were planned according to IAEA TEC-DOC 1583. All measurements were performed with a 6 MV photon beam of a linear accelerator (ELEKTA Compact) and Farmer ionization chamber. Five methods were utilized to calculate dose and were compared with measurement results as the gold standard.

**Results:**
In the homogeneity media all algorithms had good accuracy and dose difference was below 3%, but in the inhomogeneity situation dose difference increased and in some cases did not achieve agreement criteria. The Superposition algorithm overall has minimum deviations in all cases. However in some cases CCC had better accuracy. The Clarkson algorithm had maximum differences, especially when inhomogeneity correction was inactive.

**Conclusions:**
Dose calculation algorithms applied in radiotherapy treatment planning systems have different accuracy. Model-based algorithms have a better accuracy over measurement-based algorithms such as Clarkson. In the presence of large inhomogeneity, it is strongly recommended to activate manual inhomogeneity correction.

Keywords

Quality Assurance Dose Calculation Algorithms Treatment Planning Systems

### 1. Background

Radiotherapy is treatment with high-energy ionizing radiation for some kinds of cancers. In this procedure the prescribed dose of radiation was applied to damage cancer cells and stop them from growing and dividing. A high degree of accuracy regarding the delivery dose is one major requirement of radiotherapy. According to the international commission on radiation units and measurements, ICRP, the overall accuracy of the absorbed dose required to control a tumour and to avoid complications is 5% (1). Other groups proposed 3% - 4% for uncertainty in delivery dose (2-6). To achieve this level of accuracy, it is necessary to audit all steps of the dose calculation: beam data, dose calculation algorithm and so on. The IAEA has developed a practical test based on TRS430 to verify the accuracy of a dissymmetry system. In radiotherapy, monitor units (MU) or beam-on time are provided to deliver actual dose to a tumour. Error in calculation MU causes dose error and a curve in the treatment. Nowadays, a three-dimensional treatment planning system (TPS) has become common in radiotherapy and therefore verification methods for dose calculation have changed (7). Dose calculation algorithms can be broadly classified. One of the major differences between various dose calculation algorithms is the approximation of the mass stopping power (8, 9). Although Monte Carlo is the most accurate algorithm for dose calculation, it requires a long processing time (10-12). Therefore, other algorithms such as Clarkson and point kernel are used. There are uncertainties in the calculation of dose distribution especially in a heterogeneous situation, small field, or an interface region because of electronic disequilibrium. Over the last decade, pencil beam (PB) algorithms have become the standard for dose calculation algorithms in radiotherapy. Although it has acceptable accuracy and calculation speed, certain limitations of PB algorithms and higher accuracy can be addressed through some kernel-based models (8).

A convolution-superposition algorithm improved estimation of dose in an electronic equilibrium region. A pencil beam convolution (PBC) algorithm considers a narrow pencil beam and uses field intensity fluencies with the kernel that describes dose deposition around the primary photon pencil beam in water.

This study focuses on the dosimetric audit of the Isogray TPS (Dosisoft company) and investigation of the accuracy of its dose calculation algorithms.

### 2. Objectives

For this purpose an anthropomorphic phantom was used with a set of clinical test cases prepared by the IAEA that covers some typical clinical radiation techniques in three-dimensional radiotherapy.

### 3. Materials and Methods

All measurements were performed in a hospital with an accelerator (ELEKTA Compact) and photon having a low energy (6 MV) beam. Measurements were performed by an ionization chamber.

#### 3.1. Phantom

This study used a heterogeneous phantom based on a semi-anthropomorphic phantom CIRS thorax 002 LFC. This phantom is used for clinical test measurement of the chest wall and breast. The phantom was designed from tissue, bone and lung equivalent materials and it has holes for dose measurement by an ionization chamber. It consists of a 15 cm thick slice and 15 standard slices with 1 cm thickness. The 30 cm × 20 cm diameter body of Plexiglas (Perspex) has a mass density of 1.17 g/cm^{3} and relative electron density of 1.003. There is also a cylinder of Teflon as bony vertebra and special fiber as lung that is applied with a relative electron density 1.506 and 0.207, respectively. The phantom contains 10 holes to hold interchangeable rod inserts for an ionization chamber. The holes are numbered as shown in Figures 1 and 2. However the phantom has a limited number of available positions for dose measurements; it has slabs in between that accommodate a film dosimeter and measure dose in two dimensions. The phantom was scanned using a computed tomography (CT) scanner (Somatom Emotion 16-Slice, Siemens) that was calibrated using a special calibration agency. The CT images of the phantom in a DICOM format were transferred to the TPS. These images were used for planning the clinical tests.

Right, Ionization chamber measurements are performed in the 15 cm slice; left, standard slice (13).

Plugs number 1, 2, 3, 4 and 5 are tissue equivalent material, 6, 7, 8 and 9 are equivalent lung material and 10 is bone equivalent material (13).

Electron density relative to the CT number and maximum difference between CT numbers for same relative electron density was 20 Hounsfield units (Hus) (Figures 3 - 10).

#### 3.2. Clinical Test Cases

A set of clinical test cases as advised by IAEA TEC-DOC 1583 (13) were used to verify the range of basic treatment techniques applied in clinical practice. A detailed description of test cases and measurement points are given in Table 1.

^{a}

No. of Case | Description of Test | Reference Point | Measure Point(s) | Agreement Criteria, % |
---|---|---|---|---|

Single field, SSD = 100, angle Collimator and Gantry = 0°, size of field = 10 × 10 | 3 | |||

3 | 2 | |||

9 | 4 | |||

10 | 3 | |||

1 | 2 | |||

5 | 2 | |||

Oblique incidence, lack of scattering and tangential field, SSD = 100, angle Gantry = 90°, angle Collimator: depended to wedge | 1 | |||

1 | 3 | |||

Significant blocking of the field corner. SSD = 100, technique: SAD. angle Collimator = 45°,filed size 14 × 14 down to 8 × 10 | 3 | |||

3 | 3 | |||

Four field box, technique: SAD. Set-up at point: 5 | 5 | |||

5 (0°) | 2 | |||

5 (90°) | 3 | |||

5 (180°) | 3 | |||

5 ( 270°) | 3 | |||

Ʃ | ||||

6 (0°) | 4 | |||

6 (90°) | 3 | |||

6 (180°) | 4 | |||

6 (270°) | 3 | |||

Ʃ | ||||

10 (0°) | 3 | |||

10 ( 90°) | 4 | |||

10 (180°) | 3 | |||

10 (270°) | 4 | |||

Ʃ | ||||

MLC block of clinder: 10 cm, technique: SAD, set-up at point: 2, angle Gantry: 0°, angle Collimator: 0° | 2 | |||

2 | 3 | |||

7 | 4 | |||

L-shaped field with oblique incidence. Set-up at point:5, technique: SAD, angle Collimator:0°, angle Gantry: 45°, filed size 10 × 20 down to 6 × 12 | 3 | |||

3 | 3 | |||

7 | 5 | |||

10 | 5 | |||

Plan with asymmetric fields and wedge, technique: SAD. set-up at point: 3 | 5 | |||

5 (0°) | 2 | |||

5 (90°) | 4 | |||

5 (270°) | 4 | |||

Ʃ | ||||

Plan with non-coplanar field. technique: SAD, set-up at point: 5 | 5 | |||

5 (30°) | 3 | |||

5 (90°) | 3 | |||

5 (270°) | 3 | |||

Ʃ |

^{a}The case No. 5 was not performed because of lack of MLC in our radiotherapy center.

Eight cases and thirty points comprised the total number of measurements. The plan was for test cases to deliver a 2 Gy dose to the reference point in each case and then the calculated dose was obtained by different algorithms at other points. All of the measurements were performed by an ionization chamber. Due to lack of MLC, test case no. 5 was not performed.

#### 3.3. Dose Calculation Algorithms

Dose calculations were performed for some standard clinical tests with different algorithms to evaluate accuracy; the results were compared to measurements as a gold standard. The evaluated algorithms in this study were divided into two categories:

A) Measurement-based algorithm that includes Clarkson (two types: if manual heterogeneity correction is active named Clarkson* and if as default it is inactive and named Clarkson).

B) Model-based algorithms using point kernel convolution that includes fast-fourier transform (FFT), superposition and collapsed cone convolution (CCC).

These models account for density verification approximately for electrons and photons in lateral.

#### 3.4. Measurements

All measurements were carried out on an Elekta accelerator (Compact model, Elekta Oncology Systems) operating with a photon energy of 6 MV. The 6 MV photon beam was investigated in this study since higher energy photon beams have a less inhomogeneity correction effect (14).

For dose measurements in the phantom, a Farmer type ionization chamber (A19, Standard Imaging Co, USA) with electrometer (Supermax, Standard Imaging Co, USA) were employed. The chamber and electrometer were calibrated at secondary standard dosimetry laboratory (SSDL). For all measurement points, the absorbed dose was determined from ionization chamber measurements using the IAEA TRS 398 dosimetry code of practice (15). Pressure and temperature in each measurement condition were also measured.

#### 3.5. Analysis of the Results

IAEA TRS (430) determines criteria for dose measured (D_{meas}) and dose TPS calculated (D_{cal}) such that differences are normalized to the dose measured at the reference point (D_{meas, ref}) for each test case (16) (Equation 1).

### 4. Results

In the homogenous media, all of the verification doses determined by algorithms (Clarkson, Clarkson*, Superposition, CCC, FFT) approximately have a highly acceptable degree of agreement between calculations and measurements. For case 1, a single square field, all of the algorithms had acceptable accuracy and their errors were less than the agreement criteria (± 2%). For point 9 in lung equivalent material, the calculated dose was overestimated, however, it was less than the agreement criteria. For point 10 in bone equivalent material, the calculated dose was underestimated. For all algorithms the amount of MU/D was equal to 1.15.

The Clarkson* and Clarkson algorithms had more errors (than) the dose in the tangential beam (case 2). All of algorithms were in compliance with agreement criteria (± 3%). The accuracy of the model-based algorithms (CCC, Superposition and FFT) were equal, and Clarkson* and Clarkson were also equal.

#### 4.1. In Case 3

In this test, all the algorithms overestimated the dose in the field with a wedge, although all of the results in the range were acceptable. The amount of (MU/D) = 1.185 for measurement-base algorithms and for model- based algorithms was equal to 1.181.

In the four fields box test (case 4), the reference point was in the Plexiglas, point 10 in the bone equivalent and the central axis in the vertical beams. Point 6 was in the lunge equivalent material and the central axis in the lateral beams. According to the measurements the Collapsed con, superposition and FFT then Clarkson were accurate, and in some tests, Clarkson failed. The mount of (MU/D) = 1.23 for Clarkson and for point kernels that were model-based was 1.25, the amount MU for the model-base was equal, but the dose calculations in points 6, 10 were different. Accuracy for Collapsed con and Superposition and FFT were close and acceptable, but for Clarkson it was different and for some measurements were higher than the criteria agreement. For the L-shape field (case 6), all algorithms passed, in other words, the amount of error was within the criteria agreement, except for point 7 (in the lunge equivalent material) for algorithm FFT, it was above criteria agreement and failed. Three points were tested, one of them was in plastic water (point 3), where all the algorithms passed and another point was in a lung equivalent media (point 7), where FFT and Clarkson* failed (agreement criteria was 4%). The last point (point 10) was in bone equivalent media, and all of the algorithms passed. In the asymmetrical wedge field (case 7) for model-based algorithms (Superposition, CCC and FFT) the amounts of calculated and measured dose were approximately equal and the amount of error was very small. But for Clarkson, the measured dose was lower than the calculated dose and indicates that Clarkson overestimates the dose, however, Clarkson* underestimates the dose in this case.

In case 8: Deviation outside the agreement criteria was found for Clarkson, and the dose calculation overestimates. For model-based algorithms, in spite of this, the range of agreement was underestimated. For the lateral field, Clarkson and Clarkson* had error that was beyond the agreement criteria.

The result for algorithms being tested for energy of 6 MV is shown in the histogram, and the percentage of results outside of the agreement criteria for all algorithms are shown.

### 5. Discussion

#### 5.1. Phantom

This study utilized a heterogeneous phantom designed to represent lung, bone and tissue, utilizing equivalent material based on IAEA methodology and previous studies of this phantom. The phantom can be easily shared and transported between hospitals. It is easy to setup and use for testing. More details of the suitability of different phantoms for these types of tests are given (17, 18).

#### 5.2. TPS Dose Calculation Algorithms

The accuracy and quality of TPS depends on the type of algorithms utilized in every step of TPS and has been the topic of numerous studies (19). The results for a homogeneous media comply between measurement and calculated dose. The maximum difference results were: Superposition algorithm (-1.7), CCC (1.8), FFT (-2.1) and for Clarkson and Clarkson* was equal to 3.14, for lateral beam the maximum difference for Superposition was 2.9 and for CCC, FFT, Clarkson and Clarkson* was 3.1, 2.9, 4.7, 6.9 in case test 8, respectively. Except for Clarkson and Clarkson* in the lateral field, the other algorithms passed. Aarup et al. reported good agreement for five algorithms that were tested (20). Clarkson is a simple algorithm that assumes a photon scatter interaction medium. However, it does not include any change in scatter radiation and electron transport in a heterogeneous medium and it cannot consider a linear process due to insufficient accuracy in this algorithm (21). It is acceptable for clinical applications in the cranial and pelvis lesions (22). For our test, although Clarkson had the most diversions, it can be used for clinical tests. A maximum difference calculation dose was seen (3.17%) for equivalent bone material and (4.11%) for equivalent lung material, which underestimates the doses calculated. In the study by Ambolt et al. the Clarkson algorithm has an acceptable degree of agreement between calculations and measurement (23). Linthout et al. also showed good agreement between dose distribution and measured dose by Clarkson algorithm, but for a tissue-air interface the dose calculated by Clarkson overestimates by 12% (24). In test 2 that uses a tangential beam, it calculated an overestimated dose compared to dose measurement. There was a difference (-2.5%) between measurement and calculated dose. The CCC algorithm predicts a dose near to inhomogeneity exactly (25-28). The CCC algorithm as used for the Monte Carlo calculation method and the results show that it has an accuracy of 1% in complex geometries (29). The CCC algorithm had reasonable speed and accuracy in dose calculation (30). One of the limitations in the accuracy CCC is found in a lateral multiple scatter dose (8). The results show that the CCC algorithm predicts in accordance between dose calculation and measurement. Maximum difference was seen in a lung equivalent material (-3.9%) and for a bone equivalent material it was -2.9%, where both results overestimate. Asnaashari showed model-based algorithms have better compliance than measurement-based algorithms (31). Linthout et al. showed the three different calculation algorithms (Clarkson, PB and the Collapsed Cone (CC) algorithms) overestimate the dose in the PTV at the boundary with the low-density tissue, with 12, 10 and 7% for the Clarkson, the PB and the CC algorithm, respectively (24). In another study, Das et al. showed that the Collapsed cone algorithm has better accuracy than the PB algorithm in the interface of inhomogeneous media (27). Nisbet et al. showed Collapsed cone has better accuracy than the PB model in interface lung/tissue and Collapsed cone and the PB model also have inadequate accuracy in interface bone/tissue (25). The maximum deviation was 3.2% for the Sc algorithm in a lung region. Zhang et al. reported that superposition has better accuracy between FFTC and Superposition in the field edge and heterogeneous media (32). In our study FFTC had 3.8% maximum deviation for a lung equivalent media. Vicente et al. showed a 9.7% difference between measured and calculated in cork media with a 6 MV beam by the Superposition algorithm (33). Rutonjski et al. reported a mean 6% deviation with the Superposition algorithm (34). Scholz et al. showed that the Superposition algorithm had a mean 3% difference in various fields (35).

This study applied the methodology described in IAEA TECDOC1583 to perform a TPS audit and evaluate dose calculation algorithms. Maximum deviation was found for the Clarkson algorithm with active heterogeneity correction as a field of interest in (ROI) contouring, and error is reduce. Model-based algorithms had better accuracy, and collapsed cone and Superposition were better. These algorithms help to calculate dose with a higher accuracy.

### Footnotes

### References

*1.*

Shalek RJ. *Determination of Absorbed Dose in a Patient Irradiated by Beams of X or Gamma Rays in Radiotherapy Procedures: ICRU Report 24.* 1976; **24**

*2.*

Brahme A. *Accuracy Requirements and Quality Assurance of External Beam Therapy with Photons and Electrons.* 1988;

*3.*

Brenner DJ. Dose, volume, and tumor-control predictions in radiotherapy.
*Int J Radiat Oncol Biol Phys.* 1993; **26**(1) : 171
-9
[PubMed]

*4.*

Goitein M. Nonstandard deviations.
*Med Phys.* 1983; **10**(5) : 709
-11
[DOI][PubMed]

*5.*

Mijnheer BJ, Battermann JJ, Wambersie A. What degree of accuracy is required and can be achieved in photon and neutron therapy?
*Radiother Oncol.* 1987; **8**(3) : 237
-52
[PubMed]

*6.*

Brenner DJ. Reply to: Precision and accuracy in radiotherapy.
*Radiotherapy and Oncology.* 1989; **14**(2) : 163
-7
[DOI]

*7.*

Sellakumar P, Arun C, Sanjay SS, Ramesh SB. Comparison of monitor units calculated by radiotherapy treatment planning system and an independent monitor unit verification software.
*Phys Med.* 2011; **27**(1) : 21
-9
[DOI][PubMed]

*8.*

Ahnesjo A. Collapsed cone convolution of radiant energy for photon dose calculation in heterogeneous media.
*Med Phys.* 1989; **16**(4) : 577
-92
[DOI][PubMed]

*9.*

Francescon P, Cavedon C, Reccanello S, Cora S. Photon dose calculation of a three-dimensional treatment planning system compared to the Monte Carlo code BEAM.
*Med Phys.* 2000; **27**(7) : 1579
-87
[DOI][PubMed]

*10.*

Fippel M. Fast Monte Carlo dose calculation for photon beams based on the VMC electron algorithm.
*Med Phys.* 1999; **26**(8) : 1466
-75
[DOI][PubMed]

*11.*

Ahnesjo A, Aspradakis MM. Dose calculations for external photon beams in radiotherapy.
*Phys Med Biol.* 1999; **44**(11) : 99
-155
[PubMed]

*12.*

Fippel M, Haryanto F, Dohm O, Nusslin F, Kriesen S. A virtual photon energy fluence model for Monte Carlo dose calculation.
*Med Phys.* 2003; **30**(3) : 301
-11
[DOI][PubMed]

*13.*

Vatnitsky S. *Commissioning of Radiotherapy Treatment Planning Systems: Testing for Typical External Beam Treatment Techniques.* 2008;

*14.*

Papanikolaou N, Battista JJ, Boyer AL. *Tissue inhomogeneity corrections for megavoltage photon beams.* 2004; **65**

*15.*

IAEA (International Atomic Energy Agency), . *Absorbed dose determination in external beam radiotherapy: an international code of practice for dosimetry based on standards of absorbed dose to water: Technical Report Series 398* 2000;

*16.*

*Commissioning and quality assurance of computerized planning systems for radiation treatment of cancer: Technical Report Series 430* 2005;

*17.*

Knoos T,
Wieslander E,
Cozzi L,
Brink C,
Fogliata A,
Albers D,
et al.
Comparison of dose calculation algorithms for treatment planning in external photon beam therapy for clinical situations.
*Phys Med Biol.* 2006; **51**(22) : 5785
-807
[DOI][PubMed]

*18.*

Gershkevitsh E,
Schmidt R,
Velez G,
Miller D,
Korf E,
Yip F,
et al.
Dosimetric verification of radiotherapy treatment planning systems: results of IAEA pilot study.
*Radiother Oncol.* 2008; **89**(3) : 338
-46
[DOI][PubMed]

*19.*

Van Dyk J. Quality assurance of radiation therapy planning systems: current status and remaining challenges.
*Int J Radiat Oncol Biol Phys.* 2008; **71**(1 Suppl) : S23
-7
[DOI][PubMed]

*20.*

Aarup LR,
Nahum AE,
Zacharatou C,
Juhler-Nottrup T,
Knoos T,
Nystrom H,
et al.
The effect of different lung densities on the accuracy of various radiotherapy dose calculation methods: implications for tumour coverage.
*Radiother Oncol.* 2009; **91**(3) : 405
-14
[DOI][PubMed]

*21.*

Khan FM. *The physics of radiation therapy.* 1984; : 178
-80

*22.*

Miften MM, Beavis AW, Marks LB. Influence of dose calculation model on treatment plan evaluation in conformal radiotherapy: a three-case study.
*Med Dosim.* 2002; **27**(1) : 51
-7
[PubMed]

*23.*

Ambolt P,
Bowden J,
Cooper M,
Henry A,
McCormack S,
Morgan A,
et al.
Commissioning a Superposition Dose Calculation Model for 6 MV Photon Beams on a Radiotherapy Treatment Planning System.
*J Clin Oncol.* 2007; **19**(3)[DOI]

*24.*

Linthout N, Verellen D, Van Acker S, Voordeckers M, Bretz A, Storme G. Evaluation of dose calculation algorithms for dynamic arc treatments of head and neck tumors.
*Radiother Oncol.* 2002; **64**(1) : 85
-95
[PubMed]

*25.*

Nisbet A, Beange I, Vollmar HS, Irvine C, Morgan A, Thwaites DI. Dosimetric verification of a commercial collapsed cone algorithm in simulated clinical situations.
*Radiother Oncol.* 2004; **73**(1) : 79
-88
[DOI][PubMed]

*26.*

Weber L, Nilsson P. Verification of dose calculations with a clinical treatment planning system based on a point kernel dose engine.
*J Appl Clin Med Phys.* 2002; **3**(2) : 73
-87
[DOI][PubMed]

*27.*

Das IJ, Kase KR, Meigooni AS, Khan FM, Werner BL. Validity of transition-zone dosimetry at high atomic number interfaces in megavoltage photon beams.
*Med Phys.* 1990; **17**(1) : 10
-6
[DOI][PubMed]

*28.*

Irvine C, Morgan A, Crellin A, Nisbet A, Beange I. The clinical implications of the collapsed cone planning algorithm.
*Clin Oncol.* 2004; **16**(2) : 148
-54
[PubMed]

*29.*

Wang L, Yorke E, Chui CS. Monte Carlo evaluation of tissue inhomogeneity effects in the treatment of the head and neck.
*Int J Radiat Oncol Biol Phys.* 2001; **50**(5) : 1339
-49
[PubMed]

*30.*

Battista JJ, Sharpe MB. True three-dimensional dose computations for megavoltage x-ray therapy: a role for the superposition principle.
*Australas Phys Eng Sci Med.* 1992; **15**(4) : 159
-78
[PubMed]

*31.*

Asnaashari K, Nodehi MR, Mahdavi SR, Gholami S, Khosravi HR. Dosimetric comparison of different inhomogeneity correction algorithms for external photon beam dose calculations.
*J Med Phys.* 2013; **38**(2) : 74
-81
[DOI][PubMed]

*32.*

Zhang P,
Simon A,
De Crevoisier R,
Haigron P,
Nassef MH,
Li B,
et al.
A new pencil beam model for photon dose calculations in heterogeneous media.
*Phys Med.* 2014; **30**(7) : 765
-73
[DOI][PubMed]

*33.*

Garcia-Vicente F, Minambres A, Jerez I, Modolell I, Perez L, Torres J. Experimental validation tests of fast Fourier transform convolution and multigrid superposition algorithms for dose calculation in low-density media.
*Radiother Oncol.* 2003; **67**(2) : 239
-49
[DOI]

*34.*

Rutonjski L,
Petrovic B,
Baucal M,
Teodorovic M,
Cudic O,
Gershkevitsh E,
et al.
Dosimetric verification of radiotherapy treatment planning systems in Serbia: national audit.
*Radiat Oncol.* 2012; **7** : 155
-63
[DOI][PubMed]

*35.*

Scholz C, Schulze C, Oelfke U, Bortfeld T. Development and clinical application of a fast superposition algorithm in radiation therapy.
*Radiother Oncol.* 2003; **69**(1) : 79
-90
[PubMed]

## LEAVE A COMMENT HERE: