Skims/python/CosmicsPD_Skims.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("SKIM")
0004
0005 process.configurationMetadata = cms.untracked.PSet(
0006     version = cms.untracked.string('$Revision: 1.1 $'),
0007     name = cms.untracked.string('$Source: /cvs/CMSSW/CMSSW/DPGAnalysis/Skims/python/CosmicsPD_Skims.py,v $'),
0008     annotation = cms.untracked.string('CRAFT Cosmics PD skim')
0009 )
0010
0011 #
0012 #
0013 # This is for testing purposes.
0014 #
0015 #
0016 process.source = cms.Source("PoolSource",
0017                             fileNames = cms.untracked.vstring(
0018 '/store/data/Commissioning09/Cosmics/RECO/v5/000/105/755/D266D139-D871-DE11-A709-001D09F28F0C.root',
0019 '/store/data/Commissioning09/Cosmics/RECO/v5/000/105/755/CA27788D-E871-DE11-8B46-001D09F276CF.root',
0020 '/store/data/Commissioning09/Cosmics/RECO/v5/000/105/755/AC5633B2-D471-DE11-9B3A-001D09F252F3.root',
0021 '/store/data/Commissioning09/Cosmics/RECO/v5/000/105/755/9CD957E7-D071-DE11-B6AE-001D09F252F3.root',
0022 '/store/data/Commissioning09/Cosmics/RECO/v5/000/105/755/94BF68F7-D171-DE11-902B-000423D986A8.root',
0023 '/store/data/Commissioning09/Cosmics/RECO/v5/000/105/755/7838FE1E-C771-DE11-9FD5-000423D98950.root',
0024 '/store/data/Commissioning09/Cosmics/RECO/v5/000/105/755/56632803-DD71-DE11-BAF5-000423D9870C.root',
0025 '/store/data/Commissioning09/Cosmics/RECO/v5/000/105/755/42A67CB9-E971-DE11-AA86-001D09F252F3.root',
0026 '/store/data/Commissioning09/Cosmics/RECO/v5/000/105/755/407225D3-D071-DE11-809B-001D09F297EF.root',
0027 '/store/data/Commissioning09/Cosmics/RECO/v5/000/105/755/3E5E1CF0-D271-DE11-AC2B-000423D94700.root',
0028 '/store/data/Commissioning09/Cosmics/RECO/v5/000/105/755/2C57E916-D071-DE11-AF0E-001D09F24E39.root',
0029 '/store/data/Commissioning09/Cosmics/RECO/v5/000/105/755/228896A5-E571-DE11-A60B-001D09F2AF96.root'),
0030
0031 #  '/store/data/Commissioning08/Cosmics/RECO/CRAFT_ALL_V9_225-v1/0002/0A12CE23-D7F9-DD11-819E-00E081348D21.root'),
0032                             secondaryFileNames = cms.untracked.vstring(
0033 '/store/data/Commissioning09/Cosmics/RAW/v3/000/105/755/F6887FD0-9371-DE11-B69E-00304879FBB2.root',
0034 '/store/data/Commissioning09/Cosmics/RAW/v3/000/105/755/EEAF292E-9571-DE11-9A17-000423D94C68.root',
0035 '/store/data/Commissioning09/Cosmics/RAW/v3/000/105/755/ECC04DEB-9071-DE11-9F3A-001D09F23174.root',
0036 '/store/data/Commissioning09/Cosmics/RAW/v3/000/105/755/E8CE8710-9171-DE11-9211-000423D94534.root',
0037 '/store/data/Commissioning09/Cosmics/RAW/v3/000/105/755/DEB9057C-9471-DE11-BAF5-000423D94524.root',
0038 '/store/data/Commissioning09/Cosmics/RAW/v3/000/105/755/DC2A7158-A171-DE11-ACD3-001D09F24047.root',
0039 '/store/data/Commissioning09/Cosmics/RAW/v3/000/105/755/D6E68664-9271-DE11-AE97-000423D9970C.root',
0040 '/store/data/Commissioning09/Cosmics/RAW/v3/000/105/755/D2D94C9C-9171-DE11-AA96-000423D94AA8.root',
0041 '/store/data/Commissioning09/Cosmics/RAW/v3/000/105/755/B045C7CB-9371-DE11-AF39-001D09F24489.root',
0042 '/store/data/Commissioning09/Cosmics/RAW/v3/000/105/755/A2C791EF-9071-DE11-871D-001D09F2423B.root',
0043 '/store/data/Commissioning09/Cosmics/RAW/v3/000/105/755/9E1A4336-9071-DE11-81BD-001D09F251B8.root',
0044 '/store/data/Commissioning09/Cosmics/RAW/v3/000/105/755/8C28F89C-9171-DE11-B944-000423D9A2AE.root',
0045 '/store/data/Commissioning09/Cosmics/RAW/v3/000/105/755/82A37610-9371-DE11-A293-000423D98B6C.root',
0046 '/store/data/Commissioning09/Cosmics/RAW/v3/000/105/755/822D1510-9371-DE11-A329-000423D99A8E.root',
0047 '/store/data/Commissioning09/Cosmics/RAW/v3/000/105/755/76515611-9371-DE11-8BF2-001D09F24024.root',
0048 '/store/data/Commissioning09/Cosmics/RAW/v3/000/105/755/76434EA1-9171-DE11-B39F-000423D98E54.root',
0049 '/store/data/Commissioning09/Cosmics/RAW/v3/000/105/755/6E1DB87B-9471-DE11-B4E8-000423D98B28.root',
0050 '/store/data/Commissioning09/Cosmics/RAW/v3/000/105/755/6A15A1C7-9371-DE11-9ACB-0030487A18F2.root',
0051 '/store/data/Commissioning09/Cosmics/RAW/v3/000/105/755/4E95107B-9471-DE11-9E86-000423D944F0.root',
0052 '/store/data/Commissioning09/Cosmics/RAW/v3/000/105/755/42C723CA-9871-DE11-8614-000423D6B48C.root',
0053 '/store/data/Commissioning09/Cosmics/RAW/v3/000/105/755/3E934456-9271-DE11-BA19-000423D99AA2.root',
0054 '/store/data/Commissioning09/Cosmics/RAW/v3/000/105/755/3AF89B38-9371-DE11-AEE9-000423D33970.root',
0055 '/store/data/Commissioning09/Cosmics/RAW/v3/000/105/755/2C696C5B-9271-DE11-BD57-000423D94990.root')
0056
0057 #        '/store/data/Commissioning08/Cosmics/RAW/v1/000/069/578/085EFED4-E5AB-DD11-9ACA-001617C3B6FE.root')
0058 )
0059
0060 process.maxEvents = cms.untracked.PSet(
0061     input = cms.untracked.int32(500)
0062 )
0063
0064 process.options = cms.untracked.PSet(
0065     wantSummary = cms.untracked.bool(True)
0066 )
0067
0068 #------------------------------------------
0069 # Load standard sequences.
0070 #------------------------------------------
0071 process.load('Configuration/StandardSequences/MagneticField_AutoFromDBCurrent_cff')
0072 process.load('Configuration/StandardSequences/GeometryIdeal_cff')
0073
0074
0075 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0076 process.GlobalTag.globaltag = 'CRAFT09_R_V4::All'
0077
0078 process.load("Configuration/StandardSequences/RawToDigi_Data_cff")
0079 process.load("Configuration/StandardSequences/ReconstructionCosmics_cff")
0080
0081 process.load('Configuration/EventContent/EventContentCosmics_cff')
0082
0083
0084 import DPGAnalysis.Skims.CSCSkim_cfi
0085 process.extend(DPGAnalysis.Skims.CSCSkim_cfi)
0086
0087 ############################CSCSkim trial ##########################################
0088
0089 #### the path
0090
0091 process.mySkim = cms.Path(process.cscSkim)
0092
0093
0094 #### output
0095 process.outputSkim = cms.OutputModule(
0096     "PoolOutputModule",
0097     outputCommands = cms.untracked.vstring('keep *','drop *_MEtoEDMConverter_*_*'),
0098     fileName = cms.untracked.string("cscskimEvents.root"),
0099     dataset = cms.untracked.PSet(
0100       dataTier = cms.untracked.string('RAW-RECO'),
0101       filterName = cms.untracked.string('CSCSkim_trial')
0102     ),
0103     SelectEvents = cms.untracked.PSet(SelectEvents = cms.vstring('mySkim'))
0104 )
0105
0106 ##########################CSC Skim BField Studies#######################################
0107
0108 process.cscbfield=DPGAnalysis.Skims.CSCSkim_cfi.cscSkim.clone()
0109 process.cscbfield.typeOfSkim = cms.untracked.int32(9)
0110
0111 #### the path
0112
0113 process.BfieldStudySkim = cms.Path(process.cscbfield)
0114
0115
0116 #### output
0117 process.outputSkimBField = cms.OutputModule(
0118     "PoolOutputModule",
0119     outputCommands = cms.untracked.vstring('keep *','drop *_MEtoEDMConverter_*_*'),
0120     fileName = cms.untracked.string("CSCEvents_BFieldStudy.root"),
0121     dataset = cms.untracked.PSet(
0122       dataTier = cms.untracked.string('RAW-RECO'),
0123       filterName = cms.untracked.string('CSCSkim_BFieldStudies')
0124     ),
0125     SelectEvents = cms.untracked.PSet(SelectEvents = cms.vstring('BfieldStudySkim'))
0126 )
0127
0128
0129
0130 ######################################Super Pointing#################################################
0131
0132 process.cosmicMuonsBarrelOnlyFilter = cms.EDFilter("HLTMuonPointingFilter",
0133     SALabel = cms.InputTag("cosmicMuonsBarrelOnly"),
0134     PropagatorName = cms.string("SteppingHelixPropagatorAny"),
0135     radius = cms.double(10.0),
0136     maxZ = cms.double(50.0),
0137 )
0138
0139 process.cosmicMuonsFilter = cms.EDFilter("HLTMuonPointingFilter",
0140     SALabel = cms.InputTag("cosmicMuons"),
0141     PropagatorName = cms.string("SteppingHelixPropagatorAny"),
0142     radius = cms.double(10.0),
0143     maxZ = cms.double(50.0),
0144 )
0145
0146 process.cosmicMuons1LegFilter = cms.EDFilter("HLTMuonPointingFilter",
0147     SALabel = cms.InputTag("cosmicMuons1Leg"),
0148     PropagatorName = cms.string("SteppingHelixPropagatorAny"),
0149     radius = cms.double(10.0),
0150     maxZ = cms.double(50.0),
0151 )
0152
0153 process.globalCosmicMuonsBarrelOnlyFilter = cms.EDFilter("HLTMuonPointingFilter",
0154     SALabel = cms.InputTag("globalCosmicMuonsBarrelOnly"),
0155     PropagatorName = cms.string("SteppingHelixPropagatorAny"),
0156     radius = cms.double(10.0),
0157     maxZ = cms.double(50.0),
0158 )
0159
0160 process.cosmictrackfinderP5Filter = cms.EDFilter("HLTMuonPointingFilter",
0161     SALabel = cms.InputTag("cosmictrackfinderP5"),
0162     PropagatorName = cms.string("SteppingHelixPropagatorAny"),
0163     radius = cms.double(10.0),
0164     maxZ = cms.double(50.0),
0165 )
0166
0167 process.globalCosmicMuonsFilter = cms.EDFilter("HLTMuonPointingFilter",
0168     SALabel = cms.InputTag("globalCosmicMuons"),
0169     PropagatorName = cms.string("SteppingHelixPropagatorAny"),
0170     radius = cms.double(10.0),
0171     maxZ = cms.double(50.0),
0172 )
0173
0174 process.rsWithMaterialTracksP5Filter = cms.EDFilter("HLTMuonPointingFilter",
0175     SALabel = cms.InputTag("rsWithMaterialTracksP5"),
0176     PropagatorName = cms.string("SteppingHelixPropagatorAny"),
0177     radius = cms.double(10.0),
0178     maxZ = cms.double(50.0),
0179 )
0180
0181 process.globalCosmicMuons1LegFilter = cms.EDFilter("HLTMuonPointingFilter",
0182     SALabel = cms.InputTag("globalCosmicMuons1Leg"),
0183     PropagatorName = cms.string("SteppingHelixPropagatorAny"),
0184     radius = cms.double(10.0),
0185     maxZ = cms.double(50.0),
0186 )
0187
0188 process.ctfWithMaterialTracksP5Filter = cms.EDFilter("HLTMuonPointingFilter",
0189     SALabel = cms.InputTag("ctfWithMaterialTracksP5"),
0190     PropagatorName = cms.string("SteppingHelixPropagatorAny"),
0191     radius = cms.double(10.0),
0192     maxZ = cms.double(50.0),
0193 )
0194
0195
0196 process.cosmicMuonsBarrelOnlyPath = cms.Path(process.cosmicMuonsBarrelOnlyFilter)
0197 process.cosmicMuonsPath = cms.Path(process.cosmicMuonsFilter)
0198 process.cosmicMuons1LegPath = cms.Path(process.cosmicMuons1LegFilter)
0199 process.globalCosmicMuonsBarrelOnlyPath = cms.Path(process.globalCosmicMuonsBarrelOnlyFilter)
0200 process.cosmictrackfinderP5Path = cms.Path(process.cosmictrackfinderP5Filter)
0201 process.globalCosmicMuonsPath = cms.Path(process.globalCosmicMuonsFilter)
0202 process.rsWithMaterialTracksP5Path = cms.Path(process.rsWithMaterialTracksP5Filter)
0203 process.globalCosmicMuons1LegPath = cms.Path(process.globalCosmicMuons1LegFilter)
0204 process.ctfWithMaterialTracksP5Path = cms.Path(process.ctfWithMaterialTracksP5Filter)
0205
0206
0207
0208 process.outSP = cms.OutputModule("PoolOutputModule",
0209                                outputCommands = cms.untracked.vstring('keep *','drop *_MEtoEDMConverter_*_*'),
0210                                SelectEvents = cms.untracked.PSet(SelectEvents = cms.vstring('cosmicMuonsBarrelOnlyPath',
0211                                                                                             'cosmicMuonsPath',
0212                                                                                             'cosmicMuons1LegPath',
0213                                                                                             'globalCosmicMuonsBarrelOnlyPath',
0214                                                                                             'cosmictrackfinderP5Path',
0215                                                                                             'globalCosmicMuonsPath',
0216                                                                                             'rsWithMaterialTracksP5Path',
0217                                                                                             'globalCosmicMuons1LegPath',
0218                                                                                             'ctfWithMaterialTracksP5Path')),
0219                                dataset = cms.untracked.PSet(
0220                              dataTier = cms.untracked.string('RAW-RECO'),
0221                                          filterName = cms.untracked.string('SuperPointing')),
0222                                fileName = cms.untracked.string('/tmp/arizzi/superPointing.root')
0223                                )
0224
0225 #################################Tracker Pointing###############################################
0226 process.cosmictrackfinderP5TkCntFilter = cms.EDFilter("TrackCountFilter",
0227                                                       src = cms.InputTag('cosmictrackfinderP5'),
0228                                                       minNumber = cms.uint32(1)
0229                                                       )
0230
0231 process.ctfWithMaterialTracksP5TkCntFilter = cms.EDFilter("TrackCountFilter",
0232                                                           src = cms.InputTag('ctfWithMaterialTracksP5'),
0233                                                           minNumber = cms.uint32(1)
0234                                                           )
0235
0236 process.rsWithMaterialTracksP5TkCntFilter = cms.EDFilter("TrackCountFilter",
0237                                                          src = cms.InputTag('rsWithMaterialTracksP5'),
0238                                                          minNumber = cms.uint32(1)
0239                                                          )
0240
0241 process.cosmicMuonsBarrelOnlyTkFilter = cms.EDFilter("HLTMuonPointingFilter",
0242                                                      SALabel = cms.InputTag("cosmicMuonsBarrelOnly"),
0243                                                      PropagatorName = cms.string("SteppingHelixPropagatorAny"),
0244                                                      radius = cms.double(90.0),
0245                                                      maxZ = cms.double(130.0)
0246                                                      )
0247
0248 process.cosmicMuonsBarrelOnlyTkPath = cms.Path(process.cosmicMuonsBarrelOnlyTkFilter)
0249 process.cosmictrackfinderP5TkCntPath = cms.Path(process.cosmictrackfinderP5TkCntFilter)
0250 process.ctfWithMaterialTracksP5TkCntPath = cms.Path(process.ctfWithMaterialTracksP5TkCntFilter)
0251 process.rsWithMaterialTracksP5TkCntPath = cms.Path(process.rsWithMaterialTracksP5TkCntFilter)
0252
0253
0254 process.outTP = cms.OutputModule("PoolOutputModule",
0255                                outputCommands = cms.untracked.vstring('keep *','drop *_MEtoEDMConverter_*_*'),
0256                                SelectEvents = cms.untracked.PSet(SelectEvents = cms.vstring('cosmicMuonsBarrelOnlyTkPath',
0257                                                                                             'cosmictrackfinderP5TkCntPath',
0258                                                                                             'ctfWithMaterialTracksP5TkCntPath',
0259                                                                                             'rsWithMaterialTracksP5TkCntPath')),
0260                                dataset = cms.untracked.PSet(
0261                              dataTier = cms.untracked.string('RAW-RECO'),
0262                                          filterName = cms.untracked.string('TrackingPointing')),
0263                                fileName = cms.untracked.string('/tmp/arizzi/trackerPointing.root')
0264                                )
0265
0266
0267 ####################################### Multi Muon ###########################################################
0268
0269 process.multiCosmicMuonFilter = cms.EDFilter("TrackCountFilter",
0270                                              src = cms.InputTag('cosmicMuonsBarrelOnly'),
0271                                              minNumber = cms.uint32(5)
0272                                              )
0273
0274 process.multiLHCMuonFilter = cms.EDFilter("TrackCountFilter",
0275                                           src = cms.InputTag('lhcStandAloneMuonsBarrelOnly'),
0276                                           minNumber = cms.uint32(5)
0277                                           )
0278
0279 process.multiCosmicMuonPath = cms.Path(process.multiCosmicMuonFilter)
0280 process.multiLHCMuonPath = cms.Path(process.multiLHCMuonFilter)
0281
0282 process.outMM = cms.OutputModule("PoolOutputModule",
0283                                 process.RECOEventContent,
0284                                SelectEvents = cms.untracked.PSet(SelectEvents = cms.vstring('multiCosmicMuonPath',
0285                                                                                             'multiLHCMuonPath')),
0286                                dataset = cms.untracked.PSet(
0287                              dataTier = cms.untracked.string('RECO'),
0288                                          filterName = cms.untracked.string('multiCosmicMuon')),
0289                                fileName = cms.untracked.string('/tmp/arizzi/multiMuon.root')
0290                                )
0291
0292 process.outMM.outputCommands.append('drop *_MEtoEDMConverter_*_*')
0293
0294
0295
0296
0297
0298 import FWCore.ParameterSet.Config as cms
0299 import DPGAnalysis.Skims.HCALHighEnergyCombinedPath_cff
0300
0301 #################################HCAL High Energy#########################################
0302
0303 process.extend(DPGAnalysis.Skims.HCALHighEnergyCombinedPath_cff)
0304
0305 process.outHCAL = cms.OutputModule("PoolOutputModule",
0306                                outputCommands = cms.untracked.vstring('keep *','drop *_MEtoEDMConverter_*_*'),
0307                                SelectEvents = cms.untracked.PSet(
0308     SelectEvents = cms.vstring("HCALHighEnergyPath")
0309     ),
0310                                dataset = cms.untracked.PSet(
0311                                dataTier = cms.untracked.string('RAW-RECO'),
0312                                filterName = cms.untracked.string('HCALHighEnergy')),
0313                                fileName = cms.untracked.string('/tmp/arizzi/HCALHighEnergy_filter.root')
0314                                )
0315
0316 #################################ECAL High Energy#########################################
0317 process.skimming = cms.EDFilter("EcalSkim",
0318     #cosmic cluster energy threshold in GeV
0319     energyCutEB = cms.untracked.double(2.0),
0320     energyCutEE = cms.untracked.double(2.0),
0321     endcapClusterCollection = cms.InputTag("cosmicSuperClusters","CosmicEndcapSuperClusters"),
0322     barrelClusterCollection = cms.InputTag("cosmicSuperClusters","CosmicBarrelSuperClusters")
0323 )
0324 process.ecalPath = cms.Path(process.skimming)
0325
0326 process.outECAL = cms.OutputModule("PoolOutputModule",
0327     process.RECOEventContent,
0328     fileName = cms.untracked.string('/tmp/arizzi/ecalSkim.root'),
0329     dataset = cms.untracked.PSet(
0330               dataTier = cms.untracked.string('RECO'),
0331               filterName = cms.untracked.string('ecalSkim_fromRECO')),
0332     SelectEvents = cms.untracked.PSet(
0333         SelectEvents = cms.vstring('ecalPath')
0334     )
0335 )
0336
0337 process.outECAL.outputCommands.append('drop *_MEtoEDMConverter_*_*')
0338
0339 ###############################StoppedHSCP################################
0340
0341 process.outHSCP = cms.OutputModule("PoolOutputModule",
0342                                outputCommands = cms.untracked.vstring('keep *','drop *_MEtoEDMConverter_*_*'),
0343                                SelectEvents = cms.untracked.PSet(
0344     SelectEvents = cms.vstring("HLT_StoppedHSCP_*:HLT")
0345     ),
0346                                dataset = cms.untracked.PSet(
0347                                dataTier = cms.untracked.string('RAW-RECO'),
0348                                filterName = cms.untracked.string('StoppedHSCP')),
0349                                fileName = cms.untracked.string('/tmp/arizzi/StoppedHSCP_filter.root')
0350                                )
0351
0352
0353 process.outpath = cms.EndPath(process.outSP+process.outTP+process.outMM+process.outputSkim+process.outputSkimBField+process.outHCAL+process.outECAL+process.outHSCP)