Merge pull request #14043 from jdolen/puid_80x_try2

Update PUID to 80X training
cms-sw · Apr 18, 2016 · c9598ab · c9598ab
2 parents 2694fd6 + a8e7289
commit c9598ab
Show file tree

Hide file tree

Showing 3 changed files with 80 additions and 1 deletion.
diff --git a/RecoJets/JetProducers/python/PileupJetIDCutParams_cfi.py b/RecoJets/JetProducers/python/PileupJetIDCutParams_cfi.py
@@ -1,5 +1,31 @@
 import FWCore.ParameterSet.Config as cms
 
+###########################################################
+## Working points for the 80X training
+###########################################################
+full_80x_chs_wp  = cms.PSet(
+    #4 Eta Categories  0-2.5 2.5-2.75 2.75-3.0 3.0-5.0
+
+    #Tight Id            
+    Pt010_Tight    = cms.vdouble( 0.26, -0.34, -0.24, -0.26),
+    Pt1020_Tight   = cms.vdouble( 0.26, -0.34, -0.24, -0.26),
+    Pt2030_Tight   = cms.vdouble( 0.26, -0.34, -0.24, -0.26),
+    Pt3050_Tight   = cms.vdouble( 0.62, -0.21, -0.07, -0.03),
+
+    #Medium Id
+    Pt010_Medium   = cms.vdouble(-0.49, -0.53, -0.44, -0.42),
+    Pt1020_Medium  = cms.vdouble(-0.49, -0.53, -0.44, -0.42),
+    Pt2030_Medium  = cms.vdouble(-0.49, -0.53, -0.44, -0.42),
+    Pt3050_Medium  = cms.vdouble(-0.06, -0.42, -0.3 , -0.23),
+
+    #Loose Id
+    Pt010_Loose    = cms.vdouble(-0.96, -0.64, -0.56, -0.54),
+    Pt1020_Loose   = cms.vdouble(-0.96, -0.64, -0.56, -0.54),
+    Pt2030_Loose   = cms.vdouble(-0.96, -0.64, -0.56, -0.54),
+    Pt3050_Loose   = cms.vdouble(-0.92, -0.56, -0.44, -0.39)
+)
+
+
 ###########################################################
 ## Working points for the 76X training
 ###########################################################

diff --git a/RecoJets/JetProducers/python/PileupJetIDParams_cfi.py b/RecoJets/JetProducers/python/PileupJetIDParams_cfi.py
@@ -1,5 +1,57 @@
 import FWCore.ParameterSet.Config as cms
 from RecoJets.JetProducers.PileupJetIDCutParams_cfi import *
+
+
+####################################################################################################################                                                                                      
+full_80x_chs = cms.PSet(
+ impactParTkThreshold = cms.double(1.) ,
+ cutBased = cms.bool(False),
+ etaBinnedWeights = cms.bool(True),
+ tmvaWeights_jteta_0_2p5    = cms.string("RecoJets/JetProducers/data/pileupJetId_80X_Eta0to2p5_BDT.weights.xml.gz"),
+ tmvaWeights_jteta_2p5_2p75 = cms.string("RecoJets/JetProducers/data/pileupJetId_80X_Eta2p5to2p75_BDT.weights.xml.gz"),
+ tmvaWeights_jteta_2p75_3   = cms.string("RecoJets/JetProducers/data/pileupJetId_80X_Eta2p75to3_BDT.weights.xml.gz"),
+ tmvaWeights_jteta_3_5      = cms.string("RecoJets/JetProducers/data/pileupJetId_80X_Eta3to5_BDT.weights.xml.gz"),
+ tmvaMethod  = cms.string("JetIDMVAHighPt"),
+ version = cms.int32(-1),
+ tmvaVariables_jteta_0_3 = cms.vstring(
+    "nvtx",
+    "dR2Mean"     ,
+    "nParticles"     ,
+    "nCharged" ,
+    "majW" ,
+    "minW",
+    "frac01"  ,
+    "frac02"      ,
+    "frac03"   ,
+    "frac04"   ,
+    "ptD"   ,
+    "beta"   ,
+    "pull"   ,
+    "jetR"   ,
+    "jetRchg"   ,
+    ),
+ tmvaVariables_jteta_3_5 = cms.vstring(
+    "nvtx",
+    "dR2Mean"     ,
+    "nParticles"     ,
+    "majW" ,
+    "minW",
+    "frac01"  ,
+    "frac02"      ,
+    "frac03"   ,
+    "frac04"   ,
+    "ptD"   ,
+    "pull"   ,
+    "jetR"   ,
+    ),
+ tmvaSpectators = cms.vstring(
+    "jetPt"   ,
+    "jetEta"   ,
+    ),
+ JetIdParams = full_80x_chs_wp,
+ label = cms.string("full")
+ )
+
 ####################################################################################################################                                                                                      
 full_76x_chs = cms.PSet(
  impactParTkThreshold = cms.double(1.) ,

diff --git a/RecoJets/JetProducers/python/PileupJetID_cfi.py b/RecoJets/JetProducers/python/PileupJetID_cfi.py
@@ -9,8 +9,9 @@
 _chsalgos_5x = cms.VPSet(full_5x_chs,cutbased)
 _chsalgos_74x = cms.VPSet(full_74x_chs,cutbased)
 _chsalgos_76x = cms.VPSet(full_76x_chs,cutbased)
+_chsalgos_80x = cms.VPSet(full_80x_chs,cutbased)
 
-_stdalgos    = _chsalgos_76x
+_stdalgos    = _chsalgos_80x
 
 # Calculate+store variables and run MVAs
 pileupJetId = cms.EDProducer('PileupJetIdProducer',