4 from array
import array
8 import ProviderHistoHelpers
10 SystematicNames = {
'EtaIntercalibration_TotalStat' :
'EtaIntercalibration_TotalStat',\
11 'EtaIntercalibration_TotalSyst' :
'EtaIntercalibration_Modelling'}
14 'AntiKt4EMTopo' :
'AntiKt4Topo_EMJES',
15 'AntiKt6EMTopo' :
'AntiKt6Topo_EMJES',
16 'AntiKt4LCTopo' :
'AntiKt4Topo_LCJES',
17 'AntiKt6LCTopo' :
'AntiKt6Topo_LCJES'
21 if not dirName.endswith(
"/"):
22 dirName = dirName +
"/"
26 files =
sorted(glob.glob(dirName+
"*.root"))
27 for aFileName
in files:
30 for aFileDef,aJetDef
in jetDefDict.iteritems():
31 if aFileDef
in aFileName:
35 print "Failed to determine jet definition for file:",aFileName
40 inFile = TFile(aFileName,
"READ")
41 for aName
in SystematicNames.keys():
42 getsystematicName = aName +
"_" + jetDef.replace(
"Topo",
"")
43 systematicName = SystematicNames[aName] +
"_" + jetDef
45 histo = inFile.Get(getsystematicName)
47 print "Failed to get histogram:",getsystematicName
49 histo.SetName(systematicName+
"_old")
53 histoNew.SetDirectory(0)
54 histos[jetDef][SystematicNames[aName]] = histoNew