def runMVAMet(process, isMC, leptons): from RecoJets.JetProducers.ak4PFJets_cfi import ak4PFJets if not hasattr(process, "ak4PFJets"): setattr(process, "ak4PFJets", ak4PFJets.clone(src=cms.InputTag('packedPFCandidates'))) if not hasattr(process, "ak4PFL1FastL2L3Corrector"): process.load('JetMETCorrections.Configuration.JetCorrectors_cff') process.calibratedAK4PFJetsForPFMVAMEt = calibratedAK4PFJetsForPFMVAMEt.clone( correctors=cms.vstring("ak4PFL1FastL2L3")) from JetMETCorrections.Configuration.JetCorrectionServices_cff import ak4PFL1FastL2L3, ak4PFL1FastL2L3Residual, ak4PFL1Fastjet, ak4PFL2Relative, ak4PFL3Absolute, ak4PFResidual if not isMC: process.ak4PFL1FastL2L3Residual = ak4PFL1FastL2L3Residual.clone() process.ak4PFResidual = ak4PFResidual.clone() process.calibratedAK4PFJetsForPFMVAMEt.correctors = cms.vstring( "ak4PFL1FastL2L3Residual") else: process.ak4PFL1FastL2L3 = ak4PFL1FastL2L3.clone() process.ak4PFL1Fastjet = ak4PFL1Fastjet.clone() process.ak4PFL2Relative = ak4PFL2Relative.clone() process.ak4PFL3Absolute = ak4PFL3Absolute.clone() from RecoJets.JetProducers.PileupJetID_cfi import pileupJetIdEvaluator from RecoJets.JetProducers.PileupJetIDParams_cfi import JetIdParams process.puJetIdForPFMVAMEt = pileupJetIdEvaluator.clone( produceJetIds=cms.bool(True), jets=cms.InputTag("calibratedAK4PFJetsForPFMVAMEt"), vertexes=cms.InputTag("offlineSlimmedPrimaryVertices"), jec=cms.string("AK4PF")) process.mvaMET = pfMVAMEt.clone( srcPFCandidates=cms.InputTag('packedPFCandidates'), srcVertices=cms.InputTag('offlineSlimmedPrimaryVertices'), srcLeptons=cms.VInputTag(leptons), inputFileNames=cms. PSet(U=cms.FileInPath( 'RecoMET/METPUSubtraction/data/gbru_7_4_X_miniAOD_50NS_July2015.root' ), DPhi=cms.FileInPath( 'RecoMET/METPUSubtraction/data/gbrphi_7_4_X_miniAOD_50NS_July2015.root' ), CovU1=cms.FileInPath( 'RecoMET/METPUSubtraction/data/gbru1cov_7_4_X_miniAOD_50NS_July2015.root' ), CovU2=cms.FileInPath( 'RecoMET/METPUSubtraction/data/gbru2cov_7_4_X_miniAOD_50NS_July2015.root' )))
from RecoJets.JetProducers.PileupJetID_cfi import pileupJetIdEvaluator from RecoJets.JetProducers.PileupJetIDParams_cfi import JetIdParams puJetIdForPFMVAMEt = pileupJetIdEvaluator.clone( algos=cms.VPSet( cms.PSet( tmvaVariables=cms.vstring("nvtx", "jetPt", "jetEta", "jetPhi", "dZ", "beta", "betaStar", "nCharged", "nNeutrals", "dR2Mean", "ptD", "frac01", "frac02", "frac03", "frac04", "frac05"), tmvaWeights=cms.string( "RecoJets/JetProducers/data/TMVAClassificationCategory_JetID_MET_53X_Dec2012.weights.xml.gz" ), tmvaMethod=cms.string("JetID"), tmvaSpectators=cms.vstring(), JetIdParams=JetIdParams, impactParTkThreshold=cms.double(0.), version=cms.int32(-1), cutBased=cms.bool(False), label=cms.string("full"))), produceJetIds=cms.bool(True), runMvas=cms.bool(True), jets=cms.InputTag( "calibratedAK4PFJetsForPFMVAMEt"), #calibratedAK4PFJetsForPFMVAMEt applyJec=cms.bool(True), inputIsCorrected=cms.bool(True), jec=cms.string("AK4PF"), ) pfMVAMEt = cms.EDProducer( "PFMETProducerMVA",
puJetIdForPFMVAMEt = pileupJetIdEvaluator.clone( algos = cms.VPSet( cms.PSet( tmvaVariables = cms.vstring( "nvtx", "jetPt", "jetEta", "jetPhi", "dZ", "beta", "betaStar", "nCharged", "nNeutrals", "dR2Mean", "ptD", "frac01", "frac02", "frac03", "frac04", "frac05" ), etaBinnedWeights = cms.bool(False), tmvaWeights = cms.string("RecoJets/JetProducers/data/TMVAClassificationCategory_JetID_MET_53X_Dec2012.weights.xml.gz"), tmvaMethod = cms.string("JetID"), tmvaSpectators = cms.vstring(), JetIdParams = JetIdParams, impactParTkThreshold = cms.double(0.), version = cms.int32(-1), cutBased = cms.bool(False), label = cms.string("full") ) ), produceJetIds = cms.bool(True), runMvas = cms.bool(True), jets = cms.InputTag("calibratedAK4PFJetsForPFMVAMEt"),#calibratedAK4PFJetsForPFMVAMEt applyJec = cms.bool(True), inputIsCorrected = cms.bool(True), jec = cms.string("AK4PF"), )