CMS 3D CMS Logo

topSingleLeptonDQM_miniAOD_cfi.py
Go to the documentation of this file.
1 import FWCore.ParameterSet.Config as cms
2 
3 looseMuonCut = " (isGlobalMuon || isTrackerMuon) && isPFMuon"
4 looseIsoCut = "(pfIsolationR04.sumChargedHadronPt + max(0., pfIsolationR04.sumNeutralHadronEt + pfIsolationR04.sumPhotonEt - 0.5 * pfIsolationR04.sumPUPt) ) / pt < 0.2"
5 
6 tightMuonCut = " isGlobalMuon && isPFMuon && globalTrack.normalizedChi2 < 10. && globalTrack.hitPattern.numberOfValidMuonHits > 0 && " + \
7  "numberOfMatchedStations > 1 && innerTrack.hitPattern.numberOfValidPixelHits > 0 && innerTrack.hitPattern.trackerLayersWithMeasurement > 8"
8  # CB PV cut!
9 tightIsoCut = "(pfIsolationR04.sumChargedHadronPt + max(0., pfIsolationR04.sumNeutralHadronEt + pfIsolationR04.sumPhotonEt - 0.5 * pfIsolationR04.sumPUPt) ) / pt < 0.12"
10 
11 EletightIsoCut = "(pfIsolationVariables.sumChargedHadronPt + max(0., pfIsolationVariables.sumNeutralHadronEt + pfIsolationVariables.sumPhotonEt - 0.5 * pfIsolationVariables.sumPUPt) ) / pt < 0.1"
12 ElelooseIsoCut = "(pfIsolationVariables.sumChargedHadronPt + max(0., pfIsolationVariables.sumNeutralHadronEt + pfIsolationVariables.sumPhotonEt - 0.5 * pfIsolationVariables.sumPUPt) ) / pt < 0.15"
13 
14 topSingleLeptonDQM_miniAOD = cms.EDAnalyzer("TopSingleLeptonDQM_miniAOD",
15  setup = cms.PSet(
16  directory = cms.string("Physics/Top/TopSingleLeptonDQM_miniAOD/"),
17  sources = cms.PSet(
18 
19  elecs = cms.InputTag("slimmedElectrons"),
20  jets = cms.InputTag("slimmedJets"),
21  mets = cms.VInputTag("slimmedMETs", "slimmedMETsNoHF", "slimmedMETsPuppi"),
22 
23  pvs = cms.InputTag("offlineSlimmedPrimaryVertices")
24  ),
25  monitoring = cms.PSet(
26  verbosity = cms.string("DEBUG")
27  ),
28  pvExtras = cms.PSet(
29  ),
30 
31  elecExtras = cms.PSet(
32  select = cms.string("pt>15 & abs(eta)<2.5 & abs(gsfTrack.d0)<1 & abs(gsfTrack.dz)<20"),
33  isolation = cms.string(ElelooseIsoCut),
34  ),
35  muonExtras = cms.PSet(
36  select = cms.string(looseMuonCut + " && pt>10 & abs(eta)<2.4"),
37  isolation = cms.string(looseIsoCut),
38  ),
39  jetExtras = cms.PSet(
40  select = cms.string("pt>30 & abs(eta)<2.5 "),
41  ),
42  massExtras = cms.PSet(
43  lowerEdge = cms.double( 70.),
44  upperEdge = cms.double(110.)
45  ),
46  triggerExtras = cms.PSet(
47  src = cms.InputTag("TriggerResults","","HLT"),
48  paths = cms.vstring(['HLT_Mu3:HLT_QuadJet15U',
49  'HLT_Mu5:HLT_QuadJet15U',
50  'HLT_Mu7:HLT_QuadJet15U',
51  'HLT_Mu9:HLT_QuadJet15U'])
52  )
53  ),
54  preselection = cms.PSet(
55 
56  vertex = cms.PSet(
57  src = cms.InputTag("offlineSlimmedPrimaryVertices")#,
58  )
59  ),
60  selection = cms.VPSet(
61  cms.PSet(
62  label = cms.string("jets:step0"),
63  src = cms.InputTag("slimmedJets"),
64  select = cms.string("pt>30 & abs(eta)<2.5 "),
65  min = cms.int32(2),
66  ),
67  )
68 )
69 
70 topSingleMuonLooseDQM_miniAOD = cms.EDAnalyzer("TopSingleLeptonDQM_miniAOD",
71 
72  setup = cms.PSet(
73  directory = cms.string("Physics/Top/TopSingleMuonLooseDQM_miniAOD/"),
74  sources = cms.PSet(
75  muons = cms.InputTag("slimmedMuons"),
76  elecs = cms.InputTag("slimmedElectrons"),
77  jets = cms.InputTag("slimmedJets"),
78  mets = cms.VInputTag("slimmedMETs", "slimmedMETsNoHF", "slimmedMETsPuppi"),
79  pvs = cms.InputTag("offlineSlimmedPrimaryVertices")
80  ),
81  monitoring = cms.PSet(
82  verbosity = cms.string("DEBUG")
83  ),
84  pvExtras = cms.PSet(
85  ),
86  muonExtras = cms.PSet(
87  select = cms.string(looseMuonCut + " && pt > 10 & abs(eta)<2.4"),
88  isolation = cms.string(looseIsoCut)
89  ),
90  jetExtras = cms.PSet(
91  select = cms.string("pt>30 & abs(eta)<2.5"),
92  jetBTaggers = cms.PSet(
93  trackCountingEff = cms.PSet(
94  label = cms.InputTag("trackCountingHighEffBJetTags" ),
95  workingPoint = cms.double(1.25)
96  ),
97  trackCountingPur = cms.PSet(
98  label = cms.InputTag("trackCountingHighPurBJetTags" ),
99  workingPoint = cms.double(3.00)
100  ),
101  secondaryVertex = cms.PSet(
102  label = cms.InputTag("simpleSecondaryVertexHighEffBJetTags"),
103  workingPoint = cms.double(2.05)
104  ),
105  cvsVertex = cms.PSet(
106  label = cms.InputTag("combinedSecondaryVertexBJetTags"),
107  workingPoint = cms.double(0.898)
108  )
109  ),
110  ),
111  massExtras = cms.PSet(
112  lowerEdge = cms.double( 70.),
113  upperEdge = cms.double(110.)
114  ),
115 
116  triggerExtras = cms.PSet(
117  src = cms.InputTag("TriggerResults","","HLT"),
118  paths = cms.vstring(['HLT_Mu3:HLT_QuadJet15U',
119  'HLT_Mu5:HLT_QuadJet15U',
120  'HLT_Mu7:HLT_QuadJet15U',
121  'HLT_Mu9:HLT_QuadJet15U',
122  'HLT_Mu11:HLT_QuadJet15U'])
123  )
124  ),
125  preselection = cms.PSet(
126  vertex = cms.PSet(
127  src = cms.InputTag("offlineSlimmedPrimaryVertices")#,
128  )
129  ),
130  selection = cms.VPSet(
131  cms.PSet(
132  label = cms.string("muons:step0"),
133  src = cms.InputTag("slimmedMuons"),
134  select = cms.string(looseMuonCut + looseIsoCut + " && pt>10 & abs(eta)<2.4"), # CB what about iso? CD Added looseIso
135  min = cms.int32(1),
136  ),
137  cms.PSet(
138  label = cms.string("jets:step1"),
139  src = cms.InputTag("slimmedJets"),
140  select = cms.string("pt>30 & abs(eta)<2.5 "),
141  min = cms.int32(1),
142  ),
143  cms.PSet(
144  label = cms.string("jets:step2"),
145  src = cms.InputTag("slimmedJets"),
146  select = cms.string("pt>30 & abs(eta)<2.5 "),
147  min = cms.int32(2),
148  ),
149  cms.PSet(
150  label = cms.string("jets:step3"),
151  src = cms.InputTag("slimmedJets"),
152  select = cms.string("pt>30 & abs(eta)<2.5 "),
153  min = cms.int32(3),
154  ),
155  cms.PSet(
156  label = cms.string("jets:step4"),
157  src = cms.InputTag("slimmedJets"),
158  select = cms.string("pt>30 & abs(eta)<2.5 "),
159  min = cms.int32(4),
160  ),
161  )
162 )
163 topSingleMuonMediumDQM_miniAOD = cms.EDAnalyzer("TopSingleLeptonDQM_miniAOD",
164  setup = cms.PSet(
165  directory = cms.string("Physics/Top/TopSingleMuonMediumDQM_miniAOD/"),
166  sources = cms.PSet(
167  muons = cms.InputTag("slimmedMuons"),
168  elecs = cms.InputTag("slimmedElectrons"),
169  jets = cms.InputTag("slimmedJets"),
170  mets = cms.VInputTag("slimmedMETs", "slimmedMETsNoHF", "slimmedMETsPuppi"),
171  pvs = cms.InputTag("offlineSlimmedPrimaryVertices")
172 
173  ),
174  monitoring = cms.PSet(
175  verbosity = cms.string("DEBUG")
176  ),
177  pvExtras = cms.PSet(
178  ),
179  muonExtras = cms.PSet(
180  select = cms.string(looseMuonCut + " && pt>20 & abs(eta)<2.1"),
181  isolation = cms.string(looseIsoCut)
182  ),
183  jetExtras = cms.PSet(
184  select = cms.string("pt>30 & abs(eta)<2.5 "),
185  jetBTaggers = cms.PSet(
186  trackCountingEff = cms.PSet(
187  label = cms.InputTag("trackCountingHighEffBJetTags" ),
188  workingPoint = cms.double(1.25)
189  ),
190  trackCountingPur = cms.PSet(
191  label = cms.InputTag("trackCountingHighPurBJetTags" ),
192  workingPoint = cms.double(3.00)
193  ),
194  secondaryVertex = cms.PSet(
195  label = cms.InputTag("simpleSecondaryVertexHighEffBJetTags"),
196  workingPoint = cms.double(2.05)
197  ),
198  cvsVertex = cms.PSet(
199  label = cms.InputTag("combinedSecondaryVertexBJetTags"),
200  workingPoint = cms.double(0.898)
201  )
202  ),
203  ),
204  massExtras = cms.PSet(
205  lowerEdge = cms.double( 70.),
206  upperEdge = cms.double(110.)
207  ),
208  ),
209  preselection = cms.PSet(
210  vertex = cms.PSet(
211  src = cms.InputTag("offlineSlimmedPrimaryVertices")#,
212  )
213  ),
214  selection = cms.VPSet(
215  cms.PSet(
216  label = cms.string("muons:step0"),
217  src = cms.InputTag("slimmedMuons"),
218  select = cms.string(tightMuonCut +"&&"+ tightIsoCut + " && pt>20 & abs(eta)<2.1"), # CB what about iso? CD Added tightIso
219  min = cms.int32(1),
220  max = cms.int32(1),
221  ),
222  cms.PSet(
223  label = cms.string("jets:step1"),
224  src = cms.InputTag("slimmedJets"),
225  select = cms.string("pt>30 & abs(eta)<2.5 "),
226  min = cms.int32(1),
227  ),
228  cms.PSet(
229  label = cms.string("jets:step2"),
230  src = cms.InputTag("slimmedJets"),
231  select = cms.string("pt>30 & abs(eta)<2.5 "),
232  min = cms.int32(2),
233  ),
234  cms.PSet(
235  label = cms.string("jets:step3"),
236  src = cms.InputTag("slimmedJets"),
237  select = cms.string("pt>30 & abs(eta)<2.5 "),
238  min = cms.int32(3),
239  ),
240  cms.PSet(
241  label = cms.string("jets:step4"),
242  src = cms.InputTag("slimmedJets"),
243  select = cms.string("pt>30 & abs(eta)<2.5 "),
244  min = cms.int32(4),
245  ),
246  )
247 )
248 
249 topSingleElectronLooseDQM_miniAOD = cms.EDAnalyzer("TopSingleLeptonDQM_miniAOD",
250  setup = cms.PSet(
251  directory = cms.string("Physics/Top/TopSingleElectronLooseDQM_miniAOD/"),
252  sources = cms.PSet(
253  muons = cms.InputTag("slimmedMuons"),
254  elecs = cms.InputTag("slimmedElectrons"),
255  jets = cms.InputTag("slimmedJets"),
256  mets = cms.VInputTag("slimmedMETs", "slimmedMETsNoHF", "slimmedMETsPuppi"),
257  pvs = cms.InputTag("offlineSlimmedPrimaryVertices")
258 
259  ),
260  monitoring = cms.PSet(
261  verbosity = cms.string("DEBUG")
262  ),
263  pvExtras = cms.PSet(
264  ),
265  elecExtras = cms.PSet(
266  select = cms.string("pt>20 & abs(eta)<2.5"),
267  isolation = cms.string(ElelooseIsoCut),
268  ),
269  jetExtras = cms.PSet(
270  select = cms.string("pt>30 & abs(eta)<2.5 "),
271  jetBTaggers = cms.PSet(
272  trackCountingEff = cms.PSet(
273  label = cms.InputTag("trackCountingHighEffBJetTags" ),
274  workingPoint = cms.double(1.25)
275  ),
276  trackCountingPur = cms.PSet(
277  label = cms.InputTag("trackCountingHighPurBJetTags" ),
278  workingPoint = cms.double(3.00)
279  ),
280  secondaryVertex = cms.PSet(
281  label = cms.InputTag("simpleSecondaryVertexHighEffBJetTags"),
282  workingPoint = cms.double(2.05)
283  ),
284  cvsVertex = cms.PSet(
285  label = cms.InputTag("combinedSecondaryVertexBJetTags"),
286  workingPoint = cms.double(0.898)
287  )
288  ),
289  ),
290  massExtras = cms.PSet(
291  lowerEdge = cms.double( 70.),
292  upperEdge = cms.double(110.)
293  ),
294  triggerExtras = cms.PSet(
295  src = cms.InputTag("TriggerResults","","HLT"),
296  paths = cms.vstring(['HLT_Ele15_LW_L1R:HLT_QuadJetU15'])
297  )
298  ),
299  preselection = cms.PSet(
300  vertex = cms.PSet(
301  src = cms.InputTag("offlineSlimmedPrimaryVertices")#,
302  )
303  ),
304  selection = cms.VPSet(
305  cms.PSet(
306  label = cms.string("elecs:step0"),
307  src = cms.InputTag("slimmedElectrons"),
308  select = cms.string("pt>20 & abs(eta)<2.5 && "+ElelooseIsoCut),
309  min = cms.int32(1),
310  ),
311  cms.PSet(
312  label = cms.string("jets:step1"),
313  src = cms.InputTag("slimmedJets"),
314  select = cms.string("pt>30 & abs(eta)<2.5 "),
315  min = cms.int32(1),
316  ),
317  cms.PSet(
318  label = cms.string("jets:step2"),
319  src = cms.InputTag("slimmedJets"),
320  select = cms.string("pt>30 & abs(eta)<2.5 "),
321  min = cms.int32(2),
322  ),
323  cms.PSet(
324  label = cms.string("jets:step3"),
325  src = cms.InputTag("slimmedJets"),
326  select = cms.string("pt>30 & abs(eta)<2.5 "),
327  min = cms.int32(3),
328  ),
329  cms.PSet(
330  label = cms.string("jets:step4"),
331  src = cms.InputTag("slimmedJets"),
332  select = cms.string("pt>30 & abs(eta)<2.5 "),
333  min = cms.int32(4),
334  ),
335  )
336 )
337 
338 topSingleElectronMediumDQM_miniAOD = cms.EDAnalyzer("TopSingleLeptonDQM_miniAOD",
339  setup = cms.PSet(
340 
341  directory = cms.string("Physics/Top/TopSingleElectronMediumDQM_miniAOD/"),
342  sources = cms.PSet(
343  muons = cms.InputTag("slimmedMuons"),
344  elecs = cms.InputTag("slimmedElectrons"),
345  jets = cms.InputTag("slimmedJets"),
346  mets = cms.VInputTag("slimmedMETs", "slimmedMETsNoHF", "slimmedMETsPuppi"),
347  pvs = cms.InputTag("offlineSlimmedPrimaryVertices")
348 
349  ),
350  monitoring = cms.PSet(
351  verbosity = cms.string("DEBUG")
352  ),
353  pvExtras = cms.PSet(
354  ),
355  elecExtras = cms.PSet(
356  select = cms.string("pt>20 & abs(eta)<2.5"),
357  isolation = cms.string(ElelooseIsoCut),
358  ),
359  jetExtras = cms.PSet(
360  select = cms.string("pt>30 & abs(eta)<2.5 "),
361  jetBTaggers = cms.PSet(
362  trackCountingEff = cms.PSet(
363  label = cms.InputTag("trackCountingHighEffBJetTags" ),
364  workingPoint = cms.double(1.25)
365  ),
366  trackCountingPur = cms.PSet(
367  label = cms.InputTag("trackCountingHighPurBJetTags" ),
368  workingPoint = cms.double(3.00)
369  ),
370  secondaryVertex = cms.PSet(
371  label = cms.InputTag("simpleSecondaryVertexHighEffBJetTags"),
372  workingPoint = cms.double(2.05)
373  ),
374  cvsVertex = cms.PSet(
375  label = cms.InputTag("combinedSecondaryVertexBJetTags"),
376  workingPoint = cms.double(0.898)
377  )
378  ),
379  ),
380  massExtras = cms.PSet(
381  lowerEdge = cms.double( 70.),
382  upperEdge = cms.double(110.)
383  ),
384  ),
385  preselection = cms.PSet(
386  vertex = cms.PSet(
387  src = cms.InputTag("offlineSlimmedPrimaryVertices")#,
388  )
389  ),
390  selection = cms.VPSet(
391  cms.PSet(
392  label = cms.string("elecs:step0"),
393  src = cms.InputTag("slimmedElectrons"),
394  select = cms.string("pt>30 & abs(eta)<2.5 & abs(gsfTrack.d0)<0.02 & gsfTrack.hitPattern().numberOfHits('MISSING_INNER_HITS') <= 0 & (abs(superCluster.eta) <= 1.4442 || abs(superCluster.eta) >= 1.5660) & " + EletightIsoCut),
395  min = cms.int32(1),
396  max = cms.int32(1),
397  ),
398  cms.PSet(
399  label = cms.string("jets:step1"),
400  src = cms.InputTag("slimmedJets"),
401  select = cms.string("pt>30 & abs(eta)<2.5 "),
402  min = cms.int32(1),
403  ),
404  cms.PSet(
405  label = cms.string("jets:step2"),
406  src = cms.InputTag("slimmedJets"),
407  select = cms.string("pt>30 & abs(eta)<2.5 "),
408  min = cms.int32(2),
409  ),
410  cms.PSet(
411  label = cms.string("jets:step3"),
412  src = cms.InputTag("slimmedJets"),
413  select = cms.string("pt>30 & abs(eta)<2.5 "),
414  min = cms.int32(3),
415  ),
416  cms.PSet(
417  label = cms.string("jets:step4"),
418  src = cms.InputTag("slimmedJets"),
419  select = cms.string("pt>30 & abs(eta)<2.5 "),
420  min = cms.int32(4),
421  ),
422  )
423 )