% file: x1822_data.sl
% 
% Load HETG data into Event-2D for analysis...
% 
% 1/6/08 - dd
% 

variable data_dir = "../Data/obs_671/";
variable yaxis_str;

% Select a y-axis value for the data sets:
%% yaxis_str = "TG_D";
yaxis_str = "TIME";

% MEG -1
e2di=1;
e2d_load_data(data_dir+"evt2","TG_LAM",yaxis_str,[0.,0.],[2,-1],[0.3,9.0]);
e2d_setup_arf(e2di,data_dir+"meg_-1.arf");
% MEG +1
e2di=2;
e2d_load_data(data_dir+"evt2","TG_LAM",yaxis_str,[0.,0.],[2,+1],[0.3,9.0]);
e2d_setup_arf(e2di,data_dir+"meg_1.arf");
% HEG -1
e2di=3;
e2d_load_data(data_dir+"evt2","TG_LAM",yaxis_str,[0.,0.],[1,-1],[0.3,9.0]);
e2d_setup_arf(e2di,data_dir+"heg_-1.arf");
% HEG +1
e2di=4;
e2d_load_data(data_dir+"evt2","TG_LAM",yaxis_str,[0.,0.],[1,+1],[0.3,9.0]);
e2d_setup_arf(e2di,data_dir+"heg_1.arf");

% Combine all of them into one
e2di=5;
e2d_load_data(data_dir+"evt2","TG_LAM",yaxis_str,[0.,0.],[2,-1],[0.3,9.0]);
e2d_setup_arf(e2di,data_dir+"meg_-1.arf");
e2d_addin_data(5,2); 
e2d_addin_data(5,3); 
e2d_addin_data(5,4);
% and set the model events equal to the real events
% to allow plotting
e2d_data[5].mxs = e2d_data[5].dxs;
e2d_data[5].mys = e2d_data[5].dys;
e2d_meta[5].modvalid=1;

% Set the nominal viewing ranges:
if (yaxis_str == "TIME") 
{
   e2d_set_view(1,[1.0,23.0,0.2],[0.0,42.0,0.40]);
   e2d_set_view(2,[1.0,23.0,0.2],[0.0,42.0,0.40]);
   e2d_set_view(3,[1.0,15.0,0.1],[0.0,42.0,0.40]);
   e2d_set_view(4,[1.0,15.0,0.1],[0.0,42.0,0.40]);
   e2d_set_view(5,[1.0,23.0,0.1],[0.0,42.0,0.40]);
}
else
{
   % TG_D:
   e2d_set_view(1,[1.0,23.0,0.02],[-5.0,5.0,0.1]);
   e2d_set_view(2,[1.0,23.0,0.02],[-5.0,5.0,0.1]);
   e2d_set_view(3,[1.0,15.0,0.01],[-5.0,5.0,0.1]);
   e2d_set_view(4,[1.0,15.0,0.01],[-5.0,5.0,0.1]);
   e2d_set_view(5,[1.0,23.0,0.01],[-5.0,5.0,0.1]);
}


% and the zeroth order events too:
% X and Y
e2di=10;
e2d_load_data(data_dir+"evt2","X","Y",[4095.73,4169.74],[0,0],[0.3,9.0]);
e2d_setup_arf(e2di,data_dir+"HETG_0.arf");
e2d_set_view(e2di, [-5.0,5.0,0.1], [-5.0,5.0,0.1]);
% ENERGY and TIME
e2di=11;
e2d_load_data(data_dir+"evt2","ENERGY","TIME",[0.0,0.0],[0,0],[0.3,9.0]);
e2d_setup_arf(e2di,data_dir+"HETG_0.arf");
e2d_set_view(e2di, [0.3,8.0,0.1], [0.0,42.0,0.40]);


% Finish by listing the data:
e2d_list_data;