def plot_hist(filename, N=400000):
values = read_cosmos_table(filename, ["zp_best", "type"], N)
# select only galaxies
i = numpy.where(values[:, 1] == 0)
values = values[i[0], :]
# select only objects with valid photometric redshift
i = numpy.where(~numpy.isnan(values[:, 0]))
values = values[i[0], :]
print "selected %i sources" % values.shape[0]
pylab.hist(values[:, 0], bins=300, linewidth=0.1)
pylab.xlabel("z_phot")
pylab.ylabel("N(z)")
pylab.text(
0.95, 0.95, "%i sources" % values.shape[0], transform=pylab.gca().transAxes, fontsize=14, va="top", ha="right"
)
pylab.title("COSMOS photo-z distribution")
def plot_hist(filename, N=400000):
values = read_cosmos_table(filename, ['zp_best', 'type'], N)
#select only galaxies
i = numpy.where(values[:, 1] == 0)
values = values[i[0], :]
#select only objects with valid photometric redshift
i = numpy.where(~numpy.isnan(values[:, 0]))
values = values[i[0], :]
print "selected %i sources" % values.shape[0]
pylab.hist(values[:, 0], bins=300, linewidth=0.1)
pylab.xlabel('z_phot')
pylab.ylabel('N(z)')
pylab.text(0.95,
0.95,
"%i sources" % values.shape[0],
transform=pylab.gca().transAxes,
fontsize=14,
va='top',
ha='right')
pylab.title('COSMOS photo-z distribution')
def plot_by_mag(filename, bandid="i_auto", N=400000, mag_bins=numpy.arange(16, 26)):
values = read_cosmos_table(filename, ["zp_best", "type", bandid], N)
print " - %i total objects" % values.shape[0]
# select only galaxies
i = numpy.where(values[:, 1] == 0)
values = values[i[0], :]
print " - %i galaxies" % values.shape[0]
# select only objects with valid photometric redshift
i = numpy.where(~numpy.isnan(values[:, 0]))
values = values[i[0], :]
print " - %i with redshift" % values.shape[0]
# select only objects with valid <bandid> magnitude
i = numpy.where(~numpy.isnan(values[:, 2]))
values = values[i[0], :]
print " - %i with valid mag" % values.shape[0]
"""
mags = values[:,2]
mags = mags[numpy.where(mags!=-99)]
mag_min = min(mags)
mag_max = max(mags)
mag_min = int(numpy.floor(mag_min))
mag_max = int(numpy.ceil(mag_max))
mag_bins = numpy.arange(mag_min,mag_max+1)
"""
Nmag = len(mag_bins) - 1
pylab.figure(figsize=(10, 10))
for i in range(Nmag):
ind = numpy.where((values[:, 2] >= mag_bins[i]) & (values[:, 2] < mag_bins[i + 1]))
z = values[ind[0], 0]
pylab.subplot(331 + i)
print "%i sources with %.1f<%s<%.1f" % (len(z), mag_bins[i], bandid, mag_bins[i + 1])
# Nbins = len(z)/50
# Nbins = max(Nbins,30)
# Nbins = min(Nbins,300)
bins = numpy.linspace(0, 4, 201)
if bandid == "i_auto":
bandid = "I"
pylab.hist(z, bins, linewidth=0)
if i in (6, 7, 8):
pylab.xlabel("z_phot")
if i in (0, 3, 6):
pylab.ylabel("N(z)")
pylab.text(
0.95,
0.95,
"%.1f<%s<%.1f\n%i sources" % (mag_bins[i], bandid, mag_bins[i + 1], len(z)),
transform=pylab.gca().transAxes,
fontsize=12,
va="top",
ha="right",
)
pylab.xlim(0, 4)
pylab.gca().xaxis.set_major_locator(ticker.MultipleLocator(1))
def plot_by_mag(filename,
bandid='i_auto',
N=400000,
mag_bins=numpy.arange(16, 26)):
values = read_cosmos_table(filename, ['zp_best', 'type', bandid], N)
print " - %i total objects" % values.shape[0]
#select only galaxies
i = numpy.where(values[:, 1] == 0)
values = values[i[0], :]
print " - %i galaxies" % values.shape[0]
#select only objects with valid photometric redshift
i = numpy.where(~numpy.isnan(values[:, 0]))
values = values[i[0], :]
print " - %i with redshift" % values.shape[0]
#select only objects with valid <bandid> magnitude
i = numpy.where(~numpy.isnan(values[:, 2]))
values = values[i[0], :]
print " - %i with valid mag" % values.shape[0]
"""
mags = values[:,2]
mags = mags[numpy.where(mags!=-99)]
mag_min = min(mags)
mag_max = max(mags)
mag_min = int(numpy.floor(mag_min))
mag_max = int(numpy.ceil(mag_max))
mag_bins = numpy.arange(mag_min,mag_max+1)
"""
Nmag = len(mag_bins) - 1
pylab.figure(figsize=(10, 10))
for i in range(Nmag):
ind = numpy.where( (values[:,2]>=mag_bins[i])\
&(values[:,2]<mag_bins[i+1]) )
z = values[ind[0], 0]
pylab.subplot(331 + i)
print "%i sources with %.1f<%s<%.1f" % (len(z), mag_bins[i], bandid,
mag_bins[i + 1])
#Nbins = len(z)/50
#Nbins = max(Nbins,30)
#Nbins = min(Nbins,300)
bins = numpy.linspace(0, 4, 201)
if bandid == 'i_auto': bandid = 'I'
pylab.hist(z, bins, linewidth=0)
if i in (6, 7, 8):
pylab.xlabel('z_phot')
if i in (0, 3, 6):
pylab.ylabel('N(z)')
pylab.text(0.95,
0.95,
"%.1f<%s<%.1f\n%i sources" %
(mag_bins[i], bandid, mag_bins[i + 1], len(z)),
transform=pylab.gca().transAxes,
fontsize=12,
va='top',
ha='right')
pylab.xlim(0, 4)
pylab.gca().xaxis.set_major_locator(ticker.MultipleLocator(1))