text:
# General analysis of time evolution of decay spectrum in 𝐵⁰→𝜋⁺𝜋⁻ℓ⁺ℓ⁻

## 1 Introduction

CP violation has been observed in the neutral kaon systems . Great effort is devoted to the study of possible signals of CP violation in $`B`$ system, which will provide invaluable information about the origin of CP violation and is one of the most promising research area having the potential of establishing new physics beyond the standard model (SM). Started operating, two $`B`$–factories BaBar and Belle open an excited era for a comprehensive study of $`B`$ meson physics and especially its rare decays. The main physics program of these factories constitutes a detailed study of CP violation in $`B_d`$ meson and precise measurement of rare flavor–changing neutral current (FCNC) processes. Furthermore, the ultimate goal of these studies is to look for the inconsistencies within the SM (see for example ), in particular, to find indications for indications for new physics in the flavor and CP violating sectors. The factories mentioned above have both already signaled the first evidences about the CP violation in neutral $`B`$ meson decays .

One efficient way for detecting CP violation in the neutral $`B`$ meson decays is the well known Dalitz plot asymmetry . The main advantage of the Dalitz plot asymmetries from the partial rate asymmetries is that they might be present even when partial rate asymmetries vanish. It should be noted that such CP–violating asymmetries in the angular variable $`\phi `$ ($`\phi `$ is the angle between the $`\mathrm{}^+\mathrm{}^{}`$ and $`\pi ^+\pi ^{}`$ planes) in the $`K\pi ^+\pi ^{}\mathrm{}^+\mathrm{}^{}`$ decay was measured . In a recent report, NA48 Collaboration  confirms large CP–violating effect in the $`K\pi ^+\pi ^{}\mathrm{}^+\mathrm{}^{}`$ decay. Time evolution of the decay spectrum in $`K^0(\overline{K}^0\pi ^+\pi ^{}\mathrm{}^+\mathrm{}^{})`$ was investigated in framework of the SM in .

In this connection several interesting questions come into mind immediately: how does the angular variable asymmetry $`𝒜_\phi `$ evolve with time for the neutral $`B`$ meson decay? How do new physics effects change the Dalitz plot asymmetry? Can we measure this asymmetry in $`B`$–factories? This paper is devoted to answering these questions by analyzing the time dependence of the decays $`B^0,\overline{B}^0\pi ^+\pi ^{}\mathrm{}^+\mathrm{}^{}`$ in detail. It should be noted that the $`B\pi ^+\pi ^{}\mathrm{}^+\mathrm{}^{}`$ decay in SM and beyond were studied in  and , respectively.

The paper is organized as follows. In section 2 we present the general expression of the Dalitz plot asymmetries for the $`B^0(\overline{B}^0\pi ^+\pi ^{}\mathrm{}^+\mathrm{}^{})`$ decays using the most general, model independent form of the effective Hamiltonian. In section 3 we study the sensitivity of the Dalitz plot asymmetries on the new Wilson coefficients.

## 2 Time evolution of the decay spectrum of $`B^0\pi ^+\pi ^{}\mathrm{}^+\mathrm{}^{}`$

The matrix element for the $`B\rho (\pi ^+\pi ^{})\mathrm{}^+\mathrm{}^{}`$ decay is described by the $`bd\mathrm{}^+\mathrm{}^{}`$ transition at quark level. Following , the effective Hamiltonian of the $`bd\mathrm{}^+\mathrm{}^{}`$ transition can be written as the sum of the SM and new physics contribution, in a general model independent way in the following form

$`_{eff}`$ $`=`$ $`{\displaystyle \frac{G\alpha }{\sqrt{2}\pi }}V_{td}V_{tb}^{}\{C_{SL}\overline{d}i\sigma _{\mu \nu }{\displaystyle \frac{q^\nu }{q^2}}Lb\overline{\mathrm{}}\gamma ^\mu \mathrm{}+C_{BR}\overline{d}i\sigma _{\mu \nu }{\displaystyle \frac{q^\nu }{q^2}}Rb\overline{\mathrm{}}\gamma ^\mu \mathrm{}`$
$`+C_{LL}^{tot}\overline{d}_L\gamma _\mu b_L\overline{\mathrm{}}_L\gamma ^\mu \mathrm{}_L+C_{LR}^{tot}\overline{d}_L\gamma _\mu b_L\overline{\mathrm{}}_R\gamma ^\mu \mathrm{}_R+C_{RL}\overline{d}_R\gamma _\mu b_R\overline{\mathrm{}}_L\gamma ^\mu \mathrm{}_L`$
$`+C_{RR}\overline{d}_R\gamma _\mu b_R\overline{\mathrm{}}_R\gamma ^\mu \mathrm{}_R+C_{LRLR}\overline{d}_Lb_R\overline{\mathrm{}}_L\mathrm{}_R+C_{RLLR}\overline{d}_Rb_L\overline{\mathrm{}}_L\mathrm{}_R`$
$`+C_{LRRL}\overline{d}_Lb_R\overline{\mathrm{}}_R\mathrm{}_L+C_{RLRL}\overline{d}_Rb_L\overline{\mathrm{}}_R\mathrm{}_L+C_T\overline{d}\sigma _{\mu \nu }b\overline{\mathrm{}}\sigma ^{\mu \nu }\mathrm{}`$
$`+iC_{TE}ϵ^{\mu \nu \alpha \beta }\overline{d}\sigma _{\mu \nu }b\overline{\mathrm{}}\sigma _{\alpha \beta }\mathrm{}\},`$
where the chiral projection operators $`L`$ and $`R`$ in (2) are defined as

$`L={\displaystyle \frac{1\gamma _5}{2}},R={\displaystyle \frac{1+\gamma _5}{2}},`$
and $`C_X`$ are the coefficients of the four–Fermi interactions. Here $`q=p_++p_{}`$ and $`p_+`$ and $`p_{}`$ are the four momenta of $`\mathrm{}^+`$ and $`\mathrm{}^{}`$, respectively. The first two of these coefficients, $`C_{SL}`$ and $`C_{BR}`$, are the nonlocal Fermi interactions which correspond to $`2m_sC_7^{eff}`$ and $`2m_bC_7^{eff}`$ in the SM, respectively. The next four terms with coefficients $`C_{LL}`$, $`C_{LR}`$, $`C_{RL}`$ and $`C_{RR}`$ describe the vector type interactions. Two of these vector interactions containing $`C_{LL}^{tot}`$ and $`C_{LR}^{tot}`$ do already exist in the SM in combinations of the form $`(C_9^{eff}C_{10})`$ and $`(C_9^{eff}+C_{10})`$. Therefore, $`C_{LL}^{tot}`$ and $`C_{LR}^{tot}`$ describe the sum of the contributions from SM and the new physics and they are defined as

