# 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
)
