Analysis of clinical patient-specific pre-treatment quality assurance with the new helical tomotherapy platform, following the AAPM TG-218 report

Purpose

This study presents patient-specific quality assurance (QA) results from the first 395 clinical cases for the new helical TomoTherapy® platform (Radixact) coupled with dedicated Precision TPS.

Methods

The passing rate of the Gamma Index (GP%) of 395 helical QA of patient-specific tomotherapy, acquired with ArcCHECK, is presented, analysed and correlated to various parameters of the plan. Following TG-218 recommendations, the clinic specific action limit (AL_cs) and tolerance limit (TL_cs) were calculated for our clinic and monitored during the analysed period.

Results

The mean values ​​(± 1 standard deviation) of GP% (3%/2 mm) (both global and local normalization) are: 97.6% and 90.9%, respectively. The proposed AL_cs and TL_cs, after a period of two years’ process monitoring are 89.4% and 91.1% respectively.

Conclusions

The phantom measurements closely match the planned dose distributions, demonstrating that the calculation accuracy of the new Precision TPS and the delivery accuracy of the Radixact unit are adequate, with respect to international guidelines and reports. Furthermore, a first correlation with the planning parameters was made. Action and tolerance limits have been set for the new Radixact Linac.

Fluence-weighted average subfield size in helical TomoTherapy

INTRODUCTION

Helical TomoTherapy allows a highly conformal dose distribution to complex target geometries with a good protection of organs at risk. However, the small field sizes associated with this method are a possible source of dosimetrical uncertainties. The IAEA has published detector-specific field output correction factors for static fields of the TomoTherapy in the TRS483. This work investigates the average subfield size of helical TomoTherapy plans.

MATERIAL AND METHODS

A new parameter for helical TomoTherapy was defined - the fluence-weighted average subfield size. The subfield sizes were extracted from the leaf-opening time sinograms in the RT-plan files for 30 clinical prostate and head and neck plans and were put in relation to Delat4 Phantom+ measurement results. Additionally the influence of planning parameters on the subfield size was studied by varying the modulation factor, number of iterations and pitch in the dose optimization and calculation for three different clinical indications H&N, prostate and rectum cancer. Selected plans were dosimetrically verified by Delta4 measurements to examine the reliability in dependence of the average subfield size. Furthermore, the impact of the planning parameters on a) the dose distribution, with regard to the planning target volume and regions at risks, and b) machine characteristics such as delivery time, actual modulation factor and leaf-opening times were evaluated.

RESULTS

The average equivalent square subfield lengths (s¯_eq) of the two investigated indications did not differ significantly - prostate plans: 2.75±0.14cm and H&N plans: 2.70±0.16cm, both with a jaw width of 2.5cm. No correlation between field size and measured dose deviation was detected. The number of iterations and the modulation factor have a considerable influence on the average subfield size. The higher the planned modulation factor and the more iterations are used during optimization, the smaller is the subfield size. In our pilot study plans were calculated with field sizes s¯_eq between 4.2cm and 1.7cm, with a jaw width of 2.5cm. Again, the measurement results of Delta4 showed no significant deviation from the doses calculated by the TomoTherapy planning system for the whole range of subfield sizes, and no significant correlation between field sizes and dose deviations was found. As expected, the clinical dose distribution improved with increasing modulation factor and an increasing number of iterations. The compromise between an improved dose distribution and smaller s¯_eq was shown.

CONCLUSION

In this work, a method was presented to determine the average subfield size for helical TomoTherapy plans. The response of the Delta4 did not show any dependence on field size in the range of the field sizes covered by the studied plans. The influence of the subfield sizes on other dosimetry systems remains to be investigated.

Density scaling of phantom materials for a 3D dose verification system

