Use/test Dave's stuff to Get and Reprocess HETG data

Re-do E0102 data to have explicit order sorting limits
of 0.20 for '3828 and
0.25 for '120 and '968 because no tgain/cti (-110 deg).

080620 Re-do extraction centered 150 pixels in CHIPY from nominal
       to get background events, obs_3828bkg.
       Offset to use is:  
           Delta SkyX = DeltaChipY in pixels *  -1.0*sin(roll_nom)
	   Delta SkyY = DeltaChipY in pixels *  -1.0*cos(roll_nom)
- - -

0) Use Dave's sw that is in:  ~dph/h3/Analysis/TGCat/
unix>  setenv TGCat ~dph/h3/Analysis/TGCat/

0a) Make similar sub dirs here (.../Data/) and
copy whatever I use from dph to here...
unix>  mkdir bin scripts

[unix] cp $TGCat/scripts/*sl scripts/

Ia) Download:
cp $TGCat/bin/download_obsid bin/
  and also need:
cp ~dph/bin/chandra-get ~/bin/

try an E0102 HETG obs:
dd@caribou[86] bin/download_obsid 968


Ib) Configure:
copy the routine:
cp $TGCat/bin/setup_obsdir bin/
and execute it on the dir:
./bin/setup_obsdir obs_3828bkg

Ic) Reprocess:

Need to have CIAO[X] setup:

[unix] source /nfs/cxc/a1/setup/ciaox-setup 


unix> cd obs_968               % E0102 test case: wider extraction...
unix> isis
isis>       

append_to_isis_load_path("../scripts");
require("tg_repro_fun");

set_exec( 2 );
setup_pfiles;
obs_info = read_config( "evt0" );
make_obspar;
acis_process_events_te;

acis_evt_filter_ds;      %  <--- destreak used

tgdetect_single_src;
variable x,y,c;
(x,y,c) = read_src1a_pos;
x ;
y ;
c ;

% Set the center to:  approx from evt2img on evt1
% Use:   evt2img evt1 4000:4400:0.5,3800:4200:0.5 0.7,2.5
% Use:   evt2img evt1 3800:4200:0.5,3800:4200:0.5 0.7,2.5
% Use:   evt2img evt1 3900:4300:0.5,4000:4400:0.5 0.7,2.5
%
% obs_968 :
x = 4077.8;  % obsid: 120, 968, 3828: 4231.6, 4077.8, 4094.1
y = 4077.0;  %                        3964.6, 4077.0, 4166.1
% obs_120 :
x = 4231.6;  % obsid: 120, 968, 3828: 4231.6, 4077.8, 4094.1
y = 3964.6;  %                        3964.6, 4077.0, 4166.1

% obs_3828 :
x = 4094.1;  % obsid: 120, 968, 3828: 4231.6, 4077.8, 4094.1
y = 4166.1;  %                        3964.6, 4077.0, 4166.1
% obs_3828bkg :
x = 4094.1 + 137.0 ;  % obsid: 120, 968, 3828: 4231.6, 4077.8, 4094.1,  4231.1
y = 4166.1 - 61.0  ;  %                        3964.6, 4077.0, 4166.1,  4105.1

%-----------
%%%tg_create_mask;
%# For E0102 we don't use any background - so set it
%# wide enough for the whole remnant:
%# Zero radius of 60 pixels.
%# Try  +/- 30 arc seconds or 120 pixels total width.
%pset $p infile=$inf outfile=$ouf input_pos_tab=$ipos \
%     use_user_pars=$user sA_zero_x=$x0 sA_zero_y=$y0 \
%     sA_zero_rad=60 sA_width_heg=120 sA_width_meg=120 \
%     verbose=$vrb clobber=$clb
%
% instead:
rzo = 60.0;
wheg=120.0;
wmeg=120.0;
tg_create_mask(x,y,rzo,wheg,wmeg);
%-----------
% Fix the order sorting level explicitly
%  0.25 for '968 and '120;  0.20 for '3828
nps = Assoc_Type[ Any_Type ];
nps["osipfile"]="NONE";
nps["osort_lo"]=0.20;
nps["osort_hi"]=0.20;
tg_resolve_events_te(nps);
%------------
%%%tgextract;
% For E0102 extract we want to set a bunch of non-standard things a la:
%# Set extraction region:
%#   Note: +/-0.0015 deg ~ +/-9.42 pixels
%# Want +/- 30" ~ +/- 60 pix --> +/- 0.0096 deg.
%pset $p infile=$inf outfile=$ouf tg_order_list=${ords} \
%        extract_background=${exbkg} \
%        min_tg_d=-0.0096 max_tg_d=0.0096 \
%##      max_upbkg_tg_d=0.0045 min_downbkg_tg_d=-0.0045 \
%        verbose=$vrb clobber=$clb
%
% instead:
np_tgex = Assoc_Type[ Any_Type ];
np_tgex["extract_background"]= "no";
np_tgex["min_tg_d"]=-0.0096;
np_tgex["max_tg_d"]= 0.0096;
np_tgex["tg_order_list"]="-1,1";
tgextract(np_tgex);
%----------
orders=[-1,1];
make_responses(orders);
x;
y;
system("ls");



%%%  Make zeroth-order arf

***   to reprocess:
unix> cd obs_3828              % E0102 test case: wider extraction...
unix> isis
isis>       
      
require("../scripts/tg_repro_fun");
set_exec( 2 );
setup_pfiles;
obs_info = read_config( "evt0" );
make_obspar;
   exit;

%  obs_968:
x = 4077.8;  % obsid: 120, 968, 3828: 4231.6, 4077.8, 4094.1
y = 4077.0;  %                        3964.6, 4077.0, 4166.1

mkarf \
  asphistfile='./s3.asphist' \
  outfile='hetg_0.arf' \
  sourcepixelx=4077.8  \
  sourcepixely=4077.0  \
  engrid='grid(../aciss_aimpt_cy09.rmf[MATRIX][cols ENERG_LO,ENERG_HI])' \
  obsfile='./evt2' \
  pbkfile='./pbk0' \
  detsubsys='ACIS-S3' \
  grating='HETG' \
  maskfile='./msk1' \
  verbose=1 \
  clobber='yes' \
  mode='hl'

% obs_120 :
x = 4231.6;  % obsid: 120, 968, 3828: 4231.6, 4077.8, 4094.1
y = 3964.6;  %                        3964.6, 4077.0, 4166.1

mkarf \
  asphistfile='./s3.asphist' \
  outfile='hetg_0.arf' \
  sourcepixelx=4231.6  \
  sourcepixely=3964.6  \
  engrid='grid(../aciss_aimpt_cy09.rmf[MATRIX][cols ENERG_LO,ENERG_HI])' \
  obsfile='./evt2' \
  pbkfile='./pbk0' \
  detsubsys='ACIS-S3' \
  grating='HETG' \
  maskfile='./msk1' \
  verbose=1 \
  clobber='yes' \
  mode='hl'
  
%  obs_3828
x = 4094.1;  % obsid: 120, 968, 3828: 4231.6, 4077.8, 4094.1
y = 4166.1;  %                        3964.6, 4077.0, 4166.1
  
mkarf \
  asphistfile='./s3.asphist' \
  outfile='hetg_0.arf' \
  sourcepixelx=4094.1  \
  sourcepixely=4166.1  \
  engrid='grid(../aciss_aimpt_cy09.rmf[MATRIX][cols ENERG_LO,ENERG_HI])' \
  obsfile='./evt2' \
  pbkfile='./pbk0' \
  detsubsys='ACIS-S3' \
  grating='HETG' \
  maskfile='./msk1' \
  verbose=1 \
  clobber='yes' \
  mode='hl'