JetExamples/test/printJetFlavour.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
0011 process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
0012
0013 process.maxEvents = cms.untracked.PSet(
0014     input = cms.untracked.int32(5)
0015 )
0016
0017 process.source = cms.Source("PoolSource",
0018     fileNames = cms.untracked.vstring(
0019         # /TTJets_MassiveBinDECAY_TuneZ2star_8TeV-madgraph-tauola/Summer12_DR53X-PU_S10_START53_V7C-v1/AODSIM
0020         '/store/mc/Summer12_DR53X/TTJets_MassiveBinDECAY_TuneZ2star_8TeV-madgraph-tauola/AODSIM/PU_S10_START53_V7C-v1/00000/008BD264-1526-E211-897A-00266CFFA7BC.root'
0021     )
0022 )
0023
0024 process.printList = cms.EDAnalyzer("ParticleListDrawer",
0025     src = cms.InputTag("genParticles"),
0026     maxEventsToPrint  = cms.untracked.int32(1)
0027 )
0028
0029 process.myPartons = cms.EDProducer("PartonSelector",
0030     src = cms.InputTag("genParticles"),
0031     withLeptons = cms.bool(False)
0032 )
0033
0034 process.flavourByRef = cms.EDProducer("JetPartonMatcher",
0035     jets = cms.InputTag("ak5PFJets"),
0036     coneSizeToAssociate = cms.double(0.3),
0037     partons = cms.InputTag("myPartons")
0038 )
0039
0040 process.flavourByVal = cms.EDProducer("JetFlavourIdentifier",
0041     srcByReference = cms.InputTag("flavourByRef"),
0042     physicsDefinition = cms.bool(False)
0043 )
0044
0045 process.printEvent = cms.EDAnalyzer("printJetFlavour",
0046     srcSelectedPartons = cms.InputTag("myPartons"),
0047     srcByReference = cms.InputTag("flavourByRef"),
0048     srcByValue = cms.InputTag("flavourByVal")
0049 )
0050
0051 process.p = cms.Path(process.printList*process.myPartons*process.flavourByRef*process.flavourByVal*process.printEvent)
0052
0053 #process.MessageLogger.cout = dict(enable = True, threshold = 'ERROR')
0054
0055