test/python/proto_runIsolatedTracksNxNData_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("L1SKIM")
0004
0005 process.load("RecoLocalCalo.EcalRecAlgos.EcalSeverityLevelESProducer_cfi")
0006 process.load("FWCore.MessageService.MessageLogger_cfi")
0007 process.MessageLogger.cerr.FwkReport.reportEvery = 1000000
0008 process.MessageLogger.L1GtTrigReport=dict()
0009 process.MessageLogger.HLTrigReport=dict()
0010
0011 process.options = cms.untracked.PSet(
0012     wantSummary = cms.untracked.bool(True)
0013 )
0014
0015 process.source = cms.Source("PoolSource",fileNames =cms.untracked.vstring(
0016     "/store/data/Run2012A/MinimumBias/RECO/13Jul2012-v1/00000/06BDF16E-1ACF-E111-A543-00304867BED8.root"))
0017 process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(5) )
0018
0019 ##################### digi-2-raw plus L1 emulation #########################
0020
0021 process.load("Configuration.StandardSequences.Services_cff")
0022 process.load('Configuration.StandardSequences.GeometryRecoDB_cff')
0023 process.load('Configuration.StandardSequences.MagneticField_cff')
0024
0025 #################### Conditions and L1 menu ################################
0026 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0027 from Configuration.AlCa.autoCond import autoCond
0028 process.GlobalTag.globaltag=autoCond['run1_data']
0029
0030
0031 #################################################################################################
0032 process.load("Calibration.IsolatedParticles.isolatedTracksNxN_cfi")
0033 process.isolatedTracksNxN.JetSource = 'ak5CaloJets'
0034 process.TFileService = cms.Service("TFileService",
0035                                    fileName = cms.string('IsolatedTracksNxNData.root')
0036                                    )
0037 #=============================================================================
0038 # select the  events with trigger type
0039 # configure Technical Bits
0040 process.load('L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff')
0041 process.load('HLTrigger/HLTfilters/hltLevel1GTSeed_cfi')
0042 process.hltLevel1GTSeed.L1TechTriggerSeeding = cms.bool(True)
0043 process.hltLevel1GTSeed.L1SeedsLogicalExpression = cms.string('0 AND NOT (36 OR 37 OR 38 OR 39) AND NOT ((42 AND NOT 43) OR (43 AND NOT 42))')
0044 #process.hltLevel1GTSeed.L1SeedsLogicalExpression = cms.string('0 AND (40 OR 41) AND NOT (36 OR 37 OR 38 OR 39) AND NOT ((42 AND NOT 43) OR (43 AND NOT 42))')
0045
0046 # select on MinBias trigger
0047 import HLTrigger.HLTfilters.hltHighLevel_cfi
0048 process.myHLT = HLTrigger.HLTfilters.hltHighLevel_cfi.hltHighLevel.clone()
0049 process.myHLT.HLTPaths = cms.vstring("HLT_MinBiasBSC","HLT_L1Tech_BSC_minBias")
0050 process.myHLT.throw    = cms.bool(False)
0051
0052 # filter out scrapping events
0053 process.noScraping= cms.EDFilter("FilterOutScraping",
0054                                  applyfilter = cms.untracked.bool(True),
0055                                  debugOn     = cms.untracked.bool(False), ## Or 'True' to get some per-event info
0056                                  numtrack    = cms.untracked.uint32(10),
0057                                  thresh      = cms.untracked.double(0.25)
0058                                  )
0059
0060 process.primaryVertexFilter = cms.EDFilter("GoodVertexFilter",
0061                                            vertexCollection = cms.InputTag('offlinePrimaryVertices'),
0062                                            minimumNDOF      = cms.uint32(4) ,
0063                                            maxAbsZ          = cms.double(25.0),
0064                                            maxd0            = cms.double(2.0)
0065                                            )
0066
0067
0068 #=============================================================================
0069 # define an EndPath to analyze all other path results
0070 process.hltTrigReport = cms.EDAnalyzer( 'HLTrigReport',
0071       HLTriggerResults = cms.InputTag( 'TriggerResults','','HLT')
0072 )
0073
0074 process.load("L1Trigger.GlobalTriggerAnalyzer.l1GtTrigReport_cfi")
0075 #process.l1GtTrigReport.L1GtRecordInputTag = 'simGtDigis'
0076 process.l1GtTrigReport.L1GtRecordInputTag = 'gtDigis'
0077 process.l1GtTrigReport.PrintVerbosity = 1
0078 #=============================================================================
0079
0080 process.p1 = cms.Path( process.primaryVertexFilter * process.hltLevel1GTSeed * process.noScraping * process.myHLT * process.isolatedTracksNxN )
0081
0082 #process.e = cms.EndPath(process.l1GtTrigReport + process.hltTrigReport)
0083 #process.e = cms.EndPath(process.l1GtTrigReport)