Merge branch 'master' of github.com:arizzi/PisaHmm into masteronline

arizzi · Mar 10, 2020 · 7ea6faa · 7ea6faa
2 parents ccd9e63 + d47bdd2
commit 7ea6faa
Show file tree

Hide file tree

Showing 7 changed files with 165 additions and 122 deletions.
diff --git a/models2016H.py b/models2016H.py
@@ -1,37 +1,38 @@
 from samples2016 import *
 name="H"
 background={
-"DY":["DY105_2016AMCPY"], ## Alternative: "DY105_2016MGPY"
-"DYVBF":["DY105VBF_2016AMCPY"], ## Alternative: "DY105VBF_2016MGPY"
-"EWKZ":[
+#"DY":["DY105_2016AMCPY"], ## Alternative: "DY105_2016MGPY"
+#"DYVBF":["DY105VBF_2016AMCPY"], ## Alternative: "DY105VBF_2016MGPY"
+"DY":["DY105_2016AMCPY", "DY105VBF_2016AMCPY"], 
+"VBF Z":[
    #"EWKZ105FIX2_2016MGHERWIG",
    "EWKZ105FIX2_2016MGHERWIG",
    # interference with DY
-#   "EWKZint_2016MGPY"
+  "EWKZint_2016MGPY"
 ],
 "Top":[
    ### Single Top (s, t, tW channels) ###
- #  "STs_2016AMCPY",
- #  "STt_2016POW_MADSPIN_PY",
+  "STs_2016AMCPY",
+  "STt_2016POW_MADSPIN_PY",
    "STwt_2016POWPY", 
 
    ### Single Anti-Top (s, t, tW channels) ###
    #"STsbar_2016AMCPY",
-#   "STtbar_2016POW_MADSPIN_PY",
+  "STtbar_2016POW_MADSPIN_PY",
    "STwtbar_2016POWPY",
 
    ### TTbar (leptonic, semileptonic, hadronic)
    "TTlep_2016POWPY",   # 2 lept
-### MISSING ### "TTsemi_2016POWPY", # 1 lept
+   "TTsemi_2016POWPY", # 1 lept
 ### MISSING ### "TThad_2016POWPY",  # 0 lept
    ### TTbar alternatives (inclusive): "TT_2016POWPY", "TT_2016AMCPY",
 ],
 "Other":[
          ### W+jets ###
          #--- inclusive ---
-  #       "W2J_2016AMCPY",
-   #      "W1J_2016AMCPY",
-    #     "W0J_2016AMCPY", 
+        "W2J_2016AMCPY",
+        "W1J_2016AMCPY",
+        "W0J_2016AMCPY", 
 
          ### WW ###
          #--- 2 lept ---
@@ -42,21 +43,21 @@
 
          ### WZ ###
          #--- 3 lept ---
-#         "WZ3l1n_2016AMCPY", ## Alternative: #"WZ3l1n_2016POWPY"
+        "WZ3l1n_2016AMCPY", ## Alternative: #"WZ3l1n_2016POWPY"
          #--- 2 lept ---
          "WZ2l2q_2016AMC_MADSPIN_PY",
          #--- 1 lept ---
 ### MISSING ### "WZ1l1n2q_2016AMCPY",
- #        "WZ1l3n_2016AMCPY", 
+        "WZ1l3n_2016AMCPY", 
 
          ### ZZ ###
          #--- 4 lept ---
 ### MISSING ### "ZZ4l_2016POWPY",
          #--- 2 lept ---
-  #       "ZZ2l2q_2016POWPY",
+        "ZZ2l2q_2016POWPY",
 ### MISSING ### "ZZ2l2n_2016POWPY",
          #--- 0 lept ---
-### MISSING ### "ZZ2q2n_2016POWPY",
+## MISSING ### "ZZ2q2n_2016POWPY",
 
          ### Vector boson scattering ###
 ### LHE Weights broken "WWJJlnln_2016MGPY",          ## VBS W(lv)W(ln) + 2jets 
@@ -69,18 +70,26 @@
 
 
 #sorting
-backgroundSorted=["Other","Top","DY","DYVBF","EWKZ"]
+backgroundSorted=["Other","Top","DY","VBF Z"]
 backgroundSorted+=[x for x in background if x not in backgroundSorted]
 
 
 signal={
 "VBF H":["vbfHmm_2016AMCPY"], ## Alternative: "vbfHmm_2016POWPY"
 "gg H":["ggHmm_2016AMCPY"],   ## Alternative: "ggHmm_2016POWPY"
-#ZH":["zHmm_2016POWPY"],
-#WH":["WplusHmm_2016POWPY","WminusHmm_2016POWPY"],
-#ttH":["ttHmm_2016POWPY"]
+#"ZH":["zHmm_2016POWPY"],
+#"WH":["WplusHmm_2016POWPY","WminusHmm_2016POWPY"],
+"VH":["WplusHmm_2016POWPY","WminusHmm_2016POWPY", "zHmm_2016POWPY"],
+"ttH":["ttHmm_2016POWPY"]
 }
 
+#legend sorting
+backgroundSortedForLegend=["DY","VBF Z","Top", "Other"]
+backgroundSortedForLegend+=[x for x in background if x not in backgroundSortedForLegend]
+signalSortedForLegend=["VBF H","gg H"]
+signalSortedForLegend+=[x for x in signal if x not in signalSortedForLegend]
+
+
 data={
 "2016":["data2016"]
 }
@@ -89,16 +98,19 @@
 fillcolor={
 "DY": ROOT.kOrange,
 "DYVBF": ROOT.kOrange-3,
-"EWKZ": ROOT.kViolet,
-"Top": ROOT.kGreen,
-"Other" : ROOT.kGreen+1,
+"VBF Z": ROOT.kMagenta+2,
+"Top": ROOT.kGreen+1,
+"Other" : ROOT.kGreen+3,
 "VBF H":ROOT.kRed,
 "gg H":ROOT.kRed+4,
 "ZH":ROOT.kPink+4,
 "WH":ROOT.kPink+9,
+"VH":ROOT.kPink+5,
 "ttH":ROOT.kRed-4,
 }
 
+
+
 #systematicsToPlot=["MuScaleUp"]
 #systematicsToPlot=["JERUp","JERDown","JESUp","JESDown","WithJER","puWeightUp","puWeightDown"]
 #ystematicsToPlot=["JERUp","JERDown","puWeightUp","puWeightDown","LHERenUp","LHERenDown","LHEFacUp","LHEFacDown","MuScaleUp","MuScaleDown"]

diff --git a/models2016Z.py b/models2016Z.py
@@ -2,10 +2,11 @@
 from samples2016 import *
 name="Z"
 background={
-"DY0J":["DY0J_2016AMCPY"],
-"DY1J":["DY1J_2016AMCPY"],
-"DY2J":["DY2J_2016AMCPY"],
-"EWKZ":[
+#"DY0J":["DY0J_2016AMCPY"],
+#"DY1J":["DY1J_2016AMCPY"],
+#"DY2J":["DY2J_2016AMCPY"],
+"DY":["DY0J_2016AMCPY", "DY1J_2016AMCPY", "DY2J_2016AMCPY"],
+"VBF Z":[
    "EWKZ_2016MGHERWIG",  ## Alternative: "EWKZ_2016MGPY"
    "EWKZint_2016MGPY", # interference with DY
 ],
@@ -17,23 +18,23 @@
 data = models2016H.data
 
 #sorting
-backgroundSorted=["Other","Top","DY0J","DY1J","DY2J","EWKZ"]
+backgroundSorted=["Other","Top","DY","VBF Z"]
 backgroundSorted+=[x for x in background if x not in backgroundSorted]
 
+#legend sorting
+backgroundSortedForLegend=["DY","VBF Z","Top", "Other"]
+backgroundSortedForLegend+=[x for x in background if x not in backgroundSortedForLegend]
+signalSortedForLegend=["VBF H","gg H"]
+signalSortedForLegend+=[x for x in signal if x not in signalSortedForLegend]
+
+
 import ROOT
 fillcolor={
 "DY0J": ROOT.kOrange+2,
 "DY1J": ROOT.kOrange+1,
 "DY2J": ROOT.kOrange,
-"EWKZ": ROOT.kViolet,
-"Top": ROOT.kGreen,
-"Other" : ROOT.kGreen+1,
-"VBF H":ROOT.kRed,
-"gg H":ROOT.kRed+4,
-"ZH":ROOT.kPink+4,
-"WH":ROOT.kPink+9,
-"ttH":ROOT.kRed-4,
 }
+fillcolor.update(models2016H.fillcolor)
 
 #systematicsToPlot=["JERUp","JERDown","JESUp","JESDown","puWeightUp","puWeightDown"]
 #systematicsToPlot=["JERUp","JERDown","puWeightUp","puWeightDown","LHERenUp","LHERenDown","LHEFacUp","LHEFacDown","MuScaleUp","MuScaleDown"]

diff --git a/models2017H.py b/models2017H.py
@@ -1,59 +1,61 @@
 from samples2017 import *
 name="H"
 background={
-"DY":["DY105_2017AMCPY"], ## Alternative: "DY105_2017MGPY"
-"DYVBF":[ "DY105VBF_2017AMCPY"], ## Alternative: "DY105VBF_2017MGPY"
-"EWKZ":[
+#"DY":["DY105_2017AMCPY"], ## Alternative: "DY105_2017MGPY"
+#"DYVBF":[ "DY105VBF_2017AMCPY"], ## Alternative: "DY105VBF_2017MGPY"
+"DY":["DY105_2017AMCPY", "DY105VBF_2017AMCPY"], ## Alternative: "DY105VBF_2017MGPY"
+"VBF Z":[
    "EWKZ105FIX2_2017MGHERWIG",
    # interference with DY
-#   "EWKZint_2017MGPY" 
+  "EWKZint_2017MGPY" 
 ],
 "Top":[
    ### Single Top (s, t, tW channels) ###
- #  "STs_2017AMCPY",
-### MISSING ### "STt_2017POW_MADSPIN_PY",
+   "STs_2017AMCPY",
+   "STt_2017POWPY",
    "STwt_2017POWPY", 
 
    ### Single Anti-Top (s, t, tW channels) ###
 ### MISSING ### "STsbar_2017AMCPY",
-### MISSING ### "STtbar_2017POW_MADSPIN_PY",
+   #"STtbar_2017POW_MADSPIN_PY",
    "STwtbar_2017POWPY",
 
    ### TTbar (leptonic, semileptonic, hadronic)
    "TTlep_2017POWPY",   # 2 lept
-  # "TTsemi_2017POWPY",  # 1 lept
+  "TTsemi_2017POWPY",  # 1 lept
    #"TThad_2017POWPY",  # 0 lept
    ### TTbar alternatives (inclusive): "TT_2017POWPY", "TT_2017AMCPY",
 ],
 "Other":[
          ### W+jets ###
          #--- inclusive ---
-#         "W2J_2017AMCPY",
-### MISSING ### "W1J_2017AMCPY",
-   #      "W0J_2017AMCPY", 
+        "W2J_2017AMCPY",
+       ##"W1J_2017AMCPY",
+        "W0J_2017AMCPY", 
 
          ### WW ###
          #--- 2 lept ---
          "WW2l2n_2017POWPY",
+         "WWlnqq_2017POWPY",
          #--- 1 lept ---
 ### MISSING ### "WWlnqq_2017AMC_MADSPIN_PY", ## Alternative: "WWlnqq_2017POWPY",
          #--- incl ---
 
          ### WZ ###
          #--- 3 lept ---
-    #     "WZ3l1n_2017AMCPY", ## Alternative #"WZ3l1n_2017POWPY"
+        "WZ3l1n_2017AMCPY", ## Alternative #"WZ3l1n_2017POWPY"
          #--- 2 lept ---
-         "WZ2l2q_2017AMC_MADSPIN_PY",
+         ### MISSING ###"WZ2l2q_2017AMC_MADSPIN_PY",
          #--- 1 lept ---
-     #    "WZ1l1n2q_2017AMCPY",
-      #   "WZ1l3n_2017AMCPY", 
+        "WZ1l1n2q_2017AMCPY",
+        "WZ1l3n_2017AMCPY", 
 
          ### ZZ ###
          #--- 4 lept ---
-       #  "ZZ4l_2017POWPY",
+       ##  "ZZ4l_2017POWPY",
          #--- 2 lept ---
-        # "ZZ2l2q_2017POWPY",
-         #"ZZ2l2n_2017POWPY",
+        "ZZ2l2q_2017POWPY",
+         ### MISSING ###"ZZ2l2n_2017POWPY",
          #--- 0 lept ---
 ### MISSING ### "ZZ2q2n_2017POWPY",
 
@@ -68,18 +70,23 @@
 
 
 #sorting
-backgroundSorted=["Other","Top","DY","DYVBF","EWKZ"]
+backgroundSorted=["Other","Top","DY","VBF Z"]
 backgroundSorted+=[x for x in background if x not in backgroundSorted]
 
-
 signal={
 "VBF H":["vbfHmm_2017AMCPY"], ## Alternative: "vbfHmm_2017POWPY"
 "gg H":["ggHmm_2017AMCPY"],   ## Alternative: "ggHmm_2017POWPY"
-#"ZH":["zHmm_2017POWPY"],
-#"WH":["WplusHmm_2017POWPY","WminusHmm_2017POWPY"],
-#"ttH":["ttHmm_2017POWPY"]
+"VH":["WplusHmm_2017POWPY","WminusHmm_2017POWPY", "zHmm_2017POWPY"],
+"ttH":["ttHmm_2017POWPY"]
 }
 
+#legend sorting
+backgroundSortedForLegend=["DY","VBF Z","Top", "Other"]
+backgroundSortedForLegend+=[x for x in background if x not in backgroundSortedForLegend]
+signalSortedForLegend=["VBF H","gg H"]
+signalSortedForLegend+=[x for x in signal if x not in signalSortedForLegend]
+
+
 data={
 "2017":["data2017"]
 }
@@ -88,13 +95,14 @@
 fillcolor={
 "DY": ROOT.kOrange,
 "DYVBF": ROOT.kOrange-3,
-"EWKZ": ROOT.kViolet,
-"Top": ROOT.kGreen,
-"Other" : ROOT.kGreen+1,
+"VBF Z": ROOT.kMagenta+2,
+"Top": ROOT.kGreen+1,
+"Other" : ROOT.kGreen+3,
 "VBF H":ROOT.kRed,
 "gg H":ROOT.kRed+4,
 "ZH":ROOT.kPink+4,
 "WH":ROOT.kPink+9,
+"VH":ROOT.kPink+5,
 "ttH":ROOT.kRed-4,
 }
 

diff --git a/models2017Z.py b/models2017Z.py
@@ -2,10 +2,11 @@
 from samples2017 import *
 name="Z"
 background={
-"DY0J":["DY0J_2017AMCPY"],
-"DY1J":["DY1J_2017AMCPY"],
-"DY2J":["DY2J_2017AMCPY"],
-"EWKZ":[
+#"DY0J":["DY0J_2017AMCPY"],
+#"DY1J":["DY1J_2017AMCPY"],
+#"DY2J":["DY2J_2017AMCPY"],
+"DY":["DY0J_2017AMCPY", "DY1J_2017AMCPY", "DY2J_2017AMCPY"],
+"VBF Z":[
    "EWKZ_2017MGHERWIG",  ## Alternative: "EWKZ_2017MGPY"
    "EWKZint_2017MGPY", # interference with DY
 ],
@@ -17,23 +18,23 @@
 data = models2017H.data
 
 #sorting
-backgroundSorted=["Other","Top","DY0J","DY1J","DY2J","EWKZ"]
+backgroundSorted=["Other","Top","DY","VBF Z"]
 backgroundSorted+=[x for x in background if x not in backgroundSorted]
 
+#legend sorting
+backgroundSortedForLegend=["DY","VBF Z","Top", "Other"]
+backgroundSortedForLegend+=[x for x in background if x not in backgroundSortedForLegend]
+signalSortedForLegend=["VBF H","gg H"]
+signalSortedForLegend+=[x for x in signal if x not in signalSortedForLegend]
+
+
 import ROOT
 fillcolor={
 "DY0J": ROOT.kOrange+2,
 "DY1J": ROOT.kOrange+1,
 "DY2J": ROOT.kOrange,
-"EWKZ": ROOT.kViolet,
-"Top": ROOT.kGreen,
-"Other" : ROOT.kGreen+1,
-"VBF H":ROOT.kRed,
-"gg H":ROOT.kRed+4,
-"ZH":ROOT.kPink+4,
-"WH":ROOT.kPink+9,
-"ttH":ROOT.kRed-4,
 }
+fillcolor.update(models2017H.fillcolor)
 
 #systematicsToPlot=["JERUp","JERDown","JESUp","JESDown","puWeightUp","puWeightDown"]
 #systematicsToPlot=["JERUp","JERDown","puWeightUp","puWeightDown","LHERenUp","LHERenDown","LHEFacUp","LHEFacDown","MuScaleUp","MuScaleDown"]