##PURPOSE:
## To plot the value of volcano pixels with the regular pixels
##
##

from Return_CR_DbQuantities import *
execfile('/afs/slac/u/ki/lances/python/python_HxRG/HxRG_Setup.py')

PlotDir = '/nfs/slac/g/ki/ki03/lances/H1RG-022/Plots/Cosmics/'
TvFlag='null'
MinCosEneDep = 0000
QueryString=' AND COSENEDEP > '+str(MinCosEneDep)
MaxCosEneDep  = 35000.
CosEneDepInds = arange(MaxCosEneDep)

#Plot the Fermi Distribution
Et = 0.100 #meV
Ef = 1.2 
k  = 8.617e-5
Temps  = [100., 110., 120., 130., 140., 150.]
Colors = ['g','r','b','k','c','m']
Es     = frange(10000)/2000
mplot.figure(1)
mplot.clf()
for i, Temp in zip(arange(size(Temps)), Temps):
  print Temps[i]
  mplot.plot(Es, 1/(1+exp((Es-Ef)/(k*Temp))), color=Colors[i])

#Plot the fraction q_lost/d_dep
mplot.figure(0)
mplot.clf()
mplot.hold(True)
TempOff=80
Et_m_Ef=0.0254 #ev
Temp=frange(2000)
Temp=Temp+TempOff
q_frac=exp((-Et_m_Ef)/(k*Temp))*(T-80)**.5/(1+exp((-Et_m_Ef)/(k*Temp)))
mplot.plot(Temp,q_frac)
xlim(100,200)

q_frac_120 = q_frac[120-TempOff]
print 'q at 120 :' +str(q_frac_120)
q_frac_130 = q_frac[130-TempOff]
print 'q at 130 :' +str(q_frac_130)
q_frac_140 = q_frac[140-TempOff]
print 'q at 140 :' +str(q_frac_140)
q_frac_150 = q_frac[150-TempOff]
print 'q at 150 :' +str(q_frac_150)
q_frac_160 = q_frac[160-TempOff]
print 'q at 160 :' +str(q_frac_160)