##PlotFE55wMultiChannelAvg
##
##PURPOSE:
##  To plot FE55 histograms for the various number of channels being averaged.
##  These files are contained in 
##        /nfs/slac/g/ki/ki03/lances/H1RG-022/ASIC/08Jun29/
##
##SETUP:
execfile('/afs/slac/u/ki/lances/python/python_HxRG/HxRG_Setup.py')
execfile('/afs/slac/u/ki/lances/python/ThesisPlotSettings.py')
from FE55_Extract_Database import *

Bins1x1_G1, N1x1_G1, Bins2x2_G1, N2x2_G1, Bins3x3_G1, N3x3_G1, \
BinsCross_G1, NCross_G1 = \
    FE55_Extract_Database('H1RG-022','ASIC','08Jun29','1000', \
         QueryString='WHERE ADCs=1 and ADU2x2 !=0',Gain=1, NumADCs=1, \
         WindowM=1, Mode=1, UseRef=1, TvFlag=0)

Bins1x1_G2, N1x1_G2, Bins2x2_G2, N2x2_G2, Bins3x3_G2, N3x3_G2, \
BinsCross_G2, NCross_G2 = \
    FE55_Extract_Database('H1RG-022','ASIC','08Jun29','1000', \
         QueryString='WHERE ADCs=2 and ADU2x2 !=0',Gain=1, NumADCs=1, \
         WindowM=1, Mode=1, UseRef=1, TvFlag=0)

Bins1x1_G4, N1x1_G4, Bins2x2_G4, N2x2_G4, Bins3x3_G4, N3x3_G4, \
BinsCross_G4, NCross_G4 = \
    FE55_Extract_Database('H1RG-022','ASIC','08Jun29','1000', \
         QueryString='WHERE ADCs=4 and ADU2x2 !=0',Gain=1, NumADCs=1, \
         WindowM=1, Mode=1, UseRef=1, TvFlag=0)

Bins1x1_G6, N1x1_G6, Bins2x2_G6, N2x2_G6, Bins3x3_G6, N3x3_G6, \
BinsCross_G6, NCross_G6 = \
    FE55_Extract_Database('H1RG-022','ASIC','08Jun29','1000', \
         QueryString='WHERE ADCs=6 and ADU2x2 !=0',Gain=1, NumADCs=1, \
         WindowM=1, Mode=1, UseRef=1, TvFlag=0)

Bins1x1_G8, N1x1_G8, Bins2x2_G8, N2x2_G8, Bins3x3_G8, N3x3_G8, \
BinsCross_G8, NCross_G8 = \
    FE55_Extract_Database('H1RG-022','ASIC','08Jun29','1000', \
         QueryString='WHERE ADCs=8 and ADU2x2 !=0',Gain=1, NumADCs=1, \
         WindowM=1, Mode=1, UseRef=1, TvFlag=0)

mplot.figure(0)
mplot.clf()
mplot.plot(Bins2x2_G8, N2x2_G8, 'r-', label='8 Ch')
mplot.plot(Bins2x2_G6, N2x2_G6, 'b-', label='6 Ch')
mplot.plot(Bins2x2_G4, N2x2_G4, 'g-', label='4 Ch')
mplot.plot(Bins2x2_G2, N2x2_G2, 'm-', label='2 Ch')
mplot.plot(Bins2x2_G1, N2x2_G1, 'c-', label='1 Ch')
mplot.title(r'\textbf{Fe$^{55}$ with Gain = 1 and 2x2 bins}')

mplot.figure(1)
mplot.clf()
mplot.plot(Bins3x3_G8, N3x3_G8, 'r-', label='8 Ch')
mplot.plot(Bins3x3_G6, N3x3_G6, 'b-', label='6 Ch')
mplot.plot(Bins3x3_G4, N3x3_G4, 'g-', label='4 Ch')
mplot.plot(Bins3x3_G2, N3x3_G2, 'm-', label='2 Ch')
mplot.plot(Bins3x3_G1, N3x3_G1, 'c-', label='1 Ch')
mplot.title(r'\textbf{Fe$^{55}$ with Gain = 1 and 2x2 bins}')

Bins1x1_G1, N1x1_G1, Bins2x2_G1, N2x2_G1, Bins3x3_G1, N3x3_G1, \
BinsCross_G1, NCross_G1 = \
    FE55_Extract_Database('H1RG-022','ASIC','08Jun29','1000', \
         QueryString='WHERE ADCs=1 and ADU2x2 !=0',Gain=4, NumADCs=1, \
         WindowM=1, Mode=1, UseRef=1, TvFlag=0)

Bins1x1_G2, N1x1_G2, Bins2x2_G2, N2x2_G2, Bins3x3_G2, N3x3_G2, \
BinsCross_G2, NCross_G2 = \
    FE55_Extract_Database('H1RG-022','ASIC','08Jun29','1000', \
         QueryString='WHERE ADCs=2 and ADU2x2 !=0',Gain=4, NumADCs=1, \
         WindowM=1, Mode=1, UseRef=1, TvFlag=0)

