Skip to main content

Full text of "Axial-vector mesons from $τ\to APν_τ$ decays"

See other formats


Axial-vector mesons from r — APv^ decays 



G. CalderoiE 

Facultad de Ingenieria Mecdnica y Electrica, Universidad Autonoma de Coahuila, CP. 27276, Torreon, Coahuila, Mexico 

J. H. Muiio^ and C. E. Ver4!l 

Departamento de Fisica, Universidad del Tolima A. A. 546, Ibague, Colombia 

Axial-vector mesons ai(1260), /i(1285), /ii(1170), ifi(1270) and ii'i(1400) can be produced 
in semileptonic r APvt decays, where P stands for the pseudoscalar mesons tt or K. We 
calculate the branching ratios based in a meson dominance model. The exclusive channels 
t" ai(1260)"7r''i/, t~ -i> ai(1260)°7r"!/, r~ /i(1285)7r"!y and ^ /!.i(1170)7r"i/ turn out to 
beof order 0(10-*) and channels r" ^ is'i(1270)"7r°//, ^ Ki{127Qf-K~v, -> J<i(1400)"7rV 
and T~ — >■ Js'i(1400)''7r~z/ of order 0(10"^). These results indicate that the branching ratios could 
be measured in experiments. 

PACS numbers: 12.40Vv, 13.35Dx 



* |german.caIderon@ uadec.edu.mx| 
t jhmunoz@ut.edu. co] 
■f .cvera@ut.edu. CO 



2 



I. INTRODUCTION 

