-
Notifications
You must be signed in to change notification settings - Fork 0
/
caloParams_2022_v0_1_cfi.py
200 lines (147 loc) · 13.4 KB
/
caloParams_2022_v0_1_cfi.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
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
import FWCore.ParameterSet.Config as cms
from L1Trigger.L1TCalorimeter.caloParams_cfi import caloParamsSource
import L1Trigger.L1TCalorimeter.caloParams_cfi
caloStage2Params = L1Trigger.L1TCalorimeter.caloParams_cfi.caloParams.clone()
# towers
caloStage2Params.towerLsbH = cms.double(0.5)
caloStage2Params.towerLsbE = cms.double(0.5)
caloStage2Params.towerLsbSum = cms.double(0.5)
caloStage2Params.towerNBitsH = cms.int32(8)
caloStage2Params.towerNBitsE = cms.int32(8)
caloStage2Params.towerNBitsSum = cms.int32(9)
caloStage2Params.towerNBitsRatio = cms.int32(3)
caloStage2Params.towerEncoding = cms.bool(True)
# regions
caloStage2Params.regionLsb = cms.double(0.5)
caloStage2Params.regionPUSType = cms.string("None")
caloStage2Params.regionPUSParams = cms.vdouble()
# EG
caloStage2Params.egLsb = cms.double(0.5)
caloStage2Params.egEtaCut = cms.int32(28)
caloStage2Params.egSeedThreshold = cms.double(2.)
caloStage2Params.egNeighbourThreshold = cms.double(1.)
caloStage2Params.egHcalThreshold = cms.double(0.)
caloStage2Params.egTrimmingLUTFile = cms.FileInPath("L1Trigger/L1TCalorimeter/data/egTrimmingLUT_10_v16.01.19.txt")
caloStage2Params.egMaxHcalEt = cms.double(0.)
caloStage2Params.egMaxPtHOverE = cms.double(128.)
caloStage2Params.egMaxHOverELUTFile = cms.FileInPath("L1Trigger/L1TCalorimeter/data/HoverEIdentification_0.995_v15.12.23.txt")
caloStage2Params.egHOverEcutBarrel = cms.int32(3)
caloStage2Params.egHOverEcutEndcap = cms.int32(4)
caloStage2Params.egBypassExtHOverE = cms.uint32(0)
caloStage2Params.egCompressShapesLUTFile = cms.FileInPath("L1Trigger/L1TCalorimeter/data/egCompressLUT_v4.txt")
caloStage2Params.egShapeIdType = cms.string("compressed")
caloStage2Params.egShapeIdVersion = cms.uint32(0)
caloStage2Params.egShapeIdLUTFile = cms.FileInPath("L1Trigger/L1TCalorimeter/data/shapeIdentification_adapt0.99_compressedieta_compressedE_compressedshape_v15.12.08.txt")#Not used any more in the current emulator version, merged with calibration LUT
caloStage2Params.egPUSType = cms.string("None")
caloStage2Params.egIsolationType = cms.string("compressed")
caloStage2Params.egIsoLUTFile = cms.FileInPath("L1Trigger/L1TCalorimeter/data/EG_Iso_LUT_04_04_2017.2.txt")
caloStage2Params.egIsoLUTFile2 = cms.FileInPath("L1Trigger/L1TCalorimeter/data/EG_LoosestIso_2018.2.txt")
caloStage2Params.egIsoAreaNrTowersEta = cms.uint32(2)
caloStage2Params.egIsoAreaNrTowersPhi = cms.uint32(4)
caloStage2Params.egIsoVetoNrTowersPhi = cms.uint32(2)
#caloStage2Params.egIsoPUEstTowerGranularity = cms.uint32(1)
#caloStage2Params.egIsoMaxEtaAbsForTowerSum = cms.uint32(4)
#caloStage2Params.egIsoMaxEtaAbsForIsoSum = cms.uint32(27)
caloStage2Params.egPUSParams = cms.vdouble(1,4,32) #Isolation window in firmware goes up to abs(ieta)=32 for now
caloStage2Params.egCalibrationType = cms.string("compressed")
caloStage2Params.egCalibrationVersion = cms.uint32(0)
#caloStage2Params.egCalibrationLUTFile = cms.FileInPath("L1Trigger/L1TCalorimeter/data/EG_Calibration_LUT_FW_v17.04.04_shapeIdentification_adapt0.99_compressedieta_compressedE_compressedshape_v15.12.08_correct.txt")
caloStage2Params.egCalibrationLUTFile = cms.FileInPath("L1Trigger/L1TCalorimeter/data/corrections_Trimming10_compressedieta_compressedE_compressedshape_PANTELIS_v2_NEW_CALIBRATIONS_withShape_v17.04.04.txt")
# Tau
caloStage2Params.tauLsb = cms.double(0.5)
caloStage2Params.isoTauEtaMax = cms.int32(25)
caloStage2Params.tauSeedThreshold = cms.double(0.)
caloStage2Params.tauNeighbourThreshold = cms.double(0.)
caloStage2Params.tauIsoAreaNrTowersEta = cms.uint32(2)
caloStage2Params.tauIsoAreaNrTowersPhi = cms.uint32(4)
caloStage2Params.tauIsoVetoNrTowersPhi = cms.uint32(2)
caloStage2Params.tauPUSType = cms.string("None")
caloStage2Params.tauIsoLUTFile = cms.FileInPath("L1Trigger/L1TCalorimeter/data/Tau_Iso_LUT_Option_31_extrap_2018_FW_v10.0.0.txt")
caloStage2Params.tauIsoLUTFile2 = cms.FileInPath("L1Trigger/L1TCalorimeter/data/Tau_Iso_LUT_Option_31_extrap_2018_FW_v10.0.0.txt")
#caloStage2Params.tauCalibrationLUTFile = cms.FileInPath("L1Trigger/L1TCalorimeter/data/Tau_Calibration_LUT_2017_Layer1Calibration_FW_v12.0.0.txt")
caloStage2Params.tauCalibrationLUTFile = cms.FileInPath("L1Trigger/L1TCalorimeter/data/Tau_Calibration_LUT_2018_Layer1CalibrationNewHCAL_FW_v13.0.0.txt")
caloStage2Params.tauCompressLUTFile = cms.FileInPath("L1Trigger/L1TCalorimeter/data/tauCompressAllLUT_12bit_v3.txt")
caloStage2Params.tauPUSParams = cms.vdouble(1,4,32)
# jets
caloStage2Params.jetLsb = cms.double(0.5)
caloStage2Params.jetSeedThreshold = cms.double(4.0)
caloStage2Params.jetNeighbourThreshold = cms.double(0.)
caloStage2Params.jetPUSType = cms.string("ChunkyDonut")
caloStage2Params.jetBypassPUS = cms.uint32(0)
# Calibration options
caloStage2Params.jetCalibrationType = cms.string("LUT")
caloStage2Params.jetCompressPtLUTFile = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_pt_compress_2017v1.txt")
caloStage2Params.jetCompressEtaLUTFile = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_eta_compress_2017v1.txt")
caloStage2Params.jetCalibrationLUTFile = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_calib_2018v1_ECALZS_noHFJEC.txt")
# sums: 0=ET, 1=HT, 2=MET, 3=MHT
caloStage2Params.etSumLsb = cms.double(0.5)
caloStage2Params.etSumEtaMin = cms.vint32(1, 1, 1, 1, 1)
caloStage2Params.etSumEtaMax = cms.vint32(28, 26, 28, 26, 28)
caloStage2Params.etSumEtThreshold = cms.vdouble(0., 30., 0., 30., 0.) # only 2nd (HT) and 4th (MHT) values applied
caloStage2Params.etSumMetPUSType = cms.string("LUT") # et threshold from this LUT supercedes et threshold in line above
caloStage2Params.etSumEttPUSType = cms.string("None")
caloStage2Params.etSumEcalSumPUSType = cms.string("None")
caloStage2Params.etSumBypassMetPUS = cms.uint32(0)
caloStage2Params.etSumBypassEttPUS = cms.uint32(1)
caloStage2Params.etSumBypassEcalSumPUS = cms.uint32(1)
caloStage2Params.etSumXCalibrationType = cms.string("None")
caloStage2Params.etSumYCalibrationType = cms.string("None")
caloStage2Params.etSumEttCalibrationType = cms.string("None")
caloStage2Params.etSumEcalSumCalibrationType = cms.string("None")
caloStage2Params.etSumMetPUSLUTFile = cms.FileInPath("L1Trigger/L1TCalorimeter/data/newSFHBHEOnp5_METPUM_211124.txt")
caloStage2Params.etSumEttPUSLUTFile = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_towEtThresh_dummy.txt")
caloStage2Params.etSumEcalSumPUSLUTFile = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_towEtThresh_dummy.txt")
caloStage2Params.etSumXCalibrationLUTFile = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_etSumPUS_dummy.txt")
caloStage2Params.etSumYCalibrationLUTFile = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_etSumPUS_dummy.txt")
caloStage2Params.etSumEttCalibrationLUTFile = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_etSumPUS_dummy.txt")
caloStage2Params.etSumEcalSumCalibrationLUTFile = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_etSumPUS_dummy.txt")
# Layer 1 SF
caloStage2Params.layer1ECalScaleETBins = cms.vint32([3, 6, 9, 12, 15, 20, 25, 30, 35, 40, 45, 55, 70, 256])
caloStage2Params.layer1ECalScaleFactors = cms.vdouble([
1.12, 1.13, 1.13, 1.12, 1.12, 1.12, 1.13, 1.13, 1.13, 1.13, 1.13, 1.13, 1.13, 1.14, 1.14, 1.13, 1.13, 1.30, 1.16, 1.28, 1.27, 1.30, 1.30, 1.33, 1.36, 0.00, 0.00, 0.00,
1.12, 1.13, 1.13, 1.12, 1.12, 1.12, 1.13, 1.13, 1.13, 1.13, 1.13, 1.13, 1.13, 1.14, 1.14, 1.13, 1.13, 1.30, 1.16, 1.28, 1.27, 1.30, 1.30, 1.33, 1.36, 1.39, 0.00, 0.00,
1.08, 1.08, 1.09, 1.08, 1.09, 1.08, 1.08, 1.08, 1.09, 1.09, 1.15, 1.09, 1.10, 1.10, 1.10, 1.09, 1.11, 1.24, 1.11, 1.21, 1.21, 1.23, 1.24, 1.26, 1.31, 1.35, 1.22, 0.00,
1.06, 1.06, 1.06, 1.07, 1.07, 1.07, 1.07, 1.07, 1.06, 1.07, 1.07, 1.07, 1.07, 1.08, 1.07, 1.08, 1.08, 1.21, 1.09, 1.17, 1.18, 1.20, 1.22, 1.24, 1.29, 1.29, 1.16, 1.09,
1.04, 1.05, 1.05, 1.05, 1.05, 1.05, 1.04, 1.05, 1.05, 1.05, 1.06, 1.05, 1.06, 1.06, 1.06, 1.07, 1.05, 1.17, 1.09, 1.15, 1.16, 1.18, 1.19, 1.22, 1.26, 1.26, 1.15, 1.08,
1.03, 1.03, 1.04, 1.04, 1.04, 1.04, 1.05, 1.04, 1.04, 1.04, 1.04, 1.05, 1.05, 1.05, 1.06, 1.05, 1.06, 1.14, 1.07, 1.13, 1.13, 1.17, 1.19, 1.21, 1.24, 1.29, 1.14, 1.07,
1.03, 1.03, 1.03, 1.03, 1.03, 1.03, 1.04, 1.03, 1.04, 1.03, 1.04, 1.04, 1.04, 1.04, 1.04, 1.03, 1.05, 1.13, 1.06, 1.12, 1.13, 1.16, 1.18, 1.20, 1.23, 1.26, 1.09, 1.06,
1.03, 1.03, 1.03, 1.03, 1.03, 1.03, 1.03, 1.03, 1.03, 1.03, 1.02, 1.03, 1.04, 1.04, 1.04, 1.04, 1.04, 1.11, 1.06, 1.10, 1.13, 1.14, 1.16, 1.19, 1.22, 1.23, 1.07, 1.05,
1.02, 1.02, 1.02, 1.02, 1.03, 1.03, 1.03, 1.03, 1.03, 1.03, 1.03, 1.03, 1.03, 1.03, 1.03, 1.03, 1.02, 1.11, 1.05, 1.10, 1.12, 1.14, 1.15, 1.18, 1.21, 1.25, 1.07, 1.05,
1.02, 1.02, 1.02, 1.02, 1.02, 1.02, 1.02, 1.02, 1.03, 1.02, 1.03, 1.03, 1.03, 1.03, 1.03, 1.03, 1.03, 1.10, 1.05, 1.09, 1.11, 1.13, 1.15, 1.17, 1.20, 1.23, 1.06, 1.05,
1.02, 1.02, 1.02, 1.02, 1.02, 1.02, 1.02, 1.02, 1.02, 1.02, 1.02, 1.02, 1.03, 1.03, 1.03, 1.03, 1.03, 1.10, 1.05, 1.10, 1.10, 1.13, 1.14, 1.18, 1.20, 1.23, 1.09, 1.05,
1.01, 1.01, 1.01, 1.01, 1.01, 1.02, 1.02, 1.02, 1.02, 1.01, 1.02, 1.02, 1.03, 1.03, 1.03, 1.03, 1.03, 1.08, 1.05, 1.09, 1.10, 1.12, 1.14, 1.17, 1.19, 1.21, 1.05, 1.05,
1.01, 1.00, 1.00, 1.01, 1.01, 1.02, 1.01, 1.01, 1.01, 1.01, 1.01, 1.02, 1.02, 1.02, 1.02, 1.02, 1.02, 1.07, 1.05, 1.09, 1.09, 1.11, 1.13, 1.15, 1.19, 1.14, 1.00, 1.01,
1.00, 1.00, 1.00, 1.00, 1.01, 1.02, 1.01, 1.01, 1.01, 1.01, 1.01, 1.01, 1.02, 1.02, 1.02, 1.02, 1.02, 1.06, 1.04, 1.08, 1.09, 1.11, 1.13, 1.13, 1.18, 1.21, 1.00, 1.01
])
caloStage2Params.layer1HCalScaleETBins = cms.vint32([6, 9, 12, 15, 20, 25, 30, 35, 40, 45, 55, 70, 256])
caloStage2Params.layer1HCalScaleFactors = cms.vdouble([
1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00
])
caloStage2Params.layer1HFScaleETBins = cms.vint32([6, 9, 12, 15, 20, 25, 30, 35, 40, 45, 55, 70, 256])
caloStage2Params.layer1HFScaleFactors = cms.vdouble([
1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00
])