19def ReadEtaIntercalibrationHistograms(dirName):
20 if not dirName.endswith("/"):
21 dirName = dirName + "/"
22
23
24 histos = {}
25 files = sorted(glob.glob(dirName+"*.root"))
26 for aFileName in files:
27
28 jetDef = ""
29 for aFileDef,aJetDef in jetDefDict.iteritems():
30 if aFileDef in aFileName:
31 jetDef = aJetDef
32 break
33 if jetDef == "":
34 print "Failed to determine jet definition for file:",aFileName
35 return None
36 histos[jetDef] = {}
37
38
39 inFile = TFile(aFileName,"READ")
40 for aName in SystematicNames:
41 systematicName = aName + "_" + jetDef
42 histo = inFile.Get(systematicName)
43 if histo is None:
44 print "Failed to get histogram:",systematicName
45 return None
46 histo.SetName(systematicName+"_old")
47
48
50 histoNew.SetDirectory(0)
51 histos[jetDef][aName] = histoNew
52
53
54 inFile.Close()
55
56 return histos
ExtendPtRangeOfHisto(histo, histoName)