Reference Manual

Go to the documentation of this file.
00001 import FWCore.ParameterSet.Config as cms
00002 
00003 process = cms.Process("SKIM")
00004 
00005 process.configurationMetadata = cms.untracked.PSet(
00006     version = cms.untracked.string('$Revision: 1.36 $'),
00007     name = cms.untracked.string('$Source: /cvs_server/repositories/CMSSW/CMSSW/DPGAnalysis/Skims/python/MinBiasPDSkim_cfg.py,v $'),
00008     annotation = cms.untracked.string('Combined MinBias skim')
00009 )
00010 # selection eff. on 1000 events
00011 #file:/tmp/azzi/Background.root
00012 #/tmp/azzi/Background.root ( 0 events, 664191 bytes )
00013 #file:/tmp/azzi/MinBiascscskimEvents.root
00014 #/tmp/azzi/MinBiascscskimEvents.root ( 0 events, 664212 bytes )
00015 #file:/tmp/azzi/MuonDPGSkim.root
00016 #/tmp/azzi/MuonDPGSkim.root ( 95 events, 40887080 bytes )
00017 #file:/tmp/azzi/StoppedHSCP_filter.root
00018 #/tmp/azzi/StoppedHSCP_filter.root ( 124 events, 28507733 bytes )
00019 #file:/tmp/azzi/ValSkim.root
00020 #/tmp/azzi/ValSkim.root ( 22 events, 9588726 bytes )
00021 #file:/tmp/azzi/ecalrechitfilter.root
00022 #/tmp/azzi/ecalrechitfilter.root ( 23 events, 14583758 bytes )
00023 #file:/tmp/azzi/logerror_filter.root
00024 #/tmp/azzi/logerror_filter.root ( 14 events, 9389296 bytes )
00025 #file:/tmp/azzi/TPGSkim.root
00026 #/tmp/azzi/TPGSkim.root ( 46 events, 20271643 bytes )
00027 
00028 #
00029 #
00030 # This is for testing purposes.
00031 #
00032 #
00033 process.source = cms.Source("PoolSource",
00034                             fileNames = cms.untracked.vstring(
00035 # run 136066 lumi~500
00036 '/store/data/Run2010A/MinimumBias/RECO/v1/000/136/066/18F6DB82-5566-DF11-B289-0030487CAF0E.root'),
00037                            secondaryFileNames = cms.untracked.vstring(
00038 '/store/data/Run2010A/MinimumBias/RAW/v1/000/136/066/38D48BED-3C66-DF11-88A5-001D09F27003.root')
00039 )
00040 
00041 process.source.inputCommands = cms.untracked.vstring("keep *", "drop *_MEtoEDMConverter_*_*")
00042 
00043 process.maxEvents = cms.untracked.PSet(
00044     input = cms.untracked.int32(-1)
00045 )
00046 
00047 
00048 #------------------------------------------
00049 # Load standard sequences.
00050 #------------------------------------------
00051 process.load('Configuration/StandardSequences/MagneticField_AutoFromDBCurrent_cff')
00052 process.load('Configuration/StandardSequences/GeometryIdeal_cff')
00053 
00054 
00055 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
00056 process.GlobalTag.globaltag = 'GR_R_38X_V13::All' 
00057 
00058 process.load("Configuration/StandardSequences/RawToDigi_Data_cff")
00059 process.load("Configuration/StandardSequences/Reconstruction_cff")
00060 process.load('Configuration/EventContent/EventContent_cff')
00061 
00062 #drop collections created on the fly
00063 process.FEVTEventContent.outputCommands.append("drop *_MEtoEDMConverter_*_*")
00064 process.FEVTEventContent.outputCommands.append("drop *_*_*_SKIM")
00065 
00066 #
00067 #  Load common sequences
00068 #
00069 process.load('L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskAlgoTrigConfig_cff')
00070 process.load('L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff')
00071 process.load('HLTrigger/HLTfilters/hltLevel1GTSeed_cfi')
00072 
00073 ###########################################################################################
00074 #------------------------------------------
00075 # parameters for the CSCSkim module
00076 #------------------------------------------
00077 process.load("DPGAnalysis/Skims/CSCSkim_cfi")
00078 
00079 
00080 #set to minimum activity
00081 process.cscSkim.minimumSegments = 1
00082 process.cscSkim.minimumHitChambers = 1
00083 
00084 # this is for filtering on HLT path
00085 process.hltBeamHalo = cms.EDFilter("HLTHighLevel",
00086      TriggerResultsTag = cms.InputTag("TriggerResults","","HLT"),
00087      HLTPaths = cms.vstring('HLT_CSCBeamHalo','HLT_CSCBeamHaloOverlapRing1','HLT_CSCBeamHaloOverlapRing','HLT_CSCBeamHaloRing2or3'), # provide list of HLT paths (or patterns) you want
00088      eventSetupPathsKey = cms.string(''), # not empty => use read paths from AlCaRecoTriggerBitsRcd via this key
00089      andOr = cms.bool(True),             # how to deal with multiple triggers: True (OR) accept if ANY is true, False (AND) accept if ALL are true
00090      throw = cms.bool(False)    # throw exception on unknown path names
00091  )
00092 
00093 #### the path
00094 process.cscHaloSkim = cms.Path(process.hltBeamHalo+process.cscSkim)
00095 
00096 
00097 
00098 #### output 
00099 process.outputBeamHaloSkim = cms.OutputModule("PoolOutputModule",
00100     outputCommands = process.FEVTEventContent.outputCommands,
00101     fileName = cms.untracked.string("/tmp/azzi/MinBiascscskimEvents.root"),
00102     dataset = cms.untracked.PSet(
00103       dataTier = cms.untracked.string('RAW-RECO'),
00104       filterName = cms.untracked.string('CSCSkim_BeamHalo_MinBias')
00105     ),
00106     SelectEvents = cms.untracked.PSet(SelectEvents = cms.vstring('cscHaloSkim'))
00107 )
00108 
00109 
00110 ########################## Muon tracks Filter ############################
00111 process.muonSkim=cms.EDFilter("CandViewCountFilter", 
00112                  src =cms.InputTag("muons"), minNumber = cms.uint32(1))
00113 process.muonTracksSkim = cms.Path(process.muonSkim)
00114 
00115 
00116 ###########################################################################
00117 
00118 process.outputMuonSkim = cms.OutputModule("PoolOutputModule",
00119     fileName = cms.untracked.string('/tmp/azzi/MuonSkim.root'),
00120     outputCommands = cms.untracked.vstring('keep *','drop *_MEtoEDMConverter_*_*'),
00121     dataset = cms.untracked.PSet(
00122               dataTier = cms.untracked.string('RECO'),
00123               filterName = cms.untracked.string('Muon_skim')),
00124     SelectEvents = cms.untracked.PSet(
00125         SelectEvents = cms.vstring("muonTracksSkim")
00126     )
00127 )
00128 ####################################################################################
00129 
00130 ##################################stoppedHSCP############################################
00131 
00132 
00133 # this is for filtering on HLT path
00134 process.hltstoppedhscp = cms.EDFilter("HLTHighLevel",
00135      TriggerResultsTag = cms.InputTag("TriggerResults","","HLT"),
00136      HLTPaths = cms.vstring("HLT_StoppedHSCP*"), # provide list of HLT paths (or patterns) you want
00137      eventSetupPathsKey = cms.string(''), # not empty => use read paths from AlCaRecoTriggerBitsRcd via this key
00138      andOr = cms.bool(True),             # how to deal with multiple triggers: True (OR) accept if ANY is true, False (AND) accept if ALL are true
00139      throw = cms.bool(False)    # throw exception on unknown path names
00140  )
00141 
00142 process.HSCP=cms.Path(process.hltstoppedhscp)
00143 
00144 process.outHSCP = cms.OutputModule("PoolOutputModule",
00145                                outputCommands =  process.FEVTEventContent.outputCommands,
00146                                fileName = cms.untracked.string('/tmp/azzi/StoppedHSCP_filter.root'),
00147                                dataset = cms.untracked.PSet(
00148                                   dataTier = cms.untracked.string('RAW-RECO'),
00149                                   filterName = cms.untracked.string('Skim_StoppedHSCP')),
00150                                
00151                                SelectEvents = cms.untracked.PSet(
00152     SelectEvents = cms.vstring("HSCP")
00153     ))
00154 
00155 ###########################################################################################
00156 #------------------------------------------
00157 # parameters for the PFGCollisions skim3
00158 #------------------------------------------
00159 process.load('HLTrigger.special.hltPhysicsDeclared_cfi')
00160 process.hltPhysicsDeclared.L1GtReadoutRecordTag = 'gtDigis'
00161 
00162 process.hltbeamgas = cms.EDFilter("HLTHighLevel",
00163 TriggerResultsTag = cms.InputTag("TriggerResults","","HLT"),
00164   HLTPaths = cms.vstring('HLT_L1_BptxXOR_BscMinBiasOR'), # provide list of HLT paths (or patterns) you want
00165   eventSetupPathsKey = cms.string(''),
00166   andOr              = cms.bool(True),
00167   throw              = cms.bool(False)
00168 
00169 )
00170 
00171 #### the path
00172 process.pfgskim3noncross = cms.Path(process.hltPhysicsDeclared*process.hltbeamgas)
00173 
00174 #### output
00175 process.outputpfgskim3 = cms.OutputModule("PoolOutputModule",
00176  outputCommands = process.FEVTEventContent.outputCommands,
00177  fileName = cms.untracked.string("/tmp/azzi/Background.root"),
00178  dataset = cms.untracked.PSet(
00179    dataTier = cms.untracked.string('RAW-RECO'),
00180    filterName = cms.untracked.string('BEAMBKGV3')
00181  ),
00182  SelectEvents = cms.untracked.PSet(SelectEvents = cms.vstring('pfgskim3noncross'))
00183 )
00184 
00185 ###########################################################################################
00186 ###########################################################################################
00187 
00188 #===========================================================
00189 
00190 #################################logerrorharvester############################################
00191 process.load("FWCore.Modules.logErrorFilter_cfi")
00192 from Configuration.StandardSequences.RawToDigi_Data_cff import gtEvmDigis
00193 
00194 process.gtEvmDigis = gtEvmDigis.clone()
00195 process.stableBeam = cms.EDFilter("HLTBeamModeFilter",
00196                                   L1GtEvmReadoutRecordTag = cms.InputTag("gtEvmDigis"),
00197                                   AllowedBeamMode = cms.vuint32(11)
00198                                   )
00199 
00200 process.logerrorpath=cms.Path(process.gtEvmDigis+process.stableBeam+process.logErrorFilter)
00201 
00202 process.outlogerr = cms.OutputModule("PoolOutputModule",
00203                                outputCommands =  process.FEVTEventContent.outputCommands,
00204                                fileName = cms.untracked.string('/tmp/azzi/logerror_filter.root'),
00205                                dataset = cms.untracked.PSet(
00206                                   dataTier = cms.untracked.string('RAW-RECO'),
00207                                   filterName = cms.untracked.string('Skim_logerror')),
00208                                
00209                                SelectEvents = cms.untracked.PSet(
00210     SelectEvents = cms.vstring("logerrorpath")
00211     ))
00212 
00213 #===========================================================
00214 ###########################ngood event per lumi##########################################
00215 process.primaryVertexFilter = cms.EDFilter("VertexSelector",
00216    src = cms.InputTag("offlinePrimaryVertices"),
00217    cut = cms.string("!isFake && ndof > 4 && abs(z) <= 15 && position.Rho <= 2"), # tracksSize() > 3 for the older cut
00218    filter = cms.bool(True),   # otherwise it won't filter the events, just produce an empty vertex collection.
00219 )
00220 
00221 
00222 process.noscraping = cms.EDFilter("FilterOutScraping",
00223 applyfilter = cms.untracked.bool(True),
00224 debugOn = cms.untracked.bool(False),
00225 numtrack = cms.untracked.uint32(10),
00226 thresh = cms.untracked.double(0.25)
00227 )
00228 
00229 ###Tracks selection
00230 process.trackSelector  =cms.EDFilter("trackSelector",
00231                                     src = cms.InputTag("generalTracks"),
00232                                      cut = cms.string('quality("highPurity")')     
00233                                      )
00234 
00235 #process.trackSelector = cms.EDProducer("QualityFilter",
00236 #                                       TrackQuality = cms.string('highPurity'),
00237 #                                       recTracks = cms.InputTag("generalTracks")
00238 #                                       )
00239 
00240 process.trackFilter = cms.EDFilter("TrackCountFilter",
00241                                    src = cms.InputTag("trackSelector"),
00242                                    minNumber = cms.uint32(10)
00243                                    )
00244 
00245 process.nottoomanytracks = cms.EDFilter("NMaxPerLumi",
00246                                         NMaxPerLumi = cms.uint32(8)
00247                                         )
00248 process.relvaltrackskim = cms.Path(process.primaryVertexFilter+process.noscraping+
00249                                    process.trackSelector + process.trackFilter + process.nottoomanytracks )
00250 
00251 ### muon selection
00252 process.MuonSelector = cms.EDFilter("MuonSelector",
00253                                     src = cms.InputTag("muons"),
00254                                     cut = cms.string(" isGlobalMuon && isTrackerMuon && pt > 3")
00255                                     )
00256 process.muonFilter = cms.EDFilter("MuonCountFilter",
00257                                   src = cms.InputTag("MuonSelector"),
00258                                   minNumber = cms.uint32(1)
00259                                   )
00260 process.nottoomanymuons = cms.EDFilter("NMaxPerLumi",
00261                                        NMaxPerLumi = cms.uint32(2)
00262                                        )
00263 process.relvalmuonskim = cms.Path(process.primaryVertexFilter+process.noscraping+
00264                                   process.MuonSelector + process.muonFilter + process.nottoomanymuons )
00265 
00266 #### output 
00267 process.outputvalskim = cms.OutputModule("PoolOutputModule",
00268                                          outputCommands = process.FEVTEventContent.outputCommands,
00269                                          fileName = cms.untracked.string("/tmp/azzi/ValSkim.root"),
00270                                          dataset = cms.untracked.PSet(
00271     dataTier = cms.untracked.string('RAW-RECO'),
00272     filterName = cms.untracked.string('valskim')
00273     ),
00274                                          SelectEvents = cms.untracked.PSet(SelectEvents = cms.vstring('relvaltrackskim','relvalmuonskim')
00275                                                                            ))
00276 
00277 
00278 ###########################################################################
00279 ######################################TPG Performance SKIMS#####################################
00280 
00281 process.load('DPGAnalysis/Skims/singleMuonSkim_cff')
00282 process.load('DPGAnalysis/Skims/singleElectronSkim_cff')
00283 process.load('DPGAnalysis/Skims/muonTagProbeFilters_cff')
00284 process.load('DPGAnalysis/Skims/electronTagProbeFilters_cff')
00285 process.load('DPGAnalysis/Skims/singlePhotonSkim_cff')
00286 process.load('DPGAnalysis/Skims/jetSkim_cff')
00287 process.load('DPGAnalysis/Skims/METSkim_cff')
00288 process.load('DPGAnalysis/Skims/singlePfTauSkim_cff')
00289 
00290 #process.singleMuPt20SkimPath=cms.Path(process.singleMuPt20RecoQualitySeq)
00291 #process.singleMuPt15SkimPath=cms.Path(process.singleMuPt15RecoQualitySeq)
00292 #process.singleMuPt10SkimPath=cms.Path(process.singleMuPt10RecoQualitySeq)
00293 process.singleMuPt5SkimPath=cms.Path(process.singleMuPt5RecoQualitySeq)
00294 #process.singleElectronPt20SkimPath=cms.Path(process.singleElectronPt20RecoQualitySeq)
00295 #process.singleElectronPt15SkimPath=cms.Path(process.singleElectronPt15RecoQualitySeq)
00296 #process.singleElectronPt10SkimPath=cms.Path(process.singleElectronPt10RecoQualitySeq)
00297 process.singleElectronPt5SkimPath=cms.Path(process.singleElectronPt5RecoQualitySeq)
00298 #process.singlePhotonPt20SkimPath=cms.Path(process.singlePhotonPt20QualitySeq)
00299 #process.singlePhotonPt15SkimPath=cms.Path(process.singlePhotonPt15QualitySeq)
00300 #process.singlePhotonPt10SkimPath=cms.Path(process.singlePhotonPt10QualitySeq)
00301 process.singlePhotonPt5SkimPath=cms.Path(process.singlePhotonPt5QualitySeq)
00302 #process.muonZMMSkimPath=cms.Path(process.muonZMMRecoQualitySeq)
00303 process.muonJPsiMMSkimPath=cms.Path(process.muonJPsiMMRecoQualitySeq)
00304 #process.electronZEESkimPath=cms.Path(process.electronZEERecoQualitySeq)
00305 process.jetSkimPath=cms.Path(process.jetRecoQualitySeq)
00306 #process.METSkimPath=cms.Path(process.METQualitySeq)
00307 process.singlePfTauPt15SkimPath=cms.Path(process.singlePfTauPt15QualitySeq) 
00308 
00309 process.outTPGSkim = cms.OutputModule("PoolOutputModule",
00310     outputCommands = process.FEVTHLTALLEventContent.outputCommands,
00311     fileName = cms.untracked.string("/tmp/azzi/TPGSkim.root"),
00312     dataset = cms.untracked.PSet(
00313       dataTier = cms.untracked.string('USER'),
00314       filterName = cms.untracked.string('TPGSkim')
00315     ),
00316     SelectEvents = cms.untracked.PSet(SelectEvents = cms.vstring(
00317                                                                  #'singleMuPt20SkimPath',
00318                                                                  #'singleMuPt15SkimPath',
00319                                                                  #'singleMuPt10SkimPath',
00320                                                                  'singleMuPt5SkimPath',
00321                                                                  #'singleElectronPt20SkimPath',
00322                                                                  #'singleElectronPt15SkimPath',
00323                                                                  #'singleElectronPt10SkimPath',
00324                                                                  'singleElectronPt5SkimPath',
00325                                                                  #'singlePhotonPt20SkimPath',
00326                                                                  #'singlePhotonPt15SkimPath',
00327                                                                  #'singlePhotonPt10SkimPath',
00328                                                                  'singlePhotonPt5SkimPath',
00329                                                                  #'muonZMMSkimPath',
00330                                                                  'muonJPsiMMSkimPath',
00331                                                                  #'electronZEESkimPath',
00332                                                                  'jetSkimPath',
00333                                                                  #'METSkimPath',
00334                                                                  'singlePfTauPt15SkimPath'))
00335 )
00336 
00337 
00338 ###########################################################################
00339 
00340 
00341 process.options = cms.untracked.PSet(
00342  wantSummary = cms.untracked.bool(True)
00343 )
00344 
00345 #process.outpath = cms.EndPath(process.outputBeamHaloSkim+process.outputMuonSkim+process.outHSCP+process.outputpfgskim3+process.outlogerr+process.outputvalskim+process.outTPGSkim)
00346 #BeamHalo removed
00347 process.outpath = cms.EndPath(process.outputMuonSkim+process.outHSCP+process.outputpfgskim3+process.outlogerr+process.outputvalskim+process.outTPGSkim)
00348 
00349 
00350
/data/doxygen/doxygen-1.7.3/gen/CMSSW_4_2_8/src/DPGAnalysis/Skims/python/MinBiasPDSkim_cfg.py
