
We present theoretical analysis of paired $B_c$ mesons production in $e^+e^$
annihilation at different energy scales taking into account full
nexttoleading 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 pseudoscalar 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
feeddown from $\chi_b$ states to the associated $\Upsilon + c\bar c$
production, which was recently studied by the LHCb. We have found that this
feeddown 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, Bfactories and
superBfactories 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 ppcollisions.

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 colorsinglet components are
dominant. As for coloroctet 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 colorsinglet and coloroctet 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 pimeson system. These functions were obtained
from the fit of experimental data on $\tau$lepton decay and electronpositron
annihilation. Using different sets of Bc>Bs decay formfactors we present
branching fractions and distributions over the invariant mass of the final
pimeson system.

This paper is devoted to the study of the leading twist distribution
amplitudes of Pwave nonrelativistic mesons. It is shown that at the leading
order approximation in relative velocity of quarkantiquark pair inside the
mesons these distribution amplitudes can be expressed through one universal
function. As an example, the distribution amplitudes of the Pwave charmonia
mesons are considered. Within QCD sum rules the model for the universal
function of Pwave 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 Pwave 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 scalarvector 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 highenergy $\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 crosssection 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_Fdistributions of charmed
hadrons, produced in highenergy $\Sigma^$beam. Our calculation is based on
the modified mechanism of charmed quarks fragmentation as well as on the
mechanism of cquark 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 quarkdiquark approximation in the framework of Buchm\" ullerTye
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 quarkdiquark 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 nonperturbative 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 crosssections 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 crosssections 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 crosssection 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 quarkhadron duality.

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