Bins1x1_G4, N1x1_G4, Bins2x2_G4, N2x2_G4, Bins3x3_G4, N3x3_G4, \
BinsCross_G4, NCross_G4 = \
    FE55_Extract_Database('H1RG-022','ASIC','08Jun29','1000', \
         QueryString='WHERE ADCs=4 and ADU2x2 !=0',Gain=4, NumADCs=1, \
         WindowM=1, Mode=1, UseRef=1, TvFlag=0)

Bins1x1_G6, N1x1_G6, Bins2x2_G6, N2x2_G6, Bins3x3_G6, N3x3_G6, \
BinsCross_G6, NCross_G6 = \
    FE55_Extract_Database('H1RG-022','ASIC','08Jun29','1000', \
         QueryString='WHERE ADCs=6 and ADU2x2 !=0',Gain=4, NumADCs=1, \
         WindowM=1, Mode=1, UseRef=1, TvFlag=0)

Bins1x1_G8, N1x1_G8, Bins2x2_G8, N2x2_G8, Bins3x3_G8, N3x3_G8, \
BinsCross_G8, NCross_G8 = \
    FE55_Extract_Database('H1RG-022','ASIC','08Jun29','1000', \
         QueryString='WHERE ADCs=8 and ADU2x2 !=0',Gain=4, NumADCs=1, \
         WindowM=1, Mode=1, UseRef=1, TvFlag=0)

mplot.figure(2)
mplot.clf()
mplot.plot(Bins2x2_G8, N2x2_G8, 'r-', label='Gain 8')
mplot.plot(Bins2x2_G6, N2x2_G6, 'b-', label='Gain 6')
mplot.plot(Bins2x2_G4, N2x2_G4, 'g-', label='Gain 4')
mplot.plot(Bins2x2_G2, N2x2_G2, 'm-', label='Gain 2')
mplot.plot(Bins2x2_G1, N2x2_G1, 'c-', label='Gain 1')

Bins1x1_G1, N1x1_G1, Bins2x2_G1, N2x2_G1, Bins3x3_G1, N3x3_G1, \
BinsCross_G1, NCross_G1 = \
    FE55_Extract_Database('H1RG-022','ASIC','08Jun29','1000', \
         QueryString='WHERE ADCs=1 and ADUCross !=0',Gain=4, NumADCs=1, \
         WindowM=1, Mode=1, UseRef=1, TvFlag=0)

Bins1x1_G2, N1x1_G2, Bins2x2_G2, N2x2_G2, Bins3x3_G2, N3x3_G2, \
BinsCross_G2, NCross_G2 = \
    FE55_Extract_Database('H1RG-022','ASIC','08Jun29','1000', \
         QueryString='WHERE ADCs=2 and ADUCross !=0',Gain=4, NumADCs=1, \
         WindowM=1, Mode=1, UseRef=1, TvFlag=0)

Bins1x1_G4, N1x1_G4, Bins2x2_G4, N2x2_G4, Bins3x3_G4, N3x3_G4, \
BinsCross_G4, NCross_G4 = \
    FE55_Extract_Database('H1RG-022','ASIC','08Jun29','1000', \
         QueryString='WHERE ADCs=4 and ADUCross !=0',Gain=4, NumADCs=1, \
         WindowM=1, Mode=1, UseRef=1, TvFlag=0)

Bins1x1_G6, N1x1_G6, Bins2x2_G6, N2x2_G6, Bins3x3_G6, N3x3_G6, \
BinsCross_G6, NCross_G6 = \
    FE55_Extract_Database('H1RG-022','ASIC','08Jun29','1000', \
         QueryString='WHERE ADCs=6 and ADUCross !=0',Gain=4, NumADCs=1, \
         WindowM=1, Mode=1, UseRef=1, TvFlag=0)

Bins1x1_G8, N1x1_G8, Bins2x2_G8, N2x2_G8, Bins3x3_G8, N3x3_G8, \
BinsCross_G8, NCross_G8 = \
    FE55_Extract_Database('H1RG-022','ASIC','08Jun29','1000', \
         QueryString='WHERE ADCs=8 and ADUCross !=0',Gain=4, NumADCs=1, \
         WindowM=1, Mode=1, UseRef=1, TvFlag=0)

mplot.figure(3)
mplot.clf()
mplot.plot(BinsCross_G8, NCross_G8, 'r-', label='Gain 8')
mplot.plot(BinsCross_G6, NCross_G6, 'b-', label='Gain 6')
mplot.plot(BinsCross_G4, NCross_G4, 'g-', label='Gain 4')
mplot.plot(BinsCross_G2, NCross_G2, 'm-', label='Gain 2')
mplot.plot(BinsCross_G1, NCross_G1, 'c-', label='Gain 1')