Physics/test/PFJet_topDQM_production_cfg.py

0001 ## THIS IS CURRENTLY COMPATIBLE ONLY WITH SINGLE TOP MODULES
0002 ## since the b-tagging algorithms are here re-run with PFJets as input
0003
0004 import FWCore.ParameterSet.Config as cms
0005
0006 process = cms.Process('TOPDQM')
0007
0008 ## imports of standard configurations
0009 process.load('DQMOffline.Configuration.DQMOffline_cff')
0010 process.load('Configuration.StandardSequences.EndOfProcess_cff')
0011 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0012
0013
0014
0015 ## --------------------------------------------------------------------
0016 ## Frontier Conditions: (adjust accordingly!!!)
0017 ##
0018 ## For CMSSW_3_8_X MC use             ---> 'START38_V12::All'
0019 ## For Data (38X re-processing) use   ---> 'GR_R_38X_V13::All'
0020 ## For Data (38X prompt reco) use     ---> 'GR10_P_V10::All'
0021 ##
0022 ## For more details have a look at: WGuideFrontierConditions
0023 ## --------------------------------------------------------------------
0024 ##process.GlobalTag.globaltag = 'GR_R_42_V14::All'
0025 #process.GlobalTag.globaltag = 'GR_R_53_V14::All'
0026 process.GlobalTag.globaltag = 'GR_R_52_V7::All'
0027
0028
0029 ## input file(s) for testing
0030 process.source = cms.Source("PoolSource",
0031     fileNames = cms.untracked.vstring(
0032      #"/store/data/Run2012A/SingleMu/RECO/PromptReco-v1/000/193/116/0AB80D76-FA95-E111-8C46-5404A63886B9.root",
0033      #"/store/data/Run2012A/SingleMu/RECO/PromptReco-v1/000/193/116/0EDA5C6F-FA95-E111-8681-002481E0E56C.root",
0034      #'/store/data/Run2012A/SingleElectron/RECO/PromptReco-v1/000/190/456/1412AF9C-0681-E111-AF6F-003048D2BBF0.root'
0035      #'/store/relval/CMSSW_5_2_3_patch3/RelValTTbarLepton/GEN-SIM-RECO/START52_V9_special_120410-v1/0122/2C3473C4-1583-E111-8CE8-002618943870.root'
0036     "/store/data/Run2012B/SingleMu/AOD/PromptReco-v1/000/194/702/A0D074D6-E7A5-E111-A2B7-BCAEC518FF41.root"
0037      )
0038 )
0039
0040 ## number of events
0041 process.maxEvents = cms.untracked.PSet(
0042   input = cms.untracked.int32(-1)
0043 )
0044
0045 ## apply VBTF electronID (needed for the current implementation
0046 ## of topSingleElectronDQMLoose and topSingleElectronDQMMedium)
0047 process.load("Configuration.StandardSequences.GeometryRecoDB_cff")
0048 process.load("Configuration.StandardSequences.MagneticField_cff")
0049 process.load("Configuration.StandardSequences.Reconstruction_cff")
0050 process.load("DQM.Physics.topElectronID_cff")
0051
0052 # b-tagging
0053 # Re-run b-tagging with PFJets as input
0054 # b-tagging
0055 # Re-run b-tagging with PFJets as input
0056 process.load('Configuration/StandardSequences/Reconstruction_cff')
0057 process.ak4JetTracksAssociatorAtVertex.jets = cms.InputTag("ak4PFJets")
0058 process.p = cms.Path(process.ak4JetTracksAssociatorAtVertex*process.btagging)
0059
0060
0061
0062 #process.topSingleMuonLooseDQM.setup.triggerExtras.src  = cms.InputTag("TriggerResults","","REDIGI42X")
0063 #process.topSingleMuonLooseDQM.preselection.trigger.src = cms.InputTag("TriggerResults","","REDIGI42X")
0064 #process.topSingleMuonLooseDQM.preselection.trigger.select  = cms.vstring(['HLT_Mu15_v2'])
0065 #process.topSingleMuonMediumDQM.preselection.trigger.select = cms.vstring(['HLT_Mu15_v2'])
0066
0067
0068
0069 ## output
0070 process.output = cms.OutputModule("PoolOutputModule",
0071   fileName       = cms.untracked.string('topDQM_production.root'),
0072   outputCommands = cms.untracked.vstring(
0073     'drop *_*_*_*',
0074     'keep *_*_*_TOPDQM',
0075     'drop *_TriggerResults_*_TOPDQM',
0076     'drop *_simpleEleId70cIso_*_TOPDQM',
0077     'drop *_ak4JetTracksAssociatorAtVertex_*_TOPDQM',
0078     'drop *_btagging_*_TOPDQM',
0079     'drop *_jetProbabilityBJetTags_*_TOPDQM',
0080     'drop *_ghostTrackBJetTags_*_TOPDQM',
0081     'drop *_combinedSecondaryVertexMVABJetTags_*_TOPDQM',
0082     'drop *_trackCountingHighPurBJetTags_*_TOPDQM',
0083     'drop *_trackCountingHighEffBJetTags_*_TOPDQM',
0084     'drop *_simpleSecondaryVertexHighEffBJetTags_*_TOPDQM',
0085     'drop *_simpleSecondaryVertexHighPurBJetTags_*_TOPDQM',
0086     'drop *_softElectronByIP3dBJetTags_*_TOPDQM',
0087     'drop *_softElectronByPtBJetTags_*_TOPDQM',
0088     'drop *_softMuonBJetTags_*_TOPDQM',
0089     'drop *_softMuonByIP3dBJetTags_*_TOPDQM',
0090     'drop *_impactParameterTagInfos_*_TOPDQM',
0091     'drop *_combinedSecondaryVertexBJetTags_*_TOPDQM',
0092     'drop *_softMuonByPtBJetTags_*_TOPDQM',
0093     'drop *_ghostTrackVertexTagInfos_*_TOPDQM',
0094     'drop *_secondaryVertexTagInfos_*_TOPDQM',
0095     'drop *_softElectronCands_*_TOPDQM',
0096     'drop *_softMuonTagInfos_*_TOPDQM',
0097     'drop *_softElectronTagInfos_*_TOPDQM',
0098     ),
0099   splitLevel     = cms.untracked.int32(0),
0100   dataset = cms.untracked.PSet(
0101     dataTier   = cms.untracked.string(''),
0102     filterName = cms.untracked.string('')
0103   )
0104 )
0105
0106 ## load jet corrections
0107 process.load("JetMETCorrections.Configuration.JetCorrectionServicesAllAlgos_cff")
0108 process.load("JetMETCorrections.Configuration.DefaultJEC_cff")
0109 #process.prefer("ak4CaloL2L3")
0110 process.prefer("ak4PFL2L3")
0111
0112
0113 ## check the event content
0114 process.content = cms.EDAnalyzer("EventContentAnalyzer")
0115
0116 ## configure message logger
0117 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0118 process.MessageLogger.cerr.threshold = 'INFO'
0119 process.MessageLogger.cerr.TopSingleLeptonDQM    = cms.untracked.PSet(limit = cms.untracked.int32(1))
0120 process.MessageLogger.cerr.SingleTopTChannelLeptonDQM    = cms.untracked.PSet(limit = cms.untracked.int32(1))
0121 process.MessageLogger.cerr.TopDiLeptonOfflineDQM = cms.untracked.PSet(limit = cms.untracked.int32(1))
0122
0123 process.MessageLogger.cerr.FwkReport.reportEvery = 1
0124
0125
0126
0127 ## path definitions
0128 process.p      = cms.Path(
0129     process.simpleEleId70cIso          *
0130     process.singleTopMuonMediumDQM     +
0131     process.singleTopElectronMediumDQM
0132
0133 )
0134 process.endjob = cms.Path(
0135     process.endOfProcess
0136 )
0137 process.fanout = cms.EndPath(
0138     process.output
0139 )
0140
0141 ## schedule definition
0142 process.schedule = cms.Schedule(
0143     process.p,
0144     process.endjob,
0145     process.fanout
0146 )