TopTools/test/ttJetPartonMatch_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("TEST")
0004
0005 ## add message logger
0006 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0007 process.MessageLogger.cerr.threshold = 'INFO'
0008 process.MessageLogger.cerr.JetPartonMatching = cms.untracked.PSet(
0009     limit = cms.untracked.int32(-1)
0010 )
0011
0012 ## define input
0013 from TopQuarkAnalysis.TopEventProducers.tqafInputFiles_cff import relValTTbar
0014 process.source = cms.Source("PoolSource",
0015     fileNames = cms.untracked.vstring(relValTTbar)
0016 )
0017
0018 ## define maximal number of events to loop over
0019 process.maxEvents = cms.untracked.PSet(
0020     input = cms.untracked.int32(50)
0021 )
0022
0023 ## configure process options
0024 process.options = cms.untracked.PSet(
0025     wantSummary      = cms.untracked.bool(True)
0026 )
0027
0028 ## configure geometry & conditions
0029 process.load("Configuration.Geometry.GeometryIdeal_cff")
0030 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0031 from Configuration.AlCa.GlobalTag import GlobalTag
0032 process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:run2_mc')
0033 process.load("Configuration.StandardSequences.MagneticField_cff")
0034
0035 process.task = cms.Task()
0036
0037 ## std sequence for PAT
0038 process.load("PhysicsTools.PatAlgos.producersLayer1.patCandidates_cff")
0039 process.task.add(process.patCandidatesTask)
0040 process.load("PhysicsTools.PatAlgos.selectionLayer1.selectedPatCandidates_cff")
0041 process.task.add(process.selectedPatCandidatesTask)
0042
0043 ## std sequence to produce the ttGenEvt
0044 process.load("TopQuarkAnalysis.TopEventProducers.sequences.ttGenEvent_cff")
0045 process.task.add(process.makeGenEvtTask)
0046
0047 ## configure jet-parton matching
0048 process.load("TopQuarkAnalysis.TopTools.TtFullHadJetPartonMatch_cfi")
0049 process.task.add(process.ttFullHadJetPartonMatch)
0050 process.ttFullHadJetPartonMatch.verbosity  = 1              #default: 0
0051 #process.ttFullHadJetPartonMatch.algorithm  = "minSumDist"   #default: totalMinDist
0052 #process.ttFullHadJetPartonMatch.useDeltaR  = True           #default: True
0053 #process.ttFullHadJetPartonMatch.useMaxDist = True           #default: False
0054 #process.ttFullHadJetPartonMatch.maxDist    = 2.5            #default: 0.3
0055 #process.ttFullHadJetPartonMatch.maxNJets   = 7              #default: 6
0056 #process.ttFullHadJetPartonMatch.maxNComb   = 1              #default: 1
0057 process.load("TopQuarkAnalysis.TopTools.TtFullLepJetPartonMatch_cfi")
0058 process.task.add(process.ttFullLepJetPartonMatch)
0059 process.ttFullLepJetPartonMatch.verbosity  = 1              #default: 0
0060 #process.ttFullLepJetPartonMatch.algorithm  = "minSumDist"   #default: totalMinDist
0061 #process.ttFullLepJetPartonMatch.useDeltaR  = True           #default: True
0062 #process.ttFullLepJetPartonMatch.useMaxDist = True           #default: False
0063 #process.ttFullLepJetPartonMatch.maxDist    = 2.5            #default: 0.3
0064 #process.ttFullLepJetPartonMatch.maxNJets   = 3              #default: 2
0065 #process.ttFullLepJetPartonMatch.maxNComb   = 1              #default: 1
0066 process.load("TopQuarkAnalysis.TopTools.TtSemiLepJetPartonMatch_cfi")
0067 process.task.add(process.ttSemiLepJetPartonMatch)
0068 process.ttSemiLepJetPartonMatch.verbosity  = 1              #default: 0
0069 #process.ttSemiLepJetPartonMatch.algorithm  = "minSumDist"   #default: totalMinDist
0070 #process.ttSemiLepJetPartonMatch.useDeltaR  = True           #default: True
0071 #process.ttSemiLepJetPartonMatch.useMaxDist = True           #default: False
0072 #process.ttSemiLepJetPartonMatch.maxDist    = 2.5            #default: 0.3
0073 #process.ttSemiLepJetPartonMatch.maxNJets   = 5              #default: 4
0074 #process.ttSemiLepJetPartonMatch.maxNComb   = 1              #default: 1
0075
0076 ## configure output module
0077 process.out = cms.OutputModule("PoolOutputModule",
0078     fileName = cms.untracked.string('ttJetPartonMatch.root'),
0079     outputCommands = cms.untracked.vstring('drop *')
0080 )
0081 process.out.outputCommands += ['keep *_ttFullHadJetPartonMatch_*_*',
0082                                'keep *_ttFullLepJetPartonMatch_*_*',
0083                                'keep *_ttSemiLepJetPartonMatch_*_*']
0084
0085 ## output path
0086 process.outpath = cms.EndPath(process.out, process.task)