Bottomonium spectrum with a Dirac potential model in the momentum space

We study the bottomonium spectrum using a relativistic potential model in the momentum space. This model is based on a complete one gluon exchange interaction with a momentum dependent screening factor to account for the effects due to virtual pair creation that appear close to the decay thresholds. The overall model does not make use of nonrelativistic approximations. We fit well established bottomonium states below the open bottom threshold and predict the rest of the spectrum up to ≈11200\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\approx 11200$$\end{document} MeV and JPC=3--\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J^{PC}=3^{--}$$\end{document}. Uncertainties are treated rigorously and propagated in full to the parameters of the model using a Monte Carlo to identify if which deviations from experimental data can be absorbed into the statistical uncertainties of the models and which can be related to physics beyond the bb¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$b\bar{b}$$\end{document} picture, guiding future research. We get a good description of the spectrum, in particular the Belle measurement of the ηb(2S)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta _b(2S)$$\end{document} state and the Υ(10860)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varUpsilon (10860)$$\end{document} and χb(3P)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\chi _b(3P)$$\end{document} resonances.


Introduction
The heavy quark meson sector constitutes a major piece of information on the nonperturbative regime of the strong interaction. In particular, a lot of experimental information has been gathered on the bottomonium spectrum during the last years thanks to ATLAS, BaBar, Belle, BESIII, CLEO, CMS, D0, and LHCb collaborations [1][2][3][4][5][6][7][8][9][10][11][12] and further results are expected in the near future during the Belle II run [13,14] and after the CMS and LHCb upgrades [15,16]. Theory work has preceded and followed through the experimental a e-mail: djmolinab@unal.edu.co (corresponding author) b e-mail: mdesanctis@unal.edu.co c e-mail: cesar.fernandez@nucleares.unam.mx d e-mail: santopinto@ge.infn.it effort [17][18][19][20][21] in the form of Lattice QCD computations [22][23][24][25][26][27][28][29][30][31][32][33], Dyson-Schwinger-Bethe-Salpeter equations [34][35][36][37][38][39][40], effectives Hamiltonians reduced from QCD [41,42] and potential quark models [43][44][45][46][47][48][49][50][51][52][53][54][55][56]. In this paper we develop a relativistic quark model for bottomonia based on a complete one gluon exchange. The approach is completely relativistic and does not rely on nonrelativistic approximations. In this way the standard spin-orbit, spin-spin, and tensor interactions are automatically included. We also incorporate a relativistic scalar interaction and a momentum dependent screening factor to account for the effects due to virtual pair creation that appear close to the decay thresholds. All the calculations are performed in the momentum space. The same model was successfully applied to reproduce the charmonium spectrum in Ref. [57] which we refer the reader to for technicalities. We fit the model to all the known states of each J PC below the BB threshold except for the recently measured χ b1 (3P) and χ b2 (3P) which we prefer to predict in order to gain insight on their nature and the η b (2S) which we exclude of our fit owing to the disagreement between CLEO [58] and Belle [59] measurements. We perform a rigorous error estimation that allows us to assess if the inclusion of a new effect in the phenomenological model is necessary or not, and we compute the parameter correlations which provide insight on how independent are the different pieces of the model among them. A full error analyses is mandatory to identify which deviations from experimental data can be absorbed into the statistical uncertainties of the models and which can be related to physics beyond the bb picture, guiding future research.
The paper is organized as follows: in Sect. 2 we provide the relativistic quark model and the employed solution method; in Sect. 3 we describe the fitting procedure as well as the statistical method used to compute the uncertainties; in Sect. 4 we report the computed bottomonium spectrum up to J PC = 3 −− and ≈ 11200 MeV as well as the comparison to the available experimental information. We obtain a very good description of both fitted and nonfitted bottomonia and also predict many unobserved states; Sect. 5 summarizes the conclusions.

