L1TCalorimeter/python/caloStage2Params_2015_v1_0_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_5.txt")
0028 caloStage2Params.egMaxHcalEt                = cms.double(0.)
0029 caloStage2Params.egMaxPtHOverE          = cms.double(128.)
0030 caloStage2Params.egMaxHOverELUTFile         = cms.FileInPath("L1Trigger/L1TCalorimeter/data/egMaxHOverELUT_995eff.txt")
0031 caloStage2Params.egCompressShapesLUTFile    = cms.FileInPath("L1Trigger/L1TCalorimeter/data/egCompressLUT_v1.txt")
0032 caloStage2Params.egShapeIdType              = cms.string("compressed")
0033 caloStage2Params.egShapeIdVersion           = cms.uint32(0)
0034 caloStage2Params.egShapeIdLUTFile           = cms.FileInPath("L1Trigger/L1TCalorimeter/data/shapeIdentification_0.99_compressedieta_compressedE_compressedshape_v15.10.20.txt")
0035 caloStage2Params.egPUSType                  = cms.string("None")
0036 caloStage2Params.egIsoLUTFile               = cms.FileInPath("L1Trigger/L1TCalorimeter/data/egIsoLUTPU40bx25NrRings4Eff95.txt")
0037 caloStage2Params.egIsoAreaNrTowersEta       = cms.uint32(2)
0038 caloStage2Params.egIsoAreaNrTowersPhi       = cms.uint32(4)
0039 caloStage2Params.egIsoVetoNrTowersPhi       = cms.uint32(3)
0040 #caloStage2Params.egIsoPUEstTowerGranularity = cms.uint32(1)
0041 #caloStage2Params.egIsoMaxEtaAbsForTowerSum  = cms.uint32(4)
0042 #caloStage2Params.egIsoMaxEtaAbsForIsoSum    = cms.uint32(27)
0043 caloStage2Params.egPUSParams                = cms.vdouble(1,4,27)
0044 caloStage2Params.egCalibrationType          = cms.string("compressed")
0045 caloStage2Params.egCalibrationVersion       = cms.uint32(0)
0046 caloStage2Params.egCalibrationLUTFile       = cms.FileInPath("L1Trigger/L1TCalorimeter/data/corrections_Trimming5_PU40bx25_compressedieta_compressedE_compressedshape_v15.10.20.txt")
0047
0048 # Tau
0049 caloStage2Params.tauLsb                        = cms.double(0.5)
0050 caloStage2Params.tauSeedThreshold              = cms.double(0.)
0051 caloStage2Params.tauNeighbourThreshold         = cms.double(0.)
0052 caloStage2Params.tauIsoAreaNrTowersEta         = cms.uint32(2)
0053 caloStage2Params.tauIsoAreaNrTowersPhi         = cms.uint32(4)
0054 caloStage2Params.tauIsoVetoNrTowersPhi         = cms.uint32(2)
0055 caloStage2Params.tauPUSType                 = cms.string("None")
0056 caloStage2Params.tauIsoLUTFile                 = cms.FileInPath("L1Trigger/L1TCalorimeter/data/tauIsoLUTetPU.txt")
0057 caloStage2Params.tauCalibrationLUTFile         = cms.FileInPath("L1Trigger/L1TCalorimeter/data/tauCalibrationLUT.txt")
0058 caloStage2Params.tauPUSParams                  = cms.vdouble(1,4,27)
0059
0060 # jets
0061 caloStage2Params.jetLsb                = cms.double(0.5)
0062 caloStage2Params.jetSeedThreshold      = cms.double(1.5)
0063 caloStage2Params.jetNeighbourThreshold = cms.double(0.)
0064 caloStage2Params.jetPUSType            = cms.string("ChunkyDonut")
0065
0066 #Calibration options
0067 # e.g. function6PtParams22EtaBins
0068 caloStage2Params.jetCalibrationType    = cms.string("None")
0069 #caloStage2Params.jetCalibrationType = cms.string("function6PtParams22EtaBins")
0070
0071
0072 #Vector with 6 parameters for eta bin, from low eta to high
0073 # 1,0,1,0,1,1 gives no correction
0074 # must be in this form as may require > 255 arguments
0075 jetCalibParamsVector = cms.vdouble()
0076 jetCalibParamsVector.extend([
0077     1.52819258712,0,1,0,1,1, # Constant calibrations in HF bins
0078     1.32880934866,0,1,0,1,1,
0079     1.35650073562,0,1,0,1,1,
0080     1.55766262213,0,1,0,1,1,
0081     12.42641269,6.35085569,1.21845374,-14.70772743,0.00118945,-10.56931753,
0082     12.14768583,8.05320174,0.92028575,-15.48969055,0.00155426,-10.48022143,
0083     11.87301715,10.45078171,1.02636603,-16.87586292,0.00210878,-10.34892982,
0084     12.09471152,9.36873133,0.96484984,-16.09937338,0.00177759,-10.43545594,
0085     4.69610083,11.50217825,1.31071693,-23.88756705,0.00336036,-19.28130829,
0086     7.08667633,23.68511274,2.52417732,-43.09154102,0.00359229,-18.80055010,
0087     3.44975454,28.42124526,2.58428891,-137.42027137,0.00669104,-19.53448313,
0088     3.44975454,28.42124526,2.58428891,-137.42027137,0.00669104,-19.53448313,
0089     7.08667633,23.68511274,2.52417732,-43.09154102,0.00359229,-18.80055010,
0090     4.69610083,11.50217825,1.31071693,-23.88756705,0.00336036,-19.28130829,
0091     12.09471152,9.36873133,0.96484984,-16.09937338,0.00177759,-10.43545594,
0092     11.87301715,10.45078171,1.02636603,-16.87586292,0.00210878,-10.34892982,
0093     12.14768583,8.05320174,0.92028575,-15.48969055,0.00155426,-10.48022143,
0094     12.42641269,6.35085569,1.21845374,-14.70772743,0.00118945,-10.56931753,
0095     1.55766262213,0,1,0,1,1, # Constant calibrations in HF bins
0096     1.35650073562,0,1,0,1,1,
0097     1.32880934866,0,1,0,1,1,
0098     1.52819258712,0,1,0,1,1
0099 ])
0100 caloStage2Params.jetCalibrationParams  = jetCalibParamsVector
0101
0102 # sums: 0=ET, 1=HT, 2=MET, 3=MHT
0103 caloStage2Params.etSumLsb                = cms.double(0.5)
0104 caloStage2Params.etSumEtaMin             = cms.vint32(0, 0, 0, 0)
0105 caloStage2Params.etSumEtaMax             = cms.vint32(40,  36,  40,  36)
0106 caloStage2Params.etSumEtThreshold        = cms.vdouble(0.,  0.,   0.,   0.)
0107
0108 # Layer 1 LUT specification
0109 #
0110 # Et-dependent scale factors
0111 # ECal/HCal scale factors will be a 9*28 array:
0112 #   28 eta scale factors (1-28)
0113 #   in 9 ET bins (10, 15, 20, 25, 30, 35, 40, 45, Max)
0114 #  So, index = etBin*28+ieta
0115 caloStage2Params.layer1ECalScaleETBins = cms.vint32([1])
0116 caloStage2Params.layer1ECalScaleFactors = cms.vdouble([1.]*28)
0117 caloStage2Params.layer1HCalScaleETBins = cms.vint32([1])
0118 caloStage2Params.layer1HCalScaleFactors = cms.vdouble([1.]*28)
0119 caloStage2Params.layer1HFScaleETBins = cms.vint32([1])
0120 caloStage2Params.layer1HFScaleFactors = cms.vdouble([1.]*12)
0121
0122