⚛️ high-energy experiments

Dalitz plot binning optimization for the model-independent measurement of the CKM angle $β$ with $B^0 \to D π^+ π^-$ , $D \to K_S^0 π^+ π^-$ decays

This paper presents an optimized binning scheme for the $B^0 \to D \pi^+ \pi^-$ Dalitz plot, derived using Fisher information and simulated annealing, which enhances the statistical sensitivity to the CKM angle $\beta$ by approximately 10% compared to a standard equal-phase-difference baseline while maintaining robustness against experimental effects.

Original authors: Anfeng Li, Wenbin Qian, Boan Shi, Ao Xu, Liming Zhang

Published 2026-10-14

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Anfeng Li, Wenbin Qian, Boan Shi, Ao Xu, Liming Zhang

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, matter and antimatter are not perfect mirror images. While the laws of physics suggest they should behave identically, a subtle difference in how they decay allows the universe to exist as we know it, filled with stars and galaxies rather than pure energy. This difference is called CP violation, and it is governed by a specific angle in the mathematical framework that describes how quarks mix and change. Physicists call this angle beta. Measuring it with extreme precision is a primary goal of modern particle physics because it acts as a strict test for our current understanding of the universe. If the measured value does not match the prediction, it could signal the presence of new, undiscovered particles or forces.

To find this angle, scientists study the decay of particles called B mesons. One particularly useful process involves a B meson transforming into a D meson and two pions, with the D meson subsequently breaking apart into a neutral kaon and two more pions. This chain of events creates a complex map of possibilities known as a Dalitz plot, which charts the energy and momentum of the particles produced. The challenge is that the information needed to calculate the angle beta is hidden within the interference patterns of these particles. To extract this information without relying on potentially flawed theoretical models, researchers divide this map into distinct regions, or bins, and count the events in each. The way these regions are drawn is critical; if the bins are drawn poorly, valuable information is lost, and the measurement becomes less precise.

A team of researchers has now developed a method to draw these bins in the most effective way possible, significantly improving the precision of the measurement. Instead of using a standard, uniform grid, they used a sophisticated computer algorithm to rearrange the boundaries of the bins, grouping together specific areas of the map that provide the most information about the angle beta. By testing millions of different arrangements, they found a configuration that captures the essential physics more efficiently than the traditional approach. Their work shows that this optimized layout can improve the statistical sensitivity of the experiment by approximately ten percent compared to the standard method. This gain is substantial in a field where every fraction of a degree matters, allowing future experiments to reach their goals with fewer particles or in less time.

The researchers validated their findings using detailed computer simulations that mimic the behavior of real particle collisions. They tested their new binning scheme against various real-world complications, such as the fact that detectors do not see every particle equally well and that background noise from other processes can obscure the signal. They found that their optimized scheme remains robust even when these experimental imperfections are taken into account. While the improvement is most pronounced in ideal conditions, the method holds up well when decay times and detection efficiencies are included, confirming that the gain is not just a mathematical artifact but a practical advantage for real data.

The study also explored how the number of bins affects the result. They discovered that increasing the number of bins beyond a certain point yields diminishing returns, as the data in each individual bin becomes too sparse to be useful. They identified a sweet spot with eight pairs of bins that balances the need for detailed information with the statistical power of the data. This specific configuration, which groups the most informative regions of the particle map together, is now recommended for use in upcoming experiments at major particle physics facilities. By adopting this optimized map, scientists will be able to measure the fundamental properties of the universe with greater clarity, bringing us closer to understanding why matter dominates over antimatter.

Technical Summary: Dalitz Plot Binning Optimization for Model-Independent $\beta$ Measurement

Problem Statement
The paper addresses the challenge of maximizing the statistical sensitivity in a model-independent measurement of the CKM angle $\beta$ using the decay chain $B^0 \to D\pi^+\pi^-$ , with $D \to K_S^0\pi^+\pi^-$ . While time-dependent analyses of this channel allow for the determination of $\beta$ without relying on a $B$ -decay amplitude model, the method requires dividing the Dalitz plots of both the $B$ and $D$ decays into bins. The choice of binning scheme for the $B^0 \to D\pi^+\pi^-$ Dalitz plot directly impacts the statistical sensitivity to $\beta$ . Previous approaches often utilized equal-phase-difference binning, but the authors posit that optimizing this binning scheme can significantly reduce the statistical uncertainty on $\beta$ .

