Submit manuscript...
Journal of
eISSN: 2373-633X

Cancer Prevention & Current Research

Research Article Volume 13 Issue 1

Development of spreadsheet for rapid assessment of therapeutic radiation dose delivery with electron and photon beams at various energies

Md. Motiur Rahman,1 M Shamsuzzaman,2 MMH Bhuiyan,2 KA Khan,1 Sadiq R Malik3

1Department of Physics, Jagannath University, Dhaka, Bangladesh
2Bangladesh Atomic Energy Commission (BAEC), Atomic Energy Research Establishment (AERE), Bangladesh
3TMSS Cancer Center, Bangladesh

Correspondence: Md. Motiur Rahman, Department of Physics, Jagannath University, Dhaka, Bangladesh, Tel +88 01717821974

Received: December 20, 2021 | Published: January 31, 2022

Citation: Rahman MM, Shamsuzzaman M, Bhuiyan MMH, et al. Development of spreadsheet for rapid assessment of therapeutic radiation dose delivery with electron and photon beams at various energies. J Cancer Prev Curr Res. 2022;13(1):8-12. DOI: 10.15406/jcpcr.2022.13.00479

Download PDF

Abstract

Accurate radiation therapy, dose calculation and quality control is the fundamental concern and objective for Oncologists and Physicists in Cancer Treatments. Global use of Radiation Therapy Machines like, 60Co and LINAC (Old), demands an accurate and fast External Beam Radiation therapy (EBRT) machine dose evaluation, where treatment planning systems (TPS) are not available. Some of the cancer therapy centers do not have a medical physicist or a dosimetrist. For that, here we developed a spreadsheet named ”Mithu” where Radiation Oncologists, physicists, and Radiotherapy technologists can easily calculate treatment time (TT) or Monitor Unit (MU) by Using an easy Microsoft Excel spreadsheets presented here. To create a treatment planning system for 2D or Emergency plan and Dosimetry need to calculate Cobalt-60 (60Co) and LINAC Machine actual dose rate for source to surface distance (SSD) and Source to Axis Distance (SAD). Inside the spreadsheet Percentage depth dose (PDD), & Tissue Maximum Ratio (TMR), has been used from the British Journal of Radiology supplement 25 (BJR). Output factor, Wedge factor, and Tray factor have been used from machine commissioning data. The equivalent square field is another important and basic parameter for radiotherapy TT or MU calculation. Linac Machine dose rate is accurately calibrated by Quality Assurance procedures by Physicists which remains fixed, no decay. But 60Co dose rate always decrease or changes as yearly decay constant is 0.131 (ln2/HL) and source decays 1.09% per month. To develop this spreadsheet some important Excel functions like HLOOKUP, VLOOKUP, EXACT, and logical function if-Else is mostly used for radiotherapy patients treatment plan dose calculation which has shortened the man-made hand calculation error by 99%. Almost 2000 patients’ 60Co Machine-based Radiotherapy treatment plan has been calculated using this “Mithu” Spreadsheets. Here we observed that 60Co dose calculated is varied <1% where LINAC MU calculated is varying <1.5% which is an acceptable range in radiotherapy.

Keywords: EBRT, PDD, TMR, Dose, QA, MU, TT, treatment planning system, spreadsheet’s, physicist, oncologist

Introduction

