%
% file: analysis_global.txt
%

% Experiment with global fit to E0102 data...

%------------------
% get the data:
.source e0102_data

% choose ones to use, 3828 and zo for now:
e2d_ignore(1);
e2d_ignore(2);
e2d_ignore(5);
e2d_ignore(6);
e2d_ignore(12);
e2d_ignore(13);

e2d_meta[3].oversim=16.0;
e2d_meta[4].oversim=16.0;
e2d_meta[7].oversim=16.0;
e2d_meta[8].oversim=16.0;
e2d_meta[10].oversim=9.0;
e2d_meta[11].oversim=9.0;

e2d_list_data;


%------------------
% get the model
.source e0102_cal_model


% evaluate the model
e2d_pars2fom;



% Make comparison plots:

variable PSwin=open_plot("e0102_2Tcont_whetg.ps/VCPS",2,2);

% Make the 4 1d plots:
e2d_plot_1d(3,1);
e2d_plot_1d(4,1);
e2d_plot_1d(7,1);
e2d_plot_1d(8,1);

close_plot(PSwin);


variable PSwin=open_plot("e0102_2Tcont_wzeroth.ps/VCPS",1,2);

% Make the  1d plot:
e2d_plot_1d(11,2);

close_plot(PSwin);


% Output the residual 2D images
% Set them to the same wavelength range
% and use large bin sizes for more S/N.

e2d_set_view(3, [ 3 , 23 , 10.0*0.011 ] , [ -30 , 30 , 10.0*0.492 ] ); 
e2d_view_resid(3);

e2d_set_view(4, [ 3 , 23 , 10.0*0.011 ] , [ -30 , 30 , 10.0*0.492 ] ); 
e2d_view_resid(4);

e2d_set_view(7, [ 3 , 23 , 10.0*0.011 ] , [ -30 , 30 , 20.0*0.492 ] );
e2d_view_resid(7);

e2d_set_view(8, [ 3 , 23 , 10.0*0.011 ] , [ -30 , 30 , 20.0*0.492 ] );
e2d_view_resid(8);

% Do a Save to file names:  e0102_2Tcont_megm [-1.png etc.]