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 
24 
25 topSingleMuonMediumDQM_miniAOD = DQMEDAnalyzer('TopSingleLeptonDQM_miniAOD',
26  setup = cms.PSet(
27  directory = cms.string("Physics/Top/TopSingleMuonMediumDQM_miniAOD/"),
28  sources = cms.PSet(
29  muons = cms.InputTag("slimmedMuons"),
30  elecs = cms.InputTag("slimmedElectrons"),
31  jets = cms.InputTag("slimmedJets"),
32  mets = cms.VInputTag("slimmedMETs", "slimmedMETsNoHF", "slimmedMETsPuppi"),
33  pvs = cms.InputTag("offlineSlimmedPrimaryVertices")
34 
35  ),
36  monitoring = cms.PSet(
37  verbosity = cms.string("DEBUG")
38  ),
39  pvExtras = cms.PSet(
40  select = cms.string("abs(z) < 24. & position.rho < 2. & ndof > 4 & !isFake")
41  ),
42  elecExtras = cms.PSet(
43  select = cms.string(tightElecCut + "&& pt>20 & abs(eta)<2.5 & (abs(superCluster.eta) <= 1.4442 || abs(superCluster.eta) >= 1.5660)"),
44  rho = cms.InputTag("fixedGridRhoFastjetAll"),
45  #isolation = cms.string(ElelooseIsoCut),
46  ),
47 
48  muonExtras = cms.PSet(
49  select = cms.string(tightMuonCut + " && pt>20 & abs(eta)<2.4 && " + looseIsoCut),
50  isolation = cms.string(looseIsoCut)
51  ),
52  jetExtras = cms.PSet(
53  select = cms.string("pt>30 & abs(eta)<2.4 "),
54  jetBTaggers = cms.PSet(
55  trackCountingEff = cms.PSet(
56  label = cms.InputTag("trackCountingHighEffBJetTags" ),
57  workingPoint = cms.double(1.25)
58  ),
59  trackCountingPur = cms.PSet(
60  label = cms.InputTag("trackCountingHighPurBJetTags" ),
61  workingPoint = cms.double(3.00)
62  ),
63  secondaryVertex = cms.PSet(
64  label = cms.InputTag("simpleSecondaryVertexHighEffBJetTags"),
65  workingPoint = cms.double(2.05)
66  ),
67  cvsVertex = cms.PSet(
68  label = cms.InputTag("combinedSecondaryVertexBJetTags"),
69  workingPoint = cms.double(0.898)
70  )
71  ),
72  ),
73  massExtras = cms.PSet(
74  lowerEdge = cms.double( 70.),
75  upperEdge = cms.double(110.)
76  ),
77  ),
78  preselection = cms.PSet(
79  vertex = cms.PSet(
80  src = cms.InputTag("offlineSlimmedPrimaryVertices"),#,
81  select = cms.string("abs(z) < 24. & position.rho < 2. & ndof > 4 & !isFake")
82  )
83  ),
84  selection = cms.VPSet(
85  cms.PSet(
86  label = cms.string("muons:step0"),
87  src = cms.InputTag("slimmedMuons"),
88  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
89  min = cms.int32(1),
90  #max = cms.int32(1),
91  ),
92  cms.PSet(
93  label = cms.string("jets:step1"),
94  src = cms.InputTag("slimmedJets"),
95  select = cms.string("pt>30 & abs(eta)<2.4 "),
96  min = cms.int32(4),
97  ),
98  cms.PSet(
99  label = cms.string("met:step2"),
100  src = cms.InputTag("slimmedMETs"),
101  select = cms.string("pt>30"),
102  #min = cms.int32(2),
103  ),
104 # cms.PSet(
105 # label = cms.string("jets:step3"),
106 # src = cms.InputTag("slimmedJets"),
107 # select = cms.string("pt>30 & abs(eta)<2.5 "),
108 # min = cms.int32(3),
109 # ),
110 # cms.PSet(
111 # label = cms.string("jets:step4"),
112 # src = cms.InputTag("slimmedJets"),
113 # select = cms.string("pt>30 & abs(eta)<2.5 "),
114 # min = cms.int32(4),
115 # ),
116  )
117 )
118 
119 
120 
121 topSingleElectronMediumDQM_miniAOD = DQMEDAnalyzer('TopSingleLeptonDQM_miniAOD',
122  setup = cms.PSet(
123 
124  directory = cms.string("Physics/Top/TopSingleElectronMediumDQM_miniAOD/"),
125  sources = cms.PSet(
126  muons = cms.InputTag("slimmedMuons"),
127  elecs = cms.InputTag("slimmedElectrons"),
128  jets = cms.InputTag("slimmedJets"),
129  mets = cms.VInputTag("slimmedMETs", "slimmedMETsNoHF", "slimmedMETsPuppi"),
130  pvs = cms.InputTag("offlineSlimmedPrimaryVertices")
131 
132  ),
133  monitoring = cms.PSet(
134  verbosity = cms.string("DEBUG")
135  ),
136  pvExtras = cms.PSet(
137  select = cms.string("abs(z) < 24. & position.rho < 2. & ndof > 4 & !isFake")
138  ),
139  elecExtras = cms.PSet(
140  select = cms.string(tightElecCut + " && pt>20 & abs(eta)<2.5 & (abs(superCluster.eta) <= 1.4442 || abs(superCluster.eta) >= 1.5660)"),
141  #select = cms.string(looseElecCut+ "&& pt>20 & abs(eta)<2.5 & (abs(superCluster.eta) <= 1.4442 || abs(superCluster.eta) >= 1.5660)"),
142  rho = cms.InputTag("fixedGridRhoFastjetAll"),
143  #isolation = cms.string(ElelooseIsoCut),
144  ),
145  muonExtras = cms.PSet(
146  select = cms.string(tightMuonCut + " && pt>20 & abs(eta)<2.4 && " + looseIsoCut),
147  isolation = cms.string(looseIsoCut)
148  ),
149  jetExtras = cms.PSet(
150  select = cms.string("pt>30 & abs(eta)<2.4 "),
151  jetBTaggers = cms.PSet(
152  trackCountingEff = cms.PSet(
153  label = cms.InputTag("trackCountingHighEffBJetTags" ),
154  workingPoint = cms.double(1.25)
155  ),
156  trackCountingPur = cms.PSet(
157  label = cms.InputTag("trackCountingHighPurBJetTags" ),
158  workingPoint = cms.double(3.00)
159  ),
160  secondaryVertex = cms.PSet(
161  label = cms.InputTag("simpleSecondaryVertexHighEffBJetTags"),
162  workingPoint = cms.double(2.05)
163  ),
164  cvsVertex = cms.PSet(
165  label = cms.InputTag("combinedSecondaryVertexBJetTags"),
166  workingPoint = cms.double(0.898)
167  )
168  ),
169  ),
170  massExtras = cms.PSet(
171  lowerEdge = cms.double( 70.),
172  upperEdge = cms.double(110.)
173  ),
174  ),
175  preselection = cms.PSet(
176  vertex = cms.PSet(
177  src = cms.InputTag("offlineSlimmedPrimaryVertices"),#,
178  select = cms.string("abs(z) < 24. & position.rho < 2. & ndof > 4 & !isFake")
179  )
180  ),
181  selection = cms.VPSet(
182  cms.PSet(
183  label = cms.string("elecs:step0"),
184  src = cms.InputTag("slimmedElectrons"),
185  select = cms.string("pt>20 & abs(eta)<2.5 & (abs(superCluster.eta) <= 1.4442 || abs(superCluster.eta) >= 1.5660) &" + tightElecCut),
186  min = cms.int32(1),
187  max = cms.int32(1),
188  ),
189  cms.PSet(
190  label = cms.string("jets:step1"),
191  src = cms.InputTag("slimmedJets"),
192  select = cms.string("pt>30 & abs(eta)<2.4 "),
193  min = cms.int32(4),
194  ),
195  cms.PSet(
196  label = cms.string("met:step2"),
197  src = cms.InputTag("slimmedMETs"),
198  select = cms.string("pt>30"),
199  #min = cms.int32(2),
200  ),
201  #cms.PSet(
202  #label = cms.string("jets:step3"),
203  #src = cms.InputTag("slimmedJets"),
204  #select = cms.string("pt>30 & abs(eta)<2.5 "),
205  #min = cms.int32(3),
206  #),
207  #cms.PSet(
208  #label = cms.string("jets:step4"),
209  #src = cms.InputTag("slimmedJets"),
210  #select = cms.string("pt>30 & abs(eta)<2.5 "),
211  #min = cms.int32(4),
212  #),
213  )
214 )
DQMEDAnalyzer
Definition: DQMEDAnalyzer.py:1
Generated for CMSSW Reference Manual by doxygen 1.8.16