Model-agnostic likelihood for the reinterpretation of the $\boldsymbol{B^{+}\to K^{+} \nu \bar{\nu}}$ measurement at Belle II
| Abstract | We recently measured the branching fraction of the $B^{+}\rightarrow K^{+}\nu\bar{\nu}$ decay using 362~fb$^{-1}$ of on-resonance $e^+e^-$ collision data under the assumption of Standard Model kinematics, providing the first evidence for this decay. To facilitate future reinterpretations and maximize the scientific impact of this measurement, we publicly release the full analysis likelihood along with all necessary material required for reinterpretation under arbitrary theoretical models sensitive to this measurement. In this work, we demonstrate how the measurement can be reinterpreted within the framework of the Weak Effective Theory. Using a kinematic reweighting technique in combination with the published likelihood, we derive marginal posterior distributions for the Wilson coefficients, construct credible intervals, and assess the goodness of fit to the Belle~II data. For the Weak Effective Theory Wilson coefficients, the posterior mode of the magnitudes $|C_\mathrm{VL}+C_\mathrm{VR}|$, $|C_\mathrm{SL}+C_\mathrm{SR}|$, and $|C_\mathrm{TL}|$ corresponds to the point ${(11.3, 0.0, 8.2)}$. The respective 95\% credible intervals are $[1.9, 16.2]$, $[0.0, 15.4]$, and $[0.0, 11.2]$. |
|---|---|
| Tags | |
| Working group | Radiative and Electroweak Penguin |
| Principal authors | @ Lorenz Gaertner, Nikolai Krug, Thomas Kuhr, Slavomira Stefkova |
| Reference | published in PRD |
| Document | BELLE2-PUB-PH-2024-038 |
| Links | arXiv:2507.12393, |
The ITA binned joint number densities. The horizontal axis corresponds to the generated $q^2$. The vertical axis represents the binning used in the ${B^{+}\to K^{+}\nu\bar{\nu}}$ analysis~\cite{Belle-II:2023esi}. The heatmap shows the weighted signal events.
The ${B^{+}\to K^{+}\nu\bar{\nu}}$ differential branching fraction prediction from Eq. (9). Individual contributions are shown here with the combinations of vector, scalar and tensor Wilson coefficients set to unity, respectively. The uncertainties shown (bands) stem from the hadronic parameters.
The marginalized posterior for the Wilson coefficients in Eq. (9). We adopt the convention that $C_{\mathrm{VL}}+C_{\mathrm{VR}}$, $C_{\mathrm{SL}}+C_{\mathrm{SR}}$ and $C_{\mathrm{TL}}$ are real valued. Diagonal and off-diagonal panels show the 1-dimensional and 2-dimensional sample density PDFs on a linear scale, respectively. The overall scale is omitted, as all relevant information is contained in the shape of the distribution. The contours indicate $68\%$ and $95\%$ credible intervals. The dashed black lines and cross mark the SM point; the dash-dotted yellow lines and cross indicate the posterior mode; dotted red lines mark the symmetry axes used for sample symmetrization.
Observed and predicted best-fit yields in the highest-sensitivity bins of the analysis. These correspond to the $\eta(\mathrm{BDT}_2)>0.98$ region of the ITA. The signal is shown for the unconstrained ${B^{+}\to K^{+}\nu\bar{\nu}}$ SM (\textit{left}) and the WET (\textit{right}) predictions. The predicted background yields are shown individually for the neutral and charged \B-meson decays, and the summed five continuum categories. Pulls are shown in the lower panels.