Recently the BABAR Collaboration [l[ performed measurements of the branching fractions of three-prong and 
five-prong t decay modes. Specifically, the branching fraction for the — > /i(f 285)7r^i/7- decay is measured 

S{t- /i(1285)7r-zy^) = (3.79 ± 0.25) x lO^'^, (1) 

where the axial- vector meson /i(1285) is detected using the /i 27r+27r~ and /i -> 7r+7r~77 modes. We have 
averaged in quadrature the two branching ratios and added the statistical and systematic errors reported by BABAR 
Collaboration [i| . This result supersedes a value reported by BABAR Collaboration 0] . Besides CLEO Collaboration 
Q published a branching ratio of (S.S^iJ'g ± 1.8) x 10~^, which is also based on the /i — 27r+27r~ mode. 

The measurements of t — APv decays can give information about the hadronic matrix elements {AP\J^\0) at 
the moderate energy regime of t decays. This can be relevant in nonleptonic two-body D and B decays, where the 
annihilation contribution plays a relevant role in the branching ratios and CP asymmetries, see Refs. \^]^. 

The r -> APi> decays can contribute to the r decay modes with four-pseudoscalar mesons in the final state. 
The exclusive — > /i(1285)7r~i/ decay contributes as intermediate state in r decay involving five-pseudoscalar 
mesons. The former decays can be related to e+e" — >■ (47r)°. The relation between the t~ — )■ {A'k)~v decays and the 
electromagnetic annihilation e+e^ — (47r)° process appears because the components of the same current are present 
in both cases, i.e., the current (i7pM in r decay and l/\/2{u'^^u — d^y^d) in the case of e^e^ annihilation. Thus, these 
two processes are related by isospin symmetry. 

On the other hand, at theoretical level, the production of p-wave mesons in semileptonic r decays has not been 
widely studied in the literature. In effect, Ref. 8] and this paper are the only works that have considered the 
production of orbitally excited mesons in r — MPvr channels, where M can be an axial-vector meson or a tensor 
meson. The exception has been Ref. @, where the branching ratio for the channel — > /i(1285)7r^j/ estimated is 
2.91 X 10"''. 

In this paper we estimate the branching ratios of the t — > APv decays, where P is the pseudoscalar meson tt or 
K and A is the axial-vector meson ai(1260), /i(1285), ft.i(1170), Ki{l27Q) or ^1(1400). We have considered the 
kinematically allowed, the G-parity conserved and first-class currents decays. We use a meson dominance model, 
where the resonances that contribute are coming from the vector mesons p(770), p(1450) and i^*(982) and the axial- 
vector meson ai(1260). The weak decay and strong coupling constants are determined from experimental data. We 
use flavor SU{^) symmetry where there are no data to estimate the constants. 

In our estimations, we find that the branching ratios for the exclusive channels — > ai(1260)^7r°j/, — > 
ai(1260)°7r-;/, r" /i(1285)7r-i/ and t" /ii(1170)7r-j/ are of order 0(10-^). The channels r Ki{127Q)ttu 
and r -J> ifi(1400)7ri^ decays are of order 0(10"^). Finally, the modes t ai{\2m)Kv and r Ki{\210)Kv are 
kinematically suppressed. Up to now, the — ^ /i(1285)7r~i^ decay is the only mode that can be compared with 
experimental measurements. 

The paper is organized as follow. In Sec. II, we describe the amplitude in the meson dominance model and give the 
expression to estimate the width decay. The calculation of the weak decay and strong coupling constants are given 
in Sec. III. In Sec. IV, we present the estimation of the branching ratios for the modes considered in this paper. And 
finally, we conclude in Sec. V. 

II. AMPLITUDE AND WIDTH DECAY FOR SEMILEPTONIC r ^ APv DECAYS 

The decay amplitude for the T{pr) — >■ A{pa)P{Pp)i't{Pu) process, where A denotes the axial-vector meson ai(1260), 
/i(1285), /ii(1170), Ari(1270) or i^i(1400), and P is the pseudoscalar meson tt or AT, can be writing by 



G 



M{t^ APvr) = -^Vu.u{pr)r{l - 75 )w(p,) ( A(pa)P(p) | J^O) 1 0) , (2) 

where Gf is the Fermi constant, Vui or s) is the corresponding element of the Cabibbo-Kobayashi-Maskawa 

matrix and .//^(O) is the {V — A) weak current. 

The hadronic matrix element {A{pA,s)P{p)\Jfi\Q) is parametrized by 

{A{pA,e)P{p)\J^\Q) ^l{t)e^ + {e.p)[c+{t){p-pA),, + c_{t){p + pA)^,] 

+ ig(<)e^.p.e"(pA -pY{pa+pT , (3) 



3 



where the hadronic matrix element is expressed in terms of the form factors l{t), c+(t), c_(i) and q{t), which are 
Lorentz- invariant functions of squared momentum transfer t — (p^ the polarization vector describes the 

axial- vector meson A. 

To estimate the form factors, we use the meson dominance model, i.e., the decay amplitude is given by the sum of 
intermediate meson resonance contributions 



M{t ^ APi^r) ^"^Mir ^ R* fr ^ APi^r) , (4) 
i?* 

where R* is the intermediate resonance, which has the right quantum numbers for the specific processes. The sum is 
extended over all possible resonance contributions, see Fig. 1. 

The production of the resonance R* can be classified as either first- or second-class current depending on the spin J, 
parity P, and G-parity of the resonance particle. In the Standard Model, first-class current is considered to dominate. 
The second-class currents are associated with a decay constant proportional to the mass difference between an up and 
down quark and in the exact isospin limit symmetry they vanish The vector meson resonances p(770), p(1450), 
-ftr*(982) and the axial- vector meson ai(1260) are produced by first-class currents, which are the contributions included 
in this work. The vertex R* — > AP decay, which is produced by strong interaction, must conserve G-parity. 

Comparing the hadronic current, see Eq. ([3]), with the amplitude build from the meson dominance model, which 
is calculated using Feynman rules, we determine expressions for the form factors q(t) = and c_(t) = — c+(t) = c{t). 
Thus, the amplitude is expressed in terms of only two form factors l{t) and c{t). These form factors can be expressed 
by means of the Breit-Wigner (BW) function 



where BWv{t) is defined by 



c{t) 



M 
ifvgAVP 
2M 



BWvit) 



BWvit), 



AP 



IVP -t- iMT ' 



(5) 



(6) 



the parameters M and T are the mass and width decay of the intermediate meson resonant particle, where we have 
used the complex-mass prescription Mg — iMT and a constant width F. 

The unpolarized squared amplitude in terms of the form factors lit) and c(t) is the following 



pol 

where the kinematic factors are defined by 



aiit, s)|/(t)P -t- a2(t, s)|c(i)P + asit, s)Re[l*it)cit)] 



(7) 



aiit,s) 


1 - 


m\ . 


a2it,s) 


1 - 


m\ . 


a3it,s) 


to 





is - ml)is + t~ml- m^) + m^(TO^ + 2to^ - 2t - s) 
iim^ - t)^ - 2(m^ + t)m\ + m,\)ii'rr? - s)m\ + (s - rnl)is + t - ml - m^)) 



