Isospin analysis of charmless B-meson decays

TL;DR

This paper analyzes charmless B-meson decays to precisely determine the CKM angle alpha, using latest data, isospin methods, and assessing uncertainties, with implications for understanding QCD dynamics and future measurements.

Contribution

It provides a comprehensive extraction of the CKM angle alpha from multiple decay channels and evaluates uncertainties and theoretical assumptions involved.

Findings

Determined alpha with combined data: (86.2^{+4.4}_{-4.0} or 178.4^{+3.9}_{-5.1}) degrees.

Assessed uncertainties from isospin symmetry breaking and statistical choices.

Compared hadronic amplitude results with theoretical expectations.

Abstract

We discuss the determination of the CKM angle $\alpha$ using the non-leptonic two-body decays $B\to \pi\pi$, $B\to \rho\rho$ and $B\to \rho\pi$ using the latest data available. We illustrate the methods used in each case and extract the corresponding value of $\alpha$. Combining all these elements, we obtain the following determination of $\alpha_{\rm dir}=( 86.2^{+4.4}_{-4.0}~~ \cup~~ 178.4^{+3.9}_{-5.1})^\circ$. We assess the uncertainties associated to the breakdown of the isospin hypothesis and the choice of statistical framework in detail. We also determine the hadronic amplitudes (tree and penguin) describing the QCD dynamics involved in these decays, briefly comparing our results with theoretical expectations. For each observable of interest in the $B\to \pi\pi$, $B\to \rho\rho$ and $B\to \rho\pi$ systems, we perform an indirect determination based on the constraints from all the other observables available and we discuss the compatibility between indirect and direct determinations. Finally, we review the impact of future improved measurements on the determination of $\alpha$.

Tables

Table 1. Table 1: General decomposition of the amplitudes A i ​ j = ⟨ h 1 i ​ h 2 j | ℋ s | u ​ u ¯ ​ d ¯ , q ⟩ superscript 𝐴 𝑖 𝑗 quantum-operator-product superscript subscript ℎ 1 𝑖 superscript subscript ℎ 2 𝑗 subscript ℋ 𝑠 𝑢 ¯ 𝑢 ¯ 𝑑 𝑞 A^{ij}=\langle h_{1}^{i}h_{2}^{j}|{\cal H}_{s}|u\bar{u}\bar{d},q\rangle ( q = u , d 𝑞 𝑢 𝑑 q=u,d ; i , j = − , 0 , + formulae-sequence 𝑖 𝑗 0 i,j=-,0,+ ) in terms of the reduced matrix elements A Δ ​ I , I f subscript 𝐴 Δ 𝐼 subscript 𝐼 𝑓 A_{\Delta I,I_{f}} for a pair of distinguishable isovector mesons h 1 i superscript subscript ℎ 1 𝑖 h_{1}^{i} and h 2 j superscript subscript ℎ 2 𝑗 h_{2}^{j} .

$A^{ij}=⟨h_1^i{h}_2^j\left|{\mathscr{H}}_s\right|u\overline{u}\overline{d},q⟩$	$A_{\frac{5}{2},2}$	$A_{\frac{3}{2},2}$	$A_{\frac{3}{2},1}$	$A_{\frac{1}{2},1}$	$A_{\frac{1}{2},0}$
$A^{+0}=⟨h_1^{+}{h}_2^0\left|{\mathscr{H}}_s\right|u\overline{u}\overline{d},u⟩$	$-\sqrt{1/6}$	$+\sqrt{3/8}$	$-\sqrt{1/8}$	$+\sqrt{1/2}$	0
$A^{0+}=⟨h_1^0{h}_2^{+}\left|{\mathscr{H}}_s\right|u\overline{u}\overline{d},u⟩$	$-\sqrt{1/6}$	$+\sqrt{3/8}$	$+\sqrt{1/8}$	$-\sqrt{1/2}$	0
$A^{+-}=⟨h_1^{+}{h}_2^{-}\left|{\mathscr{H}}_s\right|u\overline{u}\overline{d},d⟩$	$+\sqrt{1/12}$	$+\sqrt{1/12}$	$+1/2$	$+1/2$	$+\sqrt{1/6}$
$A^{-+}=⟨h_1^{-}{h}_2^{+}\left|{\mathscr{H}}_s\right|u\overline{u}\overline{d},d⟩$	$+\sqrt{1/12}$	$+\sqrt{1/12}$	$-1/2$	$-1/2$	$+\sqrt{1/6}$
$A^{00}=⟨h_1^0{h}_2^0\left|{\mathscr{H}}_s\right|u\overline{u}\overline{d},d⟩$	$+\sqrt{1/3}$	$+\sqrt{1/3}$	0	0	$-\sqrt{1/6}$

Table 2. Table 2: Decomposition of the amplitude ⟨ h i ​ h j | ℋ s | u ​ u ¯ ​ d ¯ , q ⟩ quantum-operator-product superscript ℎ 𝑖 superscript ℎ 𝑗 subscript ℋ 𝑠 𝑢 ¯ 𝑢 ¯ 𝑑 𝑞 \langle h^{i}h^{j}|{\cal H}_{s}|u\bar{u}\bar{d},q\rangle ( q = u , d 𝑞 𝑢 𝑑 q=u,d ; i , j = − , 0 , + formulae-sequence 𝑖 𝑗 0 i,j=-,0,+ ) in terms of the isospin amplitudes A Δ ​ I , I f subscript 𝐴 Δ 𝐼 subscript 𝐼 𝑓 A_{\Delta I,I_{f}} for indistinguishable mesons in the final state ( h 1 = h 2 = h subscript ℎ 1 subscript ℎ 2 ℎ h_{1}=h_{2}=h ). A global factor 3 3 \sqrt{3} is applied to all the coefficients with respect to the general coefficients given in Tab.  1 .

$A^{ij}=⟨h^i{h}^j\left|{\mathscr{H}}_s\right|u\overline{u}\overline{d},q⟩$	$A_{\frac{5}{2},2}$	$A_{\frac{3}{2},2}$	$A_{\frac{1}{2},0}$
$A^{+0}=⟨{\pi }^{+}{\pi }^0\left|{\mathscr{H}}_s\right|u\overline{u}\overline{d},u⟩$	$-1$	$3/2$	$0$
$A^{+-}=⟨{\pi }^{+}{\pi }^{-}\left|{\mathscr{H}}_s\right|u\overline{u}\overline{d},d⟩$	$\sqrt{1/2}$	$+\sqrt{1/2}$	$+1$
$A^{00}=⟨{\pi }^0{\pi }^0\left|{\mathscr{H}}_s\right|u\overline{u}\overline{d},d⟩$	$+1$	$+1$	$-\sqrt{1/2}$

Table 3. Table 3: World averages for the relevant experimental observables in the B → π i ​ π j → 𝐵 superscript 𝜋 𝑖 superscript 𝜋 𝑗 B\rightarrow\pi^{i}\pi^{j} modes: branching fraction ℬ π ​ π i ​ j subscript superscript ℬ 𝑖 𝑗 𝜋 𝜋 {\cal B}^{ij}_{\pi\pi} , time-integrated C ​ P 𝐶 𝑃 CP asymmetry 𝒞 π ​ π i ​ j subscript superscript 𝒞 𝑖 𝑗 𝜋 𝜋 {\cal C}^{ij}_{\pi\pi} , time-dependent asymmetry 𝒮 π ​ π i ​ j subscript superscript 𝒮 𝑖 𝑗 𝜋 𝜋 {\cal S}^{ij}_{\pi\pi} and correlation ( ρ 𝜌 \rho ).

Observable	World average	References
${ℬ}_{\pi \pi }^{+-}$ $(\times {10}^6)$	$5.10\pm 0.19$	[24, 25, 26, 27, 28]
${ℬ}_{\pi \pi }^{+0}$ $(\times {10}^6)$	$5.48\pm 0.34$	[25, 26, 27, 29]
${ℬ}_{\pi \pi }^{00}$ $(\times {10}^6)$	$1.59\pm 0.18$	[24, 30]
${𝒞}_{\pi \pi }^{00}$	$-0.34\pm 0.22$	[24, 30]
${𝒞}_{\pi \pi }^{+-}$	$-0.284\pm 0.039$	[31, 32, 33]
${𝒮}_{\pi \pi }^{+-}$	$-0.672\pm 0.043$	[31, 32, 33]
$\rho (C_{\pi \pi }^{+-},S_{\pi \pi }^{+-})$	$+0.013$	[31, 32, 33]

Table 4. Table 4: World averages for the relevant experimental observables in the B → ρ i ​ ρ j → 𝐵 superscript 𝜌 𝑖 superscript 𝜌 𝑗 B\rightarrow\rho^{i}\rho^{j} modes: branching fraction ℬ ρ ​ ρ i ​ j subscript superscript ℬ 𝑖 𝑗 𝜌 𝜌 {\cal B}^{ij}_{\rho\rho} , fraction of longitudinal polarisation f L i ​ j superscript subscript 𝑓 𝐿 𝑖 𝑗 f_{L}^{ij} , time-integrated C ​ P 𝐶 𝑃 CP asymmetry 𝒞 ρ ​ ρ i ​ j subscript superscript 𝒞 𝑖 𝑗 𝜌 𝜌 {\cal C}^{ij}_{\rho\rho} , time-dependent asymmetry 𝒮 ρ ​ ρ i ​ j subscript superscript 𝒮 𝑖 𝑗 𝜌 𝜌 {\cal S}^{ij}_{\rho\rho} and correlation ( ρ 𝜌 \rho ).

Observable	World average	References
${ℬ}_{\rho \rho }^{+-}\times f_L^{+-}$ $(\times {10}^6)$	$\left(27.76\pm 1.84\right)$ $\times$ $\left(0.990\pm 0.020\right)$	[43, 44]
${ℬ}_{\rho \rho }^{+0}$ $\times$ $f_L^{+0}$ $(\times {10}^6)$	$\left(24.9\pm 1.9\right)$ $\times$ $\left(0.950\pm 0.016\right)$	[45, 46]
${ℬ}_{\rho \rho }^{00}$ $\times$ $f_L^{00}$ $(\times {10}^6)$	$\left(0.93\pm 0.14\right)$ $\times$ $\left(0.71\pm 0.06\right)$	[47, 48, 49]
${𝒞}_{{\rho }_L{\rho }_L}^{+-}$	$-0.00\pm 0.09$	[43, 44]
${𝒮}_{{\rho }_L{\rho }_L}^{+-}$	$-0.15\pm 0.13$	[43, 44]
$\rho (C_{{\rho }_L{\rho }_L}^{+-},S_{{\rho }_L{\rho }_L}^{+-})$	$+0.0002$	[43, 44]
${𝒞}_{{\rho }_L{\rho }_L}^{00}$	$0.2\pm 0.9$	[47, 48, 49]
${𝒮}_{{\rho }_L{\rho }_L}^{00}$	$0.3\pm 0.7$	[47, 48, 49]

Table 5. Table 8: World averages for the branching fractions ℬ ρ ​ π i ​ j subscript superscript ℬ 𝑖 𝑗 𝜌 𝜋 {\cal B}^{ij}_{\rho\pi} and time-integrated C ​ P 𝐶 𝑃 CP asymmetries 𝒞 ρ ​ π i ​ j subscript superscript 𝒞 𝑖 𝑗 𝜌 𝜋 {\cal C}^{ij}_{\rho\pi} for the charged (top) and neutral (bottom) B → ρ ​ π → 𝐵 𝜌 𝜋 B\rightarrow\rho\pi decay modes.

Observable	World average	References
${ℬ}_{\rho \pi }^{+0}$ $(\times {10}^6)$	$10.9\pm 1.4$	[54, 56, 57]
${𝒞}_{\rho \pi }^{+0}$	$-0.02\pm 0.11$	[54, 56, 57]
${ℬ}_{\rho \pi }^{0+}$ $(\times {10}^6)$	$8.3\pm 1.2$	[54, 58, 59]
${𝒞}_{\rho \pi }^{0+}$	$-{0.18}_{-0.09}^{+0.17}$	[54, 58, 59]
${ℬ}_{\rho \pi }^{\pm \mp }$ $(\times {10}^6)$	$23.0\pm 2.3$	[54, 52, 55]
${ℬ}_{\rho \pi }^{00}$ $(\times {10}^6)$	$2.0\pm 0.5$	[54, 53, 55]

Table 6. Table 9: World averages for the branching ratios and direct C ​ P 𝐶 𝑃 CP asymmetries of the B 0 , + → π 0 , + ​ η ( ′ ) B^{0,+}\rightarrow\pi^{0,+}\eta^{(^{\prime})} modes  [ 68 ] .

Observable	World average
${ℬ}^{{\pi }^0\eta }$ $(\times {10}^6)$	${0.41}_{-0.17}^{+0.18}$
${ℬ}^{{\pi }^0{\eta }^{\prime }}$ $(\times {10}^6)$	$1.2\pm 0.4$
${ℬ}^{{\pi }^{+}\eta }$ $(\times {10}^6)$	$4.02\pm 0.27$
${ℬ}^{{\pi }^{+}{\eta }^{\prime }}$ $(\times {10}^6)$	${2.7}_{-0.4}^{+0.5}$
${𝒞}^{{\pi }^{+}\eta }$	$0.14\pm 0.05$
${𝒞}^{{\pi }^{+}{\eta }^{\prime }}$	$-0.06\pm 0.15$

Table 7. Table 10: World averages for the C + 0 superscript 𝐶 0 C^{+0} direct C ​ P 𝐶 𝑃 CP asymmetries for the B → π + ​ π 0 → 𝐵 superscript 𝜋 superscript 𝜋 0 B\rightarrow\pi^{+}\pi^{0} decay and for the longitudinally polarised state in B + → ρ + ​ ρ 0 → superscript 𝐵 superscript 𝜌 superscript 𝜌 0 B^{+}\rightarrow\rho^{+}\rho^{0} decay.

Observable	World average	Reference
${𝒞}_{\pi \pi }^{+0}$	$\left(-2.6\pm 3.9\right){10}^{-2}$	[29, 25]
${𝒞}_{{\rho }_L{\rho }_L}^{+0}$	$\left(-5.1\pm 5.4\right){10}^{-2}$	[45, 35]

Table 8. Table 11: 68% CL intervals for the weak phase α 𝛼 \alpha for different theoretical hypotheses, statistical approaches and channels. In the “Iso” column, Y indicates that the analysis is performed assuming isospin symmetry, wheras N denotes the inclusion of isospin breaking, namely, the effect of π 0 − η − η ′ superscript 𝜋 0 𝜂 superscript 𝜂 ′ \pi^{0}-\eta-\eta^{\prime} mixing in B → π ​ π → 𝐵 𝜋 𝜋 B\rightarrow\pi\pi and the breakdown of isospin triangle relations up to | r T | , | r P | < 4 % subscript 𝑟 𝑇 subscript 𝑟 𝑃 percent 4 |r_{T}|,|r_{P}|<4\% in B → ρ ​ ρ → 𝐵 𝜌 𝜌 B\rightarrow\rho\rho , as discussed in Sec.  4.1.2 (isospin-breaking effects are neglected in the B 0 → ( ρ ​ π ) 0 → superscript 𝐵 0 superscript 𝜌 𝜋 0 B^{0}\rightarrow(\rho\pi)^{0} Dalitz analysis). In the “ ewp ” column, Y indicates that a contamination from Δ ​ I = 3 / 2 Δ 𝐼 3 2 \Delta I=3/2 electroweak penguins is included as discussed in Sec.  4.1.1 , whereas N corresponds to setting these penguin contributions to zero. In the “Stat” column, B corresponds to the bootstrap approach and W to the Wilks-based one. For each channel or combination of channels, the deviation with respect to the indirect determination α ind subscript 𝛼 ind \alpha_{\rm ind} is indicated within brackets.

​​Iso​​	​​ewp​​	​​ Stat​​	$B\to \rho \rho$	$B\to \pi \pi$	$B\to \pi \pi$ + $\rho \rho$	$B^0\to {\left(\rho \pi \right)}^0$	All combined
Y	N	W ​	${92.0}_{-4.8}^{+4.7}$ (0.03$\sigma$) ​	$93.0\pm 14.0$ (0.6$\sigma$) ​	${92.1}_{-5.5}^{+5.2}$ (0.03$\sigma$) ​	${54.1}_{-10.3}^{+7.7}$ (3.0$\sigma$) ​	${86.2}_{-4.0}^{+4.4}$ (1.3$\sigma$)
Y	N	B ​	${92.0}_{-4.0}^{+4.2}$ (0.03$\sigma$) ​	$92.5\pm 13.5$ (1.1$\sigma$) ​	${92.2}_{-5.3}^{+4.9}$ (0.03$\sigma$) ​	${54.0}_{-17.0}^{+10.0}$ (2.7$\sigma$) ​	${86.1}_{-5.0}^{+5.3}$ (1.1$\sigma$)
Y	Y	W ​	${90.1}_{-4.8}^{+4.7}$ (0.4$\sigma$) ​	$91.2\pm 14.2$ (0.5$\sigma$)​	${90.1}_{-5.6}^{+5.1}$ (0.3$\sigma$) ​	${52.9}_{-11.1}^{+8.7}$ (3.0$\sigma$) ​	${85.6}_{-4.2}^{+4.1}$ (1.5$\sigma$)
N	N	W ​	${91.4}_{-5.7}^{+4.9}$ (0.2$\sigma$) ​	$92.5\pm 15.5$ (0.4$\sigma$) ​	${91.4}_{-6.0}^{+5.4}$ (0.2$\sigma$) ​	${54.1}_{-10.3}^{+7.7}$ (3.0$\sigma$) ​	${84.4}_{-4.3}^{+5.2}$ (1.5$\sigma$)
N	Y	W ​	${89.8}_{-4.9}^{+4.6}$ (0.5$\sigma$) ​	$91.0\pm 15.0$ (0.3$\sigma$) ​	${89.8}_{-5.3}^{+4.9}$ (0.5$\sigma$) ​	${52.9}_{-11.1}^{+8.7}$ (3.0$\sigma$) ​	${83.3}_{-3.1}^{+6.1}$ (1.6$\sigma$)

Table 9. Table 12: 68% CL interval for α 𝛼 \alpha from the S ​ U ​ ( 2 ) 𝑆 𝑈 2 {SU}\mkern-2.0mu(2) isospin analysis of the B → ρ ​ ρ → 𝐵 𝜌 𝜌 B\rightarrow\rho\rho data in the case that the uncertainty of some observables is reduced by a factor 2 2 \sqrt{2} (third column) or set to zero (fourth column). The relative gain in resolution with respect to the current measurement is indicated within brackets. Only the preferred solution for α 𝛼 \alpha close to 90 ∘ superscript 90 90^{\circ} is reported.

Decay	Analysis	Improved data	${\sigma }_{𝒪}\to {\sigma }_{𝒪}^{WA}/\sqrt{2}$	${\sigma }_{𝒪}\to 0$
$B^0\to {\rho }^{+}{\rho }^{-}$	C + TD	${ℬ}^{+-},f_L^{+-},{𝒞}^{+-},{𝒮}^{+-}$	${91.8}_{-4.0}^{+4.4}$ (-12%)	$92.2\pm 3.5$ (-26%)
	C	${ℬ}^{+-},f_L^{+-}$	${91.8}_{-4.5}^{+4.9}$ (-1%)	${92.1}_{-4.7}^{+4.6}$ (-2%)
	TD	${𝒞}^{+-},S^{+-}$	${92.1}_{-4.3}^{+4.1}$ (-12%)	${92.1}_{-3.7}^{+3.5}$ (-24%)
$B^0\to {\rho }^0{\rho }^0$	C + TD	${ℬ}^{00},f_L^{00},{𝒞}^{00},{𝒮}^{00}$	${91.8}_{-4.0}^{+4.1}$ (-15%)	$91.7\pm 3.1$ (-35%)
	C	${ℬ}^{00},f_L^{00}$	${91.8}_{-4.5}^{+4.9}$ (-1%)	${92.1}_{-4.7}^{+4.6}$ (-2%)
	C + TD	${𝒞}^{00},{𝒮}^{00}$	${91.8}_{-4.0}^{+4.1}$ (-15%)	$91.6\pm 3.1$ (-35%)
$B^{+}\to {\rho }^{+}{\rho }^0$	C	${ℬ}^{+0}$	${92.2}_{-4.7}^{+4.6}$ (-2%)	${92.5}_{-4.7}^{+4.3}$ (-5%)
$B^{\pm ,0}\to {\left(\rho \rho \right)}^{\pm ,0}$	C + TI	${ℬ}^{ij},f_L^{ij},{𝒞}^{ij},{𝒮}^{ij}$	${92.0}_{-3.4}^{+3.3}$ (-29%)	-

Table 10. Table 13: 68% CL interval for α 𝛼 \alpha from the S ​ U ​ ( 2 ) 𝑆 𝑈 2 {SU}\mkern-2.0mu(2) isospin analysis of the B → π ​ π → 𝐵 𝜋 𝜋 B\rightarrow\pi\pi data in case the uncertainty of some observables is reduced by a factor 2 2 \sqrt{2} (third column) or set to zero (fourth column). The relative gain in resolution with respect to the current measurement is indicated within brackets. Only the preferred solution for α 𝛼 \alpha close to 90 ∘ superscript 90 90^{\circ} is reported.

Decay	Analysis	Improved data	${\sigma }_{𝒪}\to {\sigma }_{𝒪}^{WA}/\sqrt{2}$	${\sigma }_{𝒪}\to 0$
$B^0\to {\pi }^{+}{\pi }^{-}$	C + TD	${ℬ}^{+-},{𝒞}^{+-},{𝒮}^{+-}$	$93.0\pm 14.0$ (0%)	$93.0\pm 14.0$ (0%)
	C	${ℬ}^{+-}$	$93.0\pm 14.0$ (0%)	$93.0\pm 14.0$ (0%)
	TD	${𝒞}^{+-},S^{+-}$	$93.0\pm 14.0$ (0%)	$93.0\pm 14.0$ (0%)
$B^0\to {\pi }^0{\pi }^0$	C + TI	${ℬ}^{00},{𝒞}^{00}$	$92.6\pm 13.0$ (-8%)	${84.1}_{-3.5}^{+4.5}$ (-70%)
	C	${ℬ}^{00}$	$93.0\pm 14.0$ (0%)	$93.5\pm 13.5$ (-3.6%)
	TI	${𝒞}^{00}$	$92.0\pm 13.0$ (-8%)	${84.1}_{-5.6}^{+6.4}$ (-57%)
$B^{+}\to {\pi }^{+}{\pi }^0$	C	${ℬ}^{+0}$	$93.0\pm 14.0$ (0%)	$93.0\pm 14.0$ (0%)
$B^{\pm ,0}\to {\left(\pi \pi \right)}^{\pm ,0}$	C + TD or TI	${ℬ}^{ij},{𝒞}^{+-},{𝒮}^{+-},{𝒞}^{00}$	$92.5\pm 12.5$ (-11%)	-