Hamiltonian model
We apply to bottomonia the same model developed in [57] for charmonia. The total interaction Hamiltonian of this model is given by the sum of a vector (H (v) ) and a scalar (H (s) ) term. In the momentum space, it has the form: where p a and p b represent the three-momenta of both quark and antiquark in the center of mass of the bottomonium system. In more detail, the vector interaction is constructed in the Coulomb gauge taking into account the whole one gluon exchange contribution. We have: The effective potential function V (v) (q) will be introduced and discussed in the following Eq. (7a). Furthermore, represents the standard Dirac four-current of the quarks, σ i stands for the Pauli matrices, and γ μ i are the gamma matrices, where i = 1, 2 is the particle label. We also introduce the quark on-shell energy difference (being E(p) = p 2 + m 2 ) and the squared (positive) four momentum transfer where q = p b − p a represents the three momentum transfer.
The factors 1 ± (ΔE) 2 Q 2 give the retardation contributions in the Coulomb gauge. Analogously, the scalar interaction is defined as where I i is a scalar vertex and V (s) (q) is the scalar effective potential function of the following Eq. (7b). Finally, the vector and the scalar effective potentials functions have the following form: The first term of Eq. (7a) represents a standard vector interaction, mediated by a massless gluon; α st is the effective strong coupling constant. When inserted in Eq. (2), the complete Coulomb gauge relativistic interaction (with the color factor 4/3) is obtained. The second term of Eq. (7a) represents a regularized confining interaction corresponding to a linear term in the coordinate space, as shown in Ref. [57], in which β v corresponds to the vector confinement strength. As for the scalar potential function given by, Eq. (7b) we note that it contains a phenomenological constant term A plus a β s term which corresponds to the scalar confinement strength.
In both vector and scalar confining terms the regularization parameter b has been introduced to avoid the divergence when |q| → 0. With the definitions given above, the complete interaction Hamiltonian of Eq. (1) can be straightforwardly constructed. In any case, all the details have been given in the Appendix of Ref. [57].
As in [57] for charmonia, we use two different prescriptions for the scalar interaction: potential I → model using Eqs. (7) with β s = 0, potential II → model using Eqs. (7) with β s = 0. (8) In this way we can check if the two forms of the scalar interaction (with or without the confinement term) have the same effect on the spectrum, as in the case of the charmonium system. Besides, in order to take into account the effects of the virtual [60,61] pair creation that appear close to the decays thresholds, we include a screening momentum dependent factor. Hence, the total Hamiltonian takes the final form where the screening factors F s ( p) with p a or p b is defined as Where p s and k s are others parameter of the model. In this way, the model, with potential I and potential II, depends on seven and eight parameters, respectively (see Table 2).

Relativistic equation and solution method
The relativistic equation we use has been obtained performing a three dimensional reduction of the Bethe-Salpeter equation and keeping only the contributions of the positive energy Dirac spinors [57]. In the center of mass of the bb system, the relativistic integral equation takes the form where we have introduced the relativistic energy and M 0 represents the phenomenological zero point energy of the spectrum, M is the resonance mass (i.e. the eigenvalue of the integral equation) and Ψ (p) is the resonance wave function. The wave function Ψ n,{ν} (p) ({ν} = L , S, J ) can be written as where R n,L ( p;p) corresponds to the radial function in the momentum space with n the principal quantum number,p the variational parameter (with dimensions of momentum), Y L ,M L (p) are the spherical harmonics, and χ S,M S is the spin function. To solve Eq. (11) we use the variational method. As trial functions we use a combination of a finite subset of three dimensional harmonic oscillators. Hence, we can write the Hamiltonian matrix as The eigenvalues and the eigenstates are found through the variational method, diagonalizing and minimizing the M {ν},n b ,n a matrix in Eq. (14) [57,62]. The angular part is solved analytically and the radial part numerically. The details can be found in the Appendix of Ref. [57].

