 Research
 Open access
 Published:
Simplified sigmoidal curve fitting for a 6 MV FFF photon beam of the Halcyon to determine the field size for beam commissioning and quality assurance
Radiation Oncology volume 15, Article number: 273 (2020)
Abstract
Background
An Oring gantrytype linear accelerator (LINAC) with a 6MV flattening filterfree (FFF) photon beam, Halcyon, includes a reference beam that contains representative information such as the percent depth dose, profile and output factor for commissioning and quality assurance. However, because it does not provide information about the field size, we proposed a method to determine all field sizes according to all depths for radiation therapy using simplified sigmoidal curve fitting (SCF).
Methods
After mathematical definition of the SCF using four coefficients, the defined curves were fitted to both the reference data (RD) and the measured data (MD). For good agreement between the fitting curve and the profiles in each data set, the field sizes were determined by identifying the maximum point along the third derivative of the fitting curve. The curve fitting included the field sizes for beam profiles of 2 × 2, 4 × 4, 6 × 6, 8 × 8, 10 × 10, 20 × 20 and 28 × 28 cm^{2} as a function of depth (at 1.3, 5, 10 and 20 cm). The field size results from the RD were compared with the results from the MD using the same condition.
Results
All fitting curves show goodness of fit, R^{2}, values that are greater than 0.99. The differences in field size between the RD and the MD were within the range of 0 to 0.2 cm. The smallest difference in the field sizes at a depth of 10 cm, which is a surfacetoaxis distance, was reported.
Conclusion
Application of the SCF method has been proven to accurately capture the field size of the preconfigured RD and the measured FFF photon beam data for the Halcyon system. The current work can be useful for beam commissioning as a countercheck methodology to determine the field size from RD in the treatment planning system of a newly installed Halcyon system and for routine quality assurance to ascertain the correctness of field sizes for clinical use of the Halcyon system.
Background
The Halcyon (Varian Medical Systems Inc., Palo Alto, USA), which is a linear accelerator (LINAC) with an Oring gantry, is a radiotherapy machine with a 6MV flattening filter free (FFF) beam, such as the Tomotherapy and Cyberknife systems. In the case of the Halcyon, because the commissioning and quality assurance (QA) processes are totally different than those of the conventional LINAC process, an approach has been developed for the users. The conventional LINAC process required strict and longtime measurements for the commissioning and QA steps, while the Halcyon system provides an independent opportunity to verify the consistency of the measured data (MD) during the commissioning and QA processes with the reference data (RD) provided by the vendor. The RD includes parameters such as the percentage depth dose (PDD), dose profile, and output factor, and the new commissioning process of the Halcyon generally follows the guidelines of Association of Physicists in Medicine (AAPM) MMPG 5.a, AAPM TG51 and TG100 to satisfy international commissioning standards [1,2,3]. Therefore, many studies have checked and verified various parameters other than the basic parameters for the commissioning and QA processes of the Halcyon, and these parameters could be applied as another factor to verify the accuracy of the commissioning and QA processes of the Halcyon [4, 5]. The purpose of this study is to demonstrate the method for determining the field size of the Halcyon system using the simplified sigmoidal curve and to provide a field size parameter dataset that can improve the effectiveness of the commissioning and QA processes.
In radiotherapy, when using a LINAC, the field size of the radiation beam refers to the area of radiation delivery. For this reason, determination of accurate field sizes is a significant parameter for the delivery of radiation and an important process in commissioning and quality assurance (QA) [6, 7]. Normally, the field size can be determined during a commissioning and QA procedure [8, 9]. While conventional LINAC systems are often equipped with a flattening filter (FF) to deliver the radiation beam with a uniform dose distribution, the full width at half maximum (FWHM) method is the conventional representative methodology for determining the field size of the FF beam [10]. The field size of the FF beam is defined based on a point offaxis at a dose of 50% after dose normalization of a central axis (CAX) at 100%. This FWHM methodology is suitable for determination of the field size of the FF beam because it has a uniform region around the CAX. According to the Task Group Report #142 of AAPM, several parameters, such as flatness, symmetry and penumbra, should be considered when characterizing the FF beam [11]. Currently in radiation therapy, the accuracy of the dose is a very important factor, and the treatment time should not be overlooked. Because a long treatment time may cause patient discomfort and decrease the accuracy of treatment, a flattening filterfree (FFF) beam is used to reduce the treatment time. The FFF beam has the effect of reducing photon head scatter, leaf transmission head leakage and the peripheral dose [12, 13]. However, the FWHM method is not suitable for determining the field size of the FFF beam, which shows a specific shape of the dose profile, with a relatively higher peak at the CAX.
Several studies have defined the field sizes of an FFF beam. The most representative method to determine the field size of an FFF beam is to use the inflection point (IP) on the penumbra region of the beam’s profile [14, 15]. Nevertheless, some uncertainties remain in obtaining a correct IP from beam data measurements. To consider this uncertainty, Pönisch et al. proposed a method to identify the IP at the field edge of an FFF beam with the same level as that of an FF beam. The position of the IP can be changed according to the positional step size error to obtain the beam profile [16]. Fogliata et al. suggested a renormalization formulation to overcome the uncertainty in the IP due to the stepping size [17]. Although these two methods to define the field size of the FFF beam are based on the profiles of the FF beam and thus include a large number of values, they both exhibit a position error during measurement. The parameterized gradientbased method (PGM) that complements these two methods was proposed to determine the field size of the FFF beam using a mathematical model. Although the PGM results were effective for determining the field size, this method did not yield parameters for the specific field size. We struggled to determine the reason why some data were omitted. Some institutes may want the omitted data according to the level of an instrument or the requirements of a specific treatment case [18, 19].
In addition, because the PGM method also applies a mathematical model based on MD, uncertainty still exists in the measurements, as in the initial two methods. When users of the Halcyon use our proposed method to determine the field size, there are several convenient benefits. First, we proposed a reasonable method to determine the field size of the Halcyon beam according to all depths and all field sizes. The vendor of the Halcyon has already entered all beam data into the Eclipse TPS. This RD is used to optimize and calculate the dose, and this approach is effective for operation of the Halcyon in the Eclipse TPS. However, because the vendor does not provide any information regarding the field size in the RD, the field size should be determined by their method for all conditions. In this study, one of the methods for determining the effective value of the field size with a simple mathematical equation was employed for all depths and all field sizes.
Second, we used the RD from the vendor to determine the field size and dose without using our MD. The commissioning process of the Halcyon does not require an adjustment of the beam model by the user in the way that conventional LINAC commissioning does. Instead, the user should check the degree of correspondence between the RD included in the installed TPS and the MD to confirm the criteria of commissioning. This can reduce the possibility of errors induced by unexpected conditions during commissioning and the QA process in the conventional LINAC, and it can provide an opportunity to verify the reliability and accuracy of beam data by comparing the degree of correspondence between the RD and the MD. When we measure the data for the additional commissioning procedure or QA steps, we can check the accuracy of the field size using the reference value from the RD.
Finally, we used the methodology with only the sigmoid function without any other equations to generate the fitting curve on the beam profile. The procedure for the proposed method that is applicable to the RD has been simplified compared with the procedure of the PGM.
Methods
Preparation of data
A preconfigured reference beam dataset (RD) generated by the vendor is stored in the treatment planning system when a new Halcyon system is installed. The RD includes the lateral dose profiles for field sizes of 2 × 2, 4 × 4, 6 × 6, 8 × 8, 10 × 10, 20 × 20 and 28 × 28 cm^{2} as a function of depth at 1.3, 5, 10 and 20 cm in the water phantom. To compare the field size from the RD with that from the MD, the measurement was performed under the same conditions as those used to obtain the RD. The sourcetosurface distance (SSD) was set at 90 cm. A CC13 ionization chamber and a Blue Phantom water tank (IBA Dosimetry, Schwarzenbruck, Germany) were used to measure the relative dose profiles for field sizes > 4 × 4 cm^{2}. For field sizes ≤ 4 × 4 cm^{2}, an edge diode detector (Sun Nuclear, Melbourne, FL, USA) was used. The scanning step for acquisition of the profile on the measurement line along the offaxis position was 0.1 cm. All measurement values were processed with OmniPro Accept7 (version 7.4.24.0) software (IBA dosimetry, Schwarzenbruck, Germany).
Definition of fitting using sigmoidal curve
The sigmoidal curve originates from the sigmoid function, which is used in the field of the signal process. The shape of the sigmoidal curve is given by Eq. (1),
The coefficients α, β, γ and δ are used to determine the shape of the curve f(x). The coefficient α controls the gradient of the sigmoidal curve. The higher the value of α is, the steeper the curve gradient. The coefficient β is related to the horizontal movement of the entire sigmoidal curve. The higher the value of β is, the further the sigmoidal curve moves to the right. The coefficient γ determines the location of the only upper end of the sigmoidal curve. The higher the value of γ is, the higher the position of the upper end of the sigmoidal curve. The coefficient δ determines the vertical movement of the entire sigmoidal curve. The higher the value of δ is, the more upward the direction of the sigmoidal curve. Thus, the coefficients α and γ contribute to transforming the shape of the curve. The coefficients β and δ change the location of the sigmoidal curve.
After uploading the profile to MATLAB (2019 version, MathWorks Inc, Sherborn, MA, USA), sigmoidal curve fitting (SCF) was performed by changing each coefficient until the sigmoidal curve overlapped the profile. To fit the sigmoidal curve to the profile, four steps are required. First, the gradient level of the sigmoidal curve should be the same as the gradient of the profile. When these two gradients agree, the coefficient of α is obtained. The second step is to let the sigmoidal curve move so that it overlaps the profile through the horizontal pitch using the coefficient β. In the third step, the coefficient γ is altered to adjust the upper limitation of the sigmoidal curve so that the curve stops at the upper end of the profile. Similarly, in the last step, the coefficient δ is altered to adjust the lower limitation of the sigmoidal curve so that it stops at the lower end of the profile. In this process, we are able to change all the coefficients to edit the shape of the sigmoidal curve fitting. If the fitting curve is normal, an Sshaped curve that exactly overlaps the RD curve will be obtained.
Verification of agreement for fitting curves
To verify the accuracy of the fitting curve based on the sigmoidal curve with the profiles, the average agreement ratio (AAR) between the values in the fitting curve (f_{i}) and the values in the profiles (x_{i}) at the same step position was calculated using Eq. (2), which shows the agreement between the profile and fitting values.
In this study, if the AAR is higher than 97%, the fitting optimization terminates because a sufficient accuracy has been obtained, and the four coefficients (α, β, γ and δ) are used to define the shape of the final fitting curve. Moreover, an additional verification was performed based on the evaluation of goodness of fit, R^{2} (Eq. 3).
where \(\overline{f}\) is the mean of all f_{i} values on the fitting curve, and y_{i} is a value on the profile. The same validation procedures were applied to the MD and the RD.
Identification of specific regions and points
In this study, to describe the sigmoidal curve, three regions and two points were assigned in the definition of the halfside of the SCF (Fig. 1a). The three regions include the introductory region (IR), the growing region (GR) and the plateau region (PR). The IR is the region where the sigmoidal curve begins to increase. The GR is a continuously increasing region on the sigmoidal curve. The PR is the region where the increase slows. These regions can be identified through the second derivative of the sigmoidal curve, as shown Fig. 1a. The range between the rightmost point and the maximum point on the second derivative curve is defined as the IR. The range between the maximum point and the minimum point on the second derivative curve is defined as the GR. The region between the minimum point and the leftmost point on the second derivative curve is defined as the PR. Because of the specification of the sigmoidal shape, there are two specific points: the singular point (SP) and the IP, both of which can be identified from the third derivative curve of the sigmoidal curve. The SP is the minimum point between the range of the IR and the GR (Eq. 4). The IP is another minimum point in the range between the GR and the PR, as shown in Eq. (5). When there is no point in either the IP or the SP, the refitting process from Sect. 2.2 is performed.
Determination of the field size
Although the determined field size (DFS) can be calculated using only the first derivative, the second and third derivatives provide an opportunity to check whether or not a given DFS exists in the period between the IP and the SP and to assess the error of the fitting curve. Incorrect fitting or an insufficient fitting range for the first derivative can be defined as the wrong field size. In the extreme case, the IP and the SP cannot be found on the derivative curve. As a result, the DFS cannot be calculated. If we determine the field size using only the first derivative, the results will contain some uncertainties. Therefore, the purpose of the second derivative is to check for the presence of both the IP and the SP. Finally, the purpose of the third derivative is to determine the DFS for the clear periods between the IP and the SP. Because the right third derivative curve can show these three points at the same time, determining the field size using the third derivative curve is the most efficient approach. After the SP and the IP have been obtained, the DFS can be identified as the maximum point on the third derivative curve between the SP and the IP, as shown below in Eq. (6). Figure 1b shows the conceptual DFS on the third derivative curve, and an actual example of the DFS is presented by fitting the profile in Fig. 1b.
The factor 2 is used in Eq. (6) because only the right half of the symmetric open beam profile was used for the curve fitting.
Results
Accuracy of the fitting curve
Figure 2 shows all the final fitting curves using SCF with half of the profiles in the RD. The profiles for all field sizes (2 × 2, 4 × 4, 6 × 6, 8 × 8, 10 × 10, 20 × 20 and 28 × 28 cm^{2}) are demonstrated as a function of depth (1.3, 5, 10 and 20 cm) in (a), (b), (c) and (d), respectively. The fitting curve is indicated by the red circles on each profile. In contrast, Fig. 3 shows all the final fitting curves using the same SCF with half of the profiles in the MD. This figure also includes the profiles for all field sizes with variable depths: depths of 1.3, 5, 10 and 20 cm are shown in (a), (b), (c) and (d), respectively. The profiles in Figs. 2 and 3 along the field size are distinguished by their color and have been normalized according to the relative dose at the CAX. The Xaxis shows the offaxis position from the CAX. All the fitting curves show good agreement with each profile. The accuracy of all the fitting curves was evaluated using the AAR and R^{2} methods. The R^{2} values were all greater than 0.99.
Table 1 tabulates all coefficients to form the final sigmoidal curves for all field sizes with all depths using SCF before determination of the field size. The differences of the coefficient values at all depths were reported as a maximum difference of α of 1.3, a maximum difference of β of 4.4, a maximum difference of γ of 4.0, a maximum difference of δ of 1.7, and all coefficients had a minimum difference of 0. All IP and SP values were also tabulated. The maximum and minimum IP differences were 0.15 cm and 0 cm, respectively. The maximum and minimum SP differences were 0.18 cm and 0 cm, respectively.
Table 2 lists all DFS values from the fitting curve for all the profiles at different depths. For the field sizes, the maximum difference was 0.2 cm, and the minimum difference was 0 cm. The results for a depth of 10 cm showed the least difference between the RD and the MD.
Discussion
Previous studies have been conducted to determine the field size of the conventional FF beam using methods such as the FWHM, which is not suitable for the LINAC profile in the Halcyon due to the difference in the CAX of the dose. The limitation of the current measurement is that it is time consuming to apply a fine scanning step of 0.1 cm in the QA and commissioning processes. MSM, techniques developed by Pönisch et al. and the method of renormalization, all of which are typical methods for defining the FFF beam field size, can contain uncertainty as a result of the scanning step [14,15,16,17]. The uncertainty in the location of the IP may occur from the MD due to the use of a scanning step with a lower resolution. In contrast, although the PGM method has the advantage of not affecting the size of the scanning step, its analytical fitting procedure also obtains the parameters from the MD. In the case of PGM, the field size of the MD is largely separated as 4 × 4 cm^{2} ≤ and 4 × 4 cm^{2} > to obtain the coefficient from the MD. However, the coefficients for 2 × 2 cm^{2} and 2.5 × 2.5 cm^{2} in the calculated data, which are used to check the accuracy of the coefficients of the MD, are not presented. Although this does not rule out the possibility of equipmentspecific changes, the field size of the small field is an important factor for equipment from any vendor [18, 19]. Thus, all the approaches described above define the field size based on data measured by the user, whereas our method defines the field size using the RD, which represents the characteristics of the Halcyon. This means that the consistency of the beam data is guaranteed for all institutions using the Halcyon, and the use of mathematical fitting curves to determine Halcyon field size ensures a high accuracy and reliability; further, commissioning and QA can be conducted at each institution based on the RD [2,3,4,5]. In our study, the field size of the FFF beam was defined in a simpler and more intuitive way than the PGM method, and the process was simplified while reducing the possibility of errors in the commissioning and QA steps by using the RD. These results are meaningful because a reference criterion is provided to Halcyon users that can be used to compare results with a full dataset. Adequate parameter information and a simple methodology could be useful, especially for new Halcyon users who must validate the preconfigured reference for the beam commissioning process and perform QA because beam data can vary from machine to machine, even for the same model and vendor, as described by AAPM MMPG 5.a, TG 51, 100, and 106 [1,2,3, 20].
Conclusions
The field size for a 6MV FFF radiation beam from the Halcyon system was determined using the simple SCF method. This method covers all field sizes, including small field sizes. The coefficients for the fitting and the field sizes between the RD and the MD were in good accord. This method can be used as a repeated countercheck for users using the same LINAC model.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Abbreviations
 QA:

