Add ParticleNet W vs QCD score

Author: nurfikri89

PhysicsTools/NanoAOD/python/jetsAK8_cff.py

@@ -122,6 +122,7 @@
         particleNet_XttVsQCD = Var("bDiscriminator('pfParticleNetFromMiniAODAK8DiscriminatorsJetTags:HttvsQCD')",float,doc="ParticleNet X->tau_h tau_h vs. QCD score: Xtt/(Xtt+QCD)",precision=10),
         particleNet_XtmVsQCD = Var("bDiscriminator('pfParticleNetFromMiniAODAK8DiscriminatorsJetTags:HtmvsQCD')",float,doc="ParticleNet X->mu tau_h vs. QCD score: Xtm/(Xtm+QCD)",precision=10),
         particleNet_XteVsQCD = Var("bDiscriminator('pfParticleNetFromMiniAODAK8DiscriminatorsJetTags:HtevsQCD')",float,doc="ParticleNet X->e tau_h vs. QCD score: Xte/(Xte+QCD)",precision=10),
+        particleNet_WVsQCD = Var("bDiscriminator('pfParticleNetFromMiniAODAK8DiscriminatorsJetTags:WvsQCD')",float,doc="ParticleNet W->qq vs. QCD score: Xqq+Xcc/(Xqq+Xcc+QCD)",precision=10),
         particleNetLegacy_mass = Var("bDiscriminator('pfParticleNetMassRegressionJetTags:mass')",float,doc="ParticleNet Legacy Run-2 mass regression",precision=10),
         particleNetLegacy_Xbb = Var("bDiscriminator('pfMassDecorrelatedParticleNetJetTags:probXbb')",float,doc="Mass-decorrelated ParticleNet Legacy Run-2 tagger raw X->bb score. For X->bb vs QCD tagging, use Xbb/(Xbb+QCD)",precision=10),
         particleNetLegacy_Xcc = Var("bDiscriminator('pfMassDecorrelatedParticleNetJetTags:probXcc')",float,doc="Mass-decorrelated ParticleNet Legacy Run-2 tagger raw X->cc score. For X->cc vs QCD tagging, use Xcc/(Xcc+QCD)",precision=10),
@@ -174,6 +175,7 @@
     particleNet_XttVsQCD = None,
     particleNet_XtmVsQCD = None,
     particleNet_XteVsQCD = None,
+    particleNet_WVsQCD = None,
     # Remove for V11 and earlier versions
     chMultiplicity = None,
     neMultiplicity = None,

PhysicsTools/NanoAOD/python/nanoDQM_cfi.py

@@ -196,6 +196,7 @@
                 Plot1D('particleNet_XttVsQCD', 'particleNet_XttVsQCD', 20, 0, 1, 'ParticleNet X->tautau vs. QCD score'),
                 Plot1D('particleNet_XtmVsQCD', 'particleNet_XtmVsQCD', 20, 0, 1, 'ParticleNet X->mutau vs. QCD score'),
                 Plot1D('particleNet_XteVsQCD', 'particleNet_XteVsQCD', 20, 0, 1, 'ParticleNet X->etau vs. QCD score'),
+                Plot1D('particleNet_WVsQCD', 'particleNet_WVsQCD', 20, 0, 1, 'ParticleNet W vs. QCD score'),
                 NoPlot('electronIdx3SJ'),
                 Plot1D('eta', 'eta', 20, -4, 4, 'eta'),
                 NoPlot('genJetAK8Idx'),

RecoBTag/ONNXRuntime/python/pfParticleNetFromMiniAODAK8DiscriminatorsJetTags_cfi.py

@@ -87,5 +87,19 @@
             cms.InputTag('pfParticleNetFromMiniAODAK8JetTags', 'probQCD0hf'),
             ),
          ),
+      cms.PSet(
+         name = cms.string('WvsQCD'),
+         numerator = cms.VInputTag(
+            cms.InputTag('pfParticleNetFromMiniAODAK8JetTags', 'probHqq'),
+            cms.InputTag('pfParticleNetFromMiniAODAK8JetTags', 'probHcc'),
+            ),
+         denominator = cms.VInputTag(
+            cms.InputTag('pfParticleNetFromMiniAODAK8JetTags', 'probHqq'),
+            cms.InputTag('pfParticleNetFromMiniAODAK8JetTags', 'probHcc'),
+            cms.InputTag('pfParticleNetFromMiniAODAK8JetTags', 'probQCD2hf'),
+            cms.InputTag('pfParticleNetFromMiniAODAK8JetTags', 'probQCD1hf'),
+            cms.InputTag('pfParticleNetFromMiniAODAK8JetTags', 'probQCD0hf'),
+            ),
+         ),
       )
    )