#WindowStarPhotometryStats.py
#
#PURPOSE:
#  To extract the quantities from the WindowStarPhotometryStars database and analyze them 
#  quantitatively
#
execfile('/afs/slac/u/ki/lances/python/python_HxRG/HxRG_Setup.py')
execfile('/afs/slac/u/ki/lances/python/ThesisPlotSettings.py')

import WindowStarPhotometry_Extract_Database
reload(WindowStarPhotometry_Extract_Database)
from WindowStarPhotometry_Extract_Database import *

PixTimeStr = '10uS'
TimeStr    = '03_28_11'
FrameNums_10_1, WindowNums_10_1, XCs_10_1, YCs_10_1, XCMs_10_1, YCMs_10_1, \
   Sharps_10_1, Rounds_10_1, Residuals_10_1, Mags_10_1, FindFluxes_10_1, AperFluxes_10_1, \
   Peaks_10_1, G1s_10_1, G2s_10_1, SigXs_10_1, SigYs_10_1, FWHMXs_10_1, FWHMYs_10_1, \
   DTKs_10_1, FracSecs_10_1=\
WindowStarPhotometry_Extract_Database('H1RG-022','ASIC','SAO59392', QueryString=' WHERE FILENAME REGEXP '+"'"+PixTimeStr+"'"+ " AND FILENAME REGEXP "+"'"+TimeStr+"'" + " AND XC !=-1")

PixTimeStr = '10uS'
TimeStr    = '03_29_12'
FrameNums_10_2, WindowNums_10_2, XCs_10_2, YCs_10_2, XCMs_10_2, YCMs_10_2, \
   Sharps_10_2, Rounds_10_2, Residuals_10_2, Mags_10_2, FindFluxes_10_2, AperFluxes_10_2, \
   Peaks_10_2, G1s_10_2, G2s_10_2, SigXs_10_2, SigYs_10_2, FWHMXs_10_2, FWHMYs_10_2, \
   DTKs_10_2, FracSecs_10_2=\
WindowStarPhotometry_Extract_Database('H1RG-022','ASIC','SAO59392', QueryString=' WHERE FILENAME REGEXP '+"'"+PixTimeStr+"'"+ " AND FILENAME REGEXP "+"'"+TimeStr+"'" + " AND XC !=-1")

PixTimeStr = '8uS'
TimeStr    = '03_30_37'
FrameNums_8_1, WindowNums_8_1, XCs_8_1, YCs_8_1, XCMs_8_1, YCMs_8_1, \
   Sharps_8_1, Rounds_8_1, Residuals_8_1, Mags_8_1, FindFluxes_8_1, AperFluxes_8_1, \
   Peaks_8_1, G1s_8_1, G2s_8_1, SigXs_8_1, SigYs_8_1, FWHMXs_8_1, FWHMYs_8_1, \
   DTKs_8_1, FracSecs_8_1=\
WindowStarPhotometry_Extract_Database('H1RG-022','ASIC','SAO59392', QueryString=' WHERE FILENAME REGEXP '+"'"+PixTimeStr+"'"+ " AND FILENAME REGEXP "+"'"+TimeStr+"'" + " AND XC !=-1")

PixTimeStr = '8uS'
TimeStr    = '03_33_02'
FrameNums_8_2, WindowNums_8_2, XCs_8_2, YCs_8_2, XCMs_8_2, YCMs_8_2, \
   Sharps_8_2, Rounds_8_2, Residuals_8_2, Mags_8_2, FindFluxes_8_2, AperFluxes_8_2, \
   Peaks_8_2, G1s_8_2, G2s_8_2, SigXs_8_2, SigYs_8_2, FWHMXs_8_2, FWHMYs_8_2, \
   DTKs_8_2, FracSecs_8_2=\
