def __init__(self, pzrun, nchunk=50):
"""
Take objects from the given pzrun and select:
1) ones that are "recoverable"
2) lie within the bounds of one of our test
regions; these overlap spectroscopic regions.
"""
self.pzrun=pzrun
self.nchunk=nchunk
self.conf = zphot.cascade_config(pzrun)
self.wo = zphot.weighting.WeightedOutputs()
# load columns
self.procrun = self.conf['photo']['procrun']
self.cols = sdssgal.open_columns(self.procrun)
def plot_lens_s2n_bysample(lens_sample, pzrun, cumulative=True):
import zphot
conf=lensing.files.read_config('lcat', lens_sample)
zconf=zphot.cascade_config(pzrun)
#
# generate the expected n(z)
#
# read the weights and create weighted histogram
wstruct = zphot.weighting.read_weights(zconf['pofz']['wrun'], 2)
binsize=0.0314
zmax=1.1
hdict = eu.stat.histogram(wstruct['z'], binsize=binsize, weights=wstruct['weight'],max=zmax)
zs = hdict['center']
pzs = hdict['whist']
# now read the zl
data = lensing.files.read_original_catalog(type='lens',sample=lens_sample)
return plot_lens_s2n(data['z'], hdict['low'], hdict['high'], pzs)