The spreadsheet “Mithu” has been developed for a radiotherapy system where a treatment planning system was unavailable. Medical Physicist' time-saving in busy government and private radiotherapy centers in South Asia and the African region.1 Another goal is that Medical physicist is not available in a government hospital (Some Asian & African country) so Radiation Oncologist & technologist make a treatment dose calculation manually.2 So there have a probability to increase error in treatment dose. The manual dose calculation depends on some tables Like PDD, TMR, Output factor, etc. from the treatment machine. In BJR-25 table3-5 field size & depth are in integer form except d-max depth (field size: 4x4, 6x6, 8x8, 10x10, 12x12, 15x15 cm,......40x40 cm & Depth: 0.5, 1, 1.5, 2.4, 3, 4, 6, 8,10,12,15 cm...... 30 cm etc.). so if need to calculate a PDD for 7x7 cm field size & depth 4.3 cm its need to interpolate data from this table which is very time consuming procedure. So Spreadsheet “Mithu” will be the best solution for easy way treatment dose calculation. Another advantage of these Spreadsheets is that they can show the 25 days treatment dose schedule which is the same for LINAC all days schedule but 60Co it is very important due to the disintegration of the radioactive source.6-8 So 1st fraction dose to the 25th fraction dose will not be the same for 60Co. Every day its changes a little bit. So physicists don't need to correct the treatment dose every month for 60Co EBRT if they use it. Before apply the dose to the cancer patient’s one physicist needs to know how much dose exit from the machine at specific distance called source to surface distance. Most of the linear accelerator has SSD 100cm and Most of the 60Co machine has SSD 80 cm but SSD 100 60Co machine now a day available in market.9,10 After strike (Ƴ-ray or X-ray) a surface dose is accumulated at a point in ward from the surface is called D-Max. So it is important to know how much dose is accumulated per minute for 60Co source and per MU for LINAC in dose maximum (D-max) at fixed reference condition (field size 10x10 cm, SSD 80 or 100 cm, Chamber position 5 cm or 10 cm depth form water surface depend on Machine energy).11,12 TRS-398 & TG-51 have the details protocol for absolute dose measurement. Spreadsheet Mithu” have this good and easy solution followed by TRS-398.Also In modern radiotherapy system after installation a 60Co or Linac machine Medical physicist start beam data acquisitions work using 3D water phantom, 2D Water Phantom, Farmer type ionization chamber (volume =>0.6cc ), diode or small volume ionization chamber for small field dosimetry. After completed acceptance test its take few days for machine commissioning. So our spreadsheet “Mithu” will be the good tools /for verifying the TPS System. As per BJR Suplimentary-25 data table for LINAC Machine we know that, 6MV photon beam have 67% dose and 60Co have 56% dose at 10 cm depth for 10x10cm collimator jaw setting. So if we deliver a dose 200cGy at 10 cm depth than 298.5 cGy dose (Beam is calibrated 1 MU = 1cGy) need to deliver. So in TPS we can delivery same dose in water phantom at same depth and observed that, our TPS calculated dose near about 298.5 cGy. If the calculated dose vary more that 2% from our Spreadsheet “Mithu” than again need to check the commissioning data. Most of the duel energy Linac have Photon (X-ray) & Electron (e-) therapy treatment option.13 So electron therapy we need to know the cutout factor, energy verses collimator jaw setting, output factor and percentage depth dose which will be acquired at the commissioning time.14 Spreadsheet “Mithu” have this option to calculate dose for electron therapy. Patient’s database is the very important part in radiotherapy center. Patients Address, contact number, diagnosis and treatment dose assign doctors & physicist name need to record in database. So that easily physicist can say home much patients be treated per year or month. Can be analysis, which type of organ (site) is treated more and demographic basis which type of patients is more. Admin panel is the heart of this Spreadsheet “Mithu” After commissioning every date must be put inside the admin panel for Cobalt-60, LINAC and Electron dose calculation. This Spreadsheet “Mithu” is uploaded in “https://lnkd.in/gcrr-G-J”& “https://lnkd.in/gkF7csri” site which can be downloaded for use in any site shows in Figure 1.

Figure 1 Main panel of Spreadsheet “Mithu”.

Material and methods

