#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/H2RG-32-147/ASIC/Reduced/07Dec14/SAO116737/Astrometry/SAO116737_RefMap.txt'

if SingleFile == 1:
  #I Band
  #15 Reads
  FWHM         = 6.0
  HWHM         = FWHM/2
  DataMin      = 0.0
  Sky          = 2.0
  Thresh       = 40
  FitsFileName = '/nfs/slac/g/ki/ki04/lances/H2RG-32-147/ASIC/Reduced/07Dec14/SAO116737/SAO116737_H2RG_SIPIN_30_Reads_Dec15_2007_03_56_32_I_RPS_SlopeFlatFielded.fits[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)

  #I Band 30 Reads
  FWHM         = 6.2
  HWHM         = FWHM/2
  DataMin      = 0.0
  Sky          = 2.0
  Thresh       = 40
  FitsFileName = '/nfs/slac/g/ki/ki04/lances/H2RG-32-147/ASIC/Reduced/07Dec14/SAO116737/SAO116737_H2RG_SIPIN_30_Reads_Dec15_2007_03_58_18_I_RPS_SlopeFlatFielded.fits[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)

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)

FitsFileName = '/nfs/slac/g/ki/ki04/lances/H2RG-32-147/ASIC/Reduced/07Dec14/SAO116737/SAO116737_Dither5_20ArcSecondsWest_0ArcSecondsNorth_H2RG_SIPIN_30_Reads_Dec15_2007_02_54_09_I_RPS_SlopeFlatFielded.fits[0]'
FitsFileName = '/nfs/slac/g/ki/ki04/lances/H2RG-32-147/ASIC/Reduced/07Dec14/SAO116737/SAO116737_Dither0_0ArcSecondsEast_0ArcSecondsNorth_H2RG_SIPIN_30_Reads_Dec15_2007_02_30_47_I_RPS_SlopeFlatFielded.fits[0]'