Table 11. Table 14: Indirect S ​ U ​ ( 2 ) 𝑆 𝑈 2 {SU}\mkern-2.0mu(2) isospin determination of B → π ​ π → 𝐵 𝜋 𝜋 B\!\rightarrow\pi\pi observables compared to their experimental measurement (world average).

Observable	Experimental value	$SU\left(2\right)$ constraint (dir. measurement not in fit)		$\mathrm{Pull}$
	(world average)	68% CL	95% CL
${ℬ}_{\pi \pi }^{\text{+-}}$ $(\times {10}^6)$	$5.10\pm 0.19$	${5.9}_{-1.3}^{+2.2}$ $\cup$ ${15.4}_{-2.2}^{+1.7}$	${15.4}_{-12.0}^{+3.8}$	0.
${𝒞}_{\pi \pi }^{\text{+-}}$	$-0.304\pm 0.047$	$-{0.2}_{-0.2}^{+0.3}$ $\cup$ $0.5\pm 0.1$	$-{0.2}_{-0.3}^{+0.8}$	0.
${𝒮}_{\pi \pi }^{\text{+-}}$	$-0.662\pm 0.062$	$-{0.6}_{-0.2}^{+0.3}$ $\cup$ ${0.5}_{-0.3}^{+0.2}$	$-{0.6}_{-0.3}^{+0.5}$ $\cup$ ${0.5}_{-0.5}^{+0.3}$	0.
${ℬ}_{\pi \pi }^{\text{+0}}$ $(\times {10}^6)$	$5.48\pm 0.34$	${5.0}_{-0.9}^{+0.7}$ $\cup$ ${0.4}_{-0.1}^{+0.3}$	${5.0}_{-4.9}^{+1.3}$	0.
${ℬ}_{\pi \pi }^{\text{00}}$ $(\times {10}^6)$	$1.59\pm 0.18$	${1.3}_{-0.3}^{+1.0}$ $\cup$ ${9.9}_{-6.8}^{+4.7}$	$>0.7$	0.
${𝒞}_{\pi \pi }^{\text{00}}$	$-0.34\pm 0.22$	$-{0.49}_{-0.08}^{+0.12}$ $\cup$ ${0.82}_{-0.06}^{+0.05}$	$-{0.49}_{-0.16}^{+0.34}$ $\cup$ ${0.82}_{-0.18}^{+0.10}$	0.

Table 12. Table 15: Indirect S ​ U ​ ( 2 ) 𝑆 𝑈 2 {SU}\mkern-2.0mu(2) isospin determination of the B → ρ ​ ρ → 𝐵 𝜌 𝜌 B\!\rightarrow\rho\rho observables compared to their experimental measurement (world average).

Observable	Experimental value	$SU\left(2\right)$ constraint (dir. measurement not in fit)		$\mathrm{Pull}$
	(world average)	68% CL	95% CL
${ℬ}_{\rho \rho }^{\text{+-}}$ $(\times {10}^6)$	$27.76\pm 1.84$	${29.7}_{-4.1}^{+31.7}$	${29.7}_{-7.2}^{+37.9}$	0.39
$f_L^{\text{+-}}$	$0.990\pm 0.020$	$>0.91$	$>0.79$	0.22
${𝒞}_{{\rho }_L{\rho }_L}^{\text{+-}}$	$-0.00\pm 0.09$	$-{0.03}_{-0.17}^{+0.19}$	$-{0.03}_{-0.31}^{+0.35}$	0.14
${𝒮}_{{\rho }_L{\rho }_L}^{\text{+-}}$	$-0.15\pm 0.13$	$-0.14\pm 0.15$	$-{0.14}_{-0.33}^{+0.30}$	0.00
${ℬ}_{\rho \rho }^{\text{+0}}$ $(\times {10}^6)$	$24.9\pm 1.9$	${21.4}_{-13.0}^{+2.0}$	${21.4}_{-14.2}^{+3.6}$	0.39
$f_L^{\text{+0}}$	$0.950\pm 0.016$	$>0.34$	$>0.29$	0.39
${ℬ}_{\rho \rho }^{\text{00}}$ $(\times {10}^6)$	$0.93\pm 0.14$	$>0.6$	$>0.2$	0.37
$f_L^{\text{00}}$	$0.71\pm 0.06$	$>0.48$	$>0.20$	0.37
${𝒞}_{{\rho }_L{\rho }_L}^{\text{00}}$	$0.2\pm 0.9$	$0.02\pm 0.46$	${0.02}_{-0.81}^{+0.79}$	0.18
${𝒮}_{{\rho }_L{\rho }_L}^{\text{00}}$	$0.3\pm 0.7$	${0.21}_{-0.55}^{+0.48}$	${0.21}_{-1.10}^{+0.76}$	0.10

Table 13. Table 16: Indirect S ​ U ​ ( 2 ) 𝑆 𝑈 2 {SU}\mkern-2.0mu(2) isospin determination of the 𝒰 𝒰 {\cal U} and ℐ ℐ {\cal I} coefficients compared to their experimental measurement (world average).

Obs.	Experimental value	$SU\left(2\right)$ constraint (dir. measurement not in fit)		$\mathrm{Pull}$
	(world average)	68% CL	95% CL
${𝒰}_{\text{+}}^{\text{-}}$	$0.24\pm 0.09$	${0.98}_{-0.15}^{+0.04}$	${0.98}_{-0.60}^{+0.06}$	2.1
${ℐ}_{\text{+}}$	$0.04\pm 0.06$	${0.48}_{-0.10}^{+0.03}$	${0.48}_{-1.00}^{+0.06}$	1.6
${𝒰}_{\text{-}}^{\text{-}}$	$-0.37\pm 0.09$	$-{1.04}_{-0.16}^{+0.32}$	$-{1.04}_{-0.24}^{+2.35}$	1.8
${𝒰}_{\text{-}}^{\text{+}}$	$1.19\pm 0.07$	${0.38}_{-0.14}^{+0.24}$	${0.38}_{-0.25}^{+0.62}$	2.2
${ℐ}_{\text{-}}$	$0.03\pm 0.06$	$-{0.38}_{-0.16}^{+0.18}$	$-{0.38}_{-0.24}^{+0.34}$	2.3
${𝒰}_{\text{0}}^{\text{+}}$	$0.22\pm 0.03$	${0.09}_{-0.07}^{+0.09}$ $\cup$ ${0.31}_{-0.05}^{0.04}$	${0.09}_{-0.11}^{+0.32}$	1.3
${𝒰}_{\text{0}}^{\text{-}}$	$0.04\pm 0.06$	$-{0.13}_{-0.08}^{+0.11}$	$-{0.13}_{-0.14}^{+0.24}$ $\cup$ ${0.21}_{0.02}^{0.01}$	1.4
${ℐ}_{\text{0}}$	$-0.03\pm 0.04$	$-0.07\pm 0.03$	$-{0.07}_{-0.07}^{+0.22}$	1.0
${𝒰}_{\text{+0}}^{\text{+}ℛe}$	$0.08\pm 0.23$	$-0.49\pm 0.04$	$-0.49\pm 0.08$ $\cup$ ${0.43}_{-0.04}^{+0.02}$	2.4
${𝒰}_{\text{+0}}^{\text{+}ℐm}$	$0.27\pm 0.21$	$-{0.40}_{-0.08}^{+0.19}$	$-{0.40}_{-0.13}^{+0.88}$	1.7
${𝒰}_{\text{+0}}^{\text{-}ℛe}$	$-0.06\pm 0.47$	$-0.18\pm 0.25$ $\cup$ ${0.01}_{-0.02}^{+0.06}$	$-0.18\pm 0.10$ $\cup$ ${0.01}_{-0.05}^{+0.22}$	0.3
${𝒰}_{\text{+0}}^{\text{-}ℐm}$	$-0.12\pm 0.47$	${0.18}_{-0.12}^{+0.06}$	${0.18}_{-0.22}^{+0.12}$	1.0
${ℐ}_{\text{+0}}^{ℛe}$	$0.48\pm 0.78$	${0.06}_{-0.06}^{+0.10}$	${0.06}_{-0.55}^{+0.19}$	0.8
${ℐ}_{\text{+0}}^{ℐm}$	$0.01\pm 0.53$	$-{0.10}_{-0.04}^{+0.19}$	$-{0.10}_{-0.10}^{+0.63}$	0.6
${𝒰}_{\text{-0}}^{\text{+}ℛe}$	$0.07\pm 0.22$	$-0.50\pm 0.05$	$-0.50\pm 0.09$ $\cup$ ${0.36}_{-0.03}^{+0.05}$	2.3
${𝒰}_{\text{-0}}^{\text{+}ℐm}$	$-0.61\pm 0.28$	${0.36}_{-0.34}^{+0.12}$	${0.36}_{-0.80}^{+0.82}$ $\cup$ $-0.46\pm 0.14$	1.9
${𝒰}_{\text{-0}}^{\text{-}ℛe}$	$0.33\pm 0.47$	$-{0.13}_{-0.18}^{+0.23}$	$-{0.13}_{-0.42}^{+0.48}$	0.8
${𝒰}_{\text{-0}}^{\text{-}ℐm}$	$0.69\pm 0.60$	$0.17\pm 0.06$ $\cup$ $0.51\pm 0.05$	${0.17}_{-0.17}^{+0.11}$ $\cup$ ${0.51}_{-0.21}^{+0.09}$	0.8
${ℐ}_{\text{-0}}^{ℛe}$	$-0.46\pm 0.76$	$-{0.15}_{-0.07}^{+0.40}$	$-{0.15}_{-0.13}^{+0.55}$	0.8
${ℐ}_{\text{-0}}^{ℐm}$	$-0.57\pm 0.57$	${0.09}_{-0.07}^{+0.06}$	${0.09}_{-0.23}^{+0.11}$ $\cup$ $0.43\pm 0.08$	1.3
${𝒰}_{\text{+-}}^{\text{+}ℛe}$	$0.08\pm 0.37$	$-{1.03}_{-0.05}^{+0.10}$	$-{1.03}_{-0.10}^{+0.32}$	3.3
${𝒰}_{\text{+-}}^{\text{+}ℐm}$	$0.17\pm 0.34$	${1.04}_{-0.08}^{+0.04}$	${1.04}_{-0.29}^{+0.09}$ $\cup$ $-{1.00}_{-0.06}^{+0.14}$	2.9
${𝒰}_{\text{+-}}^{\text{-}ℛe}$	$-0.43\pm 0.81$	${0.44}_{-0.54}^{+0.35}$	${0.44}_{-1.47}^{+0.57}$	0.8
${𝒰}_{\text{+-}}^{\text{-}ℐm}$	$0.68\pm 0.84$	$-0.03\pm 0.07$	$-0.03\pm 0.14$ $\cup$ ${0.97}_{-0.35}^{+0.10}$	1.5
${ℐ}_{\text{+-}}^{ℛe}$	$-0.34\pm 1.06$	$-{0.28}_{-0.08}^{+0.11}$ $\cup$ ${0.33}_{-0.42}^{0.35}$	$-{0.28}_{-0.36}^{+1.42}$	0.4
${ℐ}_{\text{+-}}^{ℐm}$	$-0.66\pm 1.05$	${1.04}_{-0.17}^{+0.07}$	${1.04}_{-0.46}^{+0.11}$ $\cup$ $-{0.02}_{-0.25}^{+0.44}$	1.8

Table 14. Table 17: Indirect S ​ U ​ ( 2 ) 𝑆 𝑈 2 {SU}\mkern-2.0mu(2) isospin determination of the B 0 → π + ​ π − ​ π 0 → superscript 𝐵 0 superscript 𝜋 superscript 𝜋 superscript 𝜋 0 B^{0}\rightarrow\pi^{+}\pi^{-}\pi^{0} parameters derived from the combined 𝒰 𝒰 {\cal U} and ℐ ℐ {\cal I} coefficients, compared to their experimental measurement (world average).

Observable	Experimental value	$SU\left(2\right)$ constraint (dir. measurement not in fit)		$\mathrm{Pull}$
	(world average)	68% CL	95% CL
${𝒞}^{+}$	$0.24\pm 0.10$	$-{0.86}_{-0.13}^{+0.28}$	$-{0.86}_{-0.14}^{+0.70}$ $\cup$ ${0.87}_{-0.24}^{+0.12}$	2.4
${𝒮}^{+}$	$0.07\pm 0.12$	${0.95}_{-0.29}^{+0.04}$	$[-1.,+1]$	1.4
${𝒞}^{-}$	$-0.31\pm 0.08$	$-{0.99}_{-0.01}^{+0.25}$	$-{0.99}_{-0.01}^{+1.10}$	1.8
${𝒮}^{-}$	$0.04\pm 0.10$	$-{0.68}_{-0.25}^{+0.32}$	$-{0.85}_{-0.15}^{+0.77}$	2.2
${𝒞}^0$	$0.19\pm 0.25$	$-{0.62}_{-0.25}^{+0.44}$	$-{0.62}_{-0.36}^{+0.98}$	1.4
${𝒮}^0$	$-0.28\pm 0.36$	$-{0.70}_{-0.16}^{+0.96}$	$[-1.,+1]$	0.4
$f^{00}$	$0.09\pm 0.03$	${0.13}_{-0.03}^{+0.03}$	${0.13}_{-0.13}^{+0.07}$	0.7
$A^{+}$	$0.11\pm 0.06$	$-{0.93}_{-0.07}^{+0.35}$	$-{0.93}_{-0.07}^{+0.74}$	2.8
$A^{-}$	$-0.04\pm 0.09$	${0.94}_{-0.33}^{+0.06}$	${0.94}_{-0.70}^{+0.06}$	2.8
$𝒞$	$-0.04\pm 0.06$	${0.36}_{-0.16}^{+0.13}$	${0.36}_{-1.12}^{+0.49}$	1.2
$\Delta 𝒞$	$0.28\pm 0.07$	${0.96}_{-0.15}^{+0.04}$	${0.96}_{-0.42}^{+0.04}$ $\cup$ $-{0.96}_{-0.04}^{+0.10}$	2.6
$𝒮$	$0.06\pm 0.08$	${0.34}_{-0.23}^{+0.23}$	$-{0.34}_{-0.46}^{+1.30}$	1.6
$\Delta 𝒮$	$0.01\pm 0.08$	${0.82}_{-0.23}^{+0.13}$	${0.82}_{-0.52}^{+0.18}$	2.7
$A_{\rho \pi }$	$-0.09\pm 0.04$	${0.93}_{-0.33}^{+0.07}$	$-{0.93}_{-0.72}^{+0.07}$	2.8

Table 15. Table 18: Indirect S ​ U ​ ( 2 ) 𝑆 𝑈 2 {SU}\mkern-2.0mu(2) isospin determination of the branching fractions of the charged and neutral B → ρ ​ π → 𝐵 𝜌 𝜋 B\!\rightarrow\rho\pi decays and direct C ​ P 𝐶 𝑃 CP asymmetries for the charged modes compared to their experimental measurement (world average).

Observable	Experimental value	$SU\left(2\right)$ constraint (dir. measurement not in fit)		$\mathrm{Pull}$
	(world average)	68% CL	95% CL
${ℬ}_{\rho \pi }^{\pm \mp }$ $(\times {10}^6)$	$23.0\pm 2.3$	${15.5}_{-3.1}^{+3.4}$	${15.5}_{-6.1}^{+12.5}$	1.6
${ℬ}_{\rho \pi }^{\text{00}}$ $(\times {10}^6)$	$2.0\pm 0.51$	${2.3}_{-0.5}^{+0.4}$	$2.3\pm 0.9$	0.5
${ℬ}_{\rho \pi }^{\text{0+}}$ $(\times {10}^6)$	${8.3}_{-1.3}^{1.2}$	${22.3}_{-8.5}^{+67.6}$		1.9
${ℬ}_{\rho \pi }^{\text{+0}}$ $(\times {10}^6)$	${10.9}_{-1.5}^{1.3}$	${55.0}_{-36.0}^{+31.0}$	${55}_{-44}^{+47}$	2.0
${𝒞}_{\rho \pi }^{\text{0+}}$	$-{0.18}_{0.09}^{0.17}$	${0.09}_{-0.13}^{+0.10}$ $\cup$ $-{0.13}_{-0.04}^{+0.02}$	${0.09}_{-0.56}^{+0.22}$ $\cup$ ${0.80}_{-0.16}^{+0.26}$	0.8
${𝒞}_{\rho \pi }^{\text{+0}}$	$-0.02\pm 0.11$	$-{0.17}_{-0.08}^{+0.10}$	$-{0.17}_{-0.18}^{+0.08}$ $\cup$ $-{0.95}_{-0.05}^{+0.14}$	1.0

Equations

$${\alpha }_{\mathrm{dir}}=(86.2_{-4.0}^{+4.4}\cup 178.4_{-5.1}^{+3.9}{)}^{\circ }.$$

$${\alpha }_{\mathrm{ind}}=(92.5_{-1.1}^{+1.5})\text{∘}.$$

$${\mathcal{S}}_{\pi \pi }^{00}=0.65\pm \mathrm{0.13.}$$

$${\mathcal{A}}^{ij}=⟨h_1^i{h}_2^j|{\mathcal{H}}_{\mathrm{eff}}|B_d⟩=V_{ud}{V}_{ub}^{\ast }({\mathcal{T}}_u^{ij}+{\mathcal{P}}_u^{ij})+V_{cd}{V}_{cb}^{\ast }{\mathcal{P}}_c^{ij}+V_{td}{V}_{tb}^{\ast }{\mathcal{P}}_t^{ij}$$

$$\mathcal{U}-\mathrm{convention}:{\mathcal{A}}^{ij}$$

$$\mathcal{C}-\mathrm{convention}:{\mathcal{A}}^{ij}$$

$$\mathcal{T}-\mathrm{convention}:{\mathcal{A}}^{ij}$$

$${\mathcal{A}}^{ij}=V_{ud}{V}_{ub}^{\ast }{\tilde{\mathcal{T}}}^{ij}+V_{td}{V}_{tb}^{\ast }{\tilde{\mathcal{P}}}^{ij}$$

$${\mathcal{A}}^{ij}=-e^{\mathsf{i}\gamma }{T}^{ij}+e^{-\mathsf{i}\beta }{P}^{ij},$$

$$\frac{p}{q}{\bar{\mathcal{A}}}^{ij}=-e^{-\mathsf{i}\gamma }{T}^{ij}+e^{\mathsf{i}\beta }{P}^{ij},$$

$$e^{2\mathsf{i}\alpha }=\frac{e^{\mathsf{i}\alpha }{T}^{ij}}{e^{-\mathsf{i}\alpha }{T}^{ij}}=\frac{{\bar{A}}^{ij}-P^{ij}}{A^{ij}-P^{ij}}.$$

$$a_{CP}(t)=\frac{\Gamma ({\bar{B}}^0(t)\to h_1^{+}{h}_2^{-})-\Gamma (B^0(t)\to h_1^{+}{h}_2^{-})}{\Gamma ({\bar{B}}^0(t)\to h_1^{+}{h}_2^{-})+\Gamma (B^0(t)\to h_1^{+}{h}_2^{-})}={\mathcal{S}}^{+-}\sin (\Delta m_{d}t)-{\mathcal{C}}^{+-}\cos (\Delta m_{d}t),$$

$$\lambda =\frac{{\bar{A}}^{+-}}{A^{+-}},{\mathcal{S}}^{+-}=\frac{2\mathrm{Im}\lambda }{1+|\lambda {|}^2},{\mathcal{C}}^{+-}=\frac{1-|\lambda {|}^2}{1+|\lambda {|}^2}.$$

$$\lambda =|\lambda |e^{2\mathsf{i}{\alpha }_{\mathrm{eff}}},{\mathcal{S}}^{+-}=\sqrt{1-({\mathcal{C}}^{+-}{)}^2}\sin (2{\alpha }_{\mathrm{eff}}).$$

$$A^{+-}+A^{-+}+2A^{00}=\sqrt{2}(A^{+0}+A^{0+}).$$

$$P^{+-}+P^{-+}+2P^{00}=\sqrt{2}(P^{+0}+P^{0+})=0,$$

$$e^{2\mathsf{i}\alpha }=\frac{{\bar{A}}^{+0}+{\bar{A}}^{0+}}{A^{+0}+A^{0+}}=\frac{{\bar{A}}^{+-}+{\bar{A}}^{-+}+2{\bar{A}}^{00}}{A^{+-}+A^{-+}+2A^{00}}.$$

$$P^{+0}-P^{0+}=\sqrt{2}(P^{+-}-P^{-+}),$$

$$e^{2\mathsf{i}\alpha }=\frac{({\bar{A}}^{+0}-{\bar{A}}^{0+})-\sqrt{2}({\bar{A}}^{+-}-{\bar{A}}^{-+})}{(A^{+0}-A^{0+})-\sqrt{2}(A^{+-}-A^{-+})}.$$

$$P^{+0}=-P^{0+}=(P^{+-}-P^{-+})/\sqrt{2},P^{00}=-(P^{+-}+P^{-+})/2,$$

$$A_{h_1=h_2}^{+-}=\frac{A^{+-}+A^{-+}}{\sqrt{2}},A_{h_1=h_2}^{+0}=\frac{A^{+0}+A^{0+}}{\sqrt{2}},$$

$$A^{+0}=\frac{A^{+-}}{\sqrt{2}}+A^{00},$$

$$e^{-2\mathsf{i}\alpha }=\frac{{\bar{A}}^{+0}}{A^{+0}}.$$

$${\chi }^2(\vec{p})=-2\log \mathcal{L}({\vec{\mathcal{O}}}_{exp}-\vec{\mathcal{O}}(\vec{p})),$$

$$\Delta {\chi }^2(\alpha )=\underset{\vec{\mu }}{\min }[{\chi }^2(\alpha )]-{\chi }_{\min }^2$$

$$p(\alpha )=\mathrm{Prob}(\Delta {\chi }^2(\alpha ),N_{dof}=1),\mathrm{Prob}\left(\Delta {\chi }^2,N_{dof}\right)=\frac{\Gamma (N_{dof}/2,\Delta {\chi }^2/2)}{\Gamma (N_{dof}/2)}.$$

$$\mathrm{Prob}(\Delta {\chi }^2(\alpha ),N_{dof}=1)=\mathrm{Erfc}(\sqrt{\mathrm{\Delta }{\chi }^2(\alpha )/2}).$$

