execfile('/afs/slac/u/ki/lances/python/python_HxRG/HxRG_Setup.py')
execfile('/afs/slac/u/ki/lances/python/ThesisPlotSettings.py')
PlotDir = '/nfs/slac/g/ki/ki03/lances/KPNO/Latex/Thesis/ForwardAndChangeBias/Figures/'
from Return_Mean_Slope import *
MeanModeOrMedian = 2
XStart=200
XStop=500
YStart=200
YStop=500

SlopeArray_0_2 = Return_Mean_Slope('H1RG-018', 'LEACH', \
                 '/nfs/slac/g/ki/ki03/lances/H1RG-018/VSUBChanges/VSUB_0_to_2Volts.lst',\
                 XStart=XStart, XStop=XStop, YStart=YStart, YStop=YStop)
mplot.clf()
mplot.semilogy(SlopeArray_0_2[:,4], SlopeArray_0_2[:,MeanModeOrMedian], \
               label = r'\textbf{0 to 2 Volts}')

SlopeArray_2_4 = Return_Mean_Slope('H1RG-018', 'LEACH', \
                 '/nfs/slac/g/ki/ki03/lances/H1RG-018/VSUBChanges/VSUB_2_to_4Volts.lst',\
                  XStart=XStart, XStop=XStop, YStart=YStart, YStop=YStop)
mplot.semilogy(SlopeArray_2_4[:,4], SlopeArray_2_4[:,MeanModeOrMedian],\
               label = r'\textbf{2 to 4 Volts}')

SlopeArray_4_6 = Return_Mean_Slope('H1RG-018', 'LEACH', \
                 '/nfs/slac/g/ki/ki03/lances/H1RG-018/VSUBChanges/VSUB_4_to_6Volts.lst',\
                  XStart=XStart, XStop=XStop, YStart=YStart, YStop=YStop)
mplot.semilogy(SlopeArray_4_6[:,4], SlopeArray_4_6[:,MeanModeOrMedian],\
               label = r'\textbf{4 to 6 Volts}')

SlopeArray_6_8 = Return_Mean_Slope('H1RG-018', 'LEACH', \
                 '/nfs/slac/g/ki/ki03/lances/H1RG-018/VSUBChanges/VSUB_6_to_8Volts.lst',\
                  XStart=XStart, XStop=XStop, YStart=YStart, YStop=YStop)
mplot.semilogy(SlopeArray_6_8[:,4], SlopeArray_6_8[:,MeanModeOrMedian],\
               label = r'\textbf{6 to 8 Volts}')

SlopeArray_8_10 = Return_Mean_Slope('H1RG-018', 'LEACH', \
                 '/nfs/slac/g/ki/ki03/lances/H1RG-018/VSUBChanges/VSUB_8_to_10Volts.lst',\
                  XStart=XStart, XStop=XStop, YStart=YStart, YStop=YStop)
mplot.semilogy(SlopeArray_8_10[:,4], SlopeArray_8_10[:,MeanModeOrMedian],\
               label = r'\textbf{8 to 10 Volts}')
mplot.legend()
mplot.xlabel(r'\textbf{Seconds After Increasing V$_{SUB}$}')
mplot.ylabel(r'\textbf{Dark Current\ \ (e$^{-}$/sec/pix})')
mplot.figtext(.21,.730, r'\textbf{\ \ Detector: } H1RG-018')
mplot.figtext(.23,.790, (r'\textbf{Region: } [' +
                         str(int(XStart)) + ':' + str(int(XStop)) + ',' +
                         str(int(YStart)) + ':' + str(int(YStop)) + ']'))
mplot.savefig(PlotDir+'DarkCurrentVsChangeInVSUB.eps')
mplot.savefig(PlotDir+'DarkCurrentVsChangeInVSUB.png')