)
from Dead.Puppi.Puppi_PhaseII_cff_v2 import *
process.load("Dead.Puppi.Puppi_PhaseII_cff_v2")
# sequence = cms.Sequence(sequence * )
process.load("PhysicsTools.PatAlgos.selectionLayer1.leptonCountFilter_cfi")
process.count2Muons = process.countPatLeptons.clone(countElectrons=False, muonSource="selectedIsoMuons", minNumber=2)
process.count4Muons = process.countPatLeptons.clone(countElectrons=False, muonSource="selectedIsoMuons", minNumber=4)
process.load("JetCorrection_cff")
process.ak4PUPPIJets = ak5PFJets.clone(rParam=0.4, src=cms.InputTag("puppi", "Puppi"))
if doJEC:
process.ak4PUPPIL1Fastjet = cms.ESProducer(
"L1FastjetCorrectionESProducer",
era=cms.string("Summer11"),
level=cms.string("L1FastJet"),
algorithm=cms.string("AK4PUPPI"),
section=cms.string(""),
srcRho=cms.InputTag("kt6PFJets", "rho"),
useCondDB=cms.untracked.bool(True),
)
era = "PhaseII_Shashlik140PU"
# process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
# process.GlobalTag.globaltag = ""
tag=cms.string("JetCorrectorParametersCollection_" + era + "_AK4PF"),
label=cms.untracked.string("AK4PF"),
),
cms.PSet(
record=cms.string("JetCorrectionsRecord"),
tag=cms.string("JetCorrectorParametersCollection_" + era + "_AK4PFchs"),
label=cms.untracked.string("AK4PFchs"),
),
),
)
process.es_prefer_jec = cms.ESPrefer("PoolDBESSource", "jec")
process.TFileService = cms.Service("TFileService", fileName=cms.string("vbfjetanal.root"))
# apparently no CHS pf collection available?
process.ak25PFJets = ak5PFJets.clone(rParam=0.25)
process.gpfilt = cms.EDFilter(
"GenParticleFilter", GenParticleLabel=cms.InputTag("genParticles"), excludeIDlist=cms.vint32(-15, -11, 11, 15)
)
from MyEDmodules.phase2upVBFjetAnal.phase2upVBFjetAnal_cfi import *
process.vbfanaltight = vbfanaltight.clone()
process.vbfanalloose = vbfanalloose.clone()
process.vbfanaltight.benchmarkParams.isSignalDoTruthMatch = cms.bool(isSignal)
process.vbfanaltight.benchmarkParams.pfjBenchmarkDebug = cms.bool(debugPrintout)
process.vbfanaltight.RecoJetCoreLabel = cms.InputTag("ak25PFJets")
process.vbfanalloose.benchmarkParams.isSignalDoTruthMatch = cms.bool(isSignal)
process.goodOfflinePrimaryVertices = cms.EDFilter(
"PrimaryVertexObjectFilter",
filterParams=pvSelector.clone(minNdof=cms.double(4.0), maxZ=cms.double(24.0)),
src=cms.InputTag("offlinePrimaryVertices"),
)
# Compute the mean pt per unit area (rho) from the PFCHS inputs
from RecoJets.Configuration.RecoPFJets_cff import kt6PFJets, ak5PFJets
process.kt6PFJetsPFlow = kt6PFJets.clone(
src=cms.InputTag("pfNoElectron" + postfix), doAreaFastjet=cms.bool(True), doRhoFastjet=cms.bool(True)
)
process.patJetCorrFactorsPFlow.rho = cms.InputTag("kt6PFJetsPFlow", "rho")
process.ak5PFJetsPileUp = ak5PFJets.clone(src=cms.InputTag("pfPileUp" + postfix))
"""
process.pileUpN1PrimaryVertices = cms.EDFilter("RecoTauPileUpVertexSelector",
filter = cms.bool(False),
src = cms.InputTag('goodOfflinePrimaryVertices'),
minTrackSumPt = cms.double(0.)
#minTrackSumPt = cms.double(5.)
)
process.pileUpN2PrimaryVertices = process.pileUpN1PrimaryVertices.clone(
src = cms.InputTag('pileUpN1PrimaryVertices')
)
process.pileUpN3PrimaryVertices = process.pileUpN1PrimaryVertices.clone(
src = cms.InputTag('pileUpN2PrimaryVertices')
from PhysicsTools.SelectorUtils.pvSelector_cfi import pvSelector
process.goodOfflinePrimaryVertices = cms.EDFilter(
"PrimaryVertexObjectFilter",
filterParams=pvSelector.clone(minNdof=cms.double(4.0),
maxZ=cms.double(24.0)),
src=cms.InputTag('offlinePrimaryVertices'))
# Compute the mean pt per unit area (rho) from the PFCHS inputs
from RecoJets.Configuration.RecoPFJets_cff import kt6PFJets, ak5PFJets
process.kt6PFJetsPFlow = kt6PFJets.clone(src=cms.InputTag('pfNoElectron' +
postfix),
doAreaFastjet=cms.bool(True),
doRhoFastjet=cms.bool(True))
process.patJetCorrFactorsPFlow.rho = cms.InputTag("kt6PFJetsPFlow", "rho")
process.ak5PFJetsPileUp = ak5PFJets.clone(src=cms.InputTag('pfPileUp' +
postfix))
"""
process.pileUpN1PrimaryVertices = cms.EDFilter("RecoTauPileUpVertexSelector",
filter = cms.bool(False),
src = cms.InputTag('goodOfflinePrimaryVertices'),
minTrackSumPt = cms.double(0.)
#minTrackSumPt = cms.double(5.)
)
process.pileUpN2PrimaryVertices = process.pileUpN1PrimaryVertices.clone(
src = cms.InputTag('pileUpN1PrimaryVertices')
)
process.pileUpN3PrimaryVertices = process.pileUpN1PrimaryVertices.clone(
src = cms.InputTag('pileUpN2PrimaryVertices')
)
# ------------------------------------------------------------------------------
# https://twiki.cern.ch/twiki/bin/view/CMSPublic/WorkBookJetEnergyCorrections
# https://twiki.cern.ch/twiki/bin/view/CMS/JetMETAlgorithmsReconstruction
# https://twiki.cern.ch/twiki/bin/view/CMS/JetResolution
# https://twiki.cern.ch/twiki/bin/view/CMS/JECDataMC
# https://twiki.cern.ch/twiki/bin/view/CMS/JECUncertaintySources
# https://twiki.cern.ch/twiki/bin/view/CMS/JetID
# ------------------------------------------------------------------------------
from RecoJets.Configuration.RecoPFJets_cff import \
ak5PFJets, ak5PFJetsCHS, kt6PFJets, ak8PFJets, ak8PFJetsCHS, \
ak8PFJetsCHSPruned, ak8PFJetsCHSFiltered, ak8PFJetsCHSTrimmed, \
ca8PFJetsCHS, ca8PFJetsCHSPruned, ca8PFJetsCHSFiltered, ca8PFJetsCHSTrimmed, \
ca15PFJetsCHSMassDropFiltered, ca15PFJetsCHSFiltered
#setattr(process, "ak4PFJets"+postfix, ak5PFJets.clone( rParam = 0.4 ))
process.ak4PFJets = ak5PFJets.clone(rParam=0.4)
process.kt6PFJets = kt6PFJets
process.ca15PFJetsCHS = ca8PFJetsCHS.clone(rParam=1.5)
process.ca8PFJetsCHSPruned = ca8PFJetsCHSPruned.clone(
doAreaFastjet=cms.bool(True))
process.ca15PFJetsCHSFiltered = ca15PFJetsCHSFiltered.clone(
doAreaFastjet=cms.bool(True),
## Uncomment the following three lines to turn 'Filtered' into 'MassDropFiltered'
#useMassDropTagger = cms.bool(True),
#muCut = cms.double(0.667),
#yCut = cms.double(0.08),
)
if runOnMC:
# Create GenJets
from RecoJets.Configuration.RecoGenJets_cff import \