The b quark fragmentation fractions at LHCb and meson decays with heavy quark spectators

We study the current estimates of Bc→Bsπ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_c\rightarrow B_s \pi $$\end{document} to extract the fragmentation fraction fc/fs\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_c/f_s$$\end{document} at the LHCb. A rather robust estimate of Br(Bc→Bsπ)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Br(B_c\rightarrow B_s \pi )$$\end{document} based on factorization and lattice results for the form factor gives fc/fs∼0.056\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_c/f_s \sim 0.056$$\end{document} with a 16% error. We also revisit the extraction of fs/fd\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_s/f_d$$\end{document} using B→Dπ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B\rightarrow D\pi $$\end{document} instead of the theoretical cleaner but more suppressed channel B→DK\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B\rightarrow DK$$\end{document}. We also find a tension on the predictions of Br(Bc→J/ψπ)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Br(B_c\rightarrow J/\psi \pi )$$\end{document} and Br(Bc→Bsπ)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Br(B_c\rightarrow B_s\pi )$$\end{document} considering the measurements of these modes at LHCb, and find that, within a 23% uncertainty, only the lower end of the current prediction range Br(Bc→J/ψ)∼0.4%-1.7%\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Br(B_c\rightarrow J/\psi )\sim 0.4\% - 1.7\%$$\end{document} would be consistent with the LHCb measurements.


