Measurement of the branching fractions of $\chi_{cJ} \to \pi^{+}\pi^{-}\pi^{0}\pi^{0}$ via $\psi(3686) \to \gamma\chi_{cJ}$

Using a large sample of $\psi(3686)$ events collected by the BESIII detector, this study precisely measures the branching fractions of $\chi_{cJ} \to \pi^{+}\pi^{-}\pi^{0}\pi^{0}$ decays for $J=0, 1, 2$, identifying $\rho^+\rho^-$ as the dominant intermediate state and providing results that significantly improve upon previous measurements.

The Gist

The Big Picture: Catching Ghosts in a Cosmic Factory

Imagine the universe is a giant, high-speed factory. In this factory, scientists at the BESIII lab in Beijing are trying to study a very specific, very rare type of "ghost" particle called the $\chi_c$ (Chi-c).

These ghosts are made of a heavy charm quark and its anti-particle, stuck together like a tiny, vibrating dumbbell. They are unstable and fall apart almost instantly into other particles. The problem? They are incredibly hard to catch directly.

The Trick:
Instead of trying to catch the ghost directly, the scientists use a clever sleight of hand. They look at a slightly heavier, more stable cousin called the $\psi(3686)$. Think of the $\psi(3686)$ as a glowing, energetic balloon. When this balloon pops, it doesn't just disappear; it releases a flash of light (a photon, $\gamma$) and leaves behind the ghost we are looking for ($\chi_c$).

The paper reports on a massive experiment where the team collected 2.7 billion of these "balloon pops" ($\psi(3686)$ events). That is a lot of data—like watching a fireworks display for several years straight without blinking.

The Mystery: What Happens When the Ghost Dies?

Once the $\chi_c$ ghost is created, it immediately decays (dies). The scientists wanted to know: What does it turn into?

Specifically, they were looking for a very specific "four-piece puzzle" outcome:

• Two charged pions (like tiny, electrically charged marbles: $\pi^+$ and $\pi^-$)
• Two neutral pions (invisible marbles that instantly turn into light: $\pi^0$ and $\pi^0$)

The question was: How often does this specific four-piece puzzle happen? In physics, this frequency is called the branching fraction.

The Investigation: Sorting the Noise

The detector (BESIII) is like a giant, 3D camera that takes a picture of every single particle flying out of the collision. But the photos are messy.

• The Signal: The four-pion puzzle we want.
• The Noise: Millions of other things happening at the same time (other particles, random glitches, or different decay patterns).

How they cleaned up the mess:

1. The Kinematic Fit (The "Balance Scale"): The scientists used a mathematical trick called a "kinematic fit." Imagine you know exactly how much energy went into the collision. You then check if the energy of the four pions you found adds up perfectly to that starting amount. If it balances, it's likely a real event. If it's off, it's noise.
2. The Veto (The "Bouncer"): They set up strict rules to kick out imposters. For example, if the particles looked like they came from a different decay (like a $J/\psi$ particle), they were thrown out of the club.
3. The Histogram (The "Mountain Range"): After filtering, they plotted the mass of the particle groups. Real signals show up as sharp peaks (like mountains), while background noise looks like a flat, rolling hill. They found three distinct mountains corresponding to the three types of $\chi_c$ ghosts ($\chi_{c0}$, $\chi_{c1}$, and $\chi_{c2}$).

The Discovery: The "Double-Rho" Dance

Once they counted the peaks, they calculated the branching fractions (the probability of this happening). Here are their results, simplified:

• $\chi_{c0}$: About 3.1% of the time, it turns into the four-pion puzzle.
• $\chi_{c1}$: About 1.2% of the time.
• $\chi_{c2}$: About 1.9% of the time.

The "Aha!" Moment:
The scientists didn't just count the particles; they looked at how they were arranged. They discovered that these four pions don't just appear randomly. They are usually formed via an intermediate step involving $\rho$ mesons.

Think of it like this: The $\chi_c$ ghost doesn't just break into four pieces. It first splits into two pairs of dancing partners ($\rho^+$ and $\rho^-$), and then those pairs break apart into the four pions. The paper confirms that this "Double-Rho" dance is the dominant way this decay happens.

