%; Time-stamp: <2002-03-05 16:58:00 dph> 
%; File: Setup.sl
%; Author: D. Huenemoerder
%; Original version: 2002.03.04
%;====================================================================
% version: 1
%
% purpose:  set up responses, data for simulation and fit examples.
%  

add_to_isis_load_path("./Sl/");
add_to_isis_module_path("./Sl/");

%%% some required files:
%
autoload("Define_model",       "Define_model.sl");   % ad hoc model plasma
require("Define_model");
autoload("splot_data_counts",  "splot_data.sl");  % smoothed data plotting
autoload("splot_data_flux",    "splot_data.sl");
autoload("osplot_data_counts", "splot_data.sl");
autoload("osplot_data_flux",   "splot_data.sl");
autoload("shplot",             "splot_data.sl");
autoload("oshplot",            "splot_data.sl");

require("gsmooth");

evalfile("pick_el.sl");              % utility for picking elements from a list
evalfile("p_groups.sl");             % ultility for plotting different elements
evalfile("custom_plot.sl");          % ad hoc plot utility

evalfile("Trans_List_define.sl");    % Utility function: define line transitions.
evalfile("aped_fit_models.sl");      % gotta eval this to establish a variable for...
evalfile("get_line_fluxes.sl");

%
variable T_exp = 100000.0;                 % exposure time.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% Define ARF files:
%

variable fgarfs_heg = [
  "Ard/acism00006_004N000HEG_-1_garf.fits",
  "Ard/acism00006_004N000HEG_1_garf.fits"
];

variable fgarfs_meg = [
  "Ard/acism00006_004N000MEG_-1_garf.fits",
  "Ard/acism00006_004N000MEG_1_garf.fits"
];


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% Define RMF files:
%
variable frmf_heg="Ard/acisheg1D1999-07-22rmfN0004.fits";
variable frmf_meg="Ard/acismeg1D1999-07-22rmfN0004.fits";


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% Load rmfs, define indices, define number of orders
%
load_rmf(frmf_heg); variable rmf_heg=1,  nheg = length(fgarfs_heg);
load_rmf(frmf_meg); variable rmf_meg=2,  nmeg = length(fgarfs_meg);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% Load arfs, assign arf, rmf:
%
variable i;

for (i=0; i<length(fgarfs_heg); i++)
 {
 ()=load_arf(fgarfs_heg[i]);
 assign_arf(i+1,i+1);
 assign_rmf(rmf_heg, i+1);
 set_arf_exposure(i+1, T_exp);                    
 }

for (i=0; i<length(fgarfs_meg); i++)
 {
 ()=load_arf(fgarfs_meg[i]);
 assign_arf(nheg + i+1, nheg + i+1);
 assign_rmf( rmf_meg, nheg + i+1);
 set_arf_exposure(nheg + i+1, T_exp);                    
 }

%+ get data grids
%
variable   xm1 = (get_data_counts(3)).bin_lo;
variable   xm2 = (get_data_counts(4)).bin_hi;
variable   xh1 = (get_data_counts(1)).bin_lo;
variable   xh2 = (get_data_counts(2)).bin_hi;
%-  

variable p1=open_plot("/xwin"); resize(20,0.6);  erase; % default window
variable p2=open_plot("/xwin"); resize(20,0.6);  erase; % extra window
variable p3=open_plot("/xwin"); resize(20,0.6);  erase; % another extra window
window(p1) ; % default

plasma(aped);