Parameter determination
To determine the values of the parameters, the uncertainties, and the theoretical bottomonium spectrum we fit the experimental masses given in Table 1, i.e. all the known states of each J PC below the BB threshold except for the recently measured χ b1 (3P) and χ b2 (3P) which we try to predict in order to gain insight on their nature and the η b (2S) which we prefer to exclude of our fit owing to the disagreement between CLEO and Belle masses [51]. From CLEO data a mass of 9974.6 ± 2.3 ± 2.1 MeV [58] was obtained while Belle measures 9999.0 ± 3.5 +2.8 −1.9 MeV [59]. BABAR reports a range value between 9974 and 10,015 MeV [63]. The PDG favors the Belle measurement [12], therefore we show this experimental value in Table 1 and Figs. 1 and 2. To perform the fits and the error analyses we use the bootstrap technique [64][65][66][67] and proceed as follows: 1. We randomly choose values for the masses of the resonances by sampling a Gaussian distribution according to their uncertainties (systematic and statistical summed in quadrature), obtaining a resampled bottomonium spectrum; 2. We use the least-squares method to minimize the squared distance where M i are the resampled experimental bottomonia, i.e. the states 0 . The E th i represents the theoretical states calculated by solving the eigenvalue Eq. (14) with potentials I and II. The fit is performed using MINUIT [68].
This procedure is repeated 1000 times in order to obtain enough statistics to compute the expected values of the parameters as well as their uncertainties at a 1σ (68%) confidence level (CL). The expected value of the parameters ( Table 2) are computed as the mean value of the 1000 samples. The uncertainties are obtained as the the differences between the mean value and the highest and lowest masses of the best 68% of the fits. Hence, our uncertainties can be asymmetric. Once the parameters have been determined, we can compute the bottomonium spectrum and the associated uncertainties (Table 1). We find an excellent agreement between theory and fitted states within uncertainties. We note that the values of the common parameters of the two potentials are very similar. These results show that, unlike for charmonia, the scalar confinement term of the interaction does not seem to be relevant in the bottomonia description. To gain further insight on this issue we compute the correlation matrices, Tables 3 and 4, for potentials I and II, respectively.
For potential I (Table 3) we find a strong correlation between the parameters of the vector interaction (α s and β v ) and the scalar interaction parameter A, which indicates that vector and scalar interactions are physically correlated in this model. The screening parameter p s is weakly correlated with the vector interaction parameters but strongly correlated with the scalar interaction ones. For potential II, we have the additional parameter β s . In this potential, the parameters are less correlated as shown in Table 4 with one exception, the additional scalar interaction parameter β s is noticeble correlated with the vector interaction parameter β v . Consequently, we find a significant correlation between the confinement terms of the vector and the scalar interactions. The parameter p s of the screening factors is weakly correlated with the other parameters of the interactions except with the phenomenological parameter A in the scalar interaction. This sizeable correlation highlights how the screening factor impacts more on the scalar interaction.
Using the values obtained in the fitting procedure we plot the screening function F s ( p) in Fig. 3 for the two potentials.    As mentioned above we introduce the screening momenta p j 1/2 (j = I, II labels potentials I and II) which are given by F s ( p j 1/2 ) = 1/2 (we recall that F s (0) = 1). Through the fitting values, we find p I 1/2 = 3.38 GeV and p II 1/2 = 3.34 GeV. These values correspond to the screening kinetic energȳ which amount toĒ I = 11.281 GeV for potential I and E II = 11.260 GeV for potential II. This result show that the screening effect is active above the open bottom threshold as in charmonia. Nevertheless, due to the high values of E I,II , we find that the screening effect is less relevant for the low-lying part of the bottomonium spectrum than for charmonia [57].

