test
CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
singleTopDQM_miniAOD_cfi.py
Go to the documentation of this file.
1 import FWCore.ParameterSet.Config as cms
2 
3 EletightIsoCut = "(pfIsolationVariables.sumChargedHadronPt + max(0., pfIsolationVariables.sumNeutralHadronEt + pfIsolationVariables.sumPhotonEt - 0.5 * pfIsolationVariables.sumPUPt) ) / pt < 0.1"
4 ElelooseIsoCut = "(pfIsolationVariables.sumChargedHadronPt + max(0., pfIsolationVariables.sumNeutralHadronEt + pfIsolationVariables.sumPhotonEt - 0.5 * pfIsolationVariables.sumPUPt) ) / pt < 0.15"
5 
6 
7 singleTopTChannelLeptonDQM_miniAOD = cms.EDAnalyzer("SingleTopTChannelLeptonDQM_miniAOD",
8 
9  setup = cms.PSet(
10 
11  directory = cms.string("Physics/Top/SingleTopDQM_miniAOD/"),
12 
13  sources = cms.PSet(
14  muons = cms.InputTag("slimmedMuons"),
15  elecs = cms.InputTag("slimmedElectrons"),
16  jets = cms.InputTag("slimmedJets"),
17  mets = cms.VInputTag("slimmedMETs", "slimmedMETsNoHF", "slimmedMETsPuppi"),
18  pvs = cms.InputTag("offlineSlimmedPrimaryVertices")
19  ),
20 
21  monitoring = cms.PSet(
22  verbosity = cms.string("DEBUG")
23  ),
24 
25  elecExtras = cms.PSet(
26 
27  select = cms.string("pt>15 & abs(eta)<2.5 & abs(gsfTrack.d0)<1 & abs(gsfTrack.dz)<20"),
28 
29  isolation = cms.string(ElelooseIsoCut),
30  ),
31 
32  muonExtras = cms.PSet(
33 
34  select = cms.string("pt>10 & abs(eta)<2.1 & isGlobalMuon & abs(globalTrack.d0)<1 & abs(globalTrack.dz)<20"),
35 
36  ),
37 
38  jetExtras = cms.PSet(
39 
40  select = cms.string("pt>15 & abs(eta)<2.5 & emEnergyFraction>0.01"),
41  ),
42 
43  triggerExtras = cms.PSet(
44  src = cms.InputTag("TriggerResults","","HLT"),
45  paths = cms.vstring(['HLT_Mu3:HLT_QuadJet15U',
46  'HLT_Mu5:HLT_QuadJet15U',
47  'HLT_Mu7:HLT_QuadJet15U',
48  'HLT_Mu9:HLT_QuadJet15U'])
49  )
50  ),
51 
52  preselection = cms.PSet(
53 
54  ),
55 
56  selection = cms.VPSet(
57  cms.PSet(
58  label = cms.string("jets/calo:step0"),
59  src = cms.InputTag("ak4CaloJets"),
60  select = cms.string("pt>20 & abs(eta)<2.1 & 0.05<emEnergyFraction"),
61  jetID = cms.PSet(
62  label = cms.InputTag("ak4JetID"),
63  select = cms.string("fHPD < 0.98 & n90Hits>1 & restrictedEMF<1")
64  ),
65  min = cms.int32(2),
66  )
67  )
68 )
69 
70 singleTopMuonMediumDQM_miniAOD = cms.EDAnalyzer("SingleTopTChannelLeptonDQM_miniAOD",
71 
72  setup = cms.PSet(
73 
74  directory = cms.string("Physics/Top/SingleTopMuonMediumDQM_miniAOD/"),
75 
76  sources = cms.PSet(
77  muons = cms.InputTag("slimmedMuons"),
78  elecs_gsf = cms.InputTag("slimmedElectrons"),
79  elecs = cms.InputTag("slimmedElectrons"),
80  jets = cms.InputTag("slimmedJets"),
81  mets = cms.VInputTag("slimmedMETs", "slimmedMETsNoHF", "slimmedMETsPuppi"),
82  pvs = cms.InputTag("offlineSlimmedPrimaryVertices")
83  ),
84 
85  monitoring = cms.PSet(
86  verbosity = cms.string("DEBUG")
87  ),
88 
89  muonExtras = cms.PSet(
90 
91  select = cms.string("abs(eta)<2.1")
92 
93  ),
94 
95  jetExtras = cms.PSet(
96 
97  select = cms.string("pt>15 & abs(eta)<2.5"), # & neutralEmEnergyFraction >0.01 & chargedEmEnergyFraction>0.01"),
98  jetBTaggers = cms.PSet(
99  trackCountingEff = cms.PSet(
100  label = cms.InputTag("trackCountingHighEffBJetTags" ),
101  workingPoint = cms.double(1.25)
102  ),
103  trackCountingPur = cms.PSet(
104  label = cms.InputTag("trackCountingHighPurBJetTags" ),
105  workingPoint = cms.double(3.41)
106  ),
107  secondaryVertex = cms.PSet(
108  label = cms.InputTag("simpleSecondaryVertexHighEffBJetTags"),
109  workingPoint = cms.double(2.05)
110  ),
111  combinedSecondaryVertex = cms.PSet(
112  label = cms.InputTag("combinedSecondaryVertexBJetTags"),
113  workingPoint = cms.double(0.898)
114  )
115  )
116  )
117 
118  ),
119 
120  preselection = cms.PSet(
121 
122  ),
123 
124  selection = cms.VPSet(
125  cms.PSet(
126  label = cms.string("presel"),
127  src = cms.InputTag("offlineSlimmedPrimaryVertices"),
128 
129 
130  ),
131  cms.PSet(
132  label = cms.string("muons:step0"),
133  src = cms.InputTag("slimmedMuons"),
134  select = cms.string("pt>20 & abs(eta)<2.1 & isGlobalMuon & isTrackerMuon & innerTrack.numberOfValidHits>10 & globalTrack.hitPattern.numberOfValidMuonHits>0 & globalTrack.normalizedChi2<10 & innerTrack.hitPattern.pixelLayersWithMeasurement>=1 & numberOfMatches>1 & abs(innerTrack.dxy)<0.02 & (pfIsolationR04.sumChargedHadronPt + pfIsolationR04.sumNeutralHadronEt + pfIsolationR04.sumPhotonEt)/pt < 0.15"),
135 
136  min = cms.int32(1),
137  max = cms.int32(1),
138  ),
139  cms.PSet(
140  label = cms.string("jets:step1"),
141  src = cms.InputTag("slimmedJets"),
142 
143  select = cms.string(" pt>30 & abs(eta)<4.5 & numberOfDaughters>1 & ((abs(eta)>2.4) || ( chargedHadronEnergyFraction > 0 & chargedMultiplicity>0 & chargedEmEnergyFraction<0.99)) & neutralEmEnergyFraction < 0.99 & neutralHadronEnergyFraction < 0.99"),
144 
145  min = cms.int32(1),
146  max = cms.int32(1),
147  ),
148  cms.PSet(
149  label = cms.string("jets:step2"),
150  src = cms.InputTag("slimmedJets"),
151 
152  select = cms.string(" pt>30 & abs(eta)<4.5 & numberOfDaughters>1 & ((abs(eta)>2.4) || ( chargedHadronEnergyFraction > 0 & chargedMultiplicity>0 & chargedEmEnergyFraction<0.99)) & neutralEmEnergyFraction < 0.99 & neutralHadronEnergyFraction < 0.99"),
153 
154  min = cms.int32(2),
155  max = cms.int32(2),
156  )
157  )
158 )
159 
160 singleTopElectronMediumDQM_miniAOD = cms.EDAnalyzer("SingleTopTChannelLeptonDQM_miniAOD",
161 
162  setup = cms.PSet(
163 
164  directory = cms.string("Physics/Top/SingleTopElectronMediumDQM_miniAOD/"),
165 
166  sources = cms.PSet(
167  muons = cms.InputTag("slimmedMuons"),
168  elecs_gsf = cms.InputTag("slimmedElectrons"),
169  elecs = cms.InputTag("slimmedElectrons"),
170  jets = cms.InputTag("slimmedJets"),
171  mets = cms.VInputTag("slimmedMETs", "slimmedMETsNoHF", "slimmedMETsPuppi"),
172  pvs = cms.InputTag("offlineSlimmedPrimaryVertices")
173 
174  ),
175 
176  monitoring = cms.PSet(
177  verbosity = cms.string("DEBUG")
178  ),
179 
180  elecExtras = cms.PSet(
181 
182  select = cms.string("pt>25"),
183 
184 
185  ),
186 
187  jetExtras = cms.PSet(
188 
189  select = cms.string("pt>15 & abs(eta)<2.5"),
190  jetBTaggers = cms.PSet(
191  trackCountingEff = cms.PSet(
192  label = cms.InputTag("trackCountingHighEffBJetTags" ),
193  workingPoint = cms.double(1.25)
194  ),
195  trackCountingPur = cms.PSet(
196  label = cms.InputTag("trackCountingHighPurBJetTags" ),
197  workingPoint = cms.double(3.41)
198  ),
199  secondaryVertex = cms.PSet(
200  label = cms.InputTag("simpleSecondaryVertexHighEffBJetTags"),
201  workingPoint = cms.double(2.05)
202  ),
203  combinedSecondaryVertex = cms.PSet(
204  label = cms.InputTag("combinedSecondaryVertexBJetTags"),
205  workingPoint = cms.double(0.898)
206  )
207  )
208  ),
209  ),
210 
211  preselection = cms.PSet(
212 
213  ),
214 
215  selection = cms.VPSet(
216  cms.PSet(
217  label = cms.string("presel"),
218  src = cms.InputTag("offlineSlimmedPrimaryVertices"),
219  ),
220  cms.PSet(
221  label = cms.string("elecs:step0"),
222  src = cms.InputTag("slimmedElectrons"),
223  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),
224  min = cms.int32(1),
225  max = cms.int32(1),
226  ),
227  cms.PSet(
228  label = cms.string("jets:step1"),
229  src = cms.InputTag("slimmedJets"),
230  select = cms.string("pt>30 & abs(eta)<4.5 & numberOfDaughters>1 & ((abs(eta)>2.4) || ( chargedHadronEnergyFraction > 0 & chargedMultiplicity>0 & chargedEmEnergyFraction<0.99)) & neutralEmEnergyFraction < 0.99 & neutralHadronEnergyFraction < 0.99"),
231 
232  min = cms.int32(1),
233  max = cms.int32(1),
234 
235  ),
236  cms.PSet(
237  label = cms.string("jets:step2"),
238  src = cms.InputTag("slimmedJets"),
239  select = cms.string("pt>30 & abs(eta)<4.5 & numberOfDaughters>1 & ((abs(eta)>2.4) || ( chargedHadronEnergyFraction > 0 & chargedMultiplicity>0 & chargedEmEnergyFraction<0.99)) & neutralEmEnergyFraction < 0.99 & neutralHadronEnergyFraction < 0.99"),
240 
241  min = cms.int32(2),
242  max = cms.int32(2),
243 
244  ),
245  )
246 )