Introduction
The heavy meson decay B c → B s π is a remarkable case of a decay of a heavy meson where the heavy quark is a spectator of the weak interaction process, in contrast to the much studied transitions such as B → D, where it is the heavy quark that decays and the light degrees of freedom act as spectators.The latter case has been extensively studied, giving birth to the systematic treatment called heavy quark effective theory (HQET), which is a expansion in operators with coefficients in powers of Λ QCD /m b and Λ QCD /m c [1,2,3].As such, in the infinite mass limit, the results are relatively simple [4,5,6].For example, for a transition B → Dℓν at the kinematic point where B and D are at relative rest (the so-called zero recoil point), the wave functions of the two heavy mesons should be the same and the overlap of the form factors should be equivalent to their normalization (in the appropriate normalization, this means the form factor should be unity).In contrast, here we are interested in transitions such as B c → B s , where the heavy (or better said the heaviest ) quark is the spectator and the weak transition occurs in the lighter degrees of freedom.Now the wavefunctions will not coincide, not even at zero recoil, so there is no a priori known value for the form factor at any kinematic point.These decays, where the heavy quark is a spectator, therefore offer new information for hadron structure and bound states under the strong interactions.To have an intuitive picture let us consider the heavy meson as a bound state of two constituent quarks.The Bohr radius of the wavefunction goes as the inverse of the reduced mass.In decays such as B → Dℓν, in the limit m b , m c → ∞ the reduced mass is the same in both mesons, so one expects that the form factor at zero recoil, which is the overlap of the B and D wavefunctions, should be the same as their normalization (i.e. the form factor goes to a known value, usually taken as unity).However in B c → B s transitions the B c and B s reduced masses and Bohr radii are different, so the form factor at zero recoil is not a fixed norm but depends on the bound state dynamics.The first measurement (and so far the only one) of a B c → B s transition was done by LHCb [7] in the two-body decay B + c → B s π + , however the result on this branching ratio comes combined with the ratio of the fragmentation fractions of the b quarks produced in the pp collisions into the different B mesons, in this case σ(B c )/σ(B s ) ≡ f c /f s .It is therefore necessary to have an independent estimate of either the branching ratio or of f c /f s , in order to extract a value for the other.In Ref. [7] the authors exemplify the use of independent f c /f s estimates to determine the B c → B s π branching ratio.Here we will take the opposite avenue, namely to use the best theoretical estimate we find for the branching ratio, which is based on current lattice results for the form factors [8], in order to extract the value of f c /f s from the LHCb measurement.Using this estimate for f c /f s we then obtain estimates for other B c branching ratios such as B(B + c → J/ψπ + ), which in Ref. [7] was used as an input.Concerning the b fragmentation fractions at LHCb, another important ratio is f s /f d .This ratio was estimated earlier at LHCb, using the decays B s → D − s π + and B 0 → D − K + [9], and more recently doing a global fit using several other decays [10], thus reducing the uncertainty.In the former study, the reason for using B → DK instead of the more numerous B → Dπ was that its theoretical calculation contains an additional amplitude that is more difficult to estimate.Here we revisit these estimates using updated B meson branching ratios as well as more recent lattice results for the relevant form factors.In Section II we address the analysis of f c /f s and B c → B s transitions.In Section III we revisit the estimates and uncertainties of f s /f d and the non-factorizable amplitude of B → Dπ.In Section IV we revisit Br(B c → J/ψ π + ) and the extraction of f c /f d .Our conclusions are in Section V.
2 B c → B s transitions and f c /f s The decay B c → B s π was observed at LHCb [7] and its branching ratio was determined up to the factor Here f c and f s are the fragmentation fractions, i.e. the probabilities that a b quark produced in the pp collisions hadronizes into a B c or a B s meson, respectively, and Br stands for branching ratio.The ratio f c /f s appears because one does not know a priori the number of B c produced in the collisions, and so the count of B c → B s π events is given relative to the number of B s produced.Since in that analysis the B s were identified only by their decays B s → D s π and B s → J/ψ φ, then B c → B s π are identified by the events (D s π)π and (J/ψ φ)π, where the pair reconstructs the B s mass and the pair plus the bachelor pion reconstructs the B c mass.One can understand Eq.(1) by first denoting B s → D s π and B s → J/ψφ generically as B s → X; then the number of observed events B c → Xπ, where X reconstructs the B s mass, can be expressed as: where N b is the number of b quarks produced in the pp collisions and N b • f c -an unknown-is the total number of B c produced.In the same way, the total number of observed B s → X events (regardless whether B s comes from a B c decay or not) is: Dividing these two expressions one obtains Eq.( 1), where the ratio of counts n(B c → Xπ)/n(B s → X) is the result on the right hand side, measured by LHCb [7].
In order to separate Br(B c → B s π) from f c /f s in Eq. ( 1), one needs an independent input.In Ref [7] this was done by separately estimating f c /f s , expressing it as: where the denominator was previously given in Ref. [9] as f s /f d = 0.256 ± 0.020, and the numerator f c /f d was extracted from a previous LHCb result [11]: While Br(B → J/ψ K) = (1.016± 0.033) × 10 −3 has been measured [12], for Br(B c → J/ψ π) ∼ 0.06% → 0.18% these are only theoretical estimates [13].The latter clearly brought the largest uncertainty, while the next largest uncertainties are ±17% in Eq.( 2), ±15% in Eq.( 1) and ±7.8% in f s /f d according to Ref. [9].
Consequently the extracted Br(B c → B s π) from this method reproduces the large uncertainty present in Br(B c → J/ψ π): and similar uncertainties were then obtained for f c /f s ∼ 1.6% → 4.6% and f c /f d ∼ 0.4% → 1.2%.
Here we take the opposite avenue, namely to start from current lattice estimates for the relevant form factor of Br(B c → B s π) [8] to derive the value of f c /f s from Eq. (1).
The expression for the amplitude of B + c → B 0 s π + , assuming factorization, is [14], where (ūd) and (sc) denote V − A color singlet currents.Within this assumption the matrix element of the 4-quark operator is separated into a product of current matrix elements.Non-factorizable contributions, where gluons are exchanged between these currents, are neglected (for more details see Section III).Moreover, a 1 is the corresponding combination of Wilson coefficients at the renormalization scale of B mesons, µ ∼ m b [15].In Ref. [15], a value a 1 ≃ 1.02 ± 0.01 is estimated, but one can find other estimates as well, which depend also on the flavor transition.Larger values at the B scale have also been used, e.g. a 1 = 1.14 [13,16,17], or even a 1 = 1.26 [18].Several authors have used a value a 1 ≃ 1.2 for transitions c → s in beauty meson decays such as B c → B s π, and a value a 1 ≃ 1.14 for transitions b → c such as B c → J/ψπ [17,18,16,19,14,20,21,22,13].Since a 1 is a combination of Wilson coefficients that integrate QCD corrections from the electroweak scale M W down to µ ∼ m b , its estimates can be done systematically.In a detailed work that includes next-to-leading-order (NLO) calculations the authors of Ref. [23] find a 1 ≃ 1.05 with a 2% uncertainty by varying the scale µ ∼ m b /2 → 2m b .A more recent result at NNLO can be found in Ref. [24] where a 1 = 1.070 +0.010 −0.013 is obtained for B d → D + K − and very similar results for other final state mesons; here the quoted errors include the uncertainties from the variation of the scale µ of the Wilson coefficients, quark masses, Gegenbauer moments of the light meson wave functions, and the strong coupling α s (m Z ).[24] In order to separate the uncertainties in B meson decays coming from a 1 , wherever possible we will cite the results on theoretical branching ratios with a 1 exhibited as an explicit factor (see Tables I and  II).
In Eq. ( 4) the first matrix element involves the pion decay constant f π and the pion momentum q µ : The second matrix element involves two form factors [8]: The term with f + (q 2 ) vanishes when contracted with q µ , so we only need the longitudinal form factor f 0 (q 2 ) at The value of the form factor according to lattice fits in the physical-continuum limit is [8]: It is now straightforward to obtain the decay rate: where The prediction for the branching ratio using the data above is then: The uncertainty in this branching ratio is dominated by f 2 0 (m 2 π ) [see Eq.( 8)] which is near 3.5%.We have left out the factor a 2 1 as explained above, as it is usually done in the literature.Besides a 1 , there is the uncertainty of neglecting non-factorizable contributions, but these are not expected to be large for a decay such as B c → B s π + , a Class 1 decay [15,14], where only a charged meson with large momentum -in this case π + -is generated directly from the color-singlet current.

