CMS 3D CMS Logo

PileupJetIDParams_cfi.py
Go to the documentation of this file.
1 import FWCore.ParameterSet.Config as cms
3 
4 ####################################################################################################################
5 full_81x_chs = cms.PSet(
6  impactParTkThreshold = cms.double(1.),
7  cutBased = cms.bool(False),
8  etaBinnedWeights = cms.bool(True),
9  tmvaMethod = cms.string("JetIDMVAHighPt"),
10  version = cms.int32(-1),
11  nEtaBins = cms.int32(4),
12  trainings = cms.VPSet(
13  cms.PSet(
14  jEtaMin = cms.double(0.),
15  jEtaMax = cms.double(2.5),
16  tmvaWeights = cms.FileInPath("RecoJets/JetProducers/data/pileupJetId_80XvarFix_Eta0to2p5_BDT.weights.xml.gz"),
17  tmvaVariables = cms.vstring(
18  "nvtx",
19  "dR2Mean" ,
20  "nParticles" ,
21  "nCharged" ,
22  "majW" ,
23  "minW",
24  "frac01" ,
25  "frac02" ,
26  "frac03" ,
27  "frac04" ,
28  "ptD" ,
29  "beta" ,
30  "pull" ,
31  "jetR" ,
32  "jetRchg" ,
33  )
34  ),
35  cms.PSet(
36  jEtaMin = cms.double(2.5),
37  jEtaMax = cms.double(2.75),
38  tmvaWeights = cms.FileInPath("RecoJets/JetProducers/data/pileupJetId_80XvarFix_Eta2p5to2p75_BDT.weights.xml.gz"),
39  tmvaVariables = cms.vstring(
40  "nvtx",
41  "dR2Mean" ,
42  "nParticles" ,
43  "nCharged" ,
44  "majW" ,
45  "minW",
46  "frac01" ,
47  "frac02" ,
48  "frac03" ,
49  "frac04" ,
50  "ptD" ,
51  "beta" ,
52  "pull" ,
53  "jetR" ,
54  "jetRchg" ,
55  )
56  ),
57  cms.PSet(
58  jEtaMin = cms.double(2.75),
59  jEtaMax = cms.double(3.),
60  tmvaWeights = cms.FileInPath("RecoJets/JetProducers/data/pileupJetId_80XvarFix_Eta2p75to3_BDT.weights.xml.gz"),
61  tmvaVariables = cms.vstring(
62  "nvtx",
63  "dR2Mean" ,
64  "nParticles" ,
65  "nCharged" ,
66  "majW" ,
67  "minW",
68  "frac01" ,
69  "frac02" ,
70  "frac03" ,
71  "frac04" ,
72  "ptD" ,
73  "beta" ,
74  "pull" ,
75  "jetR" ,
76  "jetRchg" ,
77  )
78  ),
79  cms.PSet(
80  jEtaMin = cms.double(3.),
81  jEtaMax = cms.double(5.),
82  tmvaWeights = cms.FileInPath("RecoJets/JetProducers/data/pileupJetId_80XvarFix_Eta3to5_BDT.weights.xml.gz"),
83  tmvaVariables = cms.vstring(
84  "nvtx",
85  "dR2Mean" ,
86  "nParticles" ,
87  "majW" ,
88  "minW",
89  "frac01" ,
90  "frac02" ,
91  "frac03" ,
92  "frac04" ,
93  "ptD" ,
94  "pull" ,
95  "jetR" ,
96  )
97  ),
98  ),
99  tmvaSpectators = cms.vstring(
100  "jetPt" ,
101  "jetEta" ,
102  ),
103  JetIdParams = full_81x_chs_wp,
104  label = cms.string("full")
105 )
106 ####################################################################################################################
107 trainingVariables_102X_Eta0To3 = [
108  "nvtx" ,
109  "beta" ,
110  "dR2Mean" ,
111  "frac01" ,
112  "frac02" ,
113  "frac03" ,
114  "frac04" ,
115  "majW" ,
116  "minW" ,
117  "jetR" ,
118  "jetRchg" ,
119  "nParticles",
120  "nCharged" ,
121  "ptD" ,
122  "pull" ,
123  ]
124 trainingVariables_102X_Eta3To5 = list(trainingVariables_102X_Eta0To3)
125 trainingVariables_102X_Eta3To5.remove('beta')
126 trainingVariables_102X_Eta3To5.remove('jetRchg')
127 trainingVariables_102X_Eta3To5.remove('nCharged')
128 
129 full_102x_chs = full_81x_chs.clone(JetIdParams = full_102x_chs_wp)
130 full_102x_chs.trainings[0].tmvaWeights = "RecoJets/JetProducers/data/pileupJetId_102X_Eta0p0To2p5_chs_BDT.weights.xml.gz"
131 full_102x_chs.trainings[0].tmvaVariables = trainingVariables_102X_Eta0To3
132 full_102x_chs.trainings[1].tmvaWeights = "RecoJets/JetProducers/data/pileupJetId_102X_Eta2p5To2p75_chs_BDT.weights.xml.gz"
133 full_102x_chs.trainings[1].tmvaVariables = trainingVariables_102X_Eta0To3
134 full_102x_chs.trainings[2].tmvaWeights = "RecoJets/JetProducers/data/pileupJetId_102X_Eta2p75To3p0_chs_BDT.weights.xml.gz"
135 full_102x_chs.trainings[2].tmvaVariables = trainingVariables_102X_Eta0To3
136 full_102x_chs.trainings[3].tmvaWeights = "RecoJets/JetProducers/data/pileupJetId_102X_Eta3p0To5p0_chs_BDT.weights.xml.gz"
137 full_102x_chs.trainings[3].tmvaVariables = trainingVariables_102X_Eta3To5
138 
139 ####################################################################################################################
140 trainingVariables_94X_Eta0To3 = list(trainingVariables_102X_Eta0To3)
141 trainingVariables_94X_Eta3To5 = list(trainingVariables_102X_Eta3To5)
142 full_94x_chs = full_81x_chs.clone(JetIdParams = full_94x_chs_wp)
143 full_94x_chs.trainings[0].tmvaWeights = "RecoJets/JetProducers/data/pileupJetId_94X_Eta0p0To2p5_chs_BDT.weights.xml.gz"
144 full_94x_chs.trainings[0].tmvaVariables = trainingVariables_94X_Eta0To3
145 full_94x_chs.trainings[1].tmvaWeights = "RecoJets/JetProducers/data/pileupJetId_94X_Eta2p5To2p75_chs_BDT.weights.xml.gz"
146 full_94x_chs.trainings[1].tmvaVariables = trainingVariables_94X_Eta0To3
147 full_94x_chs.trainings[2].tmvaWeights = "RecoJets/JetProducers/data/pileupJetId_94X_Eta2p75To3p0_chs_BDT.weights.xml.gz"
148 full_94x_chs.trainings[2].tmvaVariables = trainingVariables_94X_Eta0To3
149 full_94x_chs.trainings[3].tmvaWeights = "RecoJets/JetProducers/data/pileupJetId_94X_Eta3p0To5p0_chs_BDT.weights.xml.gz"
150 full_94x_chs.trainings[3].tmvaVariables = trainingVariables_94X_Eta3To5
151 
152 
153 
154 ####################################################################################################################
155 full_80x_chs = cms.PSet(
156  impactParTkThreshold = cms.double(1.),
157  cutBased = cms.bool(False),
158  etaBinnedWeights = cms.bool(True),
159  tmvaMethod = cms.string("JetIDMVAHighPt"),
160  version = cms.int32(-1),
161  nEtaBins = cms.int32(4),
162  trainings = cms.VPSet(
163  cms.PSet(
164  jEtaMin = cms.double(0.),
165  jEtaMax = cms.double(2.5),
166  tmvaWeights = cms.FileInPath("RecoJets/JetProducers/data/pileupJetId_80X_Eta0to2p5_BDT.weights.xml.gz"),
167  tmvaVariables = cms.vstring(
168  "nvtx",
169  "dR2Mean" ,
170  "nParticles" ,
171  "nCharged" ,
172  "majW" ,
173  "minW",
174  "frac01" ,
175  "frac02" ,
176  "frac03" ,
177  "frac04" ,
178  "ptD" ,
179  "beta" ,
180  "pull" ,
181  "jetR" ,
182  "jetRchg" ,
183  )
184  ),
185  cms.PSet(
186  jEtaMin = cms.double(2.5),
187  jEtaMax = cms.double(2.75),
188  tmvaWeights = cms.FileInPath("RecoJets/JetProducers/data/pileupJetId_80X_Eta2p5to2p75_BDT.weights.xml.gz"),
189  tmvaVariables = cms.vstring(
190  "nvtx",
191  "dR2Mean" ,
192  "nParticles" ,
193  "nCharged" ,
194  "majW" ,
195  "minW",
196  "frac01" ,
197  "frac02" ,
198  "frac03" ,
199  "frac04" ,
200  "ptD" ,
201  "beta" ,
202  "pull" ,
203  "jetR" ,
204  "jetRchg" ,
205  )
206  ),
207  cms.PSet(
208  jEtaMin = cms.double(2.75),
209  jEtaMax = cms.double(3.),
210  tmvaWeights = cms.FileInPath("RecoJets/JetProducers/data/pileupJetId_80X_Eta2p75to3_BDT.weights.xml.gz"),
211  tmvaVariables = cms.vstring(
212  "nvtx",
213  "dR2Mean" ,
214  "nParticles" ,
215  "nCharged" ,
216  "majW" ,
217  "minW",
218  "frac01" ,
219  "frac02" ,
220  "frac03" ,
221  "frac04" ,
222  "ptD" ,
223  "beta" ,
224  "pull" ,
225  "jetR" ,
226  "jetRchg" ,
227  )
228  ),
229  cms.PSet(
230  jEtaMin = cms.double(3.),
231  jEtaMax = cms.double(5.),
232  tmvaWeights = cms.FileInPath("RecoJets/JetProducers/data/pileupJetId_80X_Eta3to5_BDT.weights.xml.gz"),
233  tmvaVariables = cms.vstring(
234  "nvtx",
235  "dR2Mean" ,
236  "nParticles" ,
237  "majW" ,
238  "minW",
239  "frac01" ,
240  "frac02" ,
241  "frac03" ,
242  "frac04" ,
243  "ptD" ,
244  "pull" ,
245  "jetR" ,
246  )
247  ),
248  ),
249  tmvaSpectators = cms.vstring(
250  "jetPt" ,
251  "jetEta" ,
252  ),
253  JetIdParams = full_80x_chs_wp,
254  label = cms.string("full")
255 )
256 
257 ####################################################################################################################
258 full_76x_chs = cms.PSet(
259  impactParTkThreshold = cms.double(1.) ,
260  cutBased = cms.bool(False),
261  etaBinnedWeights = cms.bool(True),
262  nEtaBins = cms.int32(4),
263  trainings = cms.VPSet(
264  cms.PSet(
265  jEtaMin = cms.double(0.),
266  jEtaMax = cms.double(2.5),
267  tmvaWeights = cms.FileInPath("RecoJets/JetProducers/data/pileupJetId_76x_Eta0to2p5_BDT.weights.xml.gz"),
268  tmvaVariables = cms.vstring(
269  "nvtx",
270  "dR2Mean" ,
271  "nParticles" ,
272  "nCharged" ,
273  "majW" ,
274  "minW",
275  "frac01" ,
276  "frac02" ,
277  "frac03" ,
278  "frac04" ,
279  "ptD" ,
280  "beta" ,
281  "pull" ,
282  "jetR" ,
283  "jetRchg" ,
284  )
285  ),
286  cms.PSet(
287  jEtaMin = cms.double(2.5),
288  jEtaMax = cms.double(2.75),
289  tmvaWeights = cms.FileInPath("RecoJets/JetProducers/data/pileupJetId_76x_Eta2p5to2p75_BDT.weights.xml.gz"),
290  tmvaVariables = cms.vstring(
291  "nvtx",
292  "dR2Mean" ,
293  "nParticles" ,
294  "nCharged" ,
295  "majW" ,
296  "minW",
297  "frac01" ,
298  "frac02" ,
299  "frac03" ,
300  "frac04" ,
301  "ptD" ,
302  "beta" ,
303  "pull" ,
304  "jetR" ,
305  "jetRchg" ,
306  )
307  ),
308  cms.PSet(
309  jEtaMin = cms.double(2.75),
310  jEtaMax = cms.double(3.),
311  tmvaWeights = cms.FileInPath("RecoJets/JetProducers/data/pileupJetId_76x_Eta2p75to3_BDT.weights.xml.gz"),
312  tmvaVariables = cms.vstring(
313  "nvtx",
314  "dR2Mean" ,
315  "nParticles" ,
316  "nCharged" ,
317  "majW" ,
318  "minW",
319  "frac01" ,
320  "frac02" ,
321  "frac03" ,
322  "frac04" ,
323  "ptD" ,
324  "beta" ,
325  "pull" ,
326  "jetR" ,
327  "jetRchg" ,
328  )
329  ),
330  cms.PSet(
331  jEtaMin = cms.double(3.),
332  jEtaMax = cms.double(5.),
333  tmvaWeights = cms.FileInPath("RecoJets/JetProducers/data/pileupJetId_76x_Eta3to5_BDT.weights.xml.gz"),
334  tmvaVariables = cms.vstring(
335  "nvtx",
336  "dR2Mean" ,
337  "nParticles" ,
338  "majW" ,
339  "minW",
340  "frac01" ,
341  "frac02" ,
342  "frac03" ,
343  "frac04" ,
344  "ptD" ,
345  "pull" ,
346  "jetR" ,
347  )
348  ),
349  ),
350  tmvaMethod = cms.string("JetIDMVAHighPt"),
351  version = cms.int32(-1),
352  tmvaSpectators = cms.vstring(
353  "jetPt" ,
354  "jetEta" ,
355  ),
356  JetIdParams = full_76x_chs_wp,
357  label = cms.string("full")
358  )
359 ####################################################################################################################
360 full_74x_chs = cms.PSet(
361  impactParTkThreshold = cms.double(1.) ,
362  cutBased = cms.bool(False),
363  etaBinnedWeights = cms.bool(True),
364  nEtaBins = cms.int32(4),
365  trainings = cms.VPSet(
366  cms.PSet(
367  jEtaMin = cms.double(0.),
368  jEtaMax = cms.double(2.),
369  tmvaWeights = cms.FileInPath("RecoJets/JetProducers/data/TMVAClassificationCategory_BDTG.weights_jteta_0_2_newNames.xml.gz"),
370  tmvaVariables = cms.vstring(
371  "dR2Mean" ,
372  "rho" ,
373  "nParticles" ,
374  "nCharged" ,
375  "majW" ,
376  "minW",
377  "frac01" ,
378  "frac02" ,
379  "frac03" ,
380  "frac04" ,
381  "ptD" ,
382  "beta" ,
383  "betaStar" ,
384  "pull" ,
385  "jetR" ,
386  "jetRchg" ,
387  )
388  ),
389  cms.PSet(
390  jEtaMin = cms.double(2.),
391  jEtaMax = cms.double(2.5),
392  tmvaWeights = cms.FileInPath("RecoJets/JetProducers/data/TMVAClassificationCategory_BDTG.weights_jteta_2_2p5_newNames.xml.gz"),
393  tmvaVariables = cms.vstring(
394  "dR2Mean" ,
395  "rho" ,
396  "nParticles" ,
397  "nCharged" ,
398  "majW" ,
399  "minW",
400  "frac01" ,
401  "frac02" ,
402  "frac03" ,
403  "frac04" ,
404  "ptD" ,
405  "beta" ,
406  "betaStar" ,
407  "pull" ,
408  "jetR" ,
409  "jetRchg" ,
410  )
411  ),
412  cms.PSet(
413  jEtaMin = cms.double(2.5),
414  jEtaMax = cms.double(3.),
415  tmvaWeights = cms.FileInPath("RecoJets/JetProducers/data/TMVAClassificationCategory_BDTG.weights_jteta_2p5_3_newNames.xml.gz"),
416  tmvaVariables = cms.vstring(
417  "dR2Mean" ,
418  "rho" ,
419  "nParticles" ,
420  "nCharged" ,
421  "majW" ,
422  "minW",
423  "frac01" ,
424  "frac02" ,
425  "frac03" ,
426  "frac04" ,
427  "ptD" ,
428  "beta" ,
429  "betaStar" ,
430  "pull" ,
431  "jetR" ,
432  "jetRchg" ,
433  )
434  ),
435  cms.PSet(
436  jEtaMin = cms.double(3.),
437  jEtaMax = cms.double(5.),
438  tmvaWeights = cms.FileInPath("RecoJets/JetProducers/data/TMVAClassificationCategory_BDTG.weights_jteta_3_5_newNames.xml.gz"),
439  tmvaVariables = cms.vstring(
440  "dR2Mean" ,
441  "rho" ,
442  "nParticles" ,
443  "majW" ,
444  "minW",
445  "frac01" ,
446  "frac02" ,
447  "frac03" ,
448  "frac04" ,
449  "ptD" ,
450  "pull" ,
451  "jetR" ,
452  )
453  ),
454  ),
455  version = cms.int32(-1),
456  tmvaSpectators = cms.vstring(
457  "jetPt" ,
458  "jetEta" ,
459  "nTrueInt" ,
460  "dRMatch" ,
461  ),
462  JetIdParams = full_74x_chs_wp,
463  label = cms.string("full")
464  )
465 ####################################################################################################################
466 full_53x = cms.PSet(
467  impactParTkThreshold = cms.double(1.) ,
468  cutBased = cms.bool(False),
469  etaBinnedWeights = cms.bool(False),
470  tmvaWeights = cms.FileInPath("CondFormats/JetMETObjects/data/TMVAClassificationCategory_JetID_53X_Dec2012.weights.xml"),
471  tmvaMethod = cms.string("JetIDMVAHighPt"),
472  version = cms.int32(-1),
473  tmvaVariables = cms.vstring(
474  "nvtx" ,
475  "dZ" ,
476  "beta" ,
477  "betaStar" ,
478  "nCharged" ,
479  "nNeutrals",
480  "dR2Mean" ,
481  "ptD" ,
482  "frac01" ,
483  "frac02" ,
484  "frac03" ,
485  "frac04" ,
486  "frac05" ,
487  ),
488  tmvaSpectators = cms.vstring(
489  "jetPt",
490  "jetEta",
491  "jetPhi"
492  ),
493  JetIdParams = full_53x_wp,
494  label = cms.string("full53x")
495  )
496 ####################################################################################################################
497 full_53x_chs = cms.PSet(
498  impactParTkThreshold = cms.double(1.) ,
499  cutBased = cms.bool(False),
500  etaBinnedWeights = cms.bool(False),
501  tmvaWeights = cms.FileInPath("CondFormats/JetMETObjects/data/TMVAClassificationCategory_JetID_53X_chs_Dec2012.weights.xml"),
502  #tmvaWeights = cms.FileInPath("RecoJets/JetProducers/data/TMVAClassificationCategory_JetID_53X_chs_Dec2012.weights.xml"),
503  tmvaMethod = cms.string("JetIDMVAHighPt"),
504  version = cms.int32(-1),
505  tmvaVariables = cms.vstring(
506  "nvtx" ,
507  "dZ" ,
508  "beta" ,
509  "betaStar" ,
510  "nCharged" ,
511  "nNeutrals",
512  "dR2Mean" ,
513  "ptD" ,
514  "frac01" ,
515  "frac02" ,
516  "frac03" ,
517  "frac04" ,
518  "frac05" ,
519  ),
520  tmvaSpectators = cms.vstring(
521  "jetPt",
522  "jetEta",
523  "jetPhi"
524  ),
525  JetIdParams = full_53x_chs_wp,
526  label = cms.string("full")
527  )
528 ####################################################################################################################
529 met_53x = cms.PSet(
530  impactParTkThreshold = cms.double(1.) ,
531  cutBased = cms.bool(False),
532  etaBinnedWeights = cms.bool(False),
533  tmvaWeights = cms.FileInPath("RecoJets/JetProducers/data/TMVAClassificationCategory_JetID_MET_53X_Dec2012.weights.xml.gz"),
534  tmvaMethod = cms.string("JetIDMVAMET"),
535  version = cms.int32(-1),
536  tmvaVariables = cms.vstring(
537  "nvtx" ,
538  "jetPt" ,
539  "jetEta" ,
540  "jetPhi" ,
541  "dZ" ,
542  "beta" ,
543  "betaStar" ,
544  "nCharged" ,
545  "nNeutrals",
546  "dR2Mean" ,
547  "ptD" ,
548  "frac01" ,
549  "frac02" ,
550  "frac03" ,
551  "frac04" ,
552  "frac05" ,
553  ),
554  tmvaSpectators = cms.vstring(),
555  JetIdParams = met_53x_wp,
556  label = cms.string("met53x")
557  )
558 ##################################################################################################################
559 full_5x = cms.PSet(
560  impactParTkThreshold = cms.double(1.) ,
561  cutBased = cms.bool(False),
562  etaBinnedWeights = cms.bool(False),
563  tmvaWeights = cms.FileInPath("RecoJets/JetProducers/data/TMVAClassificationCategory_JetID_MET_53X_Dec2012.weights.xml.gz"),
564  tmvaMethod = cms.string("BDT_fullPlusRMS"),
565  version = cms.int32(-1),
566  tmvaVariables = cms.vstring(
567  "frac01",
568  "frac02",
569  "frac03",
570  "frac04",
571  "frac05",
572  "dR2Mean",
573  "nvtx",
574  "nNeutrals",
575  "beta",
576  "betaStar",
577  "dZ",
578  "nCharged",
579  ),
580  tmvaSpectators = cms.vstring(
581  "jetPt",
582  "jetEta",
583  ),
584  JetIdParams = full_5x_wp,
585  label = cms.string("full")
586  )
587 
588 ##################################################################################################################
589 full_5x_chs = cms.PSet(
590  impactParTkThreshold = cms.double(1.) ,
591  cutBased = cms.bool(False),
592  etaBinnedWeights = cms.bool(False),
593  tmvaWeights = cms.FileInPath("RecoJets/JetProducers/data/TMVAClassification_5x_BDT_chsFullPlusRMS.weights.xml.gz"),
594  tmvaMethod = cms.string("BDT_chsFullPlusRMS"),
595  version = cms.int32(-1),
596  tmvaVariables = cms.vstring(
597  "frac01",
598  "frac02",
599  "frac03",
600  "frac04",
601  "frac05",
602  "dR2Mean",
603  "nvtx",
604  "nNeutrals",
605  "beta",
606  "betaStar",
607  "dZ",
608  "nCharged",
609  ),
610  tmvaSpectators = cms.vstring(
611  "jetPt",
612  "jetEta",
613  ),
614  JetIdParams = full_5x_chs_wp,
615  label = cms.string("full")
616  )
617 
618 ####################################################################################################################
619 cutbased = cms.PSet(
620  impactParTkThreshold = cms.double(1.),
621  cutBased = cms.bool(True),
622  JetIdParams = PuJetIdCutBased_wp,
623  label = cms.string("cutbased")
624  )
625 
626 ####################################################################################################################
627 PhilV1 = cms.PSet(
628  impactParTkThreshold = cms.double(1.) ,
629  cutBased = cms.bool(False),
630  etaBinnedWeights = cms.bool(False),
631  tmvaWeights = cms.FileInPath("RecoJets/JetProducers/data/mva_JetID_v1.weights.xml.gz"),
632  tmvaMethod = cms.string("JetID"),
633  version = cms.int32(-1),
634  tmvaVariables = cms.vstring(
635  "nvtx",
636  "jetPt",
637  "jetEta",
638  "jetPhi",
639  "dZ",
640  "d0",
641  "beta",
642  "betaStar",
643  "nCharged",
644  "nNeutrals",
645  "dRMean",
646  "frac01",
647  "frac02",
648  "frac03",
649  "frac04",
650  "frac05",
651  ),
652  tmvaSpectators = cms.vstring(),
653  JetIdParams = JetIdParams,
654  label = cms.string("philv1")
655 )
656 
How EventSelector::AcceptEvent() decides whether to accept an event for output otherwise it is excluding the probing of A single or multiple positive and the trigger will pass if any such matching triggers are PASS or EXCEPTION[A criterion thatmatches no triggers at all is detected and causes a throw.] A single negative with an expectation of appropriate bit checking in the decision and the trigger will pass if any such matching triggers are FAIL or EXCEPTION A wildcarded negative criterion that matches more than one trigger in the trigger list("!*","!HLTx*"if it matches 2 triggers or more) will accept the event if all the matching triggers are FAIL.It will reject the event if any of the triggers are PASS or EXCEPTION(this matches the behavior of"!*"before the partial wildcard feature was incorporated).Triggers which are in the READY state are completely ignored.(READY should never be returned since the trigger paths have been run