JetAnalyzers/test/runL2L3L5L7JetCorrectionExample_cfg.py

0001 # PYTHON configuration file.
0002 # Description:  Example of applying L2+L3+L5+L7 jet corrections.
0003 # Author: K. Kousouris
0004 # Date:  02 - September - 2009
0005 import FWCore.ParameterSet.Config as cms
0006
0007 process = cms.Process("Ana")
0008 process.load("FWCore.MessageService.MessageLogger_cfi")
0009 #############   Set the number of events #############
0010 process.maxEvents = cms.untracked.PSet(
0011     input = cms.untracked.int32(100)
0012 )
0013 #############   Define the source file ###############
0014 process.source = cms.Source("PoolSource",
0015     fileNames = cms.untracked.vstring('/store/relval/CMSSW_3_1_2/RelValQCD_FlatPt_15_3000/GEN-SIM-RECO/MC_31X_V3-v1/0007/9E83A122-E978-DE11-9D04-001D09F23C73.root')
0016 )
0017 #############   Define the L2 correction service #####
0018 process.L2RelativeJetCorrector = cms.ESSource("L2RelativeCorrectionService",
0019     tagName = cms.string('Summer09_L2Relative_SC5Calo'),
0020     label = cms.string('L2RelativeJetCorrector')
0021 )
0022 #############   Define the L3 correction service #####
0023 process.L3AbsoluteJetCorrector = cms.ESSource("L3AbsoluteCorrectionService",
0024     tagName = cms.string('Summer09_L3Absolute_SC5Calo'),
0025     label = cms.string('L3AbsoluteJetCorrector')
0026 )
0027 #############   Define the L5 correction service #####
0028 process.L5JetCorrector = cms.ESSource("L5FlavorCorrectionService",
0029     section = cms.string('qJ'),
0030     tagName = cms.string('L5Flavor_IC5'),# IMPORTANT: the L5 correction was derived from IC5 but it is the same for all algos
0031     label = cms.string('L5FlavorJetCorrector')
0032 )
0033 #############   Define the L7 correction service #####
0034 process.L7JetCorrector = cms.ESSource("L7PartonCorrectionService",
0035     section = cms.string('qJ'),
0036     tagName = cms.string('L7Parton_SC5'),
0037     label = cms.string('L7PartonJetCorrector')
0038 )
0039 #############   Define the chain corrector service ###
0040 process.L2L3L5L7JetCorrector = cms.ESSource("JetCorrectionServiceChain",
0041     correctors = cms.vstring('L2RelativeJetCorrector','L3AbsoluteJetCorrector','L5FlavorJetCorrector','L7PartonJetCorrector'),
0042     label = cms.string('L2L3L5L7JetCorrector')
0043 )
0044 #############   Define the chain corrector module ####
0045 process.L2L3L5L7CorJet = cms.EDProducer("CaloJetCorrectionProducer",
0046     src = cms.InputTag("sisCone5CaloJets"),
0047     correctors = cms.vstring('L2L3L5L7JetCorrector')
0048 )
0049 # set the record's IOV. Must be defined once. Choose ANY correction service. #
0050 process.prefer("L2L3L5L7JetCorrector")
0051 #############   Plots of corrected Jet collection ####
0052 process.plots = cms.EDAnalyzer("CaloJetPlotsExample",
0053     JetAlgorithm  = cms.string('L2L3L5L7CorJet'),
0054     HistoFileName = cms.string('L2L3L5L7CaloJetCorExample.root'),
0055     NJets         = cms.int32(100)
0056 )
0057 #############   Path       ###########################
0058 process.p = cms.Path(process.L2L3L5L7CorJet * process.plots)
0059 #############   Format MessageLogger #################
0060 process.MessageLogger.cerr.FwkReport.reportEvery = 10