$`C_{LL}^{tot}`$ $`=`$ $`C_9^{eff}C_{10}+C_{LL},`$
$`C_{LR}^{tot}`$ $`=`$ $`C_9^{eff}+C_{10}+C_{LR}.`$ (2)
The terms with coefficients $`C_{LRLR}`$, $`C_{RLLR}`$, $`C_{LRRL}`$ and $`C_{RLRL}`$ describe the scalar type interactions. The remaining last two terms leaded by the coefficients $`C_T`$ and $`C_{TE}`$, obviously, describe the tensor type interactions.

The matrix element for the the process $`B\pi ^+\pi ^{}\mathrm{}^+\mathrm{}^{}`$ can be obtained from the matrix element $`B\rho \mathrm{}^+\mathrm{}^{}`$, by the subsequent decay of the $`\rho `$ meson to pion pair whose interaction Hamiltonian is

$`=g_{\rho \pi \pi }(Q\epsilon ),`$ (3)
where $`\epsilon ^\mu `$ is the the polarization four vector of the $`\rho `$ meson, $`Q=p_{\pi ^+}p_\pi ^{}`$ and $`p_{\pi ^+}`$ and $`p_\pi ^{}`$ are the four momenta of $`\pi ^+`$ and $`\pi ^{}`$, respectively. Hence, in calculating the 4–body decay amplitude for the $`B\pi ^+\pi ^{}\mathrm{}^+\mathrm{}^{}`$ process, the matrix element of the $`B\rho \mathrm{}^+\mathrm{}^{}`$ decay is needed. For this purpose it is enough to sandwich the effective Hamiltonian (2) between the initial and final meson states. In other words, the following matrix elements

$`\rho \left|\overline{d}\gamma _\mu (1\pm \gamma _5)b\right|B,`$
$`\rho \left|\overline{d}(1\pm \gamma _5)b\right|B,`$
$`\rho \left|\overline{d}\sigma _{\mu \nu }(1\pm \gamma _5)b\right|B,`$
need to be calculated in order to calculate the decay amplitude for the $`B\rho \mathrm{}^+\mathrm{}^{}`$ process. These matrix elements can be written in terms of the form factors as

$`\rho (P,\epsilon )\left|\overline{d}\gamma _\mu (1\pm \gamma _5)b\right|B(p_B)=`$
$`ϵ_{\mu \nu \lambda \sigma }\epsilon ^\nu P^\lambda q^\sigma {\displaystyle \frac{2V(q^2)}{m_B+m_\rho }}\pm i\epsilon _\mu ^{}(m_B+m_\rho )A_1(q^2)`$
$`i(p_B+P)_\mu (\epsilon ^{}q){\displaystyle \frac{A_2(q^2)}{m_B+m_\rho }}iq_\mu {\displaystyle \frac{2m_\rho }{q^2}}(\epsilon ^{}q)\left[A_3(q^2)A_0(q^2)\right].`$
Multiplying both sides of Eq. (2) with $`q_\mu `$ and using equation of motion, we get

$`\rho (P,\epsilon )\left|\overline{d}(1\pm \gamma _5)b\right|B(p_B)=`$
$`{\displaystyle \frac{i}{m_b}}(\epsilon ^{}q)\left\{(m_B+m_\rho )A_1(q^2)\pm (m_B^2P^2){\displaystyle \frac{A_2(q^2)}{m_B+m_\rho }}\pm 2m_\rho \left[A_3(q^2)A_0(q^2)\right]\right\}.`$
The remaining matrix element is defined as

$`\rho (P,\epsilon )\left|\overline{d}\sigma _{\mu \nu }b\right|B(p_B)=`$
$`iϵ_{\mu \nu \lambda \sigma }\{2T_1(q^2)\epsilon _{}^{}{}_{}{}^{\lambda }(p_B+P)^\sigma +{\displaystyle \frac{2}{q^2}}(m_B^2P^2)[T_1(q^2)T_2(q^2)]\epsilon _{}^{}{}_{}{}^{\lambda }q^\sigma `$
$`{\displaystyle \frac{4}{q^2}}[T_1(q^2)T_2(q^2){\displaystyle \frac{q^2}{m_B^2P^2}}T_3(q^2)](\epsilon ^{}q)P^\lambda q^\sigma \}.`$
Moreover, the matrix element $`\rho (P,\epsilon )\left|\overline{d}\sigma _{\mu \nu }\gamma _5b\right|B(p_B)`$ can be calculated from Eq. (2) using the identity

$`\sigma _{\mu \nu }\gamma _5={\displaystyle \frac{i}{2}}ϵ_{\mu \nu \alpha \beta }\sigma ^{\alpha \beta },`$
which yields

$`\rho (P,\epsilon )\left|\overline{d}\sigma _{\mu \nu }\gamma _5b\right|B(p_B)=`$ (7)
$`2\{T_1(q^2)[\epsilon _\mu ^{}(p_B+P)_\nu \epsilon _\nu ^{}(p_B+P)_\mu ]{\displaystyle \frac{1}{q^2}}(m_B^2P^2)(\epsilon _\mu ^{}q_\nu \epsilon _\nu ^{}q_\mu )[T_1(q^2)T_2(q^2)]`$
$`+{\displaystyle \frac{2}{q^2}}[T_1(q^2)T_2(q^2){\displaystyle \frac{q^2}{m_B^2P^2}}T_3(q^2)](\epsilon ^{}q)(P^\mu q^\nu P^\nu q^\mu )\}.`$
Using Eqs. (37) for the matrix element $`B\pi ^+\pi ^{}\mathrm{}^+\mathrm{}^{}`$ decay we get