Why Does This Matter?

1. Precision: Previous measurements were like looking at a blurry photo. This new study is like switching to a 4K camera. The precision is nearly 10 times better than before.
2. Understanding the Rules: The strong force (the glue holding quarks together) is notoriously difficult to calculate. By measuring exactly how often these decays happen, scientists can test their theories about how quarks interact.
3. Superseding the Past: This paper officially replaces the best previous measurements (made by the CLEO collaboration years ago) with much more accurate numbers.

The Bottom Line

The BESIII team acted like master detectives. They sifted through 2.7 billion cosmic events, filtered out the noise, and found that the $\chi_c$ particles have a very specific habit: they love to decay into four pions, usually by first splitting into two "rho" mesons.

This isn't just about counting particles; it's about understanding the fundamental "dance steps" of the universe's building blocks. With these new, highly precise numbers, physicists can now write better rules for the game of quantum mechanics.

Technical Summary

1. Problem Statement

Charmonium decays serve as a crucial laboratory for studying quark-gluon dynamics at low energies. While the $P$-wave triplet states $\chi_{cJ}$ ($J=0, 1, 2$) are well-known, many of their hadronic decay modes remain unobserved, or existing measurements suffer from limited statistical precision.

• Specific Gap: The branching fractions for the decay $\chi_{cJ} \to \pi^+\pi^-\pi^0\pi^0$ had not been measured with high precision. Previous measurements (e.g., by CLEO) were limited by statistics.
• Theoretical Context: Understanding these decays is essential for testing theoretical models, such as the color-octet mechanism in $P$-wave charmonium decays.
• Objective: To measure the branching fractions of $\chi_{cJ} \to \pi^+\pi^-\pi^0\pi^0$ with significantly improved precision and to investigate the intermediate resonant states contributing to these decays.

2. Methodology

Data Sample and Detector

• Data Source: The analysis utilizes a dataset of $(2712.4 \pm 14.3) \times 10^6$ $\psi(3686)$ events collected by the BESIII detector at the BEPCII collider. The data corresponds to an integrated luminosity of approximately $3.9$ fb$^{-1}$ collected across three running periods (2009, 2012, 2021).
• Detector: The BESIII detector features a magnetic spectrometer with a 1.0 T superconducting solenoid, a Multi-Drift Chamber (MDC), a Time-of-Flight (TOF) system, and a CsI(Tl) Electromagnetic Calorimeter (EMC).

Event Selection

• Final State: The signal process is $\psi(3686) \to \gamma \chi_{cJ} \to \gamma \pi^+\pi^-\pi^0\pi^0$, where each $\pi^0$ decays to $\gamma\gamma$. The final state requires two charged tracks (identified as pions) and at least five photons (one from the initial radiative transition and four from the two $\pi^0$s).
• Kinematic Fit: A six-constraint (6C) kinematic fit is applied, constraining the total energy-momentum to the initial beam energy and the masses of the two $\pi^0$ candidates to the known $\pi^0$ mass.
• Optimization: The cut on the $\chi^2$ of the 6C fit ($\chi^2_{6C} < 60$) was optimized using a Figure-of-Merit (FOM) to maximize signal significance while minimizing background.
• Background Veto:
  • $\pi\pi J/\psi$ backgrounds: Vetoed by requiring the recoil mass of $\pi\pi$ pairs to be outside the $J/\psi$ mass window.
  • $\eta J/\psi$ backgrounds: Vetoed by checking the invariant mass of $\pi^+\pi^-\pi^0$ combinations.
  • $K_S^0 K_S^0$ peaking background: Vetoed using a 2D mass window on $M(\pi^+\pi^-)$ and $M(\pi^0\pi^0)$.

Signal Extraction and Branching Fraction Calculation

