Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-02-14 13:33:08

0001 import FWCore.ParameterSet.Config as cms
0002 
0003 process = cms.Process("EnergyShift")
0004 
0005 ## Declare input
0006 from PhysicsTools.PatExamples.samplesCERN_cff import *
0007 
0008 process.source = cms.Source("PoolSource",
0009   fileNames = ttbarJets
0010 )
0011 
0012 process.maxEvents = cms.untracked.PSet(
0013     input = cms.untracked.int32( -1 )
0014 )
0015 
0016 ## Message logger configuration
0017 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0018 
0019 ## Configure jet energy scaling module
0020 process.load("PhysicsTools.PatExamples.JetEnergyShift_cfi")
0021 process.scaledJets.scaleFactor =  1.0
0022 
0023 ## Select good jets
0024 from PhysicsTools.PatAlgos.selectionLayer1.jetSelector_cfi import *
0025 process.goodJets = selectedPatJets.clone(
0026     src="scaledJets:cleanPatJets",
0027     cut = 'abs(eta) < 3 & pt > 30. &'
0028     'emEnergyFraction > 0.01       &'
0029     'jetID.fHPD < 0.98             &'
0030     'jetID.n90Hits > 1'    
0031     )
0032 
0033 ## Analyze jets
0034 process.load("PhysicsTools.PatExamples.PatJetAnalyzer_cfi")
0035 process.analyzePatJets.src = 'goodJets'
0036 
0037 ## Define output file
0038 process.TFileService = cms.Service("TFileService",
0039   fileName = cms.string('analyzeEnergyShiftedJets.root')
0040 )
0041 
0042 process.p = cms.Path(
0043     process.scaledJets *
0044     process.goodJets   * 
0045     process.analyzePatJets
0046 )