$`(B\pi ^+\pi ^{}\mathrm{}^+\mathrm{}^{})={\displaystyle \frac{G\alpha }{4\sqrt{2}}}V_{tb}V_{td}^{}{\displaystyle \frac{g^{\nu \alpha }P^\nu P^\alpha /m_\rho ^2}{P^2m_\rho ^2+im_\rho \mathrm{\Gamma }}}g_{\rho \pi \pi }Q_\alpha `$ (8)
$`\times \{\overline{\mathrm{}}\gamma ^\mu (1\gamma _5)\mathrm{}[2𝒱_{L_1}ϵ_{\mu \nu \lambda \sigma }P^\lambda q^\sigma i𝒱_{L_2}g_{\mu \nu }+i𝒱_{L_3}(2P+q)_\mu q_\nu +i𝒱_{L_4}q_\mu q_\nu ]`$
$`+\overline{\mathrm{}}\gamma ^\mu (1+\gamma _5)\mathrm{}\left[2𝒱_{R_1}ϵ_{\mu \nu \lambda \sigma }P^\lambda q^\sigma i𝒱_{R_2}g_{\mu \nu }+i𝒱_{R_3}(2P+q)_\mu q_\nu +i𝒱_{R_4}q_\mu q_\nu \right]`$
$`+\overline{\mathrm{}}(1\gamma _5)\mathrm{}(iS_L)q_\nu +\overline{\mathrm{}}(1+\gamma _5)\mathrm{}q_\nu (iS_R)`$
$`+4\overline{\mathrm{}}\sigma _{\mu \beta }\mathrm{}\left(iC_Tϵ^{\mu \beta \lambda \sigma }\right)\left[2T_1g_{\lambda \nu }(2P+q)_\sigma +B_6g_{\lambda \nu }q_\sigma B_7P_\lambda q_\sigma q_\nu \right]`$
$`+16C_{TE}\overline{\mathrm{}}\sigma ^{\mu \beta }\mathrm{}[2T_1g_{\mu \nu }(2P+q)_\beta +B_6g_{\mu \nu }q_\beta B_7q_\nu P_\mu q_\beta ]\},`$
where $`𝒱_{L_i}`$ and $`𝒱_{R_i}`$ are the coefficients of left– and right–handed leptonic currents with vector structure, respectively, $`S_{L,R}`$ are the coefficients of the leptonic scalar currents and the last two terms correspond to the contribution of the leptonic tensor currents. These new coefficients are functions of the Wilson coefficients and the form factors introduced in defining the hadronic matrix elements above. Their explicit expressions are given by

$`𝒱_{L_1}`$ $`=`$ $`(C_{LL}^{tot}+C_{RL}){\displaystyle \frac{V}{m_B+m_\rho }}2(C_{BR}+C_{SL}){\displaystyle \frac{T_1}{q^2}},`$
$`𝒱_{L_2}`$ $`=`$ $`(C_{LL}^{tot}C_{RL})(m_B+m_\rho )A_12(C_{BR}C_{SL}){\displaystyle \frac{T_2}{q^2}}(m_B^2m_\rho ^2),`$
$`𝒱_{L_3}`$ $`=`$ $`{\displaystyle \frac{C_{LL}^{tot}C_{RL}}{m_B+m_\rho }}A_22(C_{BR}C_{SL}){\displaystyle \frac{1}{q^2}}\left[T_2+{\displaystyle \frac{q^2}{m_B^2m_\rho ^2}}T_3\right],`$
$`𝒱_{L_4}`$ $`=`$ $`(C_{LL}^{tot}C_{RL}){\displaystyle \frac{2m_\rho }{q^2}}(A_3A_0)2(C_{BR}C_{SL}){\displaystyle \frac{T_3}{q^2}},`$
$`𝒱_{R_1}`$ $`=`$ $`𝒱_{L_1}(C_{LL}^{tot}C_{LR}^{tot},C_{RL}C_{RR}),`$
$`𝒱_{R_2}`$ $`=`$ $`𝒱_{L_2}(C_{LL}^{tot}C_{LR}^{tot},C_{RL}C_{RR}),`$
$`𝒱_{R_3}`$ $`=`$ $`𝒱_{L_3}(C_{LL}^{tot}C_{LR}^{tot},C_{RL}C_{RR}),`$
$`𝒱_{R_4}`$ $`=`$ $`𝒱_{L_4}(C_{LL}^{tot}C_{LR}^{tot},C_{RL}C_{RR}),`$
$`𝒮_L`$ $`=`$ $`(C_{LRRL}C_{RLRL}){\displaystyle \frac{1}{m_b}}\left[(m_B+m_\rho )A_1(m_B^2P^2){\displaystyle \frac{A_2}{m_B+m_\rho }}2m_\rho (A_3A_0)\right],`$
$`𝒮_R`$ $`=`$ $`(C_{LRLR}C_{RLLR}){\displaystyle \frac{1}{m_b}}\left[(m_B+m_\rho )A_1(m_B^2P^2){\displaystyle \frac{A_2}{m_B+m_\rho }}2m_\rho (A_3A_0)\right],`$
where dependence on $`q^2`$ is implied. To obtain the amplitudes for the $`\overline{B}^0\pi ^+\pi ^{}\mathrm{}^+\mathrm{}^{}`$ decay, it is enough to use CPT invariance and make the replacements in Eq. (8)

$`𝒱_{L_1;R_1}`$ $``$ $`\overline{𝒱}_{L_1;R_1},`$
$`𝒱_{L_i;R_i}`$ $``$ $`\overline{𝒱}_{L_i;R_i},(i=2,3,4),`$
$`𝒮_{L;R}`$ $``$ $`\overline{𝒮}_{R;L},`$ (9)
$`B_{6;7;8}`$ $``$ $`\overline{B}_{6;7;8},`$
$`B_{9;10;11}`$ $``$ $`\overline{B}_{9;10;11},`$
where the meaning of the bar can be explained as follows: If any arbitrary one of the coefficients $`𝒱_{L_i;R_i}`$ is represented in the form

$`𝒱_{L_i;R_i}=\left|𝒱_{L_i;R_i}\right|e^{i\delta _{L_i;R_i}}e^{i\varphi _{L_i;R_i}},`$
where $`\delta `$ and $`\varphi `$ are the strong and the weak phases, respectively, then

$`\overline{𝒱}_{L_i;R_i}=\left|𝒱_{L_i;R_i}\right|e^{i\delta _{L_i;R_i}}e^{i\varphi _{L_i;R_i}}.`$
The starting point in analysis of the time evolution spectrum of the $`B^0(\overline{B}^0)\pi ^+\pi ^{}\mathrm{}^+\mathrm{}^{}`$ decay is Eq. (8). The time evolution of the $`B`$ mesons is given by

$`|B^0(t)`$ $`=`$ $`g_+(t)|B^0+{\displaystyle \frac{1}{\xi }}g_{}(t)|\overline{B}^0,`$
$`|\overline{B}^0(t)`$ $`=`$ $`g_+(t)|\overline{B}^0+\xi g_{}(t)|B^0,`$
with

