From ba1e7f39c137305c00dcbc6d5627cda62476cc51 Mon Sep 17 00:00:00 2001 From: Nurfikri Norjoharuddeen Date: Thu, 3 Oct 2024 22:50:32 +0200 Subject: [PATCH 1/2] Setup Pileup Jet Id for reclustered AK4 Puppi jets --- PhysicsTools/NanoAOD/python/custom_jme_cff.py | 19 ++++++++++++++++++- 1 file changed, 18 insertions(+), 1 deletion(-) diff --git a/PhysicsTools/NanoAOD/python/custom_jme_cff.py b/PhysicsTools/NanoAOD/python/custom_jme_cff.py index d82138ab7752a..87b0300ebcb26 100644 --- a/PhysicsTools/NanoAOD/python/custom_jme_cff.py +++ b/PhysicsTools/NanoAOD/python/custom_jme_cff.py @@ -662,9 +662,26 @@ def ReclusterAK4PuppiJets(proc, recoJA, runOnMC): proc.jetPuppiTable.variables.phoMultiplicity = PFJETVARS.phoMultiplicity # - # Add variables for pileup jet ID studies. + # Add Pileup Jet ID for Puppi jets # + from RecoJets.JetProducers.PileupJetID_cfi import pileupJetIdPuppi + pileupJetIdName = "pileupJetId{}".format(jetName) + setattr(proc, pileupJetIdName, pileupJetIdPuppi.clone( + jets = "updatedJetsPuppi", + srcConstituentWeights = "packedpuppi", + vertexes = "offlineSlimmedPrimaryVertices", + inputIsCorrected=True, + applyJec=False + ) + ) + proc.jetPuppiTask.add(getattr(proc, pileupJetIdName)) + proc.updatedJetsPuppiWithUserData.userFloats.puIdDisc = cms.InputTag(pileupJetIdName+':fullDiscriminant') + proc.jetPuppiTable.variables.puIdDisc.expr = "userFloat('puIdDisc')" + + # + # Add variables for pileup jet ID studies. + # proc = AddPileUpJetIDVars(proc, jetName = jetName, jetSrc = "updatedJetsPuppi", From 49af8f6914d894998541c589bfe5f4e78e4206e1 Mon Sep 17 00:00:00 2001 From: Nurfikri Norjoharuddeen Date: Fri, 4 Oct 2024 11:40:44 +0200 Subject: [PATCH 2/2] Define the puIdDisc branch for reclustered AK4 Puppi jets --- PhysicsTools/NanoAOD/python/custom_jme_cff.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/PhysicsTools/NanoAOD/python/custom_jme_cff.py b/PhysicsTools/NanoAOD/python/custom_jme_cff.py index 87b0300ebcb26..8c57d9ff5683b 100644 --- a/PhysicsTools/NanoAOD/python/custom_jme_cff.py +++ b/PhysicsTools/NanoAOD/python/custom_jme_cff.py @@ -676,7 +676,7 @@ def ReclusterAK4PuppiJets(proc, recoJA, runOnMC): ) proc.jetPuppiTask.add(getattr(proc, pileupJetIdName)) proc.updatedJetsPuppiWithUserData.userFloats.puIdDisc = cms.InputTag(pileupJetIdName+':fullDiscriminant') - proc.jetPuppiTable.variables.puIdDisc.expr = "userFloat('puIdDisc')" + proc.jetPuppiTable.variables.puIdDisc = Var("userFloat('puIdDisc')", float, doc="Pileup ID BDT discriminant with 133X Winter24 PuppiV18 training",precision=10) #