Bottomonium spectrum
Using the relativistic model interaction, with either potential I or II, we obtain the bottomonium spectrum. Through the bootstrap method, the errors in the fitted states are carried in full to the computed uncertainties in the parameters and to the spectrum. We provide the computed spectrum in Tables 1 (fitted states) and 5 (predicted states). The computed and the experimental spectra are compared in Figs. 1 (potential I) and 2 (potential II). In general, the spectrum is reproduced by the model within the experimental uncertainties. We note that the parameters obtained with both potentials are very similar, leading to closely akin spectra. This result shows that the confining part of the scalar potential does not impact the bottomonium spectrum. However, the presence of the scalar interaction is necessary for an optimum fit, i.e. the parameter A contribution in Eq. (7b). In what follows we look into the states that were not included in the fit as well as the predicted higher-lying spectrum.
As a general remark, we note that the obtained spectrum is similar to the results of nonrelativistic and semirelativistic calculations [43][44][45][46][47][48][49][50][51][52][53][54][55][56]. In this concern, it can be argued that the high value of the bottom quark mass tends to reduce the relativistic effects with respect, for example, to the charmonium case. But also, the values of the effective parameters of the nonrelativistic and semirelativistic calculations can simulate some relativistic effects. Furthermore, the different dynamical details of the various models do not allow for straightforward and definitive conclusions about this point.

Υ (4S), Υ (10860) and Υ (11020)
These resonances belong to the family with quantum numbers 1 −− . They were discovered by means of e + e − collisions in the mid-eighties [70,71] and were more recently measured by the Belle collaboration [72]. The Υ (4S) is regarded as a 4 3 S 1 state; its experimental mass is M Υ (4S) = 10579.4 ± 1.2 MeV and is not well reproduced by either potential I or II. This resonance is generally considered as a bb state, but its mass is overestimated by models that make use of different approaches: e.g., the nonrelativistic model in Ref. [43] provides M Υ ( 4S) 10630 MeV, the semirelativistic model of Ref. [51] finds M Υ (4S) = 10607 MeV, and the non-relativistic coupled channels model in Ref. [45] reports M Υ (4S) = 10603 MeV. Our computations provides approximately 10642 ± 40 MeV, with both potentials. This result is compatible with the other models, but far away from the experimental value, even when the uncertainties are taken into account. Consequently, our result combined with nonrelativistic calculations suggest that there must be beyond the qq picture effects that need to be included to properly describe the state.
The Υ (10860) resonance is generally interpreted as a Υ (5 3 S 1 ), e.g. in [43,45,49,51,52]. However, the theoretical calculations for the pion emission decay widths, to Υ (1S), Υ (2S) and Υ (3S) are two orders of magnitude [48] greater than the measurement [73] leading to different possible interpretations, such as that Υ (10860) is a mixing of a standard Υ (5S) with a P hybrid state [74], Finally, in Ref. [44] this state is interpreted as a Υ (6S), and, hence, the Υ (5S) becomes a missing resonance of the experimental spectrum In our model, this mass state can be reproduced as a Υ (5S) (5 3 S 1 ) (see Table 5; Figs. 1 and 2) or as a 4 3 D 1 state with both potentials. We do not find support the Υ (6S) interpretation. Actually, our predicted mean value mass, with Υ (5S), is only 5 (1) MeV away from the experimental value with potential I (II). Consequently, we identify this state as a bb with Υ (5S) quantum numbers. However, any final conclusion requires the explanation of the before mentioned pion emission decay widths which we leave for a future work.
Finally, the high excitation Υ (11020) state has been successfully described in standard calculations as a bb meson in a 6 3 S 1 state. We recall, in particular, Refs. [50,51]. For completeness, we also mention that in Ref. [44] it is interpreted as a 7 3 S 1 state. Our model gives a slightly higher mass value to the 6 3 S 1 state, with either potential. Our results are similar to that of the relativized calculation of Ref. [43]. At the moment, no further conclusions can be drawn for this state within the present model.

