Skip to main content

Effectiveness of Multi-Criteria Optimization-based Trade-Off exploration in combination with RapidPlan for head & neck radiotherapy planning



A new strategy is introduced combining the use of Multi-Criteria Optimization-based Trade-Off Exploration (TO) and RapidPlan™ (RP) for the selection of optimisation parameters that improve the trade-off between sparing of organs at risk (OAR) and target coverage for head and neck radiotherapy planning. Using both approaches simultaneously; three different workflows were proposed for the optimisation process of volumetric-modulated arc therapy (VMAT) plans. The generated plans were compared to the clinical plans and the plans that resulted using RP and TO individually.


Twenty clinical VMAT plans previously administered were selected. Five additional plans were created for each patient: a clinical plan further optimised with TO (Clin+TO); two plans generated by in-house built RP models, RP_1 with the model built with VMAT clinical plans and RP_TO with the model built with VMAT plans optimised by TO. Finally, these last two plans were additionally optimised with TO for the creation of the plans RP_1 + TO and RP_TO+ respectively. The TO management was standardised to maximise the sparing of the parotid glands without compromising a clinically acceptable PTV coverage. Resulting plans were inter-compared based on dose-volume parameters for OAR and PTVs, target homogeneity, conformity, and plans complexity and deliverability.


The plans optimised using TO in combination with RP showed significantly improved OAR sparing while maintaining comparable target dose coverage to the clinical plans. The largest OAR sparing compared to the clinical plans was achieved by the RP_TO+ plans, which reported a mean parotid dose average of 15.0 ± 4.6 Gy vs 22.9 ± 5.5 Gy (left) and 17.1 ± 5.0 Gy vs 24.8 ± 5.8 Gy (right). However, at the same time, RP_TO+ showed a slight dose reduction for the 99% volume of the nodal PTV and an increase for the 95% (when comparing to the clinical plans 76.0 ± 1.2 vs 77.4 ± 0.6 and 80.9 ± 0.9 vs 79.7 ± 0.4) but remained within clinical acceptance. Plans optimised with RP and TO combined, showed an increase in complexity but were proven to be deliverable.


The use of TO combined with RP during the optimisation of VMAT plans enhanced plan quality the most. For the RP_TO+ plans, acceptance of a slight deterioration in nodal PTV allowed the largest reduction in OAR dose to be achieved.


The challenge in radiotherapy planning is to manage the compromise between tumour irradiation and sparing of healthy tissue. The planning objectives involved can be inherently contradictory, which restricts an optimal result for all objectives at the same time. This situation is reflected in the irradiation of head and neck cancer (HNC) where there are many organs at risk (OAR) and their close proximity to the target volume makes it difficult to achieve clinically acceptable target coverage without radiation-induced toxicity and decreased of quality of life [1,2,3].

For the treatment of HNC, volumetric-modulated arc therapy (VMAT) has shown good dose conformity and sparing of OAR with a decrease in treatment delivery time compared to other techniques [4,5,6,7]. For this inverse planning technique, the desired dose is prescribed for the planning target volume (PTV) and for the surrounding structures the tolerated dose is defined; then through an iterative optimisation process a plan is computed. The optimisation result is the plan which best achieves the stated objectives by minimising a cost function associated with them. However, this result is often not satisfactory since the compromise between the planning goals might not be clinically acceptable and the optimisation objectives may not adequately describe the clinical situation. Different objective settings are tried until an acceptable solution is found, which is time consuming and may result in sub-optimal treatment plans. For the solution of this inverse problem, a Multi-Criteria Optimisation (MCO) approach has been proposed as a more efficient method [8,9,10,11,12,13,14,15]. Whereas single objective function optimisation problems have one best solution, in MCO a vector of objective functions is optimised resulting in many best-compromised solutions that describe a mathematical construct, an approximation of the Pareto surface. In practice this is composed of alternative plans that are created, where each of them will prioritise a specific optimisation objective over all others. Therefore, each alternative plan is optimal in that one objective can only be improved by deteriorating others [16].

Eclipse treatment planning system (TPS) Varian Medical Systems, Palo Alto, CA, version 15.5 introduces the tools for a MCO approach based on trade-offs exploration (TO). A copy of an initial plan optimisation, called the balanced plan, is used together with the chosen N optimisation objectives as the base for the subsequent generation of the 3 N + 1 alternative plans for objective. The user can then explore the trade-offs along the possible solutions and select the plan that best fulfils the treatment goals. The initial plan influences the subsequent Pareto frontier approximation, therefore, trade-offs exploration around a starting promising plan is desirable [16, 17].

