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

mplot.figure(0)
SlopeArray_40_1 = Return_Mean_Slope('H1RG-022', 'LEACH', \
                 '/nfs/slac/g/ki/ki03/lances/H1RG-022/ForwardBiasedDiodes/ForwardBiasedFollowedBy40Reads1_08Jun28_Slopes.lst',\
                 XStart=XStart, XStop=XStop, YStart=YStart, YStop=YStop)
mplot.clf()
mplot.semilogy(SlopeArray_40_1[:,4], SlopeArray_40_1[:,MeanModeOrMedian])

SlopeArray_40_2 = Return_Mean_Slope('H1RG-022', 'LEACH', \
                 '/nfs/slac/g/ki/ki03/lances/H1RG-022/ForwardBiasedDiodes/ForwardBiasedFollowedBy40Reads2_08Jun28_Slopes.lst',\
                 XStart=XStart, XStop=XStop, YStart=YStart, YStop=YStop)
mplot.semilogy(SlopeArray_40_2[:,4], SlopeArray_40_2[:,MeanModeOrMedian])

SlopeArray_50_1 = Return_Mean_Slope('H1RG-022', 'LEACH', \
                 '/nfs/slac/g/ki/ki03/lances/H1RG-022/ForwardBiasedDiodes/ForwardBiasedFollowedBy50Reads1_08Jun28_Slopes.lst',\
                 XStart=XStart, XStop=XStop, YStart=YStart, YStop=YStop)
mplot.semilogy(SlopeArray_50_1[:,4], SlopeArray_50_1[:,MeanModeOrMedian])

SlopeArray_50_2 = Return_Mean_Slope('H1RG-022', 'LEACH', \
                 '/nfs/slac/g/ki/ki03/lances/H1RG-022/ForwardBiasedDiodes/ForwardBiasedFollowedBy50Reads2_08Jun28_Slopes.lst',\
                 XStart=XStart, XStop=XStop, YStart=YStart, YStop=YStop)
mplot.semilogy(SlopeArray_50_2[:,4], SlopeArray_50_2[:,MeanModeOrMedian])

SlopeArray_10_2 = Return_Mean_Slope('H1RG-022', 'LEACH', \
                 '/nfs/slac/g/ki/ki03/lances/H1RG-022/ForwardBiasedDiodes/ForwardBiasedFollowedBy10Reads1_08Jun27_Slopes.lst',\
                 XStart=XStart, XStop=XStop, YStart=YStart, YStop=YStop)
mplot.semilogy(SlopeArray_10_2[:,4], SlopeArray_10_2[:,MeanModeOrMedian])

mplot.legend()
mplot.xlabel(r'\textbf{Seconds After Forward Bias $\rightarrow$ Reverse Bias}')
mplot.ylabel(r'\textbf{Dark Current\ \ (e$^{-}$/sec/pix})')
mplot.figtext(.21,.730, r'\textbf{\ \ Detector: } H1RG-022')
mplot.figtext(.23,.790, (r'\textbf{Region: } [' +
                         str(int(XStart)) + ':' + str(int(XStop)) + ',' +
                         str(int(YStart)) + ':' + str(int(YStop)) + ']'))
mplot.savefig(PlotDir+'DarkCurrentAfterForwardBias_All_Reads.eps')
mplot.savefig(PlotDir+'DarkCurrentAfterForwardBias_All_Reads.png')