|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
|
import FWCore.ParameterSet.Config as cms
looseMuonCut = " (isGlobalMuon || isTrackerMuon) && isPFMuon"
looseIsoCut = "(pfIsolationR04.sumChargedHadronPt + max(0., pfIsolationR04.sumNeutralHadronEt + pfIsolationR04.sumPhotonEt - 0.5 * pfIsolationR04.sumPUPt) ) / pt < 0.25"
tightMuonCut = " isGlobalMuon && isPFMuon && globalTrack.normalizedChi2 < 10. && globalTrack.hitPattern.numberOfValidMuonHits > 0 && " + \
"numberOfMatchedStations > 1 && innerTrack.hitPattern.numberOfValidPixelHits > 0 && innerTrack.hitPattern.trackerLayersWithMeasurement > 5"
# CB PV cut!
tightIsoCut = "(pfIsolationR04.sumChargedHadronPt + max(0., pfIsolationR04.sumNeutralHadronEt + pfIsolationR04.sumPhotonEt - 0.5 * pfIsolationR04.sumPUPt) ) / pt < 0.15"
EletightIsoCut = "(pfIsolationVariables.sumChargedHadronPt + max(0., pfIsolationVariables.sumNeutralHadronEt + pfIsolationVariables.sumPhotonEt - 0.5 * pfIsolationVariables.sumPUPt) ) / pt < 0.1"
ElelooseIsoCut = "(pfIsolationVariables.sumChargedHadronPt + max(0., pfIsolationVariables.sumNeutralHadronEt + pfIsolationVariables.sumPhotonEt - 0.5 * pfIsolationVariables.sumPUPt) ) / pt < 0.15"
looseElecCut = "((full5x5_sigmaIetaIeta < 0.011 && superCluster.isNonnull && superCluster.seed.isNonnull && (deltaEtaSuperClusterTrackAtVtx - superCluster.eta + superCluster.seed.eta) < 0.00477 && abs(deltaPhiSuperClusterTrackAtVtx) < 0.222 && hadronicOverEm < 0.298 && abs(1.0 - eSuperClusterOverP)*1.0/ecalEnergy < 0.241 && gsfTrack.hitPattern.numberOfHits('MISSING_INNER_HITS') <= 1 && abs(superCluster.eta) < 1.479) || (full5x5_sigmaIetaIeta() < 0.0314 && superCluster.isNonnull && superCluster.seed.isNonnull && (deltaEtaSuperClusterTrackAtVtx - superCluster.eta + superCluster.seed.eta) < 0.00868 && abs(deltaPhiSuperClusterTrackAtVtx) < 0.213 && hadronicOverEm < 0.101 && abs(1.0 - eSuperClusterOverP)*1.0/ecalEnergy < 0.14 && gsfTrack.hitPattern.numberOfHits('MISSING_INNER_HITS') <= 1 && abs(superCluster.eta) > 1.479))"
elecIPcut = "(abs(gsfTrack.d0)<0.05 & abs(gsfTrack.dz)<0.1 & abs(superCluster.eta) < 1.479)||(abs(gsfTrack.d0)<0.1 && abs(gsfTrack.dz)<0.2 && abs(superCluster.eta) > 1.479)"
tightElecCut = "((full5x5_sigmaIetaIeta < 0.00998 && superCluster.isNonnull && superCluster.seed.isNonnull && (deltaEtaSuperClusterTrackAtVtx - superCluster.eta + superCluster.seed.eta) < 0.00308 && abs(deltaPhiSuperClusterTrackAtVtx) < 0.0816 && hadronicOverEm < 0.0414 && abs(1.0 - eSuperClusterOverP)*1.0/ecalEnergy < 0.0129 && gsfTrack.hitPattern().numberOfLostHits('MISSING_INNER_HITS') <= 1 && abs(superCluster.eta) < 1.479) || (full5x5_sigmaIetaIeta() < 0.0292 && superCluster.isNonnull && superCluster.seed.isNonnull && (deltaEtaSuperClusterTrackAtVtx - superCluster.eta + superCluster.seed.eta) < 0.00605 && abs(deltaPhiSuperClusterTrackAtVtx) < 0.0394 && hadronicOverEm < 0.0641 && abs(1.0 - eSuperClusterOverP)*1.0/ecalEnergy < 0.0129 && gsfTrack.hitPattern().numberOfLostHits('MISSING_INNER_HITS') <= 1 && abs(superCluster.eta) > 1.479))"
from DQMServices.Core.DQMEDAnalyzer import DQMEDAnalyzer
topSingleMuonMediumDQM_miniAOD = DQMEDAnalyzer('TopSingleLeptonDQM_miniAOD',
setup = cms.PSet(
directory = cms.string("Physics/Top/TopSingleMuonMediumDQM_miniAOD/"),
sources = cms.PSet(
muons = cms.InputTag("slimmedMuons"),
elecs = cms.InputTag("slimmedElectrons"),
jets = cms.InputTag("slimmedJets"),
mets = cms.VInputTag("slimmedMETs", "slimmedMETsPuppi"),
pvs = cms.InputTag("offlineSlimmedPrimaryVertices")
),
monitoring = cms.PSet(
verbosity = cms.string("DEBUG")
),
pvExtras = cms.PSet(
select = cms.string("abs(z) < 24. & position.rho < 2. & ndof > 4 & !isFake")
),
elecExtras = cms.PSet(
select = cms.string(tightElecCut + "&& pt>20 & abs(eta)<2.5 & (abs(superCluster.eta) <= 1.4442 || abs(superCluster.eta) >= 1.5660)"),
rho = cms.InputTag("fixedGridRhoFastjetAll"),
#isolation = cms.string(ElelooseIsoCut),
),
muonExtras = cms.PSet(
select = cms.string(tightMuonCut + " && pt>20 & abs(eta)<2.4 && " + looseIsoCut),
isolation = cms.string(looseIsoCut)
),
jetExtras = cms.PSet(
select = cms.string("pt>30 & abs(eta)<2.4 "),
jetBTaggers = cms.PSet(
trackCountingEff = cms.PSet(
label = cms.InputTag("trackCountingHighEffBJetTags" ),
workingPoint = cms.double(1.25)
),
trackCountingPur = cms.PSet(
label = cms.InputTag("trackCountingHighPurBJetTags" ),
workingPoint = cms.double(3.00)
),
secondaryVertex = cms.PSet(
label = cms.InputTag("simpleSecondaryVertexHighEffBJetTags"),
workingPoint = cms.double(2.05)
),
cvsVertex = cms.PSet(
label = cms.InputTag("combinedSecondaryVertexBJetTags"),
workingPoint = cms.double(0.898)
)
),
),
massExtras = cms.PSet(
lowerEdge = cms.double( 70.),
upperEdge = cms.double(110.)
),
),
preselection = cms.PSet(
vertex = cms.PSet(
src = cms.InputTag("offlineSlimmedPrimaryVertices"),#,
select = cms.string("abs(z) < 24. & position.rho < 2. & ndof > 4 & !isFake")
)
),
selection = cms.VPSet(
cms.PSet(
label = cms.string("muons:step0"),
src = cms.InputTag("slimmedMuons"),
select = cms.string(tightMuonCut + " && pt>20 & abs(eta)<2.4 && " + looseIsoCut), #tightMuonCut +"&&"+ tightIsoCut + " && pt>20 & abs(eta)<2.1"), # CB what about iso? CD Added tightIso
min = cms.int32(1),
#max = cms.int32(1),
),
cms.PSet(
label = cms.string("jets:step1"),
src = cms.InputTag("slimmedJets"),
select = cms.string("pt>30 & abs(eta)<2.4 "),
min = cms.int32(4),
),
cms.PSet(
label = cms.string("met:step2"),
src = cms.InputTag("slimmedMETs"),
select = cms.string("pt>30"),
#min = cms.int32(2),
),
# cms.PSet(
# label = cms.string("jets:step3"),
# src = cms.InputTag("slimmedJets"),
# select = cms.string("pt>30 & abs(eta)<2.5 "),
# min = cms.int32(3),
# ),
# cms.PSet(
# label = cms.string("jets:step4"),
# src = cms.InputTag("slimmedJets"),
# select = cms.string("pt>30 & abs(eta)<2.5 "),
# min = cms.int32(4),
# ),
)
)
topSingleElectronMediumDQM_miniAOD = DQMEDAnalyzer('TopSingleLeptonDQM_miniAOD',
setup = cms.PSet(
directory = cms.string("Physics/Top/TopSingleElectronMediumDQM_miniAOD/"),
sources = cms.PSet(
muons = cms.InputTag("slimmedMuons"),
elecs = cms.InputTag("slimmedElectrons"),
jets = cms.InputTag("slimmedJets"),
mets = cms.VInputTag("slimmedMETs", "slimmedMETsPuppi"),
pvs = cms.InputTag("offlineSlimmedPrimaryVertices")
),
monitoring = cms.PSet(
verbosity = cms.string("DEBUG")
),
pvExtras = cms.PSet(
select = cms.string("abs(z) < 24. & position.rho < 2. & ndof > 4 & !isFake")
),
elecExtras = cms.PSet(
select = cms.string(tightElecCut + " && pt>20 & abs(eta)<2.5 & (abs(superCluster.eta) <= 1.4442 || abs(superCluster.eta) >= 1.5660)"),
#select = cms.string(looseElecCut+ "&& pt>20 & abs(eta)<2.5 & (abs(superCluster.eta) <= 1.4442 || abs(superCluster.eta) >= 1.5660)"),
rho = cms.InputTag("fixedGridRhoFastjetAll"),
#isolation = cms.string(ElelooseIsoCut),
),
muonExtras = cms.PSet(
select = cms.string(tightMuonCut + " && pt>20 & abs(eta)<2.4 && " + looseIsoCut),
isolation = cms.string(looseIsoCut)
),
jetExtras = cms.PSet(
select = cms.string("pt>30 & abs(eta)<2.4 "),
jetBTaggers = cms.PSet(
trackCountingEff = cms.PSet(
label = cms.InputTag("trackCountingHighEffBJetTags" ),
workingPoint = cms.double(1.25)
),
trackCountingPur = cms.PSet(
label = cms.InputTag("trackCountingHighPurBJetTags" ),
workingPoint = cms.double(3.00)
),
secondaryVertex = cms.PSet(
label = cms.InputTag("simpleSecondaryVertexHighEffBJetTags"),
workingPoint = cms.double(2.05)
),
cvsVertex = cms.PSet(
label = cms.InputTag("combinedSecondaryVertexBJetTags"),
workingPoint = cms.double(0.898)
)
),
),
massExtras = cms.PSet(
lowerEdge = cms.double( 70.),
upperEdge = cms.double(110.)
),
),
preselection = cms.PSet(
vertex = cms.PSet(
src = cms.InputTag("offlineSlimmedPrimaryVertices"),#,
select = cms.string("abs(z) < 24. & position.rho < 2. & ndof > 4 & !isFake")
)
),
selection = cms.VPSet(
cms.PSet(
label = cms.string("elecs:step0"),
src = cms.InputTag("slimmedElectrons"),
select = cms.string("pt>20 & abs(eta)<2.5 & (abs(superCluster.eta) <= 1.4442 || abs(superCluster.eta) >= 1.5660) &" + tightElecCut),
min = cms.int32(1),
max = cms.int32(1),
),
cms.PSet(
label = cms.string("jets:step1"),
src = cms.InputTag("slimmedJets"),
select = cms.string("pt>30 & abs(eta)<2.4 "),
min = cms.int32(4),
),
cms.PSet(
label = cms.string("met:step2"),
src = cms.InputTag("slimmedMETs"),
select = cms.string("pt>30"),
#min = cms.int32(2),
),
#cms.PSet(
#label = cms.string("jets:step3"),
#src = cms.InputTag("slimmedJets"),
#select = cms.string("pt>30 & abs(eta)<2.5 "),
#min = cms.int32(3),
#),
#cms.PSet(
#label = cms.string("jets:step4"),
#src = cms.InputTag("slimmedJets"),
#select = cms.string("pt>30 & abs(eta)<2.5 "),
#min = cms.int32(4),
#),
)
)
|