im^ — s)m\ + (m^ — i)(s — m^)(m^ -I- 



I 2 / 2 I 2.1 2 . 4 2 2 4\ 

+ mj^(s + m r + m^t — m — to to^ — m^) 



t) 



(8) 



The kinematics factors ai(<,s), a2(t, s) and a^it,s,) are functions of squared moments t = ipA — p)"^ and s ~ 
iVv and the masses of r lepton to^ , axial-vector meson toa and pseudo-scalar meson m. 

The width decay for the channel t --)■ APur is calculated integrating the unpolarized amplitude, Eq. ([T]), with 
respect to the variables t and s, and the appropriate kinematic limits. 

The vector meson resonant contribution is the only contribution for the channels studies in this work, with the ex- 
ception of the T~ — /i(1285)7r~j/ decay, where it has an axial-vector meson ai(1260) as the only present contribution. 



4 




FIG. 1. Intermediate meson dominance resonance of r -> APv decays. 

Specifically, for the channels r ai(1260)7rz^, r /ii(1170)7ri^, r Ki{l21{))Kv and r ^ i^i(1400)A'i' we have the 
contributions of vector meson p(770) and p(1450) resonances. The strange vector meson i4r*(982) contributes to the 
processes r -J> ai{l2m)Kv, t ^ Ki{1270)TTiy and r ^ /Ci(1400)7ri/. 

