search

cms-jet / JetToolbox

Python framework for configuration of jet tools via the jet toolbox.
https://twiki.cern.ch/twiki/bin/viewauth/CMS/JetToolbox
7 stars 36 forks source link

option to reuse stored puppi weights #82

Closed kpedro88 closed 2 years ago

kpedro88 commented 2 years ago

I had been under the impression that PUPPI weights weren't stored in the output miniAOD. However, they do appear to be stored in an ultra-legacy MC sample that I checked. Therefore, it can be useful to reuse the weights when clustering the existing packed PFCandidate collection, e.g. to get low-energy AK8 jets that aren't included in miniAOD by default or to get AK15 jets (this saves CPU vs. rerunning the entire PUPPI algorithm). The results vs. rerunning PUPPI are not identical, but the AK8 jet pT distributions are very similar (differences seem to be within the limit of numerical precision).

alefisico commented 2 years ago

Hi, yes I think this is also historical. If I remember correctly we were not storing PUPPI weights at some point in miniAOD and the safest way to compute PUPPI jet collections was to redo the puppi weights. I think this is a standard procedure now, so it is good to include an option to use the weights. I'll merge the PR asap.