if not config["RUNONMC"]:
process.jec.toGet[0].tag = cms.string('JetCorrectorParametersCollection_Spring16_25nsV6_DATA_AK4PFchs')
process.jec.toGet[1].tag = cms.string('JetCorrectorParametersCollection_Spring16_25nsV6_DATA_AK8PFchs')
process.jec.toGet[2].tag = cms.string('JetCorrectorParametersCollection_Spring16_25nsV6_DATA_AK8PFPuppi')
process.jec.connect = cms.string('sqlite:JEC/Spring16_25nsV6_DATA.db')
else:
sys.exit("This config file expects to run on data, but an option to run on MC was received.")
process.es_prefer_jec = cms.ESPrefer('PoolDBESSource','jec')
####### Add AK8 GenJets ##########
if config["ADDAK8GENJETS"]:
from RecoJets.Configuration.RecoGenJets_cff import ak8GenJets
process.ak8GenJets = ak8GenJets.clone(src = 'packedGenParticles')
process.NjettinessGenAK8 = cms.EDProducer("NjettinessAdder",
src=cms.InputTag("ak8GenJets"),
Njets=cms.vuint32(1,2,3,4), # compute 1-, 2-, 3-, 4- subjettiness
# variables for measure definition :
measureDefinition = cms.uint32( 0 ), # CMS default is normalized measure
beta = cms.double(1.0), # CMS default is 1
R0 = cms.double( 0.8 ), # CMS default is jet cone size
Rcutoff = cms.double( 999.0), # not used by default
# variables for axes definition :
axesDefinition = cms.uint32( 6 ), # CMS default is 1-pass KT axes
nPass = cms.int32(999), # not used by default
akAxesR0 = cms.double(-999.0) # not used by default
)
import FWCore.ParameterSet.Config as cms
from RecoJets.Configuration.GenJetParticles_cff import genParticlesForJets
from RecoJets.Configuration.RecoGenJets_cff import ak8GenJets
from RecoJets.JetProducers.SubJetParameters_cfi import SubJetParameters
ak8GenJetsPruned = ak8GenJets.clone(
SubJetParameters,
usePruning = cms.bool(True),
useExplicitGhosts = cms.bool(True),
writeCompound = cms.bool(True),
jetCollInstanceName=cms.string("SubJets")
)
ak8GenJetsPrunedLinks = cms.EDProducer("RecoJetDeltaRValueMapProducer",
src = cms.InputTag("ak8GenJets"),
matched = cms.InputTag("ak8GenJetsPruned"),
distMax = cms.double(0.8),
value = cms.string('mass')
)
from PhysicsTools.PatAlgos.producersLayer1.jetProducer_cfi import patJets as patJets
goodJets = patJets.clone(
jetSource = "ak8GenJets",
embedCaloTowers = cms.bool(False),
embedPFCandidates = cms.bool(False),
addJetCorrFactors = cms.bool(False),
addBTagInfo = cms.bool(False),
addDiscriminators = cms.bool(False),
addAssociatedTracks = cms.bool(False),
import FWCore.ParameterSet.Config as cms
from RecoJets.Configuration.GenJetParticles_cff import genParticlesForJets
from RecoJets.Configuration.RecoGenJets_cff import ak8GenJets
from RecoJets.JetProducers.SubJetParameters_cfi import SubJetParameters
ak8GenJetsPruned = ak8GenJets.clone(SubJetParameters,
usePruning=cms.bool(True),
useExplicitGhosts=cms.bool(True),
writeCompound=cms.bool(True),
jetCollInstanceName=cms.string("SubJets"))
ak8GenJetsPrunedLinks = cms.EDProducer(
"RecoJetDeltaRValueMapProducer",
src=cms.InputTag("ak8GenJets"),
matched=cms.InputTag("ak8GenJetsPruned"),
distMax=cms.double(0.8),
value=cms.string('mass'))
from PhysicsTools.PatAlgos.producersLayer1.jetProducer_cfi import patJets as patJets
goodJets = patJets.clone(jetSource="ak8GenJets",
embedCaloTowers=cms.bool(False),
embedPFCandidates=cms.bool(False),
addJetCorrFactors=cms.bool(False),
addBTagInfo=cms.bool(False),
addDiscriminators=cms.bool(False),
addAssociatedTracks=cms.bool(False),
addJetCharge=cms.bool(False),
addJetID=cms.bool(False),
addGenPartonMatch=cms.bool(False),
# ####### Logger ##########
process.load("FWCore.MessageLogger.MessageLogger_cfi")
process.MessageLogger.cerr.threshold = 'INFO'
process.MessageLogger.categories.append('Ntuple')
process.MessageLogger.cerr.INFO = cms.untracked.PSet(
limit=cms.untracked.int32(1))
process.MessageLogger.cerr.FwkReport.reportEvery = 1000
####### Redo Jet clustering sequence ##########
from RecoJets.Configuration.GenJetParticles_cff import genParticlesForJets
from RecoJets.Configuration.RecoGenJets_cff import ak8GenJets
process.genParForJets = genParticlesForJets.clone(
src=cms.InputTag("genParticles"))
process.ak8GenJets = ak8GenJets.clone(src='genParForJets')
process.NjettinessGenAK8 = cms.EDProducer(
"NjettinessAdder",
src=cms.InputTag("ak8GenJets"),
# compute 1-, 2-, 3-, 4- subjettiness
Njets=cms.vuint32(1, 2, 3, 4),
# variables for measure definition:
# CMS default is normalized measure
measureDefinition=cms.uint32(0),
# CMS default is 1
beta=cms.double(1.0),
# CMS default is jet cone size
R0=cms.double(0.8),
# not used by default
Rcutoff=cms.double(999.0),
