Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-02-14 14:21:38

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 #process.GlobalTag.globaltag = 'START38_V13::All'
0009 
0010 
0011 ##-------------------- Import the JEC services -----------------------
0012 process.load('JetMETCorrections.Configuration.DefaultJEC_cff')
0013 
0014 
0015 from CondCore.DBCommon.CondDBSetup_cfi import *
0016 process.jec = cms.ESSource("PoolDBESSource",CondDBSetup,
0017                    connect = cms.string("sqlite:JEC_Spring10.db"),
0018                    toGet =  cms.VPSet(
0019                        cms.PSet(record = cms.string("JetCorrectionsRecord"),
0020                                 tag = cms.string("JetCorrectorParametersCollection_Spring10_V6_AK5Calo"),
0021                                 label=cms.untracked.string("AK5Calo")),
0022                        cms.PSet(record = cms.string("JetCorrectionsRecord"),
0023                                 tag = cms.string("JetCorrectorParametersCollection_Spring10_V6_AK5PF"),
0024                                 label=cms.untracked.string("AK5PF")),                                   
0025                        cms.PSet(record = cms.string("JetCorrectionsRecord"),
0026                                 tag = cms.string("JetCorrectorParametersCollection_Summer10_V6_AK5JPT"),
0027                                 label=cms.untracked.string("AK5JPT")),
0028                        cms.PSet(record = cms.string("JetCorrectionsRecord"),
0029                                 tag = cms.string("JetCorrectorParametersCollection_Spring10_V6_IC5Calo"),
0030                                 label=cms.untracked.string("IC5Calo")),
0031                        cms.PSet(record = cms.string("JetCorrectionsRecord"),
0032                                 tag = cms.string("JetCorrectorParametersCollection_Spring10_V6_KT6Calo"),
0033                                 label=cms.untracked.string("KT6Calo")),
0034                        cms.PSet(record = cms.string("JetCorrectionsRecord"),
0035                                 tag = cms.string("JetCorrectorParametersCollection_Spring10_V6_KT4Calo"),
0036                                 label=cms.untracked.string("KT4Calo")),
0037                        cms.PSet(record = cms.string("JetCorrectionsRecord"),
0038                                 tag = cms.string("JetCorrectorParametersCollection_Spring10_V6_AK7Calo"),
0039                                 label=cms.untracked.string("AK7Calo")),
0040                        )
0041                    
0042                    )
0043 es_prefer_jec = cms.ESPrefer("PoolDBESSource","jec")
0044 
0045 
0046 
0047 
0048 
0049 
0050 ##-------------------- Define the source  ----------------------------
0051 process.maxEvents = cms.untracked.PSet(
0052         input = cms.untracked.int32(1)
0053         )
0054 process.source = cms.Source("EmptySource")
0055 
0056 ##-------------------- User analyzer  --------------------------------
0057 process.ak5caloDemo  = cms.EDAnalyzer('FactorizedJetCorrectorDemo',
0058     levels                   = cms.vstring( 'L2Relative', 'L3Absolute', 'L5Flavor_gJ', 'L7Parton_gJ'),
0059     UncertaintyTag           = cms.string('Uncertainty'),
0060     UncertaintyFile          = cms.string(''),
0061     PayloadName              = cms.string('AK5Calo'),
0062     NHistoPoints             = cms.int32(10000),
0063     NGraphPoints             = cms.int32(500),
0064     EtaMin                   = cms.double(-5),
0065     EtaMax                   = cms.double(5),
0066     PtMin                    = cms.double(10),
0067     PtMax                    = cms.double(1000),
0068     #--- eta values for JEC vs pt plots ----
0069     VEta                     = cms.vdouble(0.0,1.0,2.0,3.0,4.0),
0070     #--- corrected pt values for JEC vs eta plots ----
0071     VPt                      = cms.vdouble(20,30,50,100,200),
0072     Debug                    = cms.untracked.bool(True)
0073 )
0074 
0075 process.p = cms.Path(process.ak5caloDemo)
0076