|Year : 2022 | Volume
| Issue : 1 | Page : 50-56
Implementation and validation of anisotropic analytical algorithm in eclipse treatment planning system for indigenous telecobalt machine (Bhabhatron II)
Raghavendra Hajare, Md Basit Ali, KK Sreelakshmi, M Anil Kumar, Rituraj Kalita, Umesh M Mahantshetty
Department of Radiation Oncology, Homi Bhabha Cancer Hospital and Research Centre, Visakhapatnam, Andhra Pradesh, India
|Date of Submission||13-Jul-2021|
|Date of Decision||22-Dec-2021|
|Date of Acceptance||23-Dec-2021|
|Date of Web Publication||31-Mar-2022|
Ms. K K Sreelakshmi
Department of Radiation Oncology, Homi Bhabha Cancer Hospital and Research Centre, Aganampudi, Gajuwaka Mandalam, Visakhapatnam - 530 053, Andhra Pradesh
Source of Support: None, Conflict of Interest: None
| Abstract|| |
Background: The photon dose calculation model anisotropic analytical algorithm (AAA) available with eclipse integrated treatment planning system (TPS) (Varian) supports telecobalt dose calculation from Version 13.6 onward. Formerly, pencil beam convolution (PBC) was used for modeling telecobalt machines. Eclipse TPS no longer supports PBC dose calculation algorithm in v13.6 and above. The AAA dose calculation model is a three-dimensional PBC/superposition algorithm. Its configuration is based on Monte-Carlo-determined basic physical parameters that are tailored to measured clinical beam data. Aim: The study investigated the feasibility of clinical implementation of AAA in Eclipse TPS for Bhabhatron II. Materials and Methods: The indigenous telecobalt machine, Bhabhatron II, was configured as a generic machine because an inbuilt machine model for the same was not available in Varian Eclipse TPS algorithm library. In such a scenario, it was necessary to evaluate and validate dosimetric parameters of the TPS because improper tailoring would cause errors in dose calculations. Beam data measurements of the machine were carried out which were used for configuration of the algorithm. Result: After successful configuration, a variety of plans created in TPS were executed on the machine and subsequently evaluated. Conclusion: From this study, we concluded that AAA-based dose calculation in TPS is very well suited for accurate dose calculations for telecobalt machine and can be implemented for clinical use.
Keywords: Anisotropic analytical algorithm, Bhabhatron II, percentage depth dose, telecobalt, treatment planning system
|How to cite this article:|
Hajare R, Ali MB, Sreelakshmi K K, Kumar M A, Kalita R, Mahantshetty UM. Implementation and validation of anisotropic analytical algorithm in eclipse treatment planning system for indigenous telecobalt machine (Bhabhatron II). J Med Phys 2022;47:50-6
|How to cite this URL:|
Hajare R, Ali MB, Sreelakshmi K K, Kumar M A, Kalita R, Mahantshetty UM. Implementation and validation of anisotropic analytical algorithm in eclipse treatment planning system for indigenous telecobalt machine (Bhabhatron II). J Med Phys [serial online] 2022 [cited 2022 Dec 5];47:50-6. Available from: https://www.jmp.org.in/text.asp?2022/47/1/50/341437
| Introduction|| |
Telecobalt units are still widely used in many developing countries for cancer treatment and are preferred over medical linear accelerators (linacs) because of its modest cost, reduced maintenance charges, lower power requirements, ease of operation, and usually lesser down time. In most of these centers, patient treatment is done on the basis of manually calculated treatment time, depicting the dose to a point (usually the tumor center), primarily on the central axis. Visualizing the dose distribution on the patient computed tomography (CT) in a treatment planning system (TPS) can give much more information to help us make the correct choices regarding different aspects of a plan. It also helps us to predict the toxicities which the patient might encounter during or after the course of treatment.
The Bhabhatron II Tungsten, Asymmetric, motorized Wedge (TAW) is an indigenously produced, IEC (60601-2-11) compliant, affordable telecobalt machine which possesses advanced features such as asymmetric jaw, motorized wedge, and programmable control console as compared to its other conventional counterparts.,,,,, A separate vendor-specific TPS was not purchased with Bhabhatron by the institute.
Varian Eclipse anisotropic analytical dose calculation algorithm (AAA) (Version 15.6.06) is a Type B algorithm which supports the use of cobalt treatment units in external beam planning and IRREG two-dimensional (2D) planning. It also allows the configuration of blocks and standard wedges for cobalt treatment units.
The AAA model for dose calculation comprises two key components, the configuration algorithm and the dose calculation algorithm.
The configuration algorithm determines the fundamental physical parameters such as photon energy spectrum, mean radial energy, and scatter kernels, which characterizes the photon and electron fluence and their energy spectra in the treatment beam. Determining all these required parameters through measurements is not practically possible. Hence, AAA uses Monte Carlo precalculated parameters which are then modified according to measured data. These parameters are then stored after completion of configuration process, for retrieval during the dose calculation process.
The dose calculation algorithm uses a separate convolution model for primary photons, scattered extra-focal photons, and electrons scattered from the beam-limiting devices. To apply convolutions, the clinical broad beam is divided into small, finite-sized beamlets. The cross-section of the beamlet is the calculation voxel resolution. The final dose distribution is obtained by the superposition of the dose calculated with photon and electron convolutions for the individual beamlets.
All model parameters for AAA are computed in a water-equivalent medium. While performing dose distribution calculation, to account for heterogeneity in media, density scaling according to patient tissue is done.
This study aimed to evaluate accuracy and clinical implementation of AAA in Eclipse TPS for the indigenous telecobalt machine.
| Materials and Methods|| |
The modeled telecobalt machine was a Bhabhatron II TAW (Panacea Medical technologies Pvt. Ltd.). It housed a cobalt-60 source of 2.3 cm diameter and 3.7 cm length with activity 9147.8 Ci as loaded on August 02, 2019, which had an output of 290.47 cGy/min on August 22, 2019. The output is the central axis dose, at the depth of dose maxima in a 10 cm × 10 cm square field at 80 cm source-to-axis distance (SAD). The specified maximum capacity of the machine was 250 roentgen per minute at 1 m. The Bhabhatron II had asymmetric Y jaws, symmetric X jaws, and a motorized wedge of 60°. The available field sizes were from 3 cm × 3 cm to 35 cm × 35 cm at the normal treatment distance of 80 cm. A set of manual/physical wedge filters of various wedge angles 15°, 30°, 45°, and 60° were also provided.
Beam data acquisition
Different beam data measurements, as detailed in [Table 1], were performed as per requirement of Eclipse AAA algorithm.
|Table 1: Details of beam data measurements performed for configuring anisotropic analytical algorithm|
Click here to view
The measured depth dose curves, cross line profiles, diagonal profiles, and output factors (Scp) for both open and wedged fields were imported into Eclipse to calculate beam data for the anisotropic analytical algorithm (AAA Version 15.6.06 Varian Medical Systems). The grid size used for the calculation of configured data was 2.5 mm. In Eclipse, the configuration of telecobalt was done as a generic machine. This enables the algorithm to draw generic parameters from a library built from Monte Carlo simulations of the treatment head for configuring the machine. These generic parameters are then used by configuration algorithm to build a phase space which is consistent with the measured data and that will be used for clinical calculations.
For cobalt source, effective spot size is larger than linear accelerator. The spot size effects the penumbral width. These parameters were adjusted and fine-tuned so as to match calculated and measured beam profiles.
The timer setting for cobalt fields was determined based on the reference dose and time required to deliver that reference dose at the calibration date, which was thereafter corrected for decay.
A CT scan of the carbon fiber couch provided with Bhabhatron was taken and its Hounsfield Unit (HU) value at different parts was evaluated. The couch transmission factor was also measured at different parts of the treatment couch. After analyzing these data, a couch was modeled in the TPS as a support structure with a HU value of –780.
The modeled couch was inserted as a support structure in patient-specific plans and the resulting TPS dose was compared with manual calculation and measurement.
Beam data verification
To compare the measured data with configured data, in beam analysis window, beam data were calculated with 2.5 mm calculation grid size.
A variety of dose comparison tests, as detailed below, were conducted to verify the treatment planning accuracy.
Point dose verification
A virtual phantom of 40 cm × 40 cm × 40 cm (1 g/cc density) was created in the TPS for validation. Different plans were created on virtual phantom to compare treatment time provided by TPS with manually calculated treatment time from measured data. This was done at multiple depths for a normalized dose of 100 cGy at the reference depths.
A CT of combination of slab phantom (HE Solid water; Gammex, USA) (30 cm × 30 cm × 20 cm) was taken with a 0.125cc (SNC 125c™; Sun Nuclear Corp, USA) chamber placed at 10 cm depth. Different plans of multiple field sizes, symmetric, asymmetric, wedged, open, and extended SSD (90 cm and 100 cm) were created in TPS on this CT. This ensured that the basic parameters of dose computation process were evaluated for different beam configurations. All plans were calculated using AAA algorithm. All these plans were then irradiated on identical setup on the machine. The measured doses were then compared with the doses depicted in TPS.
A custom-made phantom [Figure 1] was used to perform measurements in inhomogeneous media. The nonwater equivalent inhomogeneous part was made using thermocol (mean HU value –930 HU), wax, and white cement (mean HU value of –80 HU and 1400 HU, respectively). Its dimension was about 29 cm × 24 cm × 7.8 cm. The chamber slab of SNC 125c was kept just below this unit. A backscatter of 9 cm was placed below the chamber slab. Above the thermocol unit, a build-up slab of 1 cm was kept. The physical depth of measurement was 9.8 cm and the corresponding water equivalent depth was 7.3 cm at central axis. This whole assembly was scanned in CT and different plans of multiple field sizes for both SSD (SSD = 80 cm) and SAD setup were made. Plans were irradiated on machine and the measured dose was compared to TPS predicted dose.
|Figure 1: Custom-made inhomogeneous phantom. (a) Top view of inhomogeneous part in phantom. (b) Transverse view of the complete inhomogeneous phantom setup used for measurement. (c) Axial computed tomography image slice|
Click here to view
Wedge transmission and output factors
For output factor's comparison, we created plans of different field sizes on the above-mentioned slab phantom normalized to reference point at a depth of 10 cm. The factor was then calculated by dividing the treatment time of a given field with reference field size time. This was then compared with the measured values.
To validate wedge transmission factors, plans were made on the slab phantom in TPS with and without wedges. The created TPS plans were then irradiated on the treatment machine. The measurements were compared with the factors obtained from the TPS.
Patient dose verification
Some patient-specific QA plans were made in TPS on the aforementioned slab phantom. These plans were then irradiated on the machine in the planned setup. The dose measured was then compared with the dose predicted by TPS.
Verification of fluence was performed using the SNC 2D diode array MapCHECK®3. It consisted of 1527 diode detectors (solid state) with 7.07 mm spacing in an array size of 26 cm × 32 cm. Each detector had an active area of 0.48 mm × 0.48 mm. This array was used with an accessory (MapPHAN) which provided a build-up of 5cm to MapCHECK®3. It is modeled in TPS using a synthetic CT provided by the manufacturer.
Initially, array and dose calibration of MapCHECK®3 was performed, for cobalt energy. Array calibration determined relative sensitivity differences between detectors and stored it as individual correction factor for each detector, whereas absolute dose calibration correlated the counts to the known data establishing the absolute dose calibration factor.
Different open field sizes, various wedged fields, some extended SSD (90 cm and 100 cm) plans, and some half beam block fields were irradiated on MapCHECK®3 assembly. Some fields with rotated collimator (collimator angle ≠ 0) were also irradiated. The measured dose distribution by the MapCHECK®3 device was compared with the dose distribution calculated by the TPS. The gamma analysis was done using SNC patient software (Sun Nuclear Corporation, USA). Dose distributions were analyzed using gamma criteria of 3% dose difference and 3 mm distance to agreement (DTA) as well as 2% dose difference and 2 mm DTA.
| Results|| |
Percentage depth dose and profiles
[Figure 2] shows the gamma error histogram for open fields, generated during the optimization process of AAA, for percentage depth dose (PDD) data after and before dmax and for profiles in three regions: inside field/flattened region, penumbra region, and outside field/umbra region. This compares the calculated data with the processed measured data. The γ index was computed with default setting of DTA = 3 mm and dose difference (d) = 1%. The figure shows a global agreement between measured and the optimized data.
|Figure 2: Gamma error histogram for open fields (a) PDD after Dmax; (b) PDD before Dmax; (c) Profiles inside field; (d) Profiles outside field; (e) Profiles in penumbra|
Click here to view
There was also good agreement between TPS calculated and RFA measured PDDs. The mean of difference between TPS calculated and measured PDD at 10cm for different field sizes and their corresponding standard deviation were 0.92 ± 0.3 (Open) [Supplementary Table 1]c[Additional file 1], 0.45 ± 0.2 (W15), 0.4 ± 0.1 (W30), 0.24 ± 0.2 (W45), 0.63 ± 0.1 (W60).
PDD values of different fields from TPS plans were compared to BJR data (Supplement 17) and were found to have a max difference of 0.8% [Supplementary Table 1]a and [Supplementary Table 1]b[Additional file 2].
[Figure 3] shows the configured energy spectrum. AAA configures cobalt machines with a continuous energy spectrum without its characteristic peaks at 1.17 and 1.33 MeV.
The electron contamination curve representing laterally integrated electron contamination doses at different depths was evaluated and was found to be acceptable. As shown in [Figure 4], it peaked at surface with a fall-off tail at increasing depth. The smoothing factors for (Sigma values for Gaussian) electron contamination were Sigma 0 = 75.8 and Sigma 1 = 108.29.
The Scp for field sizes of 5 cm × 5 cm to 35 cm × 35 cm varied from 0.903 to 1.126 [Supplementary Table 2]a. The Scp for wedged field sizes from 5 cm × 5 cm to 15W × 20 cm varied from 0.895–1.091 (W15) [Supplementary Table 2]c, 0.892–1.1 (W30), 0.882–1.116 (W45), and 0.893–1.106 (W60) [Supplementary Table 2]e[Additional file 3].
The mean of the percentage difference between the RFA measurement and TPS calculated Scp for square and rectangular fields in open and wedged fields were as shown in [Table 2].
|Table 2: The mean±standard deviation of percentage difference between RFA measured and calculated output factors for different field sizes of open and wedged fields|
Click here to view
The difference in couch transmission factor between measurement and TPS model was 0.0098 (1.04%).
Point dose verification
The maximum difference between treatment time predicted by TPS and the manual calculated treatment time (from measured data) for a dose of 100 cGy at different prescribed depths was 1.2 s [Supplementary Table 3][Additional file 4].
The mean ± standard deviation of the percentage difference between TPS calculated and measured point doses were 2.83 ± 0.5 [Supplementary Table 5][Additional file 5] and 3.88 ± 0.5 for square and rectangular field sizes, respectively, in SSD (SSD = 80 cm) setup, for 10 cm depth. For square field sizes in isocentric setup, it was 3.70 ± 0.7 and for point doses at 10 cm depth in different wedged field sizes, it was 2.93 ± 1.
Measurements in inhomogeneous phantom gave a mean ± standard deviation of the percentage difference between TPS calculated and measured point doses of 2.54 ± 0.5 [Supplementary Table 4][Additional file 6] and 2.74 ± 0.8 for SSD (SSD = 80 cm) and SAD setup, respectively.
Wedge transmission and output factors
The mean of percentage differences between TPS calculated and measured wedge transmission factors for different field sizes were 0.46 ± 0.6 [Supplementary Table 5], 0.11 ± 0.7, 0.64 ± 0.6 and 0.49 ± 1.1 for W15°, W30°, W45°, and W60°, respectively. For open field, output factor at 10 cm depth, it was 0.46 ± 0.57 for SSD (SSD = 80cm) setup and 0.51 ± 0.75 for SAD setup.
Patient dose verification
A maximum difference of 1.2 s was found between manually calculated and TPS calculated treatment times for a dose of 100 cGy irrespective of setup. The statistical mode of the differences was 0.6 s.
[Table 3] shows the results from MapCHECK®3 irradiation measurements. The table includes the mean of gamma pass rate for measurement of various field sizes, its standard deviation, maximum and minimum gamma pass rates. The mean gamma percentage was calculated as the averages of per-field gamma pass percentage. The tolerance for gamma evaluation was kept as dose difference of 3% and a DTA of 3 mm.
|Table 3: Mean gamma pass rate, standard deviation, maximum gamma pass rate, and minimum gamma pass rate of different field sizes for open and wedged fields|
Click here to view
| Discussion|| |
The machine was configured as a generic machine because an inbuilt machine model for Bhabhatron II was not available in Varian Eclipse TPS algorithm library. The Monte-Carlo precalculated fundamental parameters are not now specific to any machine model. Thus, an inherent possibility of discrepancy exists between the measured data and calculated output values. These differences may be either due to intrinsic approximations or limitations of the model or due to insufficient optimization/tailoring of the fundamental physical parameters in the configuration phase.
From data summarized in this paper, there is a good agreement between measured and computed data. The deviations observed are within acceptable range. During configuration, the algorithm calculates the gamma index with an inbuilt setting (1%, 3 mm) and the mean gamma was found to be within 1. The PDD value at 10 cm depth and Scp for open and wedged fields is within 1%. The point dose differences for open and wedged field widths are within 4%. A systematic under estimation of absolute point dose prediction is observed in the TPS. Although the AAPM Task group 53 in its report recommends a tolerance criterion of 1.5% for central axis doses in rectangular fields, a higher difference was accepted keeping in view the simple planning conditions prevalent for cobalt machine in the institute. It was observed that if no heterogeneity correction was applied, the dose difference between measured and TPS calculated value systematically reduces, thus implying that there is an under estimation of scatter contribution (about 2% on an average) by the algorithm for any kind of heterogeneous media. The mean gamma passing rate with 3% dose difference and 3 mm DTA is more than 95% for the most of the cases.
The minimum modeled (and clinically used) field size in TPS was restricted to 5 cm × 5 cm, as smaller field sizes, although physically possible in machine, introduce more measurement errors and modeling them using generic machine parameter would have increased the deviations to unacceptable values. The maximum modeled field size was 35 cm × 35 cm, which was the largest possible physical field size.
For each wedge angle, other than 15°, the machine has two wedges: one for smaller field sizes (5 cm × 5 cm–16 cm × 10 cm) and another for larger field sizes (used for field sizes between 16 cm × 10 cm and 20 cm × 15 cm). For such wedge angles, only the wedge pertaining to the largest field size was modeled in TPS as for a given wedge angle, multiple wedges cannot be modeled in Eclipse planning system in a single beam model.
The beam data for cobalt can also be used for IRREG planning workspace in which planning can be done on 2D images or without using any images. The algorithm calculates treatment time according to the dose per field, field size, and depth information provided. The manually calculated treatment time matched with the TPS provided in the limited trials conducted by us.
The implementation of AAA-based dose calculation in TPS for Bhabhatron was done in our institute. It has helped in visualizing and ensuring sufficient dose coverage for many palliative spine treatments. Decisions about appropriate wedge angle to be used in anterolateral head-and-neck treatments have been made easy. A case of glioblastoma with vertex (noncoplanar) beam was treated, in which the couch rotation and gantry angle were decided based on the optimum dose distribution in the TPS. Chest wall treatments are done using bitangential fields with gantry angle, field sizes, and depths decided based on dose coverage seen in TPS.
| Conclusion|| |
The Bhabhatron II TAW was successfully modeled in the Varian Eclipse TPS using AAA dose calculation algorithm. Multiple complementary as well as redundant verifications of beam data were carried out and were found satisfactory.
Our study suggests that the Varian AAA for Bhabhatron II TAW is suitable for clinical use. This will help to visualize the dose distribution and plan the treatment more efficiently which is the need of the hour of modern era. The implementation of AAA-based dose calculation in TPS is very well suited for accurate radiation therapy treatment planning and its clinical usage will minimize the uncertainties in delivery.
The resources to do this work were provided by the Dept of radiation Oncology, Homi Bhabha Cancer Hospital and Research Centre, Visakhapatnam. We thank Mr. Subhabrat Dash and Mr. Pritam Parab for assistance with data collection. We would like to thank the software support team of Varian for their assistance. The engineering team of panacea medical technologies was also helpful.
Financial support and sponsorship
Conflicts of interest
There are no conflicts of interest.
| References|| |
Ravichandran R. Has the time come for doing away with Cobalt-60 teletherapy for cancer treatments. J Med Phys 2009;34:63-5.
] [Full text]
Jayarajan K, Kar DC, Singh M. Role of Bhabhatron in Cancer Care in India. In: Proceedings of International Conference on Peaceful Uses of Atomic Energy-2009. Vol. 2; 2009.
Kumar HS, Sharma A. A Physicist Description of Indigenous Telecobalt Machine Bhabhatron-II TAW. Global Journal of Science Frontier Research: A Physics and Space Science 2015;15.
Kar DC, Jayarajan K, Sharma SD, Singh M, Subrahmanyam GV. The bhabhatron: An affordable solution for radiation therapy. Biomed Imaging Intervention J 2008;4:E50-12.
Sahani G, Kumar M, Sharma PD, Sharma DN, Chhokra K, Mishra B, et al.
Compliance of Bhabhatron-II telecobalt unit with IEC standard –
Radiation safety. J Appl Clin Med Phys 2009;10:120-30.
Kumar R, Sharma SD, Phurailapam R, Despande DD, Kannan S. Performance characteristics of indigenously developed Bhabhatron-telecobalt unit. J Med Phys 2005;30:41.
Sharma SB, Sarma G, Barthakur M, Goswami P, Yadav B, Goswami S. Installation, commissioning, and performance evaluation of Bhabhatron-II TAW – An Indian-made telecobalt unit. J Radiat Med Trop 2020;1:38. [Full text]
Varian Medical Systems. Eclipse Photon and Electron Algorithms Reference Guide. Palo Alto (CA): Varian Medical Systems; 2018.
Fogliata A, Nicolini G, Vanetti E, Clivio A, Cozzi L. Dosimetric validation of the anisotropic analytical algorithm for photon dose calculation: Fundamental characterization in water. Phys Med Biol 2006;51:1421.
Sievinen J, Ulmer W, Kaissl W. AAA photon dose calculation model in Eclipse. Vol. 118. Palo Alto (CA): Varian Medical Systems; 2005. p. 2894.
Miller D. American Association of Physicists in Medicine Radiation Therapy Committee Task Group 23: Radiation Treatment Planning Dosimetry Verification. American Institute of Physics; 1995.
Casanova Borca V, Pasquino M, Bresciani S, Catuzzo P, Ozzello F, Tofani S. Dosimetric evaluation of a commercial 3D treatment planning system using the AAPM Task Group 23 test package. Med Phys 2005;32:744-51.
Low DA, Harms WB, Mutic S, Purdy JA. A technique for the quantitative evaluation of dose distributions. Med Phys 1998;25:656-61.
Bewley DK. Central axis depth dose data for use in radiotherapy: A survey of depth doses and related data measured in water or equivalent media. Br Instit Radiol Suppl 1983;17:1-147.
Fraass B, Doppke K, Hunt M, Kutcher G, Starkschall G, Stern R, et al
. American Association of Physicists in Medicine Radiation Therapy Committee Task Group 53: Quality assurance for clinical radiotherapy treatment planning. Med Phys 1998;25:1773-829.
[Figure 1], [Figure 2], [Figure 3], [Figure 4]
[Table 1], [Table 2], [Table 3]