In this study, the optimum density scaling factors of phantom materials for a commercially available three‐dimensional (3D) dose verification system (Delta4) were investigated in order to improve the accuracy of the calculated dose distributions in the phantom materials. At field sizes of 10 × 10 and 5 × 5 cm² with the same geometry, tissue‐phantom ratios (TPRs) in water, polymethyl methacrylate (PMMA), and Plastic Water Diagnostic Therapy (PWDT) were measured, and TPRs in various density scaling factors of water were calculated by Monte Carlo simulation, Adaptive Convolve (AdC, Pinnacle³), Collapsed Cone Convolution (CCC, RayStation), and AcurosXB (AXB, Eclipse). Effective linear attenuation coefficients (μ_eff) were obtained from the TPRs. The ratios of μ_eff in phantom and water ((μ_eff)_pl,water) were compared between the measurements and calculations. For each phantom material, the density scaling factor proposed in this study (DSF) was set to be the value providing a match between the calculated and measured (μ_eff)_pl,water. The optimum density scaling factor was verified through the comparison of the dose distributions measured by Delta4 and calculated with three different density scaling factors: the nominal physical density (PD), nominal relative electron density (ED), and DSF. Three plans were used for the verifications: a static field of 10 × 10 cm² and two intensity modulated radiation therapy (IMRT) treatment plans. DSF were determined to be 1.13 for PMMA and 0.98 for PWDT. DSF for PMMA showed good agreement for AdC and CCC with 6 MV x ray, and AdC for 10 MV x ray. DSF for PWDT showed good agreement regardless of the dose calculation algorithms and x‐ray energy. DSF can be considered one of the references for the density scaling factor of Delta4 phantom materials and may help improve the accuracy of the IMRT dose verification using Delta4.

Clinical implementation of an exit detector-based dose reconstruction tool for helical tomotherapy delivery quality assurance

PURPOSE

The aim of this study was to validate the accuracy of an exit detector-based dose reconstruction tool for helical tomotherapy (HT) delivery quality assurance (DQA).

METHODS AND MATERIAL

Exit detector-based DQA tool was developed for patient-specific HT treatment verification. The tool performs a dose reconstruction on the planning image using the sinogram measured by the HT exit detector with no objects in the beam (i.e., static couch), and compares the reconstructed dose to the planned dose. Vendor supplied (three "TomoPhant") plans with a cylindrical solid water ("cheese") phantom were used for validation. Each "TomoPhant" plan was modified with intentional multileaf collimator leaf open time (MLC LOT) errors to assess the sensitivity and robustness of this tool. Four scenarios were tested; leaf 32 was "stuck open," leaf 42 was "stuck open," random leaf LOT was closed first by mean values of 2% and then 4%. A static couch DQA procedure was then run five times (once with the unmodified sinogram and four times with modified sinograms) for each of the three "TomoPhant" treatment plans. First, the original optimized delivery plan was compared with the original machine agnostic delivery plan, then the original optimized plans with a known modification applied (intentional MLC LOT error) were compared to the corresponding error plan exit detector measurements. An absolute dose comparison between calculated and ion chamber (A1SL, Standard Imaging, Inc., WI, USA) measured dose was performed for the unmodified "TomoPhant" plans. A 3D gamma evaluation (2%/2 mm global) was performed by comparing the planned dose ("original planned dose" for unmodified plans and "adjusted planned dose" for each intentional error) to exit detector-reconstructed dose for all three "Tomophant" plans. Finally, DQA for 119 clinical (treatment length <25 cm) and three cranio-spinal irradiation (CSI) plans were measured with both the ArcCHECK phantom (Sun Nuclear Corp., Melbourne, FL, USA) and the exit detector DQA tool to assess the time required for DQA and similarity between two methods.

RESULTS

