File indexing completed on 2024-04-06 12:25:28
0001
0002
0003
0004
0005
0006
0007
0008
0009 import FWCore.ParameterSet.Config as cms
0010
0011
0012
0013
0014
0015
0016
0017 isMC = True
0018
0019
0020
0021
0022
0023
0024
0025
0026
0027 NJetsToKeep = 2
0028
0029 inputFile = 'file:/uscms_data/d2/kalanand/dijet-Run2010A-JetMET-Nov4ReReco-9667events.root'
0030 if isMC:
0031 inputFile ='/store/mc/Fall10/QCD_Pt_80to120_TuneZ2_7TeV_pythia6/GEN-SIM-RECO/START38_V12-v1/0000/FEF4D100-4CCB-DF11-94CB-00E08178C12F.root'
0032
0033
0034
0035
0036
0037
0038
0039
0040 process = cms.Process("Ana")
0041
0042 process.load("FWCore.MessageService.MessageLogger_cfi")
0043 process.MessageLogger.cerr.FwkReport.reportEvery = 10
0044
0045
0046
0047
0048
0049
0050
0051
0052
0053
0054
0055 process.maxEvents = cms.untracked.PSet(
0056 input = cms.untracked.int32(100)
0057 )
0058
0059 process.source = cms.Source("PoolSource",
0060 fileNames = cms.untracked.vstring(inputFile)
0061 )
0062 process.source.inputCommands = cms.untracked.vstring("keep *","drop *_MEtoEDMConverter_*_*")
0063
0064
0065
0066
0067
0068
0069
0070
0071
0072 process.calo = cms.EDAnalyzer("CaloJetPlotsExample",
0073 JetAlgorithm = cms.string('ak4CaloJets'),
0074 HistoFileName = cms.string('CaloJetPlotsExample.root'),
0075 NJets = cms.int32(NJetsToKeep)
0076 )
0077
0078 process.pf = cms.EDAnalyzer("PFJetPlotsExample",
0079 JetAlgorithm = cms.string('ak4PFJets'),
0080 HistoFileName = cms.string('PFJetPlotsExample.root'),
0081 NJets = cms.int32(NJetsToKeep)
0082 )
0083
0084 process.jpt = cms.EDAnalyzer("JPTJetPlotsExample",
0085 JetAlgorithm = cms.string('JetPlusTrackZSPCorJetAntiKt4'),
0086 HistoFileName = cms.string('JPTJetPlotsExample.root'),
0087 NJets = cms.int32(NJetsToKeep)
0088 )
0089
0090 if isMC:
0091 process.gen = cms.EDAnalyzer("GenJetPlotsExample",
0092 JetAlgorithm = cms.string('ak4GenJets'),
0093 HistoFileName = cms.string('GenJetPlotsExample.root'),
0094 NJets = cms.int32(NJetsToKeep)
0095 )
0096
0097
0098
0099
0100
0101
0102
0103
0104
0105
0106 if isMC:
0107 process.p = cms.Path(process.calo*process.pf*process.jpt*process.gen)
0108 else:
0109 process.p = cms.Path(process.calo*process.pf*process.jpt)
0110
0111