$$\frac{1}{{\tau }_{B^{i+j}}}{\mathcal{B}}^{ij}=\frac{|A^{ij}{|}^2+|{\bar{A}}^{ij}{|}^2}{2},{\mathcal{C}}^{ij}=\frac{|A^{ij}{|}^2-|{\bar{A}}^{ij}{|}^2}{|A^{ij}{|}^2+|{\bar{A}}^{ij}{|}^2},{\mathcal{S}}^{ij}=\frac{2\mathcal{I}m({\bar{A}}^{ij}{A}^{ij}{}^{\ast })}{|A^{ij}{|}^2+|{\bar{A}}^{ij}{|}^2},$$

$$A^{+-}=T^{+-}{e}^{-\mathsf{i}\alpha }+P$$

$$\sqrt{2}{A}^{00}=T^{00}{e}^{-\mathsf{i}\alpha }-P$$

$$\sqrt{2}{A}^{+0}=(T^{+-}+T^{00})e^{-\mathsf{i}\alpha }$$

$$A^{+0}=\mu e^{\mathsf{i}(\Delta -\alpha )}$$

Full text

LPT-Orsay-17-07

Isospin analysis of charmless $B$-meson decays

J. Charlesa, O. Deschampsb ††† Contact author: Olivier Deschamps, Olivier.Deschamps@cern.ch , S. Descotes-Genonc, V. Niessb

[for the CKMfitter group]

a* CNRS, Aix Marseille Univ, Université de Toulon, CPT

13288 Marseille Cedex 9, France

b Laboratoire de Physique de Clermont, CNRS/Université Clermont Auvergne, UMR 6533,

Campus des Cézeaux, 24 Av. des Landais, F-63177 Aubière Cedex, France

c Laboratoire de Physique Théorique (UMR 8627), CNRS, Univ. Paris-Sud,

Université Paris-Saclay, 91405 Orsay Cedex, France *

We discuss the determination of the CKM angle $\alpha$ using the non-leptonic two-body decays $B\rightarrow\pi\pi$, $B\rightarrow\rho\rho$ and $B\rightarrow\rho\pi$ using the latest data available. We illustrate the methods used in each case and extract the corresponding value of $\alpha$. Combining all these elements, we obtain the determination $\alpha_{\rm dir}=(86.2^{\scriptscriptstyle+4.4}_{\scriptscriptstyle-4.0}\cup 178.4^{\scriptscriptstyle+3.9}_{\scriptscriptstyle-5.1})^{\circ}$. We assess the uncertainties associated to the breakdown of the isospin hypothesis and the choice of statistical framework in detail. We also determine the hadronic amplitudes (tree and penguin) describing the QCD dynamics involved in these decays, briefly comparing our results with theoretical expectations. For each observable of interest in the $B\rightarrow\pi\pi$, $B\rightarrow\rho\rho$ and $B\rightarrow\rho\pi$ systems, we perform an indirect determination based on the constraints from all the other observables available and we discuss the compatibility between indirect and direct determinations. Finally, we review the impact of future improved measurements on the determination of $\alpha$.

Published in Eur.Phys.J. C77 (2017) no.8, 574

1 Introduction

Over the last few decades, our understanding of $CP$ violation has made great progress, with many new constraints from BaBar, Belle and LHCb experiments among others [1, 2]. These constraints were shown to be in remarkable agreement with each other and to support the Kobayashi–Maskawa mechanism of $CP$ violation at work within the Standard Model (SM) with three generations [3, 4]. This has led to an accurate determination of the Cabibbo–Kobayashi–Maskawa matrix (CKM) encoding the pattern of $CP$ violation as well as the strength of the weak transitions among quarks of different generations [5, 6, 7, 8, 9, 10]. These constraints prove also essential in assessing the viability of New Physics models with well-motivated flavour structures [11, 12, 13, 14].

As the CKM matrix is related to quark-flavour transitions, most of these constraints are significantly affected by hadronic uncertainties due to QCD binding quarks into the observed hadrons. However, some of these constraints have the very interesting feature of being almost free from such uncertainties. This is in particular the case for the constraints on the CKM angle $\alpha$ that are derived from the isospin analysis of the charmless decay modes $B\rightarrow\pi\pi$, $B\rightarrow\rho\rho$ and $B\rightarrow\rho\pi$. Indeed, assuming the isospin symmetry and neglecting the electroweak penguin contributions, the amplitudes of the ${SU}\mkern-2.0mu(2)$-conjugated modes are related. The measured branching fractions and asymmetries in the ${\cal B}^{\pm,0}\rightarrow(\pi\pi)^{\pm,0}$ and $B^{\pm,0}\rightarrow(\rho\rho)^{\pm,0}$ modes and the bilinear form factors in the Dalitz analysis of the $B^{0}\rightarrow(\rho\pi)^{0}$ decays provide enough observables to simultaneously extract the weak phase $\beta+\gamma=\pi-\alpha$ together with the hadronic tree and penguin contributions to each mode [15, 16, 17]. Therefore, these modes probe two different corners of the SM: on one side, they yield information on $\alpha$ that is a powerful constraint on the Kobayashi–Maskawa mechanism (and the CKM matrix) involved in weak interactions, and on the other side, they provide a glimpse on the strong interaction and especially the hadronic dynamics of charmless two-body $B$-decays.

In the following, we will provide a thorough analysis of these decays, based on the data accumulated up to the conferences of Winter 2017. Combining the experimental data for the three decay modes above, we obtain the world-average value at 68% Confidence Level (CL):

[TABLE]

The solution near $90^{\circ}$ is in good agreement with the indirect determination obtained by the global fit of the flavour data performed by the CKMfitter group in Summer 2016 [18]:

[TABLE]

Considered separately, the $B\rightarrow\pi\pi$ and $B\rightarrow\rho\rho$ decays yield direct determinations of $\alpha$ in very good agreement with the indirect determination Eq. (2), whereas the $B\rightarrow\rho\pi$ decay exhibits a 3$\sigma$ discrepancy. This discrepancy affects only marginally the combination Eq. (1), which is dominated by the results from $\rho\rho$ decays, and to a lesser extent $\pi\pi$ decays, whereas $\rho\pi$ modes play only a limited role. At this level of accuracy, it proves interesting to assess uncertainties neglected up to now, namely the sources of violation of the assumptions underlying these determinations (\Delta I=\raisebox{1.72218pt}{\footnotesize3}\negthinspace/\negthinspace\raisebox{-2.58334pt}{\footnotesize2} electroweak penguins, $\pi^{0}-\eta-\eta^{\prime}$ mixing, $\rho$ width) and the role played by the statistical framework used to extract the confidence intervals. These effects may shift the central value of $\alpha_{\rm dir}$ by around $2^{\circ}$, while keeping the uncertainty around $4^{\circ}$ to $5^{\circ}$, thus remaining within the statistical uncertainty quoted in Eq. (1).

Besides the CKM angle $\alpha$, the isospin analysis of the charmless decays data provides a determination of the hadronic tree and penguin parameters for each mode. The penguin-to-tree and colour-suppression ratios in the $B\!\rightarrow\pi\pi$ and $B\!\rightarrow\rho\rho$ decay modes can be determined precisely: they show overall good agreement with theoretical expectations for $\rho\rho$ modes, whereas the ratio between colour-allowed and -suppressed tree contributions for $\pi\pi$ modes do not agree well with theoretical expectations (both for the modulus and the phase). We also perform indirect predictions of the experimental observables, using all the other available measurements to predict the value of a given observable, and we can compare our predictions to the existing measurements: a very good compatibility is observed for the $\pi\pi$ and $\rho\rho$ channels, whereas discrepancies occur in some of the observables describing the Dalitz plot for $B\rightarrow\rho\pi$ decays. Among many other quantities, the yet-to-be-measured mixing-induced $CP$ asymmetry in the $B^{0}\!\rightarrow\pi^{0}\pi^{0}$ decay is predicted at 68% CL:

[TABLE]

using as an input the indirect value of $\alpha$, see Eq. (2). Finally, we study how the determination of $\alpha$ would be affected if the accuracy of specific subsets of observables is improved through new measurements. We find that an improved accuracy for the time-dependent asymmetries in $B^{0}\rightarrow\rho^{0}\rho^{0}$ and the measurement of ${\cal S}^{00}_{\pi\pi}$ would reduce the uncertainty on $\alpha$ in a noticeable way.

The rest of this article goes as follows. In Sec. 2, we discuss the basics of isospin analysis for charmless $B$-meson decays. In Sec. 3, we provide details on the extraction of the $\alpha$ angle, focusing on the $B\rightarrow\pi\pi$, $B\rightarrow\rho\rho$ and $B\rightarrow\rho\pi$ modes in turn, before combining these extractions in a world average. In Sec. 4, we consider the uncertainties attached to the extraction of $\alpha$. First, we discuss the ${SU}\mkern-2.0mu(2)$ isospin framework underlying these analyses, considering three sources of corrections: the presence of \Delta I=\raisebox{1.72218pt}{\footnotesize3}\negthinspace/\negthinspace\raisebox{-2.58334pt}{\footnotesize2} electroweak penguins (isospin breaking due to different charges for the $u$ and $d$ quarks), $\pi^{0}-\eta-\eta^{\prime}$ mixing (isospin breaking due to different masses for the $u$ and $d$ quarks), $\rho$ width (additional amplitude to include in the isospin relations for $B\rightarrow\rho\rho$). In addition, we discuss the statistical issues related to our frequentist framework. In Sec. 5, we extract the hadronic (tree and penguin) amplitudes associated with each decay, comparing them briefly with theoretical expectations. We use our framework to perform indirect predictions of observables in each channel and discuss the compatibility with the available measurements in Sec. 6. We perform a prospective study to determine the impact of reducing the experimental uncertainties on specific observables in Sec. 7, before drawing our conclusions. Dedicated appendices gather additional numerical results for observables in three modes, separate analyses using either the BaBar or Belle inputs only, and a brief discussion of the quasi-two-body analysis of the charmless $B^{0}\rightarrow a_{1}^{\pm}\pi^{\mp}$ that may provide some further information on $\alpha$.

2 Isospin decomposition of charmless two-body $B$ decays

The decay of neutral $B_{d}$ and charged $B_{u}$ mesons into a pair of light unflavoured isovector mesons $B\rightarrow h_{1}^{i}h_{2}^{j}$ ($h$=$\pi,\rho$ and $i,j=-,0,+$) is described by the weak transition $\bar{b}\rightarrow\bar{u}u\bar{d}$, followed by the hadronisation of the $(\bar{u}u\bar{d},q)$ system, where $q=d(u)$ is the spectator quark in the neutral (charged) component of the mesons isodoublet. The weak process receives dominant contributions from both the tree-level $\bar{b}\rightarrow\bar{u}(u\bar{d})$ charged transition and the flavour-changing neutral current penguin transition, $\bar{b}\rightarrow\bar{d}(u\bar{u})$, whose topologies are shown in Fig. 1.

2.1 Penguin pollution

We can write the charmless $B_{q}\rightarrow(\bar{u}u\bar{d},q)\rightarrow h_{1}^{i}h_{2}^{j}$ transition amplitude, accounting for the CKM factors for the tree and penguin diagrams with three different up-type quark flavours $(u,c,t)$ occurring in the $W$ loop:

[TABLE]

where ${\cal H}_{\rm eff}$ is the effective Hamiltonian describing the transition. The amplitudes ${\cal T}^{ij}_{u}$ and ${\cal P}^{ij}_{u,c,t}$ represents the tree-level and $(u,c,t)$-loop mediated topologies, respectively.

The unitarity of the CKM matrix can be used to eliminate one of the three terms in Eq. (4), resulting in three conventions $\mathcal{U}$, $\mathcal{C}$ and $\mathcal{T}$ defined as:

[TABLE]

The $\mathcal{C}$-convention is adopted in the following, so that the amplitudes can be rewritten as:

[TABLE]

where the terms $\mathcal{\tilde{T}}^{ij}=\mathcal{T}^{ij}_{u}+\mathcal{P}^{ij}_{u}-\mathcal{P}^{ij}_{c}$ and $\mathcal{\tilde{P}}^{ij}=\mathcal{P}^{ij}_{t}-\mathcal{P}^{ij}_{c}$ associated with specific CKM factors will be referred to as ’tree’ and ’penguin’ amplitudes, respectively, following a conventional abuse of language (these amplitudes are not physical, in the sense that both mix contributions of distinct topologies and are not renormalisation invariant). The choice of convention is arbitrary and it has no physical implication on the determination of the weak phase affecting the transitions. However, the particular choice sets the dynamical content of the hadronic tree and penguin amplitudes: for instance, it fixes the definition of the penguin-to-tree ratio and the other hadronic quantities discussed in Sec. 5.

Pulling out the weak phases $\gamma=\arg\left(-\frac{V_{ud}V^{*}_{ub}}{V_{cd}V^{*}_{cb}}\right)$ and $\beta=\arg\left(-\frac{V_{cd}V^{*}_{cb}}{V_{td}V^{*}_{tb}}\right)$, the amplitude can be rewritten as:

[TABLE]

where the magnitude of the CKM products $R_{u}=|V_{ud}V^{*}_{ub}|$ and $R_{t}=|V_{td}V^{*}_{tb}|$ is included in the redefined amplitudes $T^{ij}=R_{u}{\cal\tilde{T}}^{ij}$ and $P^{ij}=-R_{t}{\cal\tilde{P}}^{ij}$.

Similarly, the decay amplitudes of the $CP$-conjugate isodoublet ($\bar{B}^{0}$, $B^{-}$) can be expressed as

[TABLE]

where the factor $p/q\sim e^{2\mathsf{i}\beta}$ is included to take into account the $B^{0}-\bar{B}^{0}$ mixing phase of neutral $B$-meson that arises naturally in physical observables. For consistency between all the ${SU}\mkern-2.0mu(2)$-related decay modes considered in the isospin analysis, the same phase convention has been applied to define the amplitudes for the charged $B$ meson. The $CP$ invariance of the strong interaction means that the same hadronic amplitudes $T^{ij}$ and $P^{ij}$ are involved in the $CP$-conjugate processes, whereas a complex conjugation is applied to the weak phases.

Rotating all amplitudes by the weak phase $\beta$ through the redefinition $A^{ij}=e^{\mathsf{i}\beta}\mathcal{A}^{ij}$ and ${\bar{A}}^{ij}=e^{\mathsf{i}\beta}\mathcal{\bar{A}}^{ij}$, the parameter $\alpha=\pi-\beta-\gamma$ appears as half of the phase difference between the tree contributions to the $CP$-conjugate amplitudes:

[TABLE]

In the absence of penguin contributions, $\alpha$ can be related to the relative phase of $CP$-conjugate amplitudes describing the $B^{0}$ and $\bar{B}^{0}$ mesons decaying into the same final state $h_{1}^{i}h_{2}^{j}$. In particular, the time-dependent analysis of the $B^{0}/{\bar{B}^{0}}\rightarrow h_{1}^{+}h_{2}^{-}$ decay yields the $CP$ asymmetry:

[TABLE]

where $\Delta m_{d}$ is the $B^{0}-\bar{B}^{0}$ oscillation frequency, and $t$ is either the decay time of the meson, or (in the case of $B$-factories) the time difference between the $CP$- and tag-side decays. The coefficients can be expressed as

[TABLE]

Therefore, in the absence of the penguin contributions, the measurement of ${\cal S}^{+-}$ would yield $\sin(2\alpha)$. The a priori non-negligible penguin contributions to the $B$ charmless decays modify this picture: if we introduce the effective angle corresponding to the phase of $\lambda$, we have [19]

[TABLE]

The non-negligible penguin contributions prevent us from identifying $\alpha$ and $\alpha_{\rm eff}$ obtained from the sole consideration of $B^{0}(t)\rightarrow h_{1}^{+}h_{2}^{-}$. However, since isospin is conserved during the hadronisation process, we can write useful relations among the hadronic amplitudes in ${SU}\mkern-2.0mu(2)$-related modes. These relationships will help us to determine the amount of penguin pollution from the data and thus to extract $\alpha$ in spite of hadronic effects [15, 16].

2.2 General isospin decomposition and application to the $\rho\pi$ final-state

One can factorise the decay amplitudes in two parts. First, the weak decay $\bar{b}\rightarrow u\bar{u}\bar{d}$ (common to all $B_{q}\rightarrow h_{1}^{i}h_{2}^{j}$ decay processes) corresponds to a shift of isospin $\Delta I$. The hadronisation into two light mesons can then be described as $\langle h_{1}^{i}h_{2}^{j}|{\cal H}_{s}|u\bar{u}\bar{d},q\rangle$ where ${\cal H}_{s}$ represents the isospin-conserving strong interaction Hamiltonian. The Wigner-Eckart theorem can be used to express these amplitudes in tersm of reduced matrix elements, $A_{\Delta I,I_{f}}$, identified by the shift $\Delta I$ and the final-state isospin $I_{f}$ ($I_{f}=0,1,2$)111Considering only the valence quarks, the isospin shift $\Delta I$ can only take the values \raisebox{1.72218pt}{\footnotesize3}\negthinspace/\negthinspace\raisebox{-2.58334pt}{\footnotesize2} and \raisebox{1.72218pt}{\footnotesize1}\negthinspace/\negthinspace\raisebox{-2.58334pt}{\footnotesize2} in the $\bar{b}\rightarrow(\bar{u}u\bar{d})$ transition. Tree and strong penguin topologies correspond to \Delta I=\raisebox{1.72218pt}{\footnotesize1}\negthinspace/\negthinspace\raisebox{-2.58334pt}{\footnotesize2}, whereas electroweak penguins contain both \Delta I=\raisebox{1.72218pt}{\footnotesize1}\negthinspace/\negthinspace\raisebox{-2.58334pt}{\footnotesize2} and \Delta I=\raisebox{1.72218pt}{\footnotesize3}\negthinspace/\negthinspace\raisebox{-2.58334pt}{\footnotesize2} contributions. For completeness, the possible \Delta I=\raisebox{1.72218pt}{\footnotesize5}\negthinspace/\negthinspace\raisebox{-2.58334pt}{\footnotesize2} contribution due to long-distance rescattering effects [20, 21] is also reported in Tab. 1. This contribution, suppressed by a factor $\alpha_{em}\sim 1/127$, will be neglected until Sec. 4.1.2.. Tab. 1 yields the decomposition of the $A^{ij}$ amplitudes in the general case of two distinguishable isovector mesons ($h_{1}\neq h_{2}$).

This general decomposition applies for instance to the $B\rightarrow\rho\pi$ system. Neglecting the \Delta I=\raisebox{1.72218pt}{\footnotesize5}\negthinspace/\negthinspace\raisebox{-2.58334pt}{\footnotesize2} transition, the four remaining isospin amplitudes constrain the five decays amplitudes to follow the pentagonal relation:

[TABLE]

The same identity applies to the $CP$-conjugate amplitudes. Moreover, the sum of the decay amplitudes of the charged modes $(A^{+0}+A^{0+})$ is a pure $A_{\frac{3}{2},2}$ isospin amplitude.

A usual simplifying assumption consists in neglecting the \Delta I=\raisebox{1.72218pt}{\footnotesize3}\negthinspace/\negthinspace\raisebox{-2.58334pt}{\footnotesize2} contribution from electroweak penguins, as they are expected to be small: in this limit, all penguins (gluonic and electroweak) are mediated only by the \Delta I=\raisebox{1.72218pt}{\footnotesize1}\negthinspace/\negthinspace\raisebox{-2.58334pt}{\footnotesize2} transition $\bar{b}\rightarrow\bar{d}(u\bar{u})_{I=0}$. The isospin relation for the amplitudes Eq. (13) may then be projected onto the penguin amplitudes, leading to the triangular relations:

[TABLE]

Both combinations are identical and vanish under our assumptions, as they involve only \Delta I=\raisebox{1.72218pt}{\footnotesize3}\negthinspace/\negthinspace\raisebox{-2.58334pt}{\footnotesize2} amplitudes. Under this assumption, some combinations of the decay amplitudes are free of penguin contributions. The parameter $\alpha$ can thus be identified as half of the phase difference between the sum of the $CP$-conjugate amplitudes of the charged modes, or equivalently as a function of the amplitudes of the neutral modes only:

[TABLE]

Similarly, the combination of amplitudes $(A^{+0}-A^{0+})-\sqrt{2}(A^{+-}-A^{-+})$ receives a pure $A_{\frac{3}{2},1}$ contribution and is thus free from penguin contributions, so that:

[TABLE]

and the phase $\alpha$ can also be obtained from the pure tree combination:

[TABLE]

Combining Eqs. (14) and (16), one can see that there are actually only two independent penguin amplitudes involved:

[TABLE]

illustrating the power of the ${SU}\mkern-2.0mu(2)$ isospin approach to reduce the number of independent hadronic quantities to be determined from the data in order to extract a constraint on $\alpha$.

2.3 Application to the $\pi\pi$ and $\rho\rho$ cases

In the cases where the two isovector mesons in the final state are indistinguishable from the point of view of isospin symmetry ($h_{1}=h_{2}=h$), as for the $B\rightarrow\pi\pi$ and $B\rightarrow\rho\rho$ systems, only the even amplitudes $I_{f}=0,2$ are allowed due to the Bose–Einstein statistics. Defining the symmetrised amplitudes:

[TABLE]

the isospin decomposition gets simplified as reported in Tab. 2, leading to the triangular identity:

[TABLE]

with a similar identity for the $CP$-conjugate amplitudes 222Unless otherwise stated, the subscript $h_{1}=h_{2}$ is dropped and $A^{ij}$ implicitly refers to the symmetrised amplitude in the case of a decay into two mesons of same type.. In the case of the $\rho\rho$ channel, one should consider a different set of independent amplitudes for each of the three possible polarisations (which is identical for the two $\rho$ mesons).

Under the assumption of negligible contributions from \Delta I=\raisebox{1.72218pt}{\footnotesize3}\negthinspace/\negthinspace\raisebox{-2.58334pt}{\footnotesize2} electroweak penguins, the penguin relations Eq. (18) reduce down to $P^{+0}=0$ and \raisebox{1.72218pt}{\footnotesizeP^{+-}}\negthinspace/\negthinspace\raisebox{-2.58334pt}{\footnotesize\sqrt{2}}=-P^{00}. The total amplitude of the charged modes $A^{+0}=e^{-\mathsf{i}\alpha}T^{+0}$ is free from penguin contributions, and the angle $\alpha$ can be derived from:

[TABLE]

This complex ratio cannot be determined from a single measurement, but it is possible to reconstruct the two isospin triangles corresponding to $CP$-conjugate amplitudes using branching ratios and $CP$ asymmetries for all the modes. Fig. 2 illustrates this construction, which translates the measurement of $\alpha_{\rm eff}$ into a determination of the CKM angle $\alpha$. This procedure is affected by discrete ambiguities, since there are several manners of reconstructing the two isospin triangles. This leads to a fourfold ambiguity for $\sin(2\alpha)$, i.e. an eightfold ambiguity on the solutions of $\alpha$ in $[0,180]^{\circ}$, in general. These additional solutions are called “mirror solutions”. If one or both triangles are flat, several mirror solutions become degenerate, decreasing the number of distinct solutions for $\alpha$.

