A model independent method for measurement of $B^{\pm}$ and $B^0$ meson production fractions at $\Upsilon(4S)$

This paper proposes a model-independent method for directly measuring the production fractions of $B^{\pm}$ and $B^0$ mesons at the $\Upsilon(4S)$ resonance by counting single- and double-inclusive charmed mesons, demonstrating that world-average precision can be achieved without relying on underlying theoretical assumptions.

The Gist

The Big Picture: Counting Cookies in a Bakery

Imagine a high-end bakery (the B-factory) that only sells two types of cookie boxes:

1. The "Choco-Chip" Box (Charged B mesons, $B^+$ and $B^-$).
2. The "Oatmeal-Raisin" Box (Neutral B mesons, $B^0$ and $\bar{B}^0$).

The bakery owner (the $\Upsilon(4S)$ resonance) has a very specific rule: every time a box is made, it comes in a pair. You never get just one box; you always get two. Sometimes you get two Choco-Chip boxes, sometimes two Oatmeal-Raisin boxes, and sometimes one of each.

The Problem:
Scientists want to know the exact ratio of these boxes. Are 50% of the pairs Choco-Chip and 50% Oatmeal-Raisin? Or is it 51% and 49%? This ratio is crucial because it acts like a "conversion rate" for all other physics experiments. If you don't know the ratio, you can't accurately calculate how often specific things happen inside the boxes.

The Old Way:
Previously, scientists tried to guess the ratio by looking at specific, rare decorations on the cookies (like a specific type of chocolate chip). They had to assume that the decoration appeared equally often in both types of boxes. If that assumption was slightly wrong, their whole calculation was off. It was like trying to guess the ratio of two cookie types by assuming they both use the exact same brand of chocolate, which might not be true.

