L1TCalorimeter/python/caloStage2Params_2016_v2_1_cfi.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 from L1Trigger.L1TCalorimeter.caloParams_cfi import caloParamsSource
0004 import L1Trigger.L1TCalorimeter.caloParams_cfi
0005 caloStage2Params = L1Trigger.L1TCalorimeter.caloParams_cfi.caloParams.clone()
0006
0007 # towers
0008 caloStage2Params.towerLsbH        = cms.double(0.5)
0009 caloStage2Params.towerLsbE        = cms.double(0.5)
0010 caloStage2Params.towerLsbSum      = cms.double(0.5)
0011 caloStage2Params.towerNBitsH      = cms.int32(8)
0012 caloStage2Params.towerNBitsE      = cms.int32(8)
0013 caloStage2Params.towerNBitsSum    = cms.int32(9)
0014 caloStage2Params.towerNBitsRatio  = cms.int32(3)
0015 caloStage2Params.towerEncoding    = cms.bool(True)
0016
0017 # regions
0018 caloStage2Params.regionLsb        = cms.double(0.5)
0019 caloStage2Params.regionPUSType    = cms.string("None")
0020 caloStage2Params.regionPUSParams  = cms.vdouble()
0021
0022 # EG
0023 caloStage2Params.egLsb                      = cms.double(0.5)
0024 caloStage2Params.egSeedThreshold            = cms.double(2.)
0025 caloStage2Params.egNeighbourThreshold       = cms.double(1.)
0026 caloStage2Params.egHcalThreshold            = cms.double(0.)
0027 caloStage2Params.egTrimmingLUTFile          = cms.FileInPath("L1Trigger/L1TCalorimeter/data/egTrimmingLUT_10_v16.01.19.txt")
0028 caloStage2Params.egMaxHcalEt                = cms.double(0.)
0029 caloStage2Params.egMaxPtHOverE          = cms.double(128.)
0030 caloStage2Params.egMaxHOverELUTFile         = cms.FileInPath("L1Trigger/L1TCalorimeter/data/HoverEIdentification_0.995_v15.12.23.txt")
0031 caloStage2Params.egCompressShapesLUTFile    = cms.FileInPath("L1Trigger/L1TCalorimeter/data/egCompressLUT_v4.txt")
0032 caloStage2Params.egShapeIdType              = cms.string("compressed")
0033 caloStage2Params.egShapeIdVersion           = cms.uint32(0)
0034 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
0035 caloStage2Params.egPUSType                  = cms.string("None")
0036 caloStage2Params.egIsolationType            = cms.string("compressed")
0037 caloStage2Params.egIsoLUTFile               = cms.FileInPath("L1Trigger/L1TCalorimeter/data/IsoIdentification_0.25_adapt_extrap_FW_v16.04.05.txt")
0038 caloStage2Params.egIsoAreaNrTowersEta       = cms.uint32(2)
0039 caloStage2Params.egIsoAreaNrTowersPhi       = cms.uint32(4)
0040 caloStage2Params.egIsoVetoNrTowersPhi       = cms.uint32(2)
0041 #caloStage2Params.egIsoPUEstTowerGranularity = cms.uint32(1)
0042 #caloStage2Params.egIsoMaxEtaAbsForTowerSum  = cms.uint32(4)
0043 #caloStage2Params.egIsoMaxEtaAbsForIsoSum    = cms.uint32(27)
0044 caloStage2Params.egPUSParams                = cms.vdouble(1,4,32) #Isolation window in firmware goes up to abs(ieta)=32 for now
0045 caloStage2Params.egCalibrationType          = cms.string("compressed")
0046 caloStage2Params.egCalibrationVersion       = cms.uint32(0)
0047 caloStage2Params.egCalibrationLUTFile       = cms.FileInPath("L1Trigger/L1TCalorimeter/data/corrections_Trimming10_compressedieta_compressedE_compressedshape_v16.03.14_shapeIdentification_adapt0.99_compressedieta_compressedE_compressedshape_v15.12.08.txt")
0048 #caloStage2Params.egCalibrationLUTFile       = cms.FileInPath("L1Trigger/L1TCalorimeter/data/corrections_Trimming10_compressedieta_compressedE_compressedshape_v16.03.14.txt")
0049
0050 # Tau
0051 caloStage2Params.tauLsb                        = cms.double(0.5)
0052 caloStage2Params.tauSeedThreshold              = cms.double(0.)
0053 caloStage2Params.tauNeighbourThreshold         = cms.double(0.)
0054 caloStage2Params.tauIsoAreaNrTowersEta         = cms.uint32(2)
0055 caloStage2Params.tauIsoAreaNrTowersPhi         = cms.uint32(4)
0056 caloStage2Params.tauIsoVetoNrTowersPhi         = cms.uint32(2)
0057 caloStage2Params.tauPUSType                    = cms.string("None")
0058 caloStage2Params.tauIsoLUTFile                 = cms.FileInPath("L1Trigger/L1TCalorimeter/data/Tau_Iso_LUT_Option_5_Layer1Calibration_noCompressionBlock_FW_v4.0.0.txt")
0059 caloStage2Params.tauIsoLUTFile2                = cms.FileInPath("L1Trigger/L1TCalorimeter/data/Tau_Iso_LUT_Option_5_Layer1Calibration_noCompressionBlock_FW_v4.0.0.txt")
0060 caloStage2Params.tauCalibrationLUTFile         = cms.FileInPath("L1Trigger/L1TCalorimeter/data/Tau_Calibration_LUT_Layer1Calibration_FW_v9.0.0.txt")
0061 caloStage2Params.tauCompressLUTFile            = cms.FileInPath("L1Trigger/L1TCalorimeter/data/tauCompressAllLUT_12bit_v3.txt")
0062 caloStage2Params.tauPUSParams                  = cms.vdouble(1,4,32)
0063
0064 # jets
0065 caloStage2Params.jetLsb                = cms.double(0.5)
0066 caloStage2Params.jetSeedThreshold      = cms.double(4.0)
0067 caloStage2Params.jetNeighbourThreshold = cms.double(0.)
0068 caloStage2Params.jetPUSType            = cms.string("ChunkyDonut")
0069
0070 # Calibration options
0071 # function6PtParams22EtaBins or None
0072 #caloStage2Params.jetCalibrationType    = cms.string("None")
0073 #caloStage2Params.jetCalibrationType = cms.string("function8PtParams22EtaBins")
0074 caloStage2Params.jetCalibrationType = cms.string("LUT")
0075
0076 #Vector with 6 parameters for eta bin, from low eta to high
0077 # 1,0,1,0,1,1 gives no correction
0078 # must be in this form as may require > 255 arguments
0079
0080 # Or vector with 8 parameters, which caps correction value below given pT
0081 # as 6 parameters, but last two are max correction and L1 pT below which cap is applied, respectively
0082
0083 jetCalibParamsVector = cms.vdouble()
0084 jetCalibParamsVector.extend([
0085         1,0,1,0,1,1,1.36123039014,1024,
0086         1,0,1,0,1,1,1.37830172245,1024,
0087         1,0,1,0,1,1,1.37157036457,1024,
0088         1,0,1,0,1,1,1.42460009989,1024,
0089         10.1179757811,-697.422255848,55.9767511168,599.040770412,0.00930772659892,-21.9921521313,1.77585386314,24.1202894336,
0090         12.2578170485,-736.96846599,45.3225355911,848.976802835,0.00946235693865,-21.7970133915,2.04623980351,19.6049149791,
0091         14.0198255047,-769.175319944,38.687351315,1072.9785137,0.00951954709279,-21.6277409602,2.08021511285,22.265051562,
0092         14.119589176,-766.199501821,38.7767169666,1059.63374337,0.00952979125289,-21.6477483043,2.05901166216,23.8125466978,
0093         13.7594864391,-761.860391454,39.9060363401,1019.30588542,0.00952105483129,-21.6814176696,2.03808638982,22.2127275989,
0094         10.2635352836,-466.890522023,32.5408463829,2429.03382746,0.0111274121697,-22.0890253377,2.04880080215,22.5083699943,
0095         5.46086027683,-150.888778124,18.3292242153,16968.6469599,0.0147496053457,-22.4089831889,2.08107691501,22.4129703515,
0096         5.46086027683,-150.888778124,18.3292242153,16968.6469599,0.0147496053457,-22.4089831889,2.08107691501,22.4129703515,
0097         10.2635352836,-466.890522023,32.5408463829,2429.03382746,0.0111274121697,-22.0890253377,2.04880080215,22.5083699943,
0098         13.7594864391,-761.860391454,39.9060363401,1019.30588542,0.00952105483129,-21.6814176696,2.03808638982,22.2127275989,
0099         14.119589176,-766.199501821,38.7767169666,1059.63374337,0.00952979125289,-21.6477483043,2.05901166216,23.8125466978,
0100         14.0198255047,-769.175319944,38.687351315,1072.9785137,0.00951954709279,-21.6277409602,2.08021511285,22.265051562,
0101         12.2578170485,-736.96846599,45.3225355911,848.976802835,0.00946235693865,-21.7970133915,2.04623980351,19.6049149791,
0102         10.1179757811,-697.422255848,55.9767511168,599.040770412,0.00930772659892,-21.9921521313,1.77585386314,24.1202894336,
0103         1,0,1,0,1,1,1.42460009989,1024,
0104         1,0,1,0,1,1,1.37157036457,1024,
0105         1,0,1,0,1,1,1.37830172245,1024,
0106         1,0,1,0,1,1,1.36123039014,1024
0107 ])
0108 caloStage2Params.jetCalibrationParams  = jetCalibParamsVector
0109
0110 caloStage2Params.jetCompressPtLUTFile     = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_pt_compress.txt")
0111 caloStage2Params.jetCompressEtaLUTFile    = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_eta_compress.txt")
0112 caloStage2Params.jetCalibrationLUTFile    = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_add_mult.txt")
0113
0114
0115 # sums: 0=ET, 1=HT, 2=MET, 3=MHT
0116 caloStage2Params.etSumLsb                = cms.double(0.5)
0117 caloStage2Params.etSumEtaMin             = cms.vint32(1, 1, 1, 1)
0118 caloStage2Params.etSumEtaMax             = cms.vint32(28,  28,  28,  28)
0119 caloStage2Params.etSumEtThreshold        = cms.vdouble(0.,  30.,  0.,  30.)
0120
0121 caloStage2Params.etSumXCalibrationLUTFile         = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_etSumPUS_dummy.txt")
0122 caloStage2Params.etSumYCalibrationLUTFile         = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_etSumPUS_dummy.txt")
0123 caloStage2Params.etSumEttCalibrationLUTFile       = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_etSumPUS_dummy.txt")
0124 caloStage2Params.etSumEcalSumCalibrationLUTFile   = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_etSumPUS_dummy.txt")
0125
0126
0127 # Layer 1 LUT specification
0128 #
0129 # Et-dependent scale factors
0130 # ECal/HCal scale factors will be a 9*28 array:
0131 #   28 eta scale factors (1-28)
0132 #   in 9 ET bins (10, 15, 20, 25, 30, 35, 40, 45, Max)
0133 #  So, index = etBin*28+ieta
0134 caloStage2Params.layer1ECalScaleETBins = cms.vint32([10, 15, 20, 25, 30, 35, 40, 45, 256])
0135 caloStage2Params.layer1ECalScaleFactors = cms.vdouble([
0136     1.1847, 1.16759, 1.17779, 1.19955, 1.21125, 1.214, 1.21503, 1.22515, 1.24151, 1.27836, 1.30292, 1.33526, 1.42338, 1.4931, 1.49597, 1.50405, 1.52785, 1.81552, 1.59856, 1.75692, 1.76496, 1.77562, 1.69527, 1.66827, 1.61861, 1.56645, 1, 1,
0137     1.1351, 1.12589, 1.12834, 1.13725, 1.14408, 1.1494, 1.14296, 1.14852, 1.1578, 1.17634, 1.18038, 1.19386, 1.23758, 1.27605, 1.27818, 1.28195, 1.34881, 1.71053, 1.37338, 1.52571, 1.54801, 1.53316, 1.4397, 1.40497, 1.37743, 1.33914, 1, 1,
0138     1.18043, 1.17823, 1.1751, 1.17608, 1.19152, 1.196, 1.20125, 1.2068, 1.22584, 1.22476, 1.22395, 1.22302, 1.25137, 1.28097, 1.29871, 1.2862, 1.33489, 1.60937, 1.28365, 1.41367, 1.42521, 1.42041, 1.36784, 1.34922, 1.32754, 1.29825, 1, 1,
0139     1.11664, 1.11852, 1.11861, 1.12367, 1.12405, 1.14814, 1.14304, 1.15337, 1.16607, 1.18698, 1.17048, 1.17463, 1.2185, 1.23842, 1.23214, 1.24744, 1.30047, 1.47152, 1.22868, 1.33121, 1.34841, 1.35178, 1.30048, 1.28537, 1.27012, 1.24159, 1, 1,
0140     1.08422, 1.08146, 1.08706, 1.08906, 1.08636, 1.10092, 1.10363, 1.11102, 1.1186, 1.13301, 1.12369, 1.14377, 1.16477, 1.17801, 1.18782, 1.17168, 1.24593, 1.36835, 1.20252, 1.28349, 1.29828, 1.30328, 1.26848, 1.25817, 1.2464, 1.22259, 1, 1,
0141     1.07444, 1.06774, 1.06883, 1.0707, 1.07881, 1.08859, 1.08285, 1.08747, 1.09736, 1.10678, 1.10008, 1.10717, 1.12858, 1.15383, 1.15826, 1.14855, 1.19911, 1.32567, 1.17553, 1.25976, 1.27926, 1.28459, 1.24524, 1.23706, 1.22597, 1.20006, 1, 1,
0142     1.06224, 1.05968, 1.05767, 1.06254, 1.06729, 1.0691, 1.07125, 1.07312, 1.08124, 1.08966, 1.08695, 1.08826, 1.10611, 1.13115, 1.12641, 1.13093, 1.17074, 1.28958, 1.16217, 1.22844, 1.24812, 1.25352, 1.22065, 1.21287, 1.20544, 1.18344, 1, 1,
0143     1.03589, 1.03224, 1.03229, 1.03623, 1.03979, 1.04403, 1.04574, 1.049, 1.04821, 1.06183, 1.0588, 1.06655, 1.08582, 1.10289, 1.10052, 1.10506, 1.143, 1.27373, 1.1459, 1.2156, 1.23455, 1.23968, 1.20753, 1.20127, 1.19629, 1.16809, 1, 1,
0144     1.03456, 1.02955, 1.03079, 1.03509, 1.03949, 1.0437, 1.04236, 1.04486, 1.0517, 1.05864, 1.05516, 1.06167, 1.07738, 1.0985, 1.09317, 1.09559, 1.13557, 1.26076, 1.14118, 1.20545, 1.22137, 1.22802, 1.19936, 1.19676, 1.19088, 1.16709, 1, 1,
0145     ])
0146 caloStage2Params.layer1HCalScaleETBins = cms.vint32([10, 15, 20, 25, 30, 35, 40, 45, 256])
0147 caloStage2Params.layer1HCalScaleFactors = cms.vdouble([
0148     1.511112, 1.519900, 1.499483, 1.488560, 1.528111, 1.475114, 1.476616, 1.514163, 1.515306, 1.542464, 1.511663, 1.593745, 1.493667, 1.485315, 1.419925, 1.349169, 1.312518, 1.423302, 1.478461, 1.525868, 1, 1, 1, 1, 1, 1, 1, 1,
0149     1.383350, 1.365700, 1.368470, 1.354610, 1.348480, 1.329720, 1.272250, 1.301710, 1.322210, 1.360860, 1.333850, 1.392200, 1.403060, 1.394870, 1.322050, 1.244570, 1.206910, 1.321870, 1.344160, 1.403270, 1, 1, 1, 1, 1, 1, 1, 1,
0150     1.245690, 1.238320, 1.245420, 1.234830, 1.243730, 1.249790, 1.179450, 1.213620, 1.219030, 1.252130, 1.209560, 1.250710, 1.280490, 1.262800, 1.254060, 1.186810, 1.127830, 1.260000, 1.275140, 1.305850, 1, 1, 1, 1, 1, 1, 1, 1,
0151     1.189940, 1.189120, 1.177120, 1.179690, 1.185510, 1.150590, 1.151830, 1.167860, 1.154310, 1.163190, 1.161700, 1.136100, 1.161870, 1.195050, 1.153910, 1.117900, 1.106750, 1.208120, 1.160020, 1.232800, 1, 1, 1, 1, 1, 1, 1, 1,
0152     1.122540, 1.129520, 1.125080, 1.115150, 1.118250, 1.096190, 1.108170, 1.087490, 1.109750, 1.099780, 1.081000, 1.050610, 1.078270, 1.079460, 1.047740, 1.041400, 1.041750, 1.116880, 1.097730, 1.125780, 1, 1, 1, 1, 1, 1, 1, 1,
0153     1.110470, 1.117340, 1.115980, 1.088490, 1.088260, 1.078230, 1.062720, 1.054690, 1.053270, 1.086640, 1.050620, 1.038470, 1.046440, 1.059130, 1.012240, 1.039030, 1.036040, 1.088460, 1.078880, 1.090600, 1, 1, 1, 1, 1, 1, 1, 1,
0154     1.115970, 1.111010, 1.113170, 1.079390, 1.076850, 1.063730, 1.039300, 1.049910, 1.040100, 1.025820, 1.015830, 1.015850, 1.010810, 1.014210, 0.980321, 1.023580, 1.045990, 1.073220, 1.057750, 1.059850, 1, 1, 1, 1, 1, 1, 1, 1,
0155     1.061180, 1.059770, 1.071210, 1.064420, 1.065340, 1.043070, 1.041400, 1.022680, 1.017410, 1.017690, 1.005610, 1.006360, 0.999420, 0.990866, 0.986723, 0.989036, 0.995116, 1.045620, 1.024330, 1.040660, 1, 1, 1, 1, 1, 1, 1, 1,
0156     1.083150, 1.067090, 1.083180, 1.061010, 1.075640, 1.051640, 1.038760, 1.042670, 1.010910, 1.011580, 1.006560, 0.984468, 0.986642, 0.985799, 0.968133, 1.000290, 1.011210, 1.046690, 1.016670, 1.020470, 1, 1, 1, 1, 1, 1, 1, 1,
0157     ])
0158 # HF 1x1 scale factors will be a 5*12 array:
0159 #  12 eta scale factors (30-41)
0160 #  in 5 REAL ET bins (5, 20, 30, 50, Max)
0161 #  So, index = etBin*12+ietaHF
0162 caloStage2Params.layer1HFScaleETBins = cms.vint32([5, 20, 30, 50, 256])
0163 caloStage2Params.layer1HFScaleFactors = cms.vdouble([
0164     1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
0165     1.55, 1.49, 1.35, 1.29, 1.30, 1.42, 1.49, 1.59, 1.74, 1.86, 2.18, 2.43,
0166     1.60, 1.51, 1.38, 1.38, 1.44, 1.56, 1.60, 1.67, 1.83, 2.02, 2.66, 2.79,
0167     1.56, 1.41, 1.34, 1.42, 1.44, 1.52, 1.57, 1.63, 1.73, 1.94, 2.64, 2.64,
0168     1.46, 1.30, 1.29, 1.43, 1.42, 1.49, 1.52, 1.59, 1.69, 1.87, 2.49, 2.66,
0169     ])