Skip to content

Something significantly changed in forward Lambda decays between 2025-10 and 2026-02 #1069

@DraTeots

Description

@DraTeots

I see many discrepancies in Lambdas on DD4Hep on MCParticle level between versions (using the same hepmc3 inputs). In particular:

  • We see no N + pi0 decays fully saved in MCParticles
  • We see kind of Z cut on lambda decay. Basically lambdas stop decaying before B0.

Environment: (where does this bug occur, have you tried other environments)

  • Which branch (main for latest released), 2026-02 to nightly affected. 2025-10 -better

2026-02

Image Image

2025-10.3
Image

**We compared 3 datasets now: **

  • Our own campaign 2025-10 using 2025-10.3 eic-shell container
  • First official campaign for meson-structure events of 2026-02
  • Our own campaign 2026-03 with 2026-03 image and Baptiste branch using multiple FF detectors to reconstruct lambdas (gives the same result as official 2026-02)

CLARIFICATION:

  • Lambdas we look at are decayed by Geant4/DD4Hep. And we saw the classic decay exponents at some point (decay point z).
  • I look now at both edm4hep (if available) and eicrecon outputs. In terms of MCParticles they are identical as expected. So it is on DD4Hep side.
  • I looked at "decay mode" as looking through MCParticles primary Lambda daughters via EDM4HEP and EDM4EIC C++ reader. Previously it gave 4 variants for lambdas
  1. P + Pi- (as in PDG)
  2. N + Pi0 (as in PDG)
  3. Shower (when you see like 40 daughters). Were the dominating cause is Lambda Recharge - a classic process where in matter Lambda goes Sigma, which instantly decayed back to Lambda and so on. But different hadronic showering is seen.
  4. Nothing. It just not decayed OR decay products are not written. And DD4Hep (DD4Hep governs saving to EDM4Hep) doesn't write all shower particles. So it goes as "Undecayed"

Now something is screwed up even on MCParticles level.

Since 2025-02 looking where n+pi0 gone on event by event basis we see:

"Other" lambda decay mode at evt_id=1513 daughters=1
id=5 pdg=2112

So only neutron is saved and never pi0. (even if only interacting particle is saved we should see it, as previsouly we have 20% efficiency of registration).

Steps to reproduce: (give a step by step account of how to trigger the bug)

The reproducer

https://github.com/JeffersonLab/meson-structure/tree/main/analysis/edm4hep_mcpart_lambda

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions