4 from array
import array
8 import ProviderHistoHelpers
11 SystematicNameDictionary = {
12 'RelativeNonClosure_AFII_AntiKt4EMTopo' :
'RelativeNonClosure_AFII'
16 '4EM' :
'AntiKt4Topo_EMJES',
23 if dirName.endswith(
".root") :
26 elif not dirName.endswith(
"/") :
27 dirName = dirName +
"/"
28 files =
sorted(glob.glob(dirName+
"*.root"))
32 for aFileName
in files:
35 for aFileDef,aJetDef
in jetDefDict.iteritems():
36 if aFileDef
in aFileName:
40 print "Failed to determine jet definition for file:",aFileName
44 print "Jet def is",jetDef
47 inFile = TFile(aFileName,
"READ")
48 for aName,aSystName
in SystematicNameDictionary.iteritems():
49 systematicName = aSystName +
"_" + jetDef
50 histo = inFile.Get(aName)
52 print "Failed to get histogram:",systematicName
57 for ybin
in range(histo.GetNbinsY()+2) :
60 for xbin
in range(histo.GetNbinsX()+2) :
61 inbin = histo.GetBinContent(xbin,ybin)
67 histo.SetBinContent(xbin,ybin,freezeval)
71 histoSym.SetDirectory(0)
72 histos[jetDef][aSystName] = histoSym
78 histos[jetDef][
"RelativeNonClosure_MC15"].
SetDirectory(0)