3 Determining the weak angle $\alpha$

3.1 Procedure

If we consider that each of the five amplitudes $A^{ij}$ $(ij=+-,-+,+0,0+,00)$ receives two complex contributions, tree and penguin, the hadronic contributions to the generic decay system $B^{i+j}\rightarrow h_{1}^{i}h_{2}^{j}$ can be parametrised with 20 real parameters, in addition to the weak phase $\alpha$ (one overall phase being irrelevant). In the case of a Bose-symmetric final state ($h_{1}=h_{2}$), the dimension of the parameter space reduces down to 13. Assuming the isospin relations between the amplitudes discussed in Sec. 2, the three decay systems $B\rightarrow\pi\pi$, $B\rightarrow\rho\rho$ and $B\rightarrow\rho\pi$ provide enough experimental measurements to fully constrain the parameter space, hereafter denoted $\vec{p}$=($\alpha$,$\vec{\mu}$), where $\vec{\mu}$ represents the set of independent hadronic parameters (tree and penguin amplitudes). The actual dimension of the parameter space depends on the decay system and will be discussed in the following subsections.

The constraints on the parameter $\alpha$ are determined by an exploration of the $N$-dimensional parameters space through the frequentist statistical approach discussed in detail in Refs. [5, 6, 10, 22] and in Sec. 4.2, but we find it useful to briefly summarise its main features here. The set of experimental observables, denoted $\vec{\cal O}_{exp}$, is measured in terms of likelihoods that can be used to build a $\chi^{2}$-like test statistic:

[TABLE]

where $\vec{\cal O}(\vec{p})$ represent the theoretical value of the observables in SM for fixed parameters $\vec{p}$. The test statistic $\chi^{2}$ is first minimised over the whole parameter space, letting all $N$ parameters $\vec{p}$ free to vary. The absolute minimum value of the test statistic, $\chi^{2}_{\rm min}$, quantifies the agreement of the data with the theoretical model (assuming the validity of the SM and ${SU}\mkern-2.0mu(2)$ isospin symmetry in the present case). Converting $\chi^{2}_{\rm min}$ into a $p$-value is however not trivial a priori, as one has to interpret $\chi^{2}(\vec{p})$ as a random variable distributed according to a $\chi^{2}$ law with a certain number of degrees of freedom. The actual number of degrees of freedom of the system can be ill defined in the case where the experimental observables are interdependent (see for instance the related discussion in Sec. 3.4). In the case of $M$ independent observables, the number of degrees of freedom of the system is defined as $N_{dof}=M-N$. This occurs in the Gaussian case, but it also can apply in non-Gaussian cases in the limit of large samples, under the conditions of Wilks’ theorem [23].

It is also possible to perform the metrology of specific parameters of the model, and in particular $\alpha$ [5, 6, 10, 22]. Indeed, if we consider the $\vec{\mu}$ hadronic parameters as “nuisance parameters”, we can define the test statistic from the $\chi^{2}$ difference:

[TABLE]

where $\min_{\vec{\mu}}[\chi^{2}(\alpha)]$ is the value of $\chi^{2}$, minimised with respect to the nuisance parameters for a fixed $\alpha$ value. This test statistic assesses how a given hypothesis on the true value of $\alpha$ agrees with the data, irrespective of the value of the nuisance parameters. Confidence intervals on $\alpha$ can be derived from the resulting $p$-value, which is computed assuming that $\Delta\chi^{2}(\alpha)$ is $\chi^{2}$-distributed with one degree of freedom:

[TABLE]

One can express ${\rm Prob}$ for $N_{dof}=1$ in terms of the complementary error function:

[TABLE]

Confidence intervals at a given confidence level (CL) are obtained by selecting the values of $\alpha$ with a $p$-value larger than $1-$CL. The derivation, robustness and coverage of this definition for the $p$-value will be further discussed in Sec. 4.2.

Although the relevant information on the $\alpha$ constraint is fully contained in the $p$-value function $p(\alpha)$, confidence intervals will be derived in the following subsections. It is worth noticing that the $p$-value for $\alpha$ usually presents a highly non-Gaussian profile and that the ${SU}\mkern-2.0mu(2)$ isospin analysis suffers from (pseudo-)mirror ambiguities. Therefore, the confidence intervals provided must be interpreted with particular care.

We consider here all the experimental inputs available up to Winter 2017 conferences (for each channel, we provide the corresponding references). We compare the results of the isospin-based direct determination of $\alpha$ from these inputs with the indirect result from the global CKM fit obtained in Summer 2016. This indirect determination includes all the quark-flavour constraints described in Refs. [10, 18], apart from the inputs from the $B\rightarrow\pi\pi$, $B\rightarrow\rho\rho$ and $B\rightarrow\rho\pi$ decays used for the direct determination of $\alpha$. The indirect determination, hereafter denoted $\alpha_{\rm ind}$, is illustrated in Fig. 3 and the corresponding 68% CL interval is given in Eq. (2).

3.2 Isospin analysis of the $B\rightarrow\pi\pi$ system

The three decays $B^{0,+}\rightarrow(\pi\pi)^{0,+}$ depend on 12 hadronic parameters and the weak phase $\alpha$. One can set one irrelevant global phase to zero, and one can eliminate further parameters using the complex isospin relation Eq. (20) and its $CP$ conjugate (four real constraints) as well as the absence of the penguin contribution to the charged mode amplitude (two real constraints). The remaining system of amplitudes features six degrees of freedom. The isospin-related $B\rightarrow\pi\pi$ decays (and similarly each of the helicity states of the $B\rightarrow\rho\rho$ mode) can thus be described with six real independent parameters, including the common weak phase $\alpha$.

The relevant available experimental observables and their current world-average value for the $B\rightarrow\pi\pi$ modes are summarised333The $CP$ asymmetry for the charged mode explicitly vanishes in the parametrisation of the amplitudes based on exact ${SU}\mkern-2.0mu(2)$ isospin symmetry. Therefore the corresponding observable is not included here. The current experimental measurement, fully consistent with the null hypothesis, would only affect the minimal value $\chi^{2}_{\rm min}$ of the fit, but not the metrology of the parameters. This additional observable will be useful to study isospin-breaking effects in Sec. 4.1. in Tab. 3. Six independent observables are available, allowing us to constrain the six-dimensional parameter space of the ${SU}\mkern-2.0mu(2)$ isospin analysis. These branching fractions and ${CP}$ asymmetries are related to the decay amplitudes as:

[TABLE]

where $\tau_{B^{i+j}}$ is the measured lifetime of the charged ($i+j=1$) or neutral ($i+j=0$) $B$ meson.

Taking into account the isospin relation Eq. (20), the amplitudes can be parametrised in terms of their penguin and tree contributions as

[TABLE]

where $T^{+-}$, $T^{00}$ and $P$ are three complex parameters, among which one can be taken as real to set the global phase convention.

The actual choice of the representation for the amplitudes is irrelevant from the mathematical point of view in the frequentist approach adopted here [36]. We use the following alternative representation of the amplitude system, which proves more convenient for our purposes [37]:

[TABLE]

where the parameters $a$, ${\bar{a}}$ and $\mu$ are real positive parameters related to the modulus of the decay amplitudes and ${\bar{\alpha}}$ and $\Delta$ are relative phases. $A^{+-}$ is chosen to be a real positive quantity, which sets the phase convention. The weak phase $\alpha$ clearly appears as $2\alpha=\arg({\bar{A}}^{+0}/A^{+0})$. The parameter $\bar{\alpha}$, satisfying $\bar{\alpha}=2\arg({\bar{A}}^{+-}/A^{+-})$, would coincide with $\alpha$ in the limit of a vanishing penguin contribution.444This model is a valid representation of the amplitude system for $\alpha\neq 0$. When $\alpha$ vanishes exactly, the constraints $\bar{a}=a$ and $\bar{\alpha}=0$ must be added to the system in order to satisfy the equality of the $CP$-conjugate amplitudes $A^{+-}={\bar{A}^{+-}}$.

This alternative representation provides a convenient and compact parametrisation enabling a fast and stable exploration of the six-dimensional parameter space. Besides its technical convenience, this representation has a pedagogical benefit, as it exhibits the discrete ambiguities affecting the determination of the weak phase $\alpha$ in the $B\rightarrow hh$ modes in a clear way. Indeed, the above experimental $B\rightarrow\pi\pi$ observables can be written as a function of the chosen set of parameters as

[TABLE]

where we define $c=\cos(\alpha-\Delta)$ and ${\bar{c}}=\cos(\alpha+\Delta-2{\bar{\alpha}})$. The system exhibits a fourfold trigonometric ambiguity under the phase redefinitions:

[TABLE]

There is also an additional discrete symmetry involving not only $\alpha$ and $\Delta$, but also $\bar{\alpha}$: indeed, $c$, $\bar{c}$ and ${\cal S}^{+-}$ are left invariant by the reflection:

[TABLE]

The combination of these two discrete symmetries yields an heighfold ambiguity in the determination of the $\alpha$ angle in the $[0,180]^{\circ}$ range. Geometrically speaking, the fourfold ambiguity results from the choice left concerning the position of the apex with respect to the $A^{+0}$ base for each of the two isospin triangles, and a twofold ambiguity arises in relation with the relative direction of the two $A^{+0}$ and $\bar{A}^{+0}$ bases.

In the absence of any additional input, the available $B\rightarrow\pi\pi$ observables lead thus to height strictly equivalent mirror solutions for $\alpha$. The degeneracy could be partially lifted with the measurement of the time-dependent $CP$ asymmetry in the $B^{0}\rightarrow\pi^{0}\pi^{0}$ decay, $S^{00}$. This observable can be written as:

[TABLE]

where $s=\sin(\alpha+\Delta)$ and ${\bar{s}}=\sin(\alpha-\Delta+2{\bar{\alpha}})$ are not invariant under the phase redefinition of $\alpha$ and $\Delta$ Eq. (30), leaving only a twofold ambiguity on $\alpha$. The $B^{0}\rightarrow\pi^{0}\pi^{0}$ decay is so far observed only in the 4-photon final state preventing the measurement of the time-dependent decay rate. Future high-luminosity facilities have investigated the feasibility of the ${\cal S}^{00}_{\pi\pi}$ asymmetry measurement, either by considering the rare Dalitz decays of the neutral pions or by exploiting the conversion of photons in the detector material [38, 39].

The consistency of the experimental $B\rightarrow\pi\pi$ data with the triangular isospin relation of Eq. (20) can be assessed by the means of the two-side sum:

[TABLE]

where $a^{ij}$ denotes the normalised amplitude $a^{ij}={A^{ij}}/{A^{+0}}$. The sum of the length of two sides of a triangle being greater than the third one [40], we have $t\geq 1$, the equality occurring for a flat triangle. The values

[TABLE]

are obtained for the $B\rightarrow\pi\pi$ and ${\bar{B}}\rightarrow\pi\pi$ systems and are consistent with an almost-flat and an open triangle, respectively, as illustrated in Fig. 4. Both isospin triangles display the expected mirror symmetry, therefore the $\alpha$ constraint exhibits height non-degenerate solutions in $[0,180]^{\circ}$, as illustrated in Fig. 5. The minima of the $\chi^{2}$ test statistic over the parameter space ($\chi^{2}_{\rm min}=0$) are found at $\alpha=5.9^{\circ}$, $\alpha=84.1^{\circ}$, $\alpha=100.1^{\circ}$, $\alpha=124.4^{\circ}$, $\alpha=129.6^{\circ}$, $\alpha=140.4^{\circ}$, $\alpha=125.6^{\circ}$ and $\alpha=169.9^{\circ}$. This system of six observables for six independent parameters is not over-constrained and the vanishing value of $\chi^{2}_{\rm min}$ reflects the closure of both the $\bar{B}\rightarrow\pi\pi$ and the $B\rightarrow\pi\pi$ isospin triangles. The corresponding (symmetrised) $\alpha$ intervals at 68% and 95% CL are

[TABLE]

The interval $\alpha=[20.0,71.0]^{\circ}$ is excluded at more than 3 standard deviations by the ${SU}\mkern-2.0mu(2)$ isospin analysis of the $B\rightarrow\pi\pi$ system. The solution around $90^{\circ}$ is in excellent agreement with the indirect determination of $\alpha$, Eq. (2), from the other constraints of the global CKM fit [18]. If we include $\alpha_{\rm ind}$ derived from Fig. 3 as an additional and independent constraint to the ${SU}\mkern-2.0mu(2)$ isospin fit, the minimal $\chi^{2}$ increases by 0.3, indicating a consistency at the level of 0.6 standard deviation.

It is interesting to note the isospin constraints present a discontinuity at $\alpha=0$ (modulo $\pi$) as illustrated on the right-hand side of Fig. 5. The $\alpha=0$ hypothesis is rejected with a high significance ($\chi^{2}(\alpha=0)=296$) as expected by the direct $CP$ violation observed in the $B^{0}\rightarrow\pi^{+}\pi^{-}$ decay [31, 32]. However, for any finite $\alpha$ in the vicinity of zero, the isospin relation can be satisfied if arbitrary large penguin amplitudes are allowed (the limit case $\alpha\rightarrow 0$ happens if $|T^{+-}|$, $|T^{00}|$, $|P|$ are sent to infinity, see Eq. (27) and Refs. [36, 42]). The limit $\alpha\rightarrow 0$ was also discussed in the context of the Bayesian statistical approach in Ref. [41]. External data, e.g., based on ${SU}\mkern-2.0mu(3)$ consideration, could be used to set bounds on the penguin over tree ratio and thus further constrain $\alpha$ around 0. These bounds stand beyond the pure ${SU}\mkern-2.0mu(2)$ approach adopted here and we will not attempt at including such additional information for the following reasons. On one hand, other charmless modes (such as $\rho\pi$) will provide further information on the small $\alpha$ region, and on the other hand, we will ultimately use the determination of $\alpha$ in a global fit together with other well-controlled constraints on the CKM parameters [18]: both types of constraints will rule out the region around $\alpha=0$, so that additional theoretical hypotheses on the size of the $B\rightarrow\pi\pi$ amplitudes are not necessary for our purposes.

3.3 Isospin analysis of the $B\rightarrow\rho\rho$ system

A similar ${SU}\mkern-2.0mu(2)$ isospin analysis can be performed on each polarisation state of the $B\rightarrow\rho\rho$ system. Experimentally, the decays are completely dominated by the longitudinal polarisation of the $\rho$ mesons, and in the following, we will always imply this polarisation for the $\rho\rho$ final state (which will be denoted explicitly as $\rho_{L}\rho_{L}$ when we quote experimental measurements). This final state presents several advantages with respect to $B\rightarrow\pi\pi$. The measured branching fractions ${\cal B}^{+-}$ and ${\cal B}^{+0}$ are approximately five times larger than for $B\rightarrow\pi\pi$, while ${\cal B}^{00}$ is of same order of magnitude, indicating a smaller penguin contamination in the $\rho\rho$ case. Moreover, the time-dependent $CP$ asymmetry in the $B^{0}\rightarrow\rho^{0}\rho^{0}$ decay, ${\cal S}^{00}$, is experimentally accessible for the final state with four charged pions, potentially lifting some of the discrete ambiguities affecting the determination of $\alpha$. However, the current measurement ${\cal S}^{00}=0.3\pm 0.7$ suffers from large uncertainties, leaving pseudo-mirror solutions in $\alpha$. The available experimental observables555As can be seen from Tab. 4, the longitudinal polarisation $f_{L}$ is always used as an input in combination with branching ratios ${\cal B}$, and both are defined as $CP$-averaged quantities. and their current world averages are summarised in Tab. 4. Under the ${SU}\mkern-2.0mu(2)$ isospin hypothesis, the direct $CP$ asymmetry in $B^{+}\rightarrow\rho^{+}\rho^{0}$ vanishes and we will not take into account the experimental measurement of this quantity (which is consistent with our hypothesis and will be used to test this assumption in Sec. 4.1). Seven independent observables are available for the longitudinal helicity state of the $\rho$ mesons, allowing us to over-constrain the six-dimensional parameter space of the ${SU}\mkern-2.0mu(2)$ isospin analysis.

The situation of the $B\rightarrow\rho\rho$ system is illustrated in Fig. 6. The following two-side sums are consistent with flat triangles for both $B\rightarrow\rho\rho$ and ${\bar{B}}\rightarrow\rho\rho$ systems:

[TABLE]

This further reduces the degeneracy of the (pseudo) mirror solutions in $\alpha$ leaving a twofold ambiguity; see Fig. 7. The minima of the $\chi^{2}$ test statistic over the parameter space ($\chi^{2}_{\rm min}=0.14$) are found at $\alpha=92.1^{\circ}$ and $178.0^{\circ}$. The corresponding 68 and 95% CL intervals on $\alpha$ are

[TABLE]

The solution around $90^{\circ}$ is in excellent agreement (less than 0.1 standard deviation) with the indirect determination of $\alpha$, Eq. (2), from the other constraints of the global CKM fit [18], and is more tightly constrained than in the $B\rightarrow\pi\pi$ case, in relation with the smaller penguin contamination.

The region $\alpha=[14.0,76.0]^{\circ}\cup[112.0,158.0]^{\circ}$ is excluded at more than 3 standard deviations by the ${SU}\mkern-2.0mu(2)$ isospin analysis of the $B\rightarrow\rho\rho$ system. On the right hand side of Fig. 7, a zoom around $\alpha=0$ exhibits a small discontinuity, corresponding to $\chi^{2}(\alpha=0)=1.61$ and indicating that this hypothesis is mildly disfavoured by the data, consistently with the absence of large direct CP asymmetries in $B\rightarrow\rho\rho$ decays.

3.4 Isospin analysis of the $B\rightarrow\rho\pi$ system

The three neutral and two charged $B^{0,+}\rightarrow(\rho\pi)^{0,+}$ decays can be described with 10 complex (tree and penguin) amplitudes and one weak phase $\alpha$, i.e., 21 real parameters. Assuming the pentagonal isospin relation Eq. (13) that leaves only two independent complex penguin contributions, the number of degrees of freedom of the system is reduced to 12 after setting an irrelevant global phase to zero. The dimension of the system gets down to 10 if we consider the three neutral decays only. The time-dependent Dalitz analysis of the neutral $B^{0}\rightarrow(\rho^{\pm}\pi^{\mp},\rho^{0}\pi^{0})\rightarrow\pi^{+}\pi^{-}\pi^{0}$ transitions has been shown to carry enough information to fully constrain the isospin-related $B\rightarrow\rho\pi$ system, thanks to the finite width of the intermediate $\rho$-mesons that yields a richer interference pattern of the three-body decay [15, 16]. For instance, 11 independent phenomenological observables can be defined from the flavour-tagged time-dependent Dalitz distribution of the three-pion decay, e.g., the branching fractions of the three intermediate $\rho^{i}\pi^{j}$ transitions, the corresponding direct and time-dependent $CP$ asymmetries and two relative phases between their amplitudes.

The extraction of the parameter $\alpha$ through the ${SU}\mkern-2.0mu(2)$ analysis of the neutral modes will be referred to as the “Dalitz analysis” and will be considered first in Sec. 3.4.1. An extended analysis (referred to as the “pentagonal analysis”), including the information coming from the charged decays $B^{+}\rightarrow\rho^{+}\pi^{0}$ and $B^{+}\rightarrow\rho^{0}\pi^{+}$, will be discussed in Sec. 3.4.2.

3.4.1 Dalitz analysis of the neutral $B^{0}\rightarrow(\rho\pi)^{0}$ modes

Considering the neutral modes only, the time-independent amplitude of the $\pi^{+}\pi^{-}\pi^{0}$ final state can be written as

[TABLE]

where $A^{i}$ (${\bar{A}}^{i}$) is the amplitude of the $B^{0}({\bar{B}}^{0})\rightarrow\rho^{i}\pi^{j}$ transition666With this convention the $A^{i}$ (${\bar{A}}^{i}$) amplitude carries a superscript referring to the electric charge of the $\rho^{i}$ meson in both $B^{0}$ and ${\bar{B}}^{0}$ decays, i.e., $A^{+}=A^{+-}$, $A^{-}=A^{-+}$ and ${\bar{A}}^{+}={\bar{A}}^{-+}$, ${\bar{A}}^{-}={\bar{A}}^{+-}$ where the $A^{ij}$ (${\bar{A}}^{ij}$) amplitude is defined in Sec. 2.1. and $f_{i}$ ($i=-,+,0$) is the form factor accounting for the $\rho^{i}$ line-shape. Neglecting the tiny $B^{0}$ width difference $\Delta\Gamma_{d}$, the time-dependent decay rate can be written as a function of three combinations of these amplitudes:

[TABLE]

with $\Gamma_{0}=(|A_{3\pi}|^{2}+|{\bar{A}}_{3\pi}|^{2})$, $\Gamma_{\cal C}=(|A_{3\pi}|^{2}-|{\bar{A}}_{3\pi}|^{2})$ and $\Gamma_{\cal S}={\cal I}m\left(\frac{q}{p}{\bar{A}}_{3\pi}A_{3\pi}^{*}\right)$.

These amplitudes are described phenomenologically in terms of form factors $f_{i}$ describing the decay of the $\rho$ meson into two pions, multiplied by coefficients, denoted $\cal U$ and $\cal I$, related to the $B\rightarrow\rho\pi$ dynamics:

[TABLE]

with

[TABLE]

For given form factors $f_{i}$, these 27 coefficients $\cal U$ and $\cal I$ fully parametrise the time-dependent Dalitz distribution of the $B\rightarrow\pi^{+}\pi^{-}\pi^{0}$ decay. The ${\cal U}$ parameters are related to the branching fractions and the direct $CP$ asymmetries in the $B\rightarrow\rho^{i}\pi^{j}$ intermediate states while the ${\cal I}$ coefficients parametrise the mixing-induced $CP$ asymmetries.

As long as only the neutral modes are considered, it is worth noticing that the absolute normalisation of these coefficients is irrelevant for the ${SU}\mkern-2.0mu(2)$ isospin analysis. The measured values and correlations [50, 51] of the 26 remaining $\cal U$ and $\cal I$ coefficients (normalised with respect to ${\cal U^{+}_{+}}$ that is set to unity) are collected in Tabs. 5, 6, and 7. The following compact representation of the amplitude system for the neutral transitions can be chosen:

[TABLE]

