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]")