I.A material: In the present study & development, essential tools are Microsoft Excel & Visual basic measurement from machine specific data absolute dosimetry, output factor, wedge factor, tray factor etc. PDD & TMR table used from BJR Suppl. 25(1996) both for 60Co & Linac (6 & 10 MV) beam. Some important function HLOOKUP (), VLOOKUP (), Exact () and logical function IF, Or, AND is mostly used. HLOOKUP (), VLOOKUP () we use here to extract data from the table. HLOOKUP searches for the value in a table row & VLOOKUP search for the value in a table column. HLOOKUP (), VLOOKUP () will get exactly the same result, but the table is arranged differently. EXACT function compare between two strings. IF function allows to make a logical comparison between a value like =IF (Something is True, then do something, otherwise do something else). AND, OR and NOT is an individual functions but we use it here as a combine an IF statement Example =IF (OR (Something is True, Something else is True), Value if True, Value if False). Microsoft Visual Basic is the tools which are built in with Microsoft office. Just need to active the developer option from file options menu inside the Microsoft excel. TRS-398 data and formula is also use for Absolute dosimetry. The objective of this work may, therefore, be taken as a reference for any user or, can create spreadsheets their own style for the machine treatment Time/MU verification of quality performance of the machine, while using it at any Cancer Treatment Center (CTC). Detailed information on materials and methods are compiled here.

I.B methods

Admin Panel: The admin panel (or the Administration Panel) is the primary tool for any kind of radiotherapy software. Though “Mithu” is an excel spreadsheet and we were developed a separate sheet inside the “Mithu” can say admin panel for input the basic machine-specific data before using these spreadsheets in their cancer therapy center. This admin sheet will be the interlink with others user sheets. So it will be very easy to use these spreadsheets if only once a time input the basic data Like Output factor, Wedge factor, Machine calibration value (1 MU = ?) or 60Co dose rate input . Evoke that, computer date and time is a very important and sensitive issue for the 60Co time calculation spreadsheet. Because radioactive source activity decay with the day (Figure 2).

Figure 2 Admin panel arrangement for the “Mithu” software.

Radiotherapy treatment plan: In house radiotherapy treatment planning system “Mithu” has been developed for Linac 6MV & 10MV photon energy emergency radiotherapy or 2D treatment plan for MU calculation, 60Co radiation therapy 2D treatment plan for treatment time calculation also can be easily cross check commercial treatment planning systems. Inside the treatment planning excel Spreadsheets “List” function were used from Data Menu for selectable data. Example: In our Linac Spreadsheets we can change energy and technique just press on energy cell and select energy, same for treatment technique, Tx Area, Machine, Physicist, and oncologist cells.

It’s has an advantage is that:

  1. Its can use both SSD & SAD method treatment calculation. Just change the Treatment Technique option.
  2. Its can make treatment plan schedule as per treatment fraction.
  3. No need manual data sheet table. It’s have the intelligence to find out PDD, TMR, Equivalent Square Field (ESF)3
  4. No need decay correction in patient’s treatment card for 60Co EBRT patients.
  5. All kind of correction factor can be associated with TT or MU (Figure 3,4).

Figure 3 Linac 2D treatment planning Spreadsheets which can be used for 6 & 10 MV photon dose calculation both for SSD &SAD method.  

Figure 4 60Co treatment planning Spreadsheets. Can perform treatment time or dose calculation.

Basic Equation for Manual dose calculation is used in a center is:

 For Cobalt-60 (60Co )

TT= Prescribeddose Doserate×%DDorTMR×Outputfactor MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=Mj0xXdbba91rFfpec8Eeeu0xXdbba9frFj0=OqFf ea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadsfacaWGub Gaeyypa0ZaaSaaaeaaciGGqbGaaiOCaiaadwgacaWGZbGaam4yaiaa dkhacaWGPbGaamOyaiaadwgacaWGKbGaaGPaVlaadsgacaWGVbGaam 4CaiaadwgaaeaacaWGebGaam4BaiaadohacaWGLbGaaGPaVlaadkha caWGHbGaamiDaiaadwgacqGHxdaTcaGGLaGaamiraiaadseacaaMc8 Uaam4BaiaadkhacaaMc8Uaamivaiaad2eacaWGsbGaey41aqRaam4t aiaadwhacaWG0bGaamiCaiaadwhacaWG0bGaaGPaVlaadAgacaWGHb Gaam4yaiaadshacaWGVbGaamOCaaaaaaa@6C02@  (1)

 For Linac

