Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-02-14 12:44:22

0001 import FWCore.ParameterSet.Config as cms
0002 
0003 
0004 
0005 ##
0006 ## ALCARECOTkAlMuonIsolated selection
0007 ##
0008 
0009 ## First select goodId + isolated muons
0010 TkAlGoodIdMuonSelector = cms.EDFilter("MuonSelector",
0011     src = cms.InputTag('muons'),
0012     cut = cms.string('isGlobalMuon &'
0013                      'isTrackerMuon &'
0014                      'numberOfMatches > 1 &'
0015                      'globalTrack.hitPattern.numberOfValidMuonHits > 0 &'
0016                      'abs(eta) < 2.5 &'
0017                      'globalTrack.normalizedChi2 < 20.'),
0018     filter = cms.bool(True)
0019 )
0020 TkAlRelCombIsoMuonSelector = cms.EDFilter("MuonSelector",
0021     src = cms.InputTag(''),
0022     cut = cms.string('(isolationR03().sumPt + isolationR03().emEt + isolationR03().hadEt)/pt  < 0.15'),
0023     filter = cms.bool(True)
0024 )
0025 ALCARECOTkAlMuonIsolatedGoodMuons = TkAlGoodIdMuonSelector.clone()
0026 ALCARECOTkAlMuonIsolatedRelCombIsoMuons = TkAlRelCombIsoMuonSelector.clone(src = 'ALCARECOTkAlMuonIsolatedGoodMuons')
0027 
0028 ## Then select their tracks with additional cuts
0029 import Alignment.CommonAlignmentProducer.AlignmentTrackSelector_cfi
0030 ALCARECOTkAlMuonIsolated = Alignment.CommonAlignmentProducer.AlignmentTrackSelector_cfi.AlignmentTrackSelector.clone(
0031     filter = True, ##do not store empty events
0032     applyBasicCuts = True,
0033     ptMin = 2.0, ##GeV 
0034     etaMin = -3.5,
0035     etaMax = 3.5,
0036     nHitMin = 0
0037 )
0038 ALCARECOTkAlMuonIsolated.GlobalSelector.muonSource = 'ALCARECOTkAlMuonIsolatedRelCombIsoMuons'
0039 # Isolation is shifted to the muon preselection, and then applied intrinsically if applyGlobalMuonFilter = True
0040 ALCARECOTkAlMuonIsolated.GlobalSelector.applyIsolationtest = False
0041 ALCARECOTkAlMuonIsolated.GlobalSelector.applyGlobalMuonFilter = True
0042 ALCARECOTkAlMuonIsolated.TwoBodyDecaySelector.applyMassrangeFilter = False
0043 ALCARECOTkAlMuonIsolated.TwoBodyDecaySelector.applyChargeFilter = False
0044 ALCARECOTkAlMuonIsolated.TwoBodyDecaySelector.applyAcoplanarityFilter = False
0045 
0046 ## Define ALCARECO sequence
0047 mySeqALCARECOTkAlMuonIsolated = cms.Sequence(ALCARECOTkAlMuonIsolatedGoodMuons*ALCARECOTkAlMuonIsolatedRelCombIsoMuons*ALCARECOTkAlMuonIsolated)
0048 
0049 
0050 
0051 ##
0052 ## Good Primary Vertex Selection
0053 ##
0054 goodPVs = cms.EDFilter("VertexSelector",
0055     src = cms.InputTag('offlinePrimaryVertices'),
0056     cut = cms.string('ndof>4 &'
0057                      'abs(z)<24 &'
0058                      '!isFake &'
0059                      'position.Rho<2'
0060     ),
0061 )
0062 oneGoodPVSelection = cms.EDFilter("VertexCountFilter",
0063     src = cms.InputTag('goodPVs'),
0064     minNumber = cms.uint32(1),
0065     maxNumber = cms.uint32(99999),
0066     
0067 )
0068 seqVertexSelection = cms.Sequence(goodPVs*oneGoodPVSelection)
0069 
0070 
0071 
0072 
0073 
0074 
0075