Recently, knowledge-based radiotherapy treatment planning has been incorporated into clinical practice. RapidPlan™ (RP) (Eclipse v15.5) operates with a knowledge-based treatment planning algorithm that is able to generate dose volume histogram (DVH) estimates for a new patient geometry by analysis of a database of previous plans. The DVH estimates are then translated into optimisation objectives parameters for intensity modulated radiotherapy (IMRT) or VMAT plans. The resulting plans have shown to reduce planner interaction while improving the quality and consistency [18, 19]. However, the quality of the plans generated for new patients depends highly on the quality and robustness of the plans in the training cohort [20].

In this study the key question addressed was how to maximise the performance of RP and TO for the radiotherapy treatment of HNC. We suggested a combination of both approaches, RP and TO used simultaneously during the optimisation process of a plan. Three different workflows were proposed for the VMAT planning of HNC patients and an inter-comparison of the resulting plans was performed. At the same time, these were compared with the plans optimised using RP and TO on their own and the clinical plans for reference.


Patient data

Clinical HNC plans previously administered were available from a database with tumours originating in various sites and stages, i.e. tongue, glottis, supraglottis, lip, oral cavity, nasal cavity, tonsillar region, larynx, pharynx (oro, hypo and naso), neck and parotids. Plans for the region of optical structures were not included. The clinical plans were generated with the in-house optimisation template (based on the plan constraints shown on Table 1) and expert planners modified the optimisation parameters as necessary. The planned target volumes were: the main tumour (PTV1), and neck nodes encompassed by the low-risk planning treatment volume (PTVLR). Dose levels of 65Gy and 54Gy in 30 fractions were respectively prescribed. Both PTVs were created with 3–5 mm margins from the Clinical Target Volume (CTV). Planning organ at risk volumes (PRVs) were used for spinal cord and brainstem by expanding from these structures. The planning technique was the Varian VMAT solution RapidArc™, utilising 2 full coplanar arcs with the collimator rotated to 30 and 330 degrees, and 6MV photons. Plan calculations were made using the Eclipse™ Photon Optimiser and Anisotropic Analytic Algorithm version 15.5.07. Aperture Shape Control was set to “Moderate” and all had the same Normal Tissue Objectives (NTO). The target volumes and OAR constraints followed in this study for all the sets of plans are shown on Table 1 and are based on the Phase 3 PARSPORT trial specifications [21, 22].

Table 1 Dose-volume constraints for the PTVs and OAR

Two different sets of plans were used in this work. Set_1 consisted of a collection of 54 plans that were selected randomly from the database for the creation of the RP models described in the methods section part III. These plans had a mean PTV1 volume of 325.0 cc (range 68.6–786.6 cc) and mean PTVLR volume of 182.0 cc (range 0–375.7 cc).

Set_2 included 20 different patients; these had a mean PTV1 volume of 291.8 cc (range 75.5–1069.4 cc) and a mean PTVLR volume of 274.2 cc (range 124.7–389.0 cc). Set_2 was selected from cases with no full parotid involvement (2 Tonsil, 6 Larynx, 5 Pharynx, 5 oral, and 2 face), with mean volume overlap of parotids and PTVs of 19.9%, (range 0–40.8%). Set 2 was used for the creation, comparison, and analysis of plans using TO and/or RP.

Trade-offs exploration optimisation

The inclusion of MCO in radiotherapy planning aims to allow the exploration of the trade-offs of the treatment objectives in an efficient way to then select a plan that best fulfils the prescribed clinical goals. To commence, it is required to have a starting plan (that will be the centre of the approximation of the Pareto surface). The optimisation objectives used in inverse radiotherapy planning are based on dose statistics to regions of interest, for instance:

  • Upper or lower objectives: dose-volume points that define the maximum or minimum dose that a structure may receive.

  • Mean objective: mean dose level that should not be exceeded for a structure.

  • Line objective: it is located under the lower bound of the DVH estimate range. If the OAR overlaps with the target, the highest upper objective of the target is the defining level in that region.

A starting plan will be used together with a selection of the optimisation objectives for trade-off exploration as the base for a generation of alternative plans.

The optimisation objectives selected for the trade-off exploration of this study were: right and left parotids line dose or mean dose (for the Clinical plans), PRV brainstem upper point (Dmax), PRV spinal cord combined (upper points D1 and Dmax), and larynx line dose only when the tumour did not originate in it. Upper and lower point objectives were selected for PTV1 and lower point objectives for PTVLR [23].

After the calculation of the alternative plans, a slider for each selected objective is displayed. The sliders represent the range for the objective covered by the balanced and alternative plans. A selector in the slider can be manipulated to explore the trade-offs, when it is moved to the left the associated objective improves. The manipulation of one slider automatically affects the position of other selectors depending on the trade-offs. If the planner wants to restrict the range of a specific objective a component in the slider called restrictor is used, this keeps the selector within the new range influencing how the cost is distributed during the trade-off exploration; otherwise, the algorithm aims to distribute the cost of the improvement evenly among other criteria [16]. The feedback of the exploration was constantly monitored through visual assessment of the DVHs, an axial view of the dose distribution on patient anatomy, and in the clinical protocol attached to the plan (see Fig. 5 in Appendix). Once satisfied with the result, the DVHs and the spatial distribution of the dose that resulted from the trade-offs were used as optimisation objectives parameters for the creation of a deliverable VMAT plan and the final dose was calculated.

