##PlotWellPctVsMax.py
##
##PURPOSE:
##
##To plot the decay of the persistent cores as a function of time.
##A fit of the form 
##
##f=A*exp(-t/B)+C*exp(-t/D)
##
##is applied
##
##PLOTTING SYMBOLS
##5 Reads  - Diamonds
##10 Reads - Squares
##15 Reads - Triangles
##20 Reads - Circles
##25 Reads - Pentagons
##30 Reads - Hexagons
##35 Reads - Down Triangles
PSym5='D'
PSym10='s'
PSym15='^'
PSym20='o'
PSym25='p'
PSym30='h'
PSym35='v'
PSym40='<'
PSym45='>'
PSym50='d'
PSym55='H'
PSym60='v'
PSym65='>'
PSym70='<'
PSym75='^'
##Y Band - Dark Red   #660000
##I Band - Light Red
##G Band - Blue
YCol = '#990000' 
ICol = '#009900' 
GCol = '#000099' 

execfile('/afs/slac/u/ki/lances/python/python_HxRG/HxRG_Setup.py')
execfile('/afs/slac/u/ki/lances/python/SPIEPlotSettings.py')
from Return_Persistence_DbQuantities import *
from MeanAndVarianceByBin import *
from HxRG_Class import *

import pdb

Mode=2
QueryString = " WHERE READNUM=9 AND FILENUM=1 "
PixID= 2437693

#Fitting Function for decay
FitFunc = lambda p, x: p[0]*exp(-x/p[1])+p[2]*exp(-x/p[3])
ErrFunc = lambda p, x, y: FitFunc(p,x)-y

##NGC2395:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
##5 Reads I Band
ParList = '/u/ki/lances/python/python_HxRG/PersistenceExpI/NGC956_35Reads_I.py'
execfile(ParList)

##Get some info
HxRG=HxRG_C(FitsFileName=SlopeFiles[0])
PlotDir = HxRG.DataDir+HxRG.DetStr+'/Plots/NegPersistence/'

##READ THE FILES IN FROM FILELIST
FitsFiles=ReadFileList(FileLists[0])
NumFitsFiles = len(FitsFiles)

#Get the values from the slope image
SHxRG=HxRG_C(FitsFileName=FitsFiles[1])
SHxRG.Get_Raw_Header(FitsFiles[1])

Minp1s, Maxp1s = Return_Persistence_DbQuantities(\
                 DetStr, ElecStr, Dates, NReads, NDarkReads,\
                 DataBaseName, DarkTableNames, SlopeTableNames,\
                 FileLists, NTables, QueryString, Mode, PixID=PixID, \
                 HxRG=SHxRG, TvFlag=1, Verbosity=1)