\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*} """
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 = 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')
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')
start=datetime.datetime.now() from input import * flist=[] #samples=samplespowheg #samples.remove('stop') for s in samples: for f in groups[s]['sets']: flist.append(f) flist.append("data") flist.append("ap-350") flist.append("pp-350") from extLibFull import filterByCutHDF,plotHDF,getModel model=getModel('BP4') from pylab import * print "Starting processing...", fnam="BP4-4l" #fnam="BP4-3l" varList=['mp'] nbin=[25] xl=[r'$m(\ell^\pm\ell^\pm)$ [GeV]'] yl=[r'Events / 20 GeV'] xr=[[10,510]] #xl=[r'$E^{miss}_{T}$ [GeV]'] #yl=[r'N events / 5 GeV'] #xr=[[0,100]] ll=[1]