JetExamples/test/printGenJetRatio.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("testJET")
0004
0005 process.options = cms.untracked.PSet(
0006     wantSummary = cms.untracked.bool(True)
0007 )
0008
0009 process.load('FWCore.MessageService.MessageLogger_cfi')
0010 process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
0011 process.load('IOMC.RandomEngine.IOMC_cff')
0012 process.load('Configuration.StandardSequences.Generator_cff')
0013 process.load('IOMC.EventVertexGenerators.VtxSmearedRealistic8TeVCollision_cfi')
0014
0015 process.maxEvents = cms.untracked.PSet(
0016     input = cms.untracked.int32(5)
0017 )
0018
0019 # Input source
0020 process.source = cms.Source("EmptySource")
0021
0022 from Configuration.Generator.PythiaUEZ2starSettings_cfi import *
0023 process.generator = cms.EDFilter("Pythia6GeneratorFilter",
0024     pythiaPylistVerbosity = cms.untracked.int32(0),
0025     filterEfficiency = cms.untracked.double(1.0),
0026     pythiaHepMCVerbosity = cms.untracked.bool(False),
0027     comEnergy = cms.double(8000.0),
0028     maxEventsToPrint = cms.untracked.int32(0),
0029     PythiaParameters = cms.PSet(
0030         pythiaUESettingsBlock,
0031         processParameters = cms.vstring(
0032             'MSEL=5               ! high pT process -> bb',
0033             'CKIN(3)=170.         ! minimum pt hat for hard interactions',
0034             'CKIN(4)=230.         ! maximum pt hat for hard interactions'
0035         ),
0036         # This is a vector of ParameterSet names to be read, in this order
0037         parameterSets = cms.vstring('pythiaUESettings','processParameters')
0038     )
0039 )
0040
0041 process.printList = cms.EDAnalyzer("ParticleListDrawer",
0042     src = cms.InputTag("genParticles"),
0043     maxEventsToPrint = cms.untracked.int32(1)
0044 )
0045
0046 process.BCRatio = cms.EDProducer("GenJetBCEnergyRatio",
0047     genJets = cms.InputTag("ak5GenJets")
0048 )
0049
0050 process.printEvent = cms.EDAnalyzer("printGenJetRatio",
0051     srcBratio = cms.InputTag("BCRatio","bRatioCollection"),
0052     srcCratio = cms.InputTag("BCRatio","cRatioCollection")
0053 )
0054
0055 process.p = cms.Path(process.generator*process.pgen*process.printList*process.BCRatio*process.printEvent)
0056
0057 #process.MessageLogger.cout = cms.PSet(
0058 #    enable = cms.untracked.bool(True),
0059 #    threshold = cms.untracked.string('ERROR')
0060 #)