#!/bin/sh

pd=/nfs/wiwaxia/dc/hetgs/obsid_01103-2/primary/
sd=/nfs/wiwaxia/dc/hetgs/obsid_01103-2/secondary/

irn=acisf01103_000N002
orn=acisf01103

flt="grade=0,2,3,4,6,status=000000000000xxxx0000000000000000,energy<12000"
dmcopy ${pd}/${irn}_evt1.fits.gz"[events][${flt}]" Evt/${orn}-1_evt1.fits opt=all
ln -s ${sd}/${irn}_src1a.fits Reg
ln -s ${sd}/aspect/pcad*_asol1.fits* Asp
ln -s ${sd}/${irn}_bpix1.fits Reg

fon=./Evt/${orn}_evt1.fits

punlearn acis_process_events
pset acis_process_events mode=hl verbose=0
pset acis_process_events infile=./Evt/${orn}-1_evt1.fits outfile=${fon}
pset acis_process_events acaofffile='./Asp/pcad*_asol1*'
pset acis_process_events badpixfile=./Reg/${irn}_bpix1.fits


punlearn tgdetect
pset tgdetect mode=hl verbose=0
pset tgdetect infile=${fon} outfile=./Reg/${orn}_src1a.fits


punlearn tg_create_mask
pset tg_create_mask mode=hl verbose=0 
pset tg_create_mask input_pos_tab=./Reg/${orn}_src1a.fits 
pset tg_create_mask infile=${fon} outfile=./Reg/${orn}_reg1a.fits


punlearn tg_resolve_events
pset tg_resolve_events mode=hl verbose=0
pset tg_resolve_events infile=${fon} outfile=./Evt/${orn}_evt1a.fits
pset tg_resolve_events regionfile=./Reg/${orn}_reg1a.fits
pset tg_resolve_events acaofffile='./Asp/pcad*_asol1*'


punlearn tgextgract
pset tgextract mode=hl verbose=0
pset tgextract infile=./Evt/${orn}_evt1a.fits outfile=./Pha/${orn}_pha2.fits
pset tgextract tg_order_list="-1,1"


punlearn geom
pset geom instruments=CALDB

acis_process_events
tgdetect
tg_create_mask
tg_resolve_events
tgextract


# Hey! That was fun! Let's do it again!  aHiullk...

fon=./Evt/${orn}-2_evt1.fits
pset acis_process_events infile=./Evt/${orn}-1_evt1.fits outfile=${fon}


pset tgdetect infile=${fon} outfile=./Reg/${orn}-2_src1a.fits


pset tg_create_mask input_pos_tab=./Reg/${orn}-2_src1a.fits 
pset tg_create_mask infile=${fon} outfile=./Reg/${orn}-2_reg1a.fits


pset tg_resolve_events infile=${fon} outfile=./Evt/${orn}-2_evt1a.fits
pset tg_resolve_events regionfile=./Reg/${orn}-2_reg1a.fits


pset tgextract infile=./Evt/${orn}-2_evt1a.fits outfile=./Pha/${orn}-2_pha2.fits
pset tgextract tg_order_list="-1,1"

pset geom instruments=./Ard/test_geom3b.fits

acis_process_events
tgdetect
tg_create_mask
tg_resolve_events
tgextract