# parton and gen jet matching from CommonTools.ParticleFlow.genForPF2PAT_cff import * from PhysicsTools.PatAlgos.mcMatchLayer0.jetMatch_cfi import * patJetPartonMatch.src = jetSource patJetGenJetMatch.src = jetSource patJetGenJetMatch.matched = 'ak5GenJetsNoNu' from PhysicsTools.PatAlgos.mcMatchLayer0.jetFlavourId_cff import * patJetPartonAssociation.jets = jetSource jetsPtGt1Cut = '(neutralHadronEnergy())/(correctedJet(0).pt()/pt()*energy()) < 0.99 && (neutralEmEnergy()/(correctedJet(0).pt()/pt()*energy())) < 0.99 && (nConstituents()) > 1 && ((abs(eta()) < 2.4 && chargedHadronEnergy()/(correctedJet(0).pt()/pt()*energy()) > 0 && chargedEmEnergy() /(correctedJet(0).pt()/pt()*energy()) < 0.99 && chargedMultiplicity() > 0) || abs(eta()) > 2.4) ' if not isNewerThan('CMSSW_5_2_0'): # addtl cut needed due to different MVA MET training in 44X jetsPtGt1Cut = ' && '.join([jetsPtGt1Cut,'pt()>1']) jetsPtGt1 = cmgCandSel.clone( src = 'patJets', cut = jetsPtGt1Cut ) from CMGTools.Common.miscProducers.collectionSize.candidateSize_cfi import candidateSize nJetsPtGt1 = candidateSize.clone( src = 'jetsPtGt1' ) # jet extender patJetsWithVar = cms.EDProducer('JetExtendedProducer', jets = cms.InputTag('selectedPatJets'), vertices = cms.InputTag('goodOfflinePrimaryVertices'), #debug = cms.untracked.bool(True), payload = cms.string('AK5PF') ) outPFCand = cms.EDProducer('VbfHbbPFCandOutOfJets', jets = cms.InputTag('patJetsWithVar'), vtx = cms.InputTag('goodOfflinePrimaryVertices'),
import FWCore.ParameterSet.Config as cms from CMGTools.Common.factories.cmgBaseCandMET_cfi import cmgBaseCandMET from CMGTools.Common.factories.cmgBaseMET_cfi import cmgBaseMET from CMGTools.Common.factories.cmgPFMETScaler_cfi import cmgPFMETType1Corrected from CMGTools.Common.skims.cmgCandSel_cfi import cmgCandSel #PFMET from pat::MET cmgPFMET = cmgBaseMET.clone() cmgPFMET.cfg.inputCollection = "patMETsAK5" cmgPFMETSel = cmgCandSel.clone(src='cmgPFMET') # MHT from PFJets, pt threshold 30 cmgMHTPFJet30 = cmgBaseCandMET.clone() cmgMHTPFJet30.cfg.inputCollection = "cmgPFBaseJet" cmgMHTPFJet30.cfg.ptThreshold = 30.0 cmgMHTPFJet30Sel = cmgCandSel.clone(src='cmgMHTPFJet30') # Calo MET from pat::MET cmgCaloMET = cmgBaseMET.clone() cmgCaloMET.cfg.inputCollection = "patMETs" # MHT from CaloJets, pt threshold 30 cmgMHTCaloJet30 = cmgBaseCandMET.clone() cmgMHTCaloJet30.cfg.inputCollection = "cmgCaloBaseJet" cmgMHTCaloJet30.cfg.ptThreshold = 30.0 # MET from PFCandidates # cmgMETPFCandidates = cmgBaseCandMET.clone() # cmgMETPFCandidates.cfg.inputCollection = cms.InputTag("particleFlow")
import FWCore.ParameterSet.Config as cms from CMGTools.Common.factories.cmgBaseCandMET_cfi import cmgBaseCandMET from CMGTools.Common.factories.cmgBaseMET_cfi import cmgBaseMET from CMGTools.Common.factories.cmgPFMETScaler_cfi import cmgPFMETType1Corrected from CMGTools.Common.skims.cmgCandSel_cfi import cmgCandSel #PFMET from pat::MET cmgPFMET = cmgBaseMET.clone() cmgPFMET.cfg.inputCollection = "patMETsAK5" cmgPFMETSel = cmgCandSel.clone( src = 'cmgPFMET' ) # MHT from PFJets, pt threshold 30 cmgMHTPFJet30 = cmgBaseCandMET.clone() cmgMHTPFJet30.cfg.inputCollection = "cmgPFBaseJet" cmgMHTPFJet30.cfg.ptThreshold = 30.0 cmgMHTPFJet30Sel = cmgCandSel.clone( src = 'cmgMHTPFJet30' ) # Calo MET from pat::MET cmgCaloMET = cmgBaseMET.clone() cmgCaloMET.cfg.inputCollection = "patMETs" # MHT from CaloJets, pt threshold 30 cmgMHTCaloJet30 = cmgBaseCandMET.clone() cmgMHTCaloJet30.cfg.inputCollection = "cmgCaloBaseJet" cmgMHTCaloJet30.cfg.ptThreshold = 30.0
# parton and gen jet matching from CommonTools.ParticleFlow.genForPF2PAT_cff import * from PhysicsTools.PatAlgos.mcMatchLayer0.jetMatch_cfi import * patJetPartonMatch.src = jetSource patJetGenJetMatch.src = jetSource patJetGenJetMatch.matched = 'ak5GenJetsNoNu' from PhysicsTools.PatAlgos.mcMatchLayer0.jetFlavourId_cff import * patJetPartonAssociation.jets = jetSource jetsPtGt1Cut = '(neutralHadronEnergy())/(correctedJet(0).pt()/pt()*energy()) < 0.99 && (neutralEmEnergy()/(correctedJet(0).pt()/pt()*energy())) < 0.99 && (nConstituents()) > 1 && ((abs(eta()) < 2.4 && chargedHadronEnergy()/(correctedJet(0).pt()/pt()*energy()) > 0 && chargedEmEnergy() /(correctedJet(0).pt()/pt()*energy()) < 0.99 && chargedMultiplicity() > 0) || abs(eta()) > 2.4) ' if not isNewerThan('CMSSW_5_2_0'): # addtl cut needed due to different MVA MET training in 44X jetsPtGt1Cut = ' && '.join([jetsPtGt1Cut,'pt()>1']) jetsPtGt1 = cmgCandSel.clone( src = 'patJets', cut = jetsPtGt1Cut ) from CMGTools.Common.miscProducers.collectionSize.candidateSize_cfi import candidateSize nJetsPtGt1 = candidateSize.clone( src = 'jetsPtGt1' ) # QG Tagger QGTagger = cms.EDProducer('QGTagger', srcJets = cms.InputTag('selectedPatJets'), isPatJet = cms.untracked.bool(True), useCHS = cms.untracked.bool(False), srcRho = cms.InputTag('kt6PFJets','rho'), srcRhoIso = cms.InputTag('kt6PFJetsIsoQG','rho'), ) # jet extender patJetsWithVar = cms.EDProducer('JetExtendedProducer',
"keep++ pdgId = 6",##top daughters "keep++ pdgId = -6", ) ) genParticlesFromW = cms.EDProducer( "GenParticlePruner", src = cms.InputTag("topDaughters"), select = cms.vstring( "drop * ", "keep+ pdgId = 24",##W daughters "keep+ pdgId = -24", ) ) leptonicDecay = cmgCandSel.clone( src = 'genParticlesFromW', cut = '(abs(pdgId()) == 12 || abs(pdgId()) == 14 || abs(pdgId()) == 16) ' ) countLeptonicDecay = cmgCandCount.clone( src = 'leptonicDecay', minNumber = 1 ) countDiLeptonicDecay = cmgCandCount.clone( src = 'leptonicDecay', minNumber = 2 )