We present theoretical analysis of paired $B_c$ mesons production in $e^+e^-$ annihilation at different energy scales taking into account full next-to-leading order QCD corrections. Both possible electroweak channels are considered: production via virtual photon and via virtual $Z$-boson. We study in detail the role of radiative QCD corrections, which were found to be significant especially at low energies. It is shown that the contribution from $Z$-boson is significant at high energies ($\sqrt{s} > M_Z/2$) especially in the case of paired production of pseudo-scalar and vector ($B_c + B_c^*$) mesons. Azimuthal asymmetry induced by a $P$-violating weak interaction with $Z$-boson is also analyzed.

In the present paper we study the production of $\chi_b + c\bar c$ at the LHC within single parton scattering approach. A special attention payed to the feed-down from $\chi_b$ states to the associated $\Upsilon + c\bar c$ production, which was recently studied by the LHCb. We have found that this feed-down is about percents of the total cross section seen in the experiment. It is shown that the shapes of the differential distributions are almost same as found by the LHCb except the azimuthal asymmetry, which is almost flat in the experiment. We conclude that the precise study of the single parton scattering contributions is necessary for the correct isolation of other possible production channels such as double parton scattering.

Inclusive radiative decays $B_c\to\gamma u\bar{d}$ are considered. It is shown that photon emission allows one to bypass the chiral suppression and leads the increase of the branching fraction by about four orders of magnitude. The estimates of exclusive decay branching fractions of $B_c$ meson into different sets of $\pi$ mesons are also given.

This report is devoted to the analytical calculation of heavy quarkonia production processes in modern experiments such as LHC, B-factories and superB-factories in computer. Theoretical description of heavy quarkonia is based on the factorization theorem. This theorem leads to special structure of the production amplitudes which can be used to develop computer algorithm which calculates these amplitudes automatically. This report is devoted to the description of this algorithm. As an example of its application we present the results of the calculation of double charmonia production in bottomonia decays and inclusive the $\chi_{cJ}$ mesons production in pp-collisions.

Hadronic production of $P$-wave charmonium states $\chi_{cJ}$ in hadronic interaction is considered. Using experimental results of CDF and LHCb collaborations we show, that contributions of color-singlet components are dominant. As for color-octet mechanism, we show, that contributions from $P$-wave states can also be observed, while $S$-wave states can be neglected. The best experimental observables that give information about relative importance of color-singlet and color-octet components are the ratios $\chi_{c2}/\chi_{c1}$ and $\chi_{c0}/\chi_{c1}$.

The article is devoted to Bc->Bs+n pi, Bc->Bs*+n pi decays with n=1, 2, 3, 4. In the framework of factorization theorem the branching fractions of these processes can be written as convolution of hard part, describing Bc->Bs W, Bc->Bs* W vertices, and spectral functions, that correspond to transition of virtual $W$-boson into a final pi-meson system. These functions were obtained from the fit of experimental data on $\tau$-lepton decay and electron-positron annihilation. Using different sets of Bc->Bs decay form-factors we present branching fractions and distributions over the invariant mass of the final pi-meson system.

This paper is devoted to the study of the leading twist distribution amplitudes of P-wave nonrelativistic mesons. It is shown that at the leading order approximation in relative velocity of quark-antiquark pair inside the mesons these distribution amplitudes can be expressed through one universal function. As an example, the distribution amplitudes of the P-wave charmonia mesons are considered. Within QCD sum rules the model for the universal function of P-wave charmonia mesons is built. In addition, it is found the relations between the moments of the universal function and the nonrelativistic QCD matrix elements that control relativistic corrections to any amplitude involving P-wave charmonia. Our calculation shows that characteristic size of these corrections is of order of ~30 %.

The production of the charmonium states in hadronic experiments is considered in NLO at the strong coupling constant. It is shown, that such an approach solves some significant problems, arising when only leading order processes are considered. In particular, in such a consideration distributions over the transverse momentum of the final charmonium can be obtained. There appears also a natural explanation to the existence of $\chi_{c1}$-meson in final state, that is observed experimentally and cannot be produced in leading order processes.

In this paper we show that heavy quarckonium states, similar to light mesons, form Regge trajectories in $(n,M^2)$ and $(M^2,J)$ planes and the slope of these trajectories is independent on the quantum numbers of the mesons. This fact can be useful for the prediction of the masses of heavy quarkonia and the determination of the quantum numbers of the newly discovered states.

In this paper we have analyzed the production of pair charmonium mesons in the reactions $e^+ e^- \to J/ \Psi \chi_{c0}, e^+ e^-\to \Psi(2S) \chi_{c0}$ at energy $\sqrt s = 10.6$GeV in the framework of the light cone formalism. In comparison with NRQCD the numerical results for the cross sections are in better agreement with experiment.

In the framework of scalar-vector dominance we calculate the hadronic matrix elements of scalar and tensor effective quark currents induced by virtual leptoquark interactions. Combined bounds on the product of couplings and leptoquark masses are obtained from experimental data.

We speculate on a possible interpretation of events selected by the SELEX collaboration and stress an insufficient evidence for the observation of doubly charmed baryon.

This report provides a comprehensive overview of the prospects for B physics at the Tevatron. The work was carried out during a series of workshops starting in September 1999. There were four working groups: 1) CP Violation, 2) Rare and Semileptonic Decays, 3) Mixing and Lifetimes, 4) Production, Fragmentation and Spectroscopy. The report also includes introductory chapters on theoretical and experimental tools emphasizing aspects of B physics specific to hadron colliders, as well as overviews of the CDF, D0, and BTeV detectors, and a Summary.

