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')
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()