$`g_\pm =exp\left[\left({\displaystyle \frac{\mathrm{\Gamma }_B}{2}}im_B\right)t\right]\{\mathrm{cos}\left({\displaystyle \frac{\mathrm{\Delta }m}{2}}t\right);\mathrm{sin}\left({\displaystyle \frac{\mathrm{\Delta }m}{2}}t\right)\},`$
and $`\xi =p/q`$. Here $`m_B(\mathrm{\Gamma }_B)`$ and $`\mathrm{\Delta }m_B(\mathrm{\Delta }\mathrm{\Gamma }_B)`$ are the average and the difference of the masses (widths) of the two mass eigenstates $`B_H`$ and $`B_L`$, respectively.

The matrix element of the time dependent $`B^0(t)\pi ^+\pi ^{}\mathrm{}^+\mathrm{}^{}`$ decay can be obtained from Eq. (8) by making the following replacements.

$`𝒱_{L_1;R_1}𝒱_{L_1;R_1}(t)`$ $`=`$ $`g_+(t)𝒱_{L_1;R_1}{\displaystyle \frac{1}{\xi }}g_{}(t)\overline{𝒱}_{L_1;R_1},`$
$`𝒱_{L_i;R_i}𝒱_{L_i;R_i}(t)`$ $`=`$ $`g_+(t)𝒱_{L_i;R_i}+{\displaystyle \frac{1}{\xi }}g_{}(t)\overline{𝒱}_{L_i;R_i},(i=2,3,4),`$
$`𝒮_{L;R}𝒮_{L;R}(t)`$ $`=`$ $`g_+(t)𝒮_{L;R}+{\displaystyle \frac{1}{\xi }}g_{}(t)\overline{𝒮}_{R;L},`$ (10)
$`B_{6;7;8}B_{6;7;8}(t)`$ $`=`$ $`g_+(t)B_{6;7;8}+{\displaystyle \frac{1}{\xi }}g_{}(t)\overline{B}_{6;7;8},`$
$`B_{9;10;11}B_{9;10;11}(t)`$ $`=`$ $`g_+(t)B_{9;10;11}{\displaystyle \frac{1}{\xi }}g_{}(t)\overline{B}_{9;10;11}.`$
Similarly the matrix element for the $`\overline{B}^0(t)\pi ^+\pi ^{}\mathrm{}^+\mathrm{}^{}`$ decay can be derived from the respective matrix element $`(\overline{B}^0\pi ^+\pi ^{}\mathrm{}^+\mathrm{}^{})`$ with the help of the replacements

$`\overline{𝒱}_{L_1;R_1}\overline{𝒱}_{L_1;R_1}(t)`$ $`=`$ $`g_+(t)\overline{𝒱}_{L_1;R_1}\xi g_{}(t)𝒱_{L_1;R_1},`$
$`\overline{𝒱}_{L_i;R_i}\overline{𝒱}_{L_i;R_i}(t)`$ $`=`$ $`g_+(t)\overline{𝒱}_{L_i;R_i}+\xi g_{}(t)𝒱_{L_i;R_i},(i=2,3,4),`$
$`\overline{𝒮}_{L;R}\overline{𝒮}_{L;R}(t)`$ $`=`$ $`g_+(t)\overline{𝒮}_{L;R}+\xi g_{}(t)𝒮_{R;L},`$ (11)
$`\overline{B}_{6;7;8}\overline{B}_{6;7;8}(t)`$ $`=`$ $`g_+(t)\overline{B}_{6;7;8}+\xi g_{}(t)B_{6;7;8},`$
$`\overline{B}_{9;10;11}\overline{B}_{9;10;11}(t)`$ $`=`$ $`g_+(t)\overline{B}_{9;10;11}\xi g_{}(t)B_{9;10;11}.`$
Using the formalism for the $`K_\mathrm{}_4`$ decay , from the matrix element (8) one can obtain the $`B^0\pi ^+\pi ^{}\mathrm{}^+\mathrm{}^{}`$ decay rate in terms of the following five variables: invariant mass $`s_M=P^2=(p_{\pi ^+}+p_\pi ^{})^2`$ of $`\pi ^+\pi ^{}`$ pair, invariant mass $`s_{\mathrm{}}=q^2=(p_++p_{})^2`$ of $`\mathrm{}^+\mathrm{}^{}`$ pair; the angle $`\theta `$ between $`\stackrel{}{p}_{\pi ^+}`$ and $`\stackrel{}{p}_++\stackrel{}{p}_{}`$, measured with respect to the center of mass of $`\pi ^+\pi ^{}`$ pair; the angle $`\theta _{\mathrm{}}`$ between $`\stackrel{}{p}_+`$ and $`\stackrel{}{p}_{\pi ^+}+\stackrel{}{p}_\pi ^{}`$ , measured with respect to the center of mass of $`\mathrm{}^+\mathrm{}^{}`$ pair and $`\phi `$ is the angle between the normals to the $`\pi ^+\pi ^{}`$ and $`\mathrm{}^+\mathrm{}^{}`$ planes. The final four–body phase volume in terms of the above–mentioned five variables can be written as

$`dX_{PS}={\displaystyle \frac{1}{2(4\pi )^6m_B^2}}X\beta \beta _{\mathrm{}}ds_Mds_{\mathrm{}}d(\mathrm{cos}\theta )d(\mathrm{cos}\theta _{\mathrm{}})d\phi ,`$ (12)
where

$`X`$ $`=`$ $`\left[\left(Pq\right)^2s_Ms_L\right]^{1/2},`$
$`\beta `$ $`=`$ $`{\displaystyle \frac{1}{s_M}}\lambda ^{1/2}(s_M,m_\pi ^2,m_\pi ^2),`$
$`\beta _{\mathrm{}}`$ $`=`$ $`{\displaystyle \frac{1}{s_{\mathrm{}}}}\lambda ^{1/2}(s_{\mathrm{}},m_{\mathrm{}}^2,m_{\mathrm{}}^2).`$
The bounds of the integration variables are

$`4m_\pi ^2`$ $``$ $`s_Mm_B^2,`$
$`4m_{\mathrm{}}^2`$ $``$ $`s_{\mathrm{}}\left(m_B\sqrt{s_M}\right)^2,`$
$`0`$ $``$ $`\theta ,\theta _{\mathrm{}}\pi ,`$
$`0`$ $``$ $`\phi 2\pi .`$
Different scalar quantities appearing in expression of the $`\left|\right|^2`$ can be expressed in terms of the above–mentioned five variables as