We prioritised a dosimetric reduction on the parotid glands [24,25,26]. Therefore, on the alternative plans that were generated, the management of the trade-offs was standardised as follows:

  1. 1.

    A restrictor to stop further dose increase to spinal cord and brainstem was applied; this action still allowed their initial dose to reduce as a consequence of other structures improvement.

  2. 2.

    The mean dose to right and left parotids was reduced. The limit for this was meeting the prescribed dose constraints for the PTVs.

  3. 3.

    Mean dose to larynx was reduced with the least importance over the other organs and only when its effect would not compromise the target constraints.

To begin the study, Set_2 plans were selected as the starting plans and optimised with TO to create the Clin+TO plans following the procedure described above.

RapidPlan with trade-offs optimisation

For RapidPlan knowledge-based planning a model needs to be created. The model configuration comprises two phases: data extraction and model training. During the first phase, the data related to the structure set, dose, and field geometry of each plan included is extracted; this information is used to train the model which will then estimate DVHs for a new plan based on the information obtained from the data set of plans [27].

Once the model is trained, it will be used together with the structure set of the new plan and the dose specifications for the target to provide a DVH estimation range for each of the OAR; also it will provide the optimisation objectives parameters needed to try to achieve the estimation.

Two RP models for HNC were generated using the dosimetric and geometric data from Set_1. The first model (RP_1) was trained with the 54 standard clinical VMAT plans. The second model (RP_TO) used the same VMAT plans but this time further optimisation was carried out using TO (following the process in Section II) prior to adding them to the model. Set_1 consisted of plans previously used for the head and neck RP model that is currently in use at our clinic, therefore, these plans are clinically acceptable and the model was processed for elimination of outliers. Both models were generated by training the two PTV targets and 12 OAR structures (PRV brainstem, larynx, lenses, optic nerves, orbits, parotids, PRV spinal cord, and trachea). The optimisation objectives selected are shown in Table 5 in Appendix.

The two models were used to generate a plan for each patient of Set 2 without manual intervention, RP_1 and RP_TO. The resulting plans were then used as the initial plans for a further optimisation using TO for the creation of the last two types of plans for each patient, RP_1 + TO and RP_ TO+ (following the same trade-offs exploration procedure described on section II). In this way, the trained models suggested the optimisation objectives parameters for an initial plan that would be the centre of the later approximation of the Pareto surface. Figure 1 summarizes the workflows used for the generation of the plans created for this study. All plans from this part and section II were normalized 100% to target mean, and the dosimetric data was collected.

Fig. 1
figure 1

Planning workflow of: Clinical, Clinical+TO, RP_1, RP_1 + TO, RP_TO, and RP_TO+. (See in red proposed workflows)

Data comparison and analysis

For the clinical plans and the 5 different plans generated for each patient, PTVs and OAR doses were evaluated following the dose-volume constraints specified. Additionally, the homogeneity of the PTVs was assessed by calculating the Homogeneity Index according to Eq. 1 [28].

$$ HI=\left(\frac{D_{5\%}}{D_{95\%}}\right) $$

Also, the dose delivered to normal tissue was assessed by Eq. 2 [29].

$$ \mathrm{Conformity}\kern0.5em \mathrm{Index}\kern0.5em (CI)=\frac{Vx\%}{VPTVs} $$

Where Vx% is the volume receiving at least x% of the dose and VPTVs is the sum of both target volumes PTV1 and PTVLR. This was calculated for V100%, V80%, V50%, and V10%.

To investigate plan deliverability, RadCalc version 6.3 was used to quantify several complexity parameters: the total number of monitor units (MU), the modulation factor (MF) and the average leaf pair opening (ALPO) [30, 31].

