File indexing completed on 2023-03-17 10:53:44
0001 '''
0002 Defines the selection sequence ZmmgSkimSeq for the Zmmg skim for the
0003 RAW-RECO event content. It also defines several other modules and sequences
0004 used:
0005 ZmmgHLTFilter
0006 ZmmgTrailingMuons
0007 ZmmgLeadingMuons
0008 ZmmgDimuons
0009 ZmmgDimuonFilter
0010 ZmmgDimuonSequence
0011 ZmmgMergedSuperClusters
0012 ZmmgPhotonCandidates
0013 ZmmgPhotons
0014 ZmmgPhotonSequence
0015 ZmmgCandidates
0016 ZmmgFilter
0017 ZmmgSequence
0018
0019 Jan Veverka, Caltech, 5 May 2012
0020 '''
0021
0022 import copy
0023 import FWCore.ParameterSet.Config as cms
0024
0025
0026
0027
0028
0029
0030
0031 import copy
0032 from HLTrigger.HLTfilters.hltHighLevel_cfi import *
0033 ZmmgHLTFilter = copy.deepcopy(hltHighLevel)
0034 ZmmgHLTFilter.throw = cms.bool(False)
0035 ZmmgHLTFilter.HLTPaths = ['HLT_Mu*','HLT_IsoMu*','HLT_DoubleMu*']
0036
0037
0038
0039
0040
0041
0042
0043
0044 ZmmgTrailingMuons = cms.EDFilter('MuonSelector',
0045 src = cms.InputTag('muons'),
0046 cut = cms.string('''pt > 10 &&
0047 abs(eta) < 2.4 &&
0048 isGlobalMuon = 1 &&
0049 isTrackerMuon = 1 &&
0050 abs(innerTrack().dxy)<2.0'''),
0051 filter = cms.bool(True)
0052 )
0053
0054
0055 ZmmgLeadingMuons = cms.EDFilter('MuonSelector',
0056 src = cms.InputTag('ZmmgTrailingMuons'),
0057 cut = cms.string('pt > 20'),
0058 filter = cms.bool(True)
0059 )
0060
0061
0062 ZmmgDimuons = cms.EDProducer('CandViewShallowCloneCombiner',
0063 decay = cms.string('ZmmgLeadingMuons@+ ZmmgTrailingMuons@-'),
0064 checkCharge = cms.bool(True),
0065 cut = cms.string('mass > 30'),
0066 )
0067
0068
0069 ZmmgDimuonFilter = cms.EDFilter('CandViewCountFilter',
0070 src = cms.InputTag('ZmmgDimuons'),
0071 minNumber = cms.uint32(1)
0072 )
0073
0074
0075 ZmmgDimuonSequence = cms.Sequence(
0076 ZmmgTrailingMuons *
0077 ZmmgLeadingMuons *
0078 ZmmgDimuons *
0079 ZmmgDimuonFilter
0080 )
0081
0082
0083
0084
0085
0086
0087
0088
0089 ZmmgMergedSuperClusters = cms.EDProducer('EgammaSuperClusterMerger',
0090 src = cms.VInputTag(
0091 cms.InputTag('correctedHybridSuperClusters'),
0092 cms.InputTag('correctedMulti5x5SuperClustersWithPreshower')
0093 )
0094 )
0095
0096
0097 ZmmgPhotonCandidates = cms.EDProducer('ConcreteEcalCandidateProducer',
0098 src = cms.InputTag('ZmmgMergedSuperClusters'),
0099 particleType = cms.string('gamma')
0100 )
0101
0102
0103 ZmmgPhotons = cms.EDFilter('CandViewSelector',
0104 src = cms.InputTag('ZmmgPhotonCandidates'),
0105 cut = cms.string('et > 5'),
0106 filter = cms.bool(True)
0107 )
0108
0109
0110 ZmmgPhotonSequence = cms.Sequence(
0111 ZmmgMergedSuperClusters *
0112 ZmmgPhotonCandidates *
0113 ZmmgPhotons
0114 )
0115
0116
0117
0118
0119
0120
0121
0122
0123
0124
0125
0126
0127
0128
0129
0130
0131 ZmmgCandidates = cms.EDProducer('CandViewShallowCloneCombiner',
0132 decay = cms.string('ZmmgDimuons ZmmgPhotons'),
0133 checkCharge = cms.bool(False),
0134 cut = cms.string('''
0135 daughter(0).daughter(1).pt + daughter(1).pt > 20 &
0136 min(deltaR(daughter(0).daughter(0).eta,
0137 daughter(0).daughter(0).phi,
0138 daughter(1).eta,
0139 daughter(1).phi),
0140 deltaR(daughter(0).daughter(1).eta,
0141 daughter(0).daughter(1).phi,
0142 daughter(1).eta,
0143 daughter(1).phi)) < 1.5 &
0144 mass + daughter(0).mass < 200 &
0145 mass > 40
0146 '''),
0147 )
0148
0149
0150 ZmmgFilter = cms.EDFilter('CandViewCountFilter',
0151 src = cms.InputTag('ZmmgCandidates'),
0152 minNumber = cms.uint32(1)
0153 )
0154
0155 ZmmgSequence = cms.Sequence(
0156 ZmmgCandidates *
0157 ZmmgFilter
0158 )
0159
0160
0161
0162
0163
0164
0165
0166 ZmmgSkimSeq = cms.Sequence(
0167 ZmmgHLTFilter *
0168 ZmmgDimuonSequence *
0169 ZmmgPhotonSequence *
0170 ZmmgSequence
0171 )
0172