$`Pq`$ $`=`$ $`{\displaystyle \frac{1}{2}}(m_B^2s_Ms_{\mathrm{}}),`$
$`QN`$ $`=`$ $`\beta \beta _{\mathrm{}}\left(Pq\mathrm{cos}\theta \mathrm{cos}\theta _{\mathrm{}}\sqrt{s_Ms_{\mathrm{}}}\mathrm{sin}\theta \mathrm{sin}\theta _{\mathrm{}}\mathrm{cos}\phi \right),`$
$`p_Bp_\pm `$ $`=`$ $`{\displaystyle \frac{1}{2}}\left(s_M+Pq\pm X\beta \mathrm{cos}\theta \right),`$
$`qN`$ $`=`$ $`0,`$
$`Qq`$ $`=`$ $`X\beta \mathrm{cos}\theta ,`$
$`QP`$ $`=`$ $`0,`$
$`NP`$ $`=`$ $`X\beta _{\mathrm{}}\mathrm{cos}\theta _{\mathrm{}},`$
$`Q^2`$ $`=`$ $`s_M\beta ^2,`$
$`N^2`$ $`=`$ $`s_{\mathrm{}}\beta _{\mathrm{}}^2,`$
$`ϵ_{\mu \nu \lambda \sigma }Q^\mu P^\nu N^\lambda q^\sigma `$ $`=`$ $`\sqrt{s_Ms_{\mathrm{}}}X\beta \beta _{\mathrm{}}\mathrm{sin}\theta \mathrm{sin}\theta _{\mathrm{}}\mathrm{sin}\phi ,`$
where

$`N=p_+p_{}.`$
The differential decay rate for the $`B^0(t)\pi ^+\pi ^{}\mathrm{}^+\mathrm{}^{}`$ is (in the expression below leptons are taken to be massless. However in calculation of the decay spectrum for $`B^0(t)\pi ^+\pi ^{}\tau ^+\tau ^{}`$, the $`\tau `$ lepton mass is taken into account. But this expression is quite lengthy and for this reason we do not present it in the text)

$`d\mathrm{\Gamma }`$ $`=`$ $`\left|{\displaystyle \frac{G\alpha }{4\sqrt{2}\pi }}V_{tb}V_{td}^{}\right|^2{\displaystyle \frac{g_{\rho \pi \pi }^2}{\left(s_Mm_\rho ^2\right)^2+\mathrm{\Gamma }_\rho ^2m_\rho ^2}}{\displaystyle \frac{1}{2^{14}\pi ^6m_B^3}}X\beta `$ (13)
$`\times ds_Mds_{\mathrm{}}d(\mathrm{cos}\theta )d(\mathrm{cos}\theta _{\mathrm{}})d\phi ,`$
where

$`=_1\mathrm{cos}^2\theta \mathrm{cos}^2\theta _{\mathrm{}}+_2\mathrm{cos}^2\theta \mathrm{cos}\theta _{\mathrm{}}+_3\mathrm{cos}^2\theta +_4\mathrm{sin}^2\theta \mathrm{sin}^2\theta _{\mathrm{}}\mathrm{sin}(2\phi )`$
$`+_5\mathrm{sin}(2\theta )\mathrm{sin}(2\theta _{\mathrm{}})\mathrm{cos}\phi +_6\mathrm{sin}(2\theta )\mathrm{sin}(2\theta _{\mathrm{}})\mathrm{sin}\phi +_7\mathrm{sin}(2\theta )\mathrm{sin}\theta _{\mathrm{}}\mathrm{sin}\phi `$
$`+_8\mathrm{sin}(2\theta )\mathrm{sin}\theta _{\mathrm{}}\mathrm{cos}\phi +_9\mathrm{cos}^2\theta _{\mathrm{}}+_{10}\mathrm{cos}\theta _{\mathrm{}}+_{11}\mathrm{sin}^2\theta \mathrm{sin}^2\theta _{\mathrm{}}\mathrm{cos}^2\phi +_{12},`$
and