• Fitting: Signal yields ($N_{obs}$) are extracted via an unbinned maximum likelihood fit to the invariant mass distribution $M(\pi^+\pi^-\pi^0\pi^0)$.
  • Signal PDF: Modeled using Monte Carlo (MC) shapes weighted by an energy-dependent function $w$ to account for the $E1$ transition vertex dynamics and infrared divergence, convolved with a Gaussian resolution function.
  • Background PDFs: Includes continuum background (from off-resonance data), peaking backgrounds (fixed from MC), and combinatorial backgrounds (modeled by a Gaussian for $\pi^0 h_1$ and an Argus function).
• Efficiency Correction: The detection efficiency is corrected using a multidimensional reweighting method (Eq. 5.1) to account for discrepancies between data and MC in the angular ($\cos\theta$) and momentum distributions of the pions, driven by the presence of intermediate resonances (specifically $\rho$ mesons).

Systematic Uncertainties

Systematic uncertainties are evaluated for:

• Tracking and Particle Identification (PID) efficiencies.
• $\pi^0$ reconstruction and photon detection.
• MC modeling of intermediate states (reweighting uncertainties).
• Kinematic fit performance.
• Signal yield extraction (fit range, background modeling, damping function).
• Total number of $\psi(3686)$ events and PDG branching fractions.

3. Key Contributions

1. High-Precision Measurements: The paper provides the most precise measurements to date for the branching fractions of $\chi_{c0,1,2} \to \pi^+\pi^-\pi^0\pi^0$, improving statistical precision by nearly an order of magnitude compared to previous results (e.g., CLEO).
2. Intermediate State Analysis: The analysis identifies the dominant intermediate states in these decays. The scatter plots of $M(\pi^+\pi^0)$ vs. $M(\pi^-\pi^0)$ clearly show the presence of $\rho^+\rho^-$, $\rho\pi\pi$, and $\rho(1700)\pi\pi$ contributions, with $\rho^+\rho^-$ being the dominant channel.
3. Methodological Refinement: The application of a multidimensional reweighting technique to correct MC efficiencies based on data distributions ensures that the systematic uncertainty arising from imperfect modeling of the complex multi-pion dynamics is minimized.
4. Neutral-to-Charged Ratios: The paper calculates the ratios $R_J = \mathcal{B}(\chi_{cJ} \to \pi^+\pi^-\pi^0\pi^0) / \mathcal{B}(\chi_{cJ} \to 2\pi^+2\pi^-)$, providing new constraints on the isospin structure of these decays.

4. Results

The measured branching fractions are:

• $\chi_{c0} \to \pi^+\pi^-\pi^0\pi^0$: $(3.10 \pm 0.01 \text{ (stat)} \pm 0.14 \text{ (syst)}) \times 10^{-2}$
• $\chi_{c1} \to \pi^+\pi^-\pi^0\pi^0$: $(1.16 \pm 0.01 \text{ (stat)} \pm 0.05 \text{ (syst)}) \times 10^{-2}$
• $\chi_{c2} \to \pi^+\pi^-\pi^0\pi^0$: $(1.92 \pm 0.01 \text{ (stat)} \pm 0.08 \text{ (syst)}) \times 10^{-2}$

Key Observations:

• The results supersede previous BESIII measurements and are in good agreement with CLEO data but with significantly reduced uncertainties.
• The dominant intermediate state is found to be $\rho^+\rho^-$.
• The neutral-to-charged ratios are calculated as:
  • $R_0 = 1.57 \pm 0.11$
  • $R_1 = 1.70 \pm 0.11$
  • $R_2 = 1.70 \pm 0.12$

5. Significance

• Theoretical Impact: These precise measurements provide critical input for QCD calculations involving $P$-wave charmonium decays. The dominance of the $\rho^+\rho^-$ intermediate state offers insights into the hadronization mechanisms and the role of the color-octet mechanism.
• Experimental Benchmark: The results set a new standard for precision in charmonium spectroscopy. The improved precision helps resolve discrepancies in previous datasets and constrains theoretical models of multi-pion final states, which involve complex interplay between QCD, QED, and angular momentum couplings.
• Future Physics: The methodology developed for handling complex intermediate states and correcting MC efficiencies via data-driven reweighting can be applied to other multi-body hadronic decay analyses at BESIII and future experiments.