#!/usr/bin/python 

from pylab import *
import numpy
from matplotlib.ticker import NullFormatter,MultipleLocator, FormatStrFormatter


nbin = 44
range = [-0.022,0.022]
islog=False
majorLocator   = MultipleLocator(0.01)
#majorFormatter = FormatStrFormatter('%d')
myxlim= [-0.024,0.024]
myylim= [0,88]

Rdata = loadtxt('Rresids')

ax=subplot(221)
data = loadtxt('Bresids')
hist(data,histtype='stepfilled',color='b',alpha=0.4,bins=nbin,range=range,log=islog)
hist(Rdata,histtype='stepfilled',color='k',alpha=0.1,bins=nbin,range=range,log=islog)
xlim(myxlim)
ylim(myylim)
text(0.015,70,"B",fontdict={'fontsize':20})
ax.xaxis.set_major_formatter( NullFormatter() )

ax = subplot(222)
data = loadtxt('Vresids')
hist(data,histtype='stepfilled',color='g',alpha=0.4,bins=nbin,range=range,log=islog)
hist(Rdata,histtype='stepfilled',color='k',alpha=0.1,bins=nbin,range=range,log=islog)
xlim(myxlim)
ylim(myylim)
text(0.015,70,"V",fontdict={'fontsize':20})
ax.xaxis.set_major_formatter( NullFormatter() )
ax.yaxis.set_major_formatter( NullFormatter() )


ax=subplot(223)
hist(Rdata,histtype='stepfilled',color='r',alpha=0.4,bins=nbin,range=range,log=islog)
xlim(myxlim)
ylim(myylim)
text(0.015,70,"R",fontdict={'fontsize':20})
ax.xaxis.set_major_locator( majorLocator )

ax=subplot(224)
data = loadtxt('zresids')
hist(data,histtype='stepfilled',color='k',alpha=0.4,bins=nbin,range=range,log=islog)
hist(Rdata,histtype='stepfilled',color='k',alpha=0.1,bins=nbin,range=range,log=islog)
xlim(myxlim)
ylim(myylim)
text(0.015,70,"z",fontdict={'fontsize':20})
ax.yaxis.set_major_formatter( NullFormatter() )
ax.xaxis.set_major_locator( majorLocator )
#ax.xaxis.set_major_formatter(majorFormatter)

figtext(0.3,0.02,"exposure residual (mag)",fontdict={'fontsize':18})
figtext(0.04,0.6,"# of exposures",fontdict={'fontsize':18},rotation=90)
subplots_adjust(wspace=0,hspace=0)
show()