$`_1`$ $`=`$ $`\beta ^2\{\lambda ^2\left[\right|𝒱_{L_1}|^2|𝒱_{L_3}|^2+|𝒱_{R_1}|^2\left|𝒱_{R_3}|^2\right]2\lambda (Pq)[2\text{Re}\left(𝒱_{L_2}𝒱_{L_3}^{}\right)`$
$`+`$ $`2\text{Re}\left(𝒱_{R_2}𝒱_{R_3}^{}\right)]4\left(Pq\right)^2[\lambda \left(\right|𝒱_{L_1}|^2+\left|𝒱_{R_1}|^2\right)+|𝒱_{L_2}|^2+\left|𝒱_{R_2}|^2\right]\},`$
$`_2`$ $`=`$ $`2\beta ^2\sqrt{\lambda }\left[\lambda 4\left(Pq\right)^2\right]\left[\text{Re}\left(𝒱_{L_1}𝒱_{L_2}^{}\right)+\text{Re}\left(𝒱_{R_1}𝒱_{R_2}^{}\right)\right],`$
$`_3`$ $`=`$ $`\beta ^2\lambda \left\{\lambda \left(\left|𝒱_{L_3}\right|^2+\left|𝒱_{R_3}\right|^2\right)+\left|𝒱_{R_2}\right|^24\left(Pq\right)\left[\text{Re}\left(𝒱_{L_2}𝒱_{L_3}^{}\right)+\text{Re}\left(𝒱_{R_2}𝒱_{R_3}^{}\right)\right]\right\},`$
$`_4`$ $`=`$ $`4\beta ^2\sqrt{\lambda }\left[\text{Im}\left(𝒱_{L_2}𝒱_{L_1}^{}\right)+\text{Im}\left(𝒱_{R_2}𝒱_{L_1}^{}\right)\right],`$
$`_5`$ $`=`$ $`\beta ^2\sqrt{s_Ms_{\mathrm{}}}\left\{\lambda \left[\text{Re}\left(𝒱_{L_2}𝒱_{L_3}^{}\right)+\text{Re}\left(𝒱_{R_2}𝒱_{R_3}^{}\right)\right]2\left(Pq\right)\left(\left|𝒱_{L_2}\right|^2+\left|𝒱_{R_2}\right|^2\right)\right\},`$
$`_6`$ $`=`$ $`\beta ^2\sqrt{\lambda }\{\lambda [\text{Im}\left(𝒱_{L_3}𝒱_{L_1}^{}\right)+\text{Im}\left(𝒱_{R_3}𝒱_{R_1}^{}\right)]+2\left(Pq\right)[\text{Im}\left(𝒱_{L_1}𝒱_{L_2}^{}\right)`$
$`+`$ $`\text{Im}\left(𝒱_{R_1}𝒱_{R_2}^{}\right)]\},`$
$`_7`$ $`=`$ $`2\beta ^2\lambda \sqrt{s_Ms_{\mathrm{}}}\left[\text{Im}\left(𝒱_{L_3}𝒱_{L_2}^{}\right)\text{Im}\left(𝒱_{R_3}𝒱_{R_2}^{}\right)\right],`$
$`_8`$ $`=`$ $`2\beta ^2\sqrt{\lambda s_Ms_{\mathrm{}}}\{\lambda [\text{Re}\left(𝒱_{L_1}𝒱_{L_3}^{}\right)\text{Re}\left(𝒱_{R_1}𝒱_{R_3}^{}\right)]+2\left(Pq\right)[\text{Re}\left(𝒱_{L_1}𝒱_{L_2}^{}\right)`$
$`+`$ $`\text{Re}\left(𝒱_{R_1}𝒱_{R_2}^{}\right)]\},`$
$`_9`$ $`=`$ $`4\beta ^2s_{\mathrm{}}s_M\lambda \left[\left|𝒱_{L_1}\right|^2+\left|𝒱_{R_1}\right|^2\right],`$
$`_{10}`$ $`=`$ $`8\beta ^2\sqrt{\lambda }\left[\text{Re}\left(𝒱_{L_1}𝒱_{L_2}^{}\right)+\text{Re}\left(𝒱_{R_1}𝒱_{R_2}^{}\right)\right],`$
$`_{11}`$ $`=`$ $`4\beta ^2s_{\mathrm{}}s_M\left[\lambda \left(\left|𝒱_{L_1}\right|^2+\left|𝒱_{R_1}\right|^2\right)\left|𝒱_{L_2}\right|^2\left|𝒱_{R_2}\right|^2\right],`$
$`_{12}`$ $`=`$ $`4\beta ^2s_{\mathrm{}}s_M\left[\left|𝒱_{L_2}\right|^2+\left|𝒱_{R_2}\right|^2\right].`$
Taking the narrow resonance limit of $`\rho `$ meson, i.e., using the equations

$`\mathrm{\Gamma }={\displaystyle \frac{g_{\rho \pi \pi }^2m_\rho }{48\pi }}\left(1{\displaystyle \frac{4m_\pi ^2}{m_\rho ^2}}\right)^{3/2},\text{and},`$
$`\underset{\mathrm{\Gamma }0}{lim}{\displaystyle \frac{\mathrm{\Gamma }m_\rho }{\left(s_Mm_\rho ^2\right)^2+m_\rho ^2\mathrm{\Gamma }^2}}=\pi \delta \left(s_Mm_\rho ^2\right),`$
we perform integration over $`s_M`$ easily in Eq. (2). After integrating over $`s_M`$, $`s_{\mathrm{}}`$ and $`\theta `$, we get for the angular distribution of the $`B^0(t)\pi ^+\pi ^{}\mathrm{}^+\mathrm{}^{}`$ decay

$`{\displaystyle \frac{d\mathrm{\Gamma }}{d(cos\theta _{\mathrm{}})d\phi }}=\left|{\displaystyle \frac{G\alpha }{4\sqrt{2}\pi }}V_{tb}V_{td}^{}\right|^2{\displaystyle \frac{1}{2^{14}\pi ^6m_B^3}}{\displaystyle \frac{48\pi ^2}{m_\rho ^2\left(14m_\pi ^2/m_\rho ^2\right)^{1/2}}}`$ (15)
$`\times \{({\displaystyle \frac{2}{3}}_1^{}+2_9^{})\mathrm{cos}^2\theta _{\mathrm{}}+({\displaystyle \frac{2}{3}}_2^{}+2_{10}^{})\mathrm{cos}\theta _{\mathrm{}}+{\displaystyle \frac{4}{3}}_4^{}\mathrm{sin}^2\theta _{\mathrm{}}\mathrm{sin}(2\phi )`$
$`+{\displaystyle \frac{2}{3}}_{11}^{}\mathrm{sin}^2\theta _{\mathrm{}}[1+\mathrm{cos}(2\phi )]+({\displaystyle \frac{2}{3}}_3^{}+2_{12}^{})\},`$
where we have introduced the notation $`_i^{}X_i𝑑s_{\mathrm{}}`$. CP violation manifests itself in the coefficients $`_4^{}`$ and $`_{11}^{}`$. The dominant CP violation term is contained in the coefficient $`(4/3)_4^{}`$. After integrating over $`\mathrm{cos}\theta _{\mathrm{}}`$ we obtain the $`\phi `$–distribution

$`{\displaystyle \frac{d\mathrm{\Gamma }}{d\phi }}{\displaystyle \frac{2}{3}}\left({\displaystyle \frac{2}{3}}_1^{}+2_9^{}\right)+{\displaystyle \frac{16}{9}}_4^{}\mathrm{sin}(2\phi )+{\displaystyle \frac{8}{9}}_{11}^{}\left[1+\mathrm{cos}(2\phi )\right]+2\left({\displaystyle \frac{2}{3}}_3^{}+2_{12}^{}\right).`$
As has already been noted the term $`\mathrm{sin}(2\phi )`$ is CP– and T–odd. Hence this term can be isolated by considering the following asymmetry in terms of the differential decay rates of the $`B`$ meson with respect to $`\phi `$

$`𝒜_\phi (t)={\displaystyle \frac{\left\{{\displaystyle _0^{\pi /2}}{\displaystyle _{\pi /2}^\pi }+{\displaystyle _\pi ^{3\pi /2}}{\displaystyle _{3\pi /2}^{2\pi }}\right\}{\displaystyle \frac{d\mathrm{\Gamma }}{d\phi }}d\phi }{\left\{{\displaystyle _0^{\pi /2}}+{\displaystyle _{\pi /2}^\pi }+{\displaystyle _\pi ^{3\pi /2}}+{\displaystyle _{3\pi /2}^{2\pi }}\right\}{\displaystyle \frac{d\mathrm{\Gamma }}{d\phi }}d\phi }}.`$
The main feature of such CP violating asymmetries is that they can be obtained by considering the sum of the differential decay rates of the $`B`$ and $`\overline{B}`$ rather than the difference of these rates, as usual.

