Skip to content
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

Fix gen jet variables, add new jet variables #138

Merged
merged 2 commits into from
Nov 18, 2022

Conversation

riga
Copy link

@riga riga commented Mar 11, 2022

This PR fixes the lookup of the gen jet index per pat::Jet, and adds a jet variables which seem helpful.

Fixed indices:

  • Before the index was obtained by comparing pointers, but since the gen jet lookup creates a copy, the pointers will never match any object in the actual gen jet collection.
  • The fix uses a feature of the forward reference, which already stores the index.

New variables:

  • bParticleNetAK4JetTags_probbb
  • bParticleNetAK4JetTags_probpu
  • bParticleNetAK4JetTags_probcc
  • bParticleNetAK4JetTags_probundef
  • bParticleNetAK4JetTags_probc
  • bParticleNetAK4JetTags_probb
  • bParticleNetAK4JetTags_probuds
  • bParticleNetAK4JetTags_probg
  • ak8jets_massIndependentDeepDoubleBvLJetTags_probHbb
  • ak8jets_deepDoubleBvLJetTags_probHbb
  • ak8jets_deepBoostedJetTags_probHbb
  • ak8jets_particleNetJetTags_probHbb
  • ak8jets_particleNetDiscriminatorsJetTags_HbbvsQCD

Edit: Sorry, I saw only now that my editor converted places that had mixed spaces/tabs into just spaces. I can revert that if needed.

@riga riga changed the title Fix gen jet index in HTauTauNtuplizer. Fix gen jet variables, add new jet variables Mar 11, 2022
@dzuolo
Copy link

dzuolo commented Mar 17, 2022

@riga : I assume this PR has been tested with the whole 2017 production right? ;) If there is nothing else you would like to add i propose we merge it. Do you agree @portalesHEP ?

@portalesHEP
Copy link

I was thinking of waiting for some preliminary results to be sure we'll use the new variables, but it's fine for me either way

@riga
Copy link
Author

riga commented Mar 17, 2022

@riga : I assume this PR has been tested with the whole 2017 production right? ;) If there is nothing else you would like to add i propose we merge it. Do you agree @portalesHEP ?

Yes, we ran the whole 2017 bigNtuple production with it, but we haven't looked into their distributions yet. We we're planning to do this with the small ntuples asap, and we bring them into one of the next meetings.

@portalesHEP
Copy link

Following up on this morning's discussions, since we'll very likely want to study ParticleNet for future analyses as well I guess it does not hurt to have the variables included already even if we end up not using it in the current analysis (although I guess the variables and framework may change in the meanwhile). So if nothing else should be added, I'll merge

@dzuolo
Copy link

dzuolo commented Mar 18, 2022

I agree with you Louis, if these variables are ready to be used I propose to add them!

Comment on lines +2382 to +2390
_bdiscr15.push_back(ijet->bDiscriminator("pfParticleNetAK4JetTags:probbb"));
_bdiscr16.push_back(ijet->bDiscriminator("pfParticleNetAK4JetTags:probpu"));
_bdiscr17.push_back(ijet->bDiscriminator("pfParticleNetAK4JetTags:probcc"));
_bdiscr18.push_back(ijet->bDiscriminator("pfParticleNetAK4JetTags:probundef"));
_bdiscr19.push_back(ijet->bDiscriminator("pfParticleNetAK4JetTags:probc"));
_bdiscr20.push_back(ijet->bDiscriminator("pfParticleNetAK4JetTags:probb"));
_bdiscr21.push_back(ijet->bDiscriminator("pfParticleNetAK4JetTags:probuds"));
_bdiscr22.push_back(ijet->bDiscriminator("pfParticleNetAK4JetTags:probg"));

Copy link

@francescobrivio francescobrivio Mar 18, 2022

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Just double checking: why are these called AK4 when they are used for AK8 jets?

Copy link

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I learned in the talk by Raffa at our MIB meeting that ParticleNet can tag also resolved jets ;)

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Ok cool" Thanks Davide! My question was more about the "correctness" of the names of the taggers?

Copy link

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I assume these are the taggers for AK4 jets, the one for AK8 have AK8 in the name...

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I'll let @riga confirm, but these names look correct (and are defined here)

@dzuolo
Copy link

dzuolo commented Mar 18, 2022

One question that came to my mind: for the AK8 jets (maybe there is something similar also for AK4) is the regressed mass also available? I know it has been used in a couple of boosted HH analysis.

@portalesHEP
Copy link

portalesHEP commented Mar 18, 2022

One question that came to my mind: for the AK8 jets (maybe there is something similar also for AK4) is the regressed mass also available? I know it has been used in a couple of boosted HH analysis.

To my understanding, the regressed mass was obtained by "hacking" ParticleNet (removing the output layer iirc, probably in addition to some more changes)
It'd seem that the regression was included (at least) in NanoAODs (cf. cms-sw/cmssw#33588). If it's also available in MiniAODs it could be interesting to look into as well indeed

@dzuolo
Copy link

dzuolo commented Mar 18, 2022

I have the same recollection, I was just proposing to check if something similar was added to miniAOD

@kramerto kramerto mentioned this pull request Nov 14, 2022
@riga
Copy link
Author

riga commented Nov 14, 2022

This has been a while ... 😇

I just checked again and it's does not seem to be added on mini-level.

@portalesHEP
Copy link

Indeed this was left hanging. Can it (/should it) be merged now or are there further tests being considered?

@riga
Copy link
Author

riga commented Nov 16, 2022

I re-checked the event content shipped with the miniAOD files and I think the code in PR should cover everything that we might need, so from our side the PR could be merged 👍

@portalesHEP portalesHEP merged commit c8e36ab into LLRCMS:106X_HH_UL Nov 18, 2022
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

4 participants