⚛️ high-energy experiments

Forward decay chain reconstruction with displaced vertices in intense magnetic fields

This paper presents a C++ framework for reconstructing complex particle decay chains with displaced vertices in intense, inhomogeneous magnetic fields, demonstrating significant improvements in efficiency and resolution over existing tools for precision flavour physics and spectroscopy.

Original authors: Sergio Javier Jaimes Elles, Simone Libralon, Fernando Martinez Vidal, Nicola Neri, Miguel Rebollo de Miguel, Izaac Sanderswood, Volodymyr Svintozelskyi, Giorgia Tonani, Javier Velilla Serna, Ziyi Wang

Published 2026-10-14

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Sergio Javier Jaimes Elles, Simone Libralon, Fernando Martinez Vidal, Nicola Neri, Miguel Rebollo de Miguel, Izaac Sanderswood, Volodymyr Svintozelskyi, Giorgia Tonani, Javier Velilla Serna, Ziyi Wang

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the subatomic world, particles are rarely still; they are born in collisions, travel for a fleeting instant, and then transform into other particles. To understand the fundamental rules of nature, physicists must reconstruct these transformations, tracing the path of every fragment back to its origin. This is like trying to solve a three-dimensional puzzle where the pieces are moving at nearly the speed of light and the picture is obscured by a powerful, twisting force. In many experiments, scientists rely on standard tools to map these paths, assuming that the particles travel in straight lines or simple curves. However, this assumption breaks down when particles travel through regions where the magnetic field is incredibly strong and changes rapidly, or when they travel for several meters before they are measured. In these difficult conditions, the standard tools struggle to find the true starting point of a decay, often losing precision or failing to connect the pieces correctly.

A team of researchers at CERN, working with the LHCb experiment, has developed a new software framework designed to solve exactly this problem. They created a system called ForwardTreeFitter, which is built to handle the most challenging decay chains: those where particles travel long distances through intense, uneven magnetic fields before they break apart. The researchers tested this new tool using simulated data from the Large Hadron Collider, focusing on specific particles known as b-hadrons and c-hadrons that decay into other particles several meters away from the collision point. Their work shows that by treating the particle paths with a more sophisticated mathematical approach, they can reconstruct these events with far greater accuracy than current methods allow.

The core difficulty the team addressed is how to track a particle that travels through a magnetic field that is not uniform. In the LHCb detector, the magnetic field is strongest and most irregular in a specific region several meters away from the collision point. When a charged particle enters this zone, its path curves in complex ways that simple straight-line approximations cannot describe. If a particle travels for several meters through this field before it is detected, and then decays into other particles that also travel through the field, the standard tools used by physicists often produce inaccurate results. They might guess the wrong starting point for the decay or miscalculate the energy of the particles involved. This is particularly problematic for particles that live longer than average, as they travel further into these difficult magnetic regions before they break apart.

To fix this, the researchers changed the fundamental way the computer calculates the particle's path. Instead of trying to move the measurement data from the detector back to the decay point using simple geometry, their new system keeps the measurements fixed where they were recorded. It then calculates what the particle's path must have been by running a detailed simulation of the particle moving through the magnetic field to reach those measurement points. This process is repeated many times, refining the estimate of the decay point and the particle's energy until the solution is stable. This approach allows the software to account for the exact shape of the magnetic field and the specific way the particle curves, rather than relying on a simplified guess.

The team tested their new framework against the standard tools used in the LHCb experiment, using simulated examples of particles decaying into a mix of protons, pions, and muons. They looked at two specific scenarios: one where a particle decays between three and six meters from the collision point, and another where it decays between six and eight meters. In the region between three and six meters, where the magnetic field is at its peak intensity, the new system showed a dramatic improvement. For one type of particle decay, the ability to pinpoint the exact location of the decay improved by a factor of nearly four compared to the old method. For another type of decay, the improvement was even more significant, with the new system being more than seven times better at locating the vertex. The new tool also proved much more reliable at reconstructing the total mass of the original particle, a key measurement for identifying what kind of particle was involved.

