efficiencyPlots = [] effPlots = EfficiencyPlots() effPlots.name = name effPlots.publicationDir = publicationDir effPlots.histoBaseName = "hFakeRate" effPlots.inputFileNames = [inputFileNames] effPlots.systems = systems effPlots.selectionLevels = selectionLevels effPlots.plotInfos = plotInfos effPlots.referenceLevels = referenceLevels effPlots.individualNames = names effPlots.variables = variables effPlots.variableNames = variableNames effPlots.outputFile = outputFile effPlots.plot(0., 0.5) efficiencyPlots.append(effPlots) effPlots = EfficiencyPlots() effPlots.name = name effPlots.publicationDir = publicationDir effPlots.histoBaseName = "hFakeRate" effPlots.inputFileNames = [inputFileNames] effPlots.systems = systems effPlots.selectionLevels = selectionLevels2 effPlots.plotInfos = plotInfos effPlots.referenceLevels = referenceLevels2 effPlots.individualNames = names2 effPlots.variables = variables effPlots.variableNames = variableNames effPlots.outputFile = outputFile
plotInfos = [PlotInfo()]
plotInfos[0].markerStyle = 20
plotInfos[0].yTitle = "Fake factor"
if not os.path.exists(plotDir + "/" + name):
os.makedirs(plotDir + "/" + name)
outputFile = ROOT.TFile.Open(plotDir + "/" + name + "/" + name + ".root",
"RECREATE")
efficiencyPlots = []
effPlots = EfficiencyPlots()
effPlots.name = name
effPlots.publicationDir = publicationDir
effPlots.histoBaseName = "hFakeRate"
effPlots.inputFileNames = [inputFileNames]
effPlots.systems = systems
effPlots.selectionLevels = selectionLevels2
effPlots.plotInfos = plotInfos
effPlots.referenceLevels = referenceLevels2
effPlots.individualNames = names2
effPlots.variables = variables
effPlots.variableNames = variableNames
effPlots.outputFile = outputFile
effPlots.divideOption = "pois"
effPlots.plot(0., 0.5)
efficiencyPlots.append(effPlots)
outputFile.Close()
effPlots2 = EfficiencyPlots()
effPlots2.name = name
effPlots2.publicationDir = publicationDir
effPlots2.histoBaseName = "hFakeRate"
effPlots2.inputFileNames = [[inputFileName]]
effPlots2.systems = systems
effPlots2.selectionLevels = selectionLevels
effPlots2.plotInfos = plotInfos
effPlots2.referenceLevels = referenceLevels
effPlots2.individualNames = names
effPlots2.variables = variables
effPlots2.variableNames = variableNames
effPlots2.outputFile = outputFile
effPlots2.divideOption = "pois"
effPlots2.rebin = [0.,10.,20.,30.,40.,50.,60.,70.,200.]
effPlots2.plot(0., 0.3)
efficiencyPlots.append(effPlots2)
########################################################
## Compute corrections for mT>70GeV corrections
def graphDivide(graph1, graph2):
x = []
y = []
xe = []
ye = []
for p in range(0, graph1.GetN()):
x.append(graph1.GetX()[p])
