c.Print(sys.argv[1]+".pdf")
def plotstuff3(cuts,cutname):
leg = TLegend(0.6,0.75,0.9,0.9)
events.Draw("uncVZ>>h1(150,-50,100)",makeCutString(cuts[0:2]),"colz")
leg.AddEntry(gDirectory.Get("h1"),"after base cuts")
gDirectory.Get("h1").SetLineColor(1)
gDirectory.Get("h1").SetTitle("Effect of "+cutname+" cut")
for i in range(2,len(cuts)):
events.Draw("uncVZ>>h{0}(150,-50,100)".format(i),makeCutString(cuts[0:i+1]),"colz,same")
leg.AddEntry(gDirectory.Get("h{0}".format(i)),"after {0} cut".format(cutnames[i]))
gDirectory.Get("h{0}".format(i)).SetLineColor(i)
leg.Draw()
gDirectory.Get("h1").SetTitle(cutname)
gDirectory.Get("h1").GetXaxis().SetTitle("vertex Z [mm]")
gDirectory.Get("h1").GetYaxis().SetTitle("events/mm")
c.Print(sys.argv[1]+".pdf")
c.SetLogy(1)
plotstuff2(makeCutString(cuts[0:2]),makeCutString(cuts),"vertex Z, |m(e+e-)-0.03|<0.0026 GeV")
#plotstuff(makeCutString(cuts),"allcuts")
plotstuff3(cuts,"vertex Z, |m(e+e-)-0.03|<0.0026 GeV")
for i in range(2,len(cuts)):
cutstring = flipCut(cuts,i)
cutname = cutnames[i]
plotstuff(makeCutString(cuts,i),makeCutString(cuts),flipCut(cuts,i),cutname)
c.Print(sys.argv[1]+".pdf]")
c.Print(sys.argv[1]+".pdf[")
events.Draw("triEndZ>>base1d(100,-5,100)",cuts[0],"colz,goff")
base1d = gDirectory.Get("base1d")
def plotstuff(cutstring,cutname):
c.Clear()
c.Divide(1,3)
c.cd(1)
events.Draw("uncVZ:triEndZ>>plot2d(100,-5,100,100,-50,100)",cutstring,"colz")
c.cd(2)
events.Draw("triEndZ>>plot1d(100,-5,100)",cutstring,"colz")
eff=gDirectory.Get("plot1d").Clone("effloss")
eff.Sumw2()
eff.Divide(base1d)
c.cd(3)
eff.Draw()
c.Modified()
c.Print(sys.argv[1]+".pdf","Title:"+cutname)
plotstuff(cuts[0],"basecuts")
plotstuff(makeCutString(cuts),"allcuts")
for i in range(1,len(cuts)):
cutstring = flipCut(cuts,i)
cutname = cuts[i]
plotstuff(cutstring,cutname)
c.Print(sys.argv[1]+".pdf]")
hs.Add(gDirectory.Get("h1"))
leg = TLegend(0.6,0.75,0.9,0.9)
events.Draw("eleClT-posClT>>h0(200,-10,10)",cuts[0],"goff")
leg.AddEntry(gDirectory.Get("h0"),"after base cuts")
hs.Add(gDirectory.Get("h0"))
gDirectory.Get("h0").SetLineColor(1)
#gDirectory.Get("h1").SetTitle(cutname)
for i in range(1,len(cuts)):
events.Draw("eleClT-posClT>>h{0}(100,-10,10)".format(i),makeCutString(cuts[0:i+1]),"goff")
leg.AddEntry(gDirectory.Get("h{0}".format(i)),"after {0} cut".format(cutnames[i]))
hs.Add(gDirectory.Get("h{0}".format(i)))
gDirectory.Get("h{0}".format(i)).SetLineColor(i+1)
hs.Draw("nostack")
leg.Draw()
hs.GetXaxis().SetTitle("tEle-tPos [ns]")
hs.GetYaxis().SetTitle("events/(0.1 ns)")
c.Print(sys.argv[1]+".pdf")
c.SetLogy(1)
plotstuff2(cuts[0],makeCutString(cuts),"Effect of cuts")
#plotstuff(makeCutString(cuts),"allcuts")
plotstuff3(cuts,"Effect of cuts")
for i in range(1,len(cuts)):
cutstring = flipCut(cuts,i)
cutname = cutnames[i]
plotstuff(makeCutString(cuts,i),makeCutString(cuts),flipCut(cuts,i),cutname)
c.Print(sys.argv[1]+".pdf]")
events.Draw("triEndZ>>base1d(100,-5,100)", cuts[0], "colz,goff")
base1d = gDirectory.Get("base1d")
def plotstuff(cutstring, cutname):
c.Clear()
c.Divide(1, 3)
c.cd(1)
events.Draw("uncVZ:triEndZ>>plot2d(100,-5,100,100,-50,100)", cutstring,
"colz")
c.cd(2)
events.Draw("triEndZ>>plot1d(100,-5,100)", cutstring, "colz")
eff = gDirectory.Get("plot1d").Clone("effloss")
eff.Sumw2()
eff.Divide(base1d)
c.cd(3)
eff.Draw()
c.Modified()
c.Print(sys.argv[1] + ".pdf", "Title:" + cutname)
plotstuff(cuts[0], "basecuts")
plotstuff(makeCutString(cuts), "allcuts")
for i in range(1, len(cuts)):
cutstring = flipCut(cuts, i)
cutname = cuts[i]
plotstuff(cutstring, cutname)
c.Print(sys.argv[1] + ".pdf]")