Furthermore we study the dependence of the asymmetry between the spectrum integrated rates of $`B^0(t)\pi ^+\pi ^{}\mathrm{}^+\mathrm{}^{}`$ and $`\overline{B}^0(t)\pi ^+\pi ^{}\mathrm{}^+\mathrm{}^{}`$ decays, on the new Wilson coefficients, which is defined as

$`𝒜_1(t)={\displaystyle \frac{\mathrm{\Gamma }(B^0(t)\pi ^+\pi ^{}\mathrm{}^+\mathrm{}^{})\mathrm{\Gamma }(\overline{B}^0(t)\pi ^+\pi ^{}\mathrm{}^+\mathrm{}^{})}{\mathrm{\Gamma }(B^0(t)\pi ^+\pi ^{}\mathrm{}^+\mathrm{}^{})+\mathrm{\Gamma }(\overline{B}^0(t)\pi ^+\pi ^{}\mathrm{}^+\mathrm{}^{})}}.`$ (16)
It should noted that it is also possible to construct a differential asymmetry that isolates another combination of the terms containing imaginary parts (see for example ). Such terms, for example, can be isolated by considering the difference distribution of the same hemisphere and opposite hemisphere events. The corresponding asymmetry can be defined then as follows

$`𝒜_2(t)`$ $`=`$ $`{\displaystyle \frac{1}{\mathrm{\Gamma }}}\left\{{\displaystyle _0^\pi }{\displaystyle _0^{2\pi }}\right\}d\phi \left\{{\displaystyle _1^0}{\displaystyle _0^{+1}}\right\}d(\mathrm{cos}\theta _{\mathrm{}})\left\{{\displaystyle _1^0}{\displaystyle _0^{+1}}\right\}d(\mathrm{cos}\theta )`$ (17)
$`\times `$ $`{\displaystyle \frac{d\mathrm{\Gamma }}{d\phi d(\mathrm{cos}\theta _{\mathrm{}})d(\mathrm{cos}\theta )}}.`$
## 3 Numerical analysis

In this section we study the dependence of the asymmetries $`𝒜_\phi ,𝒜_1`$ and $`𝒜_2`$ on new Wilson coefficients and time for the $`B^0(t)\pi ^+\pi ^{}\mathrm{}^+\mathrm{}^{}`$ decay. Before presenting the numerical results we want to note that, although the expressions for the asymmetries $`𝒜_\phi ,𝒜_1`$ and $`𝒜_2`$ are given for the most general case, i.e., each one of new Wilson coefficients might have new strong and weak phases. In our numerical analysis we choose all new Wilson coefficients to be real and they all vary in the region between $`4`$ and $`+4`$.

For the form factors which are needed in the course of performing numerical calculations, we have used the prediction of the light cone QCD sum rules. In our numerical analysis we will use the results of the work  (see also ) in which the form factors are described by a three–parameter fit where the radiative corrections up to leading twist contribution and SU(3)–breaking effects are taken into account. The $`q^2`$–dependence of the form factors which appear in our analysis could be parametrized as

$`F(s)={\displaystyle \frac{F(0)}{1a_Fs+b_Fs^2}},`$
where $`s=q^2/m_B^2`$ is the dilepton invariant mass in units of $`B`$ meson mass, and the parameters $`F(0)`$, $`a_F`$ and $`b_F`$ are listed in Table 1 for each form factor.

As it comes to the values of the Wilson coefficients $`C_7^{eff}(m_b)`$ and $`C_{10}(m_b)`$, whose analytical expressions in the standard model are given in , are strictly real as can be read off from Table 2. In the leading logarithmic approximation, at the scale $`𝒪(\mu =m_b)`$, we have

$`C_7^{eff}(m_b)`$ $`=`$ $`0.313,`$
$`C_{10}^{eff}(m_b)`$ $`=`$ $`4.669.`$
Although individual Wilson coefficients at $`\mu m_b`$ level are all real (see Table 2), the effective Wilson coefficient $`C_9^{eff}(m_b,\widehat{s})`$ has a finite phase, and in next–to–leading order

$`C_9^{eff}(m_b,\widehat{s})=C_9(m_b)\left[1+{\displaystyle \frac{\alpha _s(\mu )}{\pi }}\omega (\widehat{s})\right]+Y_{SD}(m_b,\widehat{s})+Y_{LD}(m_B,s),`$ (18)
where $`C_9(m_b)=4.344`$. Here $`\omega \left(\widehat{s}\right)`$ represents the $`𝒪(\alpha _s)`$ corrections coming from one–gluon exchange in the matrix element of the corresponding operator, whose explicit form can be found in . In (18) $`Y_{SD}`$ and $`Y_{LD}`$ represent, respectively, the short– and long–distance contributions of the four–quark operators $`𝒪_{i=1,\mathrm{},6}`$ . Here $`Y_{SD}`$ can be obtained by a perturbative calculation

$`Y_{SD}(m_b,\widehat{s})`$ $`=`$ $`g(\widehat{m}_c,\widehat{s})C^{(0)}{\displaystyle \frac{1}{2}}g(1,\widehat{s})\left[4C_3+4C_4+3C_5+C_6\right]`$
$``$ $`{\displaystyle \frac{1}{2}}g(0,\widehat{s})\left[C_3+3C_4\right]+{\displaystyle \frac{2}{9}}\left[3C_3+C_4+3C_5+C_6\right]`$
$``$ $`\lambda _u\left[3C_1+C_2\right]\left[g(0,\widehat{s})g(\widehat{m}_c,\widehat{s})\right],`$
where

$`C^{(0)}`$ $`=`$ $`3C_1+C_2+3C_3+C_4+3C_5+C_6,`$
$`\lambda _u`$ $`=`$ $`{\displaystyle \frac{V_{ub}V_{ud}^{}}{V_{tb}V_{td}^{}}},`$
and the loop function $`g(m_q,s)`$ stands for the loops of quarks with mass $`m_q`$ at the dilepton invariant mass $`s`$. This function develops absorptive parts for dilepton energies $`s=4m_q^2`$:

