#!/usr/bin/python
import numpy
from pylab import *
y= numpy.array([[0.214,0.133,0.144,0.161,0.170,0.085],
[0.116,0.083,0.109,0.145,0.164,0.083],
[0.066,0.061,0.095,0.139,0.163,0.083],
[0.055,0.056,0.093,0.138,0.162,0.083]])
priors = [0.1,0.05,0.02,0.01]
z=numpy.array([0.1,0.3,0.5,0.7,0.9,1.5])
for i in range(len(priors)):
plot(z,y[i,:],label='%d%%' % (priors[i]*100))
legend(title='Prior on lens mass')
xlabel('Redshift (bin centers)',fontdict={'fontsize':16})
ylabel('Uncertainty on bin fractions',fontdict={'fontsize':14})
show()