def __init__(self):
self.data = {}
self.br = helpers.readCSVXS("etc/brHmm.txt")
self.vbf = {}
self.gg = {}
self.wh = {}
self.zh = {}
for e in ["8TeV", "7TeV"]:
self.vbf[e] = helpers.readCSVXS("etc/" + e + "-vbfH.txt")
self.gg[e] = helpers.readCSVXS("etc/" + e + "-ggH.txt")
self.wh[e] = helpers.readCSVXS("etc/" + e + "-WH.txt")
self.zh[e] = helpers.readCSVXS("etc/" + e + "-ZH.txt")
def __init__(self):
self.data = {}
self.br = helpers.readCSVXS("etc/br.csv")
self.vbf = {}
self.gg = {}
self.wh = {}
self.zh = {}
for e in ["8TeV", "7TeV"]:
self.vbf[e] = helpers.readCSVXS("etc/vbfH_" + e + ".csv")
self.gg[e] = helpers.readCSVXS("etc/ggH_" + e + ".csv")
self.wh[e] = helpers.readCSVXS("etc/wH_" + e + ".csv")
self.zh[e] = helpers.readCSVXS("etc/zH_" + e + ".csv")
def __init__(self):
self.data = {}
self.br = helpers.readCSVXS("etc/br.csv")
self.vbf = {}
self.gg = {}
self.wh = {}
self.zh = {}
for e in ["8TeV","7TeV"]:
self.vbf[e] = helpers.readCSVXS("etc/vbfH_"+e+".csv")
self.gg[e] = helpers.readCSVXS("etc/ggH_"+e+".csv")
self.wh[e] = helpers.readCSVXS("etc/wH_"+e+".csv")
self.zh[e] = helpers.readCSVXS("etc/zH_"+e+".csv")
def __init__(self):
self.data = {}
self.br = helpers.readCSVXS("etc/br.csv")
self.vbf = {}
self.gg = {}
self.wh = {}
self.zh = {}
for e in ["8TeV", "7TeV"]:
self.vbf[e] = helpers.readCSVXS("etc/vbfH_" + e + ".csv")
self.gg[e] = helpers.readCSVXS("etc/ggH_" + e + ".csv")
self.wh[e] = helpers.readCSVXS("etc/wH_" + e + ".csv")
self.zh[e] = helpers.readCSVXS("etc/zH_" + e + ".csv")
self.lumi = {
#"14TeV" : 1.044, # 2012 HF-Lumi
#"8TeV" : 1.044, # 2012 HF-Lumi
"14TeV": 1.026, # 2012 Pixel-Lumi
"8TeV": 1.026, # 2012 Pixel-Lumi
"7TeV": 1.022,
}
self.JES = {
'gg': {
'7TeV': {
'Jet2CutsFailVBFGF': 1.0832,
'Jet2CutsGFPass': 1.0585,
'Jet2CutsVBFPass': 1.0759,
'Jets01FailPtG10': None,
'Jets01PassPtG10': None,
},
'8TeV': {
'Jet2CutsFailVBFGF': 1.0799,
'Jet2CutsGFPass': 1.0477,
'Jet2CutsVBFPass': 1.0606,
'Jets01FailPtG10': None,
'Jets01PassPtG10': None,
},
},
'vbf': {
'7TeV': {
'Jet2CutsFailVBFGF': 1.0572,
'Jet2CutsGFPass': 1.0381,
'Jet2CutsVBFPass': 1.0358,
'Jets01FailPtG10': -1.0321,
'Jets01PassPtG10': -1.0308,
},
'8TeV': {
'Jet2CutsFailVBFGF': 1.0482,
'Jet2CutsGFPass': 1.0249,
'Jet2CutsVBFPass': 1.0185,
'Jets01FailPtG10': -1.0278,
'Jets01PassPtG10': -1.0250,
},
},
'wh': {
'7TeV': {
'Jet2CutsFailVBFGF': None,
'Jet2CutsGFPass': None,
'Jet2CutsVBFPass': None,
'Jets01FailPtG10': None,
'Jets01PassPtG10': None,
},
'8TeV': {
'Jet2CutsFailVBFGF': None,
'Jet2CutsGFPass': None,
'Jet2CutsVBFPass': None,
'Jets01FailPtG10': None,
'Jets01PassPtG10': None,
},
},
'zh': {
'7TeV': {
'Jet2CutsFailVBFGF': None,
'Jet2CutsGFPass': None,
'Jet2CutsVBFPass': None,
'Jets01FailPtG10': None,
'Jets01PassPtG10': None,
},
'8TeV': {
'Jet2CutsFailVBFGF': None,
'Jet2CutsGFPass': None,
'Jet2CutsVBFPass': None,
'Jets01FailPtG10': None,
'Jets01PassPtG10': None,
},
},
}
self.JER = {
'gg': {
'7TeV': {
'Jet2CutsFailVBFGF': -1.0166,
'Jet2CutsGFPass': -1.0126,
'Jet2CutsVBFPass': 1.0321,
'Jets01FailPtG10': None,
'Jets01PassPtG10': None,
},
'8TeV': {
'Jet2CutsFailVBFGF': -1.0131,
'Jet2CutsGFPass': -1.0101,
'Jet2CutsVBFPass': -1.0273,
'Jets01FailPtG10': None,
'Jets01PassPtG10': None,
},
},
'vbf': {
'7TeV': {
'Jet2CutsFailVBFGF': -1.0111,
'Jet2CutsGFPass': None,
'Jet2CutsVBFPass': None,
'Jets01FailPtG10': -1.0167,
'Jets01PassPtG10': None,
},
'8TeV': {
'Jet2CutsFailVBFGF': None,
'Jet2CutsGFPass': None,
'Jet2CutsVBFPass': None,
'Jets01FailPtG10': -1.0162,
'Jets01PassPtG10': None,
},
},
'wh': {
'7TeV': {
'Jet2CutsFailVBFGF': None,
'Jet2CutsGFPass': None,
'Jet2CutsVBFPass': None,
'Jets01FailPtG10': None,
'Jets01PassPtG10': None,
},
'8TeV': {
'Jet2CutsFailVBFGF': None,
'Jet2CutsGFPass': None,
'Jet2CutsVBFPass': None,
'Jets01FailPtG10': None,
'Jets01PassPtG10': None,
},
},
'zh': {
'7TeV': {
'Jet2CutsFailVBFGF': None,
'Jet2CutsGFPass': None,
'Jet2CutsVBFPass': None,
'Jets01FailPtG10': None,
'Jets01PassPtG10': None,
},
'8TeV': {
'Jet2CutsFailVBFGF': None,
'Jet2CutsGFPass': None,
'Jet2CutsVBFPass': None,
'Jets01FailPtG10': None,
'Jets01PassPtG10': None,
},
},
}
self.PUID = {
'gg': {
'7TeV': {
'Jet2CutsFailVBFGF': 1.0109,
'Jet2CutsGFPass': 1.0158,
'Jet2CutsVBFPass': 1.0368,
'Jets01FailPtG10': None,
'Jets01PassPtG10': None,
},
'8TeV': {
'Jet2CutsFailVBFGF': 1.0142,
'Jet2CutsGFPass': 1.0174,
'Jet2CutsVBFPass': 1.0386,
'Jets01FailPtG10': None,
'Jets01PassPtG10': None,
},
},
'vbf': {
'7TeV': {
'Jet2CutsFailVBFGF': 1.0134,
'Jet2CutsGFPass': 1.0160,
'Jet2CutsVBFPass': 1.0324,
'Jets01FailPtG10': 1.0115,
'Jets01PassPtG10': 1.0108,
},
'8TeV': {
'Jet2CutsFailVBFGF': 1.0155,
'Jet2CutsGFPass': 1.0168,
'Jet2CutsVBFPass': 1.0328,
'Jets01FailPtG10': 1.0129,
'Jets01PassPtG10': 1.0128,
},
},
'wh': {
'7TeV': {
'Jet2CutsFailVBFGF': None,
'Jet2CutsGFPass': None,
'Jet2CutsVBFPass': None,
'Jets01FailPtG10': None,
'Jets01PassPtG10': None,
},
'8TeV': {
'Jet2CutsFailVBFGF': None,
'Jet2CutsGFPass': None,
'Jet2CutsVBFPass': None,
'Jets01FailPtG10': None,
'Jets01PassPtG10': None,
},
},
'zh': {
'7TeV': {
'Jet2CutsFailVBFGF': None,
'Jet2CutsGFPass': None,
'Jet2CutsVBFPass': None,
'Jets01FailPtG10': None,
'Jets01PassPtG10': None,
},
'8TeV': {
'Jet2CutsFailVBFGF': None,
'Jet2CutsGFPass': None,
'Jet2CutsVBFPass': None,
'Jets01FailPtG10': None,
'Jets01PassPtG10': None,
},
},
}
self.MCStat = {
'gg': {
'7TeV': {
'Jet2CutsFailVBFGF': 1.0301,
'Jet2CutsGFPass': 1.0572,
'Jet2CutsVBFPass': 1.1628,
'Jets01FailPtG10': 1.0154,
'Jets01PassPtG10': None,
},
'8TeV': {
'Jet2CutsFailVBFGF': 1.0308,
'Jet2CutsGFPass': 1.0577,
'Jet2CutsVBFPass': 1.1068,
'Jets01FailPtG10': 1.0147,
'Jets01PassPtG10': None,
},
},
'vbf': {
'7TeV': {
'Jet2CutsFailVBFGF': 1.0218,
'Jet2CutsGFPass': 1.0177,
'Jet2CutsVBFPass': 1.0208,
'Jets01FailPtG10': 1.0628,
'Jets01PassPtG10': None,
},
'8TeV': {
'Jet2CutsFailVBFGF': 1.0178,
'Jet2CutsGFPass': 1.0152,
'Jet2CutsVBFPass': 1.0135,
'Jets01FailPtG10': 1.0642,
'Jets01PassPtG10': None,
},
},
'wh': {
'7TeV': {
'Jet2CutsFailVBFGF': None,
'Jet2CutsGFPass': None,
'Jet2CutsVBFPass': None,
'Jets01FailPtG10': None,
'Jets01PassPtG10': None,
},
'8TeV': {
'Jet2CutsFailVBFGF': None,
'Jet2CutsGFPass': None,
'Jet2CutsVBFPass': None,
'Jets01FailPtG10': None,
'Jets01PassPtG10': None,
},
},
'zh': {
'7TeV': {
'Jet2CutsFailVBFGF': None,
'Jet2CutsGFPass': None,
'Jet2CutsVBFPass': None,
'Jets01FailPtG10': None,
'Jets01PassPtG10': None,
},
'8TeV': {
'Jet2CutsFailVBFGF': None,
'Jet2CutsGFPass': None,
'Jet2CutsVBFPass': None,
'Jets01FailPtG10': None,
'Jets01PassPtG10': None,
},
},
}
self.PDF = {
'gg': {
'8TeV': {
'Jets01PassPtG10': 1.0799,
'Jets01FailPtG10': 1.0844,
'Jet2CutsVBFPass': 1.0578,
'Jet2CutsGFPass': 1.0715,
'Jet2CutsFailVBFGF': 1.0809,
},
'7TeV': {
'Jets01PassPtG10': 1.0830,
'Jets01FailPtG10': 1.0868,
'Jet2CutsVBFPass': 1.0593,
'Jet2CutsGFPass': 1.0743,
'Jet2CutsFailVBFGF': 1.0801,
},
},
'vbf': {
'8TeV': {
'Jets01PassPtG10': 1.0197,
'Jets01FailPtG10': 1.0186,
'Jet2CutsVBFPass': 1.0378,
'Jet2CutsGFPass': 1.0212,
'Jet2CutsFailVBFGF': 1.0196,
},
'7TeV': {
'Jets01PassPtG10': 1.0220,
'Jets01FailPtG10': 1.0214,
'Jet2CutsVBFPass': 1.0448,
'Jet2CutsGFPass': 1.0299,
'Jet2CutsFailVBFGF': 1.0204,
},
},
'wh': {
'8TeV': {
'Jets01PassPtG10': None,
'Jets01FailPtG10': None,
'Jet2CutsVBFPass': None,
'Jet2CutsGFPass': None,
'Jet2CutsFailVBFGF': None,
},
'7TeV': {
'Jets01PassPtG10': None,
'Jets01FailPtG10': None,
'Jet2CutsVBFPass': None,
'Jet2CutsGFPass': None,
'Jet2CutsFailVBFGF': None,
},
},
'zh': {
'8TeV': {
'Jets01PassPtG10': None,
'Jets01FailPtG10': None,
'Jet2CutsVBFPass': None,
'Jet2CutsGFPass': None,
'Jet2CutsFailVBFGF': None,
},
'7TeV': {
'Jets01PassPtG10': None,
'Jets01FailPtG10': None,
'Jet2CutsVBFPass': None,
'Jet2CutsGFPass': None,
'Jet2CutsFailVBFGF': None,
},
},
}
self.PU = {
'gg': {
'8TeV': {
'Jets01PassPtG10': None,
'Jets01FailPtG10': None,
'Jet2CutsVBFPass': 1.0222,
'Jet2CutsGFPass': 1.0103,
'Jet2CutsFailVBFGF': 1.0112,
},
'7TeV': {
'Jets01PassPtG10': None,
'Jets01FailPtG10': None,
'Jet2CutsVBFPass': 1.0222,
'Jet2CutsGFPass': 1.0103,
'Jet2CutsFailVBFGF': 1.0112,
},
},
'vbf': {
'8TeV': {
'Jets01PassPtG10': None,
'Jets01FailPtG10': 1.0207,
'Jet2CutsVBFPass': None,
'Jet2CutsGFPass': None,
'Jet2CutsFailVBFGF': None,
},
'7TeV': {
'Jets01PassPtG10': None,
'Jets01FailPtG10': 1.0207,
'Jet2CutsVBFPass': None,
'Jet2CutsGFPass': None,
'Jet2CutsFailVBFGF': None,
},
},
'wh': {
'8TeV': {
'Jets01PassPtG10': None,
'Jets01FailPtG10': None,
'Jet2CutsVBFPass': None,
'Jet2CutsGFPass': None,
'Jet2CutsFailVBFGF': None,
},
'7TeV': {
'Jets01PassPtG10': None,
'Jets01FailPtG10': None,
'Jet2CutsVBFPass': None,
'Jet2CutsGFPass': None,
'Jet2CutsFailVBFGF': None,
},
},
'zh': {
'8TeV': {
'Jets01PassPtG10': None,
'Jets01FailPtG10': None,
'Jet2CutsVBFPass': None,
'Jet2CutsGFPass': None,
'Jet2CutsFailVBFGF': None,
},
'7TeV': {
'Jets01PassPtG10': None,
'Jets01FailPtG10': None,
'Jet2CutsVBFPass': None,
'Jet2CutsGFPass': None,
'Jet2CutsFailVBFGF': None,
},
},
}
self.QCDScale = {
'gg': {
'8TeV': {
'Jets01PassPtG10': -1.0121,
'Jets01FailPtG10': 1.0450,
'Jet2CutsVBFPass': -1.1585,
'Jet2CutsGFPass': 1.1515,
'Jet2CutsFailVBFGF': -1.0662,
},
'7TeV': {
'Jets01PassPtG10': 1.0104,
'Jets01FailPtG10': 1.0417,
'Jet2CutsVBFPass': -1.2424,
'Jet2CutsGFPass': -1.1608,
'Jet2CutsFailVBFGF': 1.0626,
},
},
'vbf': {
'8TeV': {
'Jets01PassPtG10': -1.0139,
'Jets01FailPtG10': -1.1598,
'Jet2CutsVBFPass': 1.0313,
'Jet2CutsGFPass': -1.0372,
'Jet2CutsFailVBFGF': -1.0427,
},
'7TeV': {
'Jets01PassPtG10': 1.0101,
'Jets01FailPtG10': 1.0571,
'Jet2CutsVBFPass': -1.0489,
'Jet2CutsGFPass': -1.0380,
'Jet2CutsFailVBFGF': 1.0653,
},
},
'wh': {
'8TeV': {
'Jets01PassPtG10': None,
'Jets01FailPtG10': None,
'Jet2CutsVBFPass': None,
'Jet2CutsGFPass': None,
'Jet2CutsFailVBFGF': None,
},
'7TeV': {
'Jets01PassPtG10': None,
'Jets01FailPtG10': None,
'Jet2CutsVBFPass': None,
'Jet2CutsGFPass': None,
'Jet2CutsFailVBFGF': None,
},
},
'zh': {
'8TeV': {
'Jets01PassPtG10': None,
'Jets01FailPtG10': None,
'Jet2CutsVBFPass': None,
'Jet2CutsGFPass': None,
'Jet2CutsFailVBFGF': None,
},
'7TeV': {
'Jets01PassPtG10': None,
'Jets01FailPtG10': None,
'Jet2CutsVBFPass': None,
'Jet2CutsGFPass': None,
'Jet2CutsFailVBFGF': None,
},
},
}
self.UE = {
'gg': {
'8TeV': {
'Jets01PassPtG10': -1.0618,
'Jets01FailPtG10': 1.1698,
'Jet2CutsVBFPass': 1.4228,
'Jet2CutsGFPass': 1.2317,
'Jet2CutsFailVBFGF': -1.0847,
},
'7TeV': {
'Jets01PassPtG10': -1.0658,
'Jets01FailPtG10': 1.1715,
'Jet2CutsVBFPass': 1.6606,
'Jet2CutsGFPass': 1.2570,
'Jet2CutsFailVBFGF': -1.0669,
},
},
'vbf': {
'8TeV': {
'Jets01PassPtG10': -1.0309,
'Jets01FailPtG10': -1.0794,
'Jet2CutsVBFPass': -1.1030,
'Jet2CutsGFPass': -1.0585,
'Jet2CutsFailVBFGF': 1.0432,
},
'7TeV': {
'Jets01PassPtG10': -1.0515,
'Jets01FailPtG10': 1.0953,
'Jet2CutsVBFPass': -1.1421,
'Jet2CutsGFPass': -1.1014,
'Jet2CutsFailVBFGF': 1.0763,
},
},
'wh': {
'8TeV': {
'Jets01PassPtG10': None,
'Jets01FailPtG10': None,
'Jet2CutsVBFPass': None,
'Jet2CutsGFPass': None,
'Jet2CutsFailVBFGF': None,
},
'7TeV': {
'Jets01PassPtG10': None,
'Jets01FailPtG10': None,
'Jet2CutsVBFPass': None,
'Jet2CutsGFPass': None,
'Jet2CutsFailVBFGF': None,
},
},
'zh': {
'8TeV': {
'Jets01PassPtG10': None,
'Jets01FailPtG10': None,
'Jet2CutsVBFPass': None,
'Jet2CutsGFPass': None,
'Jet2CutsFailVBFGF': None,
},
'7TeV': {
'Jets01PassPtG10': None,
'Jets01FailPtG10': None,
'Jet2CutsVBFPass': None,
'Jet2CutsGFPass': None,
'Jet2CutsFailVBFGF': None,
},
},
}
# The list of systematics to be applied
self._keys = [
"xs_ggH", "xs_vbfH", "xs_whH", "xs_zhH", "br_Hmm", "lumi", "PDF",
"JES", "JER", "PUID", "MCStat", "PU", "UE", "QCDScale"
]
# The list of systematics which are correlated between energies
self.keysEnergyCorr = [
"xs_ggH", "xs_vbfH", "xs_whH", "xs_zhH", "br_Hmm", "PDF", "JER",
"UE", "QCDScale"
]
# The list of systematics which are not correlated between energies or categories
self.keysNotCatCorr = ["MCStat"]
# The list of systematics which are not correlated between energies, but correlated w/ categories
# (everything else)
self.keysNotEnergyCorr = []
for i in self._keys:
if not i in self.keysEnergyCorr and not i in self.keysNotCatCorr:
self.keysNotEnergyCorr.append(i)
import helpers
lumiDict = {}
#lumiDict["8TeV"] = 19.79 #2012ABCD 22Jan2013 HF-Lumi
lumiDict["8TeV"] = 19.712 #2012ABCD 22Jan2013 Pixel-Lumi
lumiDict["7TeV"] = 5.05 #2011AB
#LUMI=0.807 #2012A
#LUMI=4.421 #2012B
#LUMI=0.495 #2012Cv1
#LUMI=6.311 #2012Cv2
#LUMI=2.311 #2011A
#LUMI=2.739 #2011B
brDict = helpers.readCSVXS("etc/br.csv")
ggHDict8 = helpers.readCSVXS("etc/ggH_8TeV.csv")
vbfHDict8 = helpers.readCSVXS("etc/vbfH_8TeV.csv")
wHDict8 = helpers.readCSVXS("etc/wH_8TeV.csv")
zHDict8 = helpers.readCSVXS("etc/zH_8TeV.csv")
ggHDict7 = helpers.readCSVXS("etc/ggH_7TeV.csv")
vbfHDict7 = helpers.readCSVXS("etc/vbfH_7TeV.csv")
wHDict7 = helpers.readCSVXS("etc/wH_7TeV.csv")
zHDict7 = helpers.readCSVXS("etc/zH_7TeV.csv")
class CrossSections:
def __init__(self):
self.data = {}
self.br = helpers.readCSVXS("etc/br.csv")
def __init__(self):
self.data = {}
self.br = helpers.readCSVXS("etc/brHmm.txt")
self.vbf = {}
self.gg = {}
self.wh = {}
self.zh = {}
for e in ["8TeV", "7TeV"]:
self.vbf[e] = helpers.readCSVXS("etc/" + e + "-vbfH.txt")
self.gg[e] = helpers.readCSVXS("etc/" + e + "-ggH.txt")
self.wh[e] = helpers.readCSVXS("etc/" + e + "-WH.txt")
self.zh[e] = helpers.readCSVXS("etc/" + e + "-ZH.txt")
self.lumi = {
# "14TeV" : 1.044, # 2012 HF-Lumi
# "8TeV" : 1.044, # 2012 HF-Lumi
"14TeV": 1.026, # 2012 Pixel-Lumi
"8TeV": 1.026, # 2012 Pixel-Lumi
"7TeV": 1.022,
}
self.CMS_scale_j = {
"gg": {
"7TeV": {
"Jet2CutsFailVBFGF": 1.0832,
"Jet2CutsGFPass": 1.0585,
"Jet2CutsVBFPass": 1.0759,
"Jets01FailPtG10": None,
"Jets01PassPtG10": None,
},
"8TeV": {
"Jet2CutsFailVBFGF": 1.0799,
"Jet2CutsGFPass": 1.0477,
"Jet2CutsVBFPass": 1.0606,
"Jets01FailPtG10": None,
"Jets01PassPtG10": None,
},
},
"vbf": {
"7TeV": {
"Jet2CutsFailVBFGF": 1.0572,
"Jet2CutsGFPass": 1.0381,
"Jet2CutsVBFPass": 1.0358,
"Jets01FailPtG10": -1.0321,
"Jets01PassPtG10": -1.0308,
},
"8TeV": {
"Jet2CutsFailVBFGF": 1.0482,
"Jet2CutsGFPass": 1.0249,
"Jet2CutsVBFPass": 1.0185,
"Jets01FailPtG10": -1.0278,
"Jets01PassPtG10": -1.0250,
},
},
"wh": {
"7TeV": {
"Jet2CutsFailVBFGF": None,
"Jet2CutsGFPass": None,
"Jet2CutsVBFPass": None,
"Jets01FailPtG10": None,
"Jets01PassPtG10": None,
},
"8TeV": {
"Jet2CutsFailVBFGF": None,
"Jet2CutsGFPass": None,
"Jet2CutsVBFPass": None,
"Jets01FailPtG10": None,
"Jets01PassPtG10": None,
},
},
"zh": {
"7TeV": {
"Jet2CutsFailVBFGF": None,
"Jet2CutsGFPass": None,
"Jet2CutsVBFPass": None,
"Jets01FailPtG10": None,
"Jets01PassPtG10": None,
},
"8TeV": {
"Jet2CutsFailVBFGF": None,
"Jet2CutsGFPass": None,
"Jet2CutsVBFPass": None,
"Jets01FailPtG10": None,
"Jets01PassPtG10": None,
},
},
}
self.CMS_res_j = {
"gg": {
"7TeV": {
"Jet2CutsFailVBFGF": -1.0166,
"Jet2CutsGFPass": -1.0126,
"Jet2CutsVBFPass": 1.0321,
"Jets01FailPtG10": None,
"Jets01PassPtG10": None,
},
"8TeV": {
"Jet2CutsFailVBFGF": -1.0131,
"Jet2CutsGFPass": -1.0101,
"Jet2CutsVBFPass": -1.0273,
"Jets01FailPtG10": None,
"Jets01PassPtG10": None,
},
},
"vbf": {
"7TeV": {
"Jet2CutsFailVBFGF": -1.0111,
"Jet2CutsGFPass": None,
"Jet2CutsVBFPass": None,
"Jets01FailPtG10": -1.0167,
"Jets01PassPtG10": None,
},
"8TeV": {
"Jet2CutsFailVBFGF": None,
"Jet2CutsGFPass": None,
"Jet2CutsVBFPass": None,
"Jets01FailPtG10": -1.0162,
"Jets01PassPtG10": None,
},
},
"wh": {
"7TeV": {
"Jet2CutsFailVBFGF": None,
"Jet2CutsGFPass": None,
"Jet2CutsVBFPass": None,
"Jets01FailPtG10": None,
"Jets01PassPtG10": None,
},
"8TeV": {
"Jet2CutsFailVBFGF": None,
"Jet2CutsGFPass": None,
"Jet2CutsVBFPass": None,
"Jets01FailPtG10": None,
"Jets01PassPtG10": None,
},
},
"zh": {
"7TeV": {
"Jet2CutsFailVBFGF": None,
"Jet2CutsGFPass": None,
"Jet2CutsVBFPass": None,
"Jets01FailPtG10": None,
"Jets01PassPtG10": None,
},
"8TeV": {
"Jet2CutsFailVBFGF": None,
"Jet2CutsGFPass": None,
"Jet2CutsVBFPass": None,
"Jets01FailPtG10": None,
"Jets01PassPtG10": None,
},
},
}
self.MCStat = {
"gg": {
"7TeV": {
"Jet2CutsFailVBFGF": 1.0301,
"Jet2CutsGFPass": 1.0572,
"Jet2CutsVBFPass": 1.1628,
"Jets01FailPtG10": 1.0154,
"Jets01PassPtG10": None,
},
"8TeV": {
"Jet2CutsFailVBFGF": 1.0308,
"Jet2CutsGFPass": 1.0577,
"Jet2CutsVBFPass": 1.1068,
"Jets01FailPtG10": 1.0147,
"Jets01PassPtG10": None,
},
},
"vbf": {
"7TeV": {
"Jet2CutsFailVBFGF": 1.0218,
"Jet2CutsGFPass": 1.0177,
"Jet2CutsVBFPass": 1.0208,
"Jets01FailPtG10": 1.0628,
"Jets01PassPtG10": None,
},
"8TeV": {
"Jet2CutsFailVBFGF": 1.0178,
"Jet2CutsGFPass": 1.0152,
"Jet2CutsVBFPass": 1.0135,
"Jets01FailPtG10": 1.0642,
"Jets01PassPtG10": None,
},
},
"wh": {
"7TeV": {
"Jet2CutsFailVBFGF": None,
"Jet2CutsGFPass": None,
"Jet2CutsVBFPass": None,
"Jets01FailPtG10": None,
"Jets01PassPtG10": None,
},
"8TeV": {
"Jet2CutsFailVBFGF": None,
"Jet2CutsGFPass": None,
"Jet2CutsVBFPass": None,
"Jets01FailPtG10": None,
"Jets01PassPtG10": None,
},
},
"zh": {
"7TeV": {
"Jet2CutsFailVBFGF": None,
"Jet2CutsGFPass": None,
"Jet2CutsVBFPass": None,
"Jets01FailPtG10": None,
"Jets01PassPtG10": None,
},
"8TeV": {
"Jet2CutsFailVBFGF": None,
"Jet2CutsGFPass": None,
"Jet2CutsVBFPass": None,
"Jets01FailPtG10": None,
"Jets01PassPtG10": None,
},
},
}
self.PDF = {
"gg": {
"8TeV": {
"Jets01PassPtG10": 1.0799,
"Jets01FailPtG10": 1.0844,
"Jet2CutsVBFPass": 1.0578,
"Jet2CutsGFPass": 1.0715,
"Jet2CutsFailVBFGF": 1.0809,
},
"7TeV": {
"Jets01PassPtG10": 1.0830,
"Jets01FailPtG10": 1.0868,
"Jet2CutsVBFPass": 1.0593,
"Jet2CutsGFPass": 1.0743,
"Jet2CutsFailVBFGF": 1.0801,
},
},
"vbf": {
"8TeV": {
"Jets01PassPtG10": 1.0197,
"Jets01FailPtG10": 1.0186,
"Jet2CutsVBFPass": 1.0378,
"Jet2CutsGFPass": 1.0212,
"Jet2CutsFailVBFGF": 1.0196,
},
"7TeV": {
"Jets01PassPtG10": 1.0220,
"Jets01FailPtG10": 1.0214,
"Jet2CutsVBFPass": 1.0448,
"Jet2CutsGFPass": 1.0299,
"Jet2CutsFailVBFGF": 1.0204,
},
},
"wh": {
"8TeV": {
"Jets01PassPtG10": None,
"Jets01FailPtG10": None,
"Jet2CutsVBFPass": None,
"Jet2CutsGFPass": None,
"Jet2CutsFailVBFGF": None,
},
"7TeV": {
"Jets01PassPtG10": None,
"Jets01FailPtG10": None,
"Jet2CutsVBFPass": None,
"Jet2CutsGFPass": None,
"Jet2CutsFailVBFGF": None,
},
},
"zh": {
"8TeV": {
"Jets01PassPtG10": None,
"Jets01FailPtG10": None,
"Jet2CutsVBFPass": None,
"Jet2CutsGFPass": None,
"Jet2CutsFailVBFGF": None,
},
"7TeV": {
"Jets01PassPtG10": None,
"Jets01FailPtG10": None,
"Jet2CutsVBFPass": None,
"Jet2CutsGFPass": None,
"Jet2CutsFailVBFGF": None,
},
},
}
self.QCDScale = {
"gg": {
"8TeV": {
"Jets01PassPtG10": -1.0121,
"Jets01FailPtG10": 1.0450,
"Jet2CutsVBFPass": -1.1585,
"Jet2CutsGFPass": 1.1515,
"Jet2CutsFailVBFGF": -1.0662,
},
"7TeV": {
"Jets01PassPtG10": 1.0104,
"Jets01FailPtG10": 1.0417,
"Jet2CutsVBFPass": -1.2424,
"Jet2CutsGFPass": -1.1608,
"Jet2CutsFailVBFGF": 1.0626,
},
},
"vbf": {
"8TeV": {
"Jets01PassPtG10": -1.0139,
"Jets01FailPtG10": -1.1598,
"Jet2CutsVBFPass": 1.0313,
"Jet2CutsGFPass": -1.0372,
"Jet2CutsFailVBFGF": -1.0427,
},
"7TeV": {
"Jets01PassPtG10": 1.0101,
"Jets01FailPtG10": 1.0571,
"Jet2CutsVBFPass": -1.0489,
"Jet2CutsGFPass": -1.0380,
"Jet2CutsFailVBFGF": 1.0653,
},
},
"wh": {
"8TeV": {
"Jets01PassPtG10": None,
"Jets01FailPtG10": None,
"Jet2CutsVBFPass": None,
"Jet2CutsGFPass": None,
"Jet2CutsFailVBFGF": None,
},
"7TeV": {
"Jets01PassPtG10": None,
"Jets01FailPtG10": None,
"Jet2CutsVBFPass": None,
"Jet2CutsGFPass": None,
"Jet2CutsFailVBFGF": None,
},
},
"zh": {
"8TeV": {
"Jets01PassPtG10": None,
"Jets01FailPtG10": None,
"Jet2CutsVBFPass": None,
"Jet2CutsGFPass": None,
"Jet2CutsFailVBFGF": None,
},
"7TeV": {
"Jets01PassPtG10": None,
"Jets01FailPtG10": None,
"Jet2CutsVBFPass": None,
"Jet2CutsGFPass": None,
"Jet2CutsFailVBFGF": None,
},
},
}
self.UEPS = {
"gg": {
"8TeV": {
"Jets01PassPtG10": -1.0618,
"Jets01FailPtG10": 1.1698,
"Jet2CutsVBFPass": 1.4228,
"Jet2CutsGFPass": 1.2317,
"Jet2CutsFailVBFGF": -1.0847,
},
"7TeV": {
"Jets01PassPtG10": -1.0658,
"Jets01FailPtG10": 1.1715,
"Jet2CutsVBFPass": 1.6606,
"Jet2CutsGFPass": 1.2570,
"Jet2CutsFailVBFGF": -1.0669,
},
},
"vbf": {
"8TeV": {
"Jets01PassPtG10": -1.0309,
"Jets01FailPtG10": -1.0794,
"Jet2CutsVBFPass": -1.1030,
"Jet2CutsGFPass": -1.0585,
"Jet2CutsFailVBFGF": 1.0432,
},
"7TeV": {
"Jets01PassPtG10": -1.0515,
"Jets01FailPtG10": 1.0953,
"Jet2CutsVBFPass": -1.1421,
"Jet2CutsGFPass": -1.1014,
"Jet2CutsFailVBFGF": 1.0763,
},
},
"wh": {
"8TeV": {
"Jets01PassPtG10": None,
"Jets01FailPtG10": None,
"Jet2CutsVBFPass": None,
"Jet2CutsGFPass": None,
"Jet2CutsFailVBFGF": None,
},
"7TeV": {
"Jets01PassPtG10": None,
"Jets01FailPtG10": None,
"Jet2CutsVBFPass": None,
"Jet2CutsGFPass": None,
"Jet2CutsFailVBFGF": None,
},
},
"zh": {
"8TeV": {
"Jets01PassPtG10": None,
"Jets01FailPtG10": None,
"Jet2CutsVBFPass": None,
"Jet2CutsGFPass": None,
"Jet2CutsFailVBFGF": None,
},
"7TeV": {
"Jets01PassPtG10": None,
"Jets01FailPtG10": None,
"Jet2CutsVBFPass": None,
"Jet2CutsGFPass": None,
"Jet2CutsFailVBFGF": None,
},
},
}
self.CMS_eff_j = {
"gg": {
"8TeV": {
"Jets01PassPtG10": None,
"Jets01FailPtG10": None,
"Jet2CutsFailVBFGF": 1.0181,
"Jet2CutsGFPass": 1.0202,
"Jet2CutsVBFPass": 1.0445,
},
"7TeV": {
"Jets01PassPtG10": None,
"Jets01FailPtG10": None,
"Jet2CutsFailVBFGF": 1.0156,
"Jet2CutsGFPass": 1.0189,
"Jet2CutsVBFPass": 1.0430,
},
},
"vbf": {
"8TeV": {
"Jets01PassPtG10": 1.0128,
"Jets01FailPtG10": 1.0244,
"Jet2CutsFailVBFGF": 1.0155,
"Jet2CutsGFPass": 1.0168,
"Jet2CutsVBFPass": 1.0328,
},
"7TeV": {
"Jets01PassPtG10": 1.0108,
"Jets01FailPtG10": 1.0237,
"Jet2CutsFailVBFGF": 1.0134,
"Jet2CutsGFPass": 1.0160,
"Jet2CutsVBFPass": 1.0324,
},
},
"wh": {
"8TeV": {
"Jets01PassPtG10": None,
"Jets01FailPtG10": None,
"Jet2CutsFailVBFGF": None,
"Jet2CutsGFPass": None,
"Jet2CutsVBFPass": None,
},
"7TeV": {
"Jets01PassPtG10": None,
"Jets01FailPtG10": None,
"Jet2CutsFailVBFGF": None,
"Jet2CutsGFPass": None,
"Jet2CutsVBFPass": None,
},
},
"zh": {
"8TeV": {
"Jets01PassPtG10": None,
"Jets01FailPtG10": None,
"Jet2CutsFailVBFGF": None,
"Jet2CutsGFPass": None,
"Jet2CutsVBFPass": None,
},
"7TeV": {
"Jets01PassPtG10": None,
"Jets01FailPtG10": None,
"Jet2CutsFailVBFGF": None,
"Jet2CutsGFPass": None,
"Jet2CutsVBFPass": None,
},
},
}
# Muon efficiency uncertainty on yield; same for everything (not known for VH)
self.CMS_eff_m = 1.016
# The list of systematics to be applied
self._keys = [
"QCDscale_ggH",
"QCDscale_qqH",
"QCDscale_VH",
"pdf_gg",
"pdf_qqbar",
"lumi",
"pdf_gg_ACCEPT",
"pdf_qqbar_ACCEPT",
"CMS_scale_j",
"CMS_res_j",
"MCStat",
"UEPS",
"QCDscale_ggH_ACCEPT",
"QCDscale_qqH_ACCEPT",
"CMS_eff_m",
"CMS_eff_j",
"br_Hmm",
]
# The list of systematics which are correlated between energies
self.keysEnergyCorr = [
"QCDscale_ggH",
"QCDscale_qqH",
"QCDscale_VH",
"pdf_gg",
"pdf_qqbar",
"pdf_gg_ACCEPT",
"pdf_qqbar_ACCEPT",
"CMS_scale_j",
"CMS_res_j",
"UEPS",
"QCDscale_ggH_ACCEPT",
"QCDscale_qqH_ACCEPT",
"CMS_eff_m",
"CMS_eff_j",
"br_Hmm",
]
# The list of systematics which are not correlated between energies or categories
self.keysNotCatCorr = ["MCStat"]
# The list of systematics which are not correlated between energies, but correlated w/ categories
# (everything else)
self.keysNotEnergyCorr = []
for i in self._keys:
if not i in self.keysEnergyCorr and not i in self.keysNotCatCorr:
self.keysNotEnergyCorr.append(i)
self.titleMap = {
"QCDscale_ggH": "QCD Scale",
"QCDscale_qqH": "QCD Scale",
"QCDscale_VH": "QCD Scale",
"pdf_gg": "PDF",
"pdf_qqbar": "PDF",
"lumi": "Luminosity",
"pdf_gg_ACCEPT": "PDF",
"pdf_qqbar_ACCEPT": "PDF",
"CMS_scale_j": "Jet Energy Scale",
"CMS_res_j": "Jet Energy Resolution",
"MCStat": "MC Statistics",
"UEPS": "UE/PS",
"QCDscale_ggH_ACCEPT": "QCD Scale",
"QCDscale_qqH_ACCEPT": "QCD Scale",
"CMS_eff_m": "Muon Efficiency",
"CMS_eff_j": "Jet Efficiency",
"br_Hmm": r"\BF(\hmm{})",
}
lumiDict = {}
# lumiDict["8TeV"] = 19.79 #2012ABCD 22Jan2013 HF-Lumi
lumiDict["8TeV"] = 19.712 # 2012ABCD 22Jan2013 Pixel-Lumi
lumiDict["7TeV"] = 5.05 # 2011AB
# LUMI=0.807 #2012A
# LUMI=4.421 #2012B
# LUMI=0.495 #2012Cv1
# LUMI=6.311 #2012Cv2
# LUMI=2.311 #2011A
# LUMI=2.739 #2011B
brDict = helpers.readCSVXS("etc/brHmm.txt")
ggHDict8 = helpers.readCSVXS("etc/8TeV-ggH.txt")
vbfHDict8 = helpers.readCSVXS("etc/8TeV-vbfH.txt")
wHDict8 = helpers.readCSVXS("etc/8TeV-WH.txt")
zHDict8 = helpers.readCSVXS("etc/8TeV-ZH.txt")
ggHDict7 = helpers.readCSVXS("etc/7TeV-ggH.txt")
vbfHDict7 = helpers.readCSVXS("etc/7TeV-vbfH.txt")
wHDict7 = helpers.readCSVXS("etc/7TeV-WH.txt")
zHDict7 = helpers.readCSVXS("etc/7TeV-ZH.txt")
class CrossSections:
def __init__(self):
self.data = {}
self.br = helpers.readCSVXS("etc/brHmm.txt")
def __init__(self):
self.data = {}
self.br = helpers.readCSVXS("etc/br.csv")
self.vbf = {}
self.gg = {}
self.wh = {}
self.zh = {}
for e in ["8TeV","7TeV"]:
self.vbf[e] = helpers.readCSVXS("etc/vbfH_"+e+".csv")
self.gg[e] = helpers.readCSVXS("etc/ggH_"+e+".csv")
self.wh[e] = helpers.readCSVXS("etc/wH_"+e+".csv")
self.zh[e] = helpers.readCSVXS("etc/zH_"+e+".csv")
self.lumi = {
#"14TeV" : 1.044, # 2012 HF-Lumi
#"8TeV" : 1.044, # 2012 HF-Lumi
"14TeV" : 1.026, # 2012 Pixel-Lumi
"8TeV" : 1.026, # 2012 Pixel-Lumi
"7TeV" : 1.022,
}
self.JES = {
'gg' : {
'7TeV' : {
'Jet2CutsFailVBFGF' : 1.0832,
'Jet2CutsGFPass' : 1.0585,
'Jet2CutsVBFPass' : 1.0759,
'Jets01FailPtG10' : None,
'Jets01PassPtG10' : None,
},
'8TeV' : {
'Jet2CutsFailVBFGF' : 1.0799,
'Jet2CutsGFPass' : 1.0477,
'Jet2CutsVBFPass' : 1.0606,
'Jets01FailPtG10' : None,
'Jets01PassPtG10' : None,
},
},
'vbf' : {
'7TeV' : {
'Jet2CutsFailVBFGF' : 1.0572,
'Jet2CutsGFPass' : 1.0381,
'Jet2CutsVBFPass' : 1.0358,
'Jets01FailPtG10' : -1.0321,
'Jets01PassPtG10' : -1.0308,
},
'8TeV' : {
'Jet2CutsFailVBFGF' : 1.0482,
'Jet2CutsGFPass' : 1.0249,
'Jet2CutsVBFPass' : 1.0185,
'Jets01FailPtG10' : -1.0278,
'Jets01PassPtG10' : -1.0250,
},
},
'wh' : {
'7TeV' : {
'Jet2CutsFailVBFGF' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsVBFPass' : None,
'Jets01FailPtG10' : None,
'Jets01PassPtG10' : None,
},
'8TeV' : {
'Jet2CutsFailVBFGF' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsVBFPass' : None,
'Jets01FailPtG10' : None,
'Jets01PassPtG10' : None,
},
},
'zh' : {
'7TeV' : {
'Jet2CutsFailVBFGF' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsVBFPass' : None,
'Jets01FailPtG10' : None,
'Jets01PassPtG10' : None,
},
'8TeV' : {
'Jet2CutsFailVBFGF' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsVBFPass' : None,
'Jets01FailPtG10' : None,
'Jets01PassPtG10' : None,
},
},
}
self.JER = {
'gg' : {
'7TeV' : {
'Jet2CutsFailVBFGF' : -1.0166,
'Jet2CutsGFPass' : -1.0126,
'Jet2CutsVBFPass' : 1.0321,
'Jets01FailPtG10' : None,
'Jets01PassPtG10' : None,
},
'8TeV' : {
'Jet2CutsFailVBFGF' : -1.0131,
'Jet2CutsGFPass' : -1.0101,
'Jet2CutsVBFPass' : -1.0273,
'Jets01FailPtG10' : None,
'Jets01PassPtG10' : None,
},
},
'vbf' : {
'7TeV' : {
'Jet2CutsFailVBFGF' : -1.0111,
'Jet2CutsGFPass' : None,
'Jet2CutsVBFPass' : None,
'Jets01FailPtG10' : -1.0167,
'Jets01PassPtG10' : None,
},
'8TeV' : {
'Jet2CutsFailVBFGF' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsVBFPass' : None,
'Jets01FailPtG10' : -1.0162,
'Jets01PassPtG10' : None,
},
},
'wh' : {
'7TeV' : {
'Jet2CutsFailVBFGF' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsVBFPass' : None,
'Jets01FailPtG10' : None,
'Jets01PassPtG10' : None,
},
'8TeV' : {
'Jet2CutsFailVBFGF' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsVBFPass' : None,
'Jets01FailPtG10' : None,
'Jets01PassPtG10' : None,
},
},
'zh' : {
'7TeV' : {
'Jet2CutsFailVBFGF' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsVBFPass' : None,
'Jets01FailPtG10' : None,
'Jets01PassPtG10' : None,
},
'8TeV' : {
'Jet2CutsFailVBFGF' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsVBFPass' : None,
'Jets01FailPtG10' : None,
'Jets01PassPtG10' : None,
},
},
}
self.PUID = {
'gg' : {
'7TeV' : {
'Jet2CutsFailVBFGF' : 1.0109,
'Jet2CutsGFPass' : 1.0158,
'Jet2CutsVBFPass' : 1.0368,
'Jets01FailPtG10' : None,
'Jets01PassPtG10' : None,
},
'8TeV' : {
'Jet2CutsFailVBFGF' : 1.0142,
'Jet2CutsGFPass' : 1.0174,
'Jet2CutsVBFPass' : 1.0386,
'Jets01FailPtG10' : None,
'Jets01PassPtG10' : None,
},
},
'vbf' : {
'7TeV' : {
'Jet2CutsFailVBFGF' : 1.0134,
'Jet2CutsGFPass' : 1.0160,
'Jet2CutsVBFPass' : 1.0324,
'Jets01FailPtG10' : 1.0115,
'Jets01PassPtG10' : 1.0108,
},
'8TeV' : {
'Jet2CutsFailVBFGF' : 1.0155,
'Jet2CutsGFPass' : 1.0168,
'Jet2CutsVBFPass' : 1.0328,
'Jets01FailPtG10' : 1.0129,
'Jets01PassPtG10' : 1.0128,
},
},
'wh' : {
'7TeV' : {
'Jet2CutsFailVBFGF' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsVBFPass' : None,
'Jets01FailPtG10' : None,
'Jets01PassPtG10' : None,
},
'8TeV' : {
'Jet2CutsFailVBFGF' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsVBFPass' : None,
'Jets01FailPtG10' : None,
'Jets01PassPtG10' : None,
},
},
'zh' : {
'7TeV' : {
'Jet2CutsFailVBFGF' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsVBFPass' : None,
'Jets01FailPtG10' : None,
'Jets01PassPtG10' : None,
},
'8TeV' : {
'Jet2CutsFailVBFGF' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsVBFPass' : None,
'Jets01FailPtG10' : None,
'Jets01PassPtG10' : None,
},
},
}
self.MCStat = {
'gg' : {
'7TeV' : {
'Jet2CutsFailVBFGF' : 1.0301,
'Jet2CutsGFPass' : 1.0572,
'Jet2CutsVBFPass' : 1.1628,
'Jets01FailPtG10' : 1.0154,
'Jets01PassPtG10' : None,
},
'8TeV' : {
'Jet2CutsFailVBFGF' : 1.0308,
'Jet2CutsGFPass' : 1.0577,
'Jet2CutsVBFPass' : 1.1068,
'Jets01FailPtG10' : 1.0147,
'Jets01PassPtG10' : None,
},
},
'vbf' : {
'7TeV' : {
'Jet2CutsFailVBFGF' : 1.0218,
'Jet2CutsGFPass' : 1.0177,
'Jet2CutsVBFPass' : 1.0208,
'Jets01FailPtG10' : 1.0628,
'Jets01PassPtG10' : None,
},
'8TeV' : {
'Jet2CutsFailVBFGF' : 1.0178,
'Jet2CutsGFPass' : 1.0152,
'Jet2CutsVBFPass' : 1.0135,
'Jets01FailPtG10' : 1.0642,
'Jets01PassPtG10' : None,
},
},
'wh' : {
'7TeV' : {
'Jet2CutsFailVBFGF' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsVBFPass' : None,
'Jets01FailPtG10' : None,
'Jets01PassPtG10' : None,
},
'8TeV' : {
'Jet2CutsFailVBFGF' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsVBFPass' : None,
'Jets01FailPtG10' : None,
'Jets01PassPtG10' : None,
},
},
'zh' : {
'7TeV' : {
'Jet2CutsFailVBFGF' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsVBFPass' : None,
'Jets01FailPtG10' : None,
'Jets01PassPtG10' : None,
},
'8TeV' : {
'Jet2CutsFailVBFGF' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsVBFPass' : None,
'Jets01FailPtG10' : None,
'Jets01PassPtG10' : None,
},
},
}
self.PDF = {
'gg' : {
'8TeV' : {
'Jets01PassPtG10' : 1.0799,
'Jets01FailPtG10' : 1.0844,
'Jet2CutsVBFPass' : 1.0578,
'Jet2CutsGFPass' : 1.0715,
'Jet2CutsFailVBFGF' : 1.0809,
},
'7TeV' : {
'Jets01PassPtG10' : 1.0830,
'Jets01FailPtG10' : 1.0868,
'Jet2CutsVBFPass' : 1.0593,
'Jet2CutsGFPass' : 1.0743,
'Jet2CutsFailVBFGF' : 1.0801,
},
},
'vbf' : {
'8TeV' : {
'Jets01PassPtG10' : 1.0197,
'Jets01FailPtG10' : 1.0186,
'Jet2CutsVBFPass' : 1.0378,
'Jet2CutsGFPass' : 1.0212,
'Jet2CutsFailVBFGF' : 1.0196,
},
'7TeV' : {
'Jets01PassPtG10' : 1.0220,
'Jets01FailPtG10' : 1.0214,
'Jet2CutsVBFPass' : 1.0448,
'Jet2CutsGFPass' : 1.0299,
'Jet2CutsFailVBFGF' : 1.0204,
},
},
'wh' : {
'8TeV' : {
'Jets01PassPtG10' : None,
'Jets01FailPtG10' : None,
'Jet2CutsVBFPass' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsFailVBFGF' : None,
},
'7TeV' : {
'Jets01PassPtG10' : None,
'Jets01FailPtG10' : None,
'Jet2CutsVBFPass' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsFailVBFGF' : None,
},
},
'zh' : {
'8TeV' : {
'Jets01PassPtG10' : None,
'Jets01FailPtG10' : None,
'Jet2CutsVBFPass' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsFailVBFGF' : None,
},
'7TeV' : {
'Jets01PassPtG10' : None,
'Jets01FailPtG10' : None,
'Jet2CutsVBFPass' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsFailVBFGF' : None,
},
},
}
self.PU = {
'gg' : {
'8TeV' : {
'Jets01PassPtG10' : None,
'Jets01FailPtG10' : None,
'Jet2CutsVBFPass' : 1.0222,
'Jet2CutsGFPass' : 1.0103,
'Jet2CutsFailVBFGF' : 1.0112,
},
'7TeV' : {
'Jets01PassPtG10' : None,
'Jets01FailPtG10' : None,
'Jet2CutsVBFPass' : 1.0222,
'Jet2CutsGFPass' : 1.0103,
'Jet2CutsFailVBFGF' : 1.0112,
},
},
'vbf' : {
'8TeV' : {
'Jets01PassPtG10' : None,
'Jets01FailPtG10' : 1.0207,
'Jet2CutsVBFPass' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsFailVBFGF' : None,
},
'7TeV' : {
'Jets01PassPtG10' : None,
'Jets01FailPtG10' : 1.0207,
'Jet2CutsVBFPass' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsFailVBFGF' : None,
},
},
'wh' : {
'8TeV' : {
'Jets01PassPtG10' : None,
'Jets01FailPtG10' : None,
'Jet2CutsVBFPass' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsFailVBFGF' : None,
},
'7TeV' : {
'Jets01PassPtG10' : None,
'Jets01FailPtG10' : None,
'Jet2CutsVBFPass' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsFailVBFGF' : None,
},
},
'zh' : {
'8TeV' : {
'Jets01PassPtG10' : None,
'Jets01FailPtG10' : None,
'Jet2CutsVBFPass' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsFailVBFGF' : None,
},
'7TeV' : {
'Jets01PassPtG10' : None,
'Jets01FailPtG10' : None,
'Jet2CutsVBFPass' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsFailVBFGF' : None,
},
},
}
self.QCDScale = {
'gg' : {
'8TeV' : {
'Jets01PassPtG10' : -1.0121,
'Jets01FailPtG10' : 1.0450,
'Jet2CutsVBFPass' : -1.1585,
'Jet2CutsGFPass' : 1.1515,
'Jet2CutsFailVBFGF' : -1.0662,
},
'7TeV' : {
'Jets01PassPtG10' : 1.0104,
'Jets01FailPtG10' : 1.0417,
'Jet2CutsVBFPass' : -1.2424,
'Jet2CutsGFPass' : -1.1608,
'Jet2CutsFailVBFGF' : 1.0626,
},
},
'vbf' : {
'8TeV' : {
'Jets01PassPtG10' : -1.0139,
'Jets01FailPtG10' : -1.1598,
'Jet2CutsVBFPass' : 1.0313,
'Jet2CutsGFPass' : -1.0372,
'Jet2CutsFailVBFGF' : -1.0427,
},
'7TeV' : {
'Jets01PassPtG10' : 1.0101,
'Jets01FailPtG10' : 1.0571,
'Jet2CutsVBFPass' : -1.0489,
'Jet2CutsGFPass' : -1.0380,
'Jet2CutsFailVBFGF' : 1.0653,
},
},
'wh' : {
'8TeV' : {
'Jets01PassPtG10' : None,
'Jets01FailPtG10' : None,
'Jet2CutsVBFPass' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsFailVBFGF' : None,
},
'7TeV' : {
'Jets01PassPtG10' : None,
'Jets01FailPtG10' : None,
'Jet2CutsVBFPass' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsFailVBFGF' : None,
},
},
'zh' : {
'8TeV' : {
'Jets01PassPtG10' : None,
'Jets01FailPtG10' : None,
'Jet2CutsVBFPass' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsFailVBFGF' : None,
},
'7TeV' : {
'Jets01PassPtG10' : None,
'Jets01FailPtG10' : None,
'Jet2CutsVBFPass' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsFailVBFGF' : None,
},
},
}
self.UE = {
'gg' : {
'8TeV' : {
'Jets01PassPtG10' : -1.0618,
'Jets01FailPtG10' : 1.1698,
'Jet2CutsVBFPass' : 1.4228,
'Jet2CutsGFPass' : 1.2317,
'Jet2CutsFailVBFGF' : -1.0847,
},
'7TeV' : {
'Jets01PassPtG10' : -1.0658,
'Jets01FailPtG10' : 1.1715,
'Jet2CutsVBFPass' : 1.6606,
'Jet2CutsGFPass' : 1.2570,
'Jet2CutsFailVBFGF' : -1.0669,
},
},
'vbf' : {
'8TeV' : {
'Jets01PassPtG10' : -1.0309,
'Jets01FailPtG10' : -1.0794,
'Jet2CutsVBFPass' : -1.1030,
'Jet2CutsGFPass' : -1.0585,
'Jet2CutsFailVBFGF' : 1.0432,
},
'7TeV' : {
'Jets01PassPtG10' : -1.0515,
'Jets01FailPtG10' : 1.0953,
'Jet2CutsVBFPass' : -1.1421,
'Jet2CutsGFPass' : -1.1014,
'Jet2CutsFailVBFGF' : 1.0763,
},
},
'wh' : {
'8TeV' : {
'Jets01PassPtG10' : None,
'Jets01FailPtG10' : None,
'Jet2CutsVBFPass' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsFailVBFGF' : None,
},
'7TeV' : {
'Jets01PassPtG10' : None,
'Jets01FailPtG10' : None,
'Jet2CutsVBFPass' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsFailVBFGF' : None,
},
},
'zh' : {
'8TeV' : {
'Jets01PassPtG10' : None,
'Jets01FailPtG10' : None,
'Jet2CutsVBFPass' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsFailVBFGF' : None,
},
'7TeV' : {
'Jets01PassPtG10' : None,
'Jets01FailPtG10' : None,
'Jet2CutsVBFPass' : None,
'Jet2CutsGFPass' : None,
'Jet2CutsFailVBFGF' : None,
},
},
}
# The list of systematics to be applied
self._keys = ["xs_ggH","xs_vbfH","xs_whH","xs_zhH","br_Hmm","lumi","PDF","JES","JER","PUID","MCStat","PU","UE","QCDScale"]
# The list of systematics which are correlated between energies
self.keysEnergyCorr = ["xs_ggH","xs_vbfH","xs_whH","xs_zhH","br_Hmm","PDF","JER","UE","QCDScale"]
# The list of systematics which are not correlated between energies or categories
self.keysNotCatCorr = ["MCStat"]
# The list of systematics which are not correlated between energies, but correlated w/ categories
# (everything else)
self.keysNotEnergyCorr = []
for i in self._keys:
if not i in self.keysEnergyCorr and not i in self.keysNotCatCorr:
self.keysNotEnergyCorr.append(i)
lumiDict={}
#lumiDict["8TeV"] = 19.79 #2012ABCD 22Jan2013 HF-Lumi
lumiDict["8TeV"] = 19.712 #2012ABCD 22Jan2013 Pixel-Lumi
lumiDict["7TeV"] = 5.05 #2011AB
#LUMI=0.807 #2012A
#LUMI=4.421 #2012B
#LUMI=0.495 #2012Cv1
#LUMI=6.311 #2012Cv2
#LUMI=2.311 #2011A
#LUMI=2.739 #2011B
brDict = helpers.readCSVXS("etc/br.csv")
ggHDict8 = helpers.readCSVXS("etc/ggH_8TeV.csv")
vbfHDict8 = helpers.readCSVXS("etc/vbfH_8TeV.csv")
wHDict8 = helpers.readCSVXS("etc/wH_8TeV.csv")
zHDict8 = helpers.readCSVXS("etc/zH_8TeV.csv")
ggHDict7 = helpers.readCSVXS("etc/ggH_7TeV.csv")
vbfHDict7 = helpers.readCSVXS("etc/vbfH_7TeV.csv")
wHDict7 = helpers.readCSVXS("etc/wH_7TeV.csv")
zHDict7 = helpers.readCSVXS("etc/zH_7TeV.csv")
class CrossSections:
def __init__(self):
self.data = {}
self.br = helpers.readCSVXS("etc/br.csv")
self.vbf = {}