process.MessageLogger.cerr.FwkReport.reportEvery = 100
process.MessageLogger.cerr.default.limit = 100
# jet corrections
#process.load("JetMETCorrections.Configuration.L2L3Corrections_Summer09_7TeV_ReReco332_cff")
process.load('JetMETCorrections.Configuration.DefaultJEC_cff')
#process.load("JetMETCorrections.Configuration.L2L3Corrections_900GeV_cff")
#process.load("JetMETCorrections.Configuration.L2L3Corrections_2360GeV_cff")
# summary
process.options = cms.untracked.PSet( wantSummary = cms.untracked.bool(True) )
process.load('L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff')
from HLTrigger.HLTfilters.hltLevel1GTSeed_cfi import hltLevel1GTSeed
process.bit40OR41 = hltLevel1GTSeed.clone(L1TechTriggerSeeding = cms.bool(True), L1SeedsLogicalExpression = cms.string('40 OR 41'))
from HLTrigger.HLTfilters.hltHighLevelDev_cfi import hltHighLevelDev
process.physDecl = hltHighLevelDev.clone(HLTPaths = ['HLT_PhysicsDeclared'], HLTPathsPrescales = [1])
process.promptanaTree = cms.EDAnalyzer("PromptAnaTree",
outputCommands = cms.untracked.vstring(
'drop *',
'keep *_promptanaevent_*_*',
'keep *_promptanamet_*_*',
'keep *_promptanatcmet_*_*',
'keep *_promptanapfmet_*_*',
'keep *_promptananohf_*_*',
'keep *_promptanaic5calojet_*_*',
#'keep *_promptanasc5calojet_*_*',
#'keep *_promptanakt4calojet_*_*',
import FWCore.ParameterSet.Config as cms
# ================ Physics Declared Bit ================
physDeclFilter = cms.EDFilter("PhysDecl",
applyfilter = cms.untracked.bool(False),
HLTriggerResults = cms.InputTag("TriggerResults","","HLT")
)
# ================ L1 Filters =====================
from L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff import *
from HLTrigger.HLTfilters.hltLevel1GTSeed_cfi import hltLevel1GTSeed
from L1Trigger.Skimmer.l1Filter_cfi import l1Filter
bptxAnd = hltLevel1GTSeed.clone(
L1TechTriggerSeeding = cms.bool(True),
L1SeedsLogicalExpression = cms.string('0')
)
bscOr = hltLevel1GTSeed.clone(
L1TechTriggerSeeding = cms.bool(True),
L1SeedsLogicalExpression = cms.string('34 OR 40 OR 41')
)
bscNoHalo = hltLevel1GTSeed.clone(
L1TechTriggerSeeding = cms.bool(True),
L1SeedsLogicalExpression = cms.string('NOT (36 OR 37 OR 38 OR 39)')
)
bscOrBptxOr = l1Filter.clone(
algorithms = ["L1_BscMinBiasOR_BptxPlusORMinus"]
)
PAprimaryVertexFilter = cms.EDFilter("VertexSelector",
src = cms.InputTag("offlinePrimaryVertices"),
cut = cms.string("!isFake && abs(z) <= 25 && position.Rho <= 2 && tracksSize >= 2"),
filter = cms.bool(True), # otherwise it won't filter the events
)
# Cluster-shape filter re-run offline
from RecoLocalTracker.SiPixelRecHits.SiPixelRecHits_cfi import *
from HLTrigger.special.hltPixelClusterShapeFilter_cfi import *
hltPixelClusterShapeFilter.inputTag = "siPixelRecHits"
# Reject BSC beam halo L1 technical bits
from L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff import *
from HLTrigger.HLTfilters.hltLevel1GTSeed_cfi import hltLevel1GTSeed
noBSChalo = hltLevel1GTSeed.clone(
L1TechTriggerSeeding = cms.bool(True),
L1SeedsLogicalExpression = cms.string('NOT (36 OR 37 OR 38 OR 39)')
)
#Reject beam scraping events standard pp configuration
NoScraping = cms.EDFilter("FilterOutScraping",
applyfilter = cms.untracked.bool(True),
debugOn = cms.untracked.bool(False),
numtrack = cms.untracked.uint32(10),
thresh = cms.untracked.double(0.25)
)
collisionEventSelection = cms.Sequence(noBSChalo *
hfCoincFilter3 *
primaryVertexFilter *
siPixelRecHits *
hltPixelClusterShapeFilter)
# '/store/data/Run2011A/MinimumBias/RECO/May10ReReco-v2/0002/C0CE5857-8482-E011-9BD9-002354EF3BDE.root',
# '/store/data/Run2011A/MinimumBias/RECO/May10ReReco-v2/0002/BC5BF6CB-6682-E011-B5E2-00248C0BE01E.root',
# '/store/data/Run2011A/MinimumBias/RECO/May10ReReco-v2/0002/AA6C29E6-6F82-E011-8F12-001A92971B62.root',
# '/store/data/Run2011A/MinimumBias/RECO/May10ReReco-v2/0002/A476C008-8A82-E011-A676-002354EF3BDA.root'
# )
# )
# Configure the trigger
process.load('L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff')
from HLTrigger.HLTfilters.hltLevel1GTSeed_cfi import hltLevel1GTSeed
# L1 technical triggers:
# 0 -- beam crossing
# 36-39 -- beam halo
process.VetoHardInt = hltLevel1GTSeed.clone(
L1TechTriggerSeeding = cms.bool(True),
L1SeedsLogicalExpression = cms.string('0 AND NOT (36 OR 37 OR 38 OR 39)')
)
# Remove the so-called "scraping" events
process.noScraping = cms.EDFilter("FilterOutScraping",
applyfilter = cms.untracked.bool(True),
debugOn = cms.untracked.bool(False), ## Or 'True' to get some per-event info
numtrack = cms.untracked.uint32(10),
thresh = cms.untracked.double(0.25)
)
# min bias filter
process.HLTZeroBias =cms.EDFilter("HLTHighLevel",
TriggerResultsTag = cms.InputTag("TriggerResults","","HLT"),
HLTPaths = cms.vstring('HLT_ZeroBias','HLT_ZeroBias_v1','HLT_ZeroBias_v2','HLT_ZeroBias_v3'),
eventSetupPathsKey = cms.string(''),
process.load("DQMServices.Core.DQM_cfg")
process.load("RecoBTag.Configuration.RecoBTag_cff")
process.load('Configuration.StandardSequences.Services_cff')
process.load('Configuration.StandardSequences.GeometryExtended_cff')
process.load('Configuration.StandardSequences.MagneticField_AutoFromDBCurrent_cff')
process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
#Global tag for 3_6_1
process.GlobalTag.globaltag = 'GR_R_36X_V12A::All'
process.load('L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff')
from HLTrigger.HLTfilters.hltLevel1GTSeed_cfi import hltLevel1GTSeed
#Good Bunch Crossings
process.bptxAnd = hltLevel1GTSeed.clone(L1TechTriggerSeeding = cms.bool(True), L1SeedsLogicalExpression = cms.string('0'))
#BSCNOBEAMHALO
process.bit40 = hltLevel1GTSeed.clone(L1TechTriggerSeeding = cms.bool(True), L1SeedsLogicalExpression = cms.string('(40 OR 41) AND NOT (36 OR 37 OR 38 OR 39) AND NOT ((42 AND NOT 43) OR (43 AND NOT 42))'))
#Physics-declared Bit
from HLTrigger.HLTfilters.hltHighLevelDev_cfi import hltHighLevelDev
process.physDecl = hltHighLevelDev.clone(HLTPaths = ['HLT_PhysicsDeclared'], HLTPathsPrescales = [1])
#Select events based on the HLT triggers (HLT_L1Jet6U || HLT_L1Jet10U || HLT_Jet15U)
# Use the instructions provided at:
# http://cmssw.cvs.cern.ch/cgi-bin/cmssw.cgi/CMSSW/HLTrigger/HLTfilters/python/hltHighLevel_cfi.py?hideattic=1&revision=1.5&view=markup
import HLTrigger.HLTfilters.hltHighLevel_cfi
process.singleJetHLTFilter = HLTrigger.HLTfilters.hltHighLevel_cfi.hltHighLevel.clone()
process.singleJetHLTFilter.TriggerResultsTag = cms.InputTag("TriggerResults","","HLT")
process.singleJetHLTFilter.HLTPaths = ["HLT_L1Jet6U", "HLT_L1Jet10U", "HLT_Jet15U"]
process.singleJetHLTFilter.andOr = cms.bool(True) # how to deal with multiple triggers: True (OR) accept if ANY is true, False (AND) accept if ALL are true
andOr = cms.bool(True), # how to deal with multiple triggers: True (OR) accept if ANY is true, False (AND) accept if ALL are true
throw = cms.bool(False) # throw exception on unknown path names
)
process.dtHLTSkim = cms.Path(process.HLTDT)
process.dtSkim=cms.Path(process.muonDTDigis+process.hltDTActivityFilter)
###########################################################################
from HLTrigger.HLTfilters.hltLevel1GTSeed_cfi import hltLevel1GTSeed
############################ L1 Muon bits #################################
process.l1RequestPhAlgos = hltLevel1GTSeed.clone()
# False allows to read directly from L1 instead fo candidate ObjectMap
process.l1RequestPhAlgos.L1UseL1TriggerObjectMaps = cms.bool(False)
#
# option used forL1UseL1TriggerObjectMaps = False only
# number of BxInEvent: 1: L1A=0; 3: -1, L1A=0, 1; 5: -2, -1, L1A=0, 1,
# online is used 5
process.l1RequestPhAlgos.L1NrBxInEvent = cms.int32(5)
# Request the or of the following bits: from 54 to 62 and 106-107
process.l1RequestPhAlgos.L1SeedsLogicalExpression = cms.string(
'L1_SingleMuBeamHalo 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_DoubleMuOpen OR L1_DoubleMu3')
process.l1MuBitsSkim = cms.Path(process.l1RequestPhAlgos)
process.load('Configuration.StandardSequences.Services_cff')
process.load('Configuration.StandardSequences.GeometryExtended_cff')
process.load(
'Configuration.StandardSequences.MagneticField_AutoFromDBCurrent_cff')
process.load(
'Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
#Global tag for 3_6_1
process.GlobalTag.globaltag = 'GR_R_36X_V12A::All'
process.load(
'L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff')
from HLTrigger.HLTfilters.hltLevel1GTSeed_cfi import hltLevel1GTSeed
#Good Bunch Crossings
process.bptxAnd = hltLevel1GTSeed.clone(
L1TechTriggerSeeding=cms.bool(True),
L1SeedsLogicalExpression=cms.string('0'))
#BSCNOBEAMHALO
process.bit40 = hltLevel1GTSeed.clone(
L1TechTriggerSeeding=cms.bool(True),
L1SeedsLogicalExpression=cms.string(
'(40 OR 41) AND NOT (36 OR 37 OR 38 OR 39) AND NOT ((42 AND NOT 43) OR (43 AND NOT 42))'
))
#Physics-declared Bit
from HLTrigger.HLTfilters.hltHighLevelDev_cfi import hltHighLevelDev
process.physDecl = hltHighLevelDev.clone(HLTPaths=['HLT_PhysicsDeclared'],
HLTPathsPrescales=[1])
#Select events based on the HLT triggers (HLT_L1Jet6U || HLT_L1Jet10U || HLT_Jet15U)
process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1) )
process.options = cms.untracked.PSet(wantSummary = cms.untracked.bool(True))
process.GlobalTag.globaltag = 'GR10_P_V12::All'
# =============== Trigger selection ====================
from HLTrigger.HLTfilters.hltHighLevel_cfi import hltHighLevel
process.whateverHLT = hltHighLevel.clone(
#HLTPaths = cms.vstring('HLT_HIJet90U'),
HLTPaths = cms.vstring('HLT_HIMinBiasHF'),
andOr = cms.bool(True)
)
from HLTrigger.HLTfilters.hltLevel1GTSeed_cfi import hltLevel1GTSeed
process.whateverL1 = hltLevel1GTSeed.clone(
L1TechTriggerSeeding = cms.bool(True),
L1SeedsLogicalExpression = cms.string('35')
)
#process.eventFilter_step = cms.Path(process.whateverHLT)
#process.eventFilter_step = cms.Path(process.whateverL1)
# =============== Trigger summary =======================
# L1
process.load('L1Trigger.GlobalTriggerAnalyzer.l1GtTrigReport_cfi')
process.l1GtTrigReport.L1GtRecordInputTag = cms.InputTag( "gtDigis") # simGtDigis
# HLT
process.load('HLTrigger.HLTanalyzers.hlTrigReport_cfi')
process.TriggerAnalyzerEndpath = cms.EndPath(process.l1GtTrigReport*
process.hlTrigReport)
cut=cms.string(
"!isFake && abs(z) <= 25 && position.Rho <= 2 && tracksSize >= 2"),
filter=cms.bool(True), # otherwise it won't filter the events
)
# Cluster-shape filter re-run offline
process.load("RecoLocalTracker.SiPixelRecHits.SiPixelRecHits_cfi")
process.load("HLTrigger.special.hltPixelClusterShapeFilter_cfi")
process.hltPixelClusterShapeFilter.inputTag = "siPixelRecHits"
# Reject BSC beam halo L1 technical bits
from L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff import *
from HLTrigger.HLTfilters.hltLevel1GTSeed_cfi import hltLevel1GTSeed
process.noBSChalo = hltLevel1GTSeed.clone(
L1TechTriggerSeeding=cms.bool(True),
L1SeedsLogicalExpression=cms.string('NOT (36 OR 37 OR 38 OR 39)'))
process.collisionEventSelection = cms.Sequence(
process.noBSChalo * process.hfCoincFilter3 * process.primaryVertexFilter *
process.siPixelRecHits * process.hltPixelClusterShapeFilter)
process.load("HLTrigger.HLTfilters.hltHighLevel_cfi")
### JetHI SD
process.hltJetHI = process.hltHighLevel.clone(HLTPaths=['HLT_HIJet35U'])
process.filterSdJetHI = cms.Path(
process.bscOrHfCoinc * process.collisionEventSelection * process.hltJetHI)
### PhotonHI SD
process.hltPhotonHI = process.hltHighLevel.clone(HLTPaths=['HLT_HIPhoton15'])
process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(10))
process.options = cms.untracked.PSet(wantSummary = cms.untracked.bool(True))
# ================ Physics Declared ================
process.physDeclFilter = cms.EDFilter("PhysDecl",
applyfilter = cms.untracked.bool(True),
debugOn = cms.untracked.bool(True),
HLTriggerResults = cms.InputTag("TriggerResults","","HLT")
)
# ================ bptx =======================
process.load('L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff')
from HLTrigger.HLTfilters.hltLevel1GTSeed_cfi import hltLevel1GTSeed
process.bptxAnd = hltLevel1GTSeed.clone(
L1TechTriggerSeeding = cms.bool(True),
L1SeedsLogicalExpression = cms.string('0')
)
# ================ BSC, beam halo veto =======================
# no work here, running on bsc+nohalo skimmed sample
#process.L1T1=process.hltLevel1GTSeed.clone()
#process.L1T1.L1TechTriggerSeeding = cms.bool(True)
#process.L1T1.L1SeedsLogicalExpression = cms.string('(40 OR 41) AND NOT (36 OR 37 OR 38 OR 39)')
# ================ monster event rejection =======================
process.monsterVeto = cms.EDFilter("FilterOutScraping",
applyfilter = cms.untracked.bool(True),
debugOn = cms.untracked.bool(True),
numtrack = cms.untracked.uint32(10),
### source
process.source = cms.Source("PoolSource",
fileNames = cms.untracked.vstring(
#'root://pcmssd12.cern.ch//data/gpetrucc/7TeV/hlt/MuHLT_MinBiasMC357_185_1.root',
'root://pcmssd12.cern.ch//data/gpetrucc/7TeV/jpsi/ppMuX_Spring10_REDIGI_START3X_V26_S09_GEN-SIM-RECO_C0AC7DEB-8144-DF11-A1E1-00304867D838.root'
)
)
### number of events
process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(100) )
### FILTERS for GoodCollision
process.load('L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff')
from HLTrigger.HLTfilters.hltLevel1GTSeed_cfi import hltLevel1GTSeed
hltLevel1GTSeed.L1TechTriggerSeeding = cms.bool(True)
process.bscFilter = hltLevel1GTSeed.clone(L1SeedsLogicalExpression = cms.string('(40 OR 41) AND NOT (36 OR 37 OR 38 OR 39)'))
process.bit40 = hltLevel1GTSeed.clone(L1SeedsLogicalExpression = cms.string('(40 OR 41)'))
process.haloVeto = hltLevel1GTSeed.clone(L1SeedsLogicalExpression = cms.string('NOT (36 OR 37 OR 38 OR 39)'))
process.oneGoodVertexFilter = cms.EDFilter("VertexSelector",
src = cms.InputTag("offlinePrimaryVertices"),
cut = cms.string("!isFake && ndof >= 4 && abs(z) <= 15 && position.Rho <= 2"),
filter = cms.bool(True), # otherwise it won't filter the events, just produce an empty vertex collection.
)
process.noScraping = cms.EDFilter("FilterOutScraping",
applyfilter = cms.untracked.bool(True),
debugOn = cms.untracked.bool(False), ## Or 'True' to get some per-event info
numtrack = cms.untracked.uint32(10),
thresh = cms.untracked.double(0.25)
)
process.preFilter = cms.Sequence( process.oneGoodVertexFilter * process.noScraping )
process.oneRecoMu = cms.EDFilter("CandViewCountFilter", src = cms.InputTag("muons"), minNumber = cms.uint32(1))
### source
process.source = cms.Source("PoolSource",
fileNames = cms.untracked.vstring(
#'rfio:/castor/cern.ch/user/g/gpetrucc/7TeV/DATA/Muons_With_HLTDEBUG.root'
#'file:/afs/cern.ch/user/g/gpetrucc/scratch0/mu10/clean/CMSSW_3_5_6/src/hlt.root'
'rfio:/castor/cern.ch/user/g/gpetrucc/7TeV/DATA/Muons_With_HLTDEBUG_v9_Run134542_Ls44to52.root'
)
)
### number of events
process.maxEvents = cms.untracked.PSet( output = cms.untracked.int32(100) )
from HLTrigger.HLTfilters.hltLevel1GTSeed_cfi import hltLevel1GTSeed
hltLevel1GTSeed.L1TechTriggerSeeding = cms.bool(True)
process.bptxAnd = hltLevel1GTSeed.clone(L1SeedsLogicalExpression = cms.string('0'))
process.bscFilter = hltLevel1GTSeed.clone(L1SeedsLogicalExpression = cms.string('(40 OR 41) AND NOT (36 OR 37 OR 38 OR 39)'))
process.oneGoodVertexFilter = cms.EDFilter("VertexSelector",
src = cms.InputTag("offlinePrimaryVertices"),
cut = cms.string("!isFake && ndof >= 4 && abs(z) <= 15 && position.Rho <= 2"),
filter = cms.bool(True), # otherwise it won't filter the events, just produce an empty vertex collection.
)
process.noScraping = cms.EDFilter("FilterOutScraping",
applyfilter = cms.untracked.bool(True),
debugOn = cms.untracked.bool(False), ## Or 'True' to get some per-event info
numtrack = cms.untracked.uint32(10),
thresh = cms.untracked.double(0.25)
)
#process.preFilter = cms.Sequence(process.bptxAnd * process.bscFilter * process.oneGoodVertexFilter * process.noScraping)
process.preFilter = cms.Sequence(process.noScraping * process.oneGoodVertexFilter)
process.test = cms.Path(process.preFilter)
TrackAssociatorByHits.Cut_RecoToSim = 0.5
from Tests.MaterialNtuplizer.ConversionNtuplizer_cfi import conv
from Tests.MaterialNtuplizer.ConversionNtuplizer_cfi import newConv
from Tests.MaterialNtuplizer.NuclIntNtuplizer_cfi import nucl
#convHit = conv.clone(hitassoc = cms.bool(True))
#====================================================================================
#
# New wrt PAS: correct selection of the HLT for Minimum Bias
#most similar to altTrig for MC is bit34
from L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff import *
from HLTrigger.HLTfilters.hltLevel1GTSeed_cfi import hltLevel1GTSeed
bit34 = hltLevel1GTSeed.clone(L1TechTriggerSeeding = cms.bool(True), L1SeedsLogicalExpression = cms.string('34'))
#Below not applicable for MC
###unmasking
#from L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff import *
#from HLTrigger.HLTfilters.hltLevel1GTSeed_cfi import hltLevel1GTSeed
#es_prefer_l1GtTriggerMaskTechTrig = cms.ESPrefer("L1GtTriggerMaskTechTrigTrivialProducer","l1GtTriggerMaskTechTrig")
###bit 0 selection: BPTX_AND
#bptxAnd = hltLevel1GTSeed.clone(L1TechTriggerSeeding = cms.bool(True), L1SeedsLogicalExpression = cms.string('0'))
###HLT
#from HLTrigger.HLTfilters.hltHighLevel_cfi import *
#hltMinimumBiasSelection = hltHighLevel.clone(HLTPaths = cms.vstring("HLT_L1_BscMinBiasOR_BptxPlusORMinus","HLT_L1Tech_BSC_minBias_OR"), throw = cms.bool(False))
process.main = cms.Path(process.upsilonSkimMainSequence)
## ____ _ _ _ _ _____ _ ____ _ _ _
## / ___|___ | | (_)___(_) ___ _ __ | ____|_ _____ _ __ | |_ / ___| ___| | ___ ___| |_(_) ___ _ __
## | | / _ \| | | / __| |/ _ \| '_ \ | _| \ \ / / _ \ '_ \| __| \___ \ / _ \ |/ _ \/ __| __| |/ _ \| '_ \
## | |__| (_) | | | \__ \ | (_) | | | | | |___ \ V / __/ | | | |_ ___) | __/ | __/ (__| |_| | (_) | | | |
## \____\___/|_|_|_|___/_|\___/|_| |_| |_____| \_/ \___|_| |_|\__| |____/ \___|_|\___|\___|\__|_|\___/|_| |_|
##
##
# Some event filters
process.load(
'L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff')
from HLTrigger.HLTfilters.hltLevel1GTSeed_cfi import hltLevel1GTSeed
process.bit40 = hltLevel1GTSeed.clone(
L1TechTriggerSeeding=cms.bool(True),
L1SeedsLogicalExpression=cms.string('(40 OR 41)'))
process.haloVeto = hltLevel1GTSeed.clone(
L1TechTriggerSeeding=cms.bool(True),
L1SeedsLogicalExpression=cms.string('NOT (36 OR 37 OR 38 OR 39)'))
process.bptxAnd = hltLevel1GTSeed.clone(
L1TechTriggerSeeding=cms.bool(True),
L1SeedsLogicalExpression=cms.string('0'))
from HLTrigger.HLTfilters.hltHighLevelDev_cfi import hltHighLevelDev
process.physDecl = hltHighLevelDev.clone(HLTPaths=['HLT_PhysicsDeclared'],
HLTPathsPrescales=[1])
process.oneGoodVertexFilter = cms.EDFilter(
"VertexSelector",
src=cms.InputTag("offlinePrimaryVertices"),
'root://pcmssd12.cern.ch//data/gpetrucc/Feb9Skims/MC_CollisionEvents_MuonSkim_1.root'
),
inputCommands=RECOSIMEventContent.
outputCommands, # keep only RECO out of RAW+RECO, for tests
dropDescendantsOfDroppedBranches=cms.untracked.bool(
False), # keep only RECO out of RAW+RECO, for tests
)
process.maxEvents = cms.untracked.PSet(input=cms.untracked.int32(-1))
process.options = cms.untracked.PSet(wantSummary=cms.untracked.bool(True))
process.load(
'L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff')
from HLTrigger.HLTfilters.hltLevel1GTSeed_cfi import hltLevel1GTSeed
hltLevel1GTSeed.L1TechTriggerSeeding = cms.bool(True)
process.bscFilter = hltLevel1GTSeed.clone(L1SeedsLogicalExpression=cms.string(
'(40 OR 41) AND NOT (36 OR 37 OR 38 OR 39)'))
process.oneGoodVertexFilter = cms.EDFilter(
"VertexSelector",
src=cms.InputTag("offlinePrimaryVertices"),
cut=cms.string(
"!isFake && ndof >= 4 && abs(z) <= 15 && position.Rho <= 2"),
filter=cms.bool(
True
), # otherwise it won't filter the events, just produce an empty vertex collection.
)
process.noScraping = cms.EDFilter(
"FilterOutScraping",
applyfilter=cms.untracked.bool(True),
debugOn=cms.untracked.bool(False), ## Or 'True' to get some per-event info
numtrack=cms.untracked.uint32(10),
thresh=cms.untracked.double(0.25))
),
SelectEvents = cms.untracked.PSet(
SelectEvents = cms.vstring('filter_step')
),
)
# Additional output definition
# Other statements
process.GlobalTag.globaltag = 'GR10_P_V12::All'
# filter on L1 bits
from HLTrigger.HLTfilters.hltLevel1GTSeed_cfi import hltLevel1GTSeed
process.bscSplash = hltLevel1GTSeed.clone(
L1TechTriggerSeeding = cms.bool(True),
L1SeedsLogicalExpression = cms.string('42 OR 43')
)
process.bscHighMult = hltLevel1GTSeed.clone(
L1TechTriggerSeeding = cms.bool(True),
L1SeedsLogicalExpression = cms.string('35')
)
process.bscOR = hltLevel1GTSeed.clone(
L1TechTriggerSeeding = cms.bool(True),
L1SeedsLogicalExpression = cms.string('(34) AND NOT (36 OR 37 OR 38 OR 39)')
)
process.bscAND = hltLevel1GTSeed.clone(
L1TechTriggerSeeding = cms.bool(True),
L1SeedsLogicalExpression = cms.string('40')
process.FEVT = cms.OutputModule("PoolOutputModule",
splitLevel = cms.untracked.int32(0),
outputCommands = process.RECOEventContent.outputCommands,
fileName = cms.untracked.string('rerecoOutput.root'),
dataset = cms.untracked.PSet(
dataTier = cms.untracked.string('RECO'),
filterName = cms.untracked.string('')
)
)
#======================================
# Trigger Filter
#======================================
process.load('L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff')
from HLTrigger.HLTfilters.hltLevel1GTSeed_cfi import hltLevel1GTSeed
process.bit40 = hltLevel1GTSeed.clone(
L1TechTriggerSeeding = cms.bool(True),
L1SeedsLogicalExpression = cms.string('0 AND (40 OR 41) AND NOT (36 OR 37 OR 38 OR 39)')
)
#======================================
# Services
#======================================
process.Timing = cms.Service("Timing")
#======================================
# Analyzer
#======================================
process.test = cms.EDAnalyzer("Analyzer",
VxInputTag = cms.InputTag("pixelVertices"),
AnalysisType = cms.string("Vertex"),
PrintMsg = cms.bool(False),
# '/store/data/Run2011A/MinimumBias/RECO/May10ReReco-v2/0002/BC5BF6CB-6682-E011-B5E2-00248C0BE01E.root',
# '/store/data/Run2011A/MinimumBias/RECO/May10ReReco-v2/0002/AA6C29E6-6F82-E011-8F12-001A92971B62.root',
# '/store/data/Run2011A/MinimumBias/RECO/May10ReReco-v2/0002/A476C008-8A82-E011-A676-002354EF3BDA.root'
# )
# )
# Configure the trigger
process.load(
'L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff')
from HLTrigger.HLTfilters.hltLevel1GTSeed_cfi import hltLevel1GTSeed
# L1 technical triggers:
# 0 -- beam crossing
# 36-39 -- beam halo
process.VetoHardInt = hltLevel1GTSeed.clone(
L1TechTriggerSeeding=cms.bool(True),
L1SeedsLogicalExpression=cms.string('0 AND NOT (36 OR 37 OR 38 OR 39)'))
# Remove the so-called "scraping" events
process.noScraping = cms.EDFilter(
"FilterOutScraping",
applyfilter=cms.untracked.bool(True),
debugOn=cms.untracked.bool(False), ## Or 'True' to get some per-event info
numtrack=cms.untracked.uint32(10),
thresh=cms.untracked.double(0.25))
# min bias filter
process.HLTZeroBias = cms.EDFilter(
"HLTHighLevel",
TriggerResultsTag=cms.InputTag("TriggerResults", "", "HLT"),
HLTPaths=cms.vstring('HLT_ZeroBias', 'HLT_ZeroBias_v1', 'HLT_ZeroBias_v2',
