from CMGTools.Common.factories.cmgBaseJet_cfi import cmgBaseJet
from CMGTools.Common.skims.cmgPFJetSel_cfi import cmgPFJetSel
from CMGTools.Common.skims.cmgBaseJetSel_cfi import cmgBaseJetSel
from CMGTools.Common.skims.leadingCMGBaseJetSelector_cfi import leadingCMGBaseJetSelector
from CMGTools.Common.skims.leadingCMGPFJetSelector_cfi import leadingCMGPFJetSelector
from CMGTools.External.pujetidsequence_cff import puJetMva
# PF jets -------------------------------------------------------
# this pt cut is also going to be used to select GenJets,
# see CMGTools.Common.gen.generator.genJets_cff
ptCut = 'pt()>10'
cmgPFJetSel.cut = ptCut
# Jets for Type 1 MET corrections
cmgPFJetType1MET = cmgPFJetSel.clone(src = 'cmgPFJet',cut='pt() >= 10 && abs(eta()) < 4.7')
# PU JetID
cmgPUJetMva = puJetMva.clone()
cmgPUJetMva.jets = cms.InputTag('patJetsWithVar')
from CMGTools.Common.jetId_cff import *
pfJetSequence = cms.Sequence(
cmgPUJetMva +
cmgPFJet +
cmgPFJetSel +
cmgPFJetType1MET +
jetIdSequence
)
from RecoJets.JetProducers.ak5PFJets_cfi import *
from RecoJets.JetProducers.kt4PFJets_cfi import *
from RecoJets.JetProducers.ca4PFJets_cfi import *
from CMGTools.Common.skims.cmgPFJetSel_cfi import cmgPFJetSel
##################################
#####def input jet collection ####
##################################
#########playing with PU ID
cmgPFJetSelNoPU = cmgPFJetSel.clone(
src = 'cmgPFJetSel',
cut = 'puId("full") > 3'#'puId("full") > 5'#passPuJetId("full", PileupJetIdentifier::kMedium )' #does not work
)
#########playing with gen jet matching
#produce a collection of jets NOT matched to gen jet
deltaRJetGenJet = cms.EDProducer(
"DeltaRVetoProducerPFJet",
inputCollection = cms.InputTag('cmgPFJetSelCHS'),
MatchingParams = cms.VPSet(
cms.PSet(
vetoCollection=cms.InputTag("genJetSel"),
minDeltaR=cms.double(0.4),
removeMatchedObject=cms.bool(True)
)
),
verbose = cms.untracked.bool(False)
import FWCore.ParameterSet.Config as cms
from CMGTools.Utilities.metRecoilCorrection.recoilCorrectedMETTauMu_cfi import *
from CMGTools.Utilities.metRecoilCorrection.recoilCorrectedMETTauEle_cfi import *
from CMGTools.Utilities.metRecoilCorrection.recoilCorrectedMETMuEle_cfi import *
from CMGTools.Utilities.metRecoilCorrection.recoilCorrectedMETDiTau_cfi import *
from CMGTools.Common.generator.genWorZ_cfi import *
from SimGeneral.HepPDTESSource.pythiapdt_cfi import *
from CMGTools.Common.skims.cmgPFJetSel_cfi import cmgPFJetSel
cmgPFJetForRecoil = cmgPFJetSel.clone(
cut = 'pt()>30 && abs(eta)<4.5 && getSelection("cuts_looseJetId") && passLooseFullPuJetId()',
src = 'cmgPFJetSel')
metRecoilCorrectionInputSequence = cms.Sequence(
cmgPFJetForRecoil +
genWorZ
)
metRecoilCorrectionSequence = cms.Sequence(
metRecoilCorrectionInputSequence +
recoilCorrectedMETTauMu +
recoilCorrectedMETTauEle +
recoilCorrectedMETMuEle
# Recoil corrections for DiTau case not studied at all yet
#recoilCorrectedMETDiTau
)
import FWCore.ParameterSet.Config as cms
from CMGTools.Common.skims.cmgPFJetSel_cfi import cmgPFJetSel
jetIDJet = cmgPFJetSel.clone(
cut = cms.string( "getSelection(\"cuts_looseJetId\") && pt() > 30" )
)
process.schedule = cms.Schedule(process.hcallasereventfilter2012Path,*[p for p in process.schedule])
#### no AK5 pruned
if runPATCMG:
process.p.remove(process.PATCMGJetSequenceCHSpruned)
process.p.remove(process.razorMJStructureJetSel30)
#### recalibrated jets
if recalibrateJets:
from CMGTools.Common.miscProducers.cmgPFJetCorrector_cfi import cmgPFJetCorrector
process.cmgPFJetCor = cmgPFJetCorrector.clone(src='cmgPFJetSel',
payload='AK5PF')
process.cmgPFJetCorCHS = cmgPFJetCorrector.clone(src='cmgPFJetSelCHS',
payload='AK5PFchs')
from CMGTools.Common.skims.cmgPFJetSel_cfi import cmgPFJetSel
process.cmgPFJetSel = cmgPFJetSel.clone(src='cmgPFJetCor',
cut='pt()>30')
process.cmgPFJetSelCHS = cmgPFJetSel.clone(src='cmgPFJetCorCHS',
cut='pt()>30')
if runOnMC:
process.cmgPFJetCor.levels = cms.vstring('L1FastJet','L2Relative','L3Absolute')
process.cmgPFJetCorCHS.levels = cms.vstring('L1FastJet','L2Relative','L3Absolute')
else:
process.cmgPFJetCor.levels = cms.vstring('L1FastJet','L2Relative','L3Absolute','L2L3Residual')
process.cmgPFJetCorCHS.levels = cms.vstring('L1FastJet','L2Relative','L3Absolute','L2L3Residual')
process.p.insert(0,process.cmgPFJetSelCHS)
process.p.insert(0,process.cmgPFJetSel)
process.p.insert(0,process.cmgPFJetCorCHS)
process.p.insert(0,process.cmgPFJetCor)
#### evaluate speed
import FWCore.ParameterSet.Config as cms
from CMGTools.Common.skims.cmgPFJetSel_cfi import cmgPFJetSel
from CMGTools.Common.skims.cmgCandCount_cfi import cmgCandCount
#####Jet Selection
cmgPFJetSelPt = cmgPFJetSel.clone(
src = 'cmgPFJetSel',
cut = 'pt > 80. && abs(eta)<2.4')
cmgPFJetSelQuadJet = cmgCandCount.clone(
src = 'cmgPFJetSelPt',
minNumber = 4
)
cmgPFJetSelDiJet = cmgCandCount.clone(
src = 'cmgPFJetSelPt',
minNumber = 2
)
######MET selection
cmgPFMETSelPt = cms.EDFilter(
"CandViewSelector",
src = cms.InputTag("cmgPFMET"),
cut = cms.string( "pt()>200." )
)
cmgPFMETSelCut = cmgCandCount.clone(
src = 'cmgPFMETSelPt',
minNumber = 1
)
skimmingTopTupleTrgSequence = cms.Sequence(
import FWCore.ParameterSet.Config as cms
# from CMGTools.H2TauTau.generator.metRecoilCorrection.recoilCorrectedMETTauMu_cfi import *
# from CMGTools.H2TauTau.generator.metRecoilCorrection.recoilCorrectedMETTauEle_cfi import *
# from CMGTools.H2TauTau.generator.metRecoilCorrection.recoilCorrectedMETMuEle_cfi import *
# from CMGTools.H2TauTau.generator.metRecoilCorrection.recoilCorrectedMETDiTau_cfi import *
from CMGTools.Utilities.generator.metRecoilCorrection.recoilCorrectedMETTauMu_cfi import *
from CMGTools.Utilities.generator.metRecoilCorrection.recoilCorrectedMETTauEle_cfi import *
from CMGTools.Utilities.generator.metRecoilCorrection.recoilCorrectedMETMuEle_cfi import *
from CMGTools.Utilities.generator.metRecoilCorrection.recoilCorrectedMETDiTau_cfi import *
from CMGTools.Common.generator.genWorZ_cfi import *
from SimGeneral.HepPDTESSource.pythiapdt_cfi import *
from CMGTools.Common.skims.cmgPFJetSel_cfi import cmgPFJetSel
cmgPFJetForRecoil = cmgPFJetSel.clone(
cut=
'pt()>30 && abs(eta)<4.5 && getSelection("cuts_looseJetId") && passLooseFullPuJetId()',
src='cmgPFJetSel')
metRecoilCorrectionInputSequence = cms.Sequence(cmgPFJetForRecoil + genWorZ)
metRecoilCorrectionSequence2012 = cms.Sequence(
metRecoilCorrectionInputSequence + recoilCorrectedMETTauMu2012 +
recoilCorrectedMETTauEle2012 + recoilCorrectedMETMuEle2012
# Recoil corrections for DiTau case not studied at all yet
#recoilCorrectedMETDiTau2012
)
import FWCore.ParameterSet.Config as cms
from CMGTools.Utilities.metRecoilCorrection.recoilCorrectedMETTauMu_cfi import *
from CMGTools.Utilities.metRecoilCorrection.recoilCorrectedMETTauEle_cfi import *
from CMGTools.Utilities.metRecoilCorrection.recoilCorrectedMETMuEle_cfi import *
from CMGTools.Utilities.metRecoilCorrection.recoilCorrectedMETDiTau_cfi import *
from CMGTools.Common.generator.genWorZ_cfi import *
from SimGeneral.HepPDTESSource.pythiapdt_cfi import *
from CMGTools.Common.skims.cmgPFJetSel_cfi import cmgPFJetSel
cmgPFJetForRecoilPresel = cmgPFJetSel.clone(
cut = 'pt()>30 && abs(eta)<4.7 && getSelection("cuts_looseJetId")',
src = 'cmgPFJetSel')
from CMGTools.Common.skims.cmgPFJetPUIDSel_cfi import cmgPFJetPUIDSel
cmgPFJetForRecoil = cmgPFJetPUIDSel.clone(
src = 'cmgPFJetForRecoilPresel')
metRecoilCorrectionInputSequence = cms.Sequence(
cmgPFJetForRecoilPresel +
cmgPFJetForRecoil +
genWorZ
)
metRecoilCorrectionSequence = cms.Sequence(
metRecoilCorrectionInputSequence +
recoilCorrectedMETTauMu +
recoilCorrectedMETTauEle +
recoilCorrectedMETMuEle
# Recoil corrections for DiTau case not studied at all yet
import FWCore.ParameterSet.Config as cms
from CMGTools.Common.skims.cmgPFJetSel_cfi import cmgPFJetSel
jetIDJet = cmgPFJetSel.clone(
cut=cms.string("getSelection(\"cuts_looseJetId\") && pt() > 30"))