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_EBEE_v16.04.04.txt")
0038 caloStage2Params.egIsoAreaNrTowersEta       = cms.uint32(2)
0039 caloStage2Params.egIsoAreaNrTowersPhi       = cms.uint32(4)
0040 caloStage2Params.egIsoVetoNrTowersPhi       = cms.uint32(3)
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_Trimming5_compressedieta_compressedE_compressedshape_v16.04.01.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_NoLayer1Calibration_noCompressionBlock_v4.0.0.txt")
0058 caloStage2Params.tauIsoLUTFile2                = cms.FileInPath("L1Trigger/L1TCalorimeter/data/Tau_Iso_LUT_Option_21_NoLayer1Calibration_noCompressionBlock_v4.0.0.txt")
0059 caloStage2Params.tauCalibrationLUTFile         = cms.FileInPath("L1Trigger/L1TCalorimeter/data/Tau_Calibration_LUT_NoLayer1Calibration_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(1.5)
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("function6PtParams22EtaBins")
0073 
0074 
0075 #Vector with 6 parameters for eta bin, from low eta to high
0076 # 1,0,1,0,1,1 gives no correction
0077 # must be in this form as may require > 255 arguments
0078 jetCalibParamsVector = cms.vdouble()
0079 jetCalibParamsVector.extend([
0080     1.52819258712,0,1,0,1,1, # Constant calibrations in HF bins
0081     1.32880934866,0,1,0,1,1,
0082     1.35650073562,0,1,0,1,1,
0083     1.55766262213,0,1,0,1,1,
0084     12.42641269,6.35085569,1.21845374,-14.70772743,0.00118945,-10.56931753,
0085     12.14768583,8.05320174,0.92028575,-15.48969055,0.00155426,-10.48022143,
0086     11.87301715,10.45078171,1.02636603,-16.87586292,0.00210878,-10.34892982,
0087     12.09471152,9.36873133,0.96484984,-16.09937338,0.00177759,-10.43545594,
0088     4.69610083,11.50217825,1.31071693,-23.88756705,0.00336036,-19.28130829,
0089     7.08667633,23.68511274,2.52417732,-43.09154102,0.00359229,-18.80055010,
0090     3.44975454,28.42124526,2.58428891,-137.42027137,0.00669104,-19.53448313,
0091     3.44975454,28.42124526,2.58428891,-137.42027137,0.00669104,-19.53448313,
0092     7.08667633,23.68511274,2.52417732,-43.09154102,0.00359229,-18.80055010,
0093     4.69610083,11.50217825,1.31071693,-23.88756705,0.00336036,-19.28130829,
0094     12.09471152,9.36873133,0.96484984,-16.09937338,0.00177759,-10.43545594,
0095     11.87301715,10.45078171,1.02636603,-16.87586292,0.00210878,-10.34892982,
0096     12.14768583,8.05320174,0.92028575,-15.48969055,0.00155426,-10.48022143,
0097     12.42641269,6.35085569,1.21845374,-14.70772743,0.00118945,-10.56931753,
0098     1.55766262213,0,1,0,1,1, # Constant calibrations in HF bins
0099     1.35650073562,0,1,0,1,1,
0100     1.32880934866,0,1,0,1,1,
0101     1.52819258712,0,1,0,1,1
0102 ])
0103 caloStage2Params.jetCalibrationParams  = jetCalibParamsVector 
0104 
0105 # sums: 0=ET, 1=HT, 2=MET, 3=MHT
0106 caloStage2Params.etSumLsb                = cms.double(0.5)
0107 caloStage2Params.etSumEtaMin             = cms.vint32(1, 1, 1, 1)
0108 caloStage2Params.etSumEtaMax             = cms.vint32(28,  28,  28,  28)
0109 caloStage2Params.etSumEtThreshold        = cms.vdouble(0.,  35.,  0.,  35.)
0110 
0111 # Layer 1 LUT specification
0112 #
0113 # Et-dependent scale factors
0114 # ECal/HCal scale factors will be a 9*28 array:
0115 #   28 eta scale factors (1-28)
0116 #   in 9 ET bins (10, 15, 20, 25, 30, 35, 40, 45, Max)
0117 #  So, index = etBin*28+ieta
0118 caloStage2Params.layer1ECalScaleETBins = cms.vint32([1])
0119 caloStage2Params.layer1ECalScaleFactors = cms.vdouble([1.]*28)
0120 caloStage2Params.layer1HCalScaleETBins = cms.vint32([1])
0121 caloStage2Params.layer1HCalScaleFactors = cms.vdouble([1.]*28)
0122 caloStage2Params.layer1HFScaleETBins = cms.vint32([1])
0123 caloStage2Params.layer1HFScaleFactors = cms.vdouble([1.]*12)