This modelisation explicitly embeds the isospin relations Eq. (13) as well as the normalisation choice ${\cal U^{+}_{+}}=|A^{+}|^{2}+|{\bar{A}}^{+}|^{2}=1$. We thus obtain a system depending on 9 real parameters, consisting in four relative amplitudes ($\cos(\Theta^{+})$, $\mu^{-}$, $\mu^{0}$, $\nu^{-}$) and 5 phases ($\alpha$, $\phi^{-}$, $\phi^{0}$, $\Psi^{+}$, $\Psi^{-}$) (including the weak phase $\alpha$). It is worth stressing again that the actual parametrisation of the amplitudes has no impact on the results obtained in our frequentist framework [36]. The parametrisation Eq. (3.4.1) is mostly based on technical considerations: this minimal representation ensures a fast and stable exploration of the nine-dimensional parameter space constrained by the 26 correlated observables. The isospin relations for both $CP$-conjugate modes are illustrated on Fig. 8. The resulting constraint on $\alpha$ is shown on the left-hand side of Fig. 9.

Using the 26 correlated $\cal U$ and $\cal I$ observables, the minimum of the $\chi^{2}$ test statistic over the parameter spaces ($\chi^{2}_{\rm min}=8.6$) is found at $\alpha=54.1^{\circ}$. The corresponding 68% and 95% CL intervals are

[TABLE]

The interval $\alpha=[92.0,112.0]^{\circ}$ is excluded at more than 3 standard deviations by the ${SU}\mkern-2.0mu(2)$ isospin analysis of the $B^{0}\rightarrow(\rho\pi)^{0}$ system. These results do not show a good agreement with the indirect determination of $\alpha$, Eq. (2), based on the other constraints of the global CKM fit [18]. Adding the indirect determination as an additional constraint to the ${SU}\mkern-2.0mu(2)$ isospin fit, the minimal $\chi^{2}$ increases by 9.1 units, corresponding in the Gaussian limit to a 3.0 $\sigma$ discrepancy. In front of this uncomfortable discrepancy, we should remember that this result will be combined with constraints from $\rho\rho$ and $\pi\pi$ systems. Without a dedicated statistical study, it is difficult to determine the effective number of degrees of freedom associated with the set of 26 interdependent observables for the $\rho\pi$ system, but it remains interesting to notice that the $\chi^{2}_{\rm min}$ for $\rho\pi$ appears less significant (statistically) than in the case of $\rho\rho$ (which dominates the combination) and $\pi\pi$ (which contributes significantly). The latter constraints will dominate, leading to a combined determination remaining in good agreement with the indirect value of $\alpha$ obtained from the other constraints of the global fit. A detailed study of the actual significance of this $\chi^{2}$ offset in the $\rho\pi$ case is discussed in Sec. 4.2. In addition, the $\rho\pi$ Dalitz observables have been scrutinized to understand more precisely the origin of the tension, as reported in section 6.3.

3.4.2 “Pentagonal” analysis of the $B\rightarrow\rho\pi$ system

The amplitudes of the two charged decay modes $B^{+}\rightarrow\rho^{+}\pi^{0}$ and $B^{+}\rightarrow\rho^{0}\pi^{+}$ are related to the amplitudes of the three neutral decays $B^{0}\rightarrow(\rho\pi)^{0}$ through the pentagonal relation Eq. (13). The measurement of the branching ratios and ${CP}$-asymmetries of the charged modes may provide additional constraints to the isospin system. Considering simultaneously the charged and the neutral modes, the $\cal U$ and $\cal I$ observables that describe the relative amplitudes of the neutral decays must be supplemented with an absolute normalisation. This is achieved by identifying the sum of branching ratios, ${\cal B}^{+-}_{\rho\pi}+{\cal B}^{-+}_{\rho\pi}$, with the scaled amplitude $\mu^{+}({\cal U}^{+}_{+}+{\cal U}^{+}_{-}){\tau_{B^{0}}}/{2}$, where $\mu^{+}$ is the absolute normalisation of the relative $\cal U$ and $\cal I$ coefficients. The experimental measurements of the branching fractions for the charged and neutral $B\rightarrow\rho\pi$ modes and the ${CP}$ asymmetries for the charged modes are listed in Tab. 8.

At this stage, we would like to comment on some of these inputs. The measurement of the branching fractions for the neutral intermediate states $B^{0}\rightarrow\rho^{i}\pi^{j}$ is experimentally tricky. Indeed, the three $\rho$ mesons, charged and neutral, contribute to the neutral modes, whereas only one $\rho$ state contributes to each of the charged modes. A quasi-two-body analysis was performed by the BaBar and CLEO experiments [54, 52, 53] assuming that the interferences are negligible in the Dalitz region where the measurements were performed. The Belle measurement [55] extracts the partial branching ratio ${\cal B}^{\pm\mp}_{\rho\pi}$ and ${\cal B}^{00}_{\rho\pi}$ simultaneously from a global $B^{0}\rightarrow\pi^{+}\pi^{-}\pi^{0}$ Dalitz analysis, rescaling the ${\cal U}$ coefficients to the overall three-body branching fractions. The reliability of the measurements on the quasi-two-body intermediate states, $B^{0}\rightarrow\rho^{\pm}\pi^{\mp}$ and $B^{0}\rightarrow\rho^{0}\pi^{0}$, is thus questionable and this may bias the pentagonal constraint. Since the Dalitz-based ${SU}\mkern-2.0mu(2)$ isospin analysis already provides a complete and self-consistent description of the $B^{0}\rightarrow(\rho\pi)^{0}$ amplitude system, we will use the Dalitz analysis rather than the pentagonal one in the combined direct determination of the $\alpha$ angle. However, for completeness, we briefly discuss the results of the pentagonal analysis, interpreting the ${\cal B}^{\pm\mp}_{\rho\pi}$ measurement as an independent constraint on ${\cal B}^{+-}_{\rho\pi}+{\cal B}^{-+}_{\rho\pi}$. Further discussions can be found in Sec. 6.3.3, where predictions for these observables can be found.

Including the charged decays, the amplitude model for the neutral modes discussed in Sec. 3.4.1 can be completed with two additional real parameters ($\mu^{+0}$,$\phi^{+0}$):

[TABLE]

where the second identity uses the pentagonal relation Eq. (13) explicitly. Using Eq. (17), the $CP$-conjugate amplitudes are defined as

[TABLE]

Adding the normalisation factor $\mu^{+0}$, the whole $B\rightarrow\rho\pi$ dynamics is described through a 12-parameter system that is over-constrained by the set of observables for neutral and charged modes given in Tabs. 5, 6, 7 and 8.

The minimum of the $\chi^{2}$ test statistic over the parameter space ($\chi^{2}_{\rm min}=12.1$) is found at $\alpha=49.8^{\circ}$. The following 68% and 95% CL intervals on $\alpha$ are obtained:

[TABLE]

The compatibility with the indirect $\alpha$ determination is estimated at 3.2 standard deviations, slightly larger than for the Dalitz analysis of the neutral modes. The constraint on $\alpha$ is displayed on the right-hand side of Fig. 9. The same comments as in Sec. 3.4.1 apply also here concerning the relative statistical significance of this constraint compared to $B\rightarrow\pi\pi$ and $B\rightarrow\rho\rho$ channels.

3.5 Combined result

The ${SU}\mkern-2.0mu(2)$ isospin analyses of $B\rightarrow\pi\pi$, $B\rightarrow\rho\rho$ and $B\rightarrow\rho\pi$ provide three independent constraints on $\alpha$. A combined analysis can be performed by summing the individual $\chi^{2}(\alpha)$ curves. Using the Dalitz analysis for the $B^{0}\rightarrow(\rho\pi)^{0}$ mode together with $B\rightarrow\pi\pi$ and $B\rightarrow\rho\rho$, the minimum of the $\chi^{2}$ test statistic ($\chi^{2}_{\rm min}=17.1$) is obtained for $\alpha=86.2^{\circ}$, as illustrated on Fig. 10. A slightly disfavoured second solution peaks at $\alpha=178.4^{\circ}$ with $\Delta\chi^{2}=\chi^{2}(\alpha=178.4^{\circ})-\chi^{2}_{\rm min}=0.4$. The corresponding 68% and 95% CL intervals are

[TABLE]

The preferred solution near $90^{\circ}$ is consistent with the indirect determination $\alpha_{\rm ind.}$ given in Eq. (2) at 1.3 standard deviations. As discussed earlier, the combined constraint is dominated by $B\rightarrow\rho\rho$ and to a lesser extent, by $B\rightarrow\pi\pi$, whereas $B\rightarrow\rho\pi$ plays only a limited role. The constraint resulting of the partial combination based on $B\rightarrow\pi\pi$ and $B\rightarrow\rho\rho$ only is reported in Tab. 11.

The one-dimensional constraint on $\alpha$ can be recast in a constraint on the $(\bar{\rho},\bar{\eta})$ Wolfenstein parameters of the CKM matrix representing the apex of the $B$-meson unitarity triangle [10, 18]. The following relation can be derived:

[TABLE]

so that the curves at fixed $\alpha$ consist in circles centred on the point $(\bar{\rho},\bar{\eta})=(1/2,{\rm cotan}(\alpha)/2)$. In particular, the curve at $\alpha$=$90^{\circ}$ is a circle of radius $1/2$ centered on $(\bar{\rho},\bar{\eta})=(1/2,0)$, while $\alpha$=$0^{\circ}$ amounts to a circle of an infinite radius tangent to the axis $\bar{\eta}=0$. All the curves of constant $\alpha$ meet at the points $(0,0)$ and $(1,0)$. Fig. 11 displays the constraint resulting from the direct determination of $\alpha$ in the $(\bar{\rho},\bar{\eta})$ plane.

In this section, we have used world averages for all the observables, combining information from the BaBar, Belle and LHCb experiments. Since both $B$-factory experiments, BaBar and Belle, have measured all these observables independently, it is also possible to perform separate ${SU}\mkern-2.0mu(2)$ isospin analyses for each of the three decay channels and each experiment. The corresponding constraints are discussed in App. A.

4 Additional uncertainties in the extraction of $\alpha$

In this section, we are going to test the limits of the assumptions made to extract $\alpha$ from the data: the breakdown of isospin symmetry and the statistical approach used to build confidence intervals.

4.1 Testing the ${SU}\mkern-2.0mu(2)$ isospin framework

As discussed in Sec. 2, the extraction of the weak phase $\alpha$ relies on isospin symmetry, which is used at different levels [60, 61]:

• •

The charges of the $u$ and $d$ quarks are taken as identical. The \Delta I=\raisebox{1.72218pt}{\footnotesize3}\negthinspace/\negthinspace\raisebox{-2.58334pt}{\footnotesize2} contribution induced by the electroweak penguins topology to the $B^{+}\rightarrow h^{+}h^{0}$ decay is considered negligible. As the gluonic penguins only yield a \Delta I=\raisebox{1.72218pt}{\footnotesize1}\negthinspace/\negthinspace\raisebox{-2.58334pt}{\footnotesize2} isospin contribution, the pure \Delta I=\raisebox{1.72218pt}{\footnotesize3}\negthinspace/\negthinspace\raisebox{-2.58334pt}{\footnotesize2} $B^{+}\rightarrow h^{+}h^{0}$ decay only receives tree contributions in the absence of electroweak penguins. In this limit, the weak phase $\alpha$ can be identified as half the phase difference between the amplitude of the charged mode and its $CP$ conjugate.

• •

The masses of the $u$ and $d$ quarks are taken as identical. Isospin symmetry is assumed to be exact in the strong hadronisation process following the weak transition $\bar{b}(q)\rightarrow\bar{u}u\bar{d}(q)$ where $q$ represents the light spectator quark $u$ or $d$. This assumption allows one to relate the decay amplitudes of the charged ($\bar{b}u$) meson to the decay amplitudes of its isospin-related neutral state $\bar{b}d$ according to the pentagonal (triangular) identity given in Eq. (13) (respectively Eq. (20)).

• •

The $\rho\rho$ final state is supposed to obey Bose–Einstein statistics, and thus the analysis for the $\pi\pi$ and the $\rho\rho$ systems are supposed to follow the same isospin decomposition (for a given $\rho$ polarisation). However, the $\rho$ mesons cannot be distinguished only in the limit of a vanishing width. Once the finite $\rho$ width is taken into account, additional terms (forbidden by Bose symmetry fin the limit $\Gamma_{\rho}=0$) must be taken into account in the isospin decomposition of the amplitudes.

All these hypotheses are valid a priori to a very good approximation, but the accuracy reached in the determination of $\alpha$ in Sec. 3.5 is an incentive to investigate their limitations in more detail.

4.1.1 \Delta I=\raisebox{1.72218pt}{\footnotesize3}\negthinspace/\negthinspace\raisebox{-2.58334pt}{\footnotesize2} electroweak penguins

While preserving the isospin relations between charged and neutral decay amplitudes, the electroweak penguin topology induces a \Delta I=\raisebox{1.72218pt}{\footnotesize3}\negthinspace/\negthinspace\raisebox{-2.58334pt}{\footnotesize2} contribution to the charged modes $B^{+}\rightarrow h^{+}h^{0}$. With such a contribution, the system of amplitudes in the $\pi\pi$ and $\rho\rho$ cases still obeys the triangular relation Eq. (20), but Eq. (27) has to be rewritten as

[TABLE]

where the amplitude $P^{+0}_{\textsc{ew}}$ accounts for the \Delta I=\raisebox{1.72218pt}{\footnotesize3}\negthinspace/\negthinspace\raisebox{-2.58334pt}{\footnotesize2} contribution from electroweak penguins to the charged mode while $P^{+-}$ is redefined to absorb the contributions from both gluonic and electroweak penguins to the $B^{0}\rightarrow h^{+}h^{-}$ neutral decay.

The electroweak penguin contribution can be related to the tree amplitude in a model-independent way by performing the Fierz transformation of the relevant current-current operators in the effective Hamiltonian for $B\rightarrow\pi\pi$ decays [62, 63, 64]. This leads to an estimate of the relative contribution of the electroweak penguin, $P^{+0}_{\textsc{ew}}$, compared to the tree amplitude, $T^{+0}=(T^{00}+T^{+-})$, for the charged decay. Neglecting the penguins with internal light quarks $u$ and $c$ as well as the electroweak operators ${\cal O}_{7}$ and ${\cal O}_{8}$ (suppressed by tiny Wilson coefficients), the amplitude ratio:

[TABLE]

can be computed in terms of the short-distance Wilson coefficients ${\cal C}_{i}$ associated to the effective operators ${\cal O}_{i}$ describing the dominant electroweak penguin processes ($i=9,10$) and the current-current tree processes ($i=1,2$). Following this estimate, the electroweak penguin amplitude does exhibit no strong phase difference compared to the tree amplitude, preserving the charge symmetry $|A^{+0}|=|\bar{A}^{+0}|$ in the charged $B$ decay. Consequently, the impact of the electroweak penguin can be accounted for introducing the single real parameter $r_{\textsc{ewp}}={P_{\textsc{ew}}^{+0}}/{T^{+0}}$ in the isospin analysis of the $B\rightarrow hh$ systems.

The above modification of the amplitude $A^{+0}$ affects the determination of $\alpha$, which is the argument of $A^{+0}$ in the absence of \Delta I=\raisebox{1.72218pt}{\footnotesize3}\negthinspace/\negthinspace\raisebox{-2.58334pt}{\footnotesize2} penguin contributions. The CKM phase measured in $B\rightarrow hh$ systems under this hypothesis is an effective phase $\alpha_{0}$, which can be related to the true phase $\alpha$ through the relation:

[TABLE]

and the relative electroweak contribution $r_{\textsc{ewp}}$ induces the shift:

[TABLE]

on the determination of $\alpha$. A numerical evaluation yields [6, 62, 63, 64]:

[TABLE]

using the CKM matrix parameters constrained from the global analysis of the flavour constraints [18]. A maximal shift $\delta\alpha=-1.9^{\circ}$ occurs at $\alpha_{0}=90^{\circ}$, in agreement with the numerical results shown in Fig. 12. The resulting 68% CL intervals are

[TABLE]

Conversely, neglecting all other isospin-breaking effects, we can set a limit on the electroweak penguin contamination in the $B\rightarrow\pi\pi$ and $B\rightarrow\rho\rho$ systems. For this purpose, we require the modified isospin analysis to agree with the indirect prediction of $\alpha$, Eq. (2), determined from all the other constraints used in the global CKM fit [18]. Fig. 13 shows the resulting constraints on the amplitude ratio $r_{\textsc{ewp}}$ for $B\rightarrow\pi\pi$ (left panel) and $B\rightarrow\rho\rho$ (right panel) systems. We have rejected the mirror solutions that would yield very large, unphysical, differences $\delta\alpha$ and thus exceedingly large electroweak penguin contributions (see the discussion at the end of Sec. 3.2). This is achieved by applying an additional constraint on the strong and electroweak penguin hierarchy $|P_{\textsc{ew}}^{+0}|<|P^{+-}|$ for this specific study.777While eliminating mirror solutions, this constraint does not distort the $p$-value for $\alpha$ in the vicinity of the solution compatible with the indirect determination. We stress that we do not use such a constraint anywhere else in this article. The following 68% CL intervals are obtained:

[TABLE]

showing results consistent with an electroweak penguin contribution of only a few percents, in agreement with our theoretical expectations.

The contribution of the electroweak penguin can be added in a similar way to the $B\rightarrow\rho\pi$ system [21]:

[TABLE]

where $T^{+}$ and $P_{\textsc{ew}}^{+}$ denote the penguin and tree contributions to the sum of the charged amplitudes. Considering only the Dalitz analysis of the neutral modes, the penguin triangular relation Eq. (14) is violated:

[TABLE]

The weak phase $\alpha$ derived from

[TABLE]

can still be extracted by constraining the electroweak penguin-to-tree ratio:

[TABLE]

As discussed in Refs. [65, 66], the numerical expectation given in Eq. (53) does not apply to the $\rho\pi$ system: there are additional contributions in this case compared to $\pi\pi$ and $\rho\rho$, coming from matrix elements of tree operators that do not cancel anymore since the final state does not have a symmetric wave function. Following Refs. [65, 66], we parametrise the penguin-to-tree ratio for the $B\rightarrow\rho\pi$ transition as

[TABLE]

$r_{\textsc{ewp}}$ is the ratio given by Eq. (53) for the symmetrical $B\rightarrow hh$ transition where the Bose-Einstein symmetry applies. On the other hand, $r_{a}=|r_{a}|e^{\mathsf{i}\delta_{a}}$ is a correction parameter accounting for the non-vanishing $I_{f}=1$ contribution, with the notation $\langle(\rho\pi)^{+}|=\langle\rho^{+}\pi^{0}|+\langle\rho^{0}\pi^{+}|$. We let the correction parameter free to vary up to $|r_{a}|<0.3$: this rather conventional rule of thumb seems fairly conservative compared to estimates based on naive factorisation [65, 66], but data-based extractions in the case of $B\rightarrow K^{*}\pi$ decays suggest that it is the appropriate order of magnitude [67]. The unknown phase $\delta_{a}$, which can generate a charge asymmetry in $B^{+}\rightarrow(\rho\pi)^{+}$ decays, is unconstrained in the $B^{0}\rightarrow(\rho\pi)^{0}$ Dalitz analysis considered here.

As illustrated on the left panel of Fig. 14, including the estimate of \Delta I=\raisebox{1.72218pt}{\footnotesize3}\negthinspace/\negthinspace\raisebox{-2.58334pt}{\footnotesize2} electroweak penguin in the $B^{0}\rightarrow(\rho\pi)^{0}$ isospin analysis induces a small shift of $\delta\alpha=-1.2^{\circ}$ for the preferred value for $\alpha$:

[TABLE]

indicating that even a larger bound on $|r_{a}|$ would have only a limited impact.

The right panel of Fig. 14 shows the combination of the three channels in the presence of the \Delta I=\raisebox{1.72218pt}{\footnotesize3}\negthinspace/\negthinspace\raisebox{-2.58334pt}{\footnotesize2} electroweak penguin that produces an overall shift of $\delta\alpha=-0.7^{\circ}$ on the solution near $90^{\circ}$. The corresponding 68% CL interval is:

[TABLE]

in agreement with the indirect value $\alpha_{\rm ind}$ given in Eq. (2) at 1.5 standard deviations.

Assuming the same amplitudes ratio, $r_{\textsc{ewp}}$, hold for the three decay systems (modulo the correction term included for the asymmetrical $B\rightarrow\rho\pi$ decay), the combined 68% CL interval:

[TABLE]

is obtained when $\alpha$ is constrained to the indirect determination from Eq. (2) and any other ${SU}\mkern-2.0mu(2)$-breaking effect is neglected, in good agreement with the theoretical expectation Eq. (53).

4.1.2 Isospin-breaking effects due to mixing in the $\pi\pi$ system

We have worked up to now under the assumption that isospin symmetry was exact for the hadronic part of the $B$-meson decay. Even though isospin-breaking effects are known to be tiny, due to the very small mass difference between the $u$ and $d$ quarks, it is interesting to assess more precisely how it could affect our analysis.

In the $B\rightarrow\pi\pi$ system, the breaking of isospin symmetry due to the difference of quark masses triggers a mixing between light pseudoscalar mesons. This affects in particular the $\pi^{0}$ meson. Following Refs. [60, 61], at leading order of isospin breaking, the $\pi^{0}$ state can be written as

[TABLE]

where the flavour states are defined as follows: $|\pi_{3}\rangle$ represents the $I=0$ component of the ${SU}\mkern-2.0mu(2)$ triplet and $|\eta\rangle$, $|\eta^{\prime}\rangle$ are the states resulting from the mixing of the ${SU}\mkern-2.0mu(3)$ pseudoscalar meson octet and singlet components (according to the flavour decomposition $3\otimes 3=8\oplus 1$). The triangular isospin relation that applies to the isospin amplitudes:

[TABLE]

can be translated into mass-eigenstate amplitudes upon introducing the shifts $\delta A^{i0}$:

[TABLE]

The mixing among light pseudoscalar mesons thus generates a slight modification of the triangular relation in the $B\rightarrow\pi\pi$ amplitudes:

[TABLE]

with a similar equation for the $CP$-conjugate amplitudes.

At leading order in $\epsilon_{\eta^{(^{\prime})}}$, the amplitude shifts $\delta A^{i0}$ can be written in terms of the $B\rightarrow\pi\eta^{(^{\prime})}$ amplitudes $A^{\pi\eta^{(^{\prime})}}$ as

[TABLE]

where the factor $\sqrt{2}$ in the first equation accounts for the Bose–Einstein symmetry in the symmetric $A^{00}-A^{33}$ amplitude. Both neutral and charged $B\rightarrow\pi\eta^{(^{\prime})}$ decay amplitudes can be constrained using the experimental branching ratios and ${\cal CP}$-asymmetries reported in Tab. 9. Contrary to the \Delta I=\raisebox{1.72218pt}{\footnotesize3}\negthinspace/\negthinspace\raisebox{-2.58334pt}{\footnotesize2} contribution discussed in Sec. 4.1.1, the a priori unknown strong phase affecting $\delta A^{+0}$ may generate a charge asymmetry $|A^{+0}|\neq|\bar{A}^{+0}|$ in the charged $B$ meson decay. Consequently, the amplitude system is further constrained by introducing the measured $CP$-asymmetry in the $B^{+}\rightarrow\pi^{+}\pi^{0}$ mode, see Tab. 10, which was not considered in the isospin-symmetric analysis.