The measured ion chamber dose agreed to within 1.5% of the reconstructed dose computed by the exit detector DQA tool on a cheese phantom for all unmodified "Tomophant" plans. Excellent agreement in gamma pass rate (>95%) was observed between the planned and reconstructed dose for all "Tomophant" plans considered using the tool. The gamma pass rate from 119 clinical plan DQA measurements was 94.9% ± 1.5% and 91.9% ± 4.37% for the exit detector DQA tool and ArcCHECK phantom measurements (P = 0.81), respectively. For the clinical plans (treatment length <25 cm), the average time required to perform DQA was 24.7 ± 3.5 and 39.5 ± 4.5 min using the exit detector QA tool and ArcCHECK phantom, respectively, whereas the average time required for the 3 CSI treatments was 35 ± 3.5 and 90 ± 5.2 min, respectively.

CONCLUSION

The exit detector tool has been demonstrated to be faster for performing the DQA with equivalent sensitivity for detecting MLC LOT errors relative to a conventional phantom-based QA method. In addition, comprehensive MLC performance evaluation and features of reconstructed dose provide additional insight into understanding DQA failures and the clinical relevance of DQA results.

Long-term dosimetric stability of multiple TomoTherapy delivery systems

The dosimetric stability of six TomoTherapy units was analyzed to investigate changes in performance over time and with system upgrades. Energy and output were tracked using monitor chamber signal, onboard megavoltage computed tomography (MVCT) detector profile, and external ion chamber measurements. The systems (and monitoring periods) include three Hi-Art (67, 61, and 65 mos.), two TomoHDA (31 and 26 mos.), and one Radixact unit (11 mos.), representing approximately 10 years of clinical use. The four newest systems use the Dose Control Stability (DCS) system and Fixed Target Linear Accelerator (linac) (FTL). The output stability is reported as deviation from reference monitor chamber signal for all systems and/or from an external chamber signal. The energy stability was monitored using relative (center versus off-axis) MVCT detector signal (beam profile) and/or the ratio of chamber measurements at 2 depths. The clinical TomoHDA data were used to benchmark the Radixact stability, which has the same FTL but runs at a higher dose rate. The output based on monitor chamber data of all systems is very stable. The standard deviation of daily output on the non-DCS systems was 0.94-1.52%. As expected, the DCS systems had improved standard deviation: 0.004-0.06%. The beam energy was also very stable for all units. The standard deviation in profile flatness was 0.23-0.62% for rotating target systems and 0.04-0.09% for FTL. Ion chamber output and PDD ratios supported these results. The output stability on the Radixact system during extended treatment delivery (20, 30, and 40 min) was comparable to a clinical TomoHDA system. For each system, results are consistent between different measurement tools and techniques, proving not only the dosimetric stability, but also these quality parameters can be confirmed with various metrics. The replacement history over extended time periods of the major dosimetric components of the different delivery systems (target, linac, and magnetron) is also reported.

Measurement-guided volumetric dose reconstruction for helical tomotherapy

It was previously demonstrated that dose delivered by a conventional linear accelerator using IMRT or VMAT can be reconstructed - on patient or phantom datasets - using helical diode array measurements and a technique called planned dose perturbation (PDP). This allows meaningful and intuitive analysis of the agreement between the planned and delivered dose, including direct comparison of the dose-volume histograms. While conceptually similar to modulated arc techniques, helical tomotherapy introduces significant challenges to the PDP formalism, arising primarily from TomoTherapy delivery dynamics. The temporal characteristics of the delivery are of the same order or shorter than the dosimeter's update interval (50 ms). Additionally, the prevalence of often small and complex segments, particularly with the 1 cm Y jaw setting, lead to challenges related to detector spacing. Here, we present and test a novel method of tomotherapy-PDP (TPDP) designed to meet these challenges. One of the novel techniques introduced for TPDP is organization of the subbeams into larger subunits called sectors, which assures more robust synchronization of the measurement and delivery dynamics. Another important change is the optional application of a correction based on ion chamber (IC) measurements in the phantom. The TPDP method was validated by direct comparisons to the IC and an independent, biplanar diode array dosimeter previously evaluated for tomotherapy delivery quality assurance. Nineteen plans with varying complexity were analyzed for the 2.5 cm tomotherapy jaw setting and 18 for the 1 cm opening. The dose differences between the TPDP and IC were 1.0% ± 1.1% and 1.1% ± 1.1%, for 2.5 and 1.0 cm jaw plans, respectively. Gamma analysis agreement rates between TPDP and the independent array were: 99.1%± 1.8% (using 3% global normalization/3 mm criteria) and 93.4% ± 7.1% (using 2% global/2 mm) for the 2.5 cm jaw plans; for 1 cm plans, they were 95.2% ± 6.7% (3% G/3) and 83.8% ± 12% (2% G/2). We conclude that TPDP is capable of volumetric dose reconstruction with acceptable accuracy. However, the challenges of fast tomotherapy delivery dynamics make TPDP less precise than the IMRT/VMAT PDP version, particularly for the 1 cm jaw setting.

