Exemplo n.º 1
0
 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")
Exemplo n.º 2
0
 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")
Exemplo n.º 3
0
 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")
Exemplo n.º 4
0
    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)
Exemplo n.º 5
0
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")
Exemplo n.º 6
0
    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{})",
        }
Exemplo n.º 7
0
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")
Exemplo n.º 8
0
  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)
Exemplo n.º 9
0
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 = {}