We present the calculation of the inclusive $x_F$-distributions of charmed hadrons, produced in high-energy $\Sigma^-$-beam. The calculation is based on the modified mechanism of charmed quarks fragmentation as well as on the mechanism of $c$-quark recombination with the valence quarks from initial hadrons. We predict the additional asymmetry in the production of charmed hadrons due to the different distributions of the valence s and d quarks in $\Sigma^-$-beam.

In the framework of the vector dominance model (VDM) new contribution into the cross-section of $D^*$-meson photoproduction at HERA has been estimated. This contribution is due to the interaction of the virtual vector $c\bar c$-mesons from photon with the proton. It has been shown that the mechanism under discussion plays the essential role for all values of transverse momentum of $D^*$-mesons at HERA. Taking into account for the contribution of this mechanism improves the description of the experimental cross section distribution.

We present the calculation of the inclusive x_F-distributions of charmed hadrons, produced in high-energy $\Sigma^-$-beam. Our calculation is based on the modified mechanism of charmed quarks fragmentation as well as on the mechanism of c-quark recombination with the valence quarks from initial hadrons.

We perform a detailed investigation of total lifetimes for the doubly heavy baryons $\Xi_{QQ'}$, $\Omega_{QQ'}$ in the framework of operator product expansion over the inverse heavy quark mass, whereas, to estimate matrix elements of operators obtained in OPE, approximations of nonrelativistic QCD are used.

Using the quark-diquark approximation in the framework of Buchm\" uller-Tye potential model, we investigate the spectroscopy of doubly charmed baryons: $\Xi_{cc}^{++}$ and $\Xi_{cc}^{+}$. Our results include the masses, parameters of radial wave functions of states with the different excitations of both diquark and light quark-diquark system. We calculate the values of fine and hyperfine splittings of these levels and discuss some new features, connected to the identity of heavy quarks, in the dynamics of hadronic and radiative transitions between the states of these baryons.

Estimates of lifetimes and partial branching ratios for the baryons $\Xi_{bc}^{+}$ and $\Xi_{bc}^{0}$ are presented using the inverse heavy quark mass expansion technique carried out in the Operator Product Expansion approach. We take into account both the perturbative QCD corrections to the spectator contributions and, depending on the quark contents of hadrons, the Pauli interference and weak scattering effects between the constituents, using a potential model for the evaluation of the non-perturbative parameters.

We have performed a detailed investigation of total lifetimes for the $\Xi_{cc}^{++}$ and $\Xi_{cc}^{+}$ baryons in the framework of operator product expansion over the inverse mass of charmed quark, whereas, to estimate matrix elements of operators obtained in OPE, some approximations of nonrelativistic QCD are used. This approach allows one to take into account the corrections to the spectator decays of $c$-quarks, which reflect the fact, that these quarks are bound, as well as the contributions, connected to the effects of both the Pauli interference for the $\Xi_{cc}^{++}$-baryon and the weak scattering for the $\Xi_{cc}^{+}$-baryon. The realization of such program leads to the following estimates for the total lifetimes of doubly charmed baryons: $\tau_{\Xi_{cc}^{++}}=0.43\pm 0.1~ps$ and $\tau_{\Xi_{cc}^{+}}=0.11\pm 0.01~ps$.

We propose a modification of the mechanism of charmed quarks fragmentation into $D$ mesons in hadronic collisions. It is shown that the distinction in valence quark distributions in the initial $\pi^{\pm}$ and $K^{\pm}$ mesons leads to different inclusive spectra of $D$ and $D_s$ mesons produced in $\pi^{\pm}$ and $K^{\pm}$ beams.

In the leading $O(\alpha_s^4)$ order of the perturbative QCD, the hadronic production cross-sections of $S$- and $P$-wave states of $B_c$ meson are calculated. The results for the $S$-wave levels are compared with the values given in other papers as well as in the model of the $b$ quark fragmenatation into $B_c$. In the given order, the cross-sections of the hadronic production of the $P$-wave states are calculated for the first time. Their contribution into the $B_c$ meson production is less than 10\%. There is a strong difference between the predictions of the fragmentation model and the exact perturbative calculations. These difeerences are discussed in details for the differential distributions over various kinematical quantities.

The cross-section for the production of $b{\bar b} c{\bar c}$ quarks in $e^+e^-$ annihilation, that proves to be at a level of $\sigma (e^+e^- \rightarrow b{\bar b} c{\bar c})/ \sigma (e^+e^- \rightarrow b{\bar b}) \sim 10^{-2}$ for ${\sqrt s}=M_Z$ is calculated within the frames of the QCD perturbation theory. The cross sections for the associated production of $1S$- and $2S$-wave states of $B_c$-meson in the reaction $e^+e^- \rightarrow B_c b{\bar c}$ were calculated in the nonrelativistic model of a heavy quarkonium. The fragmentation function of $b\to B_c^{(*)}$ is analysed in the scaling limit. The number of $\Lambda_{bc}$-hyperons to be expected at LEP is estimated on the basis of the assumption on quark-hadron duality.

In this paper we summarize the results of the theory working group dedicated to the analysis of $B_c$ production at hadron colliders.