Independent calculation of dose distributions for helical tomotherapy using a conventional treatment planning system

PURPOSE

The dosimetric verification of treatment plans in helical tomotherapy usually is carried out via verification measurements. In this study, a method for independent dose calculation of tomotherapy treatment plans is presented, that uses a conventional treatment planning system with a pencil kernel dose calculation algorithm for generation of verification dose distributions based on patient CT data.

METHODS

A pencil beam algorithm that directly uses measured beam data was configured for dose calculation for a tomotherapy machine. Tomotherapy treatment plans were converted into a format readable by an in-house treatment planning system by assigning each projection to one static treatment field and shifting the calculation isocenter for each field in order to account for the couch movement. The modulation of the fluence for each projection is read out of the delivery sinogram, and with the kernel-based dose calculation, this information can directly be used for dose calculation without the need for decomposition of the sinogram. The sinogram values are only corrected for leaf output and leaf latency. Using the converted treatment plans, dose was recalculated with the independent treatment planning system. Multiple treatment plans ranging from simple static fields to real patient treatment plans were calculated using the new approach and either compared to actual measurements or the 3D dose distribution calculated by the tomotherapy treatment planning system. In addition, dose-volume histograms were calculated for the patient plans.

RESULTS

Except for minor deviations at the maximum field size, the pencil beam dose calculation for static beams agreed with measurements in a water tank within 2%/2 mm. A mean deviation to point dose measurements in the cheese phantom of 0.89% ± 0.81% was found for unmodulated helical plans. A mean voxel-based deviation of -0.67% ± 1.11% for all voxels in the respective high dose region (dose values >80%), and a mean local voxel-based deviation of -2.41% ± 0.75% for all voxels with dose values >20% were found for 11 modulated plans in the cheese phantom. Averaged over nine patient plans, the deviations amounted to -0.14% ± 1.97% (voxels >80%) and -0.95% ± 2.27% (>20%, local deviations). For a lung case, mean voxel-based deviations of more than 4% were found, while for all other patient plans, all mean voxel-based deviations were within ± 2.4%.

CONCLUSIONS

The presented method is suitable for independent dose calculation for helical tomotherapy within the known limitations of the pencil beam algorithm. It can serve as verification of the primary dose calculation and thereby reduce the need for time-consuming measurements. By using the patient anatomy and generating full 3D dose data, and combined with measurements of additional machine parameters, it can substantially contribute to overall patient safety.

Cross-validation of two commercial methods for volumetric high-resolution dose reconstruction on a phantom for non-coplanar VMAT beams

BACKGROUND AND PURPOSE

Delta(4) (ScandiDos AB, Uppsala, Sweden) and ArcCHECK with 3DVH software (Sun Nuclear Corp., Melbourne, FL, USA) are commercial quasi-three-dimensional diode dosimetry arrays capable of volumetric measurement-guided dose reconstruction. A method to reconstruct dose for non-coplanar VMAT beams with 3DVH is described. The Delta(4) 3D dose reconstruction on its own phantom for VMAT delivery has not been thoroughly evaluated previously, and we do so by comparison with 3DVH.

