CMS 3D CMS Logo
singleTopDQM_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 
16 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))"
17 
18 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)"
19 
20 
21 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))"
22 
23 from DQMServices.Core.DQMEDAnalyzer import DQMEDAnalyzer
24 singleTopTChannelLeptonDQM_miniAOD = DQMEDAnalyzer('SingleTopTChannelLeptonDQM_miniAOD',
25 
26  setup = cms.PSet(
27 
28  directory = cms.string("Physics/Top/SingleTopDQM_miniAOD/"),
29 
30  sources = cms.PSet(
31  muons = cms.InputTag("slimmedMuons"),
32  elecs = cms.InputTag("slimmedElectrons"),
33  jets = cms.InputTag("slimmedJets"),
34  mets = cms.VInputTag("slimmedMETs", "slimmedMETsNoHF", "slimmedMETsPuppi"),
35  pvs = cms.InputTag("offlineSlimmedPrimaryVertices")
36  ),
37 
38  monitoring = cms.PSet(
39  verbosity = cms.string("DEBUG")
40  ),
41 
42  pvExtras = cms.PSet(
43  select = cms.string("abs(z) < 24. & position.rho < 2. & ndof > 4 & !isFake")
44  ),
45  elecExtras = cms.PSet(
46 
47  select = cms.string(looseElecCut+ "&& pt>20 & abs(eta)<2.5 & (abs(superCluster.eta) <= 1.4442 || abs(superCluster.eta) >= 1.5660)"),
48  rho = cms.InputTag("fixedGridRhoFastjetAll"),
49  #isolation = cms.string(ElelooseIsoCut),
50 
51  ),
52 
53  muonExtras = cms.PSet(
54 
55  select = cms.string(looseMuonCut + " && pt>20 & abs(eta)<2.1"),
56  isolation = cms.string(looseIsoCut),
57  ),
58 
59  jetExtras = cms.PSet(
60 
61  select = cms.string("pt>30 & abs(eta)<2.4"),
62  ),
63 
64  triggerExtras = cms.PSet(
65  src = cms.InputTag("TriggerResults","","HLT"),
66  paths = cms.vstring(['HLT_Mu3:HLT_QuadJet15U',
67  'HLT_Mu5:HLT_QuadJet15U',
68  'HLT_Mu7:HLT_QuadJet15U',
69  'HLT_Mu9:HLT_QuadJet15U'])
70  )
71  ),
72 
73  preselection = cms.PSet(
74 
75  ),
76 
77  selection = cms.VPSet(
78  cms.PSet(
79  label = cms.string("jets/calo:step0"),
80  src = cms.InputTag("ak4CaloJets"),
81  select = cms.string("pt>20 & abs(eta)<2.1 & 0.05<emEnergyFraction"),
82  jetID = cms.PSet(
83  label = cms.InputTag("ak4JetID"),
84  select = cms.string("fHPD < 0.98 & n90Hits>1 & restrictedEMF<1")
85  ),
86  min = cms.int32(2),
87  ),
88  )
89 )
90 
91 singleTopMuonMediumDQM_miniAOD = DQMEDAnalyzer('SingleTopTChannelLeptonDQM_miniAOD',
92  setup = cms.PSet(
93 
94  directory = cms.string("Physics/Top/SingleTopMuonMediumDQM_miniAOD/"),
95  sources = cms.PSet(
96  muons = cms.InputTag("slimmedMuons"),
97  elecs = cms.InputTag("slimmedElectrons"),
98  jets = cms.InputTag("slimmedJets"),
99  mets = cms.VInputTag("slimmedMETs", "slimmedMETsNoHF", "slimmedMETsPuppi"),
100  pvs = cms.InputTag("offlineSlimmedPrimaryVertices")
101  ),
102 
103  monitoring = cms.PSet(
104  verbosity = cms.string("DEBUG")
105  ),
106 
107  pvExtras = cms.PSet(
108  select = cms.string("abs(z) < 24. & position.rho < 2. & ndof > 4 & !isFake")
109  ),
110  elecExtras = cms.PSet(
111  select = cms.string(tightElecCut + "&& pt>20 & abs(eta)<2.5 & (abs(superCluster.eta) <= 1.4442 || abs(superCluster.eta) >= 1.5660)"),
112  rho = cms.InputTag("fixedGridRhoFastjetAll"),
113  #isolation = cms.string(ElelooseIsoCut),
114  ),
115 
116  muonExtras = cms.PSet(
117  select = cms.string(looseMuonCut + " && pt>20 & abs(eta)<2.1"),
118  isolation = cms.string(looseIsoCut)
119  ),
120  jetExtras = cms.PSet(
121 
122  select = cms.string("pt>30 & abs(eta)<2.4"), # & neutralEmEnergyFraction >0.01 & chargedEmEnergyFraction>0.01"),
123  jetBTaggers = cms.PSet(
124  trackCountingEff = cms.PSet(
125  label = cms.InputTag("trackCountingHighEffBJetTags" ),
126  workingPoint = cms.double(1.25)
127  ),
128  trackCountingPur = cms.PSet(
129  label = cms.InputTag("trackCountingHighPurBJetTags" ),
130  workingPoint = cms.double(3.41)
131  ),
132  secondaryVertex = cms.PSet(
133  label = cms.InputTag("simpleSecondaryVertexHighEffBJetTags"),
134  workingPoint = cms.double(2.05)
135  ),
136  combinedSecondaryVertex = cms.PSet(
137  label = cms.InputTag("combinedSecondaryVertexBJetTags"),
138  workingPoint = cms.double(0.898)
139  )
140  ),
141  ),
142  massExtras = cms.PSet(
143  lowerEdge = cms.double( 70.),
144  upperEdge = cms.double(110.)
145  ),
146  ),
147  preselection = cms.PSet(
148  vertex = cms.PSet(
149  src = cms.InputTag("offlineSlimmedPrimaryVertices"),#,
150  select = cms.string("abs(z) < 24. & position.rho < 2. & ndof > 4 & !isFake")
151  )
152  ),
153 
154  selection = cms.VPSet(
155  cms.PSet(
156  label = cms.string("muons:step0"),
157  src = cms.InputTag("slimmedMuons"),
158  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
159  min = cms.int32(1),
160  #max = cms.int32(1),
161  ),
162  cms.PSet(
163  label = cms.string("jets:step1"),
164  src = cms.InputTag("slimmedJets"),
165  select = cms.string("pt>30 & abs(eta)<2.4 "),
166  min = cms.int32(2),
167  ),
168  cms.PSet(
169  label = cms.string("met:step2"),
170  src = cms.InputTag("slimmedMETs"),
171  select = cms.string("pt>30"),
172  #min = cms.int32(2),
173  ),
174  )
175 )
176 
177 singleTopElectronMediumDQM_miniAOD = DQMEDAnalyzer('SingleTopTChannelLeptonDQM_miniAOD',
178 
179  setup = cms.PSet(
180 
181  directory = cms.string("Physics/Top/SingleTopElectronMediumDQM_miniAOD/"),
182 
183  sources = cms.PSet(
184  muons = cms.InputTag("slimmedMuons"),
185  elecs = cms.InputTag("slimmedElectrons"),
186  jets = cms.InputTag("slimmedJets"),
187  mets = cms.VInputTag("slimmedMETs", "slimmedMETsNoHF", "slimmedMETsPuppi"),
188  pvs = cms.InputTag("offlineSlimmedPrimaryVertices")
189 
190  ),
191 
192  monitoring = cms.PSet(
193  verbosity = cms.string("DEBUG")
194  ),
195 
196  pvExtras = cms.PSet(
197  select = cms.string("abs(z) < 24. & position.rho < 2. & ndof > 4 & !isFake")
198  ),
199  elecExtras = cms.PSet(
200  select = cms.string(tightElecCut + " && pt>20 & abs(eta)<2.5 & (abs(superCluster.eta) <= 1.4442 || abs(superCluster.eta) >= 1.5660)"),
201  #select = cms.string(looseElecCut+ "&& pt>20 & abs(eta)<2.5 & (abs(superCluster.eta) <= 1.4442 || abs(superCluster.eta) >= 1.5660)"),
202  rho = cms.InputTag("fixedGridRhoFastjetAll"),
203  #isolation = cms.string(ElelooseIsoCut),
204  ),
205  muonExtras = cms.PSet(
206  select = cms.string(looseMuonCut + " && pt>20 & abs(eta)<2.1"),
207  isolation = cms.string(looseIsoCut)
208  ),
209 
210  jetExtras = cms.PSet(
211 
212  select = cms.string("pt>30 & abs(eta)<2.4"),
213  jetBTaggers = cms.PSet(
214  trackCountingEff = cms.PSet(
215  label = cms.InputTag("trackCountingHighEffBJetTags" ),
216  workingPoint = cms.double(1.25)
217  ),
218  trackCountingPur = cms.PSet(
219  label = cms.InputTag("trackCountingHighPurBJetTags" ),
220  workingPoint = cms.double(3.41)
221  ),
222  secondaryVertex = cms.PSet(
223  label = cms.InputTag("simpleSecondaryVertexHighEffBJetTags"),
224  workingPoint = cms.double(2.05)
225  ),
226  combinedSecondaryVertex = cms.PSet(
227  label = cms.InputTag("combinedSecondaryVertexBJetTags"),
228  workingPoint = cms.double(0.898)
229  ),
230  ),
231  ),
232  massExtras = cms.PSet(
233  lowerEdge = cms.double( 70.),
234  upperEdge = cms.double(110.)
235  ),
236 
237  ),
238 
239  preselection = cms.PSet(
240  vertex = cms.PSet(
241  src = cms.InputTag("offlineSlimmedPrimaryVertices"),#,
242  select = cms.string("abs(z) < 24. & position.rho < 2. & ndof > 4 & !isFake")
243  )
244  ),
245 
246  selection = cms.VPSet(
247  cms.PSet(
248  label = cms.string("elecs:step0"),
249  src = cms.InputTag("slimmedElectrons"),
250  select = cms.string("pt>20 & abs(eta)<2.5 & (abs(superCluster.eta) <= 1.4442 || abs(superCluster.eta) >= 1.5660) &&" + tightElecCut),
251  min = cms.int32(1),
252  max = cms.int32(1),
253  ),
254  cms.PSet(
255  label = cms.string("jets:step1"),
256  src = cms.InputTag("slimmedJets"),
257  select = cms.string("pt>30 & abs(eta)<2.4 "),
258  min = cms.int32(2),
259  ),
260  cms.PSet(
261  label = cms.string("met:step2"),
262  src = cms.InputTag("slimmedMETs"),
263  select = cms.string("pt>30"),
264  #min = cms.int32(2),
265  ),
266  )
267 )
Generated for CMSSW Reference Manual by doxygen 1.8.11