New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Puppi pileupjetID #24287
Puppi pileupjetID #24287
Conversation
The code-checks are being triggered in jenkins. |
+code-checks Logs: https://cmssdt.cern.ch/SDT/code-checks/cms-sw-PR-24287/6020 |
A new Pull Request was created by @knash for master. It involves the following packages: RecoJets/JetProducers @perrotta, @cmsbuild, @slava77 can you please review it and eventually sign? Thanks. cms-bot commands are listed here |
@cmsbuild please test |
The tests are being triggered in jenkins. |
Comparison job queued. |
Comparison is ready Comparison Summary:
|
} | ||
if((abs(icand->pdgId()) == 1) || (abs(icand->pdgId()) == 2)) multNeut += candPuppiWeight; |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
What is the purpose of this line? PdgId 1,2 seem not to be used in PF? https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuideParticleFlow#Particle_types
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
These are HF had/gamma. The multNeut should reproduce patjet->neutralMultiplicity() exactly in the case that usePuppi is False.
So we have events like this (setting usePuppi to False)
Adding 130
Adding 130
Adding 22
Adding 22
Adding 22
Adding 22
Adding 22
Adding 22
multNeut = 8 ; patjet->neutralMultiplicity() = 8
and events like this
Adding 1
Adding 1
Adding 1
Adding 1
Adding 1
multNeut = 5 ; patjet->neutralMultiplicity() = 5
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
ok makes sense!
@cmsbuild please test |
The tests are being triggered in jenkins. |
@cmsbuild please test |
The tests are being triggered in jenkins. |
Comparison job queued. |
Comparison is ready Comparison Summary:
|
This pull request is fully signed and it will be integrated in one of the next master IBs (tests are also fine). This pull request will now be reviewed by the release team before it's merged. @davidlange6, @slava77, @smuzaffar, @fabiocos (and backports should be raised in the release meeting by the corresponding L2) |
+1 |
Update to pileupjetID algorithm to include an option to calculate the discriminator variables using Puppi weights.
There are a few effects, but one visible change is the fractional pt variables. previously these were calculated using the Puppi jet pt in the denominator and and the unweighted constituent sum in the numerator which leads to an undesirable >1 tail in the total pt fraction. See below
OFF
ON
There was no additional implementation beyond the feature addition (ie no standard workflows should change).
I think the changes in RecoJets/JetProducers/src/PileupJetIdAlgo.cc, Line 509-510 are appropriate but would like a second opinion @ahinzmann .
This was tested only using a modified version of the script:
RecoJets/JetProducers/test/testJetTools_onMiniAOD_cfg.py
There were a few small changes required for this script to run in 10_3_X which I am also pushing.