Database
Between June 2018 and October 2020 a total of 395 treatments have been planned with the Precision TPS v1.0.02 using the GPU-based optimizer VoLO™, which incorporates the Collapsed Cone Convolution Superposition dose calculation algorithm [1, 2]. The treatments were related to a large number of different sites and diseases; therefore, within the scope of this work they have been grouped in the following categories: abdomen, brain, head and neck (HN), thorax, pelvis, prostate and others. The dose per fraction (D/fr) ranged from 1.6 Gy to 5 Gy.
Plans Characteristics
To characterize the plans, the following parameters were collected: actual modulation factor (MF), pitch, field width (FW), gantry period, total delivery time, couch speed, couch travel distance and TTDF (ratio between the treatment time and the prescribed dose per fraction). Moreover, some descriptors of the leaf open time (LOT) distribution were collected: minimum (min-LOT), mean (mean-LOT), maximum (max-LOT) and LOT standard deviation (SD-LOT).
Qa Delivery And Analysis
All the patient-specific QA measurements were collected with ArcCHECK™ (Sun Nuclear, Melbourne, FL) and performed without the 15 cm diameter homogeneous PMMA plug. The plans were recomputed on the synthetic ArcCHECK™ CT provided by the vendor with an imposed density of 1.1836 g/cm3 following manufacturer's recommendations. Acquisition and analysis were performed with Sun Nuclear SNC Patient™ version 6.7. The absolute dose calibration of ArcCHECK™ was periodically controlled against a calibrated farmer-type ionisation chamber. Before each measurement session, the output of the Radixact unit was checked and a correction applied to the ArcCHECK™ analysis. QA plans were computed in high resolution mode which grants a 1.87 mm spaced dose grid.
The QA plan was computed placing the centre of the planning PTV at the machine’s isocenter. This procedure reduces the possibility for high dose areas to fall at the border of the ArcCHECK™ limiting the dependence of the gamma passing rate (GP%) on the maximum calculated dose [5].
According to the AAPM TG-218 report, the dose comparison was performed through 2D gamma analysis using a 10% dose threshold as well as 3% dose difference (DD) and 2 mm distance-to-agreement (DTA) criteria. In order to compare our results with previous works [5, 6, 7, 8], 3%/3mm gamma analysis was also collected. Both local and global normalizations were used.
Qa Program
Following the AAPM TG-218 report, the clinical deliverability of a plan was evaluated on the basis of the universal tolerance limit (TL) and action limit (AL), which are 95% and 90% respectively, on the GP% computed with (3%, 2mm), global normalization criteria.
Statistical process control techniques recommended in TG-218 have been employed to define a clinic specific TLcs and ALcs [4] as follows:
$${AL}_{cs}=100-3\sqrt{{\sigma }^{2}+{\left(\underset{\_}{x}-100\right)}^{2}}$$
$${TL}_{cs}=\underset{\_}{x}-2.660\underset{\_}{mR}$$
where \(\underset{\_}{x}\) is the GP% average over the investigation period, \(\sigma\) is its standard deviation and \(\underset{\_}{mR}\) the moving average. These values were computed on the first 40 collected measurements. The IMRT treatment process has been monitored and investigated by periodically computing these values to verify whether the process was under control. For the subsequent evaluations, periods of approximately 6 months were considered. Only clinically deliverable plans were considered in this analysis and 95% confidence intervals on TLcs and ALcs have been computed through a bootstrap approach with 10,000 samples.
Statistical analysis
To spot possible differences of GP% among the different treatment sites a Kruskal-Wallis test was performed followed by a post-hoc analysis where Bonferroni correction was applied.
N-way analysis of variance (n-way ANOVA) was performed to spot the influence of the planning parameters on GP%. Only parameters that strictly bore no correlations among themselves were included in the test (p > 0.05 and/or Pearson’s r < 0.8).
All tests were performed using MATLAB® R2020b (MathWorks, Inc., Natick, MA, USA). All the p-values reported are two-sided and p < 0.05 is considered as statistically significant.