##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])