χ b (3P) states
The χ b (3P) states have been the focus of several experimental collaborations during the last years. An estimation of the χ b (3P) system barycenter (i.e. spin-weighted mass average of the χ b0 (3P), χ b1 (3P), and χ b2 (3P) states) was reported by ATLAS [1] and D0 [2] collaborations, yielding 10,530 ± 5(stat) ± 9(syst) MeV and 10551 ± 14(stat) ± 17(syst) MeV, respectively. More recently, two out of the three state masses were measured; χ b1 (3P) by the LHCb collaboration obtaining 10515.7 +2.2 −3.9 (stat) +1.5 −2.1 (syst) MeV, and χ b1 (3P) and χ b2 (3P) by the CMS collaboration [7] yielding 10,513.42 ± 0.41(stat) ± 0.18(syst) MeV and 10,524.02±0.57(stat)±0.18(syst) MeV, respectively. Several predictions of these states are available in the literature, employing different frameworks. For example, in Ref. [60] a mass of 10,524 MeV is predicted for the χ b1 (3P) state employing a screened potential; in Ref. [45], 10,517 MeV for the same state by means of a coupled channel calculation; and 10,580 MeV in the unquenched quark model [46]. All of the results overestimate the mass of χ b1 (3P). In our calculation (which porpously does not fit this state) we obtain 10540 ± 30 MeV with both potentials whose central value also overestimates the mass of the state. When the uncertainties are taken into account, the experimental mass falls within our error bars and no indication of the need for additional physics is called for. This shows how important it is to perform a rigorous error estimation when performing a level-bylevel comparison between theory and experiment, as differences that can be accounted by the error analyses can be mistaken by physics beyond the bb picture. Regarding χ b2 (3P), 10,532.4 MeV is obtained in Ref. [45] using the coupled channels formalism and 10,578 MeV under the unquenched quark model [46]. We obtain 10554 +25 −28 and 10,557 +22 −42 with potentials I and II, respectively. The CMS value falls well within our uncertainites for potential I and slightly out of them for potential II, although certainly within 2σ uncertainties. Hence, the individually measured χ b (3P) states are well reproduced by our model. Finally, we obtain the barycenter mass 10,545 +24 −27 MeV for potential I and 10,549 +23 −41 MeV for potential II, both compatible with the previously quoted ATLAS and D0 estimations. Recalling that not all the individual states of the χ bJ (3P) system have been measured, we provide in Tables 6 (potential I) and 7 (potential II) the n = 1, 2, 3 barycenter masses, given by [75,76] M n P = M χ b0 (n P) + 3M χ b1 (n P) + 5M χ b2 (n P) 9 , along with the available experimental measurements and estimates from PDG values Given that both potentials produce similar spectra, the χ b (n P) barycenters are very similar. In summary, we find a good agreement, within errors, between the models and the experimental barycenters. Finally, we would like to mention that it has been theorized that some of the states in the χ b (3P) system could be the bottomonia counterparts of the X (3872) charmonium [47,77], i.e. states closely related to the opening if the BB, B B * , and B s B s thresholds. Our results do not support such hypothesis, as the model reproduces the χ b (3P) system within (large) uncertainties, contrary to the X (3872) case which was overestimated using the same model [57], and whose description (both mass and width) calls for additional dynamics beyond the cc picture. Along the same ideas, according to Ref. [78], the χ b1 (4P) state could significantly couple to the BB * and B * B * channels. The measurement of this particular state combined with the comparison to quark model calculations, like the one presented in this work, can provide insight on Table 5 Predicted bottomonia for potentials I and II compared to the existing experimental masses with their corresponding uncertainties. Notation as in Table 1 Name n 2S+1 L J Mass (MeV) Potential I Potential II Experiment  Table 6 Theoretical results obtained, using Potential I, for the states χ bJ (n P) compared with the available experimental data; n = 1, 2, 3 is the principal quantum number;M n stands for the barycenter of the system for each n.
The experimental states for n = 1, 2 are taken from Ref. [7].   Table 7 Same as in Table 6 with the results given by Potential II  Table 8 Differences ΔS n = n 3 S − n 1 S. We observe that these differences decrease when n is increased. All the differences reported in this