MU= Prescribeddose cGy MU ×%DDorTMR×Outputfactor MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=Mj0xXdbba91rFfpec8Eeeu0xXdbba9frFj0=OqFf ea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaad2eacaWGvb Gaeyypa0ZaaSaaaeaaciGGqbGaaiOCaiaadwgacaWGZbGaam4yaiaa dkhacaWGPbGaamOyaiaadwgacaWGKbGaaGPaVlaadsgacaWGVbGaam 4CaiaadwgaaeaadaWcaaqaaiaadogacaWGhbGaamyEaaqaaiaad2ea caWGvbaaaiabgEna0kaacwcacaWGebGaamiraiaaykW7caWGVbGaam OCaiaaykW7caWGubGaamytaiaadkfacqGHxdaTcaWGpbGaamyDaiaa dshacaWGWbGaamyDaiaadshacaaMc8UaamOzaiaadggacaWGJbGaam iDaiaad+gacaWGYbaaaaaa@6780@  (2)

Normally the dose rate for the Linear accelerator is 1cGy/MU for 10x10 cm field size defined at the iso-center, Percentage depth dose is depend on the equivalent square field (ESF) and the prescribed depth. Output factor is Total scattering factor which is measured at machine commissioning time. Some others factor we can also use like Wedge and Tray factor depends on organ shielding and beam hinge angle. In this developed software we use this equation (1) and (2) for treatment Time and Monitor Unit (MU) calculation.For electron therapy we use a different size of Applicator and cone cutout factor for different energy. Most of the Linac in cancer center have 4 to 6 electron energy step. Here we develop this spreadsheet using VARIAN machine which have an energy step 6 , 9, 12, 15, 18 & 22 MeV among of them first four energy steps (6,9,12,15 MeV) has been used most of the clinical treatment purpose rest of the two energy steps (18, 22 MeV) seldom use in clinical treatment purpose. Varian machine has an applicator cone size 6x4, 6x6, 10x10, 15x15, 20x20, 25x25cm. Collimator jaw setting depends on cone size vs machine physics setting. Here we use VARIAN common setting for iX machine data setup in Delta Hospital Ltd. Mirpur-1, Dhaka. Different field size cutout can be made in house for clinical use and output factor can be determined as per energy. All of the measured date Applicator size, Cutout factor we can put inside the admin sheet in electron spaces. HLOOKUP () is the main excel function in this spreadsheet. Its can horizontally search the data from the table as per selection electron energy, cone cutout (Table 1).

 

6x6

10x10

15x15

20x20

25x25

6

20x20

20x20

20x20

25x25

30x30

9

20x20

20x20

20x20

25x25

30x30

12

12x12

14x14

17x17

25x25

30x30

15

11x11

14x14

17x17

23x23

28x28

18

         

22

 

 

 

 

 

Table 1 Electron Applicator vs Collimator Jaw setting. From the spreadsheet admin panel anyone can change as per their machine setting

Results & discussion

Figure 5 & 6 shows the applications of in house Spreadsheets (name “Mithu”) from the year 2012 for radiotherapy patients treatment plan dose calculation which has shortened the man-made hand calculation error by 99%. Almost 2000 patients’ 60Co Machine based Radiotherapy treatment plan has been calculated using this Mithu software at Rajshahi Medical College and Hospital, Rajshahi, National Institute of Cancer research Hospital, Dhaka, Dhaka Medical College and Hospital, Dhaka and Delta Hospital Ltd., Dhaka for a precise treatment, Hospital Physicist Association,1 Central Axis Depth dose data for radiotherapy (BJR Suppliment-11) published an equivalent Square of rectangular fields table.