We must add eight complex amplitudes, $A^{\pi^{+}\eta}$, $A^{\pi^{+}\eta^{\prime}}$, $A^{\pi^{0}\eta}$, $A^{\pi^{0}\eta^{\prime}}$ and their $CP$ conjugates, to the $B\rightarrow\pi\pi$ system. We use Ref. [69] for the numerical estimates of the mixing parameters:

[TABLE]

leading to a slightly modified constraint on the weak phase $\alpha$, shown in Fig. 15:

[TABLE]

at 68% CL. We obtain thus a global shift for $\alpha$ of $\delta\alpha=-0.5^{\circ}$. The linear increase of the 68% interval, namely $\pm 1.5^{\circ}$, is mostly due to the limited experimental resolution on the $B\rightarrow\pi\eta^{(^{\prime})}$ data and it must be considered as an upper limit of the mixing-induced breaking effect on $\alpha$.

It is possible to go one step further using the admittedly naive simplification that there is no gluonic component to the $\eta^{(^{\prime})}$ states. Additional relations between the $B\rightarrow\pi\pi$ and the $B\rightarrow\pi\eta^{(^{\prime})}$ amplitudes can then be derived, unaffected by first-order ${SU}\mkern-2.0mu(3)$ breaking [70]:

[TABLE]

where $\phi_{P}$ is the mixing angle for the isospin singlets $|\eta_{q}\rangle$ = ${(|u{\bar{u}}\rangle+|d{\bar{d}}\rangle)}/{\sqrt{2}}$ and $|\eta_{s}\rangle$=$|s{\bar{s}}\rangle$ forming the $\eta^{(^{\prime})}$ states. With these additional constraints, only four additional complex parameters are needed to account for the $B\rightarrow\pi\eta^{(^{\prime})}$ amplitudes. In that case, the deviation for the triangular isospin relation Eq. (67) is given by:

[TABLE]

Using this relation and without further assumption on the $\eta$-$\eta^{\prime}$ mixing angle, the isospin-breaking effect on the $B\rightarrow\pi\pi$ determination of $\alpha$ is found to be $\delta\alpha=-1^{\circ}$ with a slightly reduced resolution $\pm 1.0^{\circ}$. Assuming the physical $\eta^{(^{\prime})}$ states are well described by the flavour combinations:

[TABLE]

corresponding to an $\eta-\eta^{\prime}$ mixing angle $\tan(\phi_{P})=1/\sqrt{2}$, a deviation $\delta\alpha=(-0.3\pm 1.4)^{\circ}$ is obtained for the $B\rightarrow\pi\pi$ system, in the same ball park as our previous estimate.

4.1.3 Additional isospin-breaking effects for the $\rho\rho$ and $\rho\pi$ systems

In the limit of a vanishing $\rho$ meson width, the isospin analysis of each $B\rightarrow\rho\rho$ polarisation state is identical to the $B\rightarrow\pi\pi$ analysis, with a dominance of the longitudinal polarisation. In parallel with the previous section, isospin breaking manifests itself though the $\rho-\omega-\phi$ mixing:

[TABLE]

where the mixing term $\epsilon_{\omega}$ is of ${\cal O}(1\%)$ and $\epsilon_{\phi}$ is negligible as there is an almost ideal mixing in the case of the $\phi$ meson. Additional sources of isospin breaking can occur in the $\rho\rightarrow\pi\pi$ decay (see Ref. [21] for a detailed review):

• •

The isospin breaking due to the $\pi^{0}-\eta-\eta^{\prime}$ mixing may affect the $\rho^{+}\rightarrow\pi^{+}\pi^{0}$ decay but turns out negligible, as the leading term in $\epsilon_{\eta}$ is suppressed by the small $\rho^{+}\rightarrow\pi^{+}\eta$ decay rate: ${\Gamma(\rho^{+}\rightarrow\pi^{+}\eta)}/{\Gamma(\rho^{+}\rightarrow\pi^{+}\pi^{0})}<0.6\%$ [71].

• •

Differences in the di-pion couplings for the neutral and charged $\rho$ mesons are experimentally limited to less than 1%: indeed, $1-({\Gamma_{\rho^{+}}}/{\Gamma_{\rho^{0}}})=(0.2\pm 0.9)\%$ [71]

• •

Isospin breaking affecting the $\rho^{0},\omega\rightarrow\pi^{-}\pi^{+}$ interference is restricted to a small window in the $\pi\pi$ mass spectrum: this effect, integrated over the whole $\pi\pi$ range, is estimated at the order of 2% [21].

• •

The dominant source of isospin breaking is actually due to the large decay width $\Gamma_{\rho}$ that makes the two final-state mesons distinguishable in the $B\rightarrow\rho_{1}(\pi\pi)\rho_{2}(\pi\pi)$ decay. This results in a residual $I_{f}=1$ amplitude contribution, forbidden by Bose–Einstein symmetry in the limit $\Gamma_{\rho}=0$, but potentially as large as $({\Gamma_{\rho}}/{m_{\rho}})^{2}\sim 4\%$. It is a slight abuse of language to call this effect an isospin-breaking contribution, as it does not vanish in the limit $m_{u}=m_{d}$ (but it does in the limit of a vanishing decay width $\Gamma_{\rho}\sim(m_{(\pi\pi)_{1}}-m_{(\pi\pi)_{2}})=0$).

The above isospin-breaking effects can be accounted for by modifying the isospin triangular relation in the same way as in Eq. (67):

[TABLE]

The unknown amplitude $\Delta^{+0}$ is parametrised in a generic way:

[TABLE]

where $r_{T}$ and $r_{P}$ are two independent complex parameters accounting for isospin breaking in the tree and penguin contributions, respectively. The arbitrary strong phase affecting these parameters can generate $CP$ violation in the charged decays. The measured $CP$ asymmetry in the $B^{+}\rightarrow\rho^{+}\rho^{0}$ mode given in Tab. 10 is thus included as an additional constraint in the presence of these isospin-breaking terms. Fig. 16 illustrates the determination of the weak phase $\alpha$ in the case of isospin-breaking contributions smaller than $|r_{T}|,|r_{P}|<4\%$ and the very conservative bound $|r_{T}|,|r_{P}|<10\%$. A small shift on the preferred $\alpha$ value near $90^{\circ}$, $\delta\alpha=-0.6^{\circ}$ $(-1.2^{\circ})$, results from isospin-breaking contributions limited to 4% (10%) respectively. The linear increase of the 68% interval, $({}^{+0.2}_{-1.7})^{\circ}$ for $|r_{T}|,|r_{P}|<4\%$, has to be understood as an upper limit of the impact of these contributions on the determination of $\alpha$.

Isospin breaking arises similarly in the $B\rightarrow\rho\pi$ system through mixing in the light pseudoscalar and vector sectors. However, as long as the weak phase determination is limited to the neutral $B^{0}\rightarrow(\rho\pi)^{0}$ Dalitz analysis, the isospin constraints reduces to the triangular penguin relation Eq. (14). Any isospin breaking in the tree amplitudes can be absorbed by a redefinition of the unconstrained charged amplitudes $T^{+0}$, $T^{0+}$. Corrections to the penguin amplitude relation Eq. (14) in the $B\rightarrow\rho\pi$ system will be suppressed by the small penguin-to-tree ratios. For instance, the impact of the $\pi^{0}-\eta-\eta^{\prime}$ mixing on the weak phase determination in the $B^{0}\rightarrow(\rho\pi)^{0}$ analysis is suppressed by factors like $\epsilon_{\eta^{(^{\prime})}}\left|{P_{\rho\eta^{(^{\prime})}}}/{T^{+}}\right|$, where $P_{\rho\eta^{(^{\prime})}}$ represents the penguin contribution to the $B^{+}\rightarrow\rho^{+}\eta^{(^{\prime})}$ transition and $T^{+}$ is defined in Eq. (4.1.1). The corresponding deviation is found to be negligible, $\delta\alpha<0.1^{\circ}$ [21]. Therefore, apart for the electroweak penguin contribution discussed in section 4.1.1, no further isospin breaking will be considered for the $B^{0}\rightarrow(\rho\pi)^{0}$ analysis.

4.2 Impact of the statistical treatment

4.2.1 $p$-values based on the bootstrap approach and Wilks’ theorem

In Sec. 3.1, we have outlined our framework to determine confidence intervals on $\alpha$, starting from a test statistic $\Delta\chi^{2}(\alpha)$, which is converted into a $p$-value following Eq. (24). This procedure can be proven to be exact in the simple case where the observables obey Gaussian probability distribution functions and they are linearly related to the parameters of interest. Following Wilks’ theorem, this can be extended to more general cases, at least asymptotically, if the data sample is large and the observables resolutions are small enough to consider the problem as locally linear [23]. In such case, the construction in Sec. 3.1 ensures exact coverage: if one repeated the determination of $\alpha$ using independent data sets from many identical experiments, the 68% CL interval for $\alpha$ would encompass the true value of $\alpha$ in 68% of the cases.

However, many effects can alter this picture. Even in the exact Gaussian case, the $p$-value can be distorted when the observables have a nonlinear dependence on the fundamental parameters of interest. An example consists in Eq. (12), where the reexpression of the $CP$-asymmetries in terms of $\alpha_{\rm eff}$ is nonlinear and implements the trigonometric boundary on $CP$-asymmetries $\sqrt{{\cal C}^{2}+{\cal S}^{2}}<1$ [6]. More generally, one may wonder if we stand close to the hypotheses of Wilks’ theorem with the current set of data. Otherwise, the construction given in Sec. 3.1 might suffer from under- or over-coverage. In this section, we will assess the finite-size errors associated with our statistical framework by considering a different construction of the $p$-value that takes into account some of the effects deemed subleading in the Wilks-based approach.

We start by recalling some elements related to the construction of $p$-values in our context, as discussed in Refs. [22, 72, 73, 74]. We want to assess how much the data is compatible with the hypothesis that the true value of the weak phase $\alpha$, denoted $\alpha_{t}$, is equal to some fixed value $\alpha$, i.e., ${\mathcal{H}}_{\alpha}:\alpha_{t}=\alpha$. This hypothesis is composite, as it sets the value of some of the theoretical parameters, but not all of them. Indeed, the hadronic parameters (tree and penguin amplitudes) are also theoretical parameters required in our isospin analysis but are not set in ${\mathcal{H}}_{\alpha}$. These hadronic parameters are nuisance parameters, denoted collectively as $\vec{\mu}$.

The test statistic $\Delta\chi^{2}(\alpha)$ defined in Eq. (23) can be seen as the maximal likelihood ratio comparing the most plausible configuration under ${\mathcal{H}}_{\alpha}$ with the most plausible one in general. It is a definite positive function chosen in a way that large values indicate that the data present evidence against ${\mathcal{H}}_{\alpha}$: in the following, we will state explicitly its dependence on the data $X$ by using the notation $\Delta\chi^{2}(X;\alpha)$. A $p$-value is built by calculating the probability to obtain a value for the test statistic at least as large as the one that was actually observed, assuming that the hypothesis ${\mathcal{H}}_{\alpha}$ is true:

[TABLE]

where the probability distribution function (PDF) $h$ of the test statistic is obtained from the PDF $g$ of the data as

[TABLE]

A small value of the $p$-value thus provides evidence against the hypothesis ${\mathcal{H}}_{\alpha}$. We notice that in general the $p$-value Eq. (78) exhibits a dependence on the nuisance parameters through the PDF $h$, even though the test statistic $\Delta\chi^{2}$ itself is independent of $\vec{\mu}$.

For linear models, in which the observables $X$ depend linearly on the parameter $\alpha_{t}$, $\Delta\chi^{2}(X;\alpha)$ is a sum of standard normal random variables, and is distributed as a $\chi^{2}$ with $N_{dof}=1$. Under the conditions of Wilks’ theorem [23], this property can be extended to non-Gaussian cases, the distribution of $\Delta\chi^{2}(X;\alpha)$ will converge to a $\chi^{2}$ law depending only on the number of parameters tested. The $p$-value becomes independent of the nuisance parameters $\vec{\mu}$ and can be still interpreted as coming from a $\chi^{2}$ law with $N_{dof}=1$. This is the rationale for the statistical framework presented in Sec. 3.1.

We compare this Wilks-based approach with the bootstrap one based on the plug-in principle [75]. It consists here in defining

[TABLE]

We perform the evaluation of the $p$-value for a given value of $\alpha$ by setting the nuisance parameters to the value minimising $\chi^{2}(\alpha,\vec{\mu})$. This approach assesses the role played by nuisance parameters by replacing them with an estimator $\hat{\vec{\mu}}(\alpha)$ that depends on $\alpha$. Other constructions could have been considered, with more conservative statistical properties (supremum, constrained supremum, etc.) [76], but this would go beyond the scope of our study. We focus here on the bootstrap approach: it is relatively simple to implement, it exhibits good coverage properties in the examples considered here and it provides a first glimpse of the role played by nuisance parameters in coverage, which is neglected in the Wilks-based approach used up to now.

4.2.2 Comparison of $p$-values for the extraction of $\alpha$

The comparison between the two approaches may be performed by examining their coverage properties. Once computed, a $p$-value can be used to determine confidence-level intervals for the parameter of interest $\alpha$. These intervals have a correct frequentist interpretation if the $p$-value exhibits exact coverage, i.e., for any $\beta$ between 0 and 1, $P(p\leq\beta)=\beta$. Over-coverage $P>\beta$ corresponding to a conservative $p$-value and too wide a CL interval, under-coverage $P<\beta$ corresponding to a liberal $p$-value and too narrow a CL interval. In general, an exact $p$-value, or if not possible, a reasonably conservative $p$-value, is desirable, at least for the confidence levels of interest. This conservative approach is generally adopted in high-energy physics in order to avoid rejecting a hypothesis (such as “the SM is true”) too hastily.

We have already discussed the extraction of the $p$-value using the Wilks-based approach. In the case of the bootstrap approach, the statistical coverage of the $\alpha$ intervals has been studied through a full frequentist exploration of the space of nuisance parameters. A complete analysis with toy Monte Carlo simulations was carried out in order to compute the PDF $h$ and thus the $p$-value for the angle $\alpha$. The individual constraints on $\alpha$ from the $B\rightarrow\pi\pi$, $B\rightarrow\rho\rho$ and $B^{0}\rightarrow\pi^{+}\pi^{-}\pi^{0}$ systems, are displayed in Fig. 17, comparing the bootstrap and Wilks-based approaches. Considering the 68% CL interval on $\alpha$, the Wilks-based approach is slightly more conservative than the bootstrap one for both $B\rightarrow\pi\pi$ and $B\rightarrow\rho\rho$ systems, whereas the situation is reversed for the $B^{0}\rightarrow(\rho\pi)^{0}$ analysis. The 68% CL intervals obtained with the bootstrap method for the three systems are

[TABLE]

For the latter system, the bootstrap approach slightly reduces the tension with the indirect determination $\alpha_{\rm ind}$ to 2.7 standard deviations.

As shown on Fig. 18, when combining the three analyses for the 68% CL interval, the over-coverage for $B\rightarrow\pi\pi$ and $B\rightarrow\rho\rho$ almost compensates the under-coverage for $B\rightarrow\rho\pi$. The 68% CL interval for the preferred solution close to $90^{\circ}$ increases by about 1∘, which can be considered as an estimate of the finite-size errors associated to the statistical framework. The second solution near $180^{\circ}$ is slightly less disfavoured, increasing the corresponding 68% CL interval by $\pm 1.4^{\circ}$ with respect to the Wilks-based approach. The overall combination results in

[TABLE]

This direct $\alpha$ measurement within the bootstrap approach is consistent with the indirect CKM determination $\alpha_{\rm ind}$ at the level of 1.1 standard deviation.

Going one step further, we can compare the coverage of the two methods. This can be done by studying the PDF of the $p$-values to check whether they are exact, conservative or liberal. Even though $p^{\mathrm{bootstrap}}$ does not depend on the true value of the nuisance parameters explicitly, its distribution does in general. For illustrative purposes, we choose here to compute the distribution of both tests assuming as true values $\alpha=92.5^{\circ}$ and $\vec{\mu}=\hat{\vec{\mu}}(\alpha=92.5^{\circ})$ (we do not attempt to investigate other values of $\alpha$ or the nuisance parameters). The computation of the PDF of $p^{\mathrm{bootstrap}}$ requires a twofold recursion of the bootstrap procedure (double bootstrap)888This computation is a CPU-consuming exercise which was carried on at the CC-IN2P3 computing farm (scoring $49\cdot 10^{3}\ \mathrm{HS06\cdot hours}$, i.e. approximatively 200 CPUs over one day).. The $p$-value for the $\Delta\chi^{2}$ test is obtained assuming that it obeys Wilks’ theorem while the bootstrap test is assumed to be a true $p$-value, i.e., uniformly distributed over $[0,1]$.

The results are shown in Fig. 19. Rather than showing the PDF of the two $p$-values as a function of $p$, we have expressed $p$ in units of $\sigma$ (assumed test significance) and expressed the cumulative distribution function (integral of the PDF from 0 to $p$) once again in units of $\sigma$ (true test significance). A $p$-value with exact coverage corresponds to a diagonal straight line, over-coverage (under-coverage) happens if the curve is above (below) this diagonal. In our particular case, one can clearly see the improved coverage from the bootstrap test with respect to the $\Delta\chi^{2}$ case, assuming Wilks’ theorem. The bootstrap distribution has a flat distribution at least up to the $2\sigma$ level, as shown by the curve in Fig. 19, which is close to a diagonal straight line (this is expected since we set the true values of the nuisance parameters to their plug-in values). On the other hand, the $\Delta\chi^{2}$ test combined with Wilks’ theorem is too aggressive (under-covers) for confidence levels above $0.6\sigma$ for the global combination of the three $B\rightarrow hh$ decays modes and above $1.2\sigma$ when considering $B\rightarrow\pi\pi$ and $B\rightarrow\rho\rho$ only. Due to the highly CPU-consuming nature of this computation, we have not tried to perform the same computation for other true values of $\alpha$ and the nuisance parameters, but one may expect that the results of these coverage tests should hold for values of $\alpha$ in the vicinity of the solution compatible with its indirect determination.

The difference in the 68% CL intervals between the bootstrap and Wilks-based approaches could be taken as a linear correction to the direct extraction of the parameter $\alpha$, i.e. focusing only on the solution around $90^{\circ}$:

[TABLE]

It is important to notice that the size of the bootstrap correction is an indication of the non-asymptotic regime only for the direct determination of the weak phase based on the current data on charmless decays. This estimate is likely to change when including additional or improved observables. In particular, the observed under-coverage for the Wilks-based approach is likely to be reduced in a global fit of the CKM matrix, where the direct $\alpha$ measurement Eq. (83) is combined with the more precise indirect determination Eq. (2). The bootstrap interval correction has been computed here for the direct measurement of $\alpha$ alone, for a particular set of observables. Therefore this correction cannot be considered as an absolute uncertainty to be included automatically when combined with another independent measurement.

4.3 Summary for the direct determination of $\alpha$

Our results for the direct determination of the weak phase $\alpha$ are summarized in Tab. 11 for the different model hypotheses and the different statistical approaches considered up to now. The 68% CL interval (only for the solution near $90^{\circ}$ 999For $B\rightarrow\pi\pi$, the bootstrap approach (second row in Tab. 11) provides two peaks that are just separated at 68% CL, as can be seen in Fig. 17. However, in Tab. 11, a single range for $\alpha$ is given for this method, which is obtained by merging the intervals corresponding to the two peaks. Even though this range is not the exact outcome of the bootstrap analysis for $B\rightarrow\pi\pi$ at 68% CL, it allows us to perform more meaningful comparisons with the extractions of $\alpha$ using the Wilks-based approach, for which the two peaks are not distinguished at this level of significance and a single range is obtained for $\alpha$.) is reported as well as the compatibility with the indirect $\alpha_{\rm ind}$ determination Eq. (2). We see that depending on the approach, the central value for the combination shifts by $2^{\circ}$ or less, remaining thus within the error quoted in Eq. (1). The uncertainty remains between $4^{\circ}$ and $5^{\circ}$ when isospin-breaking effects are allowed. On the other hand, as discussed in Sec. 4.2.2, the comparison between the bootstrap and Wilks-based approaches suggests a sizeable uncertainty attached to the statistical framework, but this uncertainty is attached to the direct extraction of $\alpha$ from the three $B\rightarrow\pi\pi$, $B\rightarrow\rho\rho$ and $B\rightarrow\rho\pi$ modes and cannot be used as such when it is combined with other constraints, e.g., within global fits. It is likely that one would get a smaller uncertainty if a similar analysis was performed with an extended set of observables leading to a more accurate determination of $\alpha$.

These various arguments lead us to keep Eq. (1) as our final answer for the direct extraction of $\alpha$ from charmless $B$ decays to be used in latter analyses.

5 Hadronic amplitudes

In addition to the CKM angle $\alpha$, our study of the $B\rightarrow\pi\pi$, $B\rightarrow\rho\pi$ and $B\rightarrow\rho\rho$ systems provides constraints on the hadronic amplitudes that cannot be computed in QCD directly. We can thus determine some features of hadronisation from the data, to be compared with the theoretical approaches proposed to describe these decays.

In order to sharpen the constraints on these parameters, we will not only consider the experimental measurements discussed previously, but we will also take as an additional constraint the indirect prediction $\alpha_{\rm ind}$ derived from Fig. 3 and Eq. (2). This constraint is obtained excluding $B\rightarrow\pi\pi$, $B\rightarrow\rho\rho$ and $B\rightarrow\rho\pi$ data from the global CKM fit, and therefore provides an independent constraint on the isospin analysis.

5.1 Penguin-to-tree ratios

Following the ${\cal C}$-convention defined in Sec. 2.1, the penguin-to-tree ratio is given by:

[TABLE]

The isospin analysis assumes that there are no penguin contributions to the charged decays in the $B\rightarrow\pi\pi$, $B\rightarrow\rho\rho$ and $B^{0}\rightarrow(\rho\pi)^{0}$ systems. This hypothesis allows us to extract $\alpha$ and the tree and penguin contributions simultaneously in each decay mode. We start with the (complex) penguin-to-tree ratio:

[TABLE]

