Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2024-04-06 12:19:16

0001 import FWCore.ParameterSet.Config as cms
0002 #!
0003 #! PROCESS
0004 #!
0005 process = cms.Process("JEC")
0006 
0007 #!
0008 #! INPUT
0009 #!
0010 process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(1000) )
0011 process.source = cms.Source(
0012     'PoolSource',
0013     fileNames = cms.untracked.vstring('/store/relval/CMSSW_7_2_0_pre7/RelValProdTTbar/GEN-SIM-RECO/PRE_STA72_V4-v1/00000/B223AEC2-B94B-E411-884B-00261894395F.root')
0014     )
0015 
0016 #!
0017 #! SERVICES
0018 #!
0019 
0020 
0021 #!
0022 #! JET CORRECTION
0023 #!
0024 
0025 
0026 
0027 process.load('Configuration.StandardSequences.Services_cff')
0028 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0029 from Configuration.AlCa.GlobalTag import GlobalTag
0030 process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:startup')
0031 process.load('JetMETCorrections.Configuration.CorrectedJetProducersDefault_cff')
0032 process.load('JetMETCorrections.Configuration.CorrectedJetProducersAllAlgos_cff')
0033 
0034 
0035 process.TFileService=cms.Service("TFileService",fileName=cms.string('histos.root'))
0036 
0037 #!
0038 #! MAKE SOME HISTOGRAMS
0039 #!
0040 jetPtHistogram = cms.PSet(min          = cms.untracked.double(     50),
0041                           max          = cms.untracked.double(    500),
0042                           nbins        = cms.untracked.int32 (     50),
0043                           name         = cms.untracked.string('JetPt'),
0044                           description  = cms.untracked.string(     ''),
0045                           plotquantity = cms.untracked.string(   'pt')
0046                           )
0047 
0048 process.ak4CaloL2L3Histos = cms.EDAnalyzer(
0049     'CandViewHistoAnalyzer',
0050     src = cms.InputTag('ak4CaloJetsL2L3'),
0051     histograms = cms.VPSet(jetPtHistogram)
0052     )
0053 process.ak7CaloL2L3Histos = process.ak4CaloL2L3Histos.clone(src = 'ak7CaloJetsL2L3')
0054 process.kt4CaloL2L3Histos = process.ak4CaloL2L3Histos.clone(src = 'kt4CaloJetsL2L3')
0055 #process.kt6CaloL2L3Histos = process.ak4CaloL2L3Histos.clone(src = 'kt6CaloJetsL2L3')
0056 
0057 process.ak4PFL2L3Histos = process.ak4CaloL2L3Histos.clone(src = 'ak4PFJetsL2L3')
0058 process.ak4PFCHSL2L3Histos = process.ak4CaloL2L3Histos.clone(src = 'ak4PFCHSJetsL2L3')
0059 process.ak8PFL2L3Histos = process.ak4CaloL2L3Histos.clone(src = 'ak8PFJetsL2L3')
0060 process.kt4PFL2L3Histos = process.ak4CaloL2L3Histos.clone(src = 'kt4PFJetsL2L3')
0061 #process.kt6PFL2L3Histos = process.ak4CaloL2L3Histos.clone(src = 'kt6PFJetsL2L3')
0062 
0063 process.ak4JPTL2L3Histos = process.ak4CaloL2L3Histos.clone(src = 'ak4JPTJetsL2L3')
0064 process.ak4TrackL2L3Histos = process.ak4CaloL2L3Histos.clone(src = 'ak4TrackJetsL2L3')
0065 
0066 #
0067 # RUN!
0068 #
0069 process.run = cms.Path(
0070 #------ create the corrected calojet collection and run the histogram module ------
0071 #process.ak4CaloL2L3CorrectorChain * process.ak4CaloJetsL2L3 * process.ak4CaloL2L3Histos *
0072 #process.ak7CaloL2L3CorrectorChain * process.ak7CaloJetsL2L3 * process.ak7CaloL2L3Histos *
0073 #process.kt4CaloL2L3CorrectorChain * process.kt4CaloJetsL2L3 * process.kt4CaloL2L3Histos *
0074 #process.kt6CaloL2L3CorrectorChain * process.kt6CaloJetsL2L3 * process.kt6CaloL2L3Histos *
0075 #------ create the corrected pfjet collection and run the histogram module --------
0076 process.ak4PFL2L3CorrectorChain * process.ak4PFJetsL2L3 * process.ak4PFL2L3Histos *
0077 process.ak4PFCHSL2L3CorrectorChain * process.ak4PFCHSJetsL2L3 * process.ak4PFCHSL2L3Histos
0078 #process.ak8PFL2L3CorrectorChain * process.ak8PFJetsL2L3 * process.ak8PFL2L3Histos *
0079 #process.kt4PFL2L3CorrectorChain * process.kt4PFJetsL2L3 * process.kt4PFL2L3Histos *
0080 #process.kt6PFL2L3CorrectorChain * process.kt6PFJetsL2L3 * process.kt6PFL2L3Histos *
0081 #------ create the corrected jptjet collection and run the histogram module -------
0082 #process.ak4JPTL2L3CorrectorChain * process.ak4JPTJetsL2L3 * process.ak4JPTL2L3Histos *
0083 #------ create the corrected trackjet collection and run the histogram module -----
0084 #process.ak4TrackL2L3CorrectorChain * process.ak4TrackJetsL2L3 * process.ak4TrackL2L3Histos
0085 )
0086 
0087 
0088 
0089 
0090 
0091 
0092 
0093