def compareMCGentoMCUnfolded(action, infileName):
if action == "herwigOnPythia" or action == "pythiaOnPythia":
unfoldingWasDoneOn = "QCD_Pt-15to3000_TuneZ2star_Flat_HFshowerLibrary_7TeV_pythia6"
elif action == "pythiaOnHerwig" or action == "herwigOnHerwig":
unfoldingWasDoneOn = "QCD_Pt-15to1000_TuneEE3C_Flat_7TeV_herwigpp"
else:
print "compareMCGentoMCUnfolded: wrong action", action, "skipping (usually you can ignore this message)"
return
# detaGen_central_jet15
fileWithUnfoldedPlotsName = optionsReg["odir"]+"/mnxsHistos_unfolded_"+action +".root"
fileWithUnfoldedPlots = ROOT.TFile(fileWithUnfoldedPlotsName)
#mnxsHistos_unfolded_pythiaOnHerwig.root
histos = getHistos(infileName)
#print histos[unfoldingWasDoneOn].keys()
todo = ["_jet15", "_dj15fb"]
#todo = ["_jet15"]
c = ROOT.TCanvas()
for t in todo:
genHisto = histos[unfoldingWasDoneOn]["detaGen_central"+t]
unfoldedHistoName = t+"/xsunfolded_central"+t
unfoldedHisto = fileWithUnfoldedPlots.Get(unfoldedHistoName)
#print unfoldedHistoName, type(unfoldedHisto), unfoldedHisto.ClassName()
#genHisto.Scale(0.5)
genHisto.Draw()
genHisto.GetXaxis().SetTitle(DrawMNPlots.xLabels()["xs"])
genHisto.GetYaxis().SetTitleOffset(1.8)
genHisto.GetYaxis().SetTitle(DrawMNPlots.yLabels()["xsAsPB"])
genHisto.SetMarkerColor(2)
genHisto.SetLineColor(2)
unfoldedHisto.Draw("SAME")
trueMax = max(genHisto.GetMaximum(), unfoldedHisto.GetMaximum())
genHisto.SetMaximum(trueMax*1.07)
c.Print(optionsReg["odir"]+"/MConMCunfoldingTest_"+action+t+".png")
c.Print(optionsReg["odir"]+"/MConMCunfoldingTest_"+action+t+".pdf")
