Treating lung cancer with dynamic conformal arc therapy: a dosimetric study

Background Lung cancer patients are often in poor physical condition, and a shorter treatment time would reduce their discomfort. Dynamic conformal arc therapy (DCAT) offers a shorter treatment time than conventional 3D conformal radiotherapy (3D CRT) and is usually available even in departments without inverse planning possibilities. We examined its suitability as a treatment modality for lung cancer patients. Methods On a cohort of 35 lung cancer patients, relevant dosimetric parameters were compared in respective DCAT and 3D CRT treatment plans. Radiochromic film dosimetry in an anthropomorphic phantom was used to compare both DCAT and 3D CRT dose distributions against their planned counterparts. Results In comparison with their 3D CRT counterparts, DCAT plans equal or exceed the agreement between the calculated dose and the dose measured using film dosimetry. In dosimetric comparison, DCAT performed significantly better than 3D CRT in dose conformity to PTV and the number of monitor units used per plan, and significantly worse in dose homogeneity, mean lung dose and lung volume exposed to 5 Gy or more (V5Gy). No significant difference was found in the V20Gy value to lung, dose to 1 cm3 of spinal cord, and the mean dose to oesophagus. Improvements in V20Gy and V5Gy were found to be negatively correlated. DCAT plans differ from 3D CRT by exhibiting a moderate negative correlation between target volume sphericity and dose homogeneity. Conclusions With respect to the agreement between the planned and the irradiated dose distribution, DCAT appears at least as reliable as 3D CRT. In specific conditions concerning the patient anatomy and treatment prescription, DCAT may yield more favourable dosimetric parameters. On average, however, conventional 3D CRT usually obtains better dosimetric parameters. We can thus only recommend DCAT as a complementary technique to the conventional 3D CRT. Electronic supplementary material The online version of this article (doi:10.1186/s13014-017-0823-y) contains supplementary material, which is available to authorized users.


Background
Lung cancer is the most commonly diagnosed cancer in males worldwide, accounting for 17% of the total new cancer cases and 24% of total cancer-related deaths [1]. In Slovenia, 1214 new lung cancer patients were diagnosed in 2012, 906 men and 308 women, which comprises 12% and 6% of all new cancer cases, respectively [2].
Radiotherapy is an important modality in the treatment of lung cancer either in curative or in palliative setting. In non-small cell lung cancer, curative radiotherapy includes patients with inoperable early stage and radiotherapy combined with concomitant chemotherapy *Correspondence: ppeterlin@onko-i.si Institute of Oncology Ljubljana, Zaloška 2, SI-1000, Ljubljana, Slovenia for patients with locally advanced stage. Curative radiotherapy in combination with chemotherapy is also indicated for limited stage small cell lung cancer. Among new patients diagnosed with lung cancer in Slovenia in 2012, 43.5% received radiotherapy as part of their first line treatment [2].
Since its inception almost two decades ago, intensitymodulated radiation therapy (IMRT) has proved to be a great improvement to radiotherapy, allowing for better sparing of healthy organs and a higher dose to the tumour. However, as the individual treatment field segments do not encompass the whole treatment volume projection, a concern has arisen that the interplay between the movement of multileaf collimator (MLC) leaves and the organ movement (e.g., due to breathing motion) may lead to either overdosage or underdosage of the treated volume and/or the healthy tissue. The concern about the interplay between the MLC and organ motion applies to every form of radiotherapy in which only part of the treatment volume is irradiated at a time. This applies to all forms of intensitymodulated therapy, including advanced techniques such as volumetric modulated arc therapy (VMAT) and tomotherapy, the "field-in-field" technique, or "forwardplanning IMRT" technique, as well as dynamic-or virtual wedges.
A statistical analysis of a simulated patient motion [3] has shown that due to organ motion, dose in any given voxel delivered in a single fraction deviates from its expected value, with its standard deviation reaching up to 10% of the expected dose value. By splitting the dose into several fractions, however, the probability density function of dose deviations quickly converges towards a Gaussian, and for a conventional fractionation (e.g., 30 fractions), the standard deviation for dose deviation drops down to 1% for a 5 mm organ motion amplitude. An experimental verification with a motor-driven phantom to simulate organ motion [4] was devised shortly afterwards. Using an ionization chamber, the authors confirmed their previous findings, and in addition found out that lowering the dose rate (i.e., 300 vs. 500 MU/min) reduces dose variation due to organ motion. Another experimental investigation which used film dosimetry has shown that underdosage as large as 6% may occur [5]. Early research in this area led to the recommendations on management of respiratory motion [6].
Recently, there has been a renewed interest in treating lung cancer with dynamic conformal arc therapy [7,8]. Conformal arc therapy avoids the interplay effect altogether, since the field shape encompasses the complete projection of the treated volume from every direction. This, along with the shorter treatment time, is considered the advantage of DCAT over 3D CRT treatment plans employing "field-in-field" technique.
The purpose of the present study is twofold. First, we want to ascertain that the dose distribution in lungs delivered by dynamic conformal arc therapy (DCAT) exhibits at least as good agreement with the planned dose distribution as the conventional 3D conformal radiotherapy (3D CRT). Second, we want to examine and compare the relevant radiotherapy parameters in both DCAT and 3D CRT treatment plans.