in the $B^{0}\rightarrow(h_{1}^{i}h_{2}^{j})^{0}$ neutral modes. The isospin analysis of the $B\rightarrow\pi\pi$ system constrains the ratio of moduli $R^{ij}_{\pi\pi}$ for both colour-allowed ($B\rightarrow\pi^{+}\pi^{-}$) and colour-suppressed ($B^{0}\rightarrow\pi^{0}\pi^{0}$) neutral modes at a precision level of around $20\%$. The two-dimensional constraint in the (complex) $\tilde{R}^{ij}$ plane for these two modes is shown in Fig. 20. A large penguin contamination is observed in both modes. The one-dimensional extraction of the modulus and the phase gives the following 68% (respectively 95%) CL intervals:

[TABLE]

We observe that values of the phases close to 0 or $\pi$ are favoured, in agreement with QCD factorisation expectations [77, 78], which also predicts value for $R^{+-}_{\pi\pi}$ around 0.5 [79, 80]. Recent fits of two-body non-leptonic $B$ decays into light pseudoscalar mesons based on ${SU}\mkern-2.0mu(3)$ flavour symmetry yield smaller values for these ratios [81].

The penguin contamination to the colour-allowed $B^{0}\rightarrow\rho^{+}\rho^{-}$ decay is significantly smaller than for its $B\rightarrow\pi\pi$ counterpart, as expected from the flatness of the $\rho\rho$ isospin triangles discussed in Sec. 3.3. The two-dimensional constraint in the (complex) $\tilde{R}^{ij}_{\rho\rho}$ plane is shown in Fig.  21. A tight constraint on the penguin contamination to the $B^{0}\rightarrow\rho^{+}\rho^{-}$ decay is obtained with a 68% (respectively 95%) CL upper limit:

[TABLE]

on the edge of the expectations from QCD factorisation, which predicts this ratio to be around 0.12 [80]. Due to the small overall amplitudes contributing to the $B^{0}\rightarrow\rho^{0}\rho^{0}$ mode, a looser constraint is found:

[TABLE]

No constraint is obtained on the phase of the ratio for either $B^{0}\rightarrow(\rho\rho)^{0}$ decays.

Another representation of the penguin hierarchy between the $B^{0}\rightarrow\pi^{+}\pi^{-}$ and $B^{0}\rightarrow\rho^{+}\rho^{-}$ decays is provided by the phase shift $\Delta\alpha=(\alpha-\alpha_{\rm eff})$ where the effective weak phase $\alpha_{\rm eff}$, introduced in Eq. (12), describes the $B^{0}\rightarrow h^{+}h^{-}$ transition. In the limit of a vanishing penguin contribution, the mixing phase $\alpha_{\rm eff}$ coincides with the weak phase $\alpha$. The following 68% CL intervals on $\Delta\alpha$, derived from Fig. 22, are obtained

[TABLE]

The phase shift is consistent with the absence of penguin in the $B^{0}\rightarrow\rho^{+}\rho^{-}$ decay while the sizeable shift obtained in $B^{0}\rightarrow\pi\pi$ excludes the $\Delta\alpha=0$ hypothesis at 1.9 $\sigma$ in this mode, confirming the difference of size of penguin contributions for the two decay channels. A similar quantity can be defined for the asymmetric $B^{0}\rightarrow\rho^{\pm}\pi^{\mp}$ decay and is further discussed in Sec. 6.3.2.

Due to the limited experimental statistics currently available, the Dalitz isospin analysis of the three-pion decay does not lead to a significant constraint on the penguin contribution in any of the three $B^{0}\rightarrow(\rho\pi)^{0}$ decays. The pentagonal approach, involving further constraints from the charged modes, only provides a weak constraint on $R^{+-}$ and $R^{-+}$, see Fig. 23. Phases close to 0 or $\pi$ are mildly favoured, in agreement with the expectations from QCD factorisation. No constraint can be derived for the $B^{0}\rightarrow\rho^{0}\pi^{0}$ mode. For comparison, the constraints on the modulus of the penguin contribution to the various decays modes can be found in Fig. 24.

5.2 Testing colour suppression

The hadronisation of the ${\bar{b}}d\rightarrow{\bar{u}}d(u\bar{d})$ final state into a pair of charged isovector mesons, $h^{+}_{1}h^{-}_{2}$, is colour allowed while the hadronisation of this tree transition into a neutral pair, $h^{0}_{1}h^{0}_{2}$, is colour suppressed. This colour suppression can be probed through the ratio of the corresponding tree amplitudes101010We identify this ratio of “tree” amplitudes to a colour-suppressed ratio from the analysis of tree diagram topologies. In the $\cal C$-convention used here, this identification is correct only in the limit of vanishing $u$- and $c$-penguin topological amplitudes, as discussed in Sec. 2.1. Moreover, the $W$-exchange topology that provides a colour suppressed contribution that can be absorbed in the tree amplitudes for both $B\rightarrow h^{+}h^{-}$ and $B\rightarrow h^{0}h^{0}$ transitions are neglected here. A similar statement holds for the penguin-to-tree ratio discussed previously.:

[TABLE]

where $R_{\cal C}$ and $\delta_{\cal C}$ represent the modulus and the phase of the complex ratio, respectively.

The two-dimensional constraint on the phase and modulus of the colour-suppression ratio is shown in Fig. 25 for $B\rightarrow\pi\pi$, $B\rightarrow\rho\rho$ and $B\rightarrow\rho\pi$ channels. As expected, no constraint is obtained for that latter due to the limited statistics. The modulus of the colour-suppression ratio is determined with a relative resolution of about $10\%$ for the $B\rightarrow\pi\pi$ and $B\rightarrow\rho\rho$ modes. The corresponding one-dimensional 68% (95%) CL intervals on the modulus and the phase are:

[TABLE]

indicating that the colour suppression is much more effective for $B\rightarrow\rho\rho$ than for $B\rightarrow\pi\pi$. In studies based on QCD factorisation [78, 79, 80], similar values for these moduli are found, apart from the modulus of the ratio between colour-allowed and colour-suppressed tree amplitudes for $\pi\pi$ modes, $R_{\cal C}({\pi\pi})$, smaller and around 0.2. Similarly, the relative phase for the latter ratio is significantly different from 0 and $\pi$, contrary to general expectations from leading-order QCD factorisation. This may explain the difficulties experienced with this approach to reproduce the $B\rightarrow\pi^{0}\pi^{0}$ branching ratio [82] and may point towards significant power-suppressed contributions for some of these modes. Recent fits of two-body non-leptonic $B$ decays into light pseudoscalar mesons based on ${SU}\mkern-2.0mu(3)$ yield similar values for $R_{\cal C}({\pi\pi})$, indicating that large phase and modulus of the colour-suppressed tree contribution $T^{00}$ are required not only in the $\pi\pi$ system, but also in other $B\rightarrow PP$ decays [81].

6 Prediction of observables

If we use the indirect determination of the weak phase $\alpha_{\rm ind}$, Eq. (2), as an additional input, the system of isospin-related amplitudes becomes over-constrained for the three charmless decay modes of interest. The isospin analysis can then be used to perform an indirect determination for each of the experimental observables introduced in Sec. 3 and to make predictions for yet unmeasured quantities. For each observable ${\cal O}_{i}$, this indirect determination is obtained by performing a fit to the experimental data excluding the measurement of the observable considered. The compatibility of the indirect determination with the measurement is quantified by comparing the minimal $\chi^{2}$ values reached by the fit when excluding and including the direct measurement. The difference between these two quantities is distributed as a $\chi^{2}$ law with $N_{dof}=1$, whose interpretation in terms of standard deviations yields the compatibility pull defined as:

[TABLE]

where $\vec{}{\cal O}=\{\vec{}{\cal O}_{!i},{\cal O}_{i}\}$ represents the full data set of experimental measurements, and $\vec{}{\cal O}_{!i}$ the subset out of which the specific observable ${\cal O}_{i}$ has been excluded. The pull quantifies the compatibility of the measurement of the observable ${\cal O}_{i}$ with all the other experimental observables under the isospin hypothesis. Let us add that the individual pulls are not in general independent from each other, so that one should refrain from giving a statistical interpretation of their distribution in each of the $B\!\rightarrow\pi\pi$, $B\!\rightarrow\rho\pi$, $B\!\rightarrow\rho\rho$ systems.

6.1 $B\!\rightarrow\pi\pi$ observables

As already indicated in Sec. 3.2, the number of independent amplitudes in the $B\!\rightarrow\pi\pi$ system matches the number of experimental measurements available. Including the indirect determination of $\alpha$ as an additional constraint, the minimal $\chi^{2}$ increases only slightly, reflecting the excellent agreement between $\alpha_{\rm ind}$ and $B\!\rightarrow\pi\pi$ data, at the level of 0.6 standard deviation.

The indirect constraints on the three branching fractions $B^{ij}_{\pi\pi}$ are shown in Fig. 26 while Fig. 27 displays the two-dimensional constraint on the $CP$ asymmetries (${\cal C}_{\pi\pi},{\cal S}_{\pi\pi}$) for the neutral modes $B^{0}\!\rightarrow\pi^{+}\pi^{-}$ and $B^{0}\!\rightarrow\pi^{0}\pi^{0}$. When excluding any of the experimental measurements, either for the branching ratios ${\cal B}^{ij}_{\pi\pi}$ or for the direct $CP$ asymmetries ${\cal C}^{ij}_{\pi\pi}$, the $B\!\rightarrow\pi\pi$ amplitudes system is no longer over-constrained. As a consequence, the corresponding $\chi^{2}_{\rm min}(\alpha_{\rm ind},\vec{}{\cal O}_{!i})$ vanishes and the pull associated to these observables saturates its maximal value: ${\rm Pull}(B^{ij}_{\pi\pi},C^{ij}_{\pi\pi})=\sqrt{\chi^{2}_{\rm min}(\alpha_{\rm ind},\vec{}{\cal O})}$ which turns out to vanish thanks to the closure of both isospin triangles with the current world-average data. The yet unmeasured time-dependent asymmetry, ${\cal S}^{00}_{\pi\pi}$, is predicted to be:

[TABLE]

with a 68% confidence level, and larger than $0.33$ at 95% confidence level. The corresponding one-dimensional projection of the constraint is displayed in Fig. 28.

The experimental values, confidence intervals and pulls for the measured observables are reported in Tab. 14 in App. B.

6.2 $B\!\rightarrow\rho\rho$ observables

As shown in Sec. 3.3, the $B\!\rightarrow\rho\rho$ data is in very good agreement with the indirect determination of the weak phase $\alpha_{\rm ind}$. The ${SU}\mkern-2.0mu(2)$ isospin constraint on the ${\cal B}^{ij}_{\rho\rho}$ branching ratios and the corresponding fractions of longitudinal polarisation are displayed in Fig. 29. As expected theoretically [80], the isospin fit to the $B\!\rightarrow\rho\rho$ data favours large fractions of longitudinal polarisation in each of the three decay modes. The most precise constraint is obtained for the neutral $B^{0}\!\rightarrow\rho^{+}\rho^{-}$ decay with an indirect determination $f_{\rho_{L}\rho_{L};{\rm ind}}^{+-}>0.91$ (0.79) at 68% (95%) confidence level, in very good agreement with the experimental measurement $f_{\rho_{L}\rho_{L};{\rm dir}}^{+-}=0.990\pm 0.020$.

Weaker constraints are obtained on the $CP$ asymmetries (${\cal C}_{\rho\rho},{\cal S}_{\rho\rho}$) for both neutral modes $B^{0}\!\rightarrow\rho^{+}\rho^{-}$ and $B^{0}\!\rightarrow\rho^{0}\rho^{0}$ as shown in Fig. 30. Numerical values, confidence intervals and pulls for the measured observables are reported on Tab. 15 in App. B.

6.3 $B\!\rightarrow\rho\pi$ observables

6.3.1 $B^{0}\!\rightarrow\pi^{+}\pi^{-}\pi^{0}$ Dalitz analysis

The normalised ${\cal U}$ and ${\cal I}$ observables provide a complete description of the relative $B^{0}\!\rightarrow\pi^{+}\pi^{-}\pi^{0}$ decay amplitudes. As discussed in Sec. 3.4, the data from $B^{0}\!\rightarrow(\rho\pi)^{0}\rightarrow\pi^{+}\pi^{-}\pi^{0}$ Dalitz analyses disagrees with the indirect determination of $\alpha$ at the level of 3.0 standard deviations. This disagreement is reflected through the indirect determination of several of the form-factor coefficients ${\cal U}$ and ${\cal I}$ as listed in Tab. 16 in App. B, in particular ${\cal U}_{\textrm{\tiny+}}^{\textrm{\tiny-}}$, ${\cal U}_{\textrm{\tiny-}}^{\textrm{\tiny+}}$, ${\cal I}_{\textrm{\tiny-}}$, {\cal U}_{\textrm{\tiny+0}}^{\textrm{\tiny+\mathcal{R}e}}, {\cal U}_{\textrm{\tiny-0}}^{\textrm{\tiny+\mathcal{R}e}} with pulls above 2 $\sigma$, and {\cal U}_{\textrm{\tiny+-}}^{\textrm{\tiny+\mathcal{R}e}}, {\cal U}_{\textrm{\tiny+-}}^{\textrm{\tiny+\mathcal{I}m}} with pulls near or above 3 $\sigma$.

We can separate the coefficients related to the dynamics of the $\rho\pi$ intermediate state, $\vec{\cal O}_{\rho\pi}=({\cal U}^{+}_{i},{\cal U}^{-}_{i},{\cal I}^{i})_{(i=+,0,-)}$ (quasi-two-body terms or Q2B), from the parameters describing the interference pattern, $\vec{\cal O}_{\rm interf}=({\cal U}^{\pm,\mathcal{R}e(\mathcal{I}m)}_{i},{\cal I}^{\mathcal{R}e(\mathcal{I}m)}_{ij})$ (interference terms), as can be seen in Eq. (41). The indirect determination of most of the Q2B coefficients, ${\cal U}^{i}_{j}$ and ${\cal I}_{j}$, related to the $B^{0}\!\rightarrow\rho^{\pm}\pi^{\mp}$ intermediate states, deviates from the direct measurement by more than 2$\sigma$. This is also the case for the corresponding interference terms ${\cal U}^{+\mathcal{R}e/\mathcal{I}m}_{ij}$, which represent the real and imaginary parts of the combination of $B^{0}\rightarrow\rho^{\pm}\pi^{\mp}$ amplitudes, $(A^{+}A^{-^{*}}+{\bar{A}}^{+}{\bar{A}}^{-^{*}})$. The largest deviation is observed for the $U^{+\mathcal{R}e}_{+-}$ coefficient which exceeds 3$\sigma$. Figs. 31 and 32 display the indirect determinations of the form-factor bilinear coefficients for quasi-two-body and interference coefficients, respectively.

In order to understand these discrepancies, we can rely on two parametrisations focusing on the dynamics of the $\rho\pi$ intermediate state through Q2B observables. A first, “theoretical”, parametrisation was proposed in Ref. [83]. The time-dependent $CP$ asymmetry given by Eq. (10) is written independently for the three neutral $B^{0}\!\rightarrow\rho^{i}\pi^{j}$ intermediate states:

[TABLE]

where ${\cal C}^{i}$ and ${\cal S}^{i}$ are direct and mixing-induced $CP$ asymmetries, respectively, and the subscript $i$ refers to the charge of the emitted $\rho^{i}$ meson.

The $\rho^{\pm}\pi^{\mp}$ final state being a $CP$ admixture, the parameters ${\cal C}^{\pm}$ are not the only measurements of direct $CP$ violation. Considering either the decay where the $\rho^{\pm}$ meson is emitted by the spectator interaction or by the $W$ exchange, we can use one of the following time-integrated flavour-independent asymmetries:

[TABLE]

to constrain the $B^{0}\rightarrow\rho^{\pm}\pi^{\mp}$ amplitude system completely.

Due to the overlapping final states in the Dalitz analysis, the separation of the three intermediate states is experimentally challenging and prevents a clean direct measurement of the above Q2B asymmetries. The previous “theoretical” set of observables can be traded for a second, “experimental”, parametrisation in order to avoid this problem. Focusing only on the region of the charged $\rho$ meson in the Dalitz plane and neglecting the interferences and the residual contribution from the colour-suppressed $B^{0}\!\rightarrow\rho^{0}\pi^{0}$ mode, a specific time-dependent measurement of the $CP$-mixed $B^{0}\!\rightarrow\rho^{\pm}\pi^{\mp}$ intermediate state was first performed at $B$-factories in 2005 [52, 84]. Introducing the flavour-integrated charge asymmetry:

[TABLE]

the $B^{0}(t)\rightarrow\rho^{\pm}\pi^{\mp}$ decay rate is written as

[TABLE]

where $q_{B}=+1(-1)$ for $B^{0}({\bar{B}}^{0})$ and $q_{\rho}$ is the electric charge of the emitted $\rho$-meson. The $CP$-violating ($CP$-conserving) parameters ${\cal C}$ and ${\cal S}$ ($\Delta{\cal C}$ and $\Delta{\cal S}$) generate the flavour-dependent asymmetries:

[TABLE]

while the flavour-independent $CP$ asymmetries defined in Eq. (95) can be rewritten in terms of the experimental parameters as:

[TABLE]

In order to remove experimental biases due to the imperfect separation of the intermediate states, one can derive the Q2B parameters from the subset $\vec{\cal O}_{\rho\pi}$ of the form-factor coefficients:

[TABLE]

Finally, we can introduce the relative contribution of the colour-suppressed $B^{0}\!\rightarrow\rho^{0}\pi^{0}$ mode to the overall $B^{0}\!\rightarrow(\rho\pi)^{0}$ branching fraction:

[TABLE]

We can then reparametrise the $B^{0}\rightarrow\pi^{+}\pi^{-}\pi^{0}$ amplitude system in terms of a selected set of ${\cal U}$, ${\cal I}$ and Q2B observables: on hand, $\vec{\cal O}_{\rm interf}$, $\vec{\cal O}_{\rho^{0}\pi^{0}}=\{{\cal C}^{0},S^{0},f^{00}\}$ and on the other hand, either the theoretically motivated set $\vec{\cal O}_{\rho^{\pm}\pi^{\mp}}=\{{\cal C}^{\pm},S^{\pm},A\}$ or the experimentally convenient choice $\vec{\cal O}_{\rho^{\pm}\pi^{\mp}}=\{C,S,\Delta C,\Delta S,A\}$ (where $A$ stands for $A^{+}$, $A^{-}$ or $A_{\rho\pi}$). The world averages of the Q2B parameters derived from the measured ${\cal U}$ and ${\cal I}$ coefficients are listed in Tab. 17 in App. B, together with their corresponding prediction in the ${SU}\mkern-2.0mu(2)$ isospin framework. These results are also shown in Figs. 33 and 34.

A reasonable agreement is observed for the parameters related to the $B^{0}\rightarrow\rho^{0}\pi^{0}$ mode. The ${SU}\mkern-2.0mu(2)$ isospin fit favours large values for both direct and mixing-induced $CP$ asymmetry parameters in the $B^{0}\rightarrow\rho^{\pm}\pi^{\mp}$ decays, in contrast with experimental data. The overall 3 $\sigma$ discrepancy is mostly reflected by the correlated flavour-independent $CP$ asymmetries, $A^{\pm}$, or equivalently the charge asymmetry, $A_{\rho\pi}$. While the predicted $CP$-violating asymmetry averages, ${\cal C}$ and ${\cal S}$, are in a reasonable agreement with the experimental data, the $CP$-conserving terms, $\Delta{\cal C}$ and $\Delta{\cal S}$, deviate by more than 2.5 $\sigma$.

6.3.2 Role of the strong phases in the $B^{0}\rightarrow(\rho\pi)^{0}$ analysis

In order to understand this discrepancy in more detail, the effective weak angles associated with the $B^{0}\rightarrow\rho^{+}\pi^{-}$ and $B^{0}\rightarrow\rho^{-}\pi^{+}$ decays,

[TABLE]

constitute further interesting quantities as they both reduce to $\alpha$ in the limit of a vanishing penguin contribution. The time-independent Q2B observables in the $B^{0}\rightarrow\rho^{\pm}\pi^{\mp}$ decays do not lead to a determination of each effective mixing angle: only the average $\alpha_{\rm eff}=(\alpha^{+}+\alpha^{-})/2$ can be measured up to a fourfold ambiguity in the $[0,90]^{\circ}$ range [61]. Indeed, the Q2B observables can be combined to measure the two shifted phases:

[TABLE]

where the phase shift:

[TABLE]

cannot be determined from the Q2B data alone but cancels in the average $\alpha_{\rm eff}$. A bound on the deviation $\Delta\alpha$ between the effective angle $\alpha_{\rm eff}$ and the weak phase $\alpha$ can be derived [6, 85, 86] from the penguin relation Eq. (14):

[TABLE]

where ${\cal B}^{00}$ and ${\cal B}^{\pm\mp}$ denote the branching fractions of the $B^{0}\rightarrow\rho^{0}\pi^{0}$ and $B^{0}\rightarrow\rho^{\pm}\pi^{\mp}$ decays, respectively. This relation originates with the fact that the sum of the penguin amplitudes $({P}^{+-}+{P}^{-+})=-{P}^{00}$ is bounded by the branching fraction of the colour-suppressed mode while the difference $({P}^{+-}-{P}^{-+})$ vanishes in relation with the charge-integrated $CP$ asymmetry $A_{\rho\pi}$ (see also the discussion of the Q2B analysis of the charmless $B^{0}\rightarrow a_{1}^{\pm}\pi^{\mp}$ decay in Appendix C).

The constraint on $\Delta\alpha$ given in Fig. 35 features a solution consistent with 0:

[TABLE]

and a mirror solution around $100^{\circ}$. Including the indirect determination $\alpha_{\rm ind}$, Eq. (2), as an additional input, a tighter constraint is obtained, |\Delta\alpha|_{\rm min}=$$(3.5^{\scriptscriptstyle+2.4}_{\scriptscriptstyle-2.6})\textrm{{}^{\circ}} , and a new mirror solution, disfavoured by the $B^{0}\rightarrow(\rho\pi)^{0}$ data, emerges around $45^{\circ}$. As shown in Fig. 36, when using only the Q2B data set $\{{\cal U}^{i}_{j},{\cal I}_{j}\}_{i,j=(+,-)}$, four solutions are found for $\alpha_{\rm eff}$ around $0^{\circ}$, $45^{\circ}$ and $90^{\circ}$ (the peak near $45^{\circ}$ combines two quasi-degenerate solutions). The solution close to $90^{\circ}$ is consistent with the indirect determination of $\alpha$ Eq. (2). However, only the solution for $\alpha_{\rm eff}$ near 45∘ remains once the Q2B observables and the interference-related constraints for $B^{0}\rightarrow\rho^{\pm}\pi^{\mp}$ are considered altogether. The solutions close to $0^{\circ}$ and 90∘ are disfavoured at 2.5$\sigma$. The rejection of these solutions is slightly reinforced with the addition of the observables related to $B^{0}\rightarrow\rho^{0}\pi^{0}$.

