#!/bin/bash

echo "Enter the Start Week: "
read week_beg
echo "Enter the End Week: "
read week_end

week=$week_beg

while [ $week -le $week_end ]
do
echo "=============> Filtering data of week" $week

wstr="w"
if [ $week -le 99 ]
then
    wstr="w0"
fi

echo level1
gtselect evclass=3 \
         infile=lat_photon_weekly_"$wstr""$week"_p202_v001.fits \
         outfile=lat_"$wstr""$week"_l1_v1.fits \
         ra=180 dec=0 rad=180 tmin=INDEF tmax=INDEF \
         emin=10000 emax=20000 zmax=100
gtmktime scfile=lat_spacecraft_weekly_"$wstr""$week"_p202_v001.fits \
	   filter="DATA_QUAL==1 && LAT_CONFIG==1 && ABS(ROCK_ANGLE)<52" \
	   roicut=no evfile=lat_"$wstr""$week"_l1_v1.fits \
	   outfile=lat_"$wstr""$week"_l1.fits

echo level2
gtselect evclass=3 \
         infile=lat_photon_weekly_"$wstr""$week"_p202_v001.fits \
         outfile=lat_"$wstr""$week"_l2_v1.fits \
         ra=180 dec=0 rad=180 tmin=INDEF tmax=INDEF \
         emin=20000 emax=30000 zmax=100
gtmktime scfile=lat_spacecraft_weekly_"$wstr""$week"_p202_v001.fits \
	   filter="DATA_QUAL==1 && LAT_CONFIG==1 && ABS(ROCK_ANGLE)<52" \
	   roicut=no evfile=lat_"$wstr""$week"_l2_v1.fits \
	   outfile=lat_"$wstr""$week"_l2.fits
echo level3
gtselect evclass=3 \
         infile=lat_photon_weekly_"$wstr""$week"_p202_v001.fits \
         outfile=lat_"$wstr""$week"_l3_v1.fits \
         ra=180 dec=0 rad=180 tmin=INDEF tmax=INDEF \
         emin=30000 emax=40000 zmax=100
gtmktime scfile=lat_spacecraft_weekly_"$wstr""$week"_p202_v001.fits \
	   filter="DATA_QUAL==1 && LAT_CONFIG==1 && ABS(ROCK_ANGLE)<52" \
	   roicut=no evfile=lat_"$wstr""$week"_l3_v1.fits \
	   outfile=lat_"$wstr""$week"_l3.fits

echo level4
gtselect evclass=3 \
         infile=lat_photon_weekly_"$wstr""$week"_p202_v001.fits \
         outfile=lat_"$wstr""$week"_l4_v1.fits \
         ra=180 dec=0 rad=180 tmin=INDEF tmax=INDEF \
         emin=40000 emax=50000 zmax=100
gtmktime scfile=lat_spacecraft_weekly_"$wstr""$week"_p202_v001.fits \
	   filter="DATA_QUAL==1 && LAT_CONFIG==1 && ABS(ROCK_ANGLE)<52" \
	   roicut=no evfile=lat_"$wstr""$week"_l4_v1.fits \
	   outfile=lat_"$wstr""$week"_l4.fits

echo level5
gtselect evclass=3 \
         infile=lat_photon_weekly_"$wstr""$week"_p202_v001.fits \
         outfile=lat_"$wstr""$week"_l5_v1.fits \
         ra=180 dec=0 rad=180 tmin=INDEF tmax=INDEF \
         emin=50000 emax=60000 zmax=100
gtmktime scfile=lat_spacecraft_weekly_"$wstr""$week"_p202_v001.fits \
	   filter="DATA_QUAL==1 && LAT_CONFIG==1 && ABS(ROCK_ANGLE)<52" \
	   roicut=no evfile=lat_"$wstr""$week"_l5_v1.fits \
	   outfile=lat_"$wstr""$week"_l5.fits

(( week++ ))

done