In studies of the ^ tt^tt^z^ channel carried out by the Belle [l^ and ALEPH [l^ Collaborations, was established 
that in order to describe correctly the spectral function it is necessary to introduce two resonances, the vector mesons 
p(770) and p(1450). We model the two contributions with a linear combination normalized by 

where the parameter /3 is determined by a fit to the hadronic spectral function obtained by the experiments. In 
our work, we model in the same manner the contribution of the two resonances p(770) and p(1450) to the channels 
considered above. We use a conservative value of /3 = ±(0.2 ± 0.1), see Ref. Q. 

III. WEAK DECAY AND STRONG COUPLING CONSTANTS 

In order to obtain the weak decay and strong coupling constants that are required in this work, we make use of the 
available experimental information. However, for some channels, we use flavor SU{S) and isospin SU{2) symmetries 
to relate them to the constants obtained by experimental data. 

We use the experimental reported branching ratios in Ref. [14] to determine the strong coupling constants from 
the decays of strange axial- vector mesons ii:i(1270) Kp, Xi(1270) ^ K*{982)n and i4:i(1400) K*{982)t: and 
from the /ii(1170) — >■ pir decay. 

The strong coupling constant qavp (in GeV^^) is determined from the amplitude of probability for the A — > VP 
process, which is defined by 

M{A{pA,eA) ^- V{pv,ev)P{p)) gAVP£A'^v[iPA-Pv)9p.^ - {pa)ij.{pv),^] ■ (10) 
The decay rate in the rest frame of the decaying axial- vector meson is 

T{A VP) = -— ^- [2mAmy + {1x1^ + my - m ) \ , (11) 

where A(x, y, z) is the Kallen function, which is defined by \{x, y, z) — + y"^ + — 2{xy + xz + yz). 

Using the measured width decay for the channels ^1(1270) Kp, ivri(1270) K*Tr, iCi(1400) ^ Kp, 
-ft'i(1400) K*Tr, and using the isospin SU{2) symmetry to relate with the charged modes, we obtain the fol- 
lowing values in (GeV^): 9x^(1270)- ROp- = 5.04± 0.67, 5Ki(i270)OK-p- = 5.04± 0.67, 5ifj(i27o)-K*-770 =0.69±0.13, 
9Ki{i27a)0K*-Ti- = 0.97 ± 0.19, gKi{i4Q0)-K0p- = 0.68 ± 0.34, 9Ki{i40Q)0K-p- = 0.68 ± 0.34, 9ki{i4oo)-k*-ttO = 
1.95 ±0.10 and gKi(iAoa)"K-'-TT- = 2.75 ±0.14. The error reported is due to the experimental uncertainty in the width 
decay. 

In order to extract the strong decay constant 9hipTn we have assumed that only the mode /ii(1170) — !■ pn contributes 
to the total width decay of the axial- vector meson /ii(1170). We obtain 9hip-7T+ = 3.85±0.21 GeV^^. The uncertainty 
is determined by the experimental error in the width decay. 



5 



There are two types of axial-vector mesons, which in spectroscopy notation ^*+^Lj are p-wave mesons ^Pi and "^Pi. 
Considering the flavor SU{3) symmetry, the -^Pi-nonet is integrated by ai(1260), /i(1285), /i(1420) and Kia and 
the ^Pi-nonet composed by 6i(1235), /ii(1170), /ii(1380) and Kib- Here, in this work is relevant consider that the 
physical strange axial- vector mesons ifi (1270) and (1400) are the result of a mixing of the flavor SU{3) eigenstates 
KiA and Kib UM- Thus, the mixing can be described by 

ii'i(1270) = KiA sin + Kib cos 6 

Ki{UOO) ^ KiACOse - KiBsine , (12) 

where 9 is the mixing angle. Restricting its value to < 61 < 90°, it is found two solutions 9 = 32° and 9 = 58°, see 
Ref. 0. 

We consider flavor SU{3) symmetry to relate the strong coupling constants gaip-w and gaiK'K with QkiapK and 
gxiAK'^, respectively. The relations are ga^p-nr = \/^gKiAK*7v and ga^K'K = ^/TJigKtApK- 

Furthermore, we use isospin SU(2) symmetry to relate the strong coupling constants for different charged states 
in the respective channels. For example, in the case of the axial- vector meson i<ri(1270), the branching ratio for 
the modes i^i(1270) ^ K*ti and ii:i(1270) ^ Kp are reported in Ref. [ll|. The isospin relations are 9 j^o „- — 

\/V^9KiKp, gK«K-p- = VV^9KiKp, Sk-K'-'k" = ^/VSgKiK'TT and gKOK—n- = \/V59KiK*tt- Similar analysis 
for strong coupling constant in decay modes for the axial- vector meson isri(1400) is done. The results obtained are: 
9a: p~.o= 5.75 ±0.7A, Kp--- =3.75 ±0.74, 5^^-^..^o = 1.17 ± 0.09, and gaOK-K- =1.64 ±0.13, in unit GeY-\ 
The reported error is coming from the experimental uncertainty in the width decays. 

Moreover we use the flavor SU{3) symmetry to obtain the relation gf-^g^-^+ = ^/ V^9a°p-TT- ■ Thus, the strong 
coupling constant has the value 5^^^-^+ = 2.65 ± 0.53 GeV~^. 

We have estimated the weak decay constant of hadron mesons from the experimental data on branching ratio 
from r- h^i/r decays. We obtain fp — 218 ±2 and /k' = 210 ±5 MeV. Also, weak decay constant for the physical 
states ii'i(1270) and 1^1(1400) are obtained. Using a mixing angle 9 = 32° and inverting Eq. ([T^ . we obtain the weak 
decay constant for the flavor Kia and Kib states. The weak decay constant for the axial- vector meson ai(1260) is 
fai = 213 ± 10 MeV, which is obtained by flavor SU{3) symmetry from the weak decay of the mesons Kia and Kib, 
see Ref. [Tsj . In the calculations above, the error in the weak decay constants are coming from the experimental error 
in branching ratios of — > h^Vr decays. 

The necessary values for the relevant parameters such as width decays and masses of the particles involved in the 
channels calculated in this work have been taken from Ref. Urn. 



IV. BRANCHING RATIOS 



The branching ratios for the exclusive r~ — >■ APvr decays estimated in this work are shown in Table I. 

The largest branching ratios estimated in this work turn out to be of order 0(10-^), corresponding to the exclusive 
modes ai(1260)~7r"iy, r~ ai(1260)"7r^i^, fi{1285)TT~u and /ii(1170)7r^i^. In these modes we 

have the contribution of two vector meson resonances p(770) and p(1450), which produce a broad form factors l{t) 
and c{t) as function of the momentum transfered t, in consequence a large branching ratio estimation. 

Up to now, the only measured branching ratio from the calculated in this work is for the channel — > /i(1285)7r-z^. 
Our estimate for this channel is (11.8 ± 4.5) x 10~^. We can compare with the recent result of BABAR Collaboration 
(3.79 ±0.25) X 10-4 dl and determine that our estimation differs by approximately one sigma error with respect to 
the experimental result. 

The channels with branching ratios of order 0(10-^) are r — Ki{127Q)ttv and r — > Ki{1AQQ)tiv decays. They have 
only the contribution of the vector meson if* (982) as intermediate resonance. 

The exclusive channels t~ — !• Ki{1AQQY'K~v and t~ — Ki{\AQQ)~K^v decays are not kinematically allowed. The 
modes r — )■ ai(1260)-firz/ and t Ki{127Q)Kv decays are of the order of 0(10"^), this is due to the small phase 
space available for their decay, i.e, they are kinematically suppressed decays. 

The channels r — > ai(1260)7ri^, — > /ii(1170)7r~i/, r Ki{127Q)ttv and r — >■ ifi(1400)7ri/ decays, with resonance 
contributions coming from the vector mesons /3(770) and p(1450), have the largest error due to the error in the 
parameter /3. In the other cases, where the meson resonance contribution is the strange vector meson if* (982) or 
the axial-vector meson ai(1260), the principal error in the branching ratio estimate is due to the uncertainty in the 
strong couphng constant gAVP, where V is the resonance contribution if*(982) or ai(1260). 

In addition, there are two other sources of error in our estimations not considered. We use flavor SU (3) and 
isospin SU (2) symmetries to relate the weak decay and strong coupling constants to the obtained by experimental 



6 



information. The error in flavor SU{3) symmetry can reach 20% in the estimation. The other source of theoretical 
uncertainty is coming from the contact term in the weak vertex of the hadronic current (AP| J^(0)|0). However, we 
do not consider this contribution because it is associated with the continuum, i.e. there is no BWyit). The contact 
term contribution is not possible to estimate in a meson dominance model. 

TABLE I. Branching ratios for r — >■ APv decays. 



AP mode Branching ratio Parameter /? 



ai(1260)"7r^' 


(20.3 ± 11.7) 


X 


10- 




= 


±0 


2 ± 0.1 


ai(1260)"7r" 


(25.4 ± 18.8) 


X 


10" 


-4 


= 


-0 


2 ± 0.1 


fix I IZDU ) TT 


/'1QQ_|_11 A\ 


X 




-4 


R — 

P — 




9-1-0 1 


ai(1260)"7r" 


(24.5 ± 18.2) 


X 


10" 


-4 


P = 


-0 


2 ±0.1 


ai(1260)"_R'" 


(1.6 ±0.3) 


X 


10" 


-9 








ai(1260)°_ft:' 


(4.6 ±0.6) 


X 


10" 


-9 








/i(1285)7r~ 


(11.8 ±4.5) 


X 


10" 


-4 








/ii(1170)7r- 


(29.1 ± 16.6) 


X 


10" 


-4 




±0 


2 ±0.1 


/ii(1170)7r- 


(43.3 ± 30.8) 


X 


10" 


-4 


13 = 


-0 


2 ±0.1 


/s:i(1270)-7r'^ 


(8.3 ±3.0) 


X 


10" 


-7 








^^1(1270)%- 


(1.6 ±0.6) 


X 


10" 


-6 










(1.3 ±0.2) 


X 


10" 


-6 








7^1(1400)%" 


(2.4 ±0.3) 


X 


10" 


-6 








Ki{1270)°K- 


(11.6 ±5.8) 


X 


10" 


-9 


P = 


±0 


2 ± 0.1 




(2.5 ± 1.3) 


X 


10" 


-9 


P = 


±0 


2 ± 0.1 



V. CONCLUSIONS 

We have estimated branching ratios of the t — > APv decays. The axial- vector mesons A considered are ai(1260), 
/i(1285), /ii(1170), Ki{\27G) and KiilAQQ). The pseudoscalar mesons P are tt and K . In total, 12 exclusive channels 
are kinematically allowed, first class-currents kind and G-parity conserved. 

To calculate the branching ratios we use a meson dominance model, where the form factors of the current {AP\J^\{)) 
are calculated from the vector meson contributions p(770), p(1450) and X*(982). With the exception of the channel 
T~ — )■ /i(1285)7r~z^, where the contribution is due to the axial-vector meson ai(1260). 

The channels t~ ai(1260)"7r°i^, -> ai(1260)"7r~i/, t~ fi{l2S^)iT~v and /ii(1170)7r^i^ decays are of 

order 0(10"''). We can compare the estimated branching ratio Br{T~ — > /i(1285)7r^z^) = (11.8 ±4.5) x 10^'' with 
the measured value (3.79 ± 0.25) x 10"'', recently reported by BABAR Collaboration ^j]. 

We observe a pattern in the branching ratios calculated in this work. The modes of order 0(10^^) are kinematically 
suppressed. The channels of order 0(10"^) are Cabibbo suppressed and they have only the contribution of the vector 
meson _fir*(982). Finally, the channels with branching ratios of order 0(10"''') are Cabibbo allowed and have two 
vector meson resonance contributions jo(770) and p(1450), with the exception of the — )■ /i(1285)7r"i^ where the 
resonance is the axial- vector meson ai(1260). 

The processes studied in this work can contribute as intermediate states in branching ratios with four and five 
pseudoscalar mesons tt and/or K in final states for the r decay modes. 

Eventually, the branching ratios for the channels calculated in this work, with the exception of the channel — >■ 
/i(1285)7r"i/ (which have been already measured), will be measured by four-prong of r decays, with the data sample 
of T lepton pair s accumulated by the B-factories BABAR and Belle experiments [!£)]. Additionally, SuperB [13] and 
Belle-II [isl Il9j| Collaborations have programs to produce t pairs and measure rare r decays. 

ACKNOWLEDGMENTS 

G. C. acknowledges financial support from Promep, Mexico. J. H. M. and C. E. V. are grateful to Comite Central 
de Investigaciones (CCI) of the University of Tolima, Colombia. 



[1] J. P. Lees et al. [BABAR Collaboration], larXiv:1209.2734| : [hep-ex]. 



7 



[2] B. Aubert et al. [BABAR Collaboration], Phys. Rev. D 72, 072001 (2005) hep-ex/050500?. 
[3] T. Bergfeld et al. [CLEO Collaboration], Phys. Rev. Lett. 79, 2406 (1997) liep-ex/9706020\ 
[4] H. -n. Li, C. -D. Lu and F. -S. Yu, Phys. Rev. D 86, 036012 (2012) arXiv: 1203.3120 [hep-ph]]. 
[5] Z. -J. Xiao, W. -F. Wang and Y. -y. Fan, Phys. Rev. D 85, 094003 (2012) arXiv:1111.6264 [hep-ph]]. 
[6] Q. Chang, X. -W. Cui, L. Han and Y. -D. Yang, Phys. Rev. D 86, 054016 (2012) arXiv: 1205.4325 [hep-ph]]. 
[7] X. Liu and Z. -J. Xiao, Phys. Rev. D 86, 074016 (2012) arXiv: 1203.6135 [hep-ph]j. 
[8] G. L. Castro and J. H. Munoz, Phys. Rev. D 83, 094016 (2011) arXiv:1103.2993 [^hep-ph]]. 
[9] B. A. Li, Phys. Rev. D 55, 1436 (1997) hep-ph/9606402 . 
[10] N. Isgur, D. Scora, B. Grinstein and M. B. Wise, Phys. Rev. D 39, 799 (1989); D. Scora and N. Isgur, Phys. Rev. D 52, 
2783 (1995). 

[11] S. Weinberg, Phys. Rev. 112, 1375 (1958). 

[12] M. Fujikawa et al. [Belle CoUaboration], Phys. Rev. D 78, 072006 (2008) larXiv:0805."3773l fhep-ex]]. 
[13] S. Schael et al. [ALEPH Collaboration], Phys. Rept. 421, 191 (2005) hep-e x/0506"072] . 
[14] J. Beringer et al. (Particle Data Group), Phys. Rev. D 86, 010001 (2012). 
[15] G. Nardulh and T. N. Pham, Phys. Lett. B 623, 65 (2005). 

[16] A. Lusiani, PoS HQL 2010, 054 (2010) arXiv:1012.3733 [hep-ex]]. 

[17] M. Bona et al. [SuperB Collaboration], Pisa, Italy: INFN (2007) 453 p. www.pi.infn.it/SuperB/?q=CDR [arXiv:0709.045T] 
[hep-ex]]. 

[18] T. Aushev et al. [BELLE II Collaboration] , ,arXiv:1002. 5012 1 [hep-ex]. 
[19] C. Z. Yuan [BELLE II Collaboration], arXiv:1208.3813l [phvsics.ins-det].