def compare_hists(f1_name, f2_name, hist_name):
f1 = ROOT.TFile.Open(f1_name, "read") # Data
f2 = ROOT.TFile.Open(f2_name, "read") # MC
h1 = f1.Get(hist_name)
#if "v7" in f1_name:
#h2.Rebin(5)
if "solation" in f2_name:
if "_extend" in f2_name:
h2 = f2.Get("h_isoshape")
else:
h2 = f2.Get("up")
elif "irreducible" in hist_name:
h2 = f2.Get("hreluncert_irreducible")
else:
h2 = f2.Get(hist_name)
print "binning: ", h1.GetNbinsX(), h2.GetNbinsX()
print h2.GetName()
sys_name = hist_name.split('_')[1]
factor = 1.0
is_log = False
if "bkg_total_gg_full" in hist_name:
is_log = True
factor = h1.Integral()/h2.Integral()
print "intergral: ", h1.Integral(), h2.Integral(), factor
h2.Scale(factor)
h1.SetMarkerSize(0)
h2.SetMarkerSize(0)
ROOT.compare_two_hists(h1, h2,
#"Isolation [GeV]",
"m_{#gamma#gamma} [GeV]",
sys_name+" from ws",
sys_name+" from original input",
is_log)
f1.Close()
f2.Close()
def compare_hist_hist(f1_name, f2_name, hist_name):
f1 = ROOT.TFile.Open(f1_name, "read") # New
f2 = ROOT.TFile.Open(f2_name, "read") # Old
h1 = f1.Get(hist_name)
if "olation" in f2_name:
h2 = f2.Get("up")
else:
h2 = f2.Get(hist_name)
print "binning: ", h1.GetNbinsX(), h2.GetNbinsX()
#h2.Rebin(5)
h2_new = extend_hist(h2, h1)
h1.SetMarkerSize(0)
h2_new.SetMarkerSize(0)
ROOT.compare_two_hists(h1, h2_new, "m_{#gamma#gamma} [GeV]",
"New", "Old", True)
f1.Close()
f2.Close()
