arxiv: v3 [hep-ph] 15 Sep 2009

Transcription

1 Determination of β in B J/ψK+ K Decay in the Preence of a K + K S-Wave Contribution Yuehong Xie, a Peter Clarke, b Greig Cowan c and Franz Muheim d arxiv:98.367v3 [hep-ph 15 Sep 9 School of Phyic and Atronomy, Univerity of Edinburgh, Mayfield Road, Edinburgh, EH9 3JZ, U.K. Abtract We preent the complete differential decay rate for the proce B J/ψK + K including S-wave and P-wave angular momentum tate for the K + K meon pair. We examine the effect of an S-wave component on the determination of the CP violating phae β. Data from the B-factorie indicate that an S-wave component of about 1% may be expected in the φ(1) reonance region. We find that if thi contribution i ignored in the analyi it could caue a bia in the meaured value of β toward zero of the order of 1%. When including the K + K S-wave component we oberve an increae in the tatitical error on β by le than 15%. We alo point out the poibility of meauring the ign of co β by uing the interference between the K + K S-wave and P-wave amplitude to reolve the trong phae ambiguity. We conclude that the S-wave component can be properly taken into account in the analyi. a Yuehong.Xie@cern.ch b peter.clarke@ed.ac.uk c g.cowan@ed.ac.uk d f.muheim@ed.ac.uk

2 Content 1 Introduction 1 Time-dependent angular ditribution in the decay B J/ψK + K including S-wave contribution 1 3 Meauring β in the preence of a K + K S-wave 4 4 Meauring co β 7 5 Concluion 8 1 Introduction The decay B J/ψφ i a golden channel for the meaurement of the B mixing phae β which i a very enitive probe of new phyic. It ha been extenively tudied [1,,3,4,5,6,7,8. In the decay B J/ψφ, followed by a two-body decay φ(1) K + K, the K + K meon pair i in an orbital P-wave amplitude. However, in the vicinity of the φ(1) ma, the K + K ytem can have contribution from other partial wave. The ame comment hold for the K + K ytem in the decay channel B K + K K S and D K + K π. The BaBar experiment howed that in thee decay the S-wave and P- wave contribution dominate in the ma range above threhold up to 1.1 GeV/c [9, 1. In both cae there i a dominant reonant φ(1) contribution. In addition an S-wave f (98) and a non-reonant contribution are found to be neceary to decribe the data. Thee reult motivated u to invetigate the effect of a poible S-wave contribution to B J/ψK + K in the φ(1) ma region. In the decay B J/ψK + K the K + K ytem can only arie from a quark pair while in B K + K K S and D K + K π it can have contribution from both and dd. Thi make it difficult to give a quantitative etimate for the S-wave component. In reference [11 the S-wave K + K contribution under the φ(1) peak i etimated to be 5 1% for decay mode in which the K + K arie from an quark pair. In thi tudy we conider an S-wave of imilar magnitude and ae it impact on the determination of the weak mixing phae β. Time-dependent angular ditribution in the decay B J/ψK + K including S-wave contribution We conider P- and S-wave amplitude in the decay B J/ψK + K where the invariant ma of the K + K meon pair i in the φ(1) ma region and the J/ψ meon decay into a 1

3 µ + µ pair. The S-wave contribution can be non-reonant or due to the f (98) reonance 1. We denote decay amplitude for the B J/ψK + K by A = (A, A, A, A S ). Here A, A and A are the three P-wave amplitude conitent with the K + K ytem decaying via the φ(1) reonance. A S i the amplitude for a poible S-wave contribution in the K + K ytem. The amplitude for the conjugate decay B J/ψK + K are denoted by Ā = (Ā, Ā, Ā, ĀS), which, in the abence of direct CP violation, are related to A by A = Ā, A = Ā, A = Ā and A S = ĀS. Note that A and A are CP-even wherea A and A S are CP-odd. The amplitude (A, A, A ) and the amplitude A S may have different dependence on the ma m K + K of the K+ K ytem. However, in ufficiently mall bin of m K + K, uch a the narrow ma region around the φ(1) reonance, the dependence of the amplitude on m K + K can be neglected. We define the total P-wave trength, A P A + A + A, the longitudinal and perpendicular polariation fraction relative to the P-wave trength R A /A P, and R A /A P, and the S-wave fraction, R S A S /(A P + A S ). The phae of thee decay amplitude are defined by A j = A j e iδ j, where j =,,, S. A only the relative trong phae difference can be meaured we adopt the convention δ =. An angular analyi i required to dientangle the different CP eigentate on a tatitical bai. The angular obervable are denoted a the helicity angle Ω = (θ l, θ K, ϕ). Here θ l i the angle between the µ + momentum and the direction oppoite to the B momentum in the J/ψ ret frame; θ K i the angle between the K + momentum and the direction oppoite to the B momentum in the ret frame of the K + K ytem; ϕ i the angle between the decay plane of the J/ψ µ + µ and the K + K pair, when going from the poitive kaon to the poitive lepton with a rotation around the oppoite direction of the B momentum in the J/ψ ret frame. The differential decay rate for a B meon produced at time t = decaying a B J/ψK + K at proper time t i given by d 4 Γ(B J/ψK + K ) dt d co θ d co ψ dϕ 1 k=1 wherea the differential decay rate for an initial B meon i given by d 4 Γ(B J/ψK + K ) dt d co θ d co ψ dϕ 1 k=1 h k (t)f k (Ω), (.1) h k (t)f k (Ω). (.) Each of the h k (t), hk (t) and f k (Ω) for k = 1 1 are defined in Table 1. In total there are four amplitude-quared term for the three polariation of the P-wave and the S-wave component plu ix interference term. The time-dependence of the ten function h k (t) for an initial B meon tate can be written a: A (t) = A e Γt[ coh co Φ inh A (t) = A e Γt[ coh co Φ inh + in Φ in( m t), (.3) + in Φ in( m t), (.4) 1 The ma dependence of the f (98) i ditorted a the central value of the reonance i below threhold.

4 k h k (t) hk (t) f k (θ l, θ K, ϕ) 1 A (t) Ā(t) 4 in θ l co θ K A (t) Ā (t) (1 + co θ l ) in θ K in θ l in θ K co ϕ 3 A (t) Ā (t) (1 + co θ l ) in θ K + in θ l in θ K co ϕ 4 I{A (t)a (t)} I{Ā (t)ā (t)} in θ l in θ K in ϕ 5 R{A (t)a (t)} R{Ā (t)ā (t)} in θ l in θ K co ϕ 6 I{A (t)a (t)} I{Ā (t)ā (t)} in θl in θ K in ϕ 7 A S (t) ĀS(t) 4 3 in θ l 8 R{A S (t)a (t)} R{Ā S (t)ā (t)} 3 6 in θl in θ K co ϕ 9 I{A S (t)a (t)} I{Ā S (t)ā (t)} 3 6 in θl in θ K in ϕ 1 R{A S (t)a (t)} R{Ā S (t)ā(t)} in θ l co θ K Table 1: Definition of the function h k (t), hk (t) and f k (θ l, θ K, ϕ) of Eq..1 and.. A (t) = A e Γt[ coh + co Φ inh I{A (t)a (t)} = A A e Γt[ ( co(δ δ ) in Φ inh ) in Φ in( m t), (.5) + in(δ δ ) co( m t) co(δ δ ) co Φ in( m t), (.6) R{A (t)a (t)} = [ A A e Γt co(δ δ ) coh co Φ inh + in Φ in( m t), (.7) I{A (t)a (t)} = A A e Γt[ co(δ δ ) in Φ inh + in(δ δ ) co( m t) co(δ δ ) co Φ in( m t), (.8) A S (t) = A S e Γt[ coh + co Φ inh in Φ in( m t), (.9) R{A S(t)A (t)} = A S A e Γt[ in(δ δ S ) in Φ inh + co(δ δ S ) co( m t) in(δ δ S ) co Φ in( m t), (.1) I{A S(t)A (t)} = [ A S A e Γt in(δ δ S ) coh + co Φ inh in Φ in( m t), (.11) R{A S(t)A (t)} = A S A e Γt[ in(δ δ S ) in Φ inh + co(δ δ S ) co( m t) in(δ δ S ) co Φ in( m t), (.1) where Φ = β, m, Γ and Γ denote the weak mixing phae, ma difference, decay width difference and average decay width of the B -B ytem, repectively. Here we have aumed that each of the decay amplitude in A i dominated by a ingle weak phae, therefore a common effective β can be ued for all CP eigentate. The time evolution function h k (t) for an initial B meon can be obtained by revering the ign of each term proportional to in( m t) or co( m t) in h k (t). 3

5 3 Meauring β in the preence of a K + K S-wave In thi ection we invetigate how the meaurement of β i affected by the preence of a poible K + K S-wave contribution. We ue Monte Carlo imulated toy data baed on the differential decay rate expreion of Section. We generate ignal decay only and ignore background underneath the B ma peak a well a all detector effect. The incluion of thee effect doe not alter the qualitative reult of thi tudy. We aume a tagging efficiency ɛ tag = 56% and a wrong tag probability ω tag = 33%, which correpond approximately to the expected flavour tagging performance for thi channel at the LHCb experiment [1. In Table we ummarize the value of the phyical parameter ued to generate the toy data et. We generate 5 data et for different cenario where we vary the value of the S- wave fraction R S and it phae δ S and the weak phae β. Each data et contain 3 ignal event correponding to approximately one quarter of a nominal LHCb year of fb 1. m Γ Γ δ δ δ R R R S δ S β Input 17.8 p 1.68 p 1.5 p vary vary vary Fit fix fix fix fix float float float float float float float R S i fixed to when the S-wave component i neglected. Table : Value of the phyical parameter ued in the generation of ignal decay and how thee parameter are treated in the fit. We perform fit to each data et where β, R S, δ S, R, δ, R, δ are free parameter and all other parameter are kept fixed. We alo perform fit where the S-wave component i preent in the generated toy data, but ignored in the fit (R S i et to ) in order to invetigate the bia in the determination of β. The reult of thee fit for the tatitical error and mean value of the weak phae β are ummarized in Table 3, 4 and 5 for everal different cenario with β =.368, β =. and β =.5, repectively. A an example, in Figure 1 we how the ditribution of the fitted value of β for R S =.1, δ S = π/ and β =.5 for both the S-wave fraction R S fixed to zero and R S left free in the fit. In Figure we how the ditribution of the fitted value of R S and the trong phae of the S-wave component δ S for the ame cae with R S left free in the fit. It can be een that when all parameter are fitted the reult are unbiaed, but when it i wrongly aumed that R S =, the reult for β acquire a bia with regard to the true input value. Figure 3 how the bia in β from neglecting an S-wave component with R S =.1 and δ S = π/ veru the value of β ued to generate the data et. A linear dependence i oberved, which demontrate that the bia in β i proportional to the true value of β. From Table 3, 4 and 5 we oberve biae for thee cenario which range from 7 17% in the meaurement of β if an S-wave component i preent, but left unaccounted for in the fit. The bia move the meaured value of β toward zero. Thi implie that the neglected CP-odd S-wave contribution ha a bigger probability to be mi-identified a the CP-even longitudinal or parallel component than a the CPodd perpendicular component. Therefore, although the bia from neglecting an S-wave contribution i unlikely to lead to fale ignal of new phyic, it will caue a lo of 4

6 enitivity to new phyic. On the other hand, including the S-wave in the fit remove the bia in the central value of β at a cot of an increae of le than 15% in the tatitical error. Float R S in fit Fix R S to in fit R S = σ(β ) =.45, Mean(β ) =.38 R S =.1, δ S = π/ σ(β ) =.48, Mean(β ) =.35 σ(β ) =.45, Mean(β ) =.3 R S =.1, δ S = σ(β ) =.54, Mean(β ) =.4 σ(β ) =.48, Mean(β ) =.36 R S =.5, δ S = π/ σ(β ) =.48, Mean(β ) =.4 σ(β ) =.45, Mean(β ) =.36 R S =.5, δ S = σ(β ) =.55, Mean(β ) =.38 σ(β ) =.47, Mean(β ) =.3 Table 3: Statitical error and mean value of β from 5 fit for different cenario with β =.368. The error on σ(β ) and mean(β ) are approximately.3 and., repectively. The ame data et are ued to obtain the reult in the econd and third column. Float R S in fit Fix R S to in fit R S = σ(β ) =.44, Mean(β ) =.198 R S =.1, δ S = π/ σ(β ) =.5, Mean(β ) =.199 σ(β ) =.47, Mean(β ) =.166 R S =.1, δ S = σ(β ) =.56, Mean(β ) =. σ(β ) =.49, Mean(β ) =.17 R S =.5, δ S = π/ σ(β ) =.49, Mean(β ) =.197 σ(β ) =.48, Mean(β ) =.18 R S =.5, δ S = σ(β ) =.53, Mean(β ) =.198 σ(β ) =.48, Mean(β ) =.18 Table 4: Statitical error and mean value of β from 5 fit for different cenario with β =.. The error on σ(β ) and mean(β ) are approximately.3 and., repectively. The ame data et are ued to obtain the reult in the econd and third column. Float R S in fit Fix R S to in fit R S = σ(β ) =.51, Mean(β ) =.51 R S =.1, δ S = π/ σ(β ) =.59, Mean(β ) =.51 σ(β ) =.53, Mean(β ) =.415 R S =.1, δ S = σ(β ) =.61, Mean(β ) =.51 σ(β ) =.5, Mean(β ) =.417 R S =.5, δ S = π/ σ(β ) =.51, Mean(β ) =.56 σ(β ) =.48, Mean(β ) =.463 R S =.5, δ S = σ(β ) =.53, Mean(β ) =.51 σ(β ) =.49, Mean(β ) =.461 Table 5: Statitical error and mean value of β from 5 fit for different cenario with β =.5. The error on σ(β ) and mean(β ) are approximately.3 and., repectively. The ame data et are ued to obtain the reult in the econd and third column. 5

7 a plot of -β Event / (.8 ) Prob 1 Contant 33.6 ±.9 Mean ±.78 Sigma.488 ±.146 a plot of -β Event / (.8 ) Prob 1 Contant 7.17 ± Mean -.56 ±.316 Sigma.5866 ± β β Figure 1: Ditribution of the fitted value of β for the cenario R S =.1, δ S = π/, β =.5. The left and right plot are obtained with or without fixing R S to in fitting the data, repectively. a plot of R S Event / (.5 ) Prob 1 Contant 13 ± Mean.111 ±.8973 Sigma.776 ±.5574 a plot of δ S Event / (. ) Prob 1 Contant ± Mean 1.56 ±.4613 Sigma.7435 ± R S δ Figure : Ditribution of the fitted value of R S and δ S for the cenario R S =.1, δ S = π/, β =.5 without fixing R S to in fitting the data. ) (-β.1.8 Prob.3174 p ±.17 p ± input -β Figure 3: The bia in β from neglecting an S-wave component with R S =.1 and δ S = π/ veru the value of β ued to generate the data et. The bia i the difference of the mean of the fitted to the generated β value. A linear fit i uperimpoed on the graph. 6

8 4 Meauring co β In Eq..1 and. one oberve that the differential decay rate are invariant under the tranformation ( δ δ, δ δ, δ S δ, β, Γ ) ( δ δ, π + δ δ, δ δ S, π ( β ), Γ ). (4.1) A a conequence the meaurement of β i ubject to a two-fold ambiguity, which i equivalent to co β tranforming into co β. A meaurement of co β including it ign would allow u to reolve thi ambiguity. If the interference between the P-wave and S-wave amplitude were to be ignificant in the φ(1) ma region, we could ue thi effect to meaure co β, in the ame way a BaBar meaured co β in B J/ψK Sπ [13. Thi require meauring δ S δ, the trong phae difference between the S-wave and the longitudinal P-wave, a a function of the K + K ma in the φ(1) ma region. When plotting thi function, two branche are expected with each correponding to a different olution for the weak phae (ee Figure 4 left). It i traightforward to chooe the phyical olution ince the phae of the P-wave Breit-Wigner amplitude i expected to rie rapidly through the φ(1) ma region (dahed red curve in Figure 4 right), while the phae of the S-wave amplitude, which can be decribed either by a coupled channel Breit-Wigner function in cae of an f contribution or by a contant term in cae of a non-reonant contribution, i expected to vary relatively lowly (dotted green curve in Figure 4 right), reulting in δ S δ rapidly falling with increaing K + K ma (olid blue curve in Figure 4). δ -δ (Rad) S trong phae (Rad) K K invariant ma (MeV) K K invariant ma (MeV) Figure 4: An example to illutrate the dependence of the trong phae of the S-wave δ S, of the trong phae of the longitudinal P-wave δ, and of their difference δ S δ, on the K + K ma. Left: the olid blue curve i the phyical olution for δ S δ and the dahed black curve how the mirror olution. Right: the dahed red, dotted green and olid blue curve are for δ, δ S, and δ S δ, repectively. Below we ue a Monte Carlo imulated toy data et to demontrate the feaibility of thi method in meauring the ign of co β. We generate 3 B J/ψK + K event in the K + K ma region between 1 and 1.5 GeV/c, roughly correponding to.5 fb 1 of integrated luminoity. The P-wave and f contribution are included coherently. The value of the parameter ued to generate the toy data et are the ame a in Table except that we et β =.368, and that the value of both R S and δ S depend on 7

9 the K + K ma. The f contribution account for about 1% of the total decay rate in the given ma region, a i hown in Figure 5. Event / (.5 ) K K invariant ma (MeV) Figure 5: The data point correpond to the K + K ma ditribution of a generated ample of B J/ψK + K event including 1% f contribution in the ma region. The dotted red curve indicate the f contribution. The data ample i divided into bin in the K + K ma. For each bin i, two parameter δ S,i and R S,i are ued to repreent the average trong phae and the fraction of the f contribution. Both in β and co β are treated a independent free parameter. Common free parameter in β, co β, R, R, δ, δ, Γ and Γ are ued for all bin. Note that we till adopt the convention δ = a only the relative phae difference in each bin can be meaured. A combined fit to the time-dependent angular ditribution of all the bin i performed to extract thee free parameter. The fitted value of the trong phae difference δ S δ veru the K + K ma are plotted in Figure 6. The two branche correpond to oppoite value of co β. Jut a expected, the branch correponding to the true olution decreae rapidly around the nominal φ(1) ma. Chooing thi branch lead to the unique olution which give the ambiguity-free reult in β =.43 ±.5, co β = 1.5 ±.8, (4.) β =.43 ±.5. (4.3) In thi example, the meaured β i eparated from π ( β ) by 13σ, therefore the dicrete ambiguity in β i completely reolved. Although the actual meaurement preciion in co β will depend on the ize of the f contribution a well a background, the poibility to reolve the ambiguity in β uing thi method i very promiing. 5 Concluion In the decay B J/ψK + K we expect that a K + K S-wave contribution in the narrow φ(1) ma region could be a large a 1%. The full differential decay rate for thi decay including the S-wave contribution have been preented. We have conidered a range 8

10 (Rad) 3 1 -δ S δ K K invariant ma (MeV) Figure 6: The fitted value of δ S δ veru K + K ma are hown in red and black data point, correponding to oppoite value of co β. The blue curve how the dependence of δ S δ on K + K ma implemented in imulation. of cenario which include S-wave component of 5% and 1%. We have hown that within thee cenario, if an S-wave component i ignored in the analyi, the meaurement of the weak phae β would be biaed by between 7% and 17% toward zero. We have demontrated that by properly allowing for thi S-wave component in the fit, an unbiaed meaurement of β may be obtained with a lightly increaed tatitical error. Finally, we have hown that the interference between the K + K S-wave and P-wave amplitude can be ued to reolve the two-fold ambiguity in the meaurement of the weak phae β. Acknowledgment The author would like to acknowledge the LHCb colleague for ueful and timulating dicuion. We particularly thank Tim Gerhon, Olivier Leroy and Guy Wilkinon for valuable uggetion. 9

11 Reference [1 M. Artuo et al., B, D and K decay, Eur. Phy. J. C 57 (8) 39 [arxiv: [ P. Ball and R. Fleicher, Probing new phyic through B mixing: tatu, benchmark and propect, Eur. Phy. J. C 48 (6) 413 [hep-ph/6449. [3 P. Clarke, C. Mclean and A. Oorio-Olivero, Senitivity tudie to β and Γ uing the full B J/ψφ angular analyi at the LHCb, CERN-LHCB [4 M. Ciuchini et al., Next-to-leading order trong interaction correction to the (F ) = effective hamiltonian in the MSSM, JHEP 9 (6) 13 [hep-ph/ [5 A. Lenz and U. Nierte, Theoretical update of B B mixing, JHEP 6 (7) 7 [hep-ph/ [6 Z. Ligeti, M. Papucci and G. Perez, Implication of the meaurement of the B B ma difference, Phy. Rev. Lett. 97 (6) 1181 [hep-ph/6411. [7 CDF collaboration, T. Aaltonen et al., Firt flavor-tagged determination of bound on mixing-induced CP violation in B J/ψφ decay, Phy. Rev. Lett. 1 (8) [8 D collaboration, V. M. Abazov et al., Meaurement of B mixing parameter from the flavor-tagged decay B J/ψφ, Phy. Rev. Lett 11 (8) 4181 [arxiv:8.55. [9 BABAR collaboration, B. Aubert et al., Meaurement of CP-violating aymmetrie in the B K + K K S dalitz plot [arxiv:88.7. [1 BABAR collaboration, B. Aubert et al., Amplitude analyi of the decay D K K + π, Phy. Rev. D 76 (7) 111 [arxiv: [11 S. Stone and L. Zhang, S-wave and the Meaurement of CP Violating Phae in B Decay, Phy. Rev. D 79 (9) 744 [arxiv: [1 M. Calvi, O. Leroy and M. Muy, Flavour tagging algorithm and performance in LHCb, CERN-LHCB [13 BABAR collaboration, B. Aubert et al., Ambiguity-free meaurement of co β: Time-integrated and time-dependent angular analye of B J/ψKπ, Phy. Rev. D 71 (5) 35 [hep-ex/