-
Notifications
You must be signed in to change notification settings - Fork 4
/
models2017Z.py
62 lines (47 loc) · 1.96 KB
/
models2017Z.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
import models2017H
from samples2017 import *
name="Z"
background={
#"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
],
}
background["Top"] = models2017H.background["Top"]
background["Other"] = models2017H.background["Other"]
signal = models2017H.signal
data = models2017H.data
#sorting
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,
}
fillcolor.update(models2017H.fillcolor)
#systematicsToPlot=["JERUp","JERDown","JESUp","JESDown","puWeightUp","puWeightDown"]
#systematicsToPlot=["JERUp","JERDown","puWeightUp","puWeightDown","LHERenUp","LHERenDown","LHEFacUp","LHEFacDown","MuScaleUp","MuScaleDown"]
systematicsToPlot=["PrefiringWeightUp","PrefiringWeightDown","LHEPdfUp","LHEPdfDown","QGLweightUp","QGLweightDown","JERUp","JERDown","puWeightUp","puWeightDown","LHERenUp","LHERenDown","LHEFacUp","LHEFacDown","MuScaleUp","MuScaleDown","AlternativeUp","AlternativeDown","PDFX1Up","PDFX1Down","PDFX0Up","PDFX0Down"]
systematicsToPlot+=["JESUp","JESDown"]
from btagvariations import btagsys
systematicsToPlot+=btagsys
from jesnames import jesnames2017
from jernames import jernames
jesList=jesnames2017
systematicsForDC=systematicsToPlot+[x[7:] for x in jesList ]+jernames
linecolor=fillcolor
markercolor=fillcolor
from rebinning import *
from systematicGrouping import *
systematicDetail = systematicGrouping(background, signal,jesList,"2017")