Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2024-04-06 12:31:14

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.TtFullHadronicEvent = 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.GeometryRecoDB_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 #Temporary customize to the unit tests that fail due to old input samples
0041 process.patTaus.skipMissingTauID = True
0042 process.load("PhysicsTools.PatAlgos.selectionLayer1.selectedPatCandidates_cff")
0043 process.task.add(process.selectedPatCandidatesTask)
0044 
0045 ## std sequence to produce the ttGenEvt
0046 process.load("TopQuarkAnalysis.TopEventProducers.sequences.ttGenEvent_cff")
0047 process.task.add(process.makeGenEvtTask)
0048 
0049 ## std sequence to produce the ttFullHadEvent
0050 process.load("TopQuarkAnalysis.TopEventProducers.sequences.ttFullHadEvtBuilder_cff")
0051 process.task.add(process.makeTtFullHadEventTask)
0052 process.ttFullHadEvent.verbosity = 1
0053 
0054 ## choose which hypotheses to produce
0055 from TopQuarkAnalysis.TopEventProducers.sequences.ttFullHadEvtBuilder_cff import *
0056 addTtFullHadHypotheses(process,
0057                        ["kKinFit"]
0058                        )
0059 
0060 #removeTtFullHadHypGenMatch(process)
0061 
0062 ## change maximum number of jets taken into account per event (default: 6)
0063 #from TopQuarkAnalysis.TopEventProducers.sequences.ttFullHadEvtBuilder_cff import *
0064 #setForAllTtFullHadHypotheses(process, "maxNJets", 8)
0065 
0066 ## configure output module
0067 process.out = cms.OutputModule("PoolOutputModule",
0068     fileName     = cms.untracked.string('ttFullHadEvtBuilder.root'),
0069     outputCommands = cms.untracked.vstring('drop *'),
0070     dropMetaData = cms.untracked.string('DROPPED')
0071 )
0072 process.outpath = cms.EndPath(process.out, process.task)
0073 
0074 ## PAT content
0075 from PhysicsTools.PatAlgos.patEventContent_cff import patEventContentNoCleaning
0076 process.out.outputCommands += patEventContentNoCleaning
0077 ## TQAF content
0078 from TopQuarkAnalysis.TopEventProducers.tqafEventContent_cff import tqafEventContent
0079 process.out.outputCommands += tqafEventContent