#Go through all of the NGC2395 Images and do astrometry
#
from Astrometry_StarFind import *
from Astrometry_CCXYMatch import *

SingleFile  = 0 #0 to do the guided images, 1 to go through all dithers in the dither file list
TvFlag      = 1
InterActive = 1
FitOrder    = 2
CoordList = '/nfs/slac/g/ki/ki04/lances/H1RG-022/LEACH/Reduced/2007Nov18/NGC2416/Astrometry/NGC2416_UCAC2_Catalog.txt'

if SingleFile == 1:
  #I Band
  #15 Reads
  FWHM         = 5.0
  HWHM         = FWHM/2
  DataMin      = 0.0
  Sky          = 2.0
  Thresh       = 20
  EHigh        = 0.4
  FitsFileName = '/nfs/slac/g/ki/ki04/lances/H1RG-022/LEACH/Reduced/2007Nov18/NGC2416/NGC2416_Ifilter_600_1_0_0_-20_VSUB15_0_I_RPS_SlopeFlatFielded.fits[0]'
 
  R=Astrometry_StarFind(FitsFileName, \
    TvFlag=TvFlag, HWHMPSF=HWHM, DataMin=DataMin, InterActive=InterActive, Thresh=Thresh, EHigh=EHigh)

  R=Astrometry_CCXYMatch(FitsFileName, CoordList, TvFlag=TvFlag, FitOrder=FitOrder, InterActive=InterActive)

  #I Band
  #15 Reads
  FWHM         = 5.0
  HWHM         = FWHM/2
  DataMin      = 0.0
  Sky          = 2.0
  Thresh       = 20
  EHigh        = 0.4
  FitsFileName = '/nfs/slac/g/ki/ki04/lances/H1RG-022/LEACH/Reduced/2007Nov18/NGC2416/NGC2416_Ifilter_600_1_0_0_0_VSUB15_0_I_RPS_SlopeFlatFielded.fits[0]'

  R=Astrometry_StarFind(FitsFileName, \
    TvFlag=TvFlag, HWHMPSF=HWHM, DataMin=DataMin, InterActive=InterActive, Thresh=Thresh, EHigh=EHigh)

  R=Astrometry_CCXYMatch(FitsFileName, CoordList, TvFlag=TvFlag, FitOrder=FitOrder, InterActive=InterActive)

  #I Band
  #15 Reads
  FWHM         = 5.0
  HWHM         = FWHM/2
  DataMin      = 0.0
  Sky          = 2.0
  Thresh       = 20
  EHigh        = 0.4
  FitsFileName= '/nfs/slac/g/ki/ki04/lances/H1RG-022/LEACH/Reduced/2007Nov18/NGC2416/NGC2416_Ifilter_600_1_0_20_0_VSUB15_0_I_RPS_SlopeFlatFielded.fits[0]'

  R=Astrometry_StarFind(FitsFileName, \
    TvFlag=TvFlag, HWHMPSF=HWHM, DataMin=DataMin, InterActive=InterActive, Thresh=Thresh, EHigh=EHigh)

  R=Astrometry_CCXYMatch(FitsFileName, CoordList, TvFlag=TvFlag, FitOrder=FitOrder, InterActive=InterActive)

  #I Band
  #15 Reads
  FWHM         = 5.0
  HWHM         = FWHM/2
  DataMin      = 0.0
  Sky          = 2.0
  Thresh       = 20
  EHigh        = 0.4
  FitsFileName= '/nfs/slac/g/ki/ki04/lances/H1RG-022/LEACH/Reduced/2007Nov18/NGC2416/NGC2416_Ifilter_600_1_0_20_-20_VSUB15_0_I_RPS_SlopeFlatFielded.fits[0]'

  R=Astrometry_StarFind(FitsFileName, \
    TvFlag=TvFlag, HWHMPSF=HWHM, DataMin=DataMin, InterActive=InterActive, Thresh=Thresh, EHigh=EHigh)

  R=Astrometry_CCXYMatch(FitsFileName, CoordList, TvFlag=TvFlag, FitOrder=FitOrder, InterActive=InterActive)
elif SingleFile == 0:

  TvFlag      = 0
  InterActive = 0
  FitOrder    = 2

  #First File of Dither Sequence
  FWHM         = 6.5
  HWHM         = FWHM/2
  DataMin      = 0.0
  Sky          = 2.0
  Thresh       = 40
  FitsFileList = '/nfs/slac/g/ki/ki04/lances/H2RG-32-147/ASIC/Reduced/07Dec14/SAO116737/SAO116737_I_Slopes.lst'
  ##Read in the Dithered Offset List
  FitsFiles, XOffs, YOffs = ReadDitherList(FitsFileList)

  for FitsFile in FitsFiles:
    FitsFileName = FitsFile+'[0]'

    R=Astrometry_StarFind(FitsFileName, \
      TvFlag=TvFlag, HWHMPSF=HWHM, DataMin=DataMin, InterActive=InterActive, Thresh=Thresh)

    R=Astrometry_CCXYMatch(FitsFileName, CoordList, TvFlag=TvFlag, FitOrder=FitOrder, InterActive=InterActive)