from posteriors import *
import pylab
x,y=load('plot_data/square_p2.dat',1,2)
#The temperature width.
w=0.13
xcal,ycal=calibration(x,y,w,steps=200)
uncal_total=tab_area(x,y,x[0],x[-1])
uncal_zero=tab_area(x,y,x[0],0.0)
cal_total=tab_area(xcal,ycal,xcal[0],xcal[-1])
cal_zero=tab_area(xcal,ycal,xcal[0],0.0)
print "uncalibrated:",uncal_zero/uncal_total
print "calibrated  :",cal_zero/cal_total
y=y/uncal_total
ycal=ycal/cal_total
pylab.plot(x,y)
pylab.plot(xcal,ycal)
pylab.show()
f=open('calibrated.dat','w')
for i in range(len(xcal)):
    f.write('%f\t%e\n'% (xcal[i],ycal[i]))
f.close()
f=open('uncalibrated.dat','w')
for i in range(len(x)):
    f.write('%f\t%e\n'% (x[i],y[i]))
f.close()