TopJetCombination/test/ttSemiLepJetCombMVAComputer_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
0009 ## define input
0010 from TopQuarkAnalysis.TopEventProducers.tqafInputFiles_cff import relValTTbar
0011 process.source = cms.Source("PoolSource",
0012     fileNames = cms.untracked.vstring(relValTTbar)
0013 )
0014
0015 ## define maximal number of events to loop over
0016 process.maxEvents = cms.untracked.PSet(
0017     input = cms.untracked.int32(100)
0018 )
0019
0020 ## configure process options
0021 process.options = cms.untracked.PSet(
0022     wantSummary      = cms.untracked.bool(True)
0023 )
0024
0025 ## configure geometry & conditions
0026 process.load("Configuration.Geometry.GeometryRecoDB_cff")
0027 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0028 from Configuration.AlCa.GlobalTag import GlobalTag
0029 process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:run2_mc')
0030 process.load("Configuration.StandardSequences.MagneticField_cff")
0031
0032 process.task = cms.Task()
0033
0034 ## std sequence for pat
0035 process.load("PhysicsTools.PatAlgos.producersLayer1.patCandidates_cff")
0036 process.task.add(process.patCandidatesTask)
0037 process.load("PhysicsTools.PatAlgos.selectionLayer1.selectedPatCandidates_cff")
0038 process.task.add(process.selectedPatCandidatesTask)
0039
0040 ## configure mva computer
0041 process.load("TopQuarkAnalysis.TopJetCombination.TtSemiLepJetCombMVAComputer_cff")
0042 process.task.add(process.findTtSemiLepJetCombMVA)
0043 ## change maximum number of jets taken into account per event (default: 4)
0044 #process.findTtSemiLepJetCombMVA.maxNJets = 5
0045
0046 ## configure output module
0047 process.out = cms.OutputModule("PoolOutputModule",
0048     fileName = cms.untracked.string('ttSemiLepJetCombMVAComputer_muons.root'),
0049     outputCommands = cms.untracked.vstring('drop *')
0050 )
0051 process.out.outputCommands += ['keep *_findTtSemiLepJetCombMVA_*_*']
0052
0053 ## output path
0054 process.outpath = cms.EndPath(process.out, process.task)