$`g(\widehat{m}_q,\widehat{s})={\displaystyle \frac{8}{9}}\mathrm{ln}\widehat{m}_q+{\displaystyle \frac{8}{27}}+{\displaystyle \frac{4}{9}}y_q{\displaystyle \frac{2}{9}}\left(2+y_q\right)\sqrt{\left|1y_q\right|}`$
$`\times \left[\mathrm{\Theta }(1y_q)\left(\mathrm{ln}{\displaystyle \frac{1+\sqrt{1y_q}}{1\sqrt{1y_q}}}i\pi \right)+\mathrm{\Theta }(y_q1)\mathrm{\hspace{0.17em}2}\mathrm{arctan}{\displaystyle \frac{1}{\sqrt{y_q1}}}\right],`$
where $`\widehat{m}_q=m_q/m_b`$ and $`y_q=4\widehat{m}_q^2/\widehat{s}`$. In addition to these perturbative contributions $`\overline{c}c`$ loops can excite low–lying charmonium states $`\psi (1s),\mathrm{},\psi (6s)`$ whose contributions are represented by $`Y_{LD}`$ . However in the present work we restrict ourselves to the consideration of short distance contributions only and therefore we find it redundant to display the explicit form of $`Y_{LD}`$.

Calculations show that the asymmetry $`𝒜_\phi `$ for the $`B^0(t)\pi ^+\pi ^{}\mu ^+\mu ^{}`$ decay is sensitive only to tensor type interaction. We observe that, practically, the existence of the new vector or scalar type interactions do not effect the SM results crucially, which is about $`13\%`$. However when tensor type interactions are taken into consideration this asymmetry is reduced noticeably and $`\left(𝒜_\phi \right)_{max}2.5\%`$. Therefore if in future experiments the results obtained results differ from the SM results, it is an unambiguous indication of the existence of new physics beyond SM.

As one considers the $`B^0(t)\pi ^+\pi ^{}\tau ^+\tau ^{}`$ decay, the maximum value of the asymmetry $`𝒜_\phi `$ is about $`5\%`$. Similarly, as in the $`B^0(t)\pi ^+\pi ^{}\mu ^+\mu ^{}`$ decay, $`𝒜_\phi `$ is not sensitive to the new vector and scalar interactions and for the contribution of the tensor interaction is less $`1\%`$.

We further studied the dependence of the asymmetry $`𝒜_1`$ on the new Wilson coefficients and time for both $`B^0(t)\pi ^+\pi ^{}\mu ^+\mu ^{}`$ and $`B^0(t)\pi ^+\pi ^{}\tau ^+\tau ^{}`$ decays. Our calculations show that the maximal possible value that $`𝒜_1`$ asymmetry gets for the $`B^0(t)\pi ^+\pi ^{}\mu ^+\mu ^{}`$ decay is about $`2\%`$ for all type of new interactions. Similar behavior takes place for the $`B^0(t)\pi ^+\pi ^{}\tau ^+\tau ^{}`$ case as well, with only one exception. For tensor type interactions which are controlled by the new Wilson coefficients $`C_T`$ and $`C_{TE}`$, $`𝒜_1`$ can arrive at large values about $`20\%`$ (see Figs. (1) and (2)). So any departure from the SM prediction in future experiments, is an indication of the presence of tensor type interactions.

Finally we have studied the dependence of the $`𝒜_2`$ asymmetry on new Wilson coefficients and time for the above mentioned decays. For the $`B^0(t)\pi ^+\pi ^{}\mu ^+\mu ^{}`$ decay the maximal value of $`𝒜_2`$ is about coefficient $`C_{LL}=4`$ causes $`𝒜_2`$ to decrease about $`50\%`$ from the SM prediction, but if $`C_{LL}=+4`$, practically the value of $`𝒜_2`$ seems to coincide with the SM prediction, which also gives us a clue in determining the sign of of the new vector interaction, as can be seen in Fig. (3). As is depicted in Fig. (4), similar situation seems to hold for the new Wilson coefficient $`C_{RL}`$. Our numerical analysis shows that $`𝒜_2`$ is quite insensitive to the existence of scalar interaction and yields almost same results with SM prediction. But an important observation is that the tensor interaction enhances the SM prediction by about $`50\%`$. Therefore it would be an unambiguous confirmation of the existence of new vector (tensor) interaction, if in future experiments the measured value of $`𝒜_2`$ appears to be less (larger) than the SM prediction.

On the other hand for all values of the new Wilson coefficients, the $`𝒜_2`$ asymmetry for the $`B^0(t)\pi ^+\pi ^{}\tau ^+\tau ^{}`$ case coincides almost with the SM prediction. For this reason $`𝒜_2`$ for $`B^0(t)\pi ^+\pi ^{}\tau ^+\tau ^{}`$ decay does not seem to provide us with any new information about new physics beyond SM.

As a final remark we should emphasize that similar analysis can be performed for the $`\overline{B}^0(t)\pi ^+\pi ^{}\mathrm{}^+\mathrm{}^{}`$ decay as well. For this aim it is enough to use Eqs. (8), (2) and (2).

In conclusion, in this work we have studied the time evolution of the decay spectrum for the $`B^0(t)\pi ^+\pi ^{}\mathrm{}^+\mathrm{}^{}`$ decay. The sensitivity of the experimentally measurable asymmetries $`𝒜_\phi `$, $`𝒜_1`$ and $`𝒜_2`$ to the new Wilson coefficients and time is studied in detail. It is observed that different asymmetries show different dependencies on different new Wilson coefficients for the $`B^0(t)\pi ^+\pi ^{}\mu ^+\mu ^{}`$ and $`B^0(t)\pi ^+\pi ^{}\tau ^+\tau ^{}`$ channels. Studying different asymmetries for the above–mentioned two channels on new Wilson coefficients can give essential information about new physics and can serve as an efficient tool in determining not only of the magnitude, but also of the sign of new Wilson coefficients.

## Figure captions

Fig. (1) The dependence of the asymmetry $`𝒜_1`$ on the new Wilson coefficient $`C_T`$ and time for $`B^0(t)\pi ^+\pi ^{}\tau ^+\tau ^{}`$ decay, where $`\tau _s`$ is the $`B^0`$ meson life time.
Fig. (2) The same as in Fig. (1), but for the new Wilson coefficient $`C_{TE}`$.
Fig. (3) The dependence of the asymmetry $`𝒜_2`$ on the new Wilson coefficient $`C_{LL}`$ and time for $`B^0(t)\pi ^+\pi ^{}\mu ^+\mu ^{}`$ decay.
Fig. (4) The same as in Fig. (3), but for the new Wilson coefficient $`C_{RL}`$.