METPUSubtraction/python/pfNoPUMET_cff.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 pfNoPUMEtTask = cms.Task()
0004
0005 from JetMETCorrections.Configuration.JetCorrectionServices_cff import *
0006 calibratedAK4PFJetsForPFNoPUMEt = cms.EDProducer('PFJetCorrectionProducer',
0007     src = cms.InputTag('ak4PFJets'),
0008     correctors = cms.vstring('ak4PFL1FastL2L3') # NOTE: use "ak4PFL1FastL2L3" for MC / "ak4PFL1FastL2L3Residual" for Data
0009 )
0010 ak4PFJetTaskForPFNoPUMEt = cms.Task(calibratedAK4PFJetsForPFNoPUMEt)
0011 pfNoPUMEtTask.add(ak4PFJetTaskForPFNoPUMEt)
0012
0013 from RecoJets.JetProducers.PileupJetID_cfi import *
0014 puJetIdForPFNoPUMEt = pileupJetId.clone(
0015     algos = cms.VPSet(
0016         full_53x,
0017         cutbased,
0018         PhilV1
0019         ),
0020 #    label = "fullId", #MM does not work for weird reasons, cannot be cloned properly
0021     produceJetIds    = True,
0022     runMvas          = True,
0023     jets             = "calibratedAK4PFJetsForPFNoPUMEt",
0024     applyJec         = False,
0025     inputIsCorrected = True,
0026     )
0027 pfNoPUMEtTask.add(puJetIdForPFNoPUMEt)
0028
0029 from JetMETCorrections.Type1MET.pfMETCorrectionType0_cfi import *
0030 pfNoPUMEtTask.add(type0PFMEtCorrection)
0031 pfCandidateToVertexAssociationForPFNoPUMEt = pfCandidateToVertexAssociation.clone(
0032     MaxNumberOfAssociations = 1,
0033     doReassociation         = False,
0034     FinalAssociation        = 1,
0035     nTrackWeight            = 0.
0036 )
0037 pfNoPUMEtTask.add(pfCandidateToVertexAssociationForPFNoPUMEt)
0038 pfMETcorrType0ForPFNoPUMEt = pfMETcorrType0.clone(
0039     srcPFCandidateToVertexAssociations = 'pfCandidateToVertexAssociationForPFNoPUMEt'
0040 )
0041 pfNoPUMEtTask.add(pfMETcorrType0ForPFNoPUMEt)
0042
0043 jvfJetIdForPFNoPUMEt = cms.EDProducer("JVFJetIdProducer",
0044     srcJets = cms.InputTag('calibratedAK4PFJetsForPFNoPUMEt'),
0045     srcPFCandidates = cms.InputTag('particleFlow'),
0046     srcPFCandToVertexAssociations = cms.InputTag('pfCandidateToVertexAssociationForPFNoPUMEt'),
0047     srcHardScatterVertex = cms.InputTag('selectedPrimaryVertexHighestPtTrackSumForPFMEtCorrType0'),
0048     minTrackPt = cms.double(1.),
0049     dZcut = cms.double(0.2), # cm
0050     JVFcut = cms.double(0.75),
0051     neutralJetOption = cms.string("noPU")
0052 )
0053 pfNoPUMEtTask.add(jvfJetIdForPFNoPUMEt)
0054
0055 import RecoMET.METProducers.METSigParams_cfi as met_config
0056
0057 pfNoPUMEt = cms.EDProducer("NoPileUpPFMEtProducer",
0058     srcMEt = cms.InputTag('pfMet'),
0059     srcMEtCov = cms.InputTag(''), # NOTE: leave empty to take MET covariance matrix from reco::PFMET object //MM 08/29/14, bypass hardcoded as this variable has never been used so far
0060     srcPUSubMETDataJet = cms.InputTag('pfNoPUMEtData','jetInfos'),
0061     srcPUSubMETDataPFCands = cms.InputTag('pfNoPUMEtData','pfCandInfos'),
0062     srcLeptons = cms.VInputTag(), # NOTE: you need to set this to collections of electrons, muons and tau-jets
0063                                   #       passing the lepton reconstruction & identification criteria applied in your analysis
0064     srcPUSubMETDataJetLeptonMatch = cms.InputTag('pfNoPUMEtData','jetInfos'),
0065     srcPUSubMETDataPFCandsLeptonMatch = cms.InputTag('pfNoPUMEtData','pfCandInfos'),
0066     srcType0Correction = cms.InputTag('pfMETcorrType0ForPFNoPUMEt'),
0067     sfNoPUjets = cms.double(1.0),
0068     sfNoPUjetOffsetEnCorr = cms.double(0.0),
0069     sfPUjets = cms.double(1.0),
0070     sfNoPUunclChargedCands = cms.double(1.0),
0071     sfPUunclChargedCands = cms.double(1.0),
0072     sfUnclNeutralCands = cms.double(0.6),
0073     sfType0Correction = cms.double(1.0),
0074     sfLeptonIsoCones = cms.double(0.6),
0075     resolution = met_config.METSignificance_params,
0076     sfMEtCovMin = cms.double(0.6),
0077     sfMEtCovMax = cms.double(1.0),
0078     saveInputs = cms.bool(True),
0079     verbosity = cms.int32(0)
0080 )
0081 pfNoPUMEtTask.add(pfNoPUMEt)
0082 pfNoPUMEtSequence = cms.Sequence(pfNoPUMEtTask)