The New Method (The Paper's Solution):
This paper proposes a new, "model-independent" way to count the boxes. Instead of looking for specific decorations, they count every single crumb that falls out of the boxes.

The Analogy: The "Crumb Counting" Strategy

Imagine the bakery owner dumps the contents of two boxes onto a table. The contents are "crumbs" (particles like D-mesons and leptons).

• Oatmeal boxes tend to drop more Oatmeal crumbs ($D^0$).
• Choco-Chip boxes tend to drop more Chocolate crumbs ($D^*$).

The scientists' new method works like this:

1. The Single Count: They count how many Oatmeal crumbs and Chocolate crumbs are on the table in total.
2. The Double Count: They look for pairs of crumbs that fell out of the same box versus pairs that fell out of different boxes.
3. The Mix-and-Match: They look at specific combinations, like "One Oatmeal crumb + One Chocolate crumb."

Why is this clever?
In the old days, they tried to isolate just the "Oatmeal" boxes. But here, they accept that the table is a messy mix of both. They use math to solve a giant puzzle.

They set up a system of equations (a giant recipe) that says:

• "If I see this many Oatmeal crumbs and that many Chocolate crumbs, and I see this many pairs of them..."
• "...then the only way the math works out is if the bakery made X% Oatmeal boxes and Y% Choco-Chip boxes."

The Secret Sauce: The "Twin" Effect

There is a tricky part to the puzzle. Sometimes, an Oatmeal box can magically turn into a Choco-Chip box (or vice versa) before it's opened. In physics, this is called mixing ($B^0$-$\bar{B}^0$ mixing).

The paper explains that by paying attention to the charge (positive or negative) of the crumbs, they can detect this "magic trick."

• If they see two "positive" crumbs together, it usually means they came from two different boxes.
• But if they see a "positive" and a "negative" crumb together, it might mean one box turned into its twin before opening.

By tracking these charge combinations, the scientists can untangle the mess and solve the puzzle without needing to guess how the cookies were baked.

Why This Matters

1. No Guessing: The old methods relied on assumptions (like "these two decay rates are equal"). This new method relies purely on counting. It's like weighing the boxes instead of guessing their weight based on the label.
2. Precision: The authors ran a computer simulation (a "virtual bakery") to test their method. They found that with enough data, they can measure the ratio with the same high precision as the best current methods, but without the risk of their assumptions being wrong.
3. Finding the Hidden Boxes: The method is so good at counting that it can also detect if the bakery is secretly making a third type of box (non-B meson events) that nobody has seen before.

The Bottom Line

This paper presents a new, smarter way to count particle pairs. Instead of trying to separate the "good" cookies from the "bad" ones and making assumptions about them, the scientists count every single crumb and use a complex but robust mathematical recipe to figure out exactly how many of each cookie type were baked.

It's a shift from guessing based on theory to knowing based on pure data, making our understanding of the universe's building blocks much more reliable.

Technical Summary

1. Problem Statement

The production fractions of charged ($f_{+-}$) and neutral ($f_{00}$) $B$-meson pairs at the $\Upsilon(4S)$ resonance are critical parameters for converting observed event yields into absolute branching fractions. These values are essential for precision tests of the Cabibbo–Kobayashi–Maskawa (CKM) matrix, flavor symmetry tests, and searches for physics beyond the Standard Model.

Current Limitations:

• Model Dependence: Existing measurements (e.g., by CLEO, BaBar, Belle) rely on assumptions such as isospin symmetry (assuming equal partial decay widths for charged and neutral modes) or heavy-quark expansion arguments.
• Systematic Uncertainties: Methods using exclusive decays (e.g., $B \to J/\psi K$) are limited by the accuracy of isospin symmetry assumptions ($\sim 1\%$). Methods using semileptonic tags often require modeling of excited charm states ($D^{**}$) or rely on theoretical estimates of inclusive widths.
• Data Inconsistency: Current world averages combine results from experiments running at slightly different center-of-mass energies, introducing unquantified systematic effects related to beam-energy spreads.

The paper aims to provide a model-independent determination of $f_{+-}$ and $f_{00}$ that relies solely on event counting without assumptions about decay dynamics, isospin symmetry, or external branching fraction inputs.

2. Methodology

The proposed method utilizes a statistical tagging framework based on inclusive production signatures of charmed mesons ($D^0, D^{*\pm}$) and charged leptons ($\ell^\pm$) in $\Upsilon(4S) \to B\bar{B}$ events.

Core Concept: Double-Tagging with Charge Resolution

Unlike previous "double-tag" methods that attempted to isolate pure $B^0$ samples (which introduced background modeling uncertainties), this method treats $B^+$ and $B^0$ contributions simultaneously.

• Observables: The method counts single-tag (one reconstructed particle) and double-tag (two reconstructed particles) yields.
• Charge Separation: A crucial innovation is the resolution of particle charges (e.g., distinguishing $D^0$ from $\bar{D}^0$, and $\ell^+$ from $\ell^-$). This breaks the algebraic degeneracy found in charge-averaged systems.
• Kinematic Suppression: To ensure that double-tagged particles originate from different $B$ mesons (and not a single $B \to D\bar{D}X$ decay), a kinematic cut is applied: $p_D > 1.05$ GeV/c in the center-of-mass frame. This exploits the fact that $B$ mesons are nearly at rest at the $\Upsilon(4S)$; two energetic $D$ mesons cannot be emitted into the same hemisphere from a single $B$ decay.

The System of Equations

The method constructs a system of 27 equations relating observed yields ($N$) to unknown parameters:

1. Production Fractions: $f_{+-}$ and $f_{00}$ (with $f_{+-} + f_{00} \approx 1$).
2. Effective Inclusive Branching Fractions: 18 parameters representing products of physical branching fractions and reconstruction efficiencies for various modes (e.g., $B \to D^*X$, $B \to \ell X$, $B \to D\ell X$).

Key Components of the System:

• Single-Daughter Yields: Depend linearly on production fractions and inclusive rates.
• Double-Daughter Yields: Depend on the square of rates and include terms for $B^0-\bar{B}^0$ mixing ($\chi_d \approx 18.6\%$). The mixing term is essential as it introduces non-factorizable contributions that break degeneracy.
• Triple/Four-Daughter Yields: Include combinations involving leptons (e.g., $D\ell$, $\ell\ell$) to provide additional constraints.
• Discrete Ambiguity Resolution: The system initially has a discrete ambiguity where mixing effects can be mimicked by enhanced wrong-sign charm production. This is resolved by applying weak external constraints (20% uncertainty) on a subset of $D^*$ branching fractions, which are well-measured, to stabilize the fit without biasing the result.

3. Key Contributions

• Model Independence: The method does not assume isospin symmetry or equal semileptonic widths. All branching fractions and efficiencies are treated as free parameters determined directly from the data fit.
• Statistical Tagging: Instead of isolating specific exclusive channels, it uses the statistical differences in inclusive production rates of $D^0$, $D^*$, and $\ell$ between $B^+$ and $B^0$ decays.
• Charge-Resolved Formalism: By distinguishing right-sign and wrong-sign combinations, the method lifts the algebraic degeneracy inherent in previous inclusive counting approaches.
• Simultaneous Extraction: It allows for the simultaneous determination of $f_{+-}$, $f_{00}$, and the fraction of non-$B\bar{B}$ $\Upsilon(4S)$ decays ($f_{\text{non-}B\bar{B}}$).

4. Results (Feasibility Study)

A toy Monte-Carlo simulation was performed using $10^9$ events (approximating the full Belle dataset).

• Input Parameters: $f_{00} = 0.486$, $f_{+-} = 0.511$.
• Fitted Results:
  • $f_{00} = 0.4845 \pm 0.0079$
  • $f_{+-} = 0.5111 \pm 0.0069$
• Precision: The statistical precision achieved ($\sim 0.6\% - 0.8\%$) is comparable to the current world-average uncertainty.
• Non-$B\bar{B}$ Fraction: The method can also extract the fraction of $\Upsilon(4S)$ decays to non-$B\bar{B}$ final states with a precision of $\sigma(f_{\text{non-}B\bar{B}}) \approx 0.0022$.
• Systematics: The study indicates that systematic uncertainties (e.g., from hemisphere assumptions or misidentification) can be validated and corrected using data-driven checks (e.g., comparing same-sign $D$ pairs in same vs. opposite hemispheres).

5. Significance

This paper proposes a paradigm shift in measuring $B$-meson production fractions:

1. Robustness: By removing reliance on isospin symmetry and theoretical decay models, the measurement becomes a direct, data-driven observable.
2. Complementarity: It offers a cross-check to existing methods, potentially resolving discrepancies in world averages caused by energy-dependent systematic effects.
3. Future Applications: The technique is applicable to current and future $B$-factory datasets (Belle II), potentially improving the precision of CKM matrix element determinations and searches for new physics by reducing a major source of systematic uncertainty.

In conclusion, the authors demonstrate that a purely counting-based, charge-resolved analysis of inclusive charm and lepton production can achieve world-class precision on $B$-meson production fractions without the theoretical baggage of previous methods.