Model
Br(Bc → Bsπ) Table 1: Predictions for the branching ratio Br(B c → B s π) from different authors.The values listed are the partial widths extracted from their respective cited article divided by the full width also cited therein.The last line is our calculation based on the lattice result for the form factor as explained in the text.
Our result in Eq. (10) for Br(B + c → B s π + ) is quite similar to several of the other theoretical predictions shown in Table 1, and almost all of them are lower than the LHCb estimate of Eq.( 3) based on the theoretical predictions for Br(B + c → J/ψ π + ).Therefore the LHCb results shown in Eqs. ( 1) and ( 2) reveal a tension in the theoretical predictions of Br(B + c → J/ψ π + ) and Br(B + c → B s π + ).If we trust the result of Eq.( 10) and most predictions in Table 1 for Br(B + c → B s π + ), we conclude that the actual value of Br(B + c → J/ψ π + ) should be smaller than the lowest estimate ∼ 0.06% considered in Ref. [7].
We will review the branching ratio Br(B + c → J/ψ π + ) in Section IV.In the present Section we focus on obtaining f c /f s based on the theoretical estimate for Br(B + c → B s π + ) in Eq. (10), and the LHCb measurement [7] shown in Eq. (1).A straightforward calculation gives the result: where a 1 ≃ 1.2 is the accepted value for c → s transitions [17,18,16,19,14,20,21,22,13].The first uncertainty is from the LHCb measurement [7] shown in Eq. ( 1), which is about ∼ 15%, and the second uncertainty is from the lattice estimate of the form factor [8], which is near 3.5%.Clearly the first uncertainty dominates.In the next Section we examine Concerning the determination of the ratio of fragmentation fractions f s /f d , LHCb has first used B s → D s π vs. B → DK in pp collisions at 7 TeV [9], and lately [10] used several non-leptonic and semileptonic decays of B s and B at 7, 8 and 13 TeV [26,27,9,28,29].In their first study, a ratio of events analogous to Eq.( 2) was obtained: where we have denoted n cor (B s → D s π) the efficiency-corrected number of B s → D s π events and a similar definition for n cor (B 0 → DK).To extract f s /f d from Eq.( 12) one needs independent estimates of the branching ratios.For that purpose, Ref. [9] uses the theoretical expressions based on factorization, where the main theoretical uncertainties are (i) the ratio of the form factors in D s |J µ |B s and D − |J µ |B 0 , and (ii) the neglect of nonfactorizable contributions.Moreover, they use B 0 → D − K + instead of the more abundant mode B 0 → D − π + precisely because the latter mode has an additional diagram with an internal W exchange, as shown in Fig. 1, for which factorization may not be a good approximation [30].Here we consider two other options in this analysis: on the one hand, today we have good experimental data for B 0 → D − π + and B 0 → D − K + so we may try to use the B 0 → D − π + data instead of the theoretical estimates to extract f s /f d from the result in Eq.( 12), provided we have good theoretical estimates for Br(B s → D − s π + ) on the numerator.On the other hand, there are theoretical estimates for the subdominant diagram of B 0 → D − π + (see Fig. 1), so we can verify the uncertainty they bring to the determination of The tree level amplitude for B s → D − s π + has a similar diagram as B c → B s π + (see Fig. 2) and therefore a similar expression as Eq.( 7): but the form factor f (s) 0 (q 2 ) corresponds to the current matrix element D s |( bc)|B s .Similarly the tree level amplitude for B 0 → D − K + is: with the form factor f Let us make a brief parenthesis here to review the significance of non factorizable contributions.All of these two-body nonleptonic weak decays can be treated by means of a low-energy effective weak interaction hamiltonian [15], which at leading order contains two 4-quark operators, O 1 and O 2 , both products of two V − A color singlet currents, the former a product of charged currents and the latter of neutral currents.As an example, the effective interaction hamiltonian for the decay B s → D − s π + is: where is the assumption that these matrix elements can be separated as the product of current matrix elements.
Looking at the flavor content of the mesons, the first matrix element could be approximated as: Since the currents are color singlets, this is expected to be a good approximation.Instead, the matrix element of O 2 requires a Fierz rearrangement, causing a color-suppressed factorized term and a non factorizable one: Notice that, in this case, the factorized matrix element is still O 1 F .This is the case for all three processes shown in Fig. 2, where only one diagram enters.As a result, the whole amplitude within factorization contains only O 1 F : where ) is small and therefore the non-factorizable term is expected to be relatively small.On the right hand side, we follow the notation of Ref. [31] and refer to this colored-allowed amplitude as T .In contrast, B 0 → D − π + (Fig. 1) has two diagrams: one of type T , dominant and similar to the case above, and the other, called E, the internal W exchange, which is subdominant but contains relatively larger non-factorizable contributions.The presence of two diagrams corresponds to the two different contractions of the operators O 1 and O 2 for this process.Using a definition for the factorized matrix element 15), one can see that, for this process: Thus the O 2 F terms add up to a second amplitude, called E [31]: where a 2 = C 2 (µ) + C1(µ) Nc .Here the non factorizable part is not suppressed relative to the term containing O 2 F .Because of this uncertainty, the mode B → DK was preferred instead of B → Dπ.Now, the determination of f s /f d using the method of Eq.( 12) requires a theoretical expression for B s → D s π, since there is no other independent measurement of this rate.Using B → DK instead of B → Dπ in the denominator has the advantage that the theoretical expression for it is similar and several factors tend to cancel, except for the form factors and uncertainties due to the non-factorizable contributions.The non-factorizable contributions are expected to be smaller in B → DK than in B → Dπ, so B → DK was preferred even when the number of events of that mode is smaller.Using Eqs. ( 9), ( 13) and ( 14), one obtains: Here Φ sK is a kinematic factor: is the ratio of form factors at the corresponding kinematic points: while N a is the ratio of non-factorizable corrections.
and N a are the main sources of theoretical uncertainty in Eq.( 19).Now, using Br(B 0 → D − π + ) instead of Br(B 0 → D − K + ) in the denominator would presumably introduce larger uncertainties from non-factorizable contributions.We could avoid them if we use the experimental result of B 0 → D − π + instead, which will be sensitive to the experimental error only, but then we will need the full theoretical estimate of the numerator, which includes the absolute value of the form factor, not just a ratio of form factors: Currently Br(B 0 → D − π + ) exp = (2.51 ± 0.08) × 10 −3 [32], i.e. it is known with a 3.2% uncertainty.Possibly Belle II may provide a measurement with even better precision in the near future.However the theoretical uncertainty in the numerator of Eq. ( 20) is presumably larger.Instead, let us then reevaluate the theoretical uncertainties using Br(B 0 → D − π + ).We will again have a cancellation of form factors (possibly better than in Eq.( 19) by being both evaluated at the same q 2 ), but we need to assess the E amplitude.Compared to Eq.( 19), in this case we have: with the kinematic factor Φ sπ ≃ 1.0179 (using now p (π) CM ≃ 2306.41MeV for B → Dπ), while the ratio of lifetimes [32] τ Bs /τ B 0 = 1.0007 ± 0.0042 is known to a precision better than 0.5%.The ratio of form factors is now: and N E is the correction to include the internal W -exchange diagram in B 0 → D − π + (see Fig. 1), namely where T and E are the color-allowed and W -exchange amplitudes respectively [Eqs.(17) and ( 18)] for the decay B 0 → D − π + .
For the estimate of N E we will use two theoretical results, which actually differ and therefore give us a caveat about the correct estimate of this factor.Let us first consider Ref. [33], which uses perturbative QCD including Sudakov factors to estimate the amplitudes of two-body non-leptonic B → D ( * ) π(ρ) decays.According to these authors, E has comparable factorized and non-factorizable contributions, but both are much smaller than T and mainly imaginary with respect to T : in a normalization relative to T ≡ 1, they find E f act = −0.018+0.057i, E nonf act = −0.0011+0.020i, which they add up to E = −0.019+0.077i.Therefore, the estimate of N E that includes only factorizable terms would be N (f act) E = 1/((1 − 0.018) 2 + (0.057) 2 ) ≃ 1.034, while if the non-factorizable estimates are included the result extracted from Ref. [33] is We could then use N E ≃ 1.033, where the uncertainty due to neglecting the non-factorizable contributions would be well below 1% if those contributions were anything close to the estimated E nonf act [33], and thus negligible compared to other uncertainties in Eq. (21).
Alternatively, let us consider Ref. [23] for the estimation of N E .Here the authors study in systematic detail the different corrections for the factorization hypothesis, including those from the internal Wexchange amplitude E (which they call "annihilation amplitude" A), as well as non-factorizable contributions to the dominant amplitude T .For the first correction they obtain E/T lead ≃ 0.04 and the latter correction T non.f act /T lead ≃ −0.03.Considering all these estimates from Ref. [23], we obtain while disregarding T non.f act the value changes in only 1%.However, the discrepancy in this factor extracted from Refs.[23] and [33] is larger.Taking both results as a range with a central value and an uncertainty, we get: i.e. a 3% uncertainty.Therefore, considering this estimate on N E , still the main source of uncertainty for the ratio of Eq.( 21) would come from the ratio of form factors, N which, according to Ref. [10] (see also [34,35]) it is N (π) F = 1.000 ± 0.042 (a 4% uncertainty).Consequently we would obtain the ratio It is interesting to notice that this result, which is a theoretical estimate of the ratio of branching ratios, slightly differs from the fit from data obtained by LHCb in Ref. [10], which is 1.18 ± 0.04, i.e. a discrepancy in the central values by 2.5 σ, if we consider σ as the quadrature combination of the errors in both results.In order to find out the origin of the discrepancy, we cannot analyze the fit by LHCb, but we can examine the sources of uncertainty of our estimate shown in Eq. (22).
As seen in Eq. ( 21), there are only two possible sources of uncertainty at the level of a few percent or more: the form factor ratio N (π) F , which is estimated to be known to 4% as cited above, and the internal W exchange correction N E which, according to Refs.[23] and [33], should contain an uncertainty around 3%.However, we should be cautious here: there are other estimates that differ from this result.For example, N E ≃ 0.966 ± 0.05 according to Ref. [31].Notice that this value is also below unity, which means that the real part of E should be positive instead of negative as in Ref. [33].But then again, the same Ref.[31] quotes a previous CDF result [36] that leads to a value N E ∼ 1.07 ± 0.03 (stat), i.e larger than unity.In summary, there is still an uncertainty in N E that needs to be resolved.Just considering the range of values N E ∼ 0.966 [31] and 1.033 [33], these differ in near 7%.Now, if we take this more conservative uncertainty of 7%, the theoretical ratio of Eq.( 22) would be known to that level of uncertainty.Then, measurements of B s → D s π + and B 0 → D − π + events at LHCb analogous to that of Eq.( 12) could be used to extract the ratio of fragmentation fractions f s /f d , namely: where n cor are the efficiency-corrected experimental number of events just as in Eq.( 12).One could then extract the value for f s /f d with a precision up to 7%.The actual values of n cor should be determined by experiment, which is of course beyond the scope of our work.
4 On f c /f d and B + c → J/ψ π + As we saw at the end of Section II, there is a tension between the theoretical predictions of Br(B c → B s π) and Br(B c → J/ψ π) if one tries to relate these processes with f c /f s , f c /f d and the experimental data.LHCb [7] determined Br(B c → B s π) by measuring the quantity we reproduce here in Eq. ( 1), and by using f c /f s as additional input.This input was in turn estimated using a previous result of f s /f d and a value of f c /f d obtained from their measurement shown in Eq.( 2) together with a broad range of theoretical estimates for Br(B c → J/ψ π) ∼ 0.06% → 0.18% [13].The resulting branching ratio Br(B c → B s π) ∼ 5% → 15%, in spite of this broad range, is in slight tension with the direct theoretical estimates of Br(B c → B s π) shown in Table I and Eq.(10): in order to have consistency, one should expect Br(B c → J/ψ π) to be lower than the aforementioned range.
Let us then derive the values of Br(B c → J/ψ π) that would be consistent with the theoretical estimates for Br(B c → B s π) given in Section II.To do that, first we take f c /f s in Eq. (11), which was derived using the theoretical estimates for Br(B c → B s π), and take f s /f d = 0.2390 ± 0.0076 obtained by LHCb at 7 TeV [10].One then finds the values for f c /f d : where the first uncertainty comes from the LHCb measurement of Eq.( 1), which is near 15% and by far the largest, the second is near 3.5% and comes from the lattice estimate of the B c → B s form factor shown in Eq.( 8), and the third is near 3.2% and comes from the uncertainty in f s /f d according to Ref. [10].Now, using the LHCb measurement [11] shown in Eq.( 2) and the current value of Br(B + → J/ψ K + ) = (1.020± 0.019) × 10 −3 [32], we get: where the first uncertainty (near 17%) comes from the LHCb measurement shown in Eq.( 2), the second uncertainty (near 15%) is from the LHCb measurement shown in Eq.( 1), and the last uncertainty (near 5%) is the combination of the uncertainties from the branching ratio of B + → J/ψK + , the form factor of Eq.( 8), and f s /f d shown above.The combined uncertainty is near 23%.Our result for Br(B + c → J/ψπ + ) in Eq.( 25) contrasts with the theoretical estimates [13] used by LHCb [7] to extract Br(B + c → B s π + ), namely: Using a 1 = 1.14 (the accepted coefficient a 1 for b → c transitions), our result in Eq.( 25) with its 23% uncertainty is still 3 σ below the lowest value in Eq. ( 26).Moreover, other theoretical models (see Table 2) predict values for Br(B + c → J/ψ π + ) in the range (3.8 → 13) × 10 −4 a 2 1 , and in particular the lowest two values in Table 2 differ in less than 2σ from the value we calculated in Eq. (25).While there are still large uncertainties in the estimation of the B c branching ratios and on the ratio of fragmentation functions, these latter comparisons tend to indicate that our results are approaching a consistent set of values.Table 2: Prediction for the branching ratio Br(B + c → J/ψπ + ) from different models, and our estimate of Eq. 25.The coefficient a 1 = 1.14 has been extracted for comparison between the different results.

