Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2023-10-25 09:54:53

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")
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_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.txt")
0048 
0049 # Tau
0050 caloStage2Params.tauLsb                        = cms.double(0.5)
0051 caloStage2Params.tauSeedThreshold              = cms.double(0.)
0052 caloStage2Params.tauNeighbourThreshold         = cms.double(0.)
0053 caloStage2Params.tauIsoAreaNrTowersEta         = cms.uint32(2)
0054 caloStage2Params.tauIsoAreaNrTowersPhi         = cms.uint32(4)
0055 caloStage2Params.tauIsoVetoNrTowersPhi         = cms.uint32(2)
0056 caloStage2Params.tauPUSType                    = cms.string("None")
0057 caloStage2Params.tauIsoLUTFile                 = cms.FileInPath("L1Trigger/L1TCalorimeter/data/Tau_Iso_LUT_Option_21_Layer1Calibration_noCompressionBlock_v4.0.0.txt")
0058 caloStage2Params.tauIsoLUTFile2                = cms.FileInPath("L1Trigger/L1TCalorimeter/data/Tau_Iso_LUT_Option_21_Layer1Calibration_noCompressionBlock_v4.0.0.txt")
0059 caloStage2Params.tauCalibrationLUTFile         = cms.FileInPath("L1Trigger/L1TCalorimeter/data/Tau_Calibration_LUT_Layer1Calibration_v9.0.0.txt")
0060 caloStage2Params.tauCompressLUTFile            = cms.FileInPath("L1Trigger/L1TCalorimeter/data/tauCompressAllLUT_12bit_v3.txt")
0061 caloStage2Params.tauPUSParams                  = cms.vdouble(1,4,32)
0062 
0063 # jets
0064 caloStage2Params.jetLsb                = cms.double(0.5)
0065 caloStage2Params.jetSeedThreshold      = cms.double(4.0)
0066 caloStage2Params.jetNeighbourThreshold = cms.double(0.)
0067 caloStage2Params.jetPUSType            = cms.string("ChunkyDonut")
0068 
0069 # Calibration options
0070 # function6PtParams22EtaBins or None
0071 # caloStage2Params.jetCalibrationType    = cms.string("None")
0072 # caloStage2Params.jetCalibrationType = cms.string("function8PtParams22EtaBins")
0073 caloStage2Params.jetCalibrationType = cms.string("functionErf11PtParams16EtaBins")
0074 # caloStage2Params.jetCalibrationType = cms.string("None")
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 # this vector corresponds to "function8PtParams22EtaBins"
0109 # caloStage2Params.jetCalibrationParams  = jetCalibParamsVector 
0110 
0111 
0112 # vector with 11 parameters for each eta bin
0113 # each eta bin represented by a new line, starting at first eta bin
0114 # the first seven parameters of a line are the 7 parameters for the correction function
0115 # the next four parameters are:
0116 # 1. the function value at the minimum pt of fit
0117 # 2. the minimum pt of the fit
0118 # 3. the function value at the maximum pt of the fit
0119 # 4. the maximum pt of the fit
0120 jetCalibParamsVector16 = cms.vdouble()
0121 jetCalibParamsVector16.extend([ 
0122         0.661201, -1.08715e+06, -2.59519e-08, -10.8044, -1.04431e-06, -1.358, 0.830229, 1.9517, 21.7232, 1.11497, 308.136, 
0123         0.595461, -242991, -5.38e-08, -31.5184, -8.65299e-06, -0.963764, 0.382539, 2.0419, 22.4169, 1.13073, 304.49, 
0124         2.55569, -1.76867e+06, 6.783e-08, -8.05711, -3.97612e-07, -2.83151, 1.1268, 2.00221, 20.7898, 1.14923, 248.034, 
0125         2.96018, -4.20761e+06, 7.63023e-08, -2.55568, -9.93322e-08, -6.08005, 1.29864, 2.04699, 17.3628, 1.20751, 193.13, 
0126         3.08393, -4.49833e+06, 8.14721e-08, -2.23442, -8.46574e-08, -7.22304, 1.32952, 2.03851, 21.509, 1.19884, 211.281, 
0127         2.1004, -2.69476e+06, 4.91209e-08, -4.5254, -2.40174e-07, -3.2305, 1.21127, 1.97181, 16.7042, 1.06621, 273.416, 
0128         2.14597, -2.8508e+06, 5.17083e-08, -4.21701, -2.28694e-07, -3.40691, 1.18177, 1.98222, 15.4658, 1.05329, 242.235, 
0129         2.35967, -3.55304e+06, 6.32171e-08, -2.95984, -1.3141e-07, -4.07551, 1.21547, 1.83235, 15.7161, 1.03703, 190.26, 
0130         0.856007, -364884, -1.29766e-08, -22.9827, -2.11044e-06, -2.01542, 1.02893, 1.76636, 18.165, 1.00783, 273.301, 
0131         3.21455, -286314, 3.16865e-07, -19.3786, -1.14632e-06, -4.6408, 1.33054, 1.46451, 21.5884, 0.987992, 246.126, 
0132         -2.88427, -21731.7, -1.44395e-06, -106.655, -1.78369e-05, -3.42763, 1.28446, 1.37317, 21.5443, 0.983075, 288.211, 
0133         -0.787852, -12400.7, -8.08595e-07, -152.002, -2.55708e-05, -3.91554, 1.28934, 1.27682, 27.3396, 0.962841, 225.501, 
0134         -1.90211, -17583.4, -1.16672e-06, -119.776, -8.63159e-06, -7.56003, 1.45489, 1.06435, 36.2542, 0.929621, 112.707, 
0135         -74.6452, -258031, -1.76366e-05, -11.9507, -3.53722e-05, -0.493146, 0.593305, 1.11432, 46.5187, 1.03823, 1024, 
0136         0.568061, -7606.13, -5.39389e-07, -150.07, 4.85322e-05, -2.42215, 2.07552, 1.11686, 27.012, 0.967464, 1024, 
0137         -0.819561, -12223.5, -7.93953e-07, -149.116, -3.22883e-05, -5.31212, 1.26463, 1.23498, 24.8899, 0.840316, 186  
0138 ])
0139 # this vector corresponds to "functionErf11PtParams16EtaBins"
0140 caloStage2Params.jetCalibrationParams  = jetCalibParamsVector16
0141 
0142 
0143 caloStage2Params.jetCompressPtLUTFile     = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_pt_compress.txt")
0144 caloStage2Params.jetCompressEtaLUTFile    = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_eta_compress.txt")
0145 caloStage2Params.jetCalibrationLUTFile    = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_add_mult.txt")
0146 
0147 
0148 # sums: 0=ET, 1=HT, 2=MET, 3=MHT
0149 caloStage2Params.etSumLsb                = cms.double(0.5)
0150 caloStage2Params.etSumEtaMin             = cms.vint32(1, 1, 1, 1, 1)
0151 caloStage2Params.etSumEtaMax             = cms.vint32(28,  28,  28,  28, 28)
0152 caloStage2Params.etSumEtThreshold        = cms.vdouble(0.,  30.,  0.,  30., 0.)
0153 
0154 caloStage2Params.etSumXCalibrationLUTFile         = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_etSumPUS_dummy.txt")
0155 caloStage2Params.etSumYCalibrationLUTFile         = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_etSumPUS_dummy.txt")
0156 caloStage2Params.etSumEttCalibrationLUTFile       = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_etSumPUS_dummy.txt")
0157 caloStage2Params.etSumEcalSumCalibrationLUTFile   = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_etSumPUS_dummy.txt")
0158 
0159 
0160 # Layer 1 LUT specification
0161 #
0162 # Et-dependent scale factors
0163 # ECal/HCal scale factors will be a 9*28 array:
0164 #   28 eta scale factors (1-28)
0165 #   in 9 ET bins (10, 15, 20, 25, 30, 35, 40, 45, Max)
0166 #  So, index = etBin*28+ieta
0167 caloStage2Params.layer1ECalScaleETBins = cms.vint32([10, 15, 20, 25, 30, 35, 40, 45, 256])
0168 caloStage2Params.layer1ECalScaleFactors = cms.vdouble([
0169     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.56645, 1.56645,
0170     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.33914, 1.33914,
0171     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.29825, 1.29825,
0172     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.24159, 1.24159,
0173     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.22259, 1.22259,
0174     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.20006, 1.20006,
0175     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.18344, 1.18344,
0176     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.16809, 1.16809,
0177     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.16709, 1.16709,
0178     ])
0179 caloStage2Params.layer1HCalScaleETBins = cms.vint32([10, 15, 20, 25, 30, 35, 40, 45, 256])
0180 caloStage2Params.layer1HCalScaleFactors = cms.vdouble([
0181     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.525868, 1.525868, 1.525868, 1.525868, 1.525868, 1.525868, 1.525868, 1.525868,
0182     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.403270, 1.403270, 1.403270, 1.403270, 1.403270, 1.403270, 1.403270, 1.403270,
0183     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.305850, 1.305850, 1.305850, 1.305850, 1.305850, 1.305850, 1.305850, 1.305850,
0184     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.232800, 1.232800, 1.232800, 1.232800, 1.232800, 1.232800, 1.232800, 1.232800,
0185     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.125780, 1.125780, 1.125780, 1.125780, 1.125780, 1.125780, 1.125780, 1.125780,
0186     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.090600, 1.090600, 1.090600, 1.090600, 1.090600, 1.090600, 1.090600, 1.090600,
0187     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.059850, 1.059850, 1.059850, 1.059850, 1.059850, 1.059850, 1.059850, 1.059850,
0188     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.040660, 1.040660, 1.040660, 1.040660, 1.040660, 1.040660, 1.040660, 1.040660,
0189     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.020470, 1.020470, 1.020470, 1.020470, 1.020470, 1.020470, 1.020470, 1.020470,
0190     ])
0191 # HF 1x1 scale factors will be a 5*12 array:
0192 #  12 eta scale factors (30-41)
0193 #  in 5 REAL ET bins (5, 20, 30, 50, Max)
0194 #  So, index = etBin*12+ietaHF
0195 caloStage2Params.layer1HFScaleETBins = cms.vint32([5, 20, 30, 50, 256])
0196 caloStage2Params.layer1HFScaleFactors = cms.vdouble([
0197     1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 
0198     1.767080, 1.767080, 1.755186, 1.769951, 1.763527, 1.791043, 1.898787, 1.982235, 2.071074, 2.193011, 2.356886, 2.403384, 
0199     2.170477, 2.170477, 2.123540, 2.019866, 1.907698, 1.963179, 1.989122, 2.035251, 2.184642, 2.436399, 2.810884, 2.923750, 
0200     1.943941, 1.943941, 1.899826, 1.813950, 1.714978, 1.736184, 1.785928, 1.834211, 1.944230, 2.153565, 2.720887, 2.749795, 
0201     1.679984, 1.679984, 1.669753, 1.601871, 1.547276, 1.577805, 1.611497, 1.670184, 1.775022, 1.937061, 2.488311, 2.618629, 
0202     ])