cmbexp = 'CMBS4'
fsky, thetab, DeltaT, iterative = bolometers[cmbexp]['fsky'], bolometers[cmbexp]['thetab'],\
bolometers[cmbexp]['DeltaT'], bolometers[cmbexp]['iterative']
band = 'BX'
setup = {'BX': {'colors': ['goldenrod', 'tan', 'y'], 'maglim': 25.5, 'decband': 'R'},\
'u': {'colors': ['darkblue', 'deepskyblue', 'b'], 'maglim': 24.6, 'decband': 'i'},\
'g': {'colors': ['darkgreen', 'limegreen', 'g'], 'maglim': 25.8, 'decband': 'i'},\
'r': {'colors': ['darkred', 'indianred', 'r'], 'maglim': 25.8, 'decband': 'z'}}
mlim = setup[band]['maglim']
pz = get_pz(band)
bz = lambda z: bz_callmodel(z, mlim)
nbar = get_shot(band, mlim)
decband = setup[band]['decband']
colors = setup[band]['colors']
peakz = _peakz(pz)
## Hard p(z) limits.
zmin = peakz - 2.00
zmax = peakz + 2.00
##
pl.clf()
band = 'u'
setup = {'BX': {'colors': ['goldenrod', 'tan', 'y'], 'spec-maglim': 24.0, 'phot-maglim': 25.5,'decband': 'R', 'file': 'BXDrop.dat'},\
'u': {'colors': ['darkblue', 'deepskyblue', 'b'], 'spec-maglim': 24.0, 'phot-maglim': 24.6,'decband': 'R', 'file': 'MalkanDrop.dat'},\
'g': {'colors': ['darkgreen', 'limegreen', 'g'], 'spec-maglim': 25.5, 'phot-maglim': 25.8,'decband': 'i', 'file': 'gDrop.dat'},\
'r': {'colors': ['darkred', 'indianred', 'r'], 'spec-maglim': 25.5, 'phot-maglim': 25.8,'decband': 'z', 'file': 'rDrop.dat'}}
file = setup[band]['file']
colors = setup[band]['colors']
spec_mlim = setup[band]['spec-maglim']
phot_mlim = setup[band]['phot-maglim']
## Spectroscopic sample p(z) and b(z).
ps = get_pz(band)
bs = lambda z: bz_callmodel(z, spec_mlim)
## Get the z percentiles for this dropout p(z).
percentiles = percentiles(ps, printit=True)
## Photometric sample p(z).
pp = get_pz(band)
peakz = _peakz(ps)
zmin = peakz - 1.0
zmax = peakz + 1.0
## Get dropout Schechter photometric sample counts for given band.
root = os.environ['LBGCMB']