Back to home page

Project CMSSW displayed by LXR

 
 

    


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

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_adapt_extrap_v16.07.29.txt")
0038 caloStage2Params.egIsoLUTFile               = cms.FileInPath("L1Trigger/L1TCalorimeter/data/IsolationLUT_v17.02.18.txt")
0039 caloStage2Params.egIsoAreaNrTowersEta       = cms.uint32(2)
0040 caloStage2Params.egIsoAreaNrTowersPhi       = cms.uint32(4)
0041 caloStage2Params.egIsoVetoNrTowersPhi       = cms.uint32(2)
0042 #caloStage2Params.egIsoPUEstTowerGranularity = cms.uint32(1)
0043 #caloStage2Params.egIsoMaxEtaAbsForTowerSum  = cms.uint32(4)
0044 #caloStage2Params.egIsoMaxEtaAbsForIsoSum    = cms.uint32(27)
0045 caloStage2Params.egPUSParams                = cms.vdouble(1,4,32) #Isolation window in firmware goes up to abs(ieta)=32 for now
0046 caloStage2Params.egCalibrationType          = cms.string("compressed")
0047 caloStage2Params.egCalibrationVersion       = cms.uint32(0)
0048 caloStage2Params.egCalibrationLUTFile       = cms.FileInPath("L1Trigger/L1TCalorimeter/data/corrections_Trimming10_compressedieta_compressedE_compressedshape_v17.02.18_shapeIdentification_adapt0.99_compressedieta_compressedE_compressedshape_v15.12.08.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_22_2017_FW_v8.0.0.txt")
0059 caloStage2Params.tauIsoLUTFile2                = cms.FileInPath("L1Trigger/L1TCalorimeter/data/Tau_Iso_LUT_Option_22_2017_FW_v8.0.0.txt")
0060 caloStage2Params.tauCalibrationLUTFile         = cms.FileInPath("L1Trigger/L1TCalorimeter/data/Tau_Calibration_LUT_2017_Layer1Calibration_FW_v11.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("functionErf11PtParams16EtaBins")
0075 # caloStage2Params.jetCalibrationType = cms.string("LUT")
0076 
0077 #Vector with 6 parameters for eta bin, from low eta to high
0078 # 1,0,1,0,1,1 gives no correction
0079 # must be in this form as may require > 255 arguments
0080 
0081 # Or vector with 8 parameters, which caps correction value below given pT
0082 # as 6 parameters, but last two are max correction and L1 pT below which cap is applied, respectively
0083 
0084 jetCalibParamsVector = cms.vdouble()
0085 jetCalibParamsVector.extend([
0086         1,0,1,0,1,1,1.36123039014,1024,
0087         1,0,1,0,1,1,1.37830172245,1024,
0088         1,0,1,0,1,1,1.37157036457,1024,
0089         1,0,1,0,1,1,1.42460009989,1024,
0090         10.1179757811,-697.422255848,55.9767511168,599.040770412,0.00930772659892,-21.9921521313,1.77585386314,24.1202894336,
0091         12.2578170485,-736.96846599,45.3225355911,848.976802835,0.00946235693865,-21.7970133915,2.04623980351,19.6049149791,
0092         14.0198255047,-769.175319944,38.687351315,1072.9785137,0.00951954709279,-21.6277409602,2.08021511285,22.265051562,
0093         14.119589176,-766.199501821,38.7767169666,1059.63374337,0.00952979125289,-21.6477483043,2.05901166216,23.8125466978,
0094         13.7594864391,-761.860391454,39.9060363401,1019.30588542,0.00952105483129,-21.6814176696,2.03808638982,22.2127275989,
0095         10.2635352836,-466.890522023,32.5408463829,2429.03382746,0.0111274121697,-22.0890253377,2.04880080215,22.5083699943,
0096         5.46086027683,-150.888778124,18.3292242153,16968.6469599,0.0147496053457,-22.4089831889,2.08107691501,22.4129703515,
0097         5.46086027683,-150.888778124,18.3292242153,16968.6469599,0.0147496053457,-22.4089831889,2.08107691501,22.4129703515,
0098         10.2635352836,-466.890522023,32.5408463829,2429.03382746,0.0111274121697,-22.0890253377,2.04880080215,22.5083699943,
0099         13.7594864391,-761.860391454,39.9060363401,1019.30588542,0.00952105483129,-21.6814176696,2.03808638982,22.2127275989,
0100         14.119589176,-766.199501821,38.7767169666,1059.63374337,0.00952979125289,-21.6477483043,2.05901166216,23.8125466978,
0101         14.0198255047,-769.175319944,38.687351315,1072.9785137,0.00951954709279,-21.6277409602,2.08021511285,22.265051562,
0102         12.2578170485,-736.96846599,45.3225355911,848.976802835,0.00946235693865,-21.7970133915,2.04623980351,19.6049149791,
0103         10.1179757811,-697.422255848,55.9767511168,599.040770412,0.00930772659892,-21.9921521313,1.77585386314,24.1202894336,
0104         1,0,1,0,1,1,1.42460009989,1024,
0105         1,0,1,0,1,1,1.37157036457,1024,
0106         1,0,1,0,1,1,1.37830172245,1024,
0107         1,0,1,0,1,1,1.36123039014,1024
0108 ])
0109 # this vector corresponds to "function8PtParams22EtaBins"
0110 # caloStage2Params.jetCalibrationParams  = jetCalibParamsVector  
0111 
0112 
0113 # vector with 11 parameters for each eta bin
0114 # each eta bin represented by a new line, starting at first eta bin
0115 # the first seven parameters of a line are the 7 parameters for the correction function
0116 # the next four parameters are:
0117 # 1. the function value at the minimum pt of fit
0118 # 2. the minimum pt of the fit
0119 # 3. the function value at the maximum pt of the fit
0120 # 4. the maximum pt of the fit
0121 jetCalibParamsVector16 = cms.vdouble()
0122 jetCalibParamsVector16.extend([ 
0123         1.2931, -1.398e+06, 6.09656e-08, 0.725755, -5.38716e-07, -3.40806, 1.31608, 2.08678, 19.587, 1.09283, 649.338, 
0124         1.31069, -1.46535e+06, 6.40139e-08, 0.903274, -5.38967e-07, -3.5372, 1.31318, 2.15636, 18.4244, 1.1101, 633.54, 
0125         1.29052, -1.23346e+06, 4.89144e-08, 0.754995, -6.997e-07, -2.71405, 1.223, 2.22599, 18.2578, 1.15074, 663.791, 
0126         1.19118, -748138, 2.46431e-09, -0.315317, -1.59196e-06, -1.74636, 1.01454, 2.22705, 25.0086, 1.18948, 646.903, 
0127         1.09348, -1.18135e+06, -6.02111e-09, -1.0158, -1.40543e-06, -1.39134, 0.792196, 2.12615, 28.6461, 1.13239, 574.203, 
0128         1.17443, -1.33589e+06, 7.75662e-08, 1.46473, -6.47983e-07, -3.13498, 1.17014, 2.04336, 23.6199, 1.03805, 437.174, 
0129         1.16636, -1.21175e+06, 6.18082e-08, 1.41557, -7.8438e-07, -2.09711, 1.05929, 2.0897, 21.5368, 1.05125, 630.619, 
0130         1.13949, -1.09162e+06, 5.5306e-08, 1.09139, -7.67904e-07, -1.93963, 1.13944, 2.05669, 17.394, 1.02673, 644.078, 
0131         1.12711, -1.08492e+06, 4.71761e-08, 0.873853, -8.86642e-07, -1.535, 1.00765, 2.07833, 18.748, 1.02498, 613.047, 
0132         0.3138, -456006, -3.49199e-06, 2.96106, -1.05403e-05, -0.307267, 0.491321, 1.74198, 22.5035, 1.06579, 276.679, 
0133         0.247104, -696099, -4.05576e-06, 3.80007, -1.61256e-05, -0.124935, -0.244893, 1.60903, 27.9737, 1.06, 277.597, 
0134         10.8027, -3.95478e+08, -2.38679e-09, 12.2808, -5.8088e-09, -0.579212, 0.950056, 1.50119, 30.322, 1.0339, 291.417, 
0135         -7.94334, -1.17259e+08, -1.23506e-10, 385.018, -7.00971e-08, 0.000459177, 34.9532, 1.16265, 61.5647, 0.926645, 209.196, 
0136         7.11338, -3.55228e+08, -1.14879e-09, 47.1441, -3.60311e-08, 0.0853058, 2.26173, 1.03403, 48.9575, 0.846677, 111.796, 
0137         7.18481, -3.9883e+08, -1.14597e-09, 47.9492, -3.80535e-08, 0.0267988, 2.70412, 1.01369, 34.0642, 0.835839, 113.372, 
0138         7.27504, -4.83419e+08, -1.14447e-09, 46.0531, -3.71605e-08, 0.0158836, 3.53306, 1.11904, 30.5118, 0.856453, 80.149,   
0139 ])
0140 # this vector corresponds to "functionErf11PtParams16EtaBins"
0141 caloStage2Params.jetCalibrationParams  = jetCalibParamsVector16
0142 
0143 
0144 caloStage2Params.jetCompressPtLUTFile     = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_pt_compress.txt")
0145 caloStage2Params.jetCompressEtaLUTFile    = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_30to40_hfHighPt_experiment2_changeLimits_eta.txt")
0146 caloStage2Params.jetCalibrationLUTFile    = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_30to40_hfHighPt_experiment2_changeLimits_add_mult.txt")
0147 
0148 
0149 # sums: 0=ET, 1=HT, 2=MET, 3=MHT
0150 caloStage2Params.etSumLsb                = cms.double(0.5)
0151 caloStage2Params.etSumEtaMin             = cms.vint32(1, 1, 1, 1, 1)
0152 caloStage2Params.etSumEtaMax             = cms.vint32(28,  28,  28,  28, 27)
0153 caloStage2Params.etSumEtThreshold        = cms.vdouble(0.,  30.,  0.,  30., 0.)
0154 
0155 caloStage2Params.etSumXCalibrationLUTFile         = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_etSumPUS_dummy.txt")
0156 caloStage2Params.etSumYCalibrationLUTFile         = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_etSumPUS_dummy.txt")
0157 caloStage2Params.etSumEttCalibrationLUTFile       = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_etSumPUS_dummy.txt")
0158 caloStage2Params.etSumEcalSumCalibrationLUTFile   = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_etSumPUS_dummy.txt")
0159 
0160 
0161 # Layer 1 LUT specification
0162 #
0163 # Et-dependent scale factors
0164 # ECal/HCal scale factors will be a 9*28 array:
0165 #   28 eta scale factors (1-28)
0166 #   in 9 ET bins (10, 15, 20, 25, 30, 35, 40, 45, Max)
0167 #  So, index = etBin*28+ieta
0168 caloStage2Params.layer1ECalScaleETBins = cms.vint32([10, 15, 20, 25, 30, 35, 40, 45, 256])
0169 caloStage2Params.layer1ECalScaleFactors = cms.vdouble([
0170  1.131385, 1.128532, 1.133259, 1.138950, 1.138118, 1.145239, 1.135720, 1.139045, 1.142051, 1.153916, 1.149703, 1.156583, 1.172763, 1.193265, 1.187774, 1.184272,
0171  1.217902, 1.437853, 1.278790, 1.409696, 1.411696, 1.440007, 1.421163, 1.464394, 1.474977, 1.466513, 1.466475, 1.419311, 1.087010, 1.081057, 1.085211, 1.088576, 
0172  1.087186, 1.095176, 1.090827, 1.092414, 1.094626, 1.104560, 1.101845, 1.110506, 1.121719, 1.143752, 1.142682, 1.147103, 1.193239, 1.574359, 1.279538, 1.379250, 
0173  1.358983, 1.400945, 1.402297, 1.427837, 1.464855, 1.443784, 1.427415, 1.422464, 1.070929, 1.060716, 1.064563, 1.065554, 1.069079, 1.075889, 1.068384, 1.070027, 
0174  1.075460, 1.082618, 1.083300, 1.083431, 1.093023, 1.108096, 1.107142, 1.112312, 1.154710, 1.605153, 1.225397, 1.314424, 1.308718, 1.332513, 1.323390, 1.354958, 
0175  1.398960, 1.384943, 1.362412, 1.378185, 1.065368, 1.064601, 1.065149, 1.061433, 1.059306, 1.070292, 1.067302, 1.061979, 1.069252, 1.075020, 1.081264, 1.076467, 
0176  1.083812, 1.093275, 1.089516, 1.103931, 1.133687, 1.593161, 1.202658, 1.259290, 1.277875, 1.299430, 1.285493, 1.319181, 1.363092, 1.348108, 1.323084, 1.344257,
0177  1.075912, 1.071438, 1.077006, 1.080259, 1.074469, 1.079985, 1.086186, 1.093308, 1.080566, 1.085392, 1.077386, 1.090488, 1.090022, 1.113544, 1.088709, 1.102750,
0178  1.145559, 1.572533, 1.173564, 1.239818, 1.257878, 1.268196, 1.270146, 1.304235, 1.343179, 1.325568, 1.294287, 1.323642, 1.069628, 1.064962, 1.068419, 1.087210, 
0179  1.078255, 1.078654, 1.072057, 1.069703, 1.080920, 1.089242, 1.100302, 1.092627, 1.094557, 1.108302, 1.103925, 1.101119, 1.132628, 1.459704, 1.162585, 1.228826, 
0180  1.223030, 1.254767, 1.250701, 1.279128, 1.320945, 1.316232, 1.277370, 1.312794,1.051953, 1.067999, 1.055347, 1.069504, 1.055176, 1.072807, 1.069857, 1.077925, 
0181  1.064145, 1.083152, 1.070389, 1.076320, 1.085833, 1.090424, 1.095051, 1.101244, 1.128931, 1.467267, 1.162999, 1.212314, 1.223013, 1.240110, 1.234485, 1.263299, 
0182  1.307635, 1.310394, 1.266581, 1.307123, 1.050234, 1.056716, 1.051719, 1.069001, 1.048164, 1.068260, 1.062359, 1.064068, 1.064559, 1.085571, 1.065183, 1.062737, 
0183  1.072117, 1.089194, 1.069072, 1.088440, 1.125319, 1.485413, 1.151196, 1.214000, 1.205797, 1.235945, 1.226787, 1.256677, 1.294848, 1.290351, 1.258192, 1.289197,
0184  1.119059, 1.118637, 1.118686, 1.119042, 1.120382, 1.128049, 1.122904, 1.123006, 1.125606, 1.131297, 1.125513, 1.127483, 1.130068, 1.138193, 1.133926, 1.137988, 
0185  1.156933, 1.421330, 1.190742, 1.215184, 1.223037, 1.230251, 1.232612, 1.256926, 1.291641, 1.280507, 1.266805, 1.319541,
0186     ])
0187 caloStage2Params.layer1HCalScaleETBins = cms.vint32([10, 15, 20, 25, 30, 35, 40, 45, 256])
0188 caloStage2Params.layer1HCalScaleFactors = cms.vdouble([
0189     0.938987, 1.305673, 1.477472, 0.765028, 1.815504, 1.798761, 0.778031, 1.236700, 1.095857, 1.378264, 1.433842, 1.167307, 1.352559, 1.690047, 0.891717, 1.854560, 1.793667, 1.553130, 1.707625, 1.916118, 1.481646, 1.541782, 1.490156, 1.479861, 1.468454, 1.495288, 1.610239, 1.588946,
0190     1.959768, 1.735311, 1.720083, 1.747394, 1.900315, 1.850038, 2.069131, 1.870382, 1.591709, 1.754415, 1.691840, 1.790095, 1.945210, 1.639001, 1.626256, 1.751506, 1.807836, 1.880716, 1.546263, 1.539627, 1.463240, 1.617592, 1.585718, 1.591975, 1.591791, 1.644489, 1.799823, 1.734512,
0191     1.797256, 1.578326, 1.692388, 1.503111, 1.672529, 1.519262, 1.582937, 1.886371, 1.778520, 1.612681, 1.398586, 1.603034, 1.454530, 1.328033, 1.471435, 1.636657, 1.527224, 1.512641, 1.410223, 1.388980, 1.421767, 1.572817, 1.578374, 1.588652, 1.621864, 1.729711, 1.918856, 1.803434,
0192     1.527994, 1.365275, 1.392104, 1.498491, 1.519979, 1.355727, 1.405808, 1.417970, 1.333893, 1.334817, 1.310798, 1.379249, 1.351143, 1.222402, 1.366848, 1.300051, 1.179617, 1.350694, 1.262175, 1.330276, 1.391770, 1.463739, 1.488979, 1.515451, 1.539866, 1.638972, 1.868056, 1.783022,
0193     1.258886, 1.373863, 1.306919, 1.356936, 1.370243, 1.385183, 1.251747, 1.269948, 1.232910, 1.298851, 1.233371, 1.269583, 1.283337, 1.207918, 1.230305, 1.229290, 1.279386, 1.274979, 1.274343, 1.273966, 1.336204, 1.425413, 1.427249, 1.446559, 1.494423, 1.585447, 1.732028, 1.837370,
0194     1.312747, 1.251113, 1.290323, 1.266130, 1.241814, 1.258689, 1.301229, 1.290575, 1.199071, 1.206718, 1.186692, 1.171748, 1.229947, 1.197701, 1.220179, 1.124275, 1.211136, 1.187607, 1.180801, 1.266214, 1.289450, 1.380122, 1.421375, 1.443506, 1.468912, 1.538188, 1.786963, 1.786963,
0195     1.261641, 1.242962, 1.289171, 1.176940, 1.187245, 1.230983, 1.142716, 1.183702, 1.166847, 1.166354, 1.134944, 1.191856, 1.156877, 1.177983, 1.108505, 1.149493, 1.141904, 1.158885, 1.151391, 1.183157, 1.277163, 1.366516, 1.373016, 1.392367, 1.425116, 1.527114, 1.527114, 1.527114,
0196     1.228298, 1.191839, 1.256205, 1.201946, 1.129291, 1.186304, 1.133396, 1.101866, 1.210155, 1.100159, 1.108861, 1.093585, 1.126787, 1.106975, 1.153106, 1.112337, 1.142352, 1.155770, 1.124314, 1.167838, 1.241431, 1.299161, 1.329571, 1.335993, 1.378090, 1.378090, 1.378090, 1.378090,
0197     1.115713, 1.101212, 1.093974, 1.094815, 1.090964, 1.072051, 1.062362, 1.035872, 1.031726, 1.023072, 1.010703, 1.015952, 1.020603, 1.024692, 1.025791, 1.156998, 1.096629, 1.110618, 1.091158, 1.122913, 1.195124, 1.271093, 1.277982, 1.321275, 1.321275, 1.321275, 1.321275, 1.321275,
0198     ])
0199 # HF 1x1 scale factors will be a 5*12 array:
0200 #  12 eta scale factors (30-41)
0201 #  in 5 REAL ET bins (5, 20, 30, 50, Max)
0202 #  So, index = etBin*12+ietaHF
0203 caloStage2Params.layer1HFScaleETBins = cms.vint32([5, 20, 30, 50, 256])
0204 caloStage2Params.layer1HFScaleFactors = cms.vdouble([
0205     1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 
0206     1.767080, 1.767080, 1.755186, 1.769951, 1.763527, 1.791043, 1.898787, 1.982235, 2.071074, 2.193011, 2.356886, 2.403384, 
0207     2.170477, 2.170477, 2.123540, 2.019866, 1.907698, 1.963179, 1.989122, 2.035251, 2.184642, 2.436399, 2.810884, 2.923750, 
0208     1.943941, 1.943941, 1.899826, 1.813950, 1.714978, 1.736184, 1.785928, 1.834211, 1.944230, 2.153565, 2.720887, 2.749795, 
0209     1.679984, 1.679984, 1.669753, 1.601871, 1.547276, 1.577805, 1.611497, 1.670184, 1.775022, 1.937061, 2.488311, 2.618629, 
0210     ])