L1TCalorimeter/python/caloParams_2022_v0_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.egEtaCut                   = cms.int32(28)
0025 caloStage2Params.egSeedThreshold            = cms.double(2.)
0026 caloStage2Params.egNeighbourThreshold       = cms.double(1.)
0027 caloStage2Params.egHcalThreshold            = cms.double(0.)
0028 caloStage2Params.egTrimmingLUTFile          = cms.FileInPath("L1Trigger/L1TCalorimeter/data/egTrimmingLUT_10_v16.01.19.txt")
0029 caloStage2Params.egMaxHcalEt                = cms.double(0.)
0030 caloStage2Params.egMaxPtHOverE          = cms.double(128.)
0031 caloStage2Params.egMaxHOverELUTFile         = cms.FileInPath("L1Trigger/L1TCalorimeter/data/HoverEIdentification_0.995_v15.12.23.txt")
0032 caloStage2Params.egHOverEcutBarrel          = cms.int32(3)
0033 caloStage2Params.egHOverEcutEndcap          = cms.int32(4)
0034 caloStage2Params.egBypassExtHOverE          = cms.uint32(0)
0035 caloStage2Params.egCompressShapesLUTFile    = cms.FileInPath("L1Trigger/L1TCalorimeter/data/egCompressLUT_v4.txt")
0036 caloStage2Params.egShapeIdType              = cms.string("compressed")
0037 caloStage2Params.egShapeIdVersion           = cms.uint32(0)
0038 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
0039
0040 caloStage2Params.egPUSType                  = cms.string("None")
0041 caloStage2Params.egIsolationType            = cms.string("compressed")
0042 caloStage2Params.egIsoLUTFile               = cms.FileInPath("L1Trigger/L1TCalorimeter/data/EG_Iso_LUT_04_04_2017.2.txt")
0043 caloStage2Params.egIsoLUTFile2               = cms.FileInPath("L1Trigger/L1TCalorimeter/data/EG_LoosestIso_2018.2.txt")
0044 caloStage2Params.egIsoAreaNrTowersEta       = cms.uint32(2)
0045 caloStage2Params.egIsoAreaNrTowersPhi       = cms.uint32(4)
0046 caloStage2Params.egIsoVetoNrTowersPhi       = cms.uint32(2)
0047 #caloStage2Params.egIsoPUEstTowerGranularity = cms.uint32(1)
0048 #caloStage2Params.egIsoMaxEtaAbsForTowerSum  = cms.uint32(4)
0049 #caloStage2Params.egIsoMaxEtaAbsForIsoSum    = cms.uint32(27)
0050 caloStage2Params.egPUSParams                = cms.vdouble(1,4,32) #Isolation window in firmware goes up to abs(ieta)=32 for now
0051 caloStage2Params.egCalibrationType          = cms.string("compressed")
0052 caloStage2Params.egCalibrationVersion       = cms.uint32(0)
0053 #caloStage2Params.egCalibrationLUTFile       = cms.FileInPath("L1Trigger/L1TCalorimeter/data/EG_Calibration_LUT_FW_v17.04.04_shapeIdentification_adapt0.99_compressedieta_compressedE_compressedshape_v15.12.08_correct.txt")
0054 caloStage2Params.egCalibrationLUTFile       = cms.FileInPath("L1Trigger/L1TCalorimeter/data/corrections_Trimming10_compressedieta_compressedE_compressedshape_PANTELIS_v2_NEW_CALIBRATIONS_withShape_v17.04.04.txt")
0055
0056 # Tau
0057 caloStage2Params.tauLsb                        = cms.double(0.5)
0058 caloStage2Params.isoTauEtaMax                  = cms.int32(25)
0059 caloStage2Params.tauSeedThreshold              = cms.double(0.)
0060 caloStage2Params.tauNeighbourThreshold         = cms.double(0.)
0061 caloStage2Params.tauIsoAreaNrTowersEta         = cms.uint32(2)
0062 caloStage2Params.tauIsoAreaNrTowersPhi         = cms.uint32(4)
0063 caloStage2Params.tauIsoVetoNrTowersPhi         = cms.uint32(2)
0064 caloStage2Params.tauPUSType                    = cms.string("None")
0065 caloStage2Params.tauIsoLUTFile                 = cms.FileInPath("L1Trigger/L1TCalorimeter/data/Tau_Iso_LUT_Option_31_extrap_2018_FW_v10.0.0.txt")
0066 caloStage2Params.tauIsoLUTFile2                = cms.FileInPath("L1Trigger/L1TCalorimeter/data/Tau_Iso_LUT_Option_31_extrap_2018_FW_v10.0.0.txt")
0067 #caloStage2Params.tauCalibrationLUTFile         = cms.FileInPath("L1Trigger/L1TCalorimeter/data/Tau_Calibration_LUT_2017_Layer1Calibration_FW_v12.0.0.txt")
0068 caloStage2Params.tauCalibrationLUTFile         = cms.FileInPath("L1Trigger/L1TCalorimeter/data/Tau_Calibration_LUT_2018_Layer1CalibrationNewHCAL_FW_v13.0.0.txt")
0069 caloStage2Params.tauCompressLUTFile            = cms.FileInPath("L1Trigger/L1TCalorimeter/data/tauCompressAllLUT_12bit_v3.txt")
0070 caloStage2Params.tauPUSParams                  = cms.vdouble(1,4,32)
0071
0072 # jets
0073 caloStage2Params.jetLsb                = cms.double(0.5)
0074 caloStage2Params.jetSeedThreshold      = cms.double(4.0)
0075 caloStage2Params.jetNeighbourThreshold = cms.double(0.)
0076 caloStage2Params.jetPUSType            = cms.string("ChunkyDonut")
0077 caloStage2Params.jetBypassPUS          = cms.uint32(0)
0078
0079 # Calibration options
0080 caloStage2Params.jetCalibrationType = cms.string("LUT")
0081
0082 caloStage2Params.jetCompressPtLUTFile     = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_pt_compress_2017v1.txt")
0083 caloStage2Params.jetCompressEtaLUTFile    = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_eta_compress_2017v1.txt")
0084 caloStage2Params.jetCalibrationLUTFile    = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_calib_2018v1_ECALZS_noHFJEC.txt")
0085
0086
0087 # sums: 0=ET, 1=HT, 2=MET, 3=MHT
0088 caloStage2Params.etSumLsb                = cms.double(0.5)
0089 caloStage2Params.etSumEtaMin             = cms.vint32(1, 1, 1, 1, 1)
0090 caloStage2Params.etSumEtaMax             = cms.vint32(28,  26, 28,  26, 28)
0091 caloStage2Params.etSumEtThreshold        = cms.vdouble(0.,  30.,  0.,  30., 0.) # only 2nd (HT) and 4th (MHT) values applied
0092 caloStage2Params.etSumMetPUSType         = cms.string("LUT") # et threshold from this LUT supercedes et threshold in line above
0093 caloStage2Params.etSumEttPUSType         = cms.string("None")
0094 caloStage2Params.etSumEcalSumPUSType     = cms.string("None")
0095 caloStage2Params.etSumBypassMetPUS       = cms.uint32(0)
0096 caloStage2Params.etSumBypassEttPUS       = cms.uint32(1)
0097 caloStage2Params.etSumBypassEcalSumPUS    = cms.uint32(1)
0098 caloStage2Params.etSumXCalibrationType    = cms.string("None")
0099 caloStage2Params.etSumYCalibrationType    = cms.string("None")
0100 caloStage2Params.etSumEttCalibrationType  = cms.string("None")
0101 caloStage2Params.etSumEcalSumCalibrationType = cms.string("None")
0102
0103 caloStage2Params.etSumMetPUSLUTFile               = cms.FileInPath("L1Trigger/L1TCalorimeter/data/newSFHBHEOnp5_METPUM_211124.txt")
0104 caloStage2Params.etSumEttPUSLUTFile               = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_towEtThresh_dummy.txt")
0105 caloStage2Params.etSumEcalSumPUSLUTFile           = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_towEtThresh_dummy.txt")
0106 caloStage2Params.etSumXCalibrationLUTFile         = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_etSumPUS_dummy.txt")
0107 caloStage2Params.etSumYCalibrationLUTFile         = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_etSumPUS_dummy.txt")
0108 caloStage2Params.etSumEttCalibrationLUTFile       = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_etSumPUS_dummy.txt")
0109 caloStage2Params.etSumEcalSumCalibrationLUTFile   = cms.FileInPath("L1Trigger/L1TCalorimeter/data/lut_etSumPUS_dummy.txt")
0110
0111
0112 # Layer 1 SF
0113 caloStage2Params.layer1ECalScaleETBins = cms.vint32([3, 6, 9, 12, 15, 20, 25, 30, 35, 40, 45, 55, 70, 256])
0114 caloStage2Params.layer1ECalScaleFactors = cms.vdouble([
0115 1.12, 1.13, 1.13, 1.12, 1.12, 1.12, 1.13, 1.13, 1.13, 1.13, 1.13, 1.13, 1.13, 1.14, 1.14, 1.13, 1.13, 1.30, 1.16, 1.28, 1.27, 1.30, 1.30, 1.33, 1.36, 0.00, 0.00, 0.00,
0116
0117 1.12, 1.13, 1.13, 1.12, 1.12, 1.12, 1.13, 1.13, 1.13, 1.13, 1.13, 1.13, 1.13, 1.14, 1.14, 1.13, 1.13, 1.30, 1.16, 1.28, 1.27, 1.30, 1.30, 1.33, 1.36, 1.39, 0.00, 0.00,
0118
0119 1.08, 1.08, 1.09, 1.08, 1.09, 1.08, 1.08, 1.08, 1.09, 1.09, 1.15, 1.09, 1.10, 1.10, 1.10, 1.09, 1.11, 1.24, 1.11, 1.21, 1.21, 1.23, 1.24, 1.26, 1.31, 1.35, 1.22, 0.00,
0120
0121 1.06, 1.06, 1.06, 1.07, 1.07, 1.07, 1.07, 1.07, 1.06, 1.07, 1.07, 1.07, 1.07, 1.08, 1.07, 1.08, 1.08, 1.21, 1.09, 1.17, 1.18, 1.20, 1.22, 1.24, 1.29, 1.29, 1.16, 1.09,
0122
0123 1.04, 1.05, 1.05, 1.05, 1.05, 1.05, 1.04, 1.05, 1.05, 1.05, 1.06, 1.05, 1.06, 1.06, 1.06, 1.07, 1.05, 1.17, 1.09, 1.15, 1.16, 1.18, 1.19, 1.22, 1.26, 1.26, 1.15, 1.08,
0124
0125 1.03, 1.03, 1.04, 1.04, 1.04, 1.04, 1.05, 1.04, 1.04, 1.04, 1.04, 1.05, 1.05, 1.05, 1.06, 1.05, 1.06, 1.14, 1.07, 1.13, 1.13, 1.17, 1.19, 1.21, 1.24, 1.29, 1.14, 1.07,
0126
0127 1.03, 1.03, 1.03, 1.03, 1.03, 1.03, 1.04, 1.03, 1.04, 1.03, 1.04, 1.04, 1.04, 1.04, 1.04, 1.03, 1.05, 1.13, 1.06, 1.12, 1.13, 1.16, 1.18, 1.20, 1.23, 1.26, 1.09, 1.06,
0128
0129 1.03, 1.03, 1.03, 1.03, 1.03, 1.03, 1.03, 1.03, 1.03, 1.03, 1.02, 1.03, 1.04, 1.04, 1.04, 1.04, 1.04, 1.11, 1.06, 1.10, 1.13, 1.14, 1.16, 1.19, 1.22, 1.23, 1.07, 1.05,
0130
0131 1.02, 1.02, 1.02, 1.02, 1.03, 1.03, 1.03, 1.03, 1.03, 1.03, 1.03, 1.03, 1.03, 1.03, 1.03, 1.03, 1.02, 1.11, 1.05, 1.10, 1.12, 1.14, 1.15, 1.18, 1.21, 1.25, 1.07, 1.05,
0132
0133 1.02, 1.02, 1.02, 1.02, 1.02, 1.02, 1.02, 1.02, 1.03, 1.02, 1.03, 1.03, 1.03, 1.03, 1.03, 1.03, 1.03, 1.10, 1.05, 1.09, 1.11, 1.13, 1.15, 1.17, 1.20, 1.23, 1.06, 1.05,
0134
0135 1.02, 1.02, 1.02, 1.02, 1.02, 1.02, 1.02, 1.02, 1.02, 1.02, 1.02, 1.02, 1.03, 1.03, 1.03, 1.03, 1.03, 1.10, 1.05, 1.10, 1.10, 1.13, 1.14, 1.18, 1.20, 1.23, 1.09, 1.05,
0136
0137 1.01, 1.01, 1.01, 1.01, 1.01, 1.02, 1.02, 1.02, 1.02, 1.01, 1.02, 1.02, 1.03, 1.03, 1.03, 1.03, 1.03, 1.08, 1.05, 1.09, 1.10, 1.12, 1.14, 1.17, 1.19, 1.21, 1.05, 1.05,
0138
0139 1.01, 1.00, 1.00, 1.01, 1.01, 1.02, 1.01, 1.01, 1.01, 1.01, 1.01, 1.02, 1.02, 1.02, 1.02, 1.02, 1.02, 1.07, 1.05, 1.09, 1.09, 1.11, 1.13, 1.15, 1.19, 1.14, 1.00, 1.01,
0140
0141 1.00, 1.00, 1.00, 1.00, 1.01, 1.02, 1.01, 1.01, 1.01, 1.01, 1.01, 1.01, 1.02, 1.02, 1.02, 1.02, 1.02, 1.06, 1.04, 1.08, 1.09, 1.11, 1.13, 1.13, 1.18, 1.21, 1.00, 1.01
0142 ])
0143
0144 caloStage2Params.layer1HCalScaleETBins = cms.vint32([6, 9, 12, 15, 20, 25, 30, 35, 40, 45, 55, 70, 256])
0145 caloStage2Params.layer1HCalScaleFactors = cms.vdouble([
0146 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
0147
0148 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
0149
0150 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
0151
0152 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
0153
0154 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
0155
0156 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
0157
0158 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
0159
0160 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
0161
0162 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
0163
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, 1.00, 1.00, 1.00, 1.00, 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
0166 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
0167
0168 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
0169
0170 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00
0171 ])
0172
0173 caloStage2Params.layer1HFScaleETBins = cms.vint32([6, 9, 12, 15, 20, 25, 30, 35, 40, 45, 55, 70, 256])
0174 caloStage2Params.layer1HFScaleFactors = cms.vdouble([
0175 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
0176
0177 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
0178
0179 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
0180
0181 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
0182
0183 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
0184
0185 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
0186
0187 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
0188
0189 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
0190
0191 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
0192
0193 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
0194
0195 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
0196
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
0199 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00
0200 ])