Missing resonances
Besides reproducing the experimentally established states, in Table 5 we provide predictions of states both above and below the open bottom thresholds (≈ 10.6 GeV). In total, we predict 38 states up to 11.3 GeV for 0, 1, 2 (with either ± combinations for P and C) and 3 −− quantum numbers. These predicted states are of interest for future analyses at LHCb [15,79,80] and Belle II [13,14,[81][82][83]. In particular, pinning down the Υ (6S) would provide further insight on bottomonium-like states [83]. The missing η b (nS) sector (n 1 S 0 states) can be studied through their relation to their angular momentum partners Υ (nS) (n 3 S 1 ) -known from experiment-, by computing the ΔS n = n 3 S − n 1 S mass splitting. This difference should decrease as n increases in the potential model context [84]. The experimental data for ΔS 1 and ΔS 2 shown in Table 8 support this theoretical results. Thereby, we consider our mass estimations for both η b (nS) and Υ (nS) reasonable.
We also provide predictions for states of the n 1,3 D 1,2,3 family, which remain undetected except for the 1 3 D 2 resonance. The predicted missing states (with uncertainties) provide useful information to guide the forthcoming spectroscopy programs in Belle II [13,14] and LHC [15,16]. However, the production rate of these states should be low, hence, difficult to detect [79].

Conclusions
We have developed a relativistic quark model in momentum space to study the bottomonium spectrum. The model closely follows the one used in Ref. [57] to study charmonium. It combines vector and scalar interactions with a momentum dependent screening factor to account for the effects due to virtual pair creation that appear close to the decay thresholds. We fitted our model to all the known states of each J PC below the BB threshold except for the recently measured χ b1 (3P) and χ b2 (3P) which we prefer to predict in order to gain insight on their nature and the η b (2S) which we exclude of our fit owing to the disagreement between CLEO and Belle measurements. Our prediction for η b (2S) mass agrees with the Belle result.
We have performed a full statistical error analyses using the bootstrap technique, that provides a rigorous treatment of the statistical uncertainties. In this way we obtain the uncertainties of the parameters and their correlations and we can propagate both to the predicted spectrum. Previous error analyses within phenomenological models have been very limited and incomplete. The rigorous error estimations allow us to assess if the inclusion of a new effect in the phenomenological model is necessary or not, and the correlations provide insight on how independent are the different pieces of the model among them. A full error analyses is mandatory to identify which deviations from experimental data can be absorbed into the statistical uncertainties of the models and which can be related to physics beyond the bb picture, guiding future research. We find that the model reproduces very well the fitted states as well as the nonfitted ones within uncertainties.
To asses the importance of a confining term in the scalar interaction, i.e. β s = 0 in Eq. (7b), we fitted the data with and without such contribution. The results obtained with the two potentials are very similar for the fitted and the predicted states, both in the low and the high parts of the spectrum. Therefore, such confining contribution to the scalar interaction can be disregarded in a bottomonium relativistic model. Even so, the correlations found among the parameters belonging to the scalar interaction and the rest of the model parameters, show that the scalar interaction A in Eq. (7b) is strictly necessary to reproduce the spectrum. The screening factor F s ( p) included in the interaction Hamiltonian begins to impact the predictions in significant way at ≈ 11200 MeV, i.e. further away from the open bottom decay thresholds. Hence, the screening effect is not particularly intense and has a slight impact on the bottomonium spectrum, contrary to what it was found for the charmonium one [57].
We have also studied the χ b (3P) resonances. In particular we have calculated the mass of each state of this system and its barycenter. The experimental mass value of the χ b1 (3P) falls into the theoretical uncertainty calculated with both potentials. Whereby, we conclude that the model is able to properly predict this state. Also, the model, with both potentials, reproduces the χ b1 (3P) state. Our result indicates that the χ b1,2 (3P) states are more likely to be bb mesons than the hypothetical X b states.
Our model overestimates the Υ (4S) mass and is consistent with results obtained by semirelativistic quark models, within errors. This is an indication of physics beyond the bb picture for this state. We identify the Υ (10860) as a 5 3 S 1 state. The model fails to reproduce accurately the Υ (11020), although it was well reproduced by nonrelativistic calculations. Further investigation, within our model, should be performed for this high excitation state.
Finally, we report some states that, up to now, have not been observed experimentally but the confirmation of their existence is part of the experimental plans at LHC B factories and Belle II.