#!/usr/bin/csh -f
# Usage:   cleanit eventlist [gti1 gti2]
#
# cleanit file.evt acism62520_000N001_e2_flt1.fits asp_glitches.gti

# Clean the data in the following ways:

# apply sets of GTI tables passed in 2nd parameter
# grade filtered g02346
# non-zero status events removed
# energy filtered BI: [200,8000]eV
#		  FI: [500,8000]eV
# various column, grade, energy combination removed
		  
echo " "

# Apply any GTI tables provided...
setenv gti1 $2
setenv gti2 $3

if (${gti1} != "") then
  apply_gti $1 ${gti1} clean_gti1.evt
else
  ln -s $1 clean_gti1.evt
endif

if (${gti2} != "") then
  apply_gti clean_gti1.evt ${gti2} clean_gti2.evt
else
  ln -s clean_gti1.evt clean_gti2.evt
endif


# Split the data by CCD and apply various filters...
dmcopy "clean_gti2.evt[EVENTS][status=0, grade=0,2,3,4,6, energy<8000]" temp.clean

dmcopy "temp.clean[EVENTS][ccd_id=7,        energy>200]" temp.ccd7

dmcopy "temp.clean[EVENTS][ccd_id=8,grade=0,energy>500]" temp.ccd8

dmcopy "temp.clean[EVENTS][ccd_id=0,1,2,3,  energy>500]" temp.ccd0123  

rm clean_gti1.evt clean_gti2.evt


fselect "temp.ccd7[EVENTS]" temp.ccd7.clean "(chipx .ne. 512||(fltgrade .ne. 2&&fltgrade .ne. 64))&&(chipx .ne. 513||(fltgrade .ne. 2&&fltgrade .ne. 64))" clobber=yes

fselect "temp.ccd8[EVENTS]" temp.ccd8.clean "(energy .gt. 1000&&(chipx .ne. 2&&chipx .ne. 412&&chipx .ne. 491&&chipx .ne. 957))||(energy .le. 1000&&(chipx .ne. 2&&chipx .ne. 171&&chipx .ne. 178&&chipx .ne. 180&&chipx .ne. 409&&chipx .ne. 412&&chipx .ne. 422&&chipx .ne. 487&&chipx .ne. 491&&chipx .ne. 508&&chipx .ne. 514&&chipx .ne. 541&&chipx .ne. 632&&chipx .ne. 645&&chipx .ne. 910&&chipx .ne. 952))" clobber=yes

fselect "temp.ccd0123[EVENTS]" temp.ccd0123.clean "((chipx .ne. 256||fltgrade .ne. 8)&&(chipx .ne. 257||fltgrade .ne. 16)&&(chipx .ne. 512||(fltgrade .ne. 8&&fltgrade .ne. 104))&&(chipx .ne. 513||(fltgrade .ne. 16&&fltgrade .ne. 208))&&(chipx .ne. 768||fltgrade .ne. 8)&&(chipx .ne. 769||fltgrade .ne. 16))&&(ccd_id .ne. 0||chipx .ne. 81)&&(ccd_id .ne. 6||(chipx .ne. 58&&chipx .ne. 585))&&(ccd_id .ne. 2||(chipx .ne. 559&&chipx .ne. 571))" clobber=yes

rm temp.ccd7 temp.ccd8 temp.ccd0123
echo "DONE"
exit