Quality assurance
 IMRT:

Intensity modulated radiation therapy
 VMAT:

Volumetric modulated arc therapy
 SCF:

Sigmoidal curve fitting
 RD:

Reference data
 MD:

Measured data
 AAR:

Average agreement ratio
 IR:

Introductory region
 GR:

Growing region
 PR:

Plateau region
 SP:

Singular point
 IP:

Inflection point
 DFS:

Determined field size
References
Smilowitz JB, Das IJ, Feygelman V, Fraass BA, Kry SF, Marshall IR, et al. AAPM medical physics practice guideline 5. A.: commissioning and QA of treatment planning dose calculations—megavoltage photon and electron beams. J Appl Clin Med Phys. 2015;16:14–34.
Lloyd SA, Lim TY, Fave X, FloresMartinez E, Atwood TF, Moiseenko V. TG51 reference dosimetry for the HalcyonTM: a clinical experience. J Appl Clin Med Phys. 2018;19:98–102.
Teo PT, Hwang MS, Shields W, Kosterin P, Jang SY, Heron DE, et al. Application of TG100 risk analysis methods to the acceptance testing and commissioning process of a Halcyon linear accelerator. Med Phys. 2019;46:1341–54.
De Roover R, Crijns W, Poels K, Michiels S, Nulens A, Vanstraelen B, et al. Validation and IMRT/VMAT delivery quality of a preconfigured fastrotating Oring linac system. Med Phys. 2019;46:328–39.
Netherton T, Li Y, Gao S, Klopp A, Balter P, Court LE, et al. Experience in commissioning the halcyon linac. Med Phys. 2019;46:4304–13.
Kry SF, Dromgoole L, Alvarez P, Leif J, Molineu A, Taylor P, et al. Radiation therapy deficiencies identified during onsite dosimetry visits by the imaging and radiation oncology core Houston quality assurance center. Int J Radiat Oncol Biol Phys. 2017;99:1094–100.
Sun B, Goddu SM, Yaddanapudi S, Noel C, Li H, Cai B, et al. Daily QA of linear accelerators using only EPID and OBI. Med Phys. 2015;42:5584–94.
Hrbacek J, Lang S, Klöck S. Commissioning of photon beams of a flattening filterfree linear accelerator and the accuracy of beam modeling using an anisotropic analytical algorithm. Int J Radiat Oncol Biol Phys. 2011;80:1228–37.
Shende R, Gupta G, Patel G, Kumar S. Commissioning of TrueBeam TM Medical Linear Accelerator: Quantitative and Qualitative Dosimetric Analysis and Comparison of Flattening Filter (FF) and Flattening Filter Free (FFF) Beam. IJMPCERO. 2016;5:51.
Chang KP, Wang ZW, Shiau AC. Determining optimization of the initial parameters in Monte Carlo simulation for linear accelerator radiotherapy. Radiat Phys Chem. 2014;95:161–5.
Klein EE, Hanley J, Bayouth J, Yin FF, Simon W, Dresser S, et al. Task Group 142 report: quality assurance of medical acceleratorsa. Med Phys. 2009;36:4197–212.
Cho W, Kielar KN, Mok E, Xing L, Park JH, Jung WG, et al. Multisource modeling of flattening filter free (FFF) beam and the optimization of model parameters. Med Phys. 2011;38:1931–42.
Stieler F, Fleckenstein J, Simeonova A, Wenz F, Lohr F. Intensity modulated radiosurgery of brain metastases with flattening filterfree beams. Radiother Oncol. 2013;109:448–51.
Muralidhar K. Derivation of equations to define inflection point and its analysis in flattening filter free photon beams based on the principle of polynomial function. Int J Cancer Ther Oncol. 2015;3:03015.
Pichandi A, Ganesh KM, Jerin A, Balaji K, Kilara G. Analysis of physical parameters and determination of inflection point for Flattening Filter Free beams in medical linear accelerator. Rep Pract Oncol Radiother. 2014;19:322–31.
Pönisch F, Titt U, Vassiliev ON, Kry SF, Mohan R. Properties of unflattened photon beams shaped by a multileaf collimator. Med Phys. 2006;33:1738–46.
Fogliata A, Garcia R, Knöös T, Nicolini G, Clivio A, Vanetti E, et al. Definition of parameters for quality assurance of flattening filter free (FFF) photon beams in radiation therapy. Med Phys. 2012;39:6455–64.
Lebron S, Lu B, Yan G, Kahler D, Li JG, Barraclough B, et al. Parameterization of photon beam dosimetry for a linear accelerator. Med Phys. 2016;43:748–60.
Lebron S, Yan G, Li J, Lu B, Liu C. A universal parameterized gradientbased method for photon beam field size determination. Med Phys. 2017;44:5627–37.
Das IJ, Cheng CW, Watts RJ, Ahnesjö A, Gibbons J, Li XA, et al. Accelerator beam data commissioning equipment and procedures: report of the TG106 of the Therapy Physics Committee of the AAPM. Med Phys. 2008;35:4186–215.
Acknowledgements
Not applicable.
Funding
This research was supported by the National Research Foundation of Korea Grant Funded by the Korean Government (2018A2B2005343).
Author information
Authors and Affiliations
Contributions
Study concept and design: MC, DY, MT, HM, and TS. Data acquisition: MT, KM, MO, and HM. Data analysis and interpretation: MC, DY, ML, MT, KM, MO, MK, SKD, HM, and TS. Computer program modeling: MC, DY, and MK. Manuscript preparation and editing: MC, ML, DY, MT, HM, and TS. Manuscript review: All authors read and approved the final manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Choi, MG., Law, M., Yoon, DK. et al. Simplified sigmoidal curve fitting for a 6 MV FFF photon beam of the Halcyon to determine the field size for beam commissioning and quality assurance. Radiat Oncol 15, 273 (2020). https://doi.org/10.1186/s1301402001709x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s1301402001709x