A robustness optimization is disclosed for a broad beam proton therapy plan. A nominal dose distribution (62) is computed, to be delivered to a volume by performing proton therapy on the volume according to a proton therapy plan (50) having parameters (52) defining a range compensator shape and a range band. The parameters of the proton therapy plan are adjusted to reduce a difference between the nominal dose distribution (62) and a perturbed dose distribution calculated to be delivered to the volume modified by an error scenario (64) by performing proton therapy on the volume modified by the error scenario in accordance with the proton therapy plan with the parameters adjusted by the adjusting operation. The adjusting may be repeated serially for each error scenario of a set of error scenarios (44) to produce a proton therapy plan that is robust for any of these error scenarios.La présente invention concerne une optimisation de la robustesse pour un plan de thérapie par protons à faisceau large. Une distribution de dose nominale (62) est calculée, destinée à être administrée à un volume par réalisation d'une thérapie par protons sur le volume selon un plan de thérapie par protons (50) ayant des paramètres (52) définissant une forme de compensateur de portée et une bande de portée. Les paramètres du plan de thérapie par protons sont réglés pour réduire une différence entre la distribution de dose nominale (62) et une distribution de dose perturbée calculée à administrer au volume modifié par un scénario d'erreur (64) en effectuant une thérapie par protons sur le volume modifié par le scénario d'erreur conformément au plan de thérapie par protons avec les paramètres réglés par l'opération de réglage. Le réglage peut être répété en série pour chaque scénario d'erreur d'un ensemble de scénarios d'erreur (44) afin de produire un plan de thérapie par protons qui est robuste pour l'un quelconque de ces scénarios d'erreur.