Methodology
The authors propose an optimization strategy based on the Fisher information, $I$ , which bounds the statistical uncertainty of the parameter $\beta$ via the Cramér–Rao inequality ( $\sigma_\beta \ge I^{-1/2}$ ). To quantify the impact of binning, they define a quality factor $Q^2 = I_{\text{binned}} / I_{\text{unbinned}}$ , representing the fraction of Fisher information retained after binning. The goal is to maximize $Q^2$ .

Two methods for calculating $Q^2$ are evaluated:

Time-Integrated (TI) Method: An approximation analogous to $\gamma$ measurements, where time-dependent phase factors are treated as constant. This method yields a factorized quality function dependent only on the $B$ and $D$ binning separately.
Time-Dependent (TD) Method: This approach utilizes the full time-dependent probability density function (PDF), incorporating $B^0$ – $\bar{B}^0$ oscillations, decay time dependence, and flavor tagging (including mistag probability). This method does not factorize and provides a more complete description of the $\beta$ analysis.

The optimization is performed using a simulated annealing algorithm. Starting from an equal-phase-difference baseline, the algorithm iteratively reassigns pairs of pixels in the discretized Dalitz plot to different bin labels to maximize the $Q^2$ objective function. The $B$ -decay hadronic parameters ( $k_j, c_j, s_j$ ) required for the calculation are derived from a publicly available amplitude model of $B^0 \to D^0\pi^+\pi^-$ . The $D$ -decay binning is fixed to the established optimal scheme ( $N=8$ , 16 bins) derived from external strong-phase measurements.

Key Results

Sensitivity Improvement: For a nominal configuration with $M=8$ (16 bins) in the $B$ Dalitz plot, the TD-optimized binning improves the sensitivity to $\beta$ by approximately 10% relative to the equal-phase-difference baseline. The TI method predicts a larger improvement (~24% by its own metric), but this translates to only ~5% when evaluated with the more realistic TD metric.
Bin Characteristics: The optimization tends to group low-sensitivity regions of the Dalitz plot into specific bins (often near the origin in the $c_j$ - $s_j$ plane) to avoid diluting the sensitivity of high-sensitivity bins, which are pushed toward the physical boundary ( $c_j^2 + s_j^2 \approx 1$ ).
Robustness Checks:
- Number of Bins ( $M$ ): Sensitivity saturates around $M=8$ ; increasing $M$ further yields diminishing returns and reduces the statistical precision of hadronic parameter determination per bin.
- Experimental Effects: The improvement is robust against decay time acceptance, Dalitz plot efficiency variations, and finite Dalitz plot resolution. Including Dalitz plot efficiency in the evaluation increases the predicted improvement from 9.7% to 11.6%. Finite resolution increases absolute uncertainties but does not significantly alter the relative gain from optimization.
- Background: Including a representative background component ( $f_{\text{sig}} = 0.8$ ) during optimization changes the predicted improvement by about two percentage points but does not qualitatively alter the optimized parameters.
Validation: Toy Monte Carlo studies (250 pseudoexperiments) confirm that the $Q^2$ predictions translate to actual gains in sensitivity. The observed improvement in toy fits (e.g., 9.3% vs. 9.7% predicted) shows good agreement with the Fisher information-based predictions.

Significance and Claims
The paper claims that the proposed TD-optimized binning scheme is the recommended approach for future model-independent measurements of $\beta$ at LHCb and Belle II. The authors emphasize that:

The Time-Dependent method is superior to the Time-Integrated method for this specific physics goal, as it correctly accounts for the time-dependent oscillation structure essential to the $\beta$ measurement.
The optimization is robust against realistic experimental conditions, including efficiency variations and resolution effects, ensuring that the theoretical gain is not lost in practice.
The methodology provides a ~10% reduction in statistical uncertainty without requiring additional data or complex theoretical inputs, offering a straightforward path to improved precision in testing the Standard Model and constraining the Unitarity Triangle.

The authors conclude that the machine-readable definitions of the recommended $M=8$ symmetric binning will be made publicly available to facilitate these future measurements.

More like this