'/store/data/Commissioning08/Cosmics/RAW-RECO/CRAFT_ALL_V9_SuperPointing_225-v3/0015/3014AE2E-6503-DE11-B093-003048767DCD.root',
'/store/data/Commissioning08/Cosmics/RAW-RECO/CRAFT_ALL_V9_SuperPointing_225-v3/0012/EA27ED04-0602-DE11-B31E-001A92971B8C.root'
),
)
process.maxEvents = cms.untracked.PSet(input=cms.untracked.int32(100))
process.load("Configuration.StandardSequences.Geometry_cff")
process.load(
"Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
process.load("Configuration.StandardSequences.MagneticField_cff")
process.GlobalTag.globaltag = cms.string('CRAFT_ALL_V9::All')
from L1Trigger.L1ExtraFromDigis.l1extraParticles_cfi import l1extraParticles
process.l1muonsAnyBX = l1extraParticles.clone(
#muonSource = cms.InputTag( "hltGtDigis" ),
produceCaloParticles=False, ### we don't have digis for these
centralBxOnly=False ### this is the important point
)
### one could also convert STA track to a Candidate instead of using the reco::Muon,
### but one needs to add the dictionary for std::vector<reco::RecoStandAloneMuonCandidate> and it's edm::Wrapper
### in DataFormats/RecoCandidate.
### Alternatively, you can make your own reco::Muons...
#process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi") # to get the muon mass
#process.lhcBasicMuons = cms.EDProducer("ConcreteStandAloneMuonCandidateProducer", src = cms.InputTag("lhcStandAloneMuonsBarrelOnly"), particleType = cms.string('mu+'))
#process.lhcBasicMuonsUAV = process.lhcBasicMuons.clone(src = cms.InputTag("lhcStandAloneMuonsBarrelOnly", "UpdatedAtVtx"))
# Make l1extraParticles for the three bunch crossing L1A=0,+1,-1
from L1Trigger.L1ExtraFromDigis.l1extraParticles_cfi import l1extraParticles
process.l1muonsAnyBX = l1extraParticles.clone(
produceCaloParticles=False, ### we don't have digis for these
## Global Muon Filter
process.globalMuFilter = cms.EDFilter("TrackCountFilter",
src = cms.InputTag("globalMuons"),
minNumber = cms.uint32(1)
)
############################################
## -0- Prefilter
process.preFilter = cms.Sequence(process.JPsiHLT *
process.noScraping *
process.primaryVertexFilter +
process.globalMuFilter)
############################################
## -1- l1muonsAnyBX - from L1extraParticles
from L1Trigger.L1ExtraFromDigis.l1extraParticles_cfi import l1extraParticles
process.l1muonsAnyBX = l1extraParticles.clone(
centralBxOnly = cms.bool(True)
)
############################################
## -2- patMuonsWithTrigger
process.load("MuonAnalysis.MuonAssociators.patMuonsWithTrigger_cff")
############################################
## -3- tagMuons - from patMuonsWithTrigger, Tight muon https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuideMuonId
process.tagMuons = cms.EDFilter("PATMuonRefSelector",
src = cms.InputTag("patMuonsWithTrigger"),
cut = cms.string('isGlobalMuon && isTrackerMuon'\
+'&& abs(dB)<0.2 && globalTrack().normalizedChi2<10'\
+'&& globalTrack().hitPattern().numberOfValidTrackerHits>10'\
+'&& globalTrack().hitPattern().numberOfValidMuonHits>6'\
+'&& numberOfMatchedStations()>1'\
+'&& globalTrack().hitPattern().numberOfValidPixelHits>0'\
+'&& (chargedHadronIso()+neutralHadronIso()+photonIso())<0.1*pt'\
fileNames = cms.untracked.vstring(
'/store/data/Commissioning08/Cosmics/RAW-RECO/CRAFT_ALL_V9_SuperPointing_225-v3/0015/3014AE2E-6503-DE11-B093-003048767DCD.root',
'/store/data/Commissioning08/Cosmics/RAW-RECO/CRAFT_ALL_V9_SuperPointing_225-v3/0012/EA27ED04-0602-DE11-B31E-001A92971B8C.root'
),
)
process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(100) )
process.load("Configuration.StandardSequences.Geometry_cff")
process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
process.load("Configuration.StandardSequences.MagneticField_cff")
process.GlobalTag.globaltag = cms.string('CRAFT_ALL_V9::All')
from L1Trigger.L1ExtraFromDigis.l1extraParticles_cfi import l1extraParticles
process.l1muonsAnyBX = l1extraParticles.clone(
#muonSource = cms.InputTag( "hltGtDigis" ),
produceCaloParticles = False, ### we don't have digis for these
centralBxOnly = False ### this is the important point
)
### one could also convert STA track to a Candidate instead of using the reco::Muon,
### Alternatively, you can make your own reco::Muons...
#process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi") # to get the muon mass
#process.lhcBasicMuons = cms.EDProducer("ConcreteChargedCandidateProducer", src = cms.InputTag("..."), particleType = cms.string('mu+'))
#process.lhcBasicMuonsUAV = process.lhcBasicMuons.clone(src = cms.InputTag("lhcStandAloneMuonsBarrelOnly", "UpdatedAtVtx"))
### It would be more correct to use 'ConcreteStandAloneMuonCandidateProducer' (for which the track appears as 'muon' instead of 'tracker',
### but there are no dictionaries for it (see comment in testCRAFT.py); anyway, when using 'cosmicPropagationHypothesis' it doesn't matter
# Make l1extraParticles for the three bunch crossing L1A=0,+1,-1
from L1Trigger.L1ExtraFromDigis.l1extraParticles_cfi import l1extraParticles
process.l1muonsAnyBX = l1extraParticles.clone(
produceCaloParticles = False, ### we don't have digis for these
