-
Notifications
You must be signed in to change notification settings - Fork 0
/
caloStage2Params_2016_v0.0_cfi.py
122 lines (109 loc) · 6.37 KB
/
caloStage2Params_2016_v0.0_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
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.egSeedThreshold = cms.double(2.)
caloStage2Params.egNeighbourThreshold = cms.double(1.)
caloStage2Params.egHcalThreshold = cms.double(0.)
caloStage2Params.egTrimmingLUTFile = cms.FileInPath("L1Trigger/L1TCalorimeter/data/egTrimmingLUT_5.txt")
caloStage2Params.egMaxHcalEt = cms.double(0.)
caloStage2Params.egMaxPtHOverE = cms.double(128.)
caloStage2Params.egMaxHOverELUTFile = cms.FileInPath("L1Trigger/L1TCalorimeter/data/egMaxHOverELUT_995eff.txt")
caloStage2Params.egCompressShapesLUTFile = cms.FileInPath("L1Trigger/L1TCalorimeter/data/egCompressLUT_v1.txt")
caloStage2Params.egShapeIdType = cms.string("compressed")
caloStage2Params.egShapeIdVersion = cms.uint32(0)
caloStage2Params.egShapeIdLUTFile = cms.FileInPath("L1Trigger/L1TCalorimeter/data/shapeIdentification_0.99_compressedieta_compressedE_compressedshape_v15.10.20.txt")
caloStage2Params.egPUSType = cms.string("None")
caloStage2Params.egIsoLUTFile = cms.FileInPath("L1Trigger/L1TCalorimeter/data/egIsoLUTPU40bx25NrRings4Eff95.txt")
caloStage2Params.egIsoAreaNrTowersEta = cms.uint32(2)
caloStage2Params.egIsoAreaNrTowersPhi = cms.uint32(4)
caloStage2Params.egIsoVetoNrTowersPhi = cms.uint32(3)
#caloStage2Params.egIsoPUEstTowerGranularity = cms.uint32(1)
#caloStage2Params.egIsoMaxEtaAbsForTowerSum = cms.uint32(4)
#caloStage2Params.egIsoMaxEtaAbsForIsoSum = cms.uint32(27)
caloStage2Params.egPUSParams = cms.vdouble(1,4,27)
caloStage2Params.egCalibrationType = cms.string("compressed")
caloStage2Params.egCalibrationVersion = cms.uint32(0)
caloStage2Params.egCalibrationLUTFile = cms.FileInPath("L1Trigger/L1TCalorimeter/data/corrections_Trimming5_PU40bx25_compressedieta_compressedE_compressedshape_v15.10.20.txt")
# Tau
caloStage2Params.tauLsb = cms.double(0.5)
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/tauIsoLUTetPU.txt")
caloStage2Params.tauCalibrationLUTFile = cms.FileInPath("L1Trigger/L1TCalorimeter/data/tauCalibrationLUT.txt")
caloStage2Params.tauPUSParams = cms.vdouble(1,4,27)
# jets
caloStage2Params.jetLsb = cms.double(0.5)
caloStage2Params.jetSeedThreshold = cms.double(1.5)
caloStage2Params.jetNeighbourThreshold = cms.double(0.)
caloStage2Params.jetPUSType = cms.string("ChunkyDonut")
#Calibration options
# e.g. function6PtParams22EtaBins
caloStage2Params.jetCalibrationType = cms.string("None")
#caloStage2Params.jetCalibrationType = cms.string("function6PtParams22EtaBins")
#Vector with 6 parameters for eta bin, from low eta to high
# 1,0,1,0,1,1 gives no correction
# must be in this form as may require > 255 arguments
jetCalibParamsVector = cms.vdouble()
jetCalibParamsVector.extend([
1.52819258712,0,1,0,1,1, # Constant calibrations in HF bins
1.32880934866,0,1,0,1,1,
1.35650073562,0,1,0,1,1,
1.55766262213,0,1,0,1,1,
12.42641269,6.35085569,1.21845374,-14.70772743,0.00118945,-10.56931753,
12.14768583,8.05320174,0.92028575,-15.48969055,0.00155426,-10.48022143,
11.87301715,10.45078171,1.02636603,-16.87586292,0.00210878,-10.34892982,
12.09471152,9.36873133,0.96484984,-16.09937338,0.00177759,-10.43545594,
4.69610083,11.50217825,1.31071693,-23.88756705,0.00336036,-19.28130829,
7.08667633,23.68511274,2.52417732,-43.09154102,0.00359229,-18.80055010,
3.44975454,28.42124526,2.58428891,-137.42027137,0.00669104,-19.53448313,
3.44975454,28.42124526,2.58428891,-137.42027137,0.00669104,-19.53448313,
7.08667633,23.68511274,2.52417732,-43.09154102,0.00359229,-18.80055010,
4.69610083,11.50217825,1.31071693,-23.88756705,0.00336036,-19.28130829,
12.09471152,9.36873133,0.96484984,-16.09937338,0.00177759,-10.43545594,
11.87301715,10.45078171,1.02636603,-16.87586292,0.00210878,-10.34892982,
12.14768583,8.05320174,0.92028575,-15.48969055,0.00155426,-10.48022143,
12.42641269,6.35085569,1.21845374,-14.70772743,0.00118945,-10.56931753,
1.55766262213,0,1,0,1,1, # Constant calibrations in HF bins
1.35650073562,0,1,0,1,1,
1.32880934866,0,1,0,1,1,
1.52819258712,0,1,0,1,1
])
caloStage2Params.jetCalibrationParams = jetCalibParamsVector
# sums: 0=ET, 1=HT, 2=MET, 3=MHT
caloStage2Params.etSumLsb = cms.double(0.5)
caloStage2Params.etSumEtaMin = cms.vint32(0, 0, 0, 0)
caloStage2Params.etSumEtaMax = cms.vint32(40, 36, 40, 36)
caloStage2Params.etSumEtThreshold = cms.vdouble(0., 0., 0., 0.)
# Layer 1 LUT specification
#
# Et-dependent scale factors
# ECal/HCal scale factors will be a 9*28 array:
# 28 eta scale factors (1-28)
# in 9 ET bins (10, 15, 20, 25, 30, 35, 40, 45, Max)
# So, index = etBin*28+ieta
caloStage2Params.layer1ECalScaleETBins = cms.vint32([1])
caloStage2Params.layer1ECalScaleFactors = cms.vdouble([1.]*28)
caloStage2Params.layer1HCalScaleETBins = cms.vint32([1])
caloStage2Params.layer1HCalScaleFactors = cms.vdouble([1.]*28)
caloStage2Params.layer1HFScaleETBins = cms.vint32([1])
caloStage2Params.layer1HFScaleFactors = cms.vdouble([1.]*12)