Y = [d_pois([R.Double(x)],pars) for x in X]
rplt.hist(h)
plot(X,Y)
xlim(0,80)
# <codecell>
%run ../python/oana_setup.py
# <codecell>
import ROOT as R
h10.SetEntryList(0)
h10.Draw('>>elist', 'top3pass && mmthreshpass && mmpi0pass && fidpass', 'entrylist')
elist = R.gDirectory().Get('elist')
h10.SetEntryList(elist)
# <codecell>
h = rh1('nphe[cc[0]-1]', 'ccseg==5 && ccsect==1', binning=(300,0,300))
# <codecell>
rplt.hist(h)
# <codecell>
h10.Draw('ccseg>>hseg(19,-0.5,18.5)', '','',1000)
# <codecell>
