JetProducers/test/runJetReco_Substructure_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("JETRECO")
0004
0005 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0006 process.options   = cms.untracked.PSet( wantSummary = cms.untracked.bool(True) )
0007
0008 # input
0009 process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(1000))
0010 process.source = cms.Source(
0011     "PoolSource",
0012     fileNames = cms.untracked.vstring(
0013         '/store/relval/CMSSW_5_0_0/RelValTTbar/GEN-SIM-RECO/START50_V8-v3/0074/182BA6EB-BD2A-E111-801C-003048678A88.root',
0014         '/store/relval/CMSSW_5_0_0/RelValTTbar/GEN-SIM-RECO/START50_V8-v3/0074/8E79DF51-B92A-E111-A9B4-00261894396E.root',
0015         '/store/relval/CMSSW_5_0_0/RelValTTbar/GEN-SIM-RECO/START50_V8-v3/0074/AE46DE36-C02A-E111-AF78-003048FF9AC6.root',
0016         '/store/relval/CMSSW_5_0_0/RelValTTbar/GEN-SIM-RECO/START50_V8-v3/0082/5893EB84-542B-E111-AFDF-002618943915.root',
0017         '/store/relval/CMSSW_5_0_0/RelValTTbar/GEN-SIM-RECO/START50_V8-v3/0082/8C111086-542B-E111-AA93-002618943983.root',
0018         '/store/relval/CMSSW_5_0_0/RelValTTbar/GEN-SIM-RECO/START50_V8-v3/0082/E0196E85-542B-E111-AA08-0030486791DC.root'
0019     )
0020     #inputCommands = cms.untracked.vstring('keep *_*_*_*','drop recoTrackExtrapolations_*_*_RECO')
0021     )
0022 # output
0023 process.load('Configuration/EventContent/EventContent_cff')
0024 process.output = cms.OutputModule(
0025     "PoolOutputModule",
0026     outputCommands = process.RecoJetsAOD.outputCommands,
0027     fileName = cms.untracked.string('testJetRecoRECO.root'),
0028     dataset = cms.untracked.PSet(
0029             dataTier = cms.untracked.string(''),
0030                     filterName = cms.untracked.string('')
0031                     )
0032 )
0033 process.output.outputCommands.append('drop *_*_*_RECO')
0034 #process.output.outputCommands.append('keep recoCaloJets_*_*_*')
0035 #process.output.outputCommands.append('keep recoPFJets_*_*_*')
0036 #process.output.outputCommands.append('keep recoGenJets_*_*_*')
0037 #process.output.outputCommands.append('keep recoBasicJets_*_*_*')
0038 process.output.outputCommands.append('keep *_*_*_JETRECO')
0039 process.output.outputCommands.append('keep recoPFCandidates_particleFlow_*_*')
0040 #process.output.outputCommands.append('keep *_trackRefsForJets_*_*')
0041 #process.output.outputCommands.append('keep *_generalTracks_*_*')
0042
0043 # jet reconstruction
0044 process.load('Configuration.StandardSequences.Services_cff')
0045 process.load('Configuration.StandardSequences.GeometryDB_cff')
0046 process.load('Configuration.StandardSequences.MagneticField_38T_cff')
0047 process.load('Configuration.StandardSequences.Reconstruction_cff')
0048 process.load('Configuration/StandardSequences/FrontierConditions_GlobalTag_cff')
0049 process.GlobalTag.globaltag = 'START50_V10::All'
0050
0051 process.load("RecoJets/Configuration/RecoPFClusterJets_cff")
0052 process.load("RecoMET/METProducers/PFClusterMET_cfi")
0053
0054 process.load("RecoJets/JetAssociationProducers/trackExtrapolator_cfi")
0055
0056 process.ca12PFJetsMassDropFiltered = process.ak5PFJetsMassDropFiltered.clone(
0057     rParam = cms.double(1.2)
0058     )
0059
0060 process.ak5PFJetsTrimmed.doAreaFastjet = True
0061 process.ak5PFJetsPruned.doAreaFastjet = True
0062 process.ca8PFJetsPruned = process.ak5PFJetsPruned.clone(
0063     jetAlgorithm = cms.string("CambridgeAachen"),
0064     rParam       = cms.double(0.8)
0065     )
0066
0067 #process.recoJets = cms.Path(process.trackExtrapolator+process.jetGlobalReco+process.CastorFullReco+process.jetHighLevelReco+process.recoPFClusterJets)
0068 process.recoJets = cms.Path(process.ak5PFJets+
0069                             process.ak5PFJetsTrimmed+
0070                             process.ak5PFJetsFiltered+
0071                             process.ak5PFJetsMassDropFiltered+
0072                             process.ca12PFJetsMassDropFiltered+
0073                             process.ak5PFJetsPruned+
0074                             process.ca8PFJetsPruned
0075                             )
0076
0077
0078
0079
0080
0081 process.out = cms.EndPath(process.output)
0082
0083 # Set the threshold for output logging to 'info'
0084 #process.MessageLogger.cerr.threshold = cms.untracked.string('INFO')
0085 process.MessageLogger.cerr.threshold = cms.untracked.string('DEBUG')
0086 #process.MessageLogger.cout.threshold = cms.untracked.string('DEBUG')
0087 process.MessageLogger.debugModules = cms.untracked.vstring('*')
0088 process.MessageLogger.cerr.INFO = cms.untracked.PSet(
0089              limit = cms.untracked.int32(0)
0090 )
0091
0092 # schedule
0093 process.schedule = cms.Schedule(process.recoJets,process.out)