Perhaps just as important as the improved precision was the new system's ability to find these events in the first place. The standard tools often failed to reconstruct decays that happened in the middle of the magnetic field, losing efficiency as the distance increased. The new framework maintained a high success rate across the entire range, finding about 96 percent of the events even in the most difficult conditions. This means that physicists will be able to study a much larger number of rare events that were previously too hard to analyze. The researchers also found that the new system reduced the number of large errors, where the calculated position was wildly different from the true position, which had been a persistent problem with the older tools.

The success of this work relies on a combination of advanced mathematics and careful engineering. The software is modular, meaning it can be used as a complete package or its parts can be mixed with other tools. It is designed to be flexible enough to work with different types of detectors, not just the one at CERN, but it was specifically tuned to handle the unique challenges of the LHCb experiment. By using a method that treats the particle's journey as a continuous, non-linear path through a complex field, the researchers have removed a major bottleneck in their ability to study the properties of matter. The results suggest that for experiments involving long-lived particles and strong magnetic fields, the old assumption of simple paths is no longer sufficient, and a more rigorous, simulation-based approach is necessary to see the truth of the subatomic world.

This development opens the door to more precise measurements of particles that were previously difficult to study. With the new tool, scientists can now look for subtle differences in how these particles decay, which could reveal new physics beyond the current understanding of the universe. The researchers have made their software available to the wider scientific community, allowing others to apply these techniques to their own experiments. The work demonstrates that even in a field as established as particle physics, there is still room for significant improvements in how we observe and understand the invisible forces that shape our reality.

Technical Summary: Forward Decay Chain Reconstruction with Displaced Vertices in Intense Magnetic Fields

Problem Statement
Precise reconstruction of particle decay chains is fundamental to flavor physics and spectroscopy, particularly for analyses involving complex topologies with multiple final-state particles and intermediate long-lived resonances. Standard vertex and kinematic fitting algorithms typically rely on local linear approximations of particle trajectories at the decay point. While effective for conventional experimental cases, these assumptions fail in scenarios involving hierarchical decay trees where particles traverse large detector volumes through regions of intense, inhomogeneous magnetic fields.

Specifically, the LHCb experiment faces challenges reconstructing $b$ - and $c$ -hadron decays containing long-lived particles (e.g., $\Lambda$ , $K^0_S$ ) that decay between 2.5 and 7.5 meters from the interaction point. In this region, the dipole magnetic field is maximal and non-uniform. Charged daughter particles often traverse several meters of uninstrumented volume before reaching the first measurement surface (the Scintillating Fibre tracker). Under these conditions, existing tools suffer from:

Numerical Instability: Repeated transport of track measurements and covariance matrices over long distances through non-linear magnetic fields becomes inaccurate.
Inadequate Modeling: Local linear approximations fail to describe trajectories where momentum resolution is limited and kinematic constraints (such as invariant mass) are tight.
Reduced Efficiency and Resolution: Current tools exhibit sharp drops in reconstruction efficiency and degraded invariant mass and vertex resolution for decays occurring deep within the magnetic field.

Methodology: ForwardTreeFitter (fTf)
The authors present ForwardTreeFitter, a standalone, modular C++ framework implementing a simultaneous Kalman filter technique to address these limitations. The core methodological innovation is a global track model that fundamentally changes how track residuals are evaluated:

