#! /usr/bin/env /nfs/cxc/a1/i686/opt/isis-devel/bin/i686/isis

% Time-stamp: <2001-10-16 13:08:42 dph> 
% MIT Directory: ~dph/libisis
% CfA Directory: ~dph/libisis
% File: hcspec_summary
% Author: D. Huenemoerder
% Original version: 2001.1015

%  *** check/edit paths to env and isis !!!!!!!
%
%  on HEAD network, isis is in /soft/isis/isis/bin/$ARCH/isis
%  on CSR network, linux isis is in  /nfs/cxc/a1/i686/opt/isis/bin/$ARCH/isis
%  on my laptop: /usr/local/bin/isis

%====================================================================

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% env wrapper for cspec_summary.sl
%
%  Make summary plot of HETGS spectrum from type II pha file.
%
% REQUIRES cspec_summary.sl to be in your slang load path. (or pwd)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

variable NARGS, fpha, fout, bfact, temps, wgts;

NARGS = length(__argv);
if (NARGS < 3)
{
  usage("hcspec_summary pha_file1,pha_file2,... file.ps [binfactor  [logT1,logT2,...,logTn [w1,w2,..,wn]]]");
}

autoload("hcspec_summary", "cspec_summary.sl");

fpha = strtok(__argv[1], ", ");       % array of input files.
if (length(fpha) == 1) fpha = fpha[0];

fout  = __argv[2];        % postscript output file name

if (NARGS > 3)
  bfact = eval(__argv[3]);  % integer binning factor
if (NARGS > 4)
{
  temps = eval(sprintf("[%s]", __argv[4]));    % temperature array;
  if (length(temps) == 1) temps = temps[0];
}
if (NARGS > 5)
{
  wgts =  eval(sprintf("[%s]", __argv[5]));    % weights array;
  if (length(wgts) == 1) temps = temps[0];
}

if (__is_initialized(&temps))  
  plasma(aped);           % load the database only if needed................

switch (NARGS)
{
  case 3:
  hcspec_summary(fpha, fout);
}  
{
  case 4:
  hcspec_summary(fpha, fout, bfact);
}  
{
  case 5:
  hcspec_summary(fpha, fout, bfact, temps, temps-temps+1);
}  
{
  case 6:
  hcspec_summary(fpha, fout, bfact, temps, wgts);
}