Patient selection
Based exclusively on dosimetric criteria, 35 patients were included in the study. They were considered suitable for DCAT treatment if either the total dose was below the value considered to cause myelopathy of the spinal cord, or if the distance between the planned treated volume (PTV) and the spinal cord was estimated as sufficient so that the beams traversing the spinal cord would not add up to a dose contribution high enough to cause myelopathy.
Target volume, critical organs, and dose prescription PTV was created by applying a 5-10 mm isotropic margin to clinical target volume (CTV), depending on the location of the tumour (apex vs. close to the diaphragm). Figure 1 shows the metrics characterizing the target size, its location and their relations with the prescribed dose and the clinical stage. PTV volumes range from 145.2 to 2686.5 cm 3 , with 70% of all cases corresponding to the 400-1000 cm 3 range (Fig. 1a).
Lungs were contoured semi-automatically, based on Hounsfield number threshold, using the algorithm provided by the Eclipse v10.0 treatment planning system (Varian Medical Systems, Palo Alto, CA, USA). Both left and right lung were contoured as a single lung structure. The target was not excluded from the lungs; instead, there was an overlap between the lung structure and the target. On average, 68.9 ± 45.5 cm 3 (27.3 ± 20.3%) of CTV and 214.4 ± 85.3 cm 3 (38.0 ± 19.6%) of PTV overlapped with the lung structure.
The correlation between the PTV volume and the dose prescribed to PTV (Fig. 1b) is not considered significant (Spearman's rank correlation coefficient ρ = −0.229, p = 0.185), although, as expected, the few cases of radical treatment doses included in the study were limited to PTV volumes below 1000 cm 3 . The standard palliative dose prescription for lung patients at our clinic is 36 Gy in 12 fractions; these cases comprise the bulk of the studied cases (Table 1).The correlation between the dose prescribed to PTV and the minimal distance between the PTV and the spinal cord ( Fig. 1c) is also not significant (ρ = 0.254, p = 0.140). On the other hand, significant correlation exists between the PTV volume and the reciprocal value of the minimal distance between the PTV and the spinal cord (ρ = 0.525, p = 0.001; Fig. 1d). Boxplots show that the PTV volume generally increases with the clinical stage (Fig. 1e), while the minimal distance between the CTV and the spinal cord decreases with the clinical stage (Fig. 1f). PTV volume was computed using the treatment planning system (see below), while the minimal distance between either the PTV or the CTV and the spinal cord was computed using an in-house script written in Python which scanned the slices of DICOM RT Structure file for every patient.

Treatment planning
Eclipse v10.0 was used as treatment planning system (TPS), using the AAA photon dose calculation algorithm.  fields. In 20 cases out of 35, one or more supplementary "field-in-field" fields (i.e., "forward-planning IMRT") were used in addition to the principal treatment fields. Treatment fields were shaped by applying an 8 mm margin to the planning target volume (PTV) in the cranio-caudal direction and a 5 mm margin to the PTV in the transversal directions. Treatment planning restrictions used in our department generally follow QUANTEC [9]. For lungs the prescriptions are V20Gy < 35%, V5Gy < 60%, and MLD < 15 Gy, for spinal cord, D max < 50 Gy (prevails over target coverage), and for oesophagus, D < 34 Gy. All DCAT plans were done identically, using a single arc with gantry running in a clockwise direction from 182°to 178°, with the control points evenly spaced every 4°, and the collimator tilted to 45°. An isotropic 5 mm margin was applied to PTV. In order to evaluate the dosimetric impact of the collimator angle, three additional settings for collimator angle were also tested: 0°, 30°, and 90°, yielding an overall number of 4×35 DCAT plans.

Dose verification with film dosimetry
8 in × 10 in Gafchromic EBT3 films (Ashland, Wayne, NJ, USA) from lots A03181301 and 04071601 were used for dose verification. The films were scanned before and 24 ± 1 h after the irradiation using Epson Expression 10000XL flatbed scanner (Seiko Epson Corp., Nagano, Japan) in transmission mode, driven by Epson Scan v3.0 software. Images were acquired in 48-bit RGB mode using 72 dpi resolution with the image-processing options turned off. The doses were obtained using the Radiochromic.com (http://radiochromic.com/) web application [10], v1.1 through v2.5.
For dose verification purposes, 10 patients were randomly selected from the group. Both their 3D conformal plans and their DCAT plans were recalculated on the CIRS Thorax phantom (model 002LFC; Computerized Imaging Reference Systems, Norfolk, VA, USA), and the calculated dose plane taken 26 mm caudally from the isocenter coinciding with the phantom center. The chosen plane corresponding to the gap between plates marked #2 and #3 on the phantom was exported in DICOM RT Dose format. Again, a sheet of film was scanned first prior to being irradiated inside the CIRS Thorax phantom, and then again 24 ± 1 h after the irradiation (Fig. 2, top row). The scanned images of the film before and after the irradiation, the exported dose plane (Fig. 2, middle row) and the daily output factor were imported into the Radiochromic.com web application to obtain the dose distribution and to conduct gamma index analysis.

Dose calculation in different media
For each patient and for each treatment modality (3D CRT, DCAT), we delineated the regions with different electron densities (lung, soft tissue, bone) on the registered images of the calculated and the measured planar dose distributions. For every pixel, we calculated the where D eval is the evaluation dose (calculated by the TPS) and D ref is the reference dose (measured using the radiochromic film). Next we averaged the dose over the region for a given treatment plan, and finally we computed the systematic error M, its standard deviation , and the random error σ [11].

Data analysis
Gamma index analysis. Gamma index analysis [12] was used for comparing the reference and the evaluation dose distributions. For every point r referring to the reference dose distribution, one can define a function at point r, D and d are the dose deviation and DTA criteria (commonly taken as 3% of the maximal dose and 3 mm), and r is a point referring to the evaluation dose distribution. Common criteria for agreement between the two distributions are the ratio of points r for which γ (r) < 1 holds, and the average valueγ .
Dosimetric comparison of treatment plans. Dosevolume histograms were exported from the treatment planning system and subsequently analysed with an inhouse script written in GNU R (http://www.R-project. org/; [13]) using RadOnc package [14].

Dose conformity. A measure of dose conformity is the conformation number defined as CN
where TV RI is the target volume covered by the reference isodose, TV is the target volume, and V RI is the volume covered by the reference isodose [15]. In this case, the reference isodose was set to 95% of the prescription isodose. A higher CN value signifies a better conformity of the therapeutic dose to the target volume.

Dose homogeneity.
A measure of dose homogeneity is the homogeneity index [16] defined as HI = (D 2 − D 98 )/D nom , where D 2 is the dose to the 2% of the target volume (i.e., 2% of the target volume receives this dose or higher), D 98 is the dose to the 98% of the target volume (i.e., 98% of the target volume receives this dose or higher), and D nom is the nominal prescribed dose. A lower HI value indicates a more homogeneous target dose. Another measure of PTV homogeneity, dose homogeneity dispersion (DHD), is defined as the standard deviation of the absorbed dose covering PTV, divided by its mean value. A lower value of this parameter indicates a more homogeneous target dose. A similar parameter is the S-index proposed by Yoon [17].
Target volume sphericity. DCAT performs worse on highly concave target shapes with many protrusions. In order to quantify target volume shape, we have introduced the target volume sphericity parameter (Additional file 1). Sphericity is a dimensionless parameter which attains values close to 1 for near-spherical shapes, and falls towards 0 as the shape departs from the sphere.

Target volume location
Our initial hypothesis was that DCAT works best with centrally located tumours, because treating peripheral targets would lead to hot spots in the areas proximal to the skin and cold spots in the distal areas, thus degrading dose homogeneity. In order to quantify tumour location, we introduced two geometrical parameters. The first one is the magnitude of the treatment field isocentre displacement from the patient origin (reference isocentre) in the transversal plane. Assuming that the patient origin is generally selected close to the centre of mass of an average cross-section, the displacement, calculated as a square root of the squares of displacements in the mediolateral (x ml ) and the antero-posterior direction (x ap ), R = x ml + x ap , is a measure of how centrally located is the target. The second parameter is the minimal distance between the CTV and the external contour (i.e., skin), which was computed by an in-house script written in Python which scanned the slices of DICOM RT Structure file for every patient.

Dose verification with film dosimetry
Gamma index analysis [12] was used to compare the planar dose distributions for the 10 randomly selected pairs (conventional 3D CRT and DCAT) of treatment plans.
The dose obtained from radiochromic film was taken as the reference distribution, and the corresponding dose plane exported from TPS was taken as the evaluation distribution. Several different values for the dose tolerance, positional tolerance, and dose threshold were used. Normalization to global D max was used in all cases.
With the dose threshold set to 10% D max and using 3% as dose tolerance and 3 mm as positional tolerance, the median value of the points with γ < 1 was 94.8% (lower quartile 91.5%, upper quartile 97.88%) for conventional 3D CRT plans, and 97.3% (96.5%, 98.6%) for DCAT plans. The corresponding values ofγ computed with the same parameters yielded a median of 0.38 (0.31, 0.42) for conventional 3D CRT plans and 0.26 (0.25, 0.30) for DCAT plans. Setting dose and positional tolerances to (2%, 2 mm), we obtain for γ < 1 a median of 88.2% (85.1%, 91.3%) for conventional 3D CRT plans and 95.5% (94.5%, 98.6%) for DCAT plans. The corresponding values ofγ (2%, 2 mm are 0.49 (0.38, 0.55) for conventional 3D CRT and 0.33 (0.30, 0.38) for DCAT. Repeating the calculation with the dose threshold set to 60% D max yielded results which are in agreement with the ones presented here, i.e., DCAT consistently resulted in slightly better results compared to conventional 3D CRT. Figure 3 shows a comparison of the gamma index analysis for assessing the agreement between the measured and the calculated dose distribution. Using global D max , 10% dose threshold, 1 mm positional tolerance, and 2, 3, 4 and 5% dose tolerance, the ratio of points passing the γ < 1 criterion was computed for conventional 3D CRT treatment plans (a) and DCAT treatment plans (b). We can see that overall, under the same conditions, DCAT treatment plans yield a higher ratio of points passing the γ < 1 criterion. Both techniques reveal a single outlier (the same case with both conventional 3D CRT treatment plans and DCAT treatment plans) which exhibits a markedly worse results than the rest. With the exception of the said outlier, with DCAT, most plans pass the γ < 1 criterion already with 1 mm positional and 3% dose tolerance.
The gamma value distribution maps (Fig. 2e, f) show that both in conventional 3D CRT plans and in DCAT plans, the areas with a high gamma index are predominantly located along the phantom boundary where the treatment beams enter the phantom. These areas only reflect the inaccuracy of registration and/or the inacurracy of dose calculation in the buildup region. In conventional 3D CRT plans (Fig. 2e), another area with high gamma index is along the lines which correspond to the boundaries of treatment fields, which can be attributed both to the error in phantom positioning as well as to a disagreement between the calculated and the measured dose in the penumbra region. In DCAT plans, dose gradients associated with the boundaries of treatment fields cannot be discerned. However, since most conventional 3D CRT treatments are fractionated, these dose gradients would be smeared out to some degree, which would likely diminish the difference in dose agreement between the two techniques even further.

Dose calculation in different media
We analysed the dose difference distribution for the treatment plans which were previously selected for verification with the CIRS Thorax Phantom. The results for dose calculation in lung, bone and soft tissue are shown in Table 2. There are several results discernible from this table. First, the dose difference in the bone exhibits a systematic negative shift for both techniques (M = 4.2% for conventional 3D CRT, 2.4% for DCAT), indicating that the treatment planning system underestimates the dose in bone regardless of the treatment modality and the individual patient geometry. Dose difference in the lung and the soft tissue does not exhibit such deviation, and their offsets are smaller than or comparable to . Only in the lung does the Wilcoxon signed rank test show that we can convincingly (p < 0.01) reject the null hypothesis that the difference between the pairs of mean dose follows a symmetric distribution centred at zero; both in bone (p = 0.56) and in soft tissue (p = 0.49), the results are inconclusive. Finally, lower values of random error σ for DCAT plans show that the 3D CRT plans generally produce a dose distribution with a larger standard deviation, which can be attributed to the regions with a high dose gradient. However, all the differences mentioned are unlikely to have a clinical significance.  x-coordinate represents the parameter value in the conventional 3D CRT plan, and its y-coordinate represents the parameter value in the DCAT plan. In addition, the two-tailed paired Student's t-test was performed in the cases in which the Shapiro-Wilk test showed normally distributed differences, and the non-parametric Wilcoxon matched pair signed rank sum test in the rest.

Dosimetric comparison of treatment plans
• Dose conformity (Fig. 4a). , oesophagus (f) and lungs, g-i) in both types of treatment plan. An individual patient case is represented as a point in the diagram; its x-coordinate represents the parameter value in the conventional 3D CRT plan, and its y-coordinate represents the parameter value in the DCAT plan. The points shown in gray correspond to cases in which neither treatment plan meets the planning restrictions difference stems from a higher TV RI /V RI ratio, or lower exposure of healthy tissue to a therapeutic dose, in the case of DCAT treatment plans. • Dose homogeneity (Fig. 4b).
The Wilcoxon signed rank test shows a statistically significant difference of -0.016 (p = 0.002, 95% CI [ −0.023, −0.0071]) in favour of 3D CRT. Dose homogeneity dispersion (DHD) is shown in Fig. 4c. The paired t-test shows a statistically significant difference of -0.0027 (p = 0.024, 95% CI [ −0.0051, −0.0004]) in favour of 3D CRT. • Monitor units (Fig. 4d) Figure 5 shows the correlation between dosimetric and geometric parameters as well as between different dosimetric parameters for lungs. Figure 5a-c show the difference between the mean dose to lung (MLD) in DCAT treatment plans and MLD in conventional 3D CRT treatment plans with respect to the volume of PTV (V PTV ), volume of lungs (V lung ) and the reciprocal value of the minimal distance between CTV and the spinal cord (1/d CTV−medulla ). Negative values of MLD DCAT − MLD 3D CRT signify an advantage of DCAT over conventional 3D CRT. Figure 5a shows a moderate and significant (Spearman's rank correlation coefficient ρ = 0.472, p = 0.005) correlation between MLD DCAT − MLD 3D CRT and V PTV , with DCAT being advantageous over conventional 3D CRT for small values of PTV volumes. Similarly, Fig. 5c shows a weak yet significant (ρ = 0.389, p = 0.022) correlation between MLD DCAT − MLD 3D CRT and 1/d CTV−medulla , with DCAT being sometimes advantageous when the minimal distance between CTV and the spinal cord is large. No significant correlation was found between MLD DCAT − MLD 3D CRT and V lung (Fig. 5b). Likewise, no significant correlation was found between either V20Gy DCAT − V20Gy 3D CRT or V5Gy DCAT − V5Gy 3D CRT on the one hand and V PTV , V lung , and 1/d CTV−medulla on the other (not shown). Figure 5d shows a moderate and significant (ρ =−0.489, p = 0.003) negative correlation between V20Gy DCAT − V20Gy 3D CRT and V5Gy DCAT − V5Gy 3D CRT . Again, negative values in Fig. 5d-f signify an advantage of DCAT over conventional 3D CRT. Dosimetric improvements of DCAT for the V20Gy value for lung seems to be incompatible with the dosimetric improvement for the V5Gy value. No significant correlation was found between V20Gy DCAT − V20Gy 3D CRT and MLD DCAT − MLD 3D CRT (Fig. 5e), while the correlation between V5Gy DCAT − V5Gy 3D CRT and MLD DCAT − MLD 3D CRT was found weak yet significant (ρ = 0.398, p = 0.018; Fig. 5f). A dosimetric improvement in V5Gy value of one technique over the other is likely to be accompanied by an improvement of MLD.

The choice of collimator angle
In this section, we examine the influence of the collimator angle on the observed dosimetric parameters. Table 3 shows a comparison for four different settings for collimator angle: 0°, 30°, 45°, and 90°. One can see that in general, the choice of collimator angle only minimally affects the value of dosimetric parameters. The largest difference between the largest and the smallest value (3.4%) occurs at V20Gy value for lung, followed by 2.5% D1cc difference for the spinal cord and 1.5% difference for mean lung dose. In all other observed parameters,  and the collimator set to 90°yielded the best value for V5Gy value for lung (24/35). Even though some choices appear to be convincingly better, we need to reiterate that the differences are very small.
Target volume sphericity Figure 6 shows the dependence of several clinical dosimetric parameters obtained by DCAT plans on the target volume sphericity. The dose homogeneity index HI exhibits a moderate and significant (ρ = −0.402, p = 0.02) negative correlation (Fig. 6a), meaning that the dose distribution is more homogenous in more spherical PTVs. The dose conformity index CN (Fig. 6b) exhibits moderate and significant (ρ = 0.463, p < 0.01) positive correlation. This is expected -when using a conformal technique to treat a highly concave target, the treated volume can be significantly larger than the target. The mean dose to oesophagus (Fig. 6c) shows a strong and significant negative correlation with (ρ = −0.676, p < 0.01). All three dosimetric parameters for lungs show a moderate and significant negative correlation with sphericity: V20Gy value ( Fig. 6d; ρ = −0.573, p < 0.01), V5Gy value ( Fig. 6e; ρ = −0.507, p < 0.01), and the mean lung dose ( Fig. 6f; ρ = −0.558, p < 0.01). On the other hand, the correlation of D1cc value for the spinal cord is weak and not statistically significant (ρ = 0.208, p = 0.23; not shown). Unlike DCAT plans, hand-crafted conventional 3D CRT plans only show a weak and not statistically significant (ρ = −0.265, p = 0.12; not shown) correlation between and dose homogeneity. In all other aspects, they behave similarly to DCAT plans: they exhibit weak and not statistically significant dependence of D1cc value for the spinal cord on , a strong and statistically significant dependence of the mean dose to oesophagus, and a moderate and statistically significant dependence of CN and all three lung parameters (not shown). Figure 7 shows the dependence of the PTV dose homogeneity (HI) and PTV dose conformity (CN) obtained by DCAT plans on either the isocentre displacement from the patient origin R or the minimal distance between the CTV and the external contour (d CTV−skin ). The correlation between homogeneity index HI and R (Fig. 7a) Fig. 7 The dependence of the PTV dose homogeneity expressed as HI (a, b) and PTV dose conformity CN (c, d) on either the isocentre displacement from the patient origin in the transversal plane and the minimal distance between CTV and external contour. All dosimetric parameters were calculated using DCAT treatment plans. An individual patient case is represented as a point in the diagram; its x-and y-coordinates represent its relevant parameter values weak and not statistically significant (ρ = 0.05, p = 0.80); the correlation between HI and d CTV−skin (Fig. 7b) is weak and also not statistically significant (ρ = −0.24, p = 0.16). The conformity index CN is weakly yet statistically significantly corretaled with R ( Fig. 7c; ρ = 0.35, p = 0.04); the correlation between CN and d CTV−skin (Fig. 7d) is very weak and not statistically significant (ρ = −0.17, p = 0.32).

Target volume location
Overall, most dosimetric parameters observed exhibit very weak to weak correlation with both R and d CTV−skin . The only exception is the mean dose to oesophagus, which exhibits moderate and statistically significant correlation with R (ρ = −0.50, p < 0.01; not shown). Similar correlation also holds for the conventional 3D CRT plans (ρ = −0.53, p < 0.01; not shown). This can be easily understood, as oesophagus usually lies close to the intersection of the superior-inferior axis and the left-right axis, so a target away from the patient origin is automatically also away from oesophagus.

Dose calculation in lungs
Being highly heterogeneous in its electron density, the thoracic region presents a challenge for dose calculation algorithms. The Anisotropic Analytical Algorithm (AAA), employed by the dose calculation engine in the Eclipse treatment planning system, is known to overestimate dose in the low density region [18,19]. This effect is more pronounced in regions with electron density lower than that of realistic lungs, small fields, and high energies. As the clinical examples studied here do not satisfy these conditions, it can be understood that our observations ( Table 2) do not confirm the forementioned results. The dose in lungs calculated by AAA does not exhibit significant deviation from its measured value. The only medium with consistent deviation observed was the bone, where the calculated dose value was consistently lower than its measured value, irrespective of the treatment technique (3D CRT vs. DCAT).

Clinical implications of dosimetric parameters
There are many challenges in delivering radiotherapy to lung cancer patients, one of the most important is ensuring low lung toxicity. Radiation induced lung injury in a form of radiation pneumonitis occurs post-treatment within six months after the treatment and affects patient morbidity and mortality to a considerable degree. Factors predictive of radiation pneumonitis include increased age, smoking history, lower lobe location of tumour, poor performance status, pulmonary disfunction, radiation dose to the lung and the volume of lung irradiated. The use of combined dose-volume metrics (such as V5Gy, V20Gy, or MLD) to evaluate plan safety is now standard practice [20,21]. The assessment of the radiation pneumonitis risk by the radiation oncologist is therefore a critical aspect of radiation therapy treatment planning. When one-third of the total lung volume is exposed to radiation doses of 20 Gy or more, the risk for symptomatic radiation induced pneumonitis is 10-15%. For V20 below 22%, the risk for radiation pneumonitis is nearly zero. Above a V20Gy of 35%, the risk for radiation pneumonitis rises precipitously. When V20Gy is increased to a value greater than 40%, the risk for radiation induced pneumonitis increases to nearly 40-50% [22,23].
Our results (Fig. 5d) show a negative correlation between improvements in V20Gy and V5Gy values for lung: if a DCAT plan brought an improvement over a conventional 3D CRT plan at V20Gy, it fared worse at V5Gy, and vice versa, the few DCAT plans which performed better than their conventional 3D CRT counterparts at V5Gy all did worse than 3D CRT at V20Gy. The relative importance of V20Gy vs. V5Gy (or MLD) as a predictor of radiation pneumonitis might serve as a factor in decision favouring conventional 3D CRT over DCAT or vice versa. However, the available clinical evidence does not yet warrant any such decision.
Overall, treatment planning restriction for lung proved to be the most difficult to meet in this study. However, we need to reiterate that the initial selection of cases eligible for the study excluded the ones in which the total dose on target was high and PTV was close to the spinal cord. The fact that all the plans met the criteria of clinical acceptability concerning the dose to spinal cord can be therefore expected. The very high ratio of treatment plans meeting the criterion for the mean dose to oesophagus is in our view affected by two factors; on the one hand, the threshold value (34 Gy) is quite high in comparison with the total doses used (Table 1), and on the other hand, while 3D CRT plans were created by a skillful dosimetrist/physicist who consciously avoided the structure, DCAT as a rule produces highly conformal dose distribution which compensates for the lack of intelligent guidance.

Interplay effect in clinical practice
One motivating factor for exploring the possible use of DCAT was the fact that, being a purely conformal technique, it avoids the interplay effect. Despite the early concern about the interplay effect between the MLC movement and intra-fraction organ motion, intensity-modulation techniques like IMRT [24], VMAT [25], helical tomotherapy [26] or even stereotactic body radiotherapy using radical hypofractionation [27] are being increasingly used for treatment of lung cancer and appear to be safe. A dosimetric study on the interplay effect during stereotactic VMAT lung treatment delivery [28] claims that the effects are well within clinically accepted tolerance levels. An overview of the interplay effect [29] agrees that the interplay effects appear to be small (1-2%) in most typical clinical cases. The dominant effect -the blurring of the dose distributionis independent of the treatment technique, and is not more pronounced in intensity-modulated than in conformal treatment techniques, although it may have a bigger impact on IMRT because of the tendency to reduce target margins when using advanced treatment techniques. All in all, even though in 20 out 35 cases in this study, field-infield technique was used in 3D CRT, we do not consider this issue as very important, in particular as the treatment fields which do not encompass the whole target projection only contribute a small fraction to the total dose.

Dosimetric impact of collimator angle
Setting the collimator angle to 45°in DCAT plans may not appear as the optimal choice with respect to the area shielded by MLC only. Indeed in our study, the average opening defined by the jaws at the collimator angle of 45°w as 27% larger than the average opening defined by the jaws at the collimator angle of 0°or 90°. Setting the collimator to either 0°or 90°however has its drawbacks too. When using a collimator angle of 0°, one and the same slice of tissue is always exposed to the interleaf leakage, while tilting the collimator results in spreading the leakage dose over a larger volume. For concave target shapes, collimator angle of 90°often results in more MLC travel than other choices. We have shown however that the dosimetric impact of the collimator angle is small. The impact of the collimator angle seems to matter more in VMAT: while the original study argued for a collimator angle set to 45° [30], in some specific cases even the collimator rotated to 90°appears to be advantegous [31].

Monitor units and treatment time
Lung cancer patients, particularly in the palliative setting, often have bad performance status, and every effort towards shortening the time they spend on the treatment table should be considered worthwhile. As a rule, DCAT treatment plans use fewer monitor units than conventional 3D CRT treatment plans (Fig. 4d). This alone means a shorter irradiation time for the patient. The total treatment time includes the overhead of patient set-up, which is equal in both the conventional 3D CRT treatment and DCAT treatment and generally exceeds the beamon time. On the other hand, the time spent due to the required gantry and/or collimator rotation between individual static treatment fields in the conventional 3D CRT increases the difference in the patient door-to-door time even further. All in all, while the ratio of MU spent for both techniques does not directly translate into the ratio of the overall treatment times, a shortening of 1-2 min is to be expected.

Conclusions
Using film dosimetry, we verified that DCAT is at least as reliable as conventional 3D CRT in terms of dosimetric accuracy. Similar tests are probably necessary at every site wishing to employ the technique, in particular if no tests specific for DCAT have been conducted during the initial machine commissioning.
In most cases, DCAT plans do not pose a dosimetric advantage over manually crafted conventional 3D CRT plans. Even though DCAT can thus not be recommended as a replacement for the conventional 3D CRT, patient anatomy and the treatment prescription may in some cases nevertheless lead to a dosimetrically favourable result with DCAT. Taking into account also a shorter overall treatment time with DCAT, we believe there is a niche for it in radiotherapy. Based on the results of this study, DCAT has been commisioned for use and introduced as a treatment option in our clinic, and is predominantly being used for palliative lung cancer patients.

Appendix: Target volume sphericity
In general, target volume is a complex 3D shape. For the purpose of this analysis, we will introduce a sphericity parameter which measures the departure of target volume shape from a sphere and does not depend on its size [32]. Sphericity of an object with a given surface area A and volume V is defined as the surface area of a sphere with a volume equal to the volume of the object, divided by the surface area of the object: Sphericity is a dimensionless parameter and can attain values from 0 to 1, with 1 being the value of a perfect sphere. The more an object departs from it, the lower its sphericity is.
While all treatment planning systems provide the volume of a defined structure as parameter, surface area is generally not easily available. We can circumvent this limitation by making use of the onion-like structure of ICRU target volumes and define: Here, V CTV and V PTV are volumes of CTV and PTV, respectively, and d is the safety margin used to construct PTV from CTV. Then, a sphericity-like parameter can be introduced as This parameter does not apply to either CTV or PTV, but to a structure in between the two. Overall,˜ behaves like sphericity, yielding high values for near-spherical shapes and low values for highly concave shapes with many protrusions. The results are the most consistent among the shapes in which the same CTV-PTV margin is applied. Unlike , however, it can exceed the value of 1 in certain cases. Despite its simplicity and limitations, we find it a useful parameter to characterize the target shape.