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