Type1MET/python/pfMETsysShiftCorrections_cfi.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 #--------------------------------------------------------------------------------
0004 # parametrization of MET x/y shift vs. sumEt
0005 pfMEtSysShiftCorrParameters_2011runAvsSumEt_data = cms.PSet(
0006     px = cms.string("-3.365e-1 + 4.801e-3*sumEt"),
0007     py = cms.string("+2.578e-1 - 6.124e-3*sumEt")
0008 )
0009
0010 pfMEtSysShiftCorrParameters_2011runAvsSumEt_mc = cms.PSet(
0011     px = cms.string("-9.389e-2 + 1.815e-4*sumEt"),
0012     py = cms.string("+1.571e-1 - 3.710e-3*sumEt")
0013 )
0014
0015 pfMEtSysShiftCorrParameters_2011runBvsSumEt_data = cms.PSet(
0016     px = cms.string("-3.265e-1 + 5.162e-3*sumEt"),
0017     py = cms.string("-1.956e-2 - 6.299e-3*sumEt")
0018 )
0019
0020 pfMEtSysShiftCorrParameters_2011runBvsSumEt_mc = cms.PSet(
0021     px = cms.string("-1.070e-1 + 9.587e-5*sumEt"),
0022     py = cms.string("-1.517e-2 - 3.357e-3*sumEt")
0023 )
0024
0025 pfMEtSysShiftCorrParameters_2011runAplusBvsSumEt_data = cms.PSet(
0026     px = cms.string("-5.65217e-01 + 5.42436e-03*sumEt"),
0027     py = cms.string("+4.54054e-01 - 6.73607e-03*sumEt")
0028 )
0029
0030 pfMEtSysShiftCorrParameters_2011runAplusBvsSumEt_mc = cms.PSet(
0031     px = cms.string("-4.53909e-02 - 2.55863e-05*sumEt"),
0032     py = cms.string("+1.27947e-01 - 3.62604e-03*sumEt")
0033 )
0034
0035 pfMEtSysShiftCorrParameters_2012runAvsSumEt_data = cms.PSet(
0036     px = cms.string("-7.67892e-01 + 5.76983e-03*sumEt"),
0037     py = cms.string("+5.54005e-01 - 2.94046e-03*sumEt")
0038 )
0039
0040 pfMEtSysShiftCorrParameters_2012runAvsSumEt_mc = cms.PSet(
0041     px = cms.string("+1.77344e-01 - 1.34333e-03*sumEt"),
0042     py = cms.string("+8.08402e-01 - 2.84264e-03*sumEt")
0043 )
0044 #--------------------------------------------------------------------------------
0045
0046 #--------------------------------------------------------------------------------
0047 # parametrization of MET x/y shift vs. Nvtx
0048 pfMEtSysShiftCorrParameters_2011runAvsNvtx_data = cms.PSet(
0049     px = cms.string("+3.87339e-1 + 2.58294e-1*Nvtx"),
0050     py = cms.string("-7.83502e-1 - 2.88899e-1*Nvtx")
0051 )
0052
0053 pfMEtSysShiftCorrParameters_2011runAvsNvtx_mc = cms.PSet(
0054     px = cms.string("-1.94451e-2 - 4.38986e-3*Nvtx"),
0055     py = cms.string("-4.31368e-1 - 1.90753e-1*Nvtx")
0056 )
0057
0058 pfMEtSysShiftCorrParameters_2011runBvsNvtx_data = cms.PSet(
0059     px = cms.string("+6.64470e-1 + 2.71292e-1*Nvtx"),
0060     py = cms.string("-1.23999e0 - 3.18661e-1*Nvtx")
0061 )
0062
0063 pfMEtSysShiftCorrParameters_2011runBvsNvtx_mc = cms.PSet(
0064     px = cms.string("-9.89706e-2 + 6.64796e-3*Nvtx"),
0065     py = cms.string("-5.32495e-1 - 1.82195e-1*Nvtx")
0066 )
0067
0068 pfMEtSysShiftCorrParameters_2011runAplusBvsNvtx_data = cms.PSet(
0069     px = cms.string("+3.64118e-01 + 2.93853e-01*Nvtx"),
0070     py = cms.string("-7.17757e-01 - 3.57309e-01*Nvtx")
0071 )
0072
0073 pfMEtSysShiftCorrParameters_2011runAplusBvsNvtx_mc = cms.PSet(
0074     px = cms.string("-4.79178e-02 + 8.62653e-04*Nvtx"),
0075     py = cms.string("-4.54408e-01 - 1.89684e-01*Nvtx")
0076 )
0077
0078 pfMEtSysShiftCorrParameters_2012runAvsNvtx_data = cms.PSet(
0079     px = cms.string("+3.54233e-01 + 2.65299e-01*Nvtx"),
0080     py = cms.string("+1.88923e-01 - 1.66425e-01*Nvtx")
0081 )
0082
0083 pfMEtSysShiftCorrParameters_2012runAvsNvtx_mc = cms.PSet(
0084     px = cms.string("-2.99576e-02 - 6.61932e-02*Nvtx"),
0085     py = cms.string("+3.70819e-01 - 1.48617e-01*Nvtx")
0086 )
0087
0088 pfMEtSysShiftCorrParameters_2012runAplusBvsNvtx_data = cms.PSet(
0089     px = cms.string("+1.68804e-01 + 3.37139e-01*Nvtx"),
0090     py = cms.string("-1.72555e-01 - 1.79594e-01*Nvtx")
0091 )
0092
0093 pfMEtSysShiftCorrParameters_2012runAplusBvsNvtx_mc = cms.PSet(
0094     px = cms.string("+2.22335e-02 - 6.59183e-02*Nvtx"),
0095     py = cms.string("+1.52720e-01 - 1.28052e-01*Nvtx")
0096 )
0097
0098 pfMEtSysShiftCorrParameters_2012runABCDvsNvtx_data = cms.PSet( # CV: ReReco data + Summer'13 JEC
0099     px = cms.string("+4.83642e-02 + 2.48870e-01*Nvtx"),
0100     py = cms.string("-1.50135e-01 - 8.27917e-02*Nvtx")
0101 )
0102
0103 pfMEtSysShiftCorrParameters_2012runABCDvsNvtx_mc = cms.PSet( # CV: Summer'12 MC + Summer'13 JEC
0104     px = cms.string("+1.62861e-01 - 2.38517e-02*Nvtx"),
0105     py = cms.string("+3.60860e-01 - 1.30335e-01*Nvtx")
0106 )
0107
0108 selectedVerticesForMEtCorr = cms.EDFilter("VertexSelector",
0109     src = cms.InputTag('offlinePrimaryVertices'),
0110     cut = cms.string("isValid & ndof >= 4 & chi2 > 0 & tracksSize > 0 & abs(z) < 24 & abs(position.Rho) < 2."),
0111     filter = cms.bool(False)
0112 )
0113 #--------------------------------------------------------------------------------
0114
0115 pfMEtSysShiftCorr = cms.EDProducer("SysShiftMETcorrInputProducer",
0116     src = cms.InputTag('pfMet'), # "raw"/uncorrected PFMEt, needed to access sumEt
0117     srcVertices = cms.InputTag('selectedVerticesForMEtCorr'),
0118     parameter = pfMEtSysShiftCorrParameters_2012runABCDvsNvtx_data
0119 )
0120
0121 pfMEtSysShiftCorrSequence = cms.Sequence(selectedVerticesForMEtCorr * pfMEtSysShiftCorr)