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