% [D_full, w_fin, Geometry, optGoal] = NLP_beamlet_optimizer;

function compare_RO_norm_DVH_colorwash(D_full)

% path_in='C:\010-work\003_localGit\WiscPlan_v2\data\PatientData\';
% path_in='C:\010-work\003_localGit\WiscPlan_v2\data\Tomo_Phantom\';
% path_in='C:\010-work\003_localGit\WiscPlan_v2\data\PatientData_dog5_3\';
% path_in='C:\010-work\003_localGit\WiscPlan_v2\data\PatientData_ozzy1';
path_in = 'C:\010-work\003_localGit\WiscPlan_v2\data\PatData_ausGli_005_16beam';




load([path_in '\matlab_files\optResults.mat'])
D_norm = optResults.dose{end};

orthoslice(D_full-D_norm, [-30, 30])

% orthoslice(D_full, [0,70])
% colorwash(Geometry.data, D_full, [500, 1500], [0,80])
% colorwash(Geometry.data, D_norm, [500, 1500], [0,80])

end