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