Fixed Measurements, Propagated States: Unlike conventional fitters that transport track measurements to the estimated decay vertex, fTf keeps track measurements fixed at their reference surfaces. Instead, the fitted track states are iteratively propagated through the magnetic field map to these surfaces to evaluate measurement residuals.
Non-Linear Trajectory Modeling: The reference implementation utilizes a full numerical trajectory model based on Runge–Kutta–Nyström integration with adaptive step sizes. This allows for precise propagation through inhomogeneous magnetic fields. The framework also supports linear, parabolic, and cubic Hermite interpolators to balance performance and accuracy.
Constraint Handling: The framework treats exact constraints (vertex, momentum conservation, invariant mass) differently from linear approximations. Rather than linearizing both the trajectory and constraints simultaneously, fTf preserves the full non-linear trajectory model and re-linearizes only the exact constraints around the updated solution at each iteration. This improves numerical stability.
Seeding Strategies: To ensure convergence in highly non-linear topologies, the framework employs a multi-stage seeding process:
- FDoca: Calculates the Distance-of-Closest-Approach (DOCA) and Point-of-Closest-Approach (POCA) using Golden Section and Downhill Simplex algorithms combined with RK or Hermite extrapolators to determine initial vertex positions and momenta.
- FVertexer: A simplified vertex Kalman fit used to generate initial covariance matrices, which are then inflated to prevent the initial seed from dominating the global fit.

Key Contributions

Framework Design: A modular, extensible C++ framework (available under MPL-2.0) that separates the track measurement from the trajectory model, allowing for flexible application to both forward-geometry and $4\pi$ -coverage detectors.
Algorithmic Innovation: The implementation of a simultaneous fit that combines a global non-linear track model with iterative re-linearization of constraints, specifically optimized for intense, inhomogeneous magnetic fields.
Robust Initialization: A dedicated seeding suite (FDoca and FVertexer) designed to provide stable initial conditions for complex, displaced decay topologies where standard linear approximations fail.

Results
Benchmark studies were performed using simulated data for $\Lambda^0_b \to J/\psi \Lambda$ and $B^0 \to J/\psi K^0_S$ decays, reconstructed within the LHCb geometry. The performance of fTf was compared against the standard LHCb DecayTreeFitter (DTF) and ParticleVertexFitter (PVF).

Invariant Mass Resolution: For candidates reconstructed from "T tracks" (reconstructed exclusively in the SciFi tracker) in the 3–6 m decay region, fTf with mass constraints improved the $\Lambda^0_b$ mass resolution by a factor of 12.5 and the $B^0$ mass resolution by a factor of 5.9 compared to DTF.
Vertex Resolution: The $z$ -vertex resolution for $\Lambda$ candidates improved by a factor of 3.8 (from 570 mm to 150 mm) in the 3–6 m region.
Reconstruction Efficiency: While DTF efficiency dropped sharply to ~64% for decays between 3.5–7 m, fTf maintained a flat efficiency of approximately 96% across the full 3–8 m range when mass constraints were applied.
Seeding Performance: The FDoca algorithm improved the significance figure-of-merit ( $S = N_{sig}/\sqrt{N_{sig} + N_{bkg}}$ ) by 30% for $\Lambda$ and 50% for $K^0_S$ decays from TT pairs compared to the standard DOCA tool.
Limitations: The authors note that a population of events with large, non-Gaussian $z$ -vertex residuals persists, particularly for $\Lambda$ decays near the interaction point. This is attributed to intrinsic information limitations in nearly collinear, low- $Q$ two-body decays, rather than algorithmic failure.

Significance
The paper claims that ForwardTreeFitter provides a robust and numerically stable solution for reconstructing tightly constrained decay chains in regions of intense and inhomogeneous magnetic fields, a regime where standard linearized tools are insufficient. By evaluating track residuals at their measurement surfaces via non-linear propagation, the framework achieves significant gains in reconstruction efficiency and resolution.

The authors emphasize that these improvements open new opportunities for precision measurements and searches involving highly displaced decays. The framework is already being progressively incorporated into the LHCb reconstruction and analysis software stack. While developed for LHCb, the modular design allows it to accommodate general-purpose collider detector geometries and standard approximations, making it a versatile tool for particle physics experiments facing similar topological challenges.

More like this