cut = cms.string("pt > %f" % (ptMin,) ), ) process.slimAOD = cms.Sequence( process.genMuons + process.goodTracks ) ### ==== Make PAT Muons ==== from PhysicsTools.PatAlgos.producersLayer1.muonProducer_cfi import allLayer1Muons process.patMuonsWithoutTrigger = allLayer1Muons.clone( # embed the tracks, so we don't have to carry them around embedTrack = True, embedCombinedMuon = True, embedStandAloneMuon = True, # then switch off some features we don't need embedPickyMuon = False, embedTpfmsMuon = False, isolation = cms.PSet(), # no extra isolation beyond what's in reco::Muon itself isoDeposits = cms.PSet(), # no heavy isodeposits addGenMatch = False, # no mc: T&P doesn't take it from here anyway. ) ### ==== Unpack trigger, and match ==== process.load("PhysicsTools.PatAlgos.triggerLayer1.triggerProducer_cfi") process.patTrigger.onlyStandAlone = True process.patTrigger.processName = triggerProcess from PhysicsTools.PatAlgos.triggerLayer1.triggerMatcher_cfi import muonTriggerMatchHLT1MuonIso process.muonMatchHLTMuX = muonTriggerMatchHLT1MuonIso.clone( src = 'patMuonsWithoutTrigger', pathNames = [ triggerPath ] )
"TrackSelector", src=cms.InputTag("generalTracks"), cut=cms.string("pt > %f" % (ptMin, )), ) process.slimAOD = cms.Sequence(process.genMuons + process.goodTracks) ### ==== Make PAT Muons ==== from PhysicsTools.PatAlgos.producersLayer1.muonProducer_cfi import allLayer1Muons process.patMuonsWithoutTrigger = allLayer1Muons.clone( # embed the tracks, so we don't have to carry them around embedTrack=True, embedCombinedMuon=True, embedStandAloneMuon=True, # then switch off some features we don't need embedPickyMuon=False, embedTpfmsMuon=False, isolation=cms.PSet( ), # no extra isolation beyond what's in reco::Muon itself isoDeposits=cms.PSet(), # no heavy isodeposits addGenMatch=False, # no mc: T&P doesn't take it from here anyway. ) ### ==== Unpack trigger, and match ==== process.load("PhysicsTools.PatAlgos.triggerLayer1.triggerProducer_cfi") process.patTrigger.onlyStandAlone = True process.patTrigger.processName = triggerProcess from PhysicsTools.PatAlgos.triggerLayer1.triggerMatcher_cfi import muonTriggerMatchHLT1MuonIso process.muonMatchHLTMuX = muonTriggerMatchHLT1MuonIso.clone( src='patMuonsWithoutTrigger', pathNames=[triggerPath])
process.muonL1Match.maxDeltaR = 9999.0 # no deltaR match process.muonL1Match.preselection = "" process.muonL1Match.writeExtraInfo = True #process.muonL1MatchUAV = process.muonL1Match.clone(src = 'lhcBasicMuonsUAV') #### Merge matches into pat::GenericParticle ## from PhysicsTools.PatAlgos.producersLayer1.genericParticleProducer_cfi import allLayer1GenericParticles ## process.lhcMuons = allLayer1GenericParticles.clone( ## src = "lhcBasicMuons", ... #### Merge matches into pat::Muon from PhysicsTools.PatAlgos.producersLayer1.muonProducer_cfi import allLayer1Muons process.lhcMuons = allLayer1Muons.clone( muonSource='lhcSTAMuonsBarrelOnly', # src = 'lhcBasicMuons' embedStandAloneMuon=True, addTrigMatch=True, isolation=cms.PSet(), isoDeposits=cms.PSet(), addGenMatch=False, addTeVRefits=False ## turn off unwanted pat::Muon features ) #process.lhcMuonsUAV = process.lhcMuons.clone(src = 'lhcBasicMuonsUAV') def addMatch(muons, matcher): muons.trigPrimMatch = cms.VInputTag( cms.InputTag(matcher), cms.InputTag(matcher, "propagatedReco")) muons.userData.userInts.src = cms.VInputTag( cms.InputTag(matcher, "bx"), cms.InputTag(matcher, "quality"), cms.InputTag(matcher, "isolated")) muons.userData.userFloats.src = cms.VInputTag( cms.InputTag(matcher, "deltaR"))
process.muonL1Match.useState = 'innermost' process.muonL1Match.maxDeltaPhi = 50.0*3.14/180.0 # 50 degrees, slide 3 of http://indico.cern.ch/getFile.py/access?contribId=2&resId=0&materialId=slides&confId=57817 process.muonL1Match.maxDeltaR = 9999.0 # no deltaR match process.muonL1Match.preselection = "" process.muonL1Match.writeExtraInfo = True process.muonL1Match.cosmicPropagationHypothesis = cms.bool(True) # must specify the type, as I'm adding a new parameter and not replacing one #### Merge matches into pat::GenericParticle ## from PhysicsTools.PatAlgos.producersLayer1.genericParticleProducer_cfi import allLayer1GenericParticles ## process.myMuons = allLayer1GenericParticles.clone( ## src = "lhcBasicMuons", ... #### Merge matches into pat::Muon from PhysicsTools.PatAlgos.producersLayer1.muonProducer_cfi import allLayer1Muons process.myMuons = allLayer1Muons.clone( muonSource = 'muons', # src = 'lhcBasicMuons' embedStandAloneMuon = True, addTrigMatch = True, isolation = cms.PSet(), isoDeposits = cms.PSet(), addGenMatch = False, addTeVRefits = False ## turn off unwanted pat::Muon features ) process.myMuons.trigPrimMatch = cms.VInputTag( cms.InputTag("muonL1Match"), cms.InputTag("muonL1Match","propagatedReco") ) process.myMuons.userData.userInts.src = cms.VInputTag( cms.InputTag("muonL1Match", "bx"), cms.InputTag("muonL1Match", "quality"), cms.InputTag("muonL1Match", "isolated") ) process.myMuons.userData.userFloats.src = cms.VInputTag( cms.InputTag("muonL1Match", "deltaR") ) process.p = cms.Path( process.l1muonsAnyBX * process.muonL1Match * process.myMuons ) process.out = cms.OutputModule("PoolOutputModule", fileName = cms.untracked.string('l1CosmicMuonsCRAFT.root'), outputCommands = cms.untracked.vstring(