WiscPlan_v2/planEvaluation/Compare_optPlans.m

function Compare_optPlans
    
    resultPath = 'C:\010-work\003_localGit\WiscPlan_v2\data\PatientData_ozzy1\matlab_files\';
    
    load([resultPath 'Geometry.mat'])
    load([resultPath 'NLP_result.mat'])
    load([resultPath 'optResults.mat'])
    
    roi_idx = 1;
    nfrac = 1;
    dose_WP = optResults.dose{end};
    dose_NLP = NLP_result.dose;

    [dvh_WP dosebins_WP] = dvhist(dose_WP,Geometry,roi_idx,nfrac);
    [dvh_NLP dosebins_NLP] = dvhist(dose_NLP,Geometry,roi_idx,nfrac);
    
    figure
    hold on
    plot(dosebins_WP, dvh_WP,'Color', [0.9,0.2,0.2],'LineStyle', '-','DisplayName', Geometry.ROIS{roi_idx}.name);
    plot(dosebins_NLP, dvh_NLP,'Color', [0.2,0.9,0.2],'LineStyle', '-','DisplayName', Geometry.ROIS{roi_idx}.name);
    hold off
    legend('Classical', 'Robust')
end