Figure 5 60Co treatment time schedule depends on the fraction. As per plan & field it can able to visualize the treatment information. Be used for 6 & 10 MV photon dose.

Figure 6 Shows Graphical representation between Hospital Physicist Association and sterling equation of deviation.

Table 2 and 3 show that the deviation of calculated MU is less than ±1.5% for 6 and 10 MV photon beam between the Commercial TPS system and “Mithu” spreadsheet. Figure 7 shows the Sterling et al.,2 equation and Hospital Physicist Association table3 comparative deviation between equivalent square field. The sterling equation is almost matched with Hospital Physicist Association table for small field size <10 cm but it can vary for >10cm. Example 12cm X 24 cm rectangular field size as per sterling equation its equivalent square field size is 16 cm but Hospital Physicist Association table says 15.4 cm so deviation is 6 mm. or 3.9%. Below Graph1: shows Graphical representation between Hospital Physicist Association and sterling equation deviation. So its deviations have little bit impact on TT/MU calculation due to change of PDD & TMR value but are negligible. Using both BJR Table and Sterling’s Equation all the equivalent Square field size output factors were not very too much.

Treatment Depth (cm)

Eclipse Planning System Data

“Mithu” Planning System Data & Result

Deviation  ±∆%

% DD

Calculated MU

% DD

Calculated MU

5

86.55

231

86.9

230

0.41

10

67

299

67.5

296

0.85

15

51.2

391

51.7

387

0.94

20

38.9

514

39.3

509

1.01

Table 2 6MV, 10x10 cm field size has been calculated at various depths using Commercial TPS system and “Mithu” spreadsheet and their deviation

Treatment Depth (cm)

Eclipse Planning System Data

“Mithu” Planning System Data & Result

Deviation  ±∆%

% DD

Calculated MU

% DD

Calculated MU

5

91.6

218

91.4

219

-0.30

10

73.5

272

73

274

-0.69

15

58.2

344

57.8

346

-0.68

20

46

435

45.6

439

-0.96

Table 3 10MV, 10x10 cm field size has been calculated at various depths using Commercial TPS system and “Mithu” spreadsheet and their deviation

Figure 7 Output Variation with fixed Y and Variable X Jaw for (3 X xi) 6MV Photon energy.

Figure 8 shows the experimental result we got the output deviation between these two methods is maximum ±0.45% for 60Co ±0.4% for 6MV & ±0.23% for 10MV which is not much impact on clinical dose calculation for radiotherapy treatment. Grah2 & 3 Shows the output variation for X & Y jaw setting. For PDD and TMR we got some deviation between BJR and Machine commissioning data which is not more than 2%. As per ICRU50 +7% and -5% is the acceptable dose for the prescribed dose.

Figure 8 Output Variation with fixed Y and Variable X Jaw for (5 X xi) 6MV Photon energy.

Conclusion

Modern radiation therapy quality assurance (QA) devices and treatment planning systems are available in the market. But in some government hospitals or Radiation therapy centers in South Asia and Some African countries use 60Co and LINAC Machine where TPS is not available and Medical Physicist is not appointed. The developed spreadsheet can be used for the purpose of quality and accurate Dose delivery in radiation therapy in addition to QA purposes. It’s also the cost-effective solution for the center where a commercial planning system is not available. Assurance of quality radiation therapy is presented here. It will be a good solution where, physicists and dosimetrist are not available. So that, error-free treatment may be delivered. Some busy centers will be helpful for saving the physicist's time and shortening patients waiting time for radiotherapy card preparation and calculation.

Compliance with ethical standards

Acknowledgments this research work did not receive grant from any organization. Nevertheless, the technical support was afforded by Department of Physics, Jagannath University, Dhaka throughout the execution of this research work. Disclosure of conflict of interest there is no conflict of interest declared on this research article.

Acknowledgments

