test/python/proto_runIsolatedTracksNxNNzsData_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("L1SKIM")
0004
0005 process.load("FWCore.MessageService.MessageLogger_cfi")
0006 process.MessageLogger.cerr.FwkReport.reportEvery = 100000
0007
0008 process.options = cms.untracked.PSet( wantSummary = cms.untracked.bool(True) )
0009
0010 ####################### configure pool source #############################
0011 process.source = cms.Source("PoolSource",
0012  fileNames =cms.untracked.vstring(
0013     '/store/data/Run2010A/MinimumBias/RECO/Apr21ReReco-v1/0000/08275F4A-5270-E011-9DC3-003048635E02.root'
0014     ),
0015  skipEvents = cms.untracked.uint32(0)
0016 )
0017 process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(100) )
0018
0019 ##################### digi-2-raw plus L1 emulation #########################
0020 process.load("Configuration.StandardSequences.Services_cff")
0021 process.load('Configuration.StandardSequences.GeometryRecoDB_cff')
0022 process.load('Configuration.StandardSequences.MagneticField_cff')
0023
0024 #################### Conditions and L1 menu ################################
0025
0026 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0027 from Configuration.AlCa.autoCond import autoCond
0028 process.GlobalTag.globaltag=autoCond['run1_data']
0029
0030 ############ Skim the events according to the L1 seeds  ####################
0031
0032 #select on HLT_HcalNZS_8E29 trigger
0033 import HLTrigger.HLTfilters.hltLevel1GTSeed_cfi
0034 process.skimL1Seeds = HLTrigger.HLTfilters.hltLevel1GTSeed_cfi.hltLevel1GTSeed.clone()
0035 process.skimL1Seeds.L1GtReadoutRecordTag = cms.InputTag("gtDigis")
0036 process.skimL1Seeds.L1GtObjectMapTag     = cms.InputTag("hltL1GtObjectMap")
0037 process.skimL1Seeds.L1CollectionsTag     = cms.InputTag("l1extraParticles")
0038 process.skimL1Seeds.L1MuonCollectionTag  = cms.InputTag("l1extraParticles")
0039 process.skimL1Seeds.L1SeedsLogicalExpression = "L1_SingleEG2 OR L1_SingleEG5 OR L1_SingleEG8 OR L1_SingleEG10 OR L1_SingleEG12 OR L1_SingleEG15 OR L1_SingleEG20 OR L1_SingleIsoEG5 OR L1_SingleIsoEG8 OR L1_SingleIsoEG10 OR L1_SingleIsoEG12 OR L1_SingleIsoEG15 OR L1_SingleJet6U OR L1_SingleJet10U OR L1_SingleJet20U OR L1_SingleJet30U OR L1_SingleJet40U OR L1_SingleJet50U OR L1_SingleJet60U OR L1_SingleTauJet10U OR L1_SingleTauJet20U OR L1_SingleTauJet30U OR L1_SingleTauJet50U OR L1_SingleMuOpen OR L1_SingleMu0 OR L1_SingleMu3 OR L1_SingleMu5 OR L1_SingleMu7 OR L1_SingleMu10 OR L1_SingleMu14 OR L1_SingleMu20 OR L1_ZeroBias"
0040
0041 # select on HLT_HcalPhiSym trigger
0042 process.load("HLTrigger.HLTfilters.hltLevel1Activity_cfi")
0043 process.hltLevel1Activity.L1GtReadoutRecordTag  = cms.InputTag('gtDigis')
0044
0045 ######################## Configure Analyzer ###############################
0046
0047 process.load("RecoLocalCalo.EcalRecAlgos.EcalSeverityLevelESProducer_cfi")
0048 process.load("Calibration.IsolatedParticles.isolatedTracksNxN_cfi")
0049 process.isolatedTracksNxN.Verbosity = cms.untracked.int32( 0 )
0050 process.isolatedTracksNxN.HBHERecHitSource = cms.InputTag("hbhereco")
0051 process.isolatedTracksNxN.L1TriggerAlgoInfo = True
0052 #process.isolatedTracksNxN.DebugL1Info = True
0053
0054 process.isolatedTracksNxN_NZS = process.isolatedTracksNxN.clone(
0055     Verbosity = cms.untracked.int32( 0 ),
0056     HBHERecHitSource = cms.InputTag("hbherecoMB"),
0057     L1TriggerAlgoInfo = True
0058     )
0059
0060 process.TFileService = cms.Service("TFileService",
0061                                    fileName = cms.string('IsolatedTracksNxNData.root')
0062                                    )
0063
0064 # configure Technical Bits to ensure collision and remove BeamHalo
0065 process.load('L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff')
0066 process.load('HLTrigger/HLTfilters/hltLevel1GTSeed_cfi')
0067 process.hltLevel1GTSeed.L1TechTriggerSeeding = cms.bool(True)
0068 process.hltLevel1GTSeed.L1SeedsLogicalExpression = cms.string('0 AND NOT (36 OR 37 OR 38 OR 39)')
0069
0070 # filter out scrapping events
0071 process.noScraping= cms.EDFilter("FilterOutScraping",
0072                                  applyfilter = cms.untracked.bool(True),
0073                                  debugOn     = cms.untracked.bool(False), ## Or 'True' to get some per-event info
0074                                  numtrack    = cms.untracked.uint32(10),
0075                                  thresh      = cms.untracked.double(0.25)
0076                                  )
0077
0078 # select on primary vertex
0079 process.primaryVertexFilter = cms.EDFilter("GoodVertexFilter",
0080                                            vertexCollection = cms.InputTag('offlinePrimaryVertices'),
0081                                            minimumNDOF      = cms.uint32(4) ,
0082                                            maxAbsZ          = cms.double(25.0),
0083                                            maxd0            = cms.double(5.0)
0084                                            )
0085
0086
0087 #=============================================================================
0088 # define an EndPath to analyze all other path results
0089 process.hltTrigReport = cms.EDAnalyzer( 'HLTrigReport',
0090       HLTriggerResults = cms.InputTag( 'TriggerResults','','HLT')
0091 )
0092
0093 process.load("L1Trigger.GlobalTriggerAnalyzer.l1GtTrigReport_cfi")
0094 process.l1GtTrigReport.L1GtRecordInputTag = 'gtDigis'
0095 process.l1GtTrigReport.PrintVerbosity = 1
0096 #=============================================================================
0097
0098 #### by Benedikt
0099 process.p1 = cms.Path(process.primaryVertexFilter * process.hltLevel1GTSeed * process.noScraping * process.skimL1Seeds       *process.isolatedTracksNxN * process.isolatedTracksNxN_NZS)
0100 process.e = cms.EndPath(process.l1GtTrigReport + process.hltTrigReport)
0101