% file: steps_e0102_bkg.txt
%
% in:  $SCI/E0102/Data
%
% Processing E0102 the 4 E0102 obsids for Cal etc analyses.
%   968, 120, 3828, 12147

%  *** Background extractions ***  do after the steps_e0102_src.



% - - -

% Include the 120 pixel offset amounts in these values:

%    120    4228.4 ,  3963.6              4227.9-24.    3970.6-117.
%
%    968    4112.5 ,  4105.5              4112.0-24.    4112.5-117.
%
%   3828    4093.75 , 4165.75             4093.25+109.   4172.75-49.
%
%  12147    4076.0 ,  4132.50    -->      4075.5+100.    4139.5+66.
%
% - - -

% - - -
% Use ciao:
% CIAO 4.4 Friday, December  2, 2011
%   bindir      : /nfs/cxc/a1/linux-x86_64/opt/packages/ciao-4.4/bin
source /nfs/cxc/a1/setup/ciao-setup -o


% Processing will be in $SCI/E0102/Data:

% *** Set the Obsid to be working on: ***
setenv OBSDIR 968

cd $SCI/E0102/Data/

% start isis (with ~dd/.isisrc) :
isis

% use tgcat:
%   will use the version in ~dd/hy3d/Isis_Extras/
%   .sl files from tgcat/scripts copied there, version:
%                   %; Time-stamp: <2011-06-17 11:51:50 dph>
%                 "tgcat.sl version 1.10.2 loaded."
require("tgcat");


% Setup custom parameters as in previous analyses:

% --- fixed zo location:                *****  copy from the BACKGROUND values above ***
src_info = set_source_detection_info("none", "pixel", 4112.0-24.,  4112.5-117.);

% --- mask:  55 pixels for the whole SNR
mask_info = set_tg_create_mask_info( 55.0, 110.0, 110.0 );

% --- extraction:
% Want +/- 55 pix
extract_info = set_tgextract_info( -55.0*0.492, 55.0*0.492; arcsec );

% --- fixed osort limits:
% * pick one:
% 968 and 120:
%%%%%   osort_info = set_tgre_osort_info( 0.25, 0.25);
% 3828 and 12147:
%%%%%   osort_info = set_tgre_osort_info( 0.20, 0.20);


% Do it:
run_pipe("./obs_"+getenv("OBSDIR")+""; detect_info=src_info, mask_info=mask_info,
                           extract_info=extract_info, osort_info=osort_info);

% check output a bit:
system("ls obs_"+getenv("OBSDIR"));
system("gv obs_"+getenv("OBSDIR")+"/summary_im-a.ps");
system("gv obs_"+getenv("OBSDIR")+"/summary_imsp.ps");
system("gv obs_"+getenv("OBSDIR")+"/summary_im_osum.ps");

quit;


% Make a new dir to save the products to:
mkdir obs_{$OBSDIR}_bkg
cp obs_{$OBSDIR}/evt2 obs_{$OBSDIR}_bkg/
cp obs_{$OBSDIR}/pha2 obs_{$OBSDIR}_bkg/
cp obs_{$OBSDIR}/*.arf obs_{$OBSDIR}_bkg/
cp obs_{$OBSDIR}/sum*.ps obs_{$OBSDIR}_bkg/

% reduce the evt2 size, filter 0.3 to 8 keV:
cd obs_{$OBSDIR}_bkg
source /nfs/cxc/a1/setup/lheasoft-setup
fselect evt2 evt2_Elt8
   Selection Expression[] ENERGY .GT. 300 .AND. ENERGY .LT. 8000
ls -lst *evt*

rm evt2

gzip evt2_Elt8
cd ..

% make a tar file to move products to Mac, etc:
tar -cvf obs_{$OBSDIR}_bkg_Files.tar obs_{$OBSDIR}_bkg/*
gzip obs_{$OBSDIR}_bkg_Files.tar

% copy to Mac and enjoy ;-)
%    (unzip the evt2_Elt8 file for faster reading in)



% - - -