PerformanceCurvePlotter Namespace Reference

Classes
class	TauValidationInfo

Functions
def	PlotPerformanceCurves

Variables
	HideMuonAndElectronRej = True
string	KinematicVar = 'pt'
int	KinematicVarMax = 50
int	KinematicVarMin = 20
string	LatexVar = "True p_{T}"
string	VarUnit = "GeV/c"

Detailed Description

  PerformanceCurvePlotter

  Author: Evan K. Friis, UC Davis (friis@physics.ucdavis.edu)

  Plot the efficiency for a given fake rate 
  for a series of cuts in a set pt range.

  With the default cuts (20 < pt < 50), this is designed
  to reproduce (independently) the TancBenchmark plots produced in
  RecoTauTag/TauTagTools/test/MVABenchmarks.py


      Instructions:

      Add the releases to as shown.  You must specify and signal 
      and backround file for each release, and a descriptive label.
      To choose the chain of discriminators to plot, select the appropriate
      EDProducer defined in Validation.RecoTau.TauTagValidation_cfi

      PreFourTanc = TauValidationInfo("Signal.root", 
                                      "Background.root",
                                      "310pre5 TaNC",
                                      Validation.RecoTau.RecoTauValidation_cfi.RunTancValidation)

      PreFourIso = TauValidationInfo("Signal.root", 
                                      "Background.root",
                                      "310pre4 Iso",
                                      Validation.RecoTau.RecoTauValidation_cfi.PFTausHighEfficiencyLeadingPionBothProngs)


      ReleasesToBuildFrom = []
      ReleasesToBuildFrom.append(PreFourTanc)
      ReleasesToBuildFrom.append(PreFourIso)

      PlotPerformanceCurves(ReleasesToBuildFrom, "output.png")

Function Documentation

def PerformanceCurvePlotter.PlotPerformanceCurves	(		ReleasesToBuildFrom,
			OutputFile
	)

Definition at line 184 of file PerformanceCurvePlotter.py.

def PlotPerformanceCurves(ReleasesToBuildFrom, OutputFile):
   # Setup
   myCanvas = TCanvas("Validation", "Validation", 800, 800)
   # Cut label on the plots
   CutLabel = TPaveText(0.6, 0.1, 0.9, 0.27, "brNDC")
   CutLabel.AddText("%0.1f%s < %s < %0.1f%s" % (KinematicVarMin, VarUnit, LatexVar, KinematicVarMax, VarUnit))
   CutLabel.SetFillStyle(0)
   CutLabel.SetBorderSize(0)

   # Build the TGraphs of sigEff versus bkgFakeRate for each release
   for aRelease in ReleasesToBuildFrom:
      aRelease.LoadHistograms()
      aRelease.ComputeEfficiencies()
      aRelease.BuildTGraphSummary()

   CurrentHistogram = 0

   for aRelease in ReleasesToBuildFrom:
      if CurrentHistogram == 0:
         aRelease.SummaryTGraph.Draw("ALP")
         CurrentHistogram = aRelease.SummaryTGraph.GetHistogram()
      else:
         aRelease.SummaryTGraph.Draw("LP")


   CurrentHistogram.SetAxisRange(0, 1)
   CurrentHistogram.GetYaxis().SetRangeUser(0.0001,1)
   CurrentHistogram.GetXaxis().SetTitle("Efficiency")
   CurrentHistogram.GetYaxis().SetTitle("Fake Rate")
   CurrentHistogram.SetTitle("Performance Points")

   gPad.SetLogy(True)

   for aRelease in ReleasesToBuildFrom:
      for aPoint in aRelease.DiscriminatorPoints:
         aPoint['PerformancePoint'].Draw("P")

   TauValidationInfo.SummaryLegend.Draw()
   TauValidationInfo.DiscriminatorLegend.Draw()
   CutLabel.Draw()

   myCanvas.SaveAs(OutputFile)

Variable Documentation

PerformanceCurvePlotter.HideMuonAndElectronRej = True

Definition at line 57 of file PerformanceCurvePlotter.py.

string PerformanceCurvePlotter.KinematicVar = 'pt'

Definition at line 50 of file PerformanceCurvePlotter.py.

int PerformanceCurvePlotter.KinematicVarMax = 50

Definition at line 53 of file PerformanceCurvePlotter.py.

int PerformanceCurvePlotter.KinematicVarMin = 20

Definition at line 52 of file PerformanceCurvePlotter.py.

string PerformanceCurvePlotter.LatexVar = "True p_{T}"

Definition at line 51 of file PerformanceCurvePlotter.py.

string PerformanceCurvePlotter.VarUnit = "GeV/c"

Definition at line 54 of file PerformanceCurvePlotter.py.