Access to jet/quark information

atishelm · July 13, 2023, 11:05am

Dear FCCAnalyses experts,

We are looking for ways to access certain variables related to jets reconstructed using the ExclusiveJetClusteringHelper wrapper. We call this in the following way:

github.com

BNL-FCCee/BNL-Analyses/blob/master/ZccH_stage1.py#L193


      
          # Mandatory: analysers funtion to define the analysers to process, please make sure you return the last dataframe, in this example it is df2
          def analysers(df):
              global jetClusteringHelper
              global jetFlavourHelper
          
              from examples.FCCee.weaver.config import collections
          
              tag = ""
          
              ## define jet clustering parameters
              jetClusteringHelper = ExclusiveJetClusteringHelper(collections["PFParticles"], njets, tag)
          
              ## run jet clustering
              df = jetClusteringHelper.define(df)
          
              ## define jet flavour tagging parameters
              jetFlavourHelper = JetFlavourHelper(
                  collections,
                  jetClusteringHelper.jets,
                  jetClusteringHelper.constituents,
                  tag,

where this eventually allows us to access lots of variables corresponding to the reconstructed jets:

github.com

HEP-FCC/FCCAnalyses/blob/2422e7e4b8cd8f99646f5a53e138b2c712dcfef7/examples/FCCee/weaver/config.py#L349-L436


      
          "jet_p": {
              "name": "jet_p",
              "title": "p_{jet} [GeV]",
              "bin": 100,
              "xmin": 0,
              "xmax": 100.0,
              "scale": "lin",
          },
          "jet_e": {
              "name": "jet_e",
              "title": "E_{jet} [GeV]",
              "bin": 100,
              "xmin": 0,
              "xmax": 100.0,
              "scale": "lin",
          },
          "jet_mass": {
              "name": "jet_mass",
              "title": "m_{jet} [GeV]",
              "bin": 100,

This file has been truncated. show original

Momentum, mass, etc.

There is however information regarding the jets that we were not able to easily find in the framework, which leads us to our questions:

How do we get a reconstructed jet’s radius? We assume this can be dynamic when running in “N-jets” mode, i.e. requiring 4 reconstructed jets in a given event.
How do we check which jet a given PFcand was assigned to?
How do we access per-event truth information for the quarks in an event?

Thank you!
-Abe

cc: @selvaggi @jaeyserm