Modules/test/JetCorrectorDemo_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("myprocess")
0004 process.TFileService=cms.Service("TFileService",fileName=cms.string('JECplots.root'))
0005 ##-------------------- Communicate with the DB -----------------------
0006 process.load('Configuration.StandardSequences.Services_cff')
0007 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0008 from Configuration.AlCa.GlobalTag import GlobalTag
0009 process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:run2_mc')
0010
0011 ##-------------------- Import the JEC services -----------------------
0012 process.load('JetMETCorrections.Configuration.JetCorrectors_cff')
0013
0014 ##-------------------- Define the source  ----------------------------
0015 process.maxEvents = cms.untracked.PSet(
0016         input = cms.untracked.int32(1)
0017 )
0018 process.source = cms.Source("PoolSource",
0019     fileNames = cms.untracked.vstring('/store/relval/CMSSW_7_2_0_pre7/RelValProdTTbar/AODSIM/PRE_STA72_V4-v1/00000/3E58BB46-BD4B-E411-B2EC-002618943856.root')
0020 )
0021
0022 ##-------------------- User analyzer  --------------------------------
0023 process.ak4pfl2l3Residual  = cms.EDAnalyzer('JetCorrectorDemo',
0024     JetCorrector             = cms.InputTag('ak4PFL2L3ResidualCorrector'),
0025     UncertaintyTag           = cms.string('Uncertainty'),
0026     UncertaintyFile          = cms.string(''),
0027     PayloadName              = cms.string('AK4PF'),
0028     NHistoPoints             = cms.int32(10000),
0029     NGraphPoints             = cms.int32(500),
0030     EtaMin                   = cms.double(-5),
0031     EtaMax                   = cms.double(5),
0032     PtMin                    = cms.double(10),
0033     PtMax                    = cms.double(1000),
0034     #--- eta values for JEC vs pt plots ----
0035     VEta                     = cms.vdouble(0.0,1.0,2.0,3.0,4.0),
0036     #--- corrected pt values for JEC vs eta plots ----
0037     VPt                      = cms.vdouble(20,30,50,100,200),
0038     Debug                    = cms.untracked.bool(False),
0039     UseCondDB                = cms.untracked.bool(True)
0040 )
0041
0042 process.p = cms.Path(process.ak4PFL2L3ResidualCorrectorChain * process.ak4pfl2l3Residual)
0043