Skimming/test/SDmaker_6SD_3CS_PDMinBias_1e28_reprocess370patch2_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003
0004 process = cms.Process("makeSD")
0005
0006 process.configurationMetadata = cms.untracked.PSet(
0007     version = cms.untracked.string('$Revision: 1.1 $'),
0008     annotation = cms.untracked.string('SD and central skims'),
0009     name = cms.untracked.string('$Source: /cvs_server/repositories/CMSSW/CMSSW/Configuration/Skimming/test/SDmaker_6SD_3CS_PDMinBias_1e28_reprocess370_cfg.py,v $')
0010 )
0011
0012
0013 process.maxEvents = cms.untracked.PSet(
0014     input = cms.untracked.int32(1000)
0015 )
0016
0017 process.load("Configuration.StandardSequences.MagneticField_38T_cff")
0018 process.load("Configuration.StandardSequences.GeometryExtended_cff")
0019 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0020 process.load('Configuration.EventContent.EventContent_cff')
0021 process.GlobalTag.globaltag = "GR_R_37X_V6A::All"
0022
0023
0024 process.source = cms.Source("PoolSource",
0025     fileNames = cms.untracked.vstring(
0026         '/store/data/Commissioning10/MinimumBias/RAW-RECO/v8/000/132/601/F85204EE-EB40-DF11-8F71-001A64789D1C.root'
0027         ),
0028                             secondaryFileNames = cms.untracked.vstring(
0029         '/store/data/Commissioning09/Cosmics/RAW/v3/000/105/755/F6887FD0-9371-DE11-B69E-00304879FBB2.root'
0030         )
0031 )
0032 process.source.inputCommands = cms.untracked.vstring("keep *", "drop *_MEtoEDMConverter_*_*")
0033
0034 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0035 process.MessageLogger.cerr.FwkReport.reportEvery = 100
0036
0037 import HLTrigger.HLTfilters.hltHighLevelDev_cfi
0038
0039
0040 ### JetMETTau SD
0041 process.JetMETTau_1e28 = HLTrigger.HLTfilters.hltHighLevelDev_cfi.hltHighLevelDev.clone(andOr = True)
0042 process.JetMETTau_1e28.HLTPaths = (
0043 "HLT_Jet15U",
0044 "HLT_DiJetAve15U_8E29",
0045 "HLT_FwdJet20U",
0046 "HLT_Jet30U",
0047 "HLT_Jet50U",
0048 "HLT_DiJetAve30U_8E29",
0049 "HLT_QuadJet15U",
0050 "HLT_MET45",
0051 "HLT_MET100",
0052 "HLT_HT100U",
0053 "HLT_SingleLooseIsoTau20",
0054 "HLT_DoubleLooseIsoTau15",
0055 "HLT_DoubleJet15U_ForwardBackward",
0056 "HLT_BTagMu_Jet10U",
0057 "HLT_BTagIP_Jet50U",
0058 "HLT_StoppedHSCP_8E29"
0059 )
0060 process.JetMETTau_1e28.HLTPathsPrescales  = cms.vuint32(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1)
0061 process.JetMETTau_1e28.HLTOverallPrescale = cms.uint32(1)
0062 process.JetMETTau_1e28.throw = False
0063 process.JetMETTau_1e28.andOr = True
0064
0065 process.filterSdJetMETTau_1e28 = cms.Path(process.JetMETTau_1e28)
0066
0067
0068
0069 ### JetMETTauMonitor SD
0070 process.JetMETTauMonitor_1e28 = HLTrigger.HLTfilters.hltHighLevelDev_cfi.hltHighLevelDev.clone(andOr = True)
0071 process.JetMETTauMonitor_1e28.HLTPaths = (
0072 "HLT_L1Jet6U",
0073 "HLT_L1MET20",
0074 "HLT_L1SingleCenJet",
0075 "HLT_L1SingleForJet",
0076 "HLT_L1SingleTauJet",
0077 "HLT_L1Jet10U",
0078 "HLT_L1Jet10U_NoBPTX",
0079 "HLT_L1Jet6U_NoBPTX",
0080 "HLT_L1SingleCenJet_NoBPTX",
0081 "HLT_L1SingleForJet_NoBPTX",
0082 "HLT_L1SingleTauJet_NoBPTX"
0083 )
0084 process.JetMETTauMonitor_1e28.HLTPathsPrescales  = cms.vuint32(1,1,1,1,1,1,1,1,1,1,1)
0085 process.JetMETTauMonitor_1e28.HLTOverallPrescale = cms.uint32(1)
0086 process.JetMETTauMonitor_1e28.throw = False
0087 process.JetMETTauMonitor_1e28.andOr = True
0088
0089 process.filterSdJetMETTauMonitor_1e28 = cms.Path(process.JetMETTauMonitor_1e28)
0090
0091 ### MuMonitor SD
0092 process.MuMonitor_1e28 = HLTrigger.HLTfilters.hltHighLevelDev_cfi.hltHighLevelDev.clone(andOr = True)
0093 process.MuMonitor_1e28.HLTPaths = (
0094 "HLT_L1MuOpen",
0095 "HLT_L1Mu"
0096 )
0097 process.MuMonitor_1e28.HLTPathsPrescales  = cms.vuint32(1,1)
0098 process.MuMonitor_1e28.HLTOverallPrescale = cms.uint32(1)
0099 process.MuMonitor_1e28.throw = False
0100 process.MuMonitor_1e28.andOr = True
0101
0102 process.filterSdMuMonitor_1e28 = cms.Path(process.MuMonitor_1e28)
0103
0104
0105
0106 ### Mu SD
0107 process.Mu_1e28 = HLTrigger.HLTfilters.hltHighLevelDev_cfi.hltHighLevelDev.clone(andOr = True)
0108 process.Mu_1e28.HLTPaths = (
0109 "HLT_L2Mu0",
0110 "HLT_L2Mu3",
0111 #"HLT_L2Mu5",
0112 "HLT_L1Mu20",
0113 "HLT_L2Mu9",
0114 "HLT_L2Mu11",
0115 "HLT_L1Mu14_L1SingleEG10",
0116 "HLT_L1Mu14_L1SingleJet6U",
0117 "HLT_L1Mu14_L1ETM30",
0118 "HLT_L2DoubleMu0",
0119 "HLT_L1DoubleMuOpen",
0120 "HLT_DoubleMu0",
0121 "HLT_DoubleMu3",
0122 "HLT_Mu3",
0123 "HLT_Mu5",
0124 "HLT_Mu9",
0125 "HLT_IsoMu3",
0126 "HLT_Mu0_L1MuOpen",
0127 "HLT_Mu0_Track0_Jpsi",
0128 "HLT_Mu3_L1MuOpen",
0129 "HLT_Mu3_Track0_Jpsi",
0130 "HLT_Mu5_L1MuOpen",
0131 "HLT_Mu5_Track0_Jpsi",
0132 "HLT_Mu0_L2Mu0",
0133 "HLT_Mu3_L2Mu0",
0134 "HLT_Mu5_L2Mu0"
0135 )
0136 process.Mu_1e28.HLTPathsPrescales  = cms.vuint32(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1)
0137 process.Mu_1e28.HLTOverallPrescale = cms.uint32(1)
0138 process.Mu_1e28.throw = False
0139 process.Mu_1e28.andOr = True
0140
0141 process.filterSdMu_1e28 = cms.Path(process.Mu_1e28)
0142
0143
0144 ### EGMonitor SD
0145 process.EGMonitor_1e28 = HLTrigger.HLTfilters.hltHighLevelDev_cfi.hltHighLevelDev.clone(andOr = True)
0146 process.EGMonitor_1e28.HLTPaths = (
0147 "HLT_L1SingleEG2",
0148 "HLT_L1SingleEG5",
0149 "HLT_L1SingleEG8",
0150 "HLT_L1DoubleEG5",
0151 "HLT_EgammaSuperClusterOnly_L1R",
0152 "HLT_L1SingleEG20_NoBPTX",
0153 "HLT_L1SingleEG2_NoBPTX",
0154 "HLT_L1SingleEG5_NoBPTX"
0155 )
0156 process.EGMonitor_1e28.HLTPathsPrescales  = cms.vuint32(1,1,1,1,1,1,1,1)
0157 process.EGMonitor_1e28.HLTOverallPrescale = cms.uint32(1)
0158 process.EGMonitor_1e28.throw = False
0159 process.EGMonitor_1e28.andOr = True
0160
0161 process.filterSdEGMonitor_1e28 = cms.Path(process.EGMonitor_1e28)
0162
0163
0164
0165 ### EG SD
0166 process.EG_1e28 = HLTrigger.HLTfilters.hltHighLevelDev_cfi.hltHighLevelDev.clone(andOr = True)
0167 process.EG_1e28.HLTPaths = (
0168 "HLT_Photon10_L1R",
0169 "HLT_Photon15_L1R",
0170 "HLT_Photon15_LooseEcalIso_L1R",
0171 "HLT_Photon20_L1R",
0172 "HLT_Photon30_L1R_8E29",
0173 "HLT_DoublePhoton4_Jpsi_L1R",
0174 "HLT_DoublePhoton4_Upsilon_L1R",
0175 "HLT_DoublePhoton4_eeRes_L1R",
0176 "HLT_DoublePhoton5_eeRes_L1R", #added to match the /cdaq/physics/firstCollisions10/v2.0/HLT_7TeV/V5 table
0177 "HLT_DoublePhoton5_Jpsi_L1R",
0178 "HLT_DoublePhoton5_Upsilon_L1R",
0179 "HLT_DoublePhoton5_L1R",
0180 "HLT_DoublePhoton10_L1R",
0181 "HLT_DoubleEle5_SW_L1R",
0182 "HLT_Ele20_LW_L1R",
0183 "HLT_Ele15_SiStrip_L1R",
0184 "HLT_Ele15_SC10_LW_L1R",
0185 "HLT_Ele15_LW_L1R",
0186 "HLT_Ele10_LW_EleId_L1R",
0187 "HLT_Ele10_LW_L1R",
0188 "HLT_Photon15_TrackIso_L1R"
0189 )
0190 process.EG_1e28.HLTPathsPrescales  = cms.vuint32(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1)
0191 process.EG_1e28.HLTOverallPrescale = cms.uint32(1)
0192 process.EG_1e28.throw = False
0193 process.EG_1e28.andOr = True
0194
0195 process.filterSdEG_1e28 = cms.Path(process.EG_1e28)
0196
0197
0198
0199
0200 ### JetMET AOD CS
0201 process.DiJetAve_1e29 = HLTrigger.HLTfilters.hltHighLevelDev_cfi.hltHighLevelDev.clone(andOr = True)
0202 process.DiJetAve_1e29.HLTPaths = ("HLT_DiJetAve15U_8E29","HLT_DiJetAve30U_8E29")
0203 process.DiJetAve_1e29.HLTPathsPrescales  = cms.vuint32(1,1)
0204 process.DiJetAve_1e29.HLTOverallPrescale = cms.uint32(1)
0205 process.DiJetAve_1e29.andOr = True
0206
0207 process.filterCsDiJetAve_1e29 = cms.Path(process.DiJetAve_1e29)
0208
0209
0210
0211 ### Onia skim CS
0212 process.goodMuons = cms.EDFilter("MuonRefSelector",
0213     src = cms.InputTag("muons"),
0214     cut = cms.string("isGlobalMuon || (isTrackerMuon && numberOfMatches('SegmentAndTrackArbitration')>0)"),
0215 )
0216 process.diMuons = cms.EDProducer("CandViewShallowCloneCombiner",
0217     decay       = cms.string("goodMuons goodMuons"),
0218     checkCharge = cms.bool(False),
0219     cut         = cms.string("mass > 2"),
0220 )
0221 process.diMuonFilter = cms.EDFilter("CandViewCountFilter",
0222     src       = cms.InputTag("diMuons"),
0223     minNumber = cms.uint32(1),
0224 )
0225 process.Skim_diMuons = cms.Path(
0226     process.goodMuons    *
0227     process.diMuons      *
0228     process.diMuonFilter
0229 )
0230
0231
0232 ### Tau skim CS
0233 process.load("CondCore.DBCommon.CondDBSetup_cfi")
0234 process.load("TrackingTools/TransientTrack/TransientTrackBuilder_cfi")
0235
0236 process.load("L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff")
0237 process.load("HLTrigger/HLTfilters/hltLevel1GTSeed_cfi")
0238 process.hltLevel1GTSeed.L1TechTriggerSeeding = cms.bool(True)
0239 process.hltLevel1GTSeed.L1SeedsLogicalExpression = cms.string('(0 AND (40 OR 41) AND NOT (36 OR 37 OR 38 OR 39))')
0240
0241 process.scrapping = cms.EDFilter("FilterOutScraping",
0242         applyfilter = cms.untracked.bool(True),
0243         debugOn = cms.untracked.bool(False),
0244         numtrack = cms.untracked.uint32(10),
0245         thresh = cms.untracked.double(0.25)
0246 )
0247
0248 process.PFTausSelected = cms.EDFilter("PFTauSelector",
0249     src = cms.InputTag("shrinkingConePFTauProducer"),
0250     discriminators = cms.VPSet(
0251     cms.PSet( discriminator=cms.InputTag("shrinkingConePFTauDiscriminationByIsolation"),
0252            selectionCut=cms.double(0.5)
0253     ),
0254     ),
0255     cut = cms.string('et > 15. && abs(eta) < 2.5')
0256 )
0257
0258 process.PFTauSkimmed = cms.EDFilter("CandViewCountFilter",
0259   src = cms.InputTag('PFTausSelected'),
0260   minNumber = cms.uint32(1)
0261 )
0262
0263
0264 process.tauFilter = cms.Path(
0265     process.hltLevel1GTSeed *
0266     process.scrapping *
0267     process.PFTausSelected *
0268     process.PFTauSkimmed
0269 )
0270
0271
0272
0273 process.outputSdJetMETTau = cms.OutputModule("PoolOutputModule",
0274                                           SelectEvents = cms.untracked.PSet(SelectEvents = cms.vstring('filterSdJetMETTau_1e28')),
0275                                           dataset = cms.untracked.PSet(
0276         dataTier = cms.untracked.string('RECO'),
0277         filterName = cms.untracked.string('SD_JetMETTau')),
0278                                           outputCommands = process.RECOEventContent.outputCommands,
0279                                           fileName = cms.untracked.string('SD_JetMETTau_1e28.root')
0280                                           )
0281
0282 process.outputSdJetMETTauMonitor = cms.OutputModule("PoolOutputModule",
0283                                           SelectEvents = cms.untracked.PSet(SelectEvents = cms.vstring('filterSdJetMETTauMonitor_1e28')),
0284                                           dataset = cms.untracked.PSet(
0285         dataTier = cms.untracked.string('RECO'),
0286         filterName = cms.untracked.string('SD_JetMETTauMonitor')),
0287                                           outputCommands = process.RECOEventContent.outputCommands,
0288                                           fileName = cms.untracked.string('SD_JetMETTauMonitor_1e28.root')
0289                                           )
0290 process.outputSdMuMonitor = cms.OutputModule("PoolOutputModule",
0291                                           SelectEvents = cms.untracked.PSet(SelectEvents = cms.vstring('filterSdMuMonitor_1e28')),
0292                                           dataset = cms.untracked.PSet(
0293         dataTier = cms.untracked.string('RECO'),
0294         filterName = cms.untracked.string('SD_MuMonitor')),
0295                                           outputCommands = process.RECOEventContent.outputCommands,
0296                                           fileName = cms.untracked.string('SD_MuMonitor_1e28.root')
0297                                           )
0298
0299 process.outputSdMu = cms.OutputModule("PoolOutputModule",
0300                                           SelectEvents = cms.untracked.PSet(SelectEvents = cms.vstring('filterSdMu_1e28')),
0301                                           dataset = cms.untracked.PSet(
0302         dataTier = cms.untracked.string('RECO'),
0303         filterName = cms.untracked.string('SD_Mu')),
0304                                           outputCommands = process.RECOEventContent.outputCommands,
0305                                           fileName = cms.untracked.string('SD_Mu_1e28.root')
0306                                           )
0307 process.outputSdEGMonitor = cms.OutputModule("PoolOutputModule",
0308                                           SelectEvents = cms.untracked.PSet(SelectEvents = cms.vstring('filterSdEGMonitor_1e28')),
0309                                           dataset = cms.untracked.PSet(
0310         dataTier = cms.untracked.string('RECO'),
0311         filterName = cms.untracked.string('SD_EGMonitor')),
0312                                           outputCommands = process.RECOEventContent.outputCommands,
0313                                           fileName = cms.untracked.string('SD_EGMonitor_1e28.root')
0314                                           )
0315 process.outputSdEG = cms.OutputModule("PoolOutputModule",
0316                                           SelectEvents = cms.untracked.PSet(SelectEvents = cms.vstring('filterSdEG_1e28')),
0317                                           dataset = cms.untracked.PSet(
0318         dataTier = cms.untracked.string('RECO'),
0319         filterName = cms.untracked.string('SD_EG')),
0320                                           outputCommands = process.RECOEventContent.outputCommands,
0321                                           fileName = cms.untracked.string('SD_EG_1e28.root')
0322                                           )
0323
0324
0325 process.outputCsDiJet = cms.OutputModule("PoolOutputModule",
0326                                          dataset = cms.untracked.PSet(
0327         dataTier = cms.untracked.string('USER'),
0328         filterName = cms.untracked.string('CS_DiJetAve')),
0329                                          outputCommands = cms.untracked.vstring(
0330         'drop *',
0331         #------- CaloJet collections ------
0332         'keep recoCaloJets_kt4CaloJets_*_*',
0333         'keep recoCaloJets_kt6CaloJets_*_*',
0334         'keep recoCaloJets_ak5CaloJets_*_*',
0335         'keep recoCaloJets_ak7CaloJets_*_*',
0336         'keep recoCaloJets_iterativeCone5CaloJets_*_*',
0337         #------- CaloJet ID ---------------
0338         'keep *_kt4JetID_*_*',
0339         'keep *_kt6JetID_*_*',
0340         'keep *_ak5JetID_*_*',
0341         'keep *_ak7JetID_*_*',
0342         'keep *_ic5JetID_*_*',
0343         #------- PFJet collections ------
0344         'keep recoPFJets_kt4PFJets_*_*',
0345         'keep recoPFJets_kt6PFJets_*_*',
0346         'keep recoPFJets_ak5PFJets_*_*',
0347         'keep recoPFJets_ak7PFJets_*_*',
0348         'keep recoPFJets_iterativeCone5PFJets_*_*',
0349         #------- JPTJet collections ------
0350         'keep *_JetPlusTrackZSPCorJetAntiKt5_*_*',
0351         #'keep *_ak4JPTJets_*_*',
0352         #'keep *_iterativeCone5JPTJets_*_*',
0353         #------- Trigger collections ------
0354         'keep edmTriggerResults_TriggerResults_*_*',
0355         'keep *_hltTriggerSummaryAOD_*_*',
0356         'keep L1GlobalTriggerObjectMapRecord_*_*_*',
0357         'keep L1GlobalTriggerReadoutRecord_*_*_*',
0358         #------- Tracks collection --------
0359         'keep recoTracks_generalTracks_*_*',
0360         #------- CaloTower collection -----
0361         'keep *_towerMaker_*_*',
0362         #------- Various collections ------
0363         'keep *_EventAuxilary_*_*',
0364         'keep *_offlinePrimaryVertices_*_*',
0365         'keep *_hcalnoise_*_*',
0366         #------- MET collections ----------
0367         'keep *_metHO_*_*',
0368         'keep *_metNoHF_*_*',
0369         'keep *_metNoHFHO_*_*',
0370         'keep *_met_*_*'),
0371                                          SelectEvents = cms.untracked.PSet(SelectEvents = cms.vstring('filterCsDiJetAve_1e29')),
0372                                          fileName = cms.untracked.string('CS_JetAOD_DiJetAve_1e28.root')
0373                                          )
0374
0375 process.outputCsOnia = cms.OutputModule("PoolOutputModule",
0376                                         dataset = cms.untracked.PSet(
0377         dataTier = cms.untracked.string('RECO'),
0378         filterName = cms.untracked.string('CS_Onia')),
0379                                         SelectEvents = cms.untracked.PSet(SelectEvents = cms.vstring('Skim_diMuons')),
0380                                         outputCommands = process.RECOEventContent.outputCommands,
0381                                         fileName = cms.untracked.string('CS_Onia_1e28.root')
0382                                         )
0383
0384 process.outputCsTau = cms.OutputModule("PoolOutputModule",
0385                                         dataset = cms.untracked.PSet(
0386         dataTier = cms.untracked.string('RAW-RECO'),
0387         filterName = cms.untracked.string('CS_Tau')),
0388                                         SelectEvents = cms.untracked.PSet(SelectEvents = cms.vstring('tauFilter')),
0389                                         outputCommands = process.FEVTEventContent.outputCommands,
0390                                         fileName = cms.untracked.string('CS_Tau_1e28.root')
0391                                         )
0392
0393
0394
0395 process.this_is_the_end = cms.EndPath(
0396 process.outputSdJetMETTau        +
0397 process.outputSdJetMETTauMonitor +
0398 process.outputSdMuMonitor        +
0399 process.outputSdMu               +
0400 process.outputSdEGMonitor        +
0401 process.outputSdEG               +
0402 process.outputCsDiJet            +
0403 process.outputCsOnia             +
0404 process.outputCsTau
0405 )