#!/bin/sh

# 2001.06.29
# test new geom files.

irn=hrcf01248_000N002
orn=hrcf01248-1

tmpdir=/nfs/wiwaxia/dd/dphtmp/
fon=${tmpdir}${orn}_evt1.fits

punlearn hrc_process_events
pset hrc_process_events mode=hl verbose=0
pset hrc_process_events infile=./Evt/${irn}_evt1.fits outfile=${fon}
pset hrc_process_events acaofffile='./Asp/pcadf058542036N002_asol1.fits'
pset hrc_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=${tmpdir}${orn}_evt1a.fits
pset tg_resolve_events regionfile=./Reg/${orn}_reg1a.fits
pset tg_resolve_events acaofffile=./Asp/pcadf058542036N002_asol1.fits
pset tg_resolve_events rmfile="NONE" eventdef=")stdlev1_HRC"


punlearn tgextract
pset tgextract mode=hl verbose=0
#pset tgextract inregion_file=CALDB
pset tgextract inregion_file=/nfs/cxc/a1/solaris/opt/CIAO_2.1_solaris/CALDB/data/chandra/hrc/bcf/tgmask2/letgD1999-07-22regN0001.fits
pset tgextract infile=${tmpdir}${orn}_evt1a.fits outfile=./Pha/${orn}_pha2.fits
pset tgextract tg_order_list="-1,1"


punlearn geom
pset geom instruments=CALDB

hrc_process_events
tgdetect
tg_create_mask
tg_resolve_events
tgextract


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

orn=hrcf01248-2
fon=${tmpdir}${orn}_evt1.fits
pset hrc_process_events infile=./Evt/${irn}_evt1.fits outfile=${fon}


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


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


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


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

pset geom instruments=./Ard/telD1999-07-23geomN0004B.fits

hrc_process_events
tgdetect
tg_create_mask
tg_resolve_events
tgextract