Our sincere thanks to Rajshahi Medical College & Hospital, Delta Hospital Ltd., in providing the Machine-specific data and use these spreadsheets for their clinical use purpose. Thanks to Prof. Dr. Sadiq R Malik, Radiation Oncology Physicist from USA for his encouragement and association respectively, to work on 60Co & Linac machines described here. Thanks are also due to Professor Dr. Parvin A. Banu, M.D., Sr. Consultant Oncologist, Labaid cancer center and Professor Dr. AKM Ahsan Habib, DMRT, Head of the Oncology Department, TMSS Medical College and Hospital, for their support in the clinical use of the Teletherapy Units, for cancer treatments.

Conflicts of interest

Authors declare that there is no conflict of interest.

References

  1. Akulapalli Sudhakar. History of Cancer, Ancient and Modern Treatment Methods. J Cancer Sci Ther. 2009;1(2):1–4.
  2. Baskar R, Lee KA, Yeo R, et al. Cancer and Radiation Therapy: Current Advances and Future Directions. Int J Med Sci. 2012;9(3):193–199.
  3. Central Axis Depth Dose Data for Use in Radiotherapy. A survey of depth doses and related data measured in water or equivalent media. Br J Radiol Suppl. 1983;17:1–147.
  4. Sadiq RM, Rahman M, Mia M, et al. Comprehensive Benchmark Procedures for Commissioning, Quality Assurance and Treatment Delivery for Quality Cancer Therapy Using Co60 Teletherapy   Bangladesh Journal of Medical Physics. 2013;6(1):1–20.  
  5. Motiur Rahman, Shamsuzzaman M, Sarker M, et al. Dosimetric characterization of medical linear accelerator Photon and Electron beams for the treatment accuracy of cancer patients. World Journal of Advanced Engineering Technology and Sciences. 2021;3(1):41–59.
  6. IAEA TRS-398: Absorbed dose determination in external beam radiotherapy. An International code of practice for dosimeter; 2000. 242 p.
  7. Das IJ, Ding GX, Ahnesjö A. Small fields: Non-equilibrium radiation dosimetry. Med Phys. 2008;35(1):206–215.
  8. Scott AJ, Nahum AE, Fenwick JD. Using a Monte Carlo model to predict dosimetric properties of small radiotherapy photon fields. Med Phys. 2008;35(10):4671–4684.
  9. Henry FG, Ravikumar M, Sathiyan S, et al. Analysis of small field percent depth dose and profiles: Comparison of measurements with various detectors and effects of detector orientation with different jaw settings. Journal of Medical Physics. 2016;41(1):12–20.
  10. Tiru S Subramanian, Gibbons JP, William RH. Moderator Linear accelerators used for IMRT should be designed as small field, high intensity, intermediate energy unit. Med Phys. 2002;29(11):2526–2528.
  11. Bethesda MD. International Commission on radiation Units and Measurements (ICRU), Determination of Absorbed Dose in a Patient Irradiated by Beams of X or Gamma rays in Radiotherapy Procedures; 1976. 24 p.
  12. Gerald J Kutcher, Lawrence Coia, Michael Gillin, et al. Comprehensive QA for radiation oncology: Report of AAMP radiation therapy committee. Task Group-40 (TG-40) protocol of American Association of Physicists in Medicine (AAPM). Med Phys. 1994;21(4):41.
  13. Eric E Klein, Joseph Hanley, John Bayouth, et al. Quality assurance of medical accelerators. Task Group-142 (TG-142) protocol of American Association of Physicists in Medicine (AAPM). Med Phys. 2009;36(9).
  14. Sahoo SK, Rath AK, Mukharjee RN, et al. Commissioning of a Modern LINAC for Clinical Treatment and Material Research. International Journal of Trends in Interdisciplinary Studies. 2012;10.
Creative Commons Attribution License

©2022 Rahman, et al. This is an open access article distributed under the terms of the, which permits unrestricted use, distribution, and build upon your work non-commercially.