\hline"""
print "Cut & ",
for s in samples:
print groups[s]['mplabel']+' & ',
print 'Total & Data & Signal \\\\'
print '\\hline'
for il in [0,1]:
nl=il+3
cut="(mp>"+str(mmin)+") & (mp<"+str(mmax)+")"
if il:
cut+=" & (mn>"+str(mmin)+") & (mn<"+str(mmax)+")"
if mod in extras.keys():
cut+=" & "+extras[mod][il]
for c in range(len(cuts[cMod][il])):
cut+=" & "+cuts[cMod][il][c].substitute(mass=m)
filt=filterByCutHDF(flist+['ap-'+str(m),'pp-'+str(m)],cut,["mp"],getModel(mod))
yld=getYieldHDF(filt,samples,getModel(mod),addErr=False)
print str(nl)+'l-'+str(c+1),'&',
for s in samples:
print "$"+str(round(yld[s],2))+"\pm"+str(round(yld[s+'-err'],2))+'$ &',
print "$"+str(round(yld['bg'],2))+'\pm'+str(round(yld['bg-err'],2))+'$ & $'+str(yld['data'])+'\pm'+str(round(yld['data-err'],2))+'$ & $'+str(round(yld['signal'],2))+'\pm'+str(round(yld['signal-err'],2))+'$ \\\\'
print '\\hline'
print """
\end{tabular}
}
\end{center}
\end{table*}
"""
vlist=var3
x=xl3
y=yl
nb=nb3
xr=xr3
ll=1
lloc=1
if c: lloc=3
sc='linear'
vlistFilt=["mp"]
if nl==4:
vlist=var4
vlistFilt=var4
x=xl4
y=yl4
nb=nb4
xr=xr4
ll=[1,1]
lloc=[1,1]
if c: lloc=[3,3]
sc=['linear','linear']
print il,c,cut
filt=filterByCutHDF(flist+['ap-'+str(m),'pp-'+str(m)],cut,vlistFilt,getModel(mod))
fname=mod+'-m'+str(m)+'-'+str(nl)+'l-cut'+str(c)
plotHDF(filt,samples,'appp',m,getModel(mod),vlist,x,y,nb,xr,fname+'.pdf',lloc)
savefig(fname+'.png')
draw()
fname=mod+'-m'+str(m)+'-'+str(nl)+'l-cut'+str(c)+'-lin'
plotHDF(filt,samples,'appp',m,getModel(mod),vlist,x,y,nb,xr,fname+'.pdf',ll,sc)
savefig(fname+'.png')
cutCR="(ntau=="+str(nt)+") & ((mp<"+str(mmin)+") | (mp>"+str(mmax)+"))"
cutSR="(ntau=="+str(nt)+") & (mp>"+str(mmin)+") & (mp<"+str(mmax)+")"
if il:
cutCR+=" & ((mn<"+str(mmin)+") | (mn>"+str(mmax)+"))"
cutSR+=" & (mn>"+str(mmin)+") & (mn<"+str(mmax)+")"
if mod in extras.keys():
cutCR+=" & "+extras[mod][il]
cutSR+=" & "+extras[mod][il]
cutCR+=" & "+cuts[cMod][il][0].substitute(mass=mass)
for c in range(mx):
cutSR+=" & "+cuts[cMod][il][c].substitute(mass=mass)
print "Debug:",il,nt,cutCR,cutSR
crdata=filterByCutHDF(flist2,cutCR,["mp"],getModel(mod))
srdata=filterByCutHDF(flist2,cutSR,["mp"],getModel(mod))
yldCR=getYieldHDF(crdata,samples,getModel(mod),addErr=True)
yldSR=getYieldHDF(srdata,samples,getModel(mod),addErr=True)
nDataCR[il][nt]=yldCR['data']
nDataSR[il][nt]=yldSR['data']
nBGCR[il][nt]=yldCR['bg']
nBGSR[il][nt]=yldSR['bg']
#for s in samples:
# nBGCR[il][nt]+=yldCR[s]
# nBGSR[il][nt]+=yldSR[s]
nSigPP[il][nt]+=yldSR['pp-'+str(mass)]
nSigAP[il][nt]+=yldSR['ap-'+str(mass)]
xr=[[10,510]]
#xl=[r'$E^{miss}_{T}$ [GeV]']
#yl=[r'N events / 5 GeV']
#xr=[[0,100]]
ll=[1]
#cut1='~lowMassRes & (nlep>3) & (isoworst<0.35) & (sipworst<4) & (tWorstId>1)'
cut1='~lowMassRes & (nlep>3) & (ntau<3) & (isoworst<0.35) & (sipworst<4) & (tWorstId>1)' # 4-lep
#cut1='~lowMassRes & (nlep==3) & (ntau<2) & (isoworst<0.35) & (sipworst<4) & (tWorstId>1)' # 3-lep
#cut1='(nlep==3) & (isoworst<0.35) & (sipworst<4) & (sumpt>0.85*275+125) & (dz>80) & (phi1<0.005*275+1.15) & (met>20)'
#cut1='(nlep==3) & (ntau==0) & (isoworst<0.35) & (sipworst<4) & (sumpt>0.85*275+125) & (dz>80) & (phi1<0.005*275+1.15) & (met>20) & ~lowMassRes & (tWorstId>1)'
#cut2='(nlep==3) & (ntau==1) & (isoworst<0.35) & (sipworst<4) & (sumpt>0.85*275+125) & (dz>80) & (phi1<0.005*275+1.15) & (met>20) & ~lowMassRes & (tWorstId>1)'
#cut1='recoPass & (nlep==3) & (isoworst<0.35) & (sipworst<4) & (sumpt>0.85*130+125) & (dz>80) & (phi1<0.005*130+1.15) & (met>20) & ~lowMassRes'
#cut1='recoPass & (nlep==3) & (ntau==0) & (isoworst<0.35) & (sipworst<4) & (sumpt>0.85*130+125) & (dz>80) & (phi1<0.005*130+1.15) & (met>20) & ~lowMassRes'
#cut1='recoPass & (nlep==3) & (ntau==0) & (isoworst<0.35) & (sipworst<4) & (sumpt>0.85*130+125) & (dz>80) & (phi1<0.005*130+1.15) & (met>20) & ~lowMassRes'
#cut2='recoPass & (nlep==3) & (ntau==1) & (isoworst<0.35) & (sipworst<4) & (sumpt>0.85*130+125) & (dz>80) & (phi1<0.005*130+1.15) & (met>20) & ~lowMassRes'
#cut1='(nlep==3) & (ntau<2) & (isoworst<0.35) & (sipworst<4) & ~lowMassRes & recoPass & (tWorstId>1) & (sumpt>0.85*250+125) & (dz>80) & (phi1<0.005*250+1.15) & (met>20)'
#cut1='recoPass & (nlep==3) & (ntau==0) & (isoworst<0.35) & (sipworst<4)'
#cut1='recoPass & ~lowMassRes & (nlep==3) & (isoworst<0.35) & (sipworst<4) & (tWorstId>1) & (dz<30) & (sumpt>150)'
#cut1='recoPass & ~lowMassRes & (nlep>2) & (isoworst<0.35) & (sipworst<4) & (tWorstId>1)'
#cut1='recoPass & ~lowMassRes & (nlep==3) & (isoworst<0.35) & (sipworst<4) & (dz<80) & (sumpt>0.85*275+125) & (phi1<0.005*275+1.15) & (met>20)'
filt1=filterByCutHDF(flist,cut1,varList,model)
#filt2=filterByCutHDF(flist,cut2,varList,model)
#fnam="BP4-4l"
plotHDF(filt1,samples,"appp",350,getModel("BP4"),varList,xl,yl,nbin,xr,fnam+'.pdf',ll,showRatio=False)
savefig(fnam+'.png')
draw()
#fnam="BP4-3l-1tau-medTau"
#plotHDF(filt2,samples,"appp",275,getModel("BP4"),varList,xl,yl,nbin,xr,fnam+'.png')
#draw()