Conclusions
We have studied the consistency of current determinations and estimates of the fragmentation fractions f c , f s and f d , which are the probabilities that a b quark produced at the LHCb hadronizes into a B c , B s and B d meson, respectively, together with the predictions of the branching ratios of these mesons into the specific two-body non-leptonic modes that are measured nowadays by LHCb.We compare the theoretical estimates of the branching ratios obtained by different authors using different methods, and in our own estimates we use Lattice results for the corresponding form factors which enter in the theoretical expressions of these decays within the factorization approximation.In the analysis we include the uncertainties due to non-factorizable contributions as well.
In particular, f c is estimated from the theoretical prediction of B + c → B s π + based on lattice results, and measurements by LHCb.This decay is a remarkable case, where the heavy quark is a spectator of the weak interaction process, unlike the much studied decays of heavy hadrons where it is the heavy quark that goes through the weak transition.Current theoretical estimates of Br(B + c → B s π + ) range from 0.01 to 0.05.Our estimate is Br(B + c → B s π + ) ∼ 0.03, based on Lattice results, with an uncertainty of about 4%, in which case f c /f s ∼ 0.08 could be determined from the LHCb measurement with 15% uncertainty.
Concerning the fraction f s /f d , we study the possibility obtaining better results by measuring the ratio Br(B s → D s π)/Br(D 0 → Dπ) instead of the proposed ratio Br(B s → D s π)/Br(D 0 → DK).The latter was proposed because B → Dπ contains non factorizable contributions, thus inducing possibly large uncertainties.However, B → DK has the disadvantage of being a more suppressed mode, so data is more limited.We therefore investigated the effect of non factorizable contributions in B → Dπ to see whether it could compensate the uncertainties by being a more abundant mode.In order to estimate Br(B 0 → Dπ) we have discussed various methods, such as using the experimental value from PDG [32], PQCD factorization predictions [23,33] and the various fittings [31,36].We find that the measurements of B s → D s π and B 0 → Dπ at LHCb can then provide f s /f d up to 7% uncertainty, due to the corresponding non-factorizable contributions.
In summary, we predict Br(B c → B s π) ≃ 0.04 (using a 1 = 1.2 for this c → s transition), within 3.5% uncertainty from the Lattice results for the form factor, which then gives f c /f s ≃ 0.056 with 16% uncertainty using the LHCb measurements of this decay.Concerning f s /f d , we propose to use again the data on B s → D s π and B → Dπ considering our estimate of Br(B s → D s π)/Br(B 0 → D − π + ) ≃ 1.02 ± 0.05, where the uncertainties include those from neglecting non factorizable contributions.Finally, we obtain the ratio f c /f d ≃ 0.014 with a 16% error based on f c /f s and current f s /f d results, and the estimate Br(B c → J/ψπ) ∼ 4 × 10 −4 which is consistent with Br(B c → B s π) and the current data on the fragmentation fractions, with an uncertainty of 23%.This consistent value is at the lower end of the current theoretical estimates, which still span a rather broad range, from 0.04% to 0.17%.form factors and unitarity bounds beyond the SU (3) F limit.Eur.Phys.J. C, 80(4):347, 2020.

2 )Figure 2 :
Figure 2: Tree level diagrams for (a) B + c → B s π + , (b) B s → D − s π + and (c) B 0 → D − K + .All of them are color allowed diagrams in the factorization hypothesis.

O 1 ≡
(ūd)( bc) and O 2 ≡ (ūc)( bd) (factors in parenthesis represent color-singlet V − A currents).C i (µ) are Wilson coefficients at the renormalization scale µ ∼ m b .The matrix elements of these 4-quark operators, D − s π + |O 1 |B s and D − s π + |O 2 |B s , are difficult to calculate because of the strong interactions.Factorization