MATERIALS AND METHODS

Reconstructed volumetric doses for VMAT plans delivered with different table angles were compared between the Delta(4) and 3DVH using gamma analysis.

RESULTS

The average γ (2% local dose-error normalization/2mm) passing rate comparing the directly measured Delta(4) diode dose with 3DVH was 98.2 ± 1.6% (1SD). The average passing rate for the full volumetric comparison of the reconstructed doses on a homogeneous cylindrical phantom was 95.6 ± 1.5%. No dependence on the table angle was observed.

CONCLUSIONS

Modified 3DVH algorithm is capable of 3D VMAT dose reconstruction on an arbitrary volume for the full range of table angles. Our comparison results between different dosimeters make a compelling case for the use of electronic arrays with high-resolution 3D dose reconstruction as primary means of evaluating spatial dose distributions during IMRT/VMAT verification.

Radical pleurectomy/decortication followed by high dose of radiation therapy for malignant pleural mesothelioma. Final results with long-term follow-up

PURPOSE

We have previously shown the feasibility of delivering high doses of radiotherapy in malignant pleural mesothelioma (MPM) patients who underwent radical pleurectomy/decortication (P/D) or surgical biopsy. In this report, we present the long-term results of MPM patients treated with radical P/D followed by high doses of radiotherapy.

METHODS AND MATERIALS

Twenty consecutive MPM patients were enrolled in this prospective study and underwent radical P/D followed by high dose radiotherapy. The clinical target volume was defined as the entire hemithorax excluding the intact lung. The dose prescribed was 50 Gy in 25 fractions. Any FDG-avid areas or regions of particular concern for residual disease were given a simultaneous boost to 60 Gy. Nineteen patients received cisplatin/pemetrexed chemotherapy. Kaplan-Meier analysis was used to calculate rates of overall survival (OS), progression-free survival (PFS), and loco-regional control (LRC).

RESULTS

The median follow-up was of 27 months. The median OS and PFS were 33 and 29 months, respectively. The median LRC was not reached. The Kaplan-Meier estimates of OS at 2 and 3 years were 70% and 49%, respectively. The estimates of PFS at 2 and 3 years were 65% and 46%, respectively. The estimates of LRC at 2 and 3 years were 68% and 59%, respectively. The predominant pattern of failure was distant: 7 patients developed distant metastases as the first site of relapse, whereas only 3 patients experienced an isolated loco-regional recurrence. No fatal toxicity was reported. Five Grades 2-3 pneumonitis were documented.

CONCLUSIONS

High dose radiation therapy following radical P/D led to excellent loco-regional control and survival results in MPM patients. A median OS of 33 months and a 3-year OS rate of 49% are among the best observed in recent studies, supporting the idea that this approach represents a concrete therapeutic option for malignant pleural mesothelioma.

[A practical procedure to improve the accuracy of radiochromic film dosimetry: a integration with a correction method of uniformity correction and a red/blue correction method]

It has been reported that the light scattering could worsen the accuracy of dose distribution measurement using a radiochromic film. The purpose of this study was to investigate the accuracy of two different films, EDR2 and EBT2, as film dosimetry tools. The effectiveness of a correction method for the non-uniformity caused from EBT2 film and the light scattering was also evaluated. In addition the efficacy of this correction method integrated with the red/blue correction method was assessed. EDR2 and EBT2 films were read using a flatbed charge-coupled device scanner (EPSON 10000G). Dose differences on the axis perpendicular to the scanner lamp movement axis were within 1% with EDR2, but exceeded 3% (Maximum: +8%) with EBT2. The non-uniformity correction method, after a single film exposure, was applied to the readout of the films. A corrected dose distribution data was subsequently created. The correction method showed more than 10%-better pass ratios in dose difference evaluation than when the correction method was not applied. The red/blue correction method resulted in 5%-improvement compared with the standard procedure that employed red color only. The correction method with EBT2 proved to be able to rapidly correct non-uniformity, and has potential for routine clinical IMRT dose verification if the accuracy of EBT2 is required to be similar to that of EDR2. The use of red/blue correction method may improve the accuracy, but we recommend we should use the red/blue correction method carefully and understand the characteristics of EBT2 for red color only and the red/blue correction method.

