ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/CMGWPrimeGroup/python/selectors_cff.py
Revision: 1.29
Committed: Mon Aug 5 15:47:42 2013 UTC (11 years, 8 months ago) by fantasia
Content type: text/x-python
Branch: MAIN
CVS Tags: V00-04-10, HEAD
Changes since 1.28: +11 -1 lines
Log Message:
approved wz

File Contents

# Content
1 import FWCore.ParameterSet.Config as cms
2
3 #####################
4 ###### Muons ######
5 #####################
6
7 #Note: These cuts are <= and >=
8 muonSelectors = cms.PSet(
9 VBTF = cms.PSet(
10 minPt = cms.untracked.double(10.),
11 maxEta = cms.untracked.double(2.4),
12 minIsGlobal = cms.untracked.int32(1),
13 maxDxy = cms.untracked.double(0.2),
14 maxNormalizedChi2 = cms.untracked.double(10.0),
15 minNTrackerHits = cms.untracked.int32(11),
16 minNPixelHits = cms.untracked.int32(1),
17 minNMuonHits = cms.untracked.int32(1),
18 minNMatches = cms.untracked.int32(2),
19 ),
20 PFIsoHighPtBoostedZRelaxed = cms.PSet(
21 minPt = cms.untracked.double(10.),
22 maxEta = cms.untracked.double(2.4),
23 minIsGblOrTrk = cms.untracked.int32(1),
24 maxDxy = cms.untracked.double(0.2),
25 maxDz = cms.untracked.double(0.5),
26 maxRelPtErr = cms.untracked.double(0.3),
27 minNPixelHits = cms.untracked.int32(1),
28 minNTrackerLayers = cms.untracked.int32(6),
29 minNMatches = cms.untracked.int32(2),
30 #maxPFIso = cms.untracked.double(???),#Added below in clones
31 ),
32 PFIsoHighPtRelaxed = cms.PSet(
33 minPt = cms.untracked.double(10.),
34 maxEta = cms.untracked.double(2.4),
35 minIsGlobal = cms.untracked.int32(1),
36 maxDxy = cms.untracked.double(0.2),
37 maxDz = cms.untracked.double(0.5),
38 maxRelPtErr = cms.untracked.double(0.3),
39 minNPixelHits = cms.untracked.int32(1),
40 minNMuonHits = cms.untracked.int32(1),
41 minNTrackerLayers = cms.untracked.int32(6),
42 minNMatches = cms.untracked.int32(2),
43 #maxPFIso = cms.untracked.double(???),#Added below in clones
44 ),
45 PFIsoRelaxed = cms.PSet(
46 minPt = cms.untracked.double(10.),
47 maxEta = cms.untracked.double(2.4),
48 minIsGlobal = cms.untracked.int32(1),
49 minIsPF = cms.untracked.int32(1),
50 maxDxy = cms.untracked.double(0.2),
51 maxDz = cms.untracked.double(0.5),
52 maxNormalizedChi2 = cms.untracked.double(10.0),
53 maxRelPtErr = cms.untracked.double(0.3),
54 minNPixelHits = cms.untracked.int32(1),
55 minNMuonHits = cms.untracked.int32(1),
56 minNTrackerLayers = cms.untracked.int32(6),
57 minNMatches = cms.untracked.int32(2),
58 #maxPFIso = cms.untracked.double(???),#Added below in clones
59 ),
60 Reco = cms.PSet(
61 minPt = cms.untracked.double(10.),
62 ),
63 exotica = cms.PSet(
64 maxEta = cms.untracked.double(2.4),
65 minIsGlobal = cms.untracked.int32(1),
66 minIsTracker = cms.untracked.int32(1),
67 maxDxy = cms.untracked.double(0.2),
68 # maxNormalizedChi2 = cms.untracked.double(10.), not needed for TeV muons
69 # minNTrackerHits = cms.untracked.int32(11), deprecated
70 minNPixelHits = cms.untracked.int32(1),
71 minNMuonHits = cms.untracked.int32(1),
72 minNMatches = cms.untracked.int32(2),
73 minNTrackerLayers = cms.untracked.int32(9),
74 # minTrackerValidFrac = cms.untracked.double(),
75 ),
76 )
77 ###WprimeWZ
78 muonSelectors.WZRelaxed = muonSelectors.VBTF.clone(
79 )
80 muonSelectors.WZRelaxedPt20 = muonSelectors.WZRelaxed.clone(
81 minPt = 20.
82 )
83 muonSelectors.WZLoose = muonSelectors.WZRelaxed.clone(
84 maxIso03 = cms.untracked.double(0.15)
85 )
86 muonSelectors.WZTight = muonSelectors.WZRelaxed.clone(
87 maxIso03 = cms.untracked.double(0.1)
88 )
89 muonSelectors.WZTightPt20 = muonSelectors.WZTight.clone(
90 minPt = 20.,
91 )
92 ###EWKWZ
93 muonSelectors.EWKWZRelaxed = muonSelectors.PFIsoRelaxed.clone(
94 )
95 muonSelectors.EWKWZRelaxedPt20 = muonSelectors.EWKWZRelaxed.clone(
96 minPt = 20.
97 )
98 muonSelectors.EWKWZLoose = muonSelectors.EWKWZRelaxed.clone(
99 maxPFIso = cms.untracked.double(0.2)
100 )
101 muonSelectors.EWKWZLoosePt20 = muonSelectors.EWKWZLoose.clone(
102 minPt = 20.,
103 )
104 muonSelectors.EWKWZTight = muonSelectors.EWKWZRelaxed.clone(
105 maxPFIso = cms.untracked.double(0.12)
106 )
107 muonSelectors.EWKWZTightPt20 = muonSelectors.EWKWZTight.clone(
108 minPt = 20.,
109 )
110
111 ###High Pt PfIso Boosted Z
112 muonSelectors.PFIsoHighPtBoostedZRelaxedPt20 = muonSelectors.PFIsoHighPtBoostedZRelaxed.clone(
113 minPt = 20.
114 )
115 muonSelectors.PFIsoHighPtBoostedZLoose = muonSelectors.PFIsoHighPtBoostedZRelaxed.clone(
116 maxPFIso = cms.untracked.double(0.2)
117 )
118 muonSelectors.PFIsoHighPtBoostedZLoosePt20 = muonSelectors.PFIsoHighPtBoostedZLoose.clone(
119 minPt = 20.,
120 )
121 muonSelectors.PFIsoHighPtBoostedZTight = muonSelectors.PFIsoHighPtBoostedZRelaxed.clone(
122 maxPFIso = cms.untracked.double(0.12)
123 )
124 muonSelectors.PFIsoHighPtBoostedZTightPt20 = muonSelectors.PFIsoHighPtBoostedZTight.clone(
125 minPt = 20.,
126 )
127
128 #####High Pt PfIso
129 muonSelectors.PFIsoHighPtRelaxedPt20 = muonSelectors.PFIsoHighPtRelaxed.clone(
130 minPt = 20.
131 )
132 muonSelectors.PFIsoHighPtLoose = muonSelectors.PFIsoHighPtRelaxed.clone(
133 maxPFIso = cms.untracked.double(0.2)
134 )
135 muonSelectors.PFIsoHighPtLoosePt20 = muonSelectors.PFIsoHighPtLoose.clone(
136 minPt = 20.,
137 )
138 muonSelectors.PFIsoHighPtTight = muonSelectors.PFIsoHighPtRelaxed.clone(
139 maxPFIso = cms.untracked.double(0.12)
140 )
141 muonSelectors.PFIsoHighPtTightPt20 = muonSelectors.PFIsoHighPtTight.clone(
142 minPt = 20.,
143 )
144
145
146 ###HadVZ
147 muonSelectors.HadVZLoose = muonSelectors.exotica.clone(
148 minPt = cms.untracked.double(20.),
149 # maxIso03 = cms.untracked.double(0.1)
150 )
151 muonSelectors.HadVZTight = muonSelectors.exotica.clone(
152 minPt = cms.untracked.double(35.),
153 # maxIso03 = cms.untracked.double(0.1)
154 )
155
156 #####################
157 #### Electrons ####
158 #####################
159
160 cutsMissingHits = [0, 0, 0, 0, 0, 0]
161 cutsConvDist = [0., 0., 0.02, 0.02, 0.02, 0.02]
162 cutsConvDcot = [0., 0., 0.02, 0.02, 0.02, 0.02]
163 cutsEBCombRelIso = [0.150, 0.085, 0.053, 0.070, 0.030, 0.016]#2010
164 #cutsEBCombRelIso = [0.150, 0.085, 0.053, 0.040, 0.030, 0.016]#2011
165 cutsEBSigmaIEtaIEta = [0.012, 0.01, 0.01, 0.01, 0.01, 0.01]
166 cutsEBDeltaPhi = [0.800, 0.071, 0.039, 0.027, 0.020, 0.020]
167 cutsEBDeltaEta = [0.007, 0.007, 0.005, 0.005, 0.004, 0.004]
168 #cutsEECombRelIso = [0.100, 0.051, 0.042, 0.033, 0.016, 0.008]#2011
169 cutsEECombRelIso = [0.100, 0.051, 0.042, 0.060, 0.016, 0.008]#2010
170 cutsEESigmaIEtaIEta = [0.031, 0.031, 0.031, 0.031, 0.031, 0.031]
171 cutsEEDeltaPhi = [0.7, 0.047, 0.028, 0.021, 0.021, 0.021]
172 cutsEEDeltaEta = [0.011, 0.011, 0.007, 0.006, 0.005, 0.004]
173
174 electronSelectors = cms.PSet(
175 WZLoose = cms.PSet(),
176 WZRelaxed95 = cms.PSet(),
177 WZRelaxed80 = cms.PSet(),
178 WZTight = cms.PSet(),
179 exotica = cms.PSet(
180 joint = cms.PSet(
181 minPt = cms.untracked.double(10.),
182 minIsEcalDriven = cms.untracked.int32(1),
183 maxDeltaPhi = cms.untracked.double(0.06),
184 maxHoverE = cms.untracked.double(0.05),
185 #maxTrackIso = cms.untracked.double(5.),
186 maxPFIso = cms.untracked.double(0.15),
187 maxMissingHits = cms.untracked.int32(1),
188 ),
189 barrel = cms.PSet(
190 maxDeltaEta = cms.untracked.double(0.005),
191 #maxSigmaIEtaIEta = cms.untracked.double(),
192 minPassEX5overE55 = cms.untracked.int32(1),
193 #minPassEMHadDepth1Iso = cms.untracked.int32(1),
194 maxd0 = cms.untracked.double(0.02),
195 ),
196 endcap = cms.PSet(
197 maxDeltaEta = cms.untracked.double(0.007),
198 maxSigmaIEtaIEta = cms.untracked.double(0.03),
199 minPassEX5overE55 = cms.untracked.int32(1),
200 #minPassEMHadDepth1Iso = cms.untracked.int32(1),
201 maxd0 = cms.untracked.double(0.05),
202 ),
203 ),
204 mva = cms.PSet(
205 joint = cms.PSet(
206 minPt = cms.untracked.double(10.),
207 minpassMVAPresel = cms.untracked.int32(1),
208 minpassMVATrig = cms.untracked.int32(1),
209 #maxPFIso = cms.untracked.double(???),#Added below in clones
210 ),
211 barrel = cms.PSet(
212 ),
213 endcap = cms.PSet(
214 ),
215 ),
216 EWKWZLoose = cms.PSet(),
217 EWKWZTight = cms.PSet(),
218
219 CiC2012VetoRelaxed = cms.PSet(
220 joint = cms.PSet(
221 minPt = cms.untracked.double(10.),
222 maxIsGap = cms.untracked.int32(0),
223 minPassConvVeto = cms.untracked.int32(1),
224 maxIsSCGap = cms.untracked.int32(0),
225 ),
226 barrel = cms.PSet(
227 #maxMissingHits = cms.untracked.int32(1),
228 #minVtxFitProb = cms.untracked.double(1e-6),
229 maxSigmaIEtaIEta = cms.untracked.double(0.01),
230 maxDeltaPhi = cms.untracked.double(0.8),
231 maxDeltaEta = cms.untracked.double(0.007),
232 maxHoverE = cms.untracked.double(0.15),
233 maxd0 = cms.untracked.double(0.04),
234 maxdz = cms.untracked.double(0.2),
235 #maxfabsdiffEp = cms.untracked.double(0.05),
236 ),
237 endcap = cms.PSet(
238 #maxMissingHits = cms.untracked.int32(1),
239 #minVtxFitProb = cms.untracked.double(1e-6),
240 maxSigmaIEtaIEta = cms.untracked.double(0.03),
241 maxDeltaPhi = cms.untracked.double(0.7),
242 maxDeltaEta = cms.untracked.double(0.01),
243 #maxHoverE = cms.untracked.double(0.10),
244 maxd0 = cms.untracked.double(0.04),
245 maxdz = cms.untracked.double(0.2),
246 #maxfabsdiffEp = cms.untracked.double(0.05),
247 ),
248 ),
249 CiC2012LooseRelaxed = cms.PSet(
250 joint = cms.PSet(
251 minPt = cms.untracked.double(10.),
252 maxIsGap = cms.untracked.int32(0),
253 minPassConvVeto = cms.untracked.int32(1),
254 maxIsSCGap = cms.untracked.int32(0),
255 ),
256 barrel = cms.PSet(
257 maxMissingHits = cms.untracked.int32(1),
258 minVtxFitProb = cms.untracked.double(1e-6),
259 maxSigmaIEtaIEta = cms.untracked.double(0.01),
260 maxDeltaPhi = cms.untracked.double(0.15),
261 maxDeltaEta = cms.untracked.double(0.007),
262 maxHoverE = cms.untracked.double(0.12),
263 maxd0 = cms.untracked.double(0.02),
264 maxdz = cms.untracked.double(0.2),
265 maxfabsdiffEp = cms.untracked.double(0.05),
266 ),
267 endcap = cms.PSet(
268 maxMissingHits = cms.untracked.int32(1),
269 minVtxFitProb = cms.untracked.double(1e-6),
270 maxSigmaIEtaIEta = cms.untracked.double(0.03),
271 maxDeltaPhi = cms.untracked.double(0.10),
272 maxDeltaEta = cms.untracked.double(0.009),
273 maxHoverE = cms.untracked.double(0.10),
274 maxd0 = cms.untracked.double(0.02),
275 maxdz = cms.untracked.double(0.2),
276 maxfabsdiffEp = cms.untracked.double(0.05),
277 ),
278 ),
279 Reco = cms.PSet(
280 joint = cms.PSet(
281 minPt = cms.untracked.double(10.),
282 maxIsGap = cms.untracked.int32(0),
283 minPassConvVeto = cms.untracked.int32(1),
284 maxIsSCGap = cms.untracked.int32(0),
285 ),
286 barrel = cms.PSet(),
287 endcap = cms.PSet(),
288 ),
289 CiC2012Loose = cms.PSet(),
290 CiC2012LoosePt20 = cms.PSet(),
291 CiC2012MediumRelaxed = cms.PSet(
292 joint = cms.PSet(
293 minPt = cms.untracked.double(10.),
294 maxIsGap = cms.untracked.int32(0),
295 minPassConvVeto = cms.untracked.int32(1),
296 maxIsSCGap = cms.untracked.int32(0),
297 ),
298 barrel = cms.PSet(
299 maxMissingHits = cms.untracked.int32(1),
300 minVtxFitProb = cms.untracked.double(1e-6),
301 maxSigmaIEtaIEta = cms.untracked.double(0.01),
302 maxDeltaPhi = cms.untracked.double(0.06),
303 maxDeltaEta = cms.untracked.double(0.004),
304 maxHoverE = cms.untracked.double(0.12),
305 maxd0 = cms.untracked.double(0.02),
306 maxdz = cms.untracked.double(0.1),
307 maxfabsdiffEp = cms.untracked.double(0.05),
308 ),
309 endcap = cms.PSet(
310 maxMissingHits = cms.untracked.int32(1),
311 minVtxFitProb = cms.untracked.double(1e-6),
312 maxSigmaIEtaIEta = cms.untracked.double(0.03),
313 maxDeltaPhi = cms.untracked.double(0.03),
314 maxDeltaEta = cms.untracked.double(0.007),
315 maxHoverE = cms.untracked.double(0.10),
316 maxd0 = cms.untracked.double(0.02),
317 maxdz = cms.untracked.double(0.1),
318 maxfabsdiffEp = cms.untracked.double(0.05),
319 ),
320 ),
321 CiC2012Medium = cms.PSet(),
322 CiC2012MediumPt20 = cms.PSet(),
323 )
324
325 for i, s in enumerate(["wp95", "wp90", "wp85", "wp80", "wp70", "wp60"]):
326 pset = cms.PSet(
327 joint = cms.PSet(
328 minPt = cms.untracked.double(10.),
329 maxIsSCGap = cms.untracked.int32(0),
330 ),
331 barrel = cms.PSet(
332 maxMissingHits = cms.untracked.int32(cutsMissingHits[i]),
333 minConv = cms.untracked.double(cutsConvDist[i]),#Hack bc we need an OR of these two cuts below
334 maxCombRelIso = cms.untracked.double(cutsEBCombRelIso[i]),
335 maxSigmaIEtaIEta = cms.untracked.double(cutsEBSigmaIEtaIEta[i]),
336 maxDeltaPhi = cms.untracked.double(cutsEBDeltaPhi[i]),
337 maxDeltaEta = cms.untracked.double(cutsEBDeltaEta[i]),
338 ),
339 endcap = cms.PSet(
340 maxMissingHits = cms.untracked.int32(cutsMissingHits[i]),
341 minConv = cms.untracked.double(cutsConvDist[i]),#Hack bc we need an or of these two cuts below
342 maxCombRelIso = cms.untracked.double(cutsEECombRelIso[i]),
343 maxSigmaIEtaIEta = cms.untracked.double(cutsEESigmaIEtaIEta[i]),
344 maxDeltaPhi = cms.untracked.double(cutsEEDeltaPhi[i]),
345 maxDeltaEta = cms.untracked.double(cutsEEDeltaEta[i]),
346 ),
347 )
348 setattr(electronSelectors, s, pset)
349
350 electronSelectors.WZLoose = electronSelectors.wp95.clone()
351
352 electronSelectors.WZRelaxed95 = electronSelectors.wp95.clone()
353 electronSelectors.WZRelaxed95.barrel.maxCombRelIso = cms.untracked.double(999999999.)
354 electronSelectors.WZRelaxed95.endcap.maxCombRelIso = cms.untracked.double(999999999.)
355
356 electronSelectors.WZRelaxed80 = electronSelectors.wp80.clone()
357 electronSelectors.WZRelaxed80.barrel.maxCombRelIso = cms.untracked.double(999999999.)
358 electronSelectors.WZRelaxed80.endcap.maxCombRelIso = cms.untracked.double(999999999.)
359
360 electronSelectors.WZRelaxed80Pt20 = electronSelectors.wp80.clone()
361 electronSelectors.WZRelaxed80Pt20.joint.minPt = cms.untracked.double(20.)
362 electronSelectors.WZRelaxed80Pt20.barrel.maxCombRelIso = cms.untracked.double(999999999.)
363 electronSelectors.WZRelaxed80Pt20.endcap.maxCombRelIso = cms.untracked.double(999999999.)
364
365 electronSelectors.WZTight = electronSelectors.wp80.clone()
366
367 electronSelectors.WZTightPt20 = electronSelectors.wp80.clone()
368 electronSelectors.WZTightPt20.joint.minPt = cms.untracked.double(20.)
369
370 ##EWKWZ
371 electronSelectors.EWKWZLoose = electronSelectors.mva.clone()
372 electronSelectors.EWKWZLoose.joint.maxPFIso = cms.untracked.double(0.1)
373
374 electronSelectors.EWKWZRelaxed = electronSelectors.mva.clone()
375
376 electronSelectors.EWKWZRelaxedPt20 = electronSelectors.mva.clone()
377 electronSelectors.EWKWZRelaxedPt20.joint.minPt = 20.
378
379 electronSelectors.EWKWZTight = electronSelectors.mva.clone()
380 electronSelectors.EWKWZTight.joint.maxPFIso = cms.untracked.double(0.1)
381
382 electronSelectors.EWKWZTightPt20 = electronSelectors.mva.clone()
383 electronSelectors.EWKWZTightPt20.joint.minPt = 20.
384 electronSelectors.EWKWZTightPt20.joint.maxPFIso = cms.untracked.double(0.1)
385
386 #CiC2012
387 electronSelectors.CiC2012Veto = electronSelectors.CiC2012VetoRelaxed.clone()
388 electronSelectors.CiC2012Veto.joint.maxPFIso = cms.untracked.double(0.15)
389
390 electronSelectors.CiC2012VetoPt20 = electronSelectors.CiC2012Veto.clone()
391 electronSelectors.CiC2012VetoPt20.joint.minPt = 20.
392
393 electronSelectors.CiC2012Loose = electronSelectors.CiC2012LooseRelaxed.clone()
394 electronSelectors.CiC2012Loose.joint.maxPFIso = cms.untracked.double(0.15)
395
396 electronSelectors.CiC2012LoosePt20 = electronSelectors.CiC2012Loose.clone()
397 electronSelectors.CiC2012LoosePt20.joint.minPt = 20.
398
399 electronSelectors.CiC2012MediumRelaxedPt20 = electronSelectors.CiC2012MediumRelaxed.clone()
400 electronSelectors.CiC2012MediumRelaxedPt20.joint.minPt = 20.
401
402 electronSelectors.CiC2012Medium = electronSelectors.CiC2012MediumRelaxed.clone()
403 electronSelectors.CiC2012Medium.joint.maxPFIso = cms.untracked.double(0.15)
404 #electronSelectors.CiC2012Medium.joint.maxTrackIso = cms.untracked.double(2)
405
406 electronSelectors.CiC2012MediumPt20 = electronSelectors.CiC2012Medium.clone()
407 electronSelectors.CiC2012MediumPt20.joint.minPt = 20.
408
409 ##Exotica
410 electronSelectors.exoticaPt20 = electronSelectors.exotica.clone()
411 electronSelectors.exoticaPt20.joint.minPt = cms.untracked.double(20.)
412
413 ####################
414 ##### Jets #######
415 ####################
416 jetSelectors = cms.PSet(
417 Base = cms.PSet(
418 minPt = cms.untracked.double(30.),
419 maxEta = cms.untracked.double(2.4),
420 maxNHF = cms.untracked.double(0.99),
421 maxNEF = cms.untracked.double(0.99),
422 minNDaughters = cms.untracked.int32(1),
423 minCHF = cms.untracked.double(0.0),
424 maxCEF = cms.untracked.double(0.99),
425 minCMult = cms.untracked.int32(0),
426 ),
427 Pat = cms.PSet(
428 minPt = cms.untracked.double(30.),
429 maxEta = cms.untracked.double(3.0),
430 ),
431 )
432
433 ####################
434 #### Photons #######
435 ####################
436 photonSelectors = cms.PSet(
437 Base = cms.PSet(#Made up values
438 minPt = cms.untracked.double(10.),
439 maxEta = cms.untracked.double(5),
440 maxECalIso = cms.untracked.double(0.1),
441 maxHCalIso = cms.untracked.double(0.1),
442 maxTrkIso = cms.untracked.double(0.1),
443 maxHoE = cms.untracked.double(0.05),
444 maxSigmaee = cms.untracked.double(9e9),
445 minHasSeed = cms.untracked.bool(True),
446 ),
447 )
448
449 #CiC2012Relaxed