Jets and Heavy Flavor

• Mission Statement
  • The responsibilities of the Jets and Heavy Flavor working group are to evaluate the performance of the ePIC detector as it relates to jet and heavy flavor observables utilizing a set of performance and physics benchmarks. We also work on developing the analysis techniques and software tools needed to implement the benchmarks used in detector evaluation.
• Conveners
  • Rongrong Ma (marr@bnl.gov)
  • Shyam Kumar (shyam.kumar@ba.infn.it)
• Contact and administrative info
  • Mailing list: eic-projdet-jethf-l@lists.bnl.gov.
  • To subscribe, visit: https://lists.bnl.gov/mailman/listinfo/eic-projdet-jethf-l.
  • Indico: https://indico.bnl.gov/category/420/
  • Mattermost: https://chat.epic-eic.org/main/channels/phys-jets-hf
  • Wiki: https://wiki.bnl.gov/EPIC/index.php?title=JetsHF
• Meetings
  • Tuesdays at 11:30 am ET (bi-weekly)

  ---

Simulations

• Repository for event generation
  • https://github.com/eic/HFsim-PYTHIA
  • https://github.com/eic/HFsim-BeAGLE
• D⁰ samples
  • e+p: PYTHIA v8.306
    • 18x275
      • Each event contains at least one D⁰ that can decay through any channels
      • Location: epic_craterlake/SIDIS/D0_ABCONV/pythia8.306-1.0/18x275/hiDiv
    • 10x100
      • Each event contains at least one D⁰ that decays into π+K pair, and the decay daughters are within pseudorapidity of ±3.5
      • Location:
        • Q² > 1: epic_craterlake/SIDIS/D0_ABCONV/pythia8.306-1.1/10x100/q2_1/hiDiv
        • Q² > 100: epic_craterlake/SIDIS/D0_ABCONV/pythia8.306-1.1/10x100/q2_100/hiDiv
    • 5x41
      • Each event contains at least one D⁰ that can decays into π+K pair
      • Location: epic_craterlake/SIDIS/D0_ABCONV/HFsim-PYTHIA/pythia8.306-1.2/ep/5x41/q2_1to10000/hiDiv
  • e+Au: BeAGLE v1.03.01
    • Each event contains at least one D⁰ that decays into π+K pair
    • Q² > 1
    • Location:
      • 10x100: epic_craterlake_without_zdc/SIDIS/D0_ABCONV/HFsim-BeAGLE/BeAGLE1.03.01-1.0/eAu/10x100/q2_1to10000
      • 5x41: epic_craterlake_without_zdc/SIDIS/D0_ABCONV/HFsim-BeAGLE/BeAGLE1.03.01-1.0/eAu/5x41/q2_1to10000
• L_c samples
  • e+p
    • 18x275
      • Each event contains at least one L_c that can decay through any channels
      • Location: epic_craterlake/SIDIS/Lambda_ABCONV/pythia8.306-1.0/18x275/hiDiv
    • 10x100
      • Each event contains at least one L_c that decays to a π+K+p channel
      • Location: epic_craterlake/SIDIS/Lc_ABCONV/HFsim-PYTHIA/pythia8.306-1.2/ep/10x100/q2_1to10000/hiDiv
  • e+Au: BeAGLE v1.03.01
    • Each event contains at least one L_c that decays through π+K+p channel
    • Q² > 1
    • Location:
      • 10x100: epic_craterlake_without_zdc/SIDIS/Lc_ABCONV/HFsim-BeAGLE/BeAGLE1.03.01-1.0/eAu/10x100/q2_1to10000

Analysis macros

• Jet reader: https://github.com/eic/snippets/tree/main/JetsAndHF/jetReaderExamples
• D⁰ reconstruction based on Helix swimming: https://github.com/eic/snippets/tree/main/JetsAndHF/HF/D0_helix

EICrecon Jet Reconstruction Documentation

There are currently six jet collections in the default output of EICrecon:

• GeneratedJets, which are reconstructed using charged and neutral particles from the GeneratedParticles collection(final-state primaries) using the anti-kT algorithm with R = 1.0;
• GeneratedChargedJets, which are reconstructed using charged particles from the GeneratedParticles collection using the anti-kT algorithm with R = 1.0;
• GeneratedCentauroJets, which are reconstructed using charged and neutral particles from the GeneratedParticles collection using the Centauro algorithm with R = 0.8;
• ReconstructedJets, which are reconstructed using charged and neutral particles (tracks and EMCal clusters) from the ReconstructedParticles collection using the anti-kT algorithm with R = 1.0;
• ReconstructedChargedJets, which are reconstructed using charged particles (tracks) from the ReconstructedChargedParticles collection using the anti-kT algorithm with R = 1.0;
• ReconstructedCentauroJets, which are reconstructed using charged and neutral particles from the ReconstructedParticles collection using the Centauro algorithm with R = 0.8;

Jet reconstruction is handled by FastJet3. User-configurable parameters and their default values are listed in the table below. Jets are stored as edm4eic::ReconstructedParticle objects, and since this data type was intended to describe particles rather than jets, additional information like the jet area currently is not stored..

These parameters can be adjusted at runtime, illustrated by the code snippet below. To use a contributed fastjet algorithm, the user should set jetAlgo to plugin_algorithm.

eicrecon \
  -Preco:ReconstructedChargedJets:jetAlgo=ee_genkt_algorithm \
  -Preco:ReconstructedChargedJets:rJet=0.8 \
  -Preco:ReconstructedChargedJets:pJet=-1.0 \
  <input edm4hep file>