For the plans that resulted in lowest MF, plan specific quality assurance (QA) measurements were performed at the treatment room for further evaluation of the delivered dose. A Varian TrueBeam Slim linear accelerator was used and the phantom Octavius 729 detector array together with the Octavius 4D rotating unit. A comparison between the delivered and the planned dose was performed with a 3D gamma evaluation (global and local) using the verification software VeriSoft 6.1 ( and an acceptance of 95% points passing the criteria of 3 mm for the distance to agreement (DTA) and a dose difference tolerance level of 3% [32, 33]. For all above comparative analysis, the threshold for inclusion was set at 10% of the maximum dose.

Significant differences between the 6 final sets of plans were assessed for all the investigated parameters by the Wilcoxon signed-rank test at a level of 0.05.


Dosimetric evaluation

All plans resulting from implementing TO and/or RP were deemed clinically acceptable according to the criteria of PTV coverage and OAR doses (Table 1). While maintaining comparable target coverage, the sparing of the OAR improved. Table 2 lists the dosimetric parameter values (average ± standard deviation) obtained for the PTVs and OAR for the HNC plan sets following each of the 5 planning procedures described, and for the clinical plans; the letters indicate significant difference at a level of 0.05 when inter-comparing these techniques. The HI, CI, and the number of parameters that failed to meet the constraints for each procedure, and the number of plans that correspond to these failures, are also displayed (for further detail see Table 6 in Appendix). Additionally, doses were investigated for larynx and trachea and are reported in Table 7 in Appendix.

Table 2 Summary of average ± s.d. doses to PTVs and OAR for the 20 patients in Set_2

For the plans created, it is shown in Table 2 that PTV1 reported no statistically significant dose reduction compared to the clinical plans for its 99 and 95% volumes (D99, D95) and the homogeneity index remained comparable across techniques as well. For PTVLR, the dose homogeneity was shown to be slightly affected and D95 increased significantly but a D99 significant reduction resulted only for RP_TO+.

Figure 2 presents the dose differences obtained between the clinical plan and the other five plans created for each of the 20 patients. It is shown that the largest OAR reductions were achieved when RP and TO were implemented together, with the largest sparing for brainstem, parotid glands, and spinal cord being achieved by the RP_TO+ plans. For OAR, RP_1 + TO and RP_TO reported no significant difference between them. Furthermore, the combined use of RP and TO showed improvements in CI when comparing with the solely use of TO (Clin+TO plans) however whereas RP_TO plans are comparable with the clinical plans, for RP_1 + TO and RP_TO+ the results showed an increase. The dose distribution of a representative case is illustrated in Fig. 3, which displays all plans at the parotid level. Figure 4 shows the DVHs corresponding to the patient plans shown in Fig. 3.

Fig. 2
figure 2

Box plot of the difference of dosimetric parameters (PTVs-above and OAR-below) between clinical and generated plans. The mean value is indicated by the symbol, the median value by the central horizontal line, the interquartile range is represented by the box, and the outliers are indicated by the asterisks

Fig. 3
figure 3

Dose distribution at the parotids level for a representative case resulting from each planning optimisation procedure

Fig. 4
figure 4

DVHs corresponding to representative case shown in Fig. 3. Clinical (suns), Clinical+TO (hearts), RP_1 (circles), RP_1 + TO (rhombus), RP_TO (triangles), RP_TO+ (squares)

Complexity parameters evaluation

Table 3 displays the average values (and range) of the complexity parameters assessed for each group of plans. The clinical, RP_1 and Clin+TO plans, when compared, showed no significant differences in the number of MU or the MF. However, the use of TO and RP together in the three remaining types of plans led to an increase in MU and a decrease in the MF, indicating a higher complexity of the plan. As a result of this observation the two plans with the lowest MF from each of these plan categories (RP_1 + TO, RP_TO, and RP_ TO+) were selected for plan specific QA. Table 4 presents the analysis of the plan QA including the percentage of points passing the gamma criteria for each of them and also for their corresponding clinical plan for reference. All the plans subject to QA were reported deliverable according to the acceptance criteria given above.

Table 3 Complexity parameters for each group of plans (average and range)
Table 4 Percentage of points passing the 3%/3 mm gamma criteria


In this study, the combined use of RP and TO was proposed focusing on HNC radiotherapy VMAT planning. Three workflows were followed using RP together with TO during the optimisation process: the use of an in- house RP model built with TO VMAT plans (RP_TO); the use of an in-house RP model, built with clinical VMAT plans, to obtain the starting optimisation objectives parameters and then further TO (RP_1 + TO); and the use of RP_TO to generate the initial objectives parameters and further TO (RP_TO+). The dosimetric impact was assessed on the resulting plans for a set of 20 patients. In addition, for each patient, two more plans were created using TO and RP on its own during the optimisation stage (Clin+TO, RP_1). For the five optimisation strategies and the clinical plans, a dosimetric comparison across the resulting plans was carried out. Additionally, investigation of complexity and deliverability of the plans was performed.

Although the use of RP and TO separately presents advantages over the clinical plans [11, 19], it was shown that for this treatment the combined use of both optimisation methods improves the relationship between balancing OAR doses and PTV coverage, i.e. while maintaining comparable clinically acceptable target coverage, the resulting plans reduced the OAR doses. It was demonstrated that the RP models are able to provide a semi-automatic suggestion of a plan that is a good start for TO. The OAR doses of the plans generated by the RP_TO model, built with previously TO optimised plans, reflect its benefit over the RP_1 model. Therefore, this approach suggests being convenient for centres that does not have TO but that can implement shared RP models. When comparing RP_TO and RP_1 + TO plans no significant difference was reported on OAR sparing, but the first one reduced active planner interaction. Furthermore, when RP_TO is used as the initial plan with subsequent use of TO (RP_TO+) this led to the largest dose reduction for all OAR. This reveals that for a new patient, the knowledge-based model that was trained with the plans that were already optimised by TO helped to obtain optimisation parameters that are closer to the desired trade-offs, providing the best starting point for TO. The results are in agreement with Wang J et al. [17] who showed for nasopharyngeal treatments that choosing better optimisation parameters for a plan with the aid of automated planning prior to use MCO led to improved plans.

RP_ TO+ plans achieved the lowest average OAR doses for all the dose-volume parameters investigated. Most notably, a significant reduction in mean dose of 7.9 Gy and 7.7 Gy for the left and right parotids respectively was accomplished with respect to the clinical plans. Despite the dose improvement for the OAR, the CI results showed that the dose to normal tissue was slightly compromised. Further work with NTO and the inclusion of a structure, considering normal tissue not being recognised as OAR, as part of the trade-offs exploration objectives is of our interest to investigate improvement [19]. On the contrary, OAR doses with RP_TO plans are higher than with RP_TO+ but this plans provided an OAR dose improvement respect to clinical plans while maintaining comparable CI and HI to them, being this influenced by the upper PTVLR objectives used in the RP model. The results from both workflows point to the capacity of TO to decide on the trade-offs pursued for this cohort of patients and agree with ongoing HNC trials that have shown decreased toxicity to the surrounding organs as a result of dose de-escalation to nodal PTV. The implementation of RP and TO for the treatment of this cohort of patients could lead to an improved patient quality of life [34, 35].

The results displayed in Table 2 support the value of TO as a tool for attaining further individualised plan improvements. It is shown that the application of TO as the final step in planning (Clin+TO, RP_1 + TO and RP_TO+) resulted in a consistent reduced number of plans failing to meet the OAR constraints. The reason three plan parameters failed with RP_TO+ as opposed to two with Clin+TO is that the additional failing plan had the largest PTV1 in the set (1069.0 cc) and this lay outside the range of the plans that constitute the RP models used; therefore, following our trade-offs management, PTV1 D99 and PTVLR D99 were marginally below the constraint (94.9 and 74.7% respectively). This patient was re-planned with different trade-offs exploration management and the constraints of the PTVs were met at the expense of a 2.4Gy increase in the left parotid mean dose. This example points to the fact that the success of a RP model is dependent upon the quality and robustness of the model itself [20]. Similarly, in the use of TO there is an inherent element of user influence.

The time required to create these plans is equally of relevance. The estimated time in our clinics to plan a HNC case, including optimisation and dose calculation until the final plan is produced, ranges between 50 and 300 min depending on the complexity and the experience of the planner. The use of RP and MCO in clinical practice has proven to reduce planning time [10, 36]. Following the proposed workflows, with the combined use of both approaches, the approximate time taken remained always under 80 min, being affected mostly by the speed of the calculations and the number of optimisation objectives selected rather than by the human time needed or the planner experience. Therefore, using both RP and TO at the same time did not increase our human time at the console, and overall is expected to reduce it as less interaction to reach a clinically acceptable solution is needed.

Kyroudi, et al. [37] reported dosimetric differences between the exploration of the alternative plans generated with MCO and the deliverable plans that are created after this, which take into account physical treatment unit restrictions. In this study, small discrepancies were noticed as well but they were followed by a repetition of the trade-off exploration only in a few cases when we were already in the limit of the dose required to meet the PTVs constraints. Furthermore, an evaluation of the dosimetric delivery parameters was performed to warranty plan deliverability. It was found that the complexity of the plans increased when using RP and TO together. This could be expected as the sparing of the OAR has been improved providing higher quality plans [38]. The highest mean MU value and lowest mean MF and ALPO values were observed for the RP_TO plans suggesting that these plans are more complex than the RP_1 + TO and RP_TO+. Although RP_TO+ showed better relationship between OAR sparing and PTV coverage, it can be observed that the homogeneity of the targets is slightly better when using RP_TO which has proved to have relation with an increase in complexity [39]. Despite the higher complexity of the plans that used RP and TO, the observed values were within the range of the clinical plans. Further, for the extreme values of the three proposed workflows, in most cases the results from the gamma evaluation showed improvement on the agreement of the delivered and planned dose to the OAR, and concluded that the plans are deliverable clinically.

Future work will focus on the implementation of the RP_TO and RP_TO+ workflows in clinical practice and analysis of the physicians’ preference on the resulting trade-offs. The use of the proposed optimisation approaches for other treatment sites is of interest and the opportunity to use RP and TO in combination to escalate the dose delivered to the tumour while maintaining or reducing OAR doses.


The use of RP combined with TO improved OAR sparing on VMAT HNC radiotherapy plans while maintaining clinically acceptable target coverage. The simultaneous use of both approaches reflects in the resulting plans the individual benefit that they provide. The best OAR sparing was obtained when planning starting by generating a plan with a RP model built with plans optimised with TO and then further individualised TO optimisation. For this workflow, the trade-off between a minor deterioration of the nodal PTV, while maintaining clinically acceptable PTVs coverage, allowed the most significant reduction in OAR doses. RP_TO plans, created with the RP model optimised with TO plans, also provided an OAR dose improvement respect to clinical plans while maintaining comparable conformity and target homogeneity. Plans optimised with RP and TO were proven to be deliverable despite an increase in plan modulation.



Average leaf pair opening


Clinical target volume


Distance to agreement


Dose-volume histogram


Homogeneity index


Head and neck cancer


Intensity modulated radiotherapy


Multi-criteria optimisation


Modulation factor


Monitor units


Organs at risk


Planning organ at risk volume


Planning target volume


Main planning target volume


Low risk planning target volume


Quality assurance




Multi-Criteria Optimisation-based Trade-Off Exploration


Treatment planning system


Volumetric modulated arc therapy


  1. Cooper JS, Fu K, Marks J, Silverman S. Late effects of radiation therapy in the head and neck region. Int J Radiat Oncol Biol Phys. 1995;31:1141–64.

    CAS  Article  Google Scholar 

  2. Langendijk JA, Doornaert P, Verdonck-de Leeuw IM, Leemans CR, Aaronson NK, Slotman BJ. Impact of late treatment-related toxicity on quality of life among patients with head and neck cancer treated with radiotherapy. J Clin Oncol. 2008;26:3770–6.

    Article  Google Scholar 

  3. Deasy JO, Moiseenko V, Marks L, Chao KC, Nam J, Eisbruch A. Radiotherapy dose–volume effects on salivary gland function. Int J Radiat Oncol Biol Phys. 2010;76:S58–63.

    Article  Google Scholar 

  4. Holt A, Van Gestel D, Arends MP, Korevaar EW, Schuring D, Kunze Busch MC, et al. Multi-institutional comparison of volumetric modulated arc therapy vs. intensity-modulated radiation therapy for head-and-neck cancer: a planning study. Radiat Oncol. 2013;8:26.

    Article  Google Scholar 

  5. Verbakel WF, Cuijpers JP, Hoffmans D, Bieker M, Slotman BJ, Senan S. Volumetric intensity-modulated arc therapy vs. conventional IMRT in head-and-neck cancer: a comparative planning and dosimetric study. Int J Radiat Oncol Biol Phys. 2009;74:252–9.

    Article  Google Scholar 

  6. Vanetti E, Clivio A, Nicolini G, Fogliata A, Ghosh-Laskar S, Agarwal JP. Et at. Volumetric modulated arc radiotherapy for carcinomas of the oro-pharynx, hypo-pharynx and larynx: a treatment planning comparison with fixed field IMRT. Radiother Oncol. 2009;92:111–7.

    Article  Google Scholar 

  7. Clemente S, Wu B, Sanguineti G, Fusco V, Ricchetti F, Wong J, et al. SmartArc-based volumetric modulated arc therapy for oropharyngeal cancer: a dosimetric comparison with both intensity-modulated radiation therapy and helical tomotherapy. Int J Radiat Oncol Biol Phys. 2011;80:1248–55.

    Article  Google Scholar 

  8. Craft DL, Halabi TF, Shih HA, Bortfeld TR. Approximating convex Pareto surfaces in multiobjective radiotherapy planning. Med Phys. 2006;33:3399–407.

    Article  Google Scholar 

  9. Craft D, McQuaid D, Wala J, Chen W, Salari E, Bortfeld T. Multicriteria VMAT optimization. Med Phys. 2012;39:686–96.

    Article  Google Scholar 

  10. Craft DL, Hong TS, Shih HA, Bortfeld TR. Improved planning time and plan quality through multicriteria optimization for intensity-modulated radiotherapy. Int J Radiat Oncol Biol Phys. 2012;82:e83–90.

    Article  Google Scholar 

  11. Kierkels RG, Visser R, Bijl HP, Langendijk JA, Veld v‘t, AA SRJ, et al. Multicriteria optimization enables less experienced planners to efficiently produce high quality treatment plans in head and neck cancer radiotherapy. Radiat Oncol. 2015;10:87.

    Article  Google Scholar 

  12. Wala J, Craft D, Paly J, Zietman A, Efstathiou J. Maximizing dosimetric benefits of IMRT in the treatment of localized prostate cancer through multicriteria optimization planning. Med Dosimetry. 2013;38:298–303.

    Article  Google Scholar 

  13. Kamran SC, Mueller BS, Paetzold P, Dunlap J, Niemierko A, Bortfeld T, et al. Multi-criteria optimization achieves superior normal tissue sparing in a planning study of intensity-modulated radiation therapy for RTOG 1308-eligible non-small cell lung cancer patients. Radiother Oncol. 2016;118:515–20.

    Article  Google Scholar 

  14. Müller BS, Shih HA, Efstathiou JA, Bortfeld T, Craft D. Multicriteria plan optimization in the hands of physicians: a pilot study in prostate cancer and brain tumors. Radiat Oncol. 2017;12:168.

    Article  Google Scholar 

  15. Thieke C, Küfer KH, Monz M, Scherrer A, Alonso F, Oelfke U, et al. A new concept for interactive radiotherapy planning with multicriteria optimization: first clinical evaluation. Radiother Oncol. 2007;85:292–8.

    Article  Google Scholar 

  16. Varian Medical Systems, Inc. Eclipse photon and Electron algorithms reference guide. Palo Alto CA; 2017. p.221–31.

  17. Wang J, Chen Z, Li W, Qian W, Wang X, Hu W. A new strategy for volumetric-modulated arc therapy planning using AutoPlanning based multicriteria optimization for nasopharyngeal carcinoma. Radiat Oncol. 2018;13:94.

    Article  Google Scholar 

  18. Tol JP, Delaney AR, Dahele M, Slotman BJ, Verbakel WF. Evaluation of a knowledge-based planning solution for head and neck cancer. Int J Radiat Oncol Biol Phys. 2015;91:612–20.

    Article  Google Scholar 

  19. Fogliata A, Reggiori G, Stravato A, Lobefalo F, Franzese C, Franceschini D, et al. RapidPlan head and neck model: the objectives and possible clinical benefit. Radiat Oncol. 2017;12:73.

    CAS  Article  Google Scholar 

  20. Hussein M, South CP, Barry MA, Adams EJ, Jordan TJ, Stewart AJ, et al. Clinical validation and benchmarking of knowledge-based IMRT and VMAT treatment planning in pelvic anatomy. Radiother Oncol. 2016;120:473–9.

    Article  Google Scholar 

  21. Nutting CM, Morden JP, Harrington KJ, Urbano TG, Bhide SA, Clark C, et al. Parotid-sparing intensity modulated versus conventional radiotherapy in head and neck cancer (PARSPORT): a phase 3 multicentre randomised controlled trial. Lancet Oncol. 2011;12:127–36.

    Article  Google Scholar 

  22. Urbano MG, Clark CH, Kong C, Miles E, Dearnaley DP, Harrington KJ, et al. Target volume definition for head and neck intensity modulated radiotherapy: pre-clinical evaluation of PARSPORT trial guidelines. Clinical Oncol. 2007;19:604–13.

    Article  Google Scholar 

  23. Varian Medical Systems, Inc. TPS New Features Workbook 15.5: users manual for Varian treatment planning software. Palo Alto: Varian Medical Systems, Inc; 2017. p. 104–15.

  24. Bjordal K, Kaasa S, Mastekaasa A. Quality of life in patients treated for head and neck cancer: a follow-up study 7 to 11 years after radiotherapy. Int J Radiat Oncol Biol Phys. 1994;28:847–56.

    CAS  Article  Google Scholar 

  25. Harrison LB, Zelefsky MJ, Pfister DG, Carper E, Raben A, Kraus DH, et al. Detailed quality of life assessment in patients treated with primary radiotherapy for squamous cell cancer of the base of the tongue. Head & Neck. 1997;19:169–75.

    CAS  Article  Google Scholar 

  26. Malouf JG, Aragon C, Henson BS, Eisbruch A, Ship JA. Influence of parotid-sparing radiotherapy on xerostomia in head and neck cancer patients. Cancer Detect Prev. 2003;27:305–10.

    Article  Google Scholar 

  27. Fogliata A, Belosi F, Clivio A, Navarria P, Nicolini G, Scorsetti M, et al. On the pre-clinical validation of a commercial model-based optimisation engine: application to volumetric modulated arc therapy for patients with lung or prostate cancer. Radiother Oncol. 2014;113:385–91.

    Article  Google Scholar 

  28. Kataria T, Sharma K, Subramani V, Karrthick KP, Bisht SS. Homogeneity index: an objective tool for assessment of conformal radiation treatments. J Med Phys/Association of Medical Physicists of India. 2012;37:207.

    Google Scholar 

  29. Feuvret L, Noël G, Mazeron JJ, Bey P. Conformity index: a review. Int J Radiat Oncol Biol Phys. 2006;64:333–42.

    Article  Google Scholar 

  30. Svensson E. Evaluation of complexity and deliverability of IMRT-treatment plans (Master thesis). Göteborg: The Sahlgrenska Academy; 2011.

    Google Scholar 

  31. Olofsson N. Evaluation of IMRT beam complexity metrics to be used in the IMRT QA process (Master thesis). Göteborg: The Sahlgrenska Academy; 2011.

    Google Scholar 

  32. Syamkumar SA, Padmanabhan S, Sukumar P, Nagarajan V. Characterization of responses of 2d array seven29 detector and its combined use with octavius phantom for the patient-specific quality assurance in rapidarc treatment delivery. Med Dosimetry. 2012;37:53–60.

    CAS  Article  Google Scholar 

  33. Alber M, Broggi S, De Wagter C, Eichwurzel I, Engström P, Fiorino C, et al. Guidelines for the verification of IMRT. ESTRO booklet. 2008.

  34. Nuyts S, Lambrecht M, Duprez F, Daisne JF, Van Gestel D, Van den Weyngaert D, et al. Reduction of the dose to the elective neck in head and neck squamous cell carcinoma, a randomized clinical trial using intensity modulated radiotherapy (IMRT). Dosimetrical analysis and effect on acute toxicity. Radiother Oncol. 2013;109:323–9.

    Article  Google Scholar 

  35. Mirghani H, Amen F, Blanchard P, Moreau F, Guigay J, Hartl DM, et al. Treatment de-escalation in HPV-positive oropharyngeal carcinoma: ongoing trials, critical issues and perspectives. Int J Cancer. 2014;136:1494–503.

    Article  Google Scholar 

  36. Chang AT, Hung AW, Cheung FW, Lee MC, Chan OS, Philips H, Cheng YT, Ng WT. Comparison of planning quality and efficiency between conventional and knowledge-based algorithms in nasopharyngeal cancer patients using intensity modulated radiation therapy. Int J Radiat Oncol Biol Phys. 2016;95:981–90.

    Article  Google Scholar 

  37. Kyroudi A, Petersson K, Ghandour S, Pachoud M, Matzinger O, Ozsahin M, et al. Discrepancies between selected Pareto optimal plans and final deliverable plans in radiotherapy multi-criteria optimization. Radiother Oncol. 2016;120:346–8.

    Article  Google Scholar 

  38. Craft D, Süss P, Bortfeld T. The tradeoff between treatment plan quality and required number of monitor units in intensity-modulated radiotherapy. Int Journal Radiat Oncol Biol Phys. 2007;67:1596–605.

    Article  Google Scholar 

  39. Shen L, Chen S, Zhu X, Han C, Zheng X, Deng Z, et al. Multidimensional correlation among plan complexity, quality and deliverability parameters for volumetric-modulated arc therapy using canonical correlation analysis. J Radiat Research. 2018;59:207–15.

    Article  Google Scholar 

Download references


The authors would like to thank Peter Houston and Susan Morris for the valuable comments to this work.


No funding

Availability of data and materials

The datasets used and analysed during the current study are available from the corresponding author on reasonable request.

Author information

Authors and Affiliations



Study design: GC, SC. Treatment cases planning: EM, GC, SC. Knowledge-based models construction: AJ, EM, GC, SC. Dosimetry measurements: AJ, EM. Data acquisition, statistical analysis and manuscript drafting: EM. Interpretation of data: AJ, EM, GC, SC. Manuscript editing: AJ, EM, GC, SC. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Eliane Miguel-Chumacero.

Ethics declarations

Ethics approval and consent to participate

Consent for the use of patients’ data for retrospective and scientific investigation.

Consent for publication

Consent for the use of patients’ data for retrospective and scientific investigation.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.



Fig. 5
figure 5

Eclipse real time plan navigation screen view during TO for one of the head and neck plans. 1-Yellow box: Trade-Offs exploration with slider bars for each selected objective. 2-Blue box: Graphical feedback of the modifications instantly monitored in the DVHs and the dose distribution view displayed. 3-Red box: Plan objectives displaying the protocol constraints for dose-volume monitoring

Table 5 Structures and Objectives used for the RapidPlan models. The rest of the OAR (Lenses, Optic Nerves, Orbits, Optic Chiasm, Larynx, Parotids and Trachea) have only a Line (preferring target) objective with volume, dose and priority Generated
Table 6 Number of parameter violations occurred in each group of head and neck plans
Table 7 Mean doses received for larynx and trachea by each planning method when PTV did not originate in these sites

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Miguel-Chumacero, E., Currie, G., Johnston, A. et al. Effectiveness of Multi-Criteria Optimization-based Trade-Off exploration in combination with RapidPlan for head & neck radiotherapy planning. Radiat Oncol 13, 229 (2018).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Multi-criteria optimisation
  • Trade-offs exploration
  • Head and neck cancer
  • Knowledge-based treatment planning
  • RapidPlan
  • Volumetric-modulated arc therapy