Quantitative analysis of geometric information from an end-to-end examination of IMRT and VMAT using the optimal selection method

PURPOSE

Gamma index, distance-to-agreement, and dose difference (DD) are commonly used to evaluate planar dose distributions. In this evaluation, the agreement between calculated and measured dose distributions can be susceptible to steep dose gradients along another axis perpendicular to the evaluation plane. Visual registration of the measured dose distribution may be performed to achieve better agreement, although doing so might lose geometric information related to beam targeting in an end-to-end test of intensity modulated radiation therapy (IMRT) and volumetric modulated arc therapy (VMAT). The optimal selection (OS) method was developed to take into consideration a dose distribution in three-dimensions, and also to quantitatively analyze geometric information along with better agreement.

METHODS

The OS method was composed of two steps. These steps were based on two algorithms, the gamma index and DD, to (1) find the best-matched plane, which is parallel to the planar measured dose distribution and is reconstructed by a volumetric dose distribution calculated by a treatment planning system; and (2) to get shifts and rotation along with better agreement between the calculated and measured dose distribution, compared with the planar dose distribution from the test. The OS method computes shifts and rotation against a user-defined coregistered location for the measured dose distribution. Thirteen prostate IMRT plans (two planes per plan for a total of 26 planes) were analyzed retrospectively to compare the pass ratios of DD and gamma index evaluations with and without the OS method. The computed shifts and rotations were evaluated.

RESULTS

Compared with the method without OS, the average pass ratios of DD and gamma index with the OS method increased by 8.2% and 5.7%, respectively, in the dose region from 30% to 100%. A particular result from one of the planes showed an increase of 43.5% and 32.5% in the pass ratios of DD and gamma, respectively, with the OS method in the same dose region. The shifts in the x-, y-, z-axes and rotation, which were computed using the OS method, were 0.5 ± 0.6, 0.3 ± 0.5, 1.0 ± 1.1 mm, and 0.3 ± 0.3°, respectively. In terms of the comparatively large difference between the z-shift and the x- and y-shifts, an additional geometric test was performed. A systematic error of 0.7 mm in the z-axis was found at the location of the film placed in the phantom that we used.

CONCLUSIONS

The OS method improved the quality of the end-to-end test of IMRT and VMAT plans by providing additional information regarding shifts and rotation, which were calculated and found to be in better agreement.

Evaluation of PTW Seven29 for tomotherapy patient-specific quality assurance and comparison with ScandiDos Delta(4)

For routine quality assurance of helical tomotherapy plans, an alternative method, as opposed to the TomoTherapy suggested cylindrical solid water phantom with film and ionization chamber, is proposed using the PTW Seven29 2D-ARRAY inserted in a dedicated octagonal phantom, called Octavius. First, the sensitivity of the array to pitch was studied by varying the pitch during planning to 0.287, 0.433, 1.0, and 2.0. For each pitch selected, the dependence on field size was investigated by generating plans with field widths (FWs) of 1.06 cm, 2.49 cm, and 5.02 cm, for a total of 12 plans. Secondly, a total of 15 patient QA plans were delivered using helical tomotherapy with the Delta(4) and Seven29/Octavius for comparison. Using the clinical gamma criteria, 3% and 3 mm, all FW and pitch plans had a passing percentage of >90%. For patient QA plans, the average gamma pass percentage was 97.0% (94.4-99.8%) for the Delta(4) and 97.6% (92.5-100.0%) for the Seven29/Octavius. Both the Seven29/Octavius and Delta(4) performed to a high standard of measurement accuracy and had a 90% or greater gamma percent for all plans and were considered clinically acceptable.