Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2024-02-10 03:03:33

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