WindowStarPhotometry_Extract_Database('H1RG-022','ASIC','SAO59392', QueryString=' WHERE FILENAME REGEXP '+"'"+PixTimeStr+"'"+ " AND FILENAME REGEXP "+"'"+TimeStr+"'" + " AND XC !=-1")

PixTimeStr = '8uS'
TimeStr    = '03_34_04'
FrameNums_8_3, WindowNums_8_3, XCs_8_3, YCs_8_3, XCMs_8_3, YCMs_8_3, \
   Sharps_8_3, Rounds_8_3, Residuals_8_3, Mags_8_3, FindFluxes_8_3, AperFluxes_8_3, \
   Peaks_8_3, G1s_8_3, G2s_8_3, SigXs_8_3, SigYs_8_3, FWHMXs_8_3, FWHMYs_8_3, \
   DTKs_8_3, FracSecs_8_3=\
WindowStarPhotometry_Extract_Database('H1RG-022','ASIC','SAO59392', QueryString=' WHERE FILENAME REGEXP '+"'"+PixTimeStr+"'"+ " AND FILENAME REGEXP "+"'"+TimeStr+"'" + " AND XC !=-1")

PixTimeStr = '5.92uS'
TimeStr    = '03_36_31'
FrameNums_6_1, WindowNums_6_1, XCs_6_1, YCs_6_1, XCMs_6_1, YCMs_6_1, \
   Sharps_6_1, Rounds_6_1, Residuals_6_1, Mags_6_1, FindFluxes_6_1, AperFluxes_6_1, \
   Peaks_6_1, G1s_6_1, G2s_6_1, SigXs_6_1, SigYs_6_1, FWHMXs_6_1, FWHMYs_6_1, \
   DTKs_6_1, FracSecs_6_1=\
WindowStarPhotometry_Extract_Database('H1RG-022','ASIC','SAO59392', QueryString=' WHERE FILENAME REGEXP '+"'"+PixTimeStr+"'"+ " AND FILENAME REGEXP "+"'"+TimeStr+"'" + " AND XC !=-1")

PixTimeStr = '5.92uS'
TimeStr    = '03_35_46'
FrameNums_6_2, WindowNums_6_2, XCs_6_2, YCs_6_2, XCMs_6_2, YCMs_6_2, \
   Sharps_6_2, Rounds_6_2, Residuals_6_2, Mags_6_2, FindFluxes_6_2, AperFluxes_6_2, \
   Peaks_6_2, G1s_6_2, G2s_6_2, SigXs_6_2, SigYs_6_2, FWHMXs_6_2, FWHMYs_6_2, \
   DTKs_6_2, FracSecs_6_2=\
WindowStarPhotometry_Extract_Database('H1RG-022','ASIC','SAO59392', QueryString=' WHERE FILENAME REGEXP '+"'"+PixTimeStr+"'"+ " AND FILENAME REGEXP "+"'"+TimeStr+"'" + " AND XC !=-1")

mplot.figure(0)
mplot.clf()
mplot.plot(FrameNums_10_1, FindFluxes_10_1, 'ro')
mplot.plot(FrameNums_10_2, FindFluxes_10_2, 'rv')
mplot.plot(FrameNums_8_1, FindFluxes_8_1, 'go')
mplot.plot(FrameNums_8_1, FindFluxes_8_1, 'gv') 
mplot.plot(FrameNums_6_1, FindFluxes_6_1, 'bv')
mplot.plot(FrameNums_6_2, FindFluxes_6_2, 'bo')

mplot.figure(1)
mplot.clf()
mplot.plot(FrameNums_10_1, Peaks_10_1, 'ro')
mplot.plot(FrameNums_10_2, Peaks_10_2, 'rv')
mplot.plot(FrameNums_8_1, Peaks_8_1, 'go')
mplot.plot(FrameNums_8_1, Peaks_8_1, 'gv')
mplot.plot(FrameNums_6_1, Peaks_6_1, 'bv')
mplot.plot(FrameNums_6_2, Peaks_6_2, 'bo')