CMS 3D CMS Logo

topDiLeptonOfflineDQM_cfi.py
Go to the documentation of this file.
1 import FWCore.ParameterSet.Config as cms
2 
3 looseMuonCut = "muonRef.isNonnull && (muonRef.isGlobalMuon || muonRef.isTrackerMuon) && muonRef.isPFMuon"
4 looseIsoCut = "(muonRef.pfIsolationR04.sumChargedHadronPt + max(0., muonRef.pfIsolationR04.sumNeutralHadronEt + muonRef.pfIsolationR04.sumPhotonEt - 0.5 * muonRef.pfIsolationR04.sumPUPt) ) / muonRef.pt < 0.2"
5 ElelooseIsoCut = "(gsfElectronRef.pfIsolationVariables.sumChargedHadronPt + max(0., gsfElectronRef.pfIsolationVariables.sumNeutralHadronEt + gsfElectronRef.pfIsolationVariables.sumPhotonEt - 0.5 * gsfElectronRef.pfIsolationVariables.sumPUPt) ) / gsfElectronRef.pt < 0.15"
6 EletightIsoCut = "(gsfElectronRef.pfIsolationVariables.sumChargedHadronPt + max(0., gsfElectronRef.pfIsolationVariables.sumNeutralHadronEt + gsfElectronRef.pfIsolationVariables.sumPhotonEt - 0.5 * gsfElectronRef.pfIsolationVariables.sumPUPt) ) / gsfElectronRef.pt < 0.1"
7 
8 
9 from DQMServices.Core.DQMEDAnalyzer import DQMEDAnalyzer
10 topDiLeptonOfflineDQM = DQMEDAnalyzer('TopDiLeptonOfflineDQM',
11 
17  setup = cms.PSet(
18 
21  directory = cms.string("Physics/Top/TopDiLeptonDQM/"),
22 
23 
24  sources = cms.PSet(
25  muons = cms.InputTag("pfIsolatedMuonsEI"),
26  elecs = cms.InputTag("pfIsolatedElectronsEI"),
27  jets = cms.InputTag("ak4PFJetsCHS"),
28  mets = cms.VInputTag("caloMet", "tcMet", "pfMet")
29  ),
30 
31  monitoring = cms.PSet(
32  verbosity = cms.string("DEBUG")
33  ),
34 
36  elecExtras = cms.PSet(
37 
41  select = cms.string("pt>20. && abs(eta)<2.5"),
42 
44  isolation = cms.string(ElelooseIsoCut),
45  ),
46 
48  muonExtras = cms.PSet(
49 
51  select = cms.string(looseMuonCut + " && muonRef.pt > 10. && abs(muonRef.eta)<2.4"),
52 
54  isolation = cms.string(looseIsoCut),
55  ),
56 
58  jetExtras = cms.PSet(
59 
61  jetCorrector = cms.string("topDQMak5PFCHSL2L3"),
62 
71  select = cms.string("pt>30. & abs(eta)<2.4 "),
72  ),
73 
75  massExtras = cms.PSet(
76  lowerEdge = cms.double( 76.0),
77  upperEdge = cms.double(106.0)
78  ),
79 
93  ),
94 
95 
102  preselection = cms.PSet(
103 
109  vertex = cms.PSet(
110  src = cms.InputTag("offlinePrimaryVertices"),
111  select = cms.string('abs(x)<1. && abs(y)<1. && abs(z)<20. && tracksSize>3 && !isFake')
112  )
113  ),
114 
115 
124  selection = cms.VPSet(
125  #cms.PSet(
126 
131  cms.PSet(
132  label = cms.string("muons:step0"),
133  src = cms.InputTag("pfIsolatedMuonsEI"),
134  select = cms.string(looseMuonCut +" && "+ looseIsoCut + " && muonRef.pt > 20. && abs(muonRef.eta)<2.4"), # CB what to do with iso? CD Added looseIso
135  min = cms.int32(2),
136  max = cms.int32(2),
137  ),
138  cms.PSet(
139  label = cms.string("jets/pf:step1"),
140  src = cms.InputTag("ak4PFJetsCHS"),
141  jetCorrector = cms.string("topDQMak5PFCHSL2L3"),
142 # select = cms.string("pt>30. & abs(eta)<2.4 & emEnergyFraction>0.01"),
143 # jetID = cms.PSet(
144 # label = cms.InputTag("ak5JetID"),
145 # select = cms.string("fHPD < 0.98 & n90Hits>1 & restrictedEMF<1")
146 # ),
147  min = cms.int32(2),
148  #max = cms.int32(2),
149  )
150  )
151 )
152 
153 
154 
155 DiMuonDQM = DQMEDAnalyzer('TopDiLeptonOfflineDQM',
156 
162  setup = cms.PSet(
163 
166  directory = cms.string("Physics/Top/TopDiMuonDQM/"),
167 
168 
169  sources = cms.PSet(
170  muons = cms.InputTag("pfIsolatedMuonsEI"),
171  elecs = cms.InputTag("pfIsolatedElectronsEI"),
172  jets = cms.InputTag("ak4PFJetsCHS"),
173  mets = cms.VInputTag("caloMet", "tcMet", "pfMet")
174  ),
175 
176  monitoring = cms.PSet(
177  verbosity = cms.string("DEBUG")
178  ),
179 
181  elecExtras = cms.PSet(
182 
186  select = cms.string("pt>20. && abs(eta)<2.5"),
187 
189  isolation = cms.string(ElelooseIsoCut),
190  ),
191 
193  muonExtras = cms.PSet(
194 
196  select = cms.string(looseMuonCut + " && muonRef.pt > 20. && abs(muonRef.eta)<2.4"),
197 
199  isolation = cms.string(looseIsoCut),
200  ),
201 
203  jetExtras = cms.PSet(
204 
206  jetCorrector = cms.string("topDQMak5PFCHSL2L3"),
207 
216  select = cms.string("pt>30. & abs(eta)<2.4 "),
217  ),
218 
220  massExtras = cms.PSet(
221  lowerEdge = cms.double( 76.0),
222  upperEdge = cms.double(106.0)
223  ),
224 
238  ),
239 
240 
247  preselection = cms.PSet(
248 
254  vertex = cms.PSet(
255  src = cms.InputTag("offlinePrimaryVertices"),
256  select = cms.string('abs(x)<1. && abs(y)<1. && abs(z)<20. && tracksSize>3 && !isFake')
257  )
258  ),
259 
260 
269  selection = cms.VPSet(
270  #cms.PSet(
271 
276  cms.PSet(
277  label = cms.string("muons:step0"),
278  src = cms.InputTag("pfIsolatedMuonsEI"),
279  select = cms.string(looseMuonCut + " && muonRef.pt > 20. && abs(muonRef.eta)<2.4"), # CB what to do with iso?
280  min = cms.int32(2),
281  max = cms.int32(2),
282  ),
283  cms.PSet(
284  label = cms.string("jets/pf:step1"),
285  src = cms.InputTag("ak4PFJetsCHS"),
286  jetCorrector = cms.string("topDQMak5PFCHSL2L3"),
287  select = cms.string("pt>30. & abs(eta)<2.4 "),
288 # jetID = cms.PSet(
289 # label = cms.InputTag("ak5JetID"),
290 # select = cms.string("fHPD < 0.98 & n90Hits>1 & restrictedEMF<1")
291 # ),
292  min = cms.int32(2),
293  #max = cms.int32(2),
294  ),
295  ),
296 )
297 
298 DiElectronDQM = DQMEDAnalyzer('TopDiLeptonOfflineDQM',
299 
305  setup = cms.PSet(
306 
309  directory = cms.string("Physics/Top/TopDiElectronDQM/"),
310 
311 
312  sources = cms.PSet(
313  muons = cms.InputTag("pfIsolatedMuonsEI"),
314  elecs = cms.InputTag("pfIsolatedElectronsEI"),
315  jets = cms.InputTag("ak4PFJetsCHS"),
316  mets = cms.VInputTag("caloMet", "tcMet", "pfMet")
317  ),
318 
319  monitoring = cms.PSet(
320  verbosity = cms.string("DEBUG")
321  ),
322 
324  elecExtras = cms.PSet(
325 
329  select = cms.string("pt>20. && abs(eta)<2.5"),
330 
332  isolation = cms.string(ElelooseIsoCut),
333  ),
334 
336  muonExtras = cms.PSet(
337 
339  select = cms.string(looseMuonCut + " && muonRef.pt > 20. && abs(muonRef.eta)<2.4"),
340 
342  isolation = cms.string(looseIsoCut),
343  ),
344 
346  jetExtras = cms.PSet(
347 
349  jetCorrector = cms.string("topDQMak5PFCHSL2L3"),
350 
359  select = cms.string("pt>30. & abs(eta)<2.4 "),
360  ),
361 
363  massExtras = cms.PSet(
364  lowerEdge = cms.double( 76.0),
365  upperEdge = cms.double(106.0)
366  ),
367 
381  ),
382 
383 
390  preselection = cms.PSet(
391 
397  vertex = cms.PSet(
398  src = cms.InputTag("offlinePrimaryVertices"),
399  select = cms.string('abs(x)<1. && abs(y)<1. && abs(z)<20. && tracksSize>3 && !isFake')
400  )
401  ),
402 
403 
412  selection = cms.VPSet(
413  #cms.PSet(
414 
419  cms.PSet(
420  label = cms.string("elecs:step0"),
421  src = cms.InputTag("pfIsolatedElectronsEI"),
422 
423  select = cms.string("pt>20 & abs(eta)<2.5 && gsfElectronRef.gsfTrack.hitPattern().numberOfLostHits('MISSING_INNER_HITS') <= 0 && " + ElelooseIsoCut),
424  #abs(gsfElectronRef.gsfTrack.d0)<0.04
425  min = cms.int32(2),
426  max = cms.int32(2),
427  ),
428  cms.PSet(
429  label = cms.string("jets/pf:step1"),
430  src = cms.InputTag("ak4PFJetsCHS"),
431  jetCorrector = cms.string("topDQMak5PFCHSL2L3"),
432  select = cms.string("pt>30. & abs(eta)<2.4"),
433 # jetID = cms.PSet(
434 # label = cms.InputTag("ak5JetID"),
435 # select = cms.string("fHPD < 0.98 & n90Hits>1 & restrictedEMF<1")
436 # ),
437  min = cms.int32(2),
438  #max = cms.int32(2),
439  ),
440  ),
441 )
442 
443 ElecMuonDQM = DQMEDAnalyzer('TopDiLeptonOfflineDQM',
444 
450  setup = cms.PSet(
451 
454  directory = cms.string("Physics/Top/TopElecMuonDQM/"),
455 
456 
457  sources = cms.PSet(
458  muons = cms.InputTag("pfIsolatedMuonsEI"),
459  elecs = cms.InputTag("pfIsolatedElectronsEI"),
460  jets = cms.InputTag("ak4PFJetsCHS"),
461  mets = cms.VInputTag("caloMet", "tcMet", "pfMet")
462  ),
463 
464  monitoring = cms.PSet(
465  verbosity = cms.string("DEBUG")
466  ),
467 
469  elecExtras = cms.PSet(
470 
474  select = cms.string("pt>10. && abs(eta)<2.4 && abs(gsfElectronRef.gsfTrack.d0)<1. && abs(gsfElectronRef.gsfTrack.dz)<20."),
475 
477  isolation = cms.string(ElelooseIsoCut),
478  ),
479 
481  muonExtras = cms.PSet(
482 
484  select = cms.string(looseMuonCut + " && muonRef.pt > 10. && abs(muonRef.eta)<2.4"),
485 
487  isolation = cms.string(looseIsoCut),
488  ),
489 
491  jetExtras = cms.PSet(
492 
494  jetCorrector = cms.string("topDQMak5PFCHSL2L3"),
495 
504  select = cms.string("pt>30. & abs(eta)<2.4 "),
505  ),
506 
508  massExtras = cms.PSet(
509  lowerEdge = cms.double( 76.0),
510  upperEdge = cms.double(106.0)
511  ),
512 
526  ),
527 
528 
535  preselection = cms.PSet(
536 
542  vertex = cms.PSet(
543  src = cms.InputTag("offlinePrimaryVertices"),
544  select = cms.string('abs(x)<1. && abs(y)<1. && abs(z)<20. && tracksSize>3 && !isFake')
545  )
546  ),
547 
548 
557  selection = cms.VPSet(
558  #cms.PSet(
559 
564  cms.PSet(
565  label = cms.string("muons:step0"),
566  src = cms.InputTag("pfIsolatedMuonsEI"),
567  select = cms.string(looseMuonCut + " && " + looseIsoCut + " && muonRef.pt > 20. && abs(muonRef.eta)<2.4"), # CB what to do with iso? CD Added looseIsoCut
568  min = cms.int32(1),
569  max = cms.int32(1),
570  ),
571  cms.PSet(
572  label = cms.string("elecs:step1"),
573  src = cms.InputTag("pfIsolatedElectronsEI"),
574 
575  select = cms.string("pt>20 & abs(eta)<2.5 && "+ElelooseIsoCut),
576  min = cms.int32(1),
577  max = cms.int32(1),
578  ),
579  cms.PSet(
580  label = cms.string("jets/pf:step2"),
581  src = cms.InputTag("ak4PFJetsCHS"),
582  jetCorrector = cms.string("topDQMak5PFCHSL2L3"),
583  select = cms.string("pt>30. & abs(eta)<2.4 "),
584 # jetID = cms.PSet(
585 # label = cms.InputTag("ak5JetID"),
586 # select = cms.string("fHPD < 0.98 & n90Hits>1 & restrictedEMF<1")
587 # ),
588  min = cms.int32(2),
589  #max = cms.int32(2),
590  ),
591  ),
592 )
593 
DQMEDAnalyzer
Definition: DQMEDAnalyzer.py:1