As can be expected, the interference observables provide further constraints on the relative phases of the decay amplitudes. In particular, one can easily show from Eq. (41) that the phase shift $2\delta$, which vanishes in the absence of the penguin contribution, can be determined from the interference-related coefficients ${\cal U}^{\pm\mathcal{I}m(\mathcal{R}e)}_{+-}$:

[TABLE]

Simultaneously, the related phase between the conjugate modes:

[TABLE]

can be extracted from the same coefficients:

[TABLE]

meaning that the phase average $\delta_{\rm eff}=({\bar{\delta}}+\delta)/2$ can be extracted independently using either the Q2B coefficients or the interference-related coefficients for $B^{0}\rightarrow\rho^{\pm}\pi^{\mp}$:

[TABLE]

The Q2B coefficients provide also a constraint on the average mixing angle $\alpha_{\rm eff}$:

[TABLE]

Fig. 37 compares the two-dimensional constraints in the ($\alpha_{\rm eff}$,$\delta_{\rm eff}$) plane using either subset of observables. Clearly, only the solutions for $\alpha_{\rm eff}$ near $45^{\circ}$ are favoured by both subsets.

As discussed previously, the two observables ${\cal U}^{+\mathcal{R}e(\mathcal{I}m)}_{+-}$ play an important role in the discrepancy between the $B^{0}\rightarrow(\rho\pi)^{0}$ Dalitz data and the isospin hypothesis. It is interesting to show how removing these two observables impacts the analysis: the $\alpha_{\rm eff}$ solutions near $0^{\circ}$ and $90^{\circ}$ get favoured, and the direct determination of $\alpha$ agrees much better with the indirect value, as illustrated in Fig. 38.

From the experimental side, the observables {\cal U}_{\textrm{\tiny+-}}^{\textrm{\tiny\pm$$\mathcal{R}e$($\mathcal{I}m)}} and {\cal I}_{\textrm{\tiny\pm}}^{\textrm{\tiny\mathcal{R}e$($\mathcal{I}m)}} are the coefficients of the mixed form-factors bilinear ($f^{+}f^{-^{*}}$). Therefore, they are the only observables sensitive to a possible phase between the functionals describing the line-shape of the $\rho$-mesons of opposite charge. It would be very interesting to study the role played by the isobar approximation and to determine the possible re-interpretations of the experimental data, but this lies clearly beyond the scope of the isospin analysis presented here, and we will refrain from drawing further conclusions on this discrepancy.

6.3.3 $B\!\rightarrow\rho\pi$ pentagonal analysis

The absolute scale of the neutral $B^{0}\rightarrow\pi^{-}\pi^{+}\pi^{0}$ amplitudes can be determined by adding the measured branching fractions for the two charged modes $B^{+}\rightarrow\rho^{0}\pi^{+}$ and $B^{+}\rightarrow\rho^{+}\pi^{0}$, related to the neutral modes via the pentagonal relation Eq. (13). Once the data listed in Tab. 8 is added to the fit, individual $B\rightarrow\rho\pi$ branching fractions and direct $CP$ asymmetries in the charged modes can be determined. The corresponding pulls are given in Tab. 18 in App. B.

While the predicted branching fraction for the neutral $B^{0}\rightarrow\rho^{0}\pi^{0}$ mode is in very good agreement and competitive with the experimental measurement, a discrepancy is observed in the balance of branching fractions for the charged modes and the mixed $B^{0}\rightarrow\rho^{\pm}\pi^{\mp}$ decays. The branching fraction of both charged modes are predicted to be larger than their measured value, which is compensated by a lower $B^{0}\rightarrow\rho^{\pm}\pi^{\mp}$ branching ratio. Figs. 39 and 40 display the ${SU}\mkern-2.0mu(2)$ isospin constraints for the neutral and charged modes, respectively.

7 Prospective study

In this section, we discuss how improved measurements of some $B\rightarrow\pi\pi$, $B\rightarrow\rho\rho$, $B\rightarrow\rho\pi$ observables can affect the accuracy on the CKM angle $\alpha$. So far, the experimental data is statistically limited for all the $B\rightarrow hh$ charmless modes considered here. Our prospective study aims at identifying the specific decay channels worth measuring more accurately to improve the resolution on $\alpha$ significantly.

For simplicity, we adopt a systematic approach rather than relying on the expected performance of current or forthcoming flavour experiments such as LHCb or Belle II. We consider the subset of observables related to a specific $B\rightarrow h^{i}h^{j}$ decay, and we determine how much the accuracy on $\alpha$ is improved if we reduce the uncertainties for this subset of observables by the (arbitrary) factor $\sqrt{2}$, or if we take the radical limit of setting these uncertainties to zero. The central value of the current world-average measurements and correlation coefficients are kept unchanged. For the neutral modes $B^{0}\rightarrow h^{+}h^{-}$ and $B^{0}\rightarrow h^{0}h^{0}$, we distinguish the case of a counting analysis (C) that gives access to the branching ratio only (e.g., the LHCb measurement of $B^{0}\rightarrow\rho^{0}\rho^{0}$ [49]), and the case of flavour-tagged analyses: indeed, a time-integrated (TI) analysis extracts only the direct $CP$ asymmetry whereas a time-dependent (TD) analysis yields both the direct and the mixing-induced $CP$-asymmetry parameters (e.g., the LHCb contribution to the study of $B^{0}\rightarrow\pi^{+}\pi^{-}$ [33]).

The global combination of the decay-specific determinations of $\alpha$ is so far dominated by the $B\rightarrow\rho\rho$ data that provides a constraint on $\alpha$ with a relative uncertainty at the level of 5%, i.e., \alpha_{\rho\rho}=$$(92.1^{\scriptscriptstyle+4.6}_{\scriptscriptstyle-4.9})\textrm{{}^{\circ}} for the solution near $90^{\circ}$ given in Eq. (37). As shown in Tab. 12, if all other observables remain unchanged, improving the accuracy of the branching ratio of the charged mode $B^{+}\rightarrow\rho^{+}\rho^{0}$ would improve the resolution on $\alpha$ only marginally, even in the case of a vanishing resolution (indicating that this observable has only a limited impact on the accuracy for $\alpha$). Improving the measurements for the neutral modes, in particular the colour-suppressed $B^{0}\rightarrow\rho^{0}\rho^{0}$ decay, has a larger impact, essentially driven by the $CP$-asymmetries parameters. Improving the time-dependent asymmetries in the $B^{0}\rightarrow\rho^{0}\rho^{0}$ is also worth investigating, e.g., in the second run of LHCb data taking. Reducing by $\sqrt{2}$ all the $B\rightarrow\rho\rho$ uncertainties would reduce the 68% CL interval for $\alpha$ by more than 1 degree.

The subleading contribution to the combined $\alpha$ determination is provided by the $B\rightarrow\pi\pi$ system, $\alpha_{\pi\pi}=(93.0\pm 14.0)^{\circ}$, see Eq. (35). As shown on Tab. 13, any sizable improvement on $\alpha$ is driven by the increased accuracy in the measurement of the direct $CP$ asymmetry in the colour-suppressed decay $B^{0}\rightarrow\pi^{0}\pi^{0}$. A reduction by $\sqrt{2}$ of the uncertainty of this single observable would reduce the 68% CL range for $\alpha$ by $1^{\circ}$. A similar improvement of all the measured $B\rightarrow\pi\pi$ observables would reduce by about $1.5^{\circ}$ the uncertainty on $\alpha$.

No time-dependent analysis has become available so far for the colour-suppressed mode due to the dominant di-photon decay of the neutral pions that hinders the measurement of the decay time, so that ${\cal S}^{00}_{\pi\pi}$ is not measured yet. However, the low branching fraction of the Dalitz decay $\pi^{0}\rightarrow\gamma e^{+}e^{-}$ or the photon conversion in the detector material may allow one to measure the $B^{0}$ decay time in the future high-statistics flavour experiments (LHCb upgrade, Belle II). Fig. 41 shows the $p$-value for $\alpha$ when adding the $S^{00}_{\pi\pi}$ observable under different hypotheses concerning the experimental resolution, and setting the central value as predicted in Eq. (93). In addition to a significant improvement on $\alpha$, the measurement of ${\cal S}^{00}_{\pi\pi}$ would reduce the number of mirror solutions in the $B\rightarrow\pi\pi$ isospin analysis (see the discussion around Eq. (32)).

The same prospective exercise is more delicate for the $B^{0}\rightarrow\pi^{+}\pi^{-}\pi^{0}$ system due to the discrepancy between the direct measurement $\alpha_{\rho\pi}$=(54.1^{\scriptscriptstyle+7.7}_{\scriptscriptstyle-10.3})\textrm{{}^{\circ}} $\cup$ (141.8^{\scriptscriptstyle+4.8}_{\scriptscriptstyle-5.4})\textrm{{}^{\circ}} (Eq. (43)) and the indirect determination from the global CKM fit (Eq. (2)). New measurements of this decay would certainly aim first at a better understanding of this discrepancy, rather than improving the accuracy on $\alpha$ extracted from this channel. Assuming the central values given by the individual predictions for each ${\cal U}$ and ${\cal I}$ coefficient, as listed in Tab. 16, and keeping unchanged the current experimental resolution and correlations, we obtain a mild constraint on $\alpha$, consistent with $\alpha_{\rm ind}$, namely $\alpha_{\rho\pi}^{{\rm fit}~{}{\cal U},{\cal I}}$=(82^{\scriptscriptstyle-33}_{\scriptscriptstyle-48})\textrm{{}^{\circ}} . In this case, the reduction of the observable uncertainty by $\sqrt{2}$ results in a reduction by approximately $9^{\circ}$ for the 68% CL interval for $\alpha$. On the other hand, keeping the current world-average central values and reducing the uncertainty on all the ${\cal U}$ and ${\cal I}$ observables by a factor $\sqrt{2}$ leads to a more stringent 68% CL interval $\alpha_{\rho\pi}^{\sigma\scriptscriptstyle/\sqrt{2}}$=(54.1^{\scriptscriptstyle+5.6}_{\scriptscriptstyle-6.9})\textrm{{}^{\circ}} $\cup$ (141.0^{\scriptscriptstyle+1.9}_{\scriptscriptstyle-1.9})\textrm{{}^{\circ}} , which remains difficult to interpret due to the discrepancies already discussed between the current data and our theoretical framework based on isospin symmetry. Additional measurements should hopefully provide a clearer and more consistent picture for the $B\rightarrow\rho\pi$ sector before discussing any improvement in the extraction of $\alpha$ from these modes.

8 Conclusion

Quark flavour transitions provide particularly stringent tests of the Standard Model, both through rare decays and ${\cal CP}$-violating processes. An accurate knowledge of the Cabibbo–Kobayashi–Maskawa matrix is essential for these studies and it requires the combination of many precise constraints. We have focussed on the determination of the $\alpha$ angle, which can be extracted with a high accuracy from two-body charmless $B$-meson decays extensively studied at $B$-factories and LHCb.

We recalled that this extraction can be done from $B\rightarrow\pi\pi$, $B\rightarrow\rho\pi$ and $B\rightarrow\rho\rho$ decays but it is affected by the presence of penguin contributions. We explained how ${SU}\mkern-2.0mu(2)$ isospin symmetry can be used to constrain the structure of hadronic penguin and tree amplitudes, enabling the extraction of $\alpha$ from branching ratios and ${\cal CP}$ asymmetries. We gave details on the analyses of $B\rightarrow\pi\pi$, $B\rightarrow\rho\pi$ and $B\rightarrow\rho\rho$ systems separately, before combining these results to reach an accuracy around 4∘ on the direct determination of $\alpha$. The $B\rightarrow\pi\pi$ and $B\rightarrow\rho\rho$ systems dominate the combination and they favour solutions in good agreement with the indirect determination of $\alpha$ from a global CKM fit analysis, Eq. (2), it is not the case for the $B\rightarrow\rho\pi$ system which favours different ranges of values for $\alpha$ with a discrepancy at the level of 3 $\sigma$ compared to the indirect determination. The combination of the three channels is dominated by $B\rightarrow\rho\rho$, and to a lesser extent $B\rightarrow\pi\pi$, resulting in the 68% CL confidence intervals given in Eq. (1).

We have then studied several uncertainties that may affect this extraction. We have tested the hypotheses underlying the ${SU}\mkern-2.0mu(2)$ isospin analysis: setting the \Delta I=\raisebox{1.72218pt}{\footnotesize3}\negthinspace/\negthinspace\raisebox{-2.58334pt}{\footnotesize2} electroweak penguins to zero, neglecting the difference of light-quark masses generating $\pi^{0}-\eta-\eta^{\prime}$ mixing, setting the $\rho$ width to zero to cancel $\Delta I=1$ contributions thanks to Bose–Einstein symmetry. These effects may shift the central value of $\alpha_{\rm dir}$ by around $2^{\circ}$, while keeping the uncertainty around $4^{\circ}$ to $5^{\circ}$, thus remaining within the statistical uncertainty quoted in Eq. (1). In addition, we have discussed a few aspects concerning the statistical treatment used in order to extract the $p$-value, comparing two statistical approaches to test the impact of hadronic nuisance parameters on coverage. The approach based on Wilks’ theorem and mainly used here proves to be more conservative than the bootstrap method for $B\rightarrow\pi\pi$ and $B\rightarrow\rho\rho$, but less conservative in the case of $B\rightarrow\rho\pi$. The comparison of the coverage properties of the two approaches leads to a further uncertainty of around $1^{\circ}$. We stress that this uncertainty is only attached to the combination of the three direct determinations of $\alpha$ with the current data: it is likely to be reduced if one combines the direct determination of $\alpha$ with other observables, leading to a more accurate determination of $\alpha$ (which is in particular the typical case of global CKM fits).

Assuming the validity of the Standard Model and taking as an input the indirect determination of $\alpha$ from the global CKM fit, we used the observables in $B\rightarrow\pi\pi$, $B\rightarrow\rho\pi$ and $B\rightarrow\rho\rho$ decays to extract information on ratios of hadronic amplitudes (penguin-to-tree and colour-suppressed), finding results in broad agreement with expectations from QCD factorisation, apart from the ratio of the colour-suppressed to colour-allowed tree $\pi\pi$ contributions: indeed both the phase and the modulus of this ratio do not agree well with theoretical expectations. It would be interesting to widen this discussion and see how various theoretical approaches to non-leptonic two-body decays can reproduce the patterns of hadronic amplitudes that we have extracted from the data.

Under the same hypotheses, we have also performed the indirect determinations of observables of interest (using all the other measurements available), comparing the pulls for those already observed and predicting the values of the remaining ones. The compatibility between direct measurements and indirect determinations is very good for the observables in the $B\rightarrow\pi\pi$ and $B\rightarrow\rho\rho$ systems, whereas we could identify a subset of $B\rightarrow\rho\pi$ observables likely to be responsible for the discrepancies observed with respect to the Standard Model expectations in these modes. Among many other quantities, we have predicted the yet-to-be-measured mixing-induced ${\cal CP}$ asymmetry in the $B^{0}\!\rightarrow\pi^{0}\pi^{0}$ decay; see Eq. (32).

Finally, we have performed a prospective study to analyse how improved measurements for some subsets of observables can improve the uncertainty of $\alpha$. In particular, we have noticed that an improved accuracy for the time-dependent asymmetries in $B^{0}\rightarrow\rho^{0}\rho^{0}$ and the measurement of ${\cal S}^{00}_{\pi\pi}$ would reduce the uncertainty on $\alpha$ in a noticeable way.

We have seen that the extraction of $\alpha$ is now possible to a high accuracy, using many different channels and experimental sources. It would be very interesting to measure the remaining observables that we can accurately predict using the data already available. The current accuracy reached by $\alpha$ makes it a particularly useful constraint both for the Standard Model and for searches of New Physics, although the theoretical computations for these transitions remain challenging. Some of these channels will be improved by the LHCb experiment. In addition, the advent of the Belle-II experiment will certainly lead to new and improved measurements for a large set of branching ratios and ${\cal CP}$ asymmetries, providing an opportunity to understand better the results obtained for the $B\rightarrow\rho\pi$ system and allowing us to extract hadronic parameters with a higher accuracy. Both avenues should be highly beneficial for the upcoming studies of flavour physics in the quark sector.

Acknowledgements

We would like to thank all our collaborators from the CKMfitter group for many useful discussions on the statistical issues covered in this article, and in particular K. Trabelsi and L. Vale Silva for a careful reading of the manuscript. SDG acknowledges partial support from Contract FPA2014-61478-EXP. This project has received funding from the European Union’s Horizon 2020 research and innovation programme under Grant agreements nos 690575, 674896 and 692194.

Appendix A Results from the BaBar and Belle experiments

Both BaBar and Belle experiments have measured the relevant observables for an extraction of the weak phase $\alpha$ in each of the three charmless decay systems. We have provided an analysis of their average in Sec. 3, and we briefly discuss the separate results from each $B$-factory in this appendix.

A.1 $B\rightarrow\pi\pi$ and $B\rightarrow\rho\rho$ analyses

The individual measurements of the $B\rightarrow\pi\pi$ and $B\rightarrow\rho\rho$ observables are shown in Figs. 42 and 43, and Figs. 44 and 45, respectively. For each observable, a good agreement is observed between the different sources of experimental data (BaBar, Belle, and when available, other experiments).

The individual determinations of $\alpha$ based on the BaBar and Belle data separately are shown in Fig. 46. The corresponding 68% CL intervals are:

[TABLE]

A.2 $B\rightarrow\rho\pi$ analysis

Both B-factories have performed a full Dalitz analysis of the $B^{0}\rightarrow\pi^{+}\pi^{-}\pi^{0}$ decay using the same ${\cal U}$ and ${\cal I}$ observables. The measurements of the Q2B- and interference-related coefficients are summarised in Fig. 48 and Fig. 47, respectively. Good agreement between the two experiments is observed. The correlated average gives $\chi^{2}/n_{\rm dof}=18.1/26$.

The individual determinations of $\alpha_{\rho\pi}$ based on BaBar and Belle data separately are shown in Fig. 49 (left panel). The corresponding 68% intervals are

[TABLE]

The $B^{0}\rightarrow\pi^{+}\pi^{-}\pi^{0}$ data from Belle and BaBar are consistent with the indirect determination $\alpha_{\rm ind}$ at the level of 1.3 and 2.6 standard deviations, respectively. Being the projection of the combination of constraints in a multi-dimensional parameter space, the constraint on $\alpha$ that results from the averaged data cannot be interpreted as the direct combination of the one-dimensional constraints. A more explicit picture is obtained when representing the constraint in the ($\alpha_{\rm eff}$,$\delta_{\rm eff}$) plane, as shown in Fig. 49 (right panel). The preferred $\alpha$ solution near $50^{\circ}$ clearly appears as the solution favoured by both Belle and BaBar experiments.

A.3 Combined analysis

The combined constraints on $\alpha$ using BaBar and Belle data separately are shown in Fig. 50. The corresponding 68% CL intervals are

[TABLE]

The agreement with the indirect $\alpha_{\rm ind}$ determination is 0.9$\sigma$ and 1.5$\sigma$ for BaBar and Belle data, respectively. One can also notice that BaBar and Belle data do not favour the same peak around $0^{\circ}$ or $90^{\circ}$, even though both intervals in $\alpha$ are acceptable at 95% CL.

Appendix B Numerical tables for the prediction of observables

In the following tables, we present the results for the indirect determination of various observables discussed in Sec. 6. For each quantity, we present the world average of the available measurements, the indirect determination from the ${SU}\mkern-2.0mu(2)$ isospin analysis without including this observable, and the compatibility pull between the two determinations, obtained by comparing the minimum value of the $\chi^{2}$ function with and without the experimental measurement of this quantity, see Eq. (92) (Tabs. 14, 15, 16, 17, 18).

Appendix C Quasi-two-body analysis of $B^{0}\rightarrow a_{1}^{\pm}\pi^{\mp}$

As discussed in Sec. 6.3.2, the quasi-two-body (Q2B) analysis of the neutral $B^{0}\rightarrow\rho^{\pm}\pi^{\mp}$ decay provides enough information to extract both the average mixing angle $\alpha_{\rm eff}$, which coincides with the CKM angle $\alpha$ in the limit of a vanishing penguin contribution, and the phase shift $\delta_{\rm eff}$, corresponding to the phase between $\rho^{+}\pi^{-}$ and $\rho^{-}\pi^{+}$ decay amplitudes. A similar analysis can be performed in the case of the charmless $B^{0}\rightarrow a_{1}^{\pm}\pi^{\mp}$ decay. Fig. 51 displays the two-dimensional constraint in the ($\alpha_{\rm eff}$,$\delta_{\rm eff}$) plane (left panel) and the corresponding one-dimensional projection on $\alpha_{\rm eff}$ (right panel) when using the world-average $B^{0}\rightarrow a_{1}^{\pm}\pi^{\mp}$ Q2B parameters collected in Ref. [68]. The 68% CL intervals on $\alpha_{\rm eff}$ are:

[TABLE]

The two mirror solutions close to $0^{\circ}$ and $90^{\circ}$ are consistent with a vanishing average phase shift $\delta_{\rm eff}$.

The current experimental limit [68] on the branching fraction of the colour-suppressed mode, ${\cal B}(B^{0}\rightarrow a_{1}^{0}\pi^{0})<1.1\times 10^{-3}$ at 90% CL, does not allow us to derive a ${SU}\mkern-2.0mu(2)$ bound on the difference $\Delta\alpha=\alpha-\alpha_{\rm eff}$. Considering the flavour-related modes $B\rightarrow a_{1}K$, a bound based on the ${SU}\mkern-2.0mu(3)$ symmetry can, however, be defined [87]. Such an analysis, performed by the BaBar collaboration [88], leads to the constraint $|\Delta\alpha|=|\alpha-\alpha_{\rm eff}({a_{1}\pi})|<13^{\circ}$ at 90% CL, consistent with the value from Eq. (106) obtained for the $B^{0}\rightarrow(\rho\pi)^{0}$ decay. Assuming for simplicity that the latter deviation $(\Delta\alpha)_{\rho\pi}$ has the same value in the case of the $B\rightarrow a_{1}\pi$ decay, the constraint on $\alpha$ illustrated by Fig. 52 gives :

[TABLE]

The solution near $90^{\circ}$ is consistent with the indirect $\alpha$ determination given by Eq. (2). However, this constraint is based on a hypothesis with a poor theoretical motivation, and it will not be included in our combination of direct $\alpha$ determinations (whereas we include information from the $\pi\pi$, $\rho\rho$ and $\rho\pi$ modes).