arxiv: v2 [hep-ex] 16 Oct 2015

Transcription

1 Study of B Kππ Decays arxiv: v [hep-ex] 16 Oct 15 representing the collaboration Dept. of Physics, University of Bergen, Bergen, Norway gerald.eigen@ift.uib.no Using B B decays recorded with the detector at the PEP-II e e storage ring, we present the time-dependent CP asymmetry measurement in the radiative penguin decay mode B KS ρ(77, yielding S K S ρ =.17 ± The result is extracted from the time-dependent CP asymmetry parameters S K S π π =.14 ± and C KS π π =.39±±.±.5 measured in the neutral decay B KS π π after correcting for the dilution of K (89π in Kρ. The dilution factor D K S ρ = is determined from a study of the charged mode B K π π, which produces more signal events and is related to the neutral mode by isospin. We need a detailed knowledge of the resonance structure in the K π π mass spectrum and measure branching fractions of different resonances to Kπ and ππ final states. We also measure the branching fractions B(B K π π = (7. ± 1. ± and B(B K S π π = (4. ± The European Physical Society Conference on High Energy Physics 9 July 15 Vienna, Austria Speaker. This work is supported by the Norwegian Research Council. c Copyright owned by the author(s under the terms of the Creative Commons Attribution-NonCommercial-NoDerivatives 4. International License (CC BY-NC-ND 4..

2 1. Introduction The V-A structure of the Standard Model (SM weak interaction produces predominantly left-handed photons in b s decays. Thus apart from m s /m b effects [1], a B ( B meson decays predominantly to a right-handed (left-handed photon. In the SM, the mixing-induced CP asymmetry in B f CP decays is expected to be small in the SM where f CP is a CP eigenstate. However, new physics processes in which opposite-helicity photons are involved may alter the SM prediction [1,, 3, 4]. Inclusive and exclusive radiative decays have been studied by [5, 6] and Belle [7, 8] in several channels. For example, the inclusive branching fraction of B(B X s = (3.4 ± [9] is in good agreement with the SM prediction of B(B X s = (3.15 ± [1]. has studied the exclusive decay B Kππ exploring the resonance structure of the Kππ system. The data sample consists of B B events recorded with the detector at the PEP II e e storage ring at SLAC corresponding to an integrated luminosity of 46 fb 1. The goal consists of measuring the mixing-induced CP asymmetry parameter, S K ρ, in the B radiative decay to the CP eigenstate K S ρ, which is sensitive to right-handed photons. To accomplish this we measure the time-dependent CP asymmetry parameters S K ππ and C K ππ in the neutral decay B KS π π since S K S ρ = S K S π π D. The dilution factor D K K S ρ S ρ the amplitudes of the two-body decays ρ(77 KS, K (89 π and (Kπ π [11]: depends on [ AρK A D K S ρ = S K π A (Kπ π Re(A A ρks K π Re(A A ρks (Kπ π ] dm [ AρK A S K π A (Kπ π Re(A A ρks K π ] Re(A A ρks (Kπ π dm (1.1 Thus, we need to measure all two-body amplitudes to determine D K S ρ. For this task, we use the B K π π decay since due to higher signal yields we achieve more precise measurements and the charged mode is related to B KS π π by isospin. The Kππ final state is produced by several kaonic resonances. We extract these from a fit to the K π π mass spectrum.. Determination of the dilution factor D K S ρ using the decay B K π π We combine a high-energy photon (1.5 < E < 3.5 GeV with two charged pions and a charged kaon. The signal event selection is based on two kinematic observables, the beam-energyconstrained mass m ES = 1 4 E CM p B and the energy difference E = E B 1 E CM where E CM is the center-of-mass energy and p B,E B are momentum and energy of the B meson in the B rest frame. In addition, six event shape observables that are combined into a Fisher discriminant F [1] to separate signal from e e q q continuum background (q = u,d,s,c. We further use a likelihood ratio to discriminate against photons from π and η decays. Requiring L π (L η <.86 (.957 retains 93% (95% signal and removes 83% (87% q q and 63% (1% B B backgrounds. Using an extended unbinned maximum likelihood fit to m ES, E and F, we extract a B K π π signal yield of 441 ± events for m Kππ < 1.8 GeV/c. Using the s Plot technique [13], we extract the m Kππ, m Kπ and m ππ invariant-mass spectra. We model the m Kππ invariantmass spectrum with five resonances using coherent sums of same-spin terms ( K 1 (17K 1 (14,

3 Table 1: Fit results for amplitudes and phases for different K π π resonances obtained from the maximum likelihood fit to the m Kππ spectrum and computed branching fractions after efficiency correction and division by secondary branching fractions [9]. The first uncertainty is statistical, the second is systematic, and the third, when present, comes from the error of secondary branching fraction. Mode Amplitude [Phase (rad] B(B mode B(B mode PDG B(K res Kππ [1 6 ] [1 6 ] [9] [1 6 ] B K π π 7. ± 1. ± ±. B K 1 (17 1. (fixed [. (fixed] ± 13 B K 1 ( [ ] < CL B K ( ] B K ( [. (fixed] ± 4 B K ( [. (fixed] < CL K (141K (168, K (143, each parameterized by a relativistic Breit-Wigner line shape. The fit has eight free parameters, the magnitudes of the K 1 (14, K (141, K (168 and K (143, two relative phases and the widths of the K 1 (17 and K (168. We determine the resonance fit fractions and the interference fit fractions from an 8-bin maximum likelihood fit to the m Kππ spectrum. Table 1 lists the branching fractions of the individual kaonic resonances. Figure 1 (left shows the m Kππ s Plot spectrum with the fit overlaid. We measure the B K π π branching fraction to be B(B K π π = (7. ± 1. ± Events/(.13 GeV/c Norm. Residual preliminary (GeV/c m Kππ Events/(.19 GeV/c Norm. Residual (GeV/c m Kπ preliminary Figure 1: The m Kππ (left and m Kπ (right spectra for correctly-reconstructed B K π π signal FIG. 4. Distribution of mkππ for correctly-reconstructed B K π π signal events FIG. (splot, 6. Distribution extractedofrom mkπthe for maximum correctly-reconstructed B K π π signal events (splot, extracted from the maximum likelihood fit to mes, E, and the Fisher discriminant. Points with error bars give the sum likelihood of sweights. fit to mes, The E, blueand solid thecurve Fisher discriminant. Points with error bars give the sum of sweights. The blue solid the result of the fit performed directly to this mkππ distribution to extract the different curve contributions correspondsfrom to the kaonic total PDF resonances fit projection. The small-dashed red, medium-dashed green and dotted magenta curves events extracted from the maximum likelihood decaying to K π π correspond fittotothemk. Bins with negative entries is a feature of the splot technique used for background ES (89, E subtraction., ρ(77 and and F using (Kπ contributions, the Below s respectively. Plot technique. The dashed-dotted Points gray curve corresponds to the interference between the two P-wave components, i.e. the K each bin are shown the residuals, normalized in error units. The parallel dotted and full lines mark the one and two standard (89 and the ρ(77, and the dashed-triple-dotted light with error bars represent the sum of sweights. blue curve Blue corresponds solid to curves the interference showbetween the the fit (Kπ to the andmthe ρ(77. Bins with negative entries is a feature of the deviation levels, respectively. Kππ spectrum and the splot technique used for background subtraction. Below each bin are shown the residuals, normalized in error units, where the parallel dotted and full lines mark the one and two standard deviation levels, respectively. total PDF fit projection to m Kπ, respectively. The red dashed, green long-dashed and magenta dotted lines represent the K (89, ρ(77 and (Kπ contributions, respectively. The gray dash-dotted and the light blue dash-triple-dotted lines show the interference between K (89 and ρ(77 and that between (Kπ TABLE V. Results of the fit to the correctly-reconstructed signal splot of mkππ. The first uncertainty is statistical and second is systematic. The andlatter ρ(77 are discussed, respectively. in Sec. IV D 1. Interferences Note that fornegative both J P = TABLE 1 entries andviii. 1 resonances Branching are a feature fractions are destructive. of ofthethe different s Plot resonances technique decaying todue Kπ and toππ, background extracted from the fit to the mkπ spectrum. The quoted numbers are averaged over charge-conjugate states. They are obtained using the fit fraction of subtraction. J P Kres The Magnitude twoαlowerphase plots φ (rad. show residuals Fit each fraction component in units and the of corresponding standardefficiency. deviations. R denotes an The intermediate blackresonant dotted stateand hred stands for a final state hadron: a charged pion or kaon. To correct for the secondary branching fractions we used the values from Ref. [17] and solidk1(17 lines indicate 1. (fixed the one and. (fixed two standardb(k.61 deviation.8 (89.5 K levels, π =. The 3 respectively. first uncertainty is statistical, the second is systematic, and the third (when applicable is 1 due.5 to.5 the uncertainties on the secondary branching fractions. Systematics uncertainties are discussed in Sec. IVD. The last two rows of the table are obtained by separating the contributions from the resonant and the nonresonant part of the K1(14.7 ± ± LASS.5parametrization..3 Integrating separately the resonant part, the nonresonant part, and the coherent sum we find that the nonresonant part accounts for 95.6%, the resonant contribution 7.9%, and the destructive interference 3.5%. When the We further perform a 9-bin maximum K ( ± ± symbol.8.8 likelihood is quoted, it meansfit thatto the corresponding the m branching fraction was not previously reported..7.4 Kπ spectrum after correcting for the K (168 weighted.7 ±.8 efficiency in (fixed each bin..4 In.5.9 the fit B(B Mode.4.7 Mode we include the K (89 B(B 6 Previous world and Mode a 1(K π non- B(R hπ 1 6 average [17] ( 1 6 K (143.9 ± (fixed Sum of fit fractions interference J P = 1 : {K1(17 K1(14} J P = 1 : {K (141 K (168} Line-shape parameters Kres Mean ( GeV/c Width ( GeV/c K1( (fixed.98 ±.6 ±.5 3 K (89 π 17.3 ± K ρ(77 9. ± ±1.3 ±. < at 9% CL (Kπ π 11.4 ± (Kπ π (NR < 9. at 9% CL K (143 π.91 ± ± ±.14

4 Table : Fit results for amplitudes and phases and computed branching fractions for different resonances (R decaying to hπ (h = K,π obtained from the maximum likelihood fit to the m Kπ spectrum after efficiency correction and division by secondary branching fractions [9]. The first uncertainty is statistical, the second is systematic, and the third (when present comes from the error of the secondary branching fraction. The last two rows show resonant and non-resonant parts of the LASS parametrization. Mode Amplitude [Phase (rad] B(B Mode [1 6 ] PDG [1 6 ] [9] B K π π 7. ± 1. ± ±. K (89 π 1. (fixed [. (fixed] ± K ρ( [ ] ± 1.3 ±. CL (Kπ π [ ] 11.3 ± (Kπ π (NR < CL K (168 π 1.44 ± ±.14 resonant S-wave contribution. Furthermore, we include a contribution from the ρ that yields a broad structure in the K π mass spectrum. We model the K (89 with a relativistic Breit- Wigner line shape, the ρ(77 with a Gounaris-Sakurai line shape [14] and the (K π with the LASS parameterization [15]. For each resonance, we account for line shape distortions above the pole mass caused by low-mass Kππ resonances. We include interference between Kπ and ππ P-wave as well as interference between Kπ S-wave and ππ P-wave. Figure 1 (right shows the m Kπ s Plot spectrum with the fit components superimposed. Table lists the branching fraction of the different resonances decaying to K π and π π. This is the first observation of the decay B K ρ and the B (Kπ π S-wave contribution. From the (K π S-wave component, we separate the K (143 resonant contribution by performing separate integration of the resonant part, the non-resonant part and the coherent sum. We find 95.6% non-resonant, 7.9% resonant and -3.51% destructive interference contributions. For the resonant contribution, we extract a branching fraction of (1.44 ± ± where the first uncertainty is statistical, the second is systematic and the third results from secondary branching fractions. Since in the present analysis the K (143 contribution is modeled to come exclusively from B K 1 (17 K (143π, we determine the branching fraction to be B(K 1 (17 K (143π = %. This value is in good agreement with the Belle measurement [16]. Inserting the measured two-body amplitudes into Eqn (1, we obtain a dilution factor of D K S ρ = using the same mass constraints as in the B KS π π analysis, m Kππ < 1.8 GeV/c,.6 < m ππ <.9 GeV/c, m Kπ <.845 GeV/c and m Kπ >.945 GeV/c. 3. Time-dependent analysis of B K S π π The selection of the KS π π mode is the same as that for the K π π mode, except for the replacement of the K with a KS requiring m ππ m K < 11 MeV/c, a lifetime significance S of more than five standard deviations and p K l flight / p S KS K l flight >.995 where p S KS K is the K S S momentum and l flight is the vector connecting the B and K KS S decay vertices. We apply the same m Kππ, m Kπ and m ππ mass selections as for the charged mode. For most backgrounds, we use an 4

5 Argus PDF [17] for m ES, Chebychev polymomials for E and Gaussian or exponential functions for F. We optimize the selection on F to minimize statistical errors of CP parameters. We perform an extended unbinned maximum likelihood fit to extract the B KS π π signal yield along with time-dependent CP asymmetry parameters S K S ππ and C KS ππ. The likelihood function for event i is a sum over class j (signal and backgrounds in which each PDF depends on m ES, E, F and t. For most classes we can factorize the likelihood in the following way: P i j(m ES, E,F, t,σ t,q tag,c = P i j(m ES P i j( EP i j(f P i j( t,σ t ;q tag,c (3.1 where q tag = 1( 1 for B tag = B ( B and c represents the tagging category. We use six mutually exclusive tagging categories and collect all non-tagged events in a seventh category. We parametrize the proper time distribution for B K s ρ events for tagging category c by: exp[ t Psig( t,σ i t ;q tag,c = 4τ B τ B ] [ 1 1 q tag D c q tag < D > c ( Ssin( mb d t C cos( m B d t ] 8 Events/(.3 GeV/c 1 R c sig( t,σ t (3. where 8 < D > c is the average B B tagging dilution for category c and D c is the difference in D c between B and B 8 tags [18, 19]. We use the tagging algorithm and < D > c, D c values for the 6 six tagging categories from the B (c ck ( analysis 6 [19]. The seventh category of untagged events 4 is useful for the determination of the direct CP4asymmetry []. We add t background PDFs for charged B decays, B decays to flavor eigenstates, B decays to CP eigenstates as well as those for q q contributions. Figure (left shows the m ES projection of the fit. For t > ps and σ t <.5 ps, the maximum likelihood fit yields 46± signal events from which we compute a branching fraction of B(B 5.7 m (GeV/c KS π π. = (4. ± ES E (GeV/c Events/(.4 GeV/c 1 Events/(.1 1 Events/(.3 GeV/c Data Correctly-reconstructed signal Mis-reconstructed signal Continuum ± KS π ± KS π - K π Fisher Discriminant (GeV/c m ES Data Correctly-reconstructed signal Mis-reconstructed signal Continuum ± preliminary KS π ± KS π - K π Events/(.533 ps preliminary t (ps FIG. 8. Distributions of m ES (top left, E (top right, the Fisher-discriminant output (bottom left, and t (bottom right, showingfigure the fit results : Theonm the B KSπ π ES (left and t (right dataprojections sample. The ofdistributions the fit to thehave B their K signal/background ratio enhanced by S π π sample after requiring means of the following requirements:.15 E.1 GeV (m ES; ES GeV/c ( E; m ES > 5.7 GeV/c, < E <.1 GeV for the m E.1 GeV (Fisher and t. Points with ES and t spectra plus m error bars give the data. ES > 5.7 GeV/c for the t spectrum. Points The projection of the fit result is represented by the stackedwith histogram, error bars where show the shaded the data areas andrepresent stacked histograms the contributions show from the fit different projections. backgrounds, Different as shadings described represent in the legend. Some ofdifferent the contributions backgrounds are hardly as described visible due in the to their legend. small fractions. Note that the same order for the various contributions is used in both the stacked histogram and the corresponding legend. Figure (right shows the t distribution, which looks rather symmetric indicating small CP with m Kππ < 1.8 GeV/c,.6 < m ππ <.9 GeV/c 1377 violation as expected in the SM. The time-dependent and branching fractions of the ρ K, K π and (Kπ π 1399 CP parameters are measured to be S K with m Kπ <.845 GeV/c or m Kπ >.945 GeV/c S π π = 1378, we 14 components are also measured obtain the CP -violating parameters S K S π π =.14 ± and C 138 KS π π =.39 ±. ±.5, where the first quoted uncertainties are statistical and the sec- 141 ACKNOWLEDGMENTS 138 ond are systematic. We extract from this measurement the time-dependent CP asymmetry related to the B 1383 KSρ decay and obtain S K S ρ =.17 ± A 14 We thank E. Kou for help throughout the analysis, especially with the computation of the dilution factor, 1385 larger data sample is needed in order to draw conclusions 144 and A. Le Yaouanc for valuable discussions. We are 1386 regarding the contribution from right-handed photons in 145 grateful for the extraordinary contributions of our PEP b s transitions. This measurement of time dependent 146 II colleagues in achieving the excellent luminosity and

6 .14 ± and C KS π π =.39 ±. ±.5. After correcting for D KS ρ, we measure a mixing-induced CP parameter of S K S ρ =.17 ± Our result is consistent with the Belle result [1]. Both results are consistent with zero and thus agree with the SM prediction. 4. Conclusion We observed the decays B K π π and B KS π π and measured their branching fractions. In the charged mode, we observed five kaonic resonances decaying to Kππ for which we measured fit fractions and branching fractions. We found first evidence for B K 1 (14, B K (141 and B K (14. We further determined the dilution factor D Kρ after measuring amplitudes and phases of B K (89 π, B K ρ(77 and B (K π π modes. The latter two results are first measurements. In the neutral mode, we measured the timedependent CP asymmetry parameters S K S π π and C KS π π and in turn the CP asymmetry parameter S K S ρ for the K S ρ(77 CP eigenstate. Presently, experimental uncertainties are too large to set meaningful limits on right-handed photons. This remains a task for Belle II. References [1] W. Atwood, M. Gronau and A. Soni, Phys. Rev. Lett. 79, 185 (1997. [] K. Fujikawa and A. Yamada, Phys. Rev. D49, 589 (1994. [3] K.S. Babu, K. Fujikawa and A. Yamada, Phys. Lett. B333, 196 (1994. [4] P. L. Cho and M. Misiak, Phys. Rev. D49, 5894 (1994. [5] J. P. Lees et al. (, Phys.Rev.D86, 118 (1; ibid, Phys.Rev.Lett. 19, (1. [6] B. Aubert (, Phys. Rev. Lett. 88, [7] A. Limosani et al. (Belle, Phys. Rev. Lett. 13, 4181 (9; T. Saito et al. (Belle, Phys.Rev. D91, 5, 54 (14. [8] N. Nakao et al. (Belle, Phys.Rev. D69, 111 (4. [9] K. Olive et al. (Particle Data Group, Chin. Phys. C38, 91 (14. [1] M. Misiak et al., Phys. Rev. Lett. 98, (7. [11] J. Hebinger et al., LAL-15-75; note-v3.pdf. [1] R. A. Fisher, Annals Eugen. 7, 179 (1936. [13] M. Pivk and F.R. Le Diberder, Nucl. Instrum. Meth. A555, 356 (5. [14] G.J. Gounaris and J.J. Sakurai, Phys. Rev. Lett. 1, 44 (1968. [15] D. Aston et al., Nucl. Phys. B96, 493 (1988. [16] H. Guler et al. (Belle, Phys. Rev. D83, 35 (11. [17] H. Albrecht (ARGUS, Z. Phys. C48, 543 (199. [18] B. Aubert (, Phys. Rev. Lett. 94, (5. [19] B. Aubert (, Phys. Rev. Lett. 99, (7. [] S. Gardner and J. Tandean, Phys. Rev. D69, 3411 (4. [1] J. Li et al. (Belle, Phys. Rev. Lett. 11, 5161 (8. 6