
The cross section of associated production of a Z boson with heavy flavor
jets in $pp$ collisions is calculated using the SHERPA Monte Carlo generator
and the analytical combined QCD approach based on ktfactorization at small x
and conventional collinear QCD at large x. A satisfactory description of the
ATLAS and CMS data on the $p_T$ spectra of Z bosons and cjets in the whole
rapidity, y, region is shown. Searching for the intrinsic charm (IC)
contribution in these processes, which could be visible at large y > 1.5, we
study observables very sensitive to nonzero IC contributions and less affected
by theoretical QCD scale uncertainties. One of such observables is the
socalled double ratio: the ratio of the differential cross section of Z + c
production in the central region of y < 1.5 and in the forward region 1.5 <
y < 2.5, divided by the same ratio for Z + b production. These observables
could be more promising for the search of IC at LHC as compared to the
observables considered earlier.

Constraints on the intrinsic charm probability $w = P_{c\bar c / p}$ in the
proton are obtained for the first time from LHC measurements. ATLAS
Collaboration data for the production of prompt photons, accompanied by a
charmquark jet in pp collisions at $\sqrt s = 8 $ TeV are used. Fits to the
ATLAS data result in a central estimate: $w = 1.14\%$. Upper limits: $w <
2.74\%$ ($4.32\%$) are obtained at the $68\% (95\%) $ C.L. The results,
however, do not rule out zero intrinsic charm at one standard deviation. We
suggest methods which can reduce the uncertainties for determining the
intrinsic charm probability, such as the elimination of theoretical scale
ambiguities and tests of $x_T = 2 E_T^\gamma / \sqrt{s}$  scaling at fixed
rapidity, as well as measurements of the $x_T$ distribution at different LHC
energies. Other implications of intrinsic heavy quarks in the proton for QCD
studies at the LHC are also discussed.

We consider the production of $J/\psi$ and $\psi(2S)$ mesons originating from
the decays of $b$flavored hadrons at the LHC using the $k_T$factorization
approach. Our analysis covers both inclusive charmonia production and
production of $J/\psi$ mesons in association with $Z$ bosons. We apply the
transverse momentum dependent (or unintegrated) gluon density in a proton
derived from CataniCiafaloniFioraniMarchesini (CCFM) evolution equation and
adopt fragmentation functions based on the nonrelativistic QCD factorization
to describe the inclusive $b$hadron decays into the different charmonium
states. Our predictions agree well with latest experimental data taken by the
CMS, ATLAS and LHCb Collaborations at $\sqrt s = 7$, $8$ and $13$~TeV. The
contributions from double parton scattering to the associated nonprompt
$J/\psi + Z$ production are estimated and found to be small.

We consider the production of $Z$ bosons associated with beauty quarks at the
LHC in the framework of $k_T$factorization approach. Our consideration is
based on the offshell gluongluon fusion subprocess $g^* g^* \to Z Q\bar Q$ at
the leading order ${\cal O}(\alpha\alpha_s^2)$ (where the $Z$ boson further
decays into a lepton pair) and several subleading ${\cal O}(\alpha \alpha_s^2
)$ and ${\cal O}(\alpha \alpha_s^3 )$ subprocesses involving quarkantiquark
and quarkgluon interactions, taken into account in the conventional
(collinear) QCD factorization. The contributions from double parton scattering
are discussed as well. The transverse momentum dependent (or unintegrated)
gluon densities in a proton are derived from
CataniCiafaloniFioraniMarchesini (CCFM) evolution equation. We achieve
reasonably good agreement with the latest data taken by CMS and ATLAS
Collaborations. The comparison of our results with nexttoleadingorder pQCD
predictions, obtained in the collinear QCD factorization, is presented. We
discuss the uncertainties of our calculations and demonstrate the importance of
subleading quark involving contributions in describing the LHC data in the
whole kinematic region.

We investigate the inclusive Higgs boson production in protonproton
collisions at the CERN LHC conditions using the $k_T$factorization approach.
Our analysis is based on the dominant offshell gluongluon fusion subprocess
(where the transverse momenta of initial gluons are taken into account) and
covers $H \to \gamma \gamma$, $H \to ZZ^* \to 4l$ (where $l = e,\mu$) and $H
\to W^+ W^ \to e^\pm \mu^\mp \nu \bar \nu$ decay channels. The transverse
momentum dependent (or unintegrated) gluon densities in a proton were derived
from CiafaloniCataniFioraniMarchesini equation or, alternatively, were
chosen in accordance with KimberMartinRyskin prescription. We estimate the
theoretical uncertainties of our calculations and compare our results with
nexttonexttoleadingorder plus nexttonexttoleadinglogarithmic ones
obtained using collinear QCD factorization. Our predictions agree well with the
latest experimental data taken by the CMS and ATLAS Collaborations at $\sqrt s
= 8$ and $13$~TeV.

We investigate the inclusive Higgs boson production in protonproton
collisions at high energies in the framework of ktfactorization QCD approach.
The attention is focused on the dominant offshell gluongluon fusion
subprocess, where the transverse momentum of incoming gluons are taken into
account. The transverse momentum dependent (or unintegrated) gluon densities of
the proton are determined using the CCFM evolution equation as well as the
KimberMartinRyskin prescription. We study the theoretical uncertainties of
our calculations and perform the comparison with the results of traditional
pQCD evaluations. Our predictions agree well with the first experimental data
taken by the ATLAS collaboration at the LHC. We argue that further studies of
the Higgs boson production are capable of constraining the unintegrated gluon
densities of the proton.

In the framework of $k_T$factorization approach, the production and
polarization of prompt $J/\psi$ mesons at the LHC energies is studied. Our
consideration is based on the nonrelativistic QCD formalism for bound states
and offshell amplitudes for hard partonic subprocesses. Both the direct
production mechanism and feeddown contributions from $\chi_c$ and $\psi(2S)$
decays are taken into account. The transverse momentum dependent (or
unintegrated) gluon densities in a proton were derived from
CiafaloniCataniFioraniMarchesini evolution equation or, alternatively, were
chosen in accordance with KimberMartinRyskin prescription. The
nonperturbative coloroctet matrix elements were first deduced from the fits
to the latest CMS data on $J/\psi$ transverse momentum distributions and then
applied to describe the ATLAS and LHCb data on $J/\psi$ production and
polarization at $\sqrt s = 7$, $8$ and $13$ TeV. We perform an estimation of
polarization parameters $\lambda_\theta$, $\lambda_\phi$ and $\lambda_{\theta
\phi}$ which determine $J/\psi$ spin density matrix and demonstrate that
treating the soft gluon emission as a series of explicit colorelectric dipole
transitions within NRQCD leads to unpolarized $J/\psi$ production at high
transverse momenta, that is in qualitative agreement with the LHC data.

We consider an observable very sensitive to the nonzero intrinsic charm (IC)
contribution to the proton density. It is the ratio between the differential
cross sections of the photon or $Z$boson and $c$jet production in the $pp$
collision, $\gamma(Z) + c$, and the $\gamma(Z)$ and the $b$jet production. It
is shown that this ratio can be approximately flat or increasing at large
$\gamma(Z)$ transverse momenta $p_T$ and their pseudorapidities $1.5 < \eta <
2.4$ if the IC contribution is taken into account. On the contrary, in the
absence of the IC this ratio decreases as $p_T$ grows. We also present the
ratios of the cross sections integrated over $p_T$ as a function of the IC
probability $w$. It is shown that these ratios are mostly independent on the
theoretical uncertainties, and such predictions could therefore be much more
promising for the search for the intrinsic charm signal at the LHC compared to
the predictions for $p_T$spectra, which significantly depend on these
uncertainties.

We consider the production of $b\bar b$ quarks and DrellYan lepton pairs at
LHC conditions focusing attention on the total transverse momentum of the
produced pair and on the azimuthal angle between the momenta of the outgoing
particles. Plotting the corresponding distributions in bins of the final state
invariant mass, one can reconstruct the full map of the transverse momentum
dependent parton densities in a proton. We give examples of how can these
distributions can look like at the LHC energies.

We reconsider the inclusive production of isolated prompt photons in $pp$
collisions at the LHC energies in the framework of $k_T$factorization
approach. Our analysis is based on the ${\cal O}(\alpha \alpha_s)$ offshell
(depending on the transverse momenta of initial quarks and gluons) production
amplitudes of $q^* g^* \to \gamma q$ and $q^* \bar q^* \to \gamma g$ partonic
subprocesses and transverse momentum dependent (or unintegrated) quark and
gluon densities in a proton, which are chosen in accordance with the
KimberMartinRyskin prescription. We show that the subleading highorder
${\cal O}(\alpha \alpha_s^2)$ contributions, not covered by the noncollinear
evolution of parton densities, are important to describe latest LHC data.

The production of weak gauge bosons in association with heavy flavored mesons
at the LHCb conditions is considered, and a detailed study of the different
contributing processes is presented including single (SPS) and double (DPS)
parton scattering mechanisms. We find that the usual DPS factorization formula
needs to be corrected for the limited partonic phase space, and that including
the relevant corrections reduces discrepancies in the associated $ZD$
production. We conclude finally that double parton scattering dominates the
production of samesign $W^\pm D^\pm$ states, as well as the production of
$W^$ bosons associated with $B$mesons. The latter processes can thus be
regarded as new useful DPS indicators.

We study the role of the nonperturbative input to the transverse momentum
dependent (TMD) gluon density in hard processes at the LHC. We derive the input
TMD gluon distribution at low scale mu0^2 ~ 1 GeV^2 from the fit of the
inclusive hadron spectra measured at low transverse momenta in pp collisions at
the LHC and demonstrate that the best description of these spectra for larger
hadron transverse momenta can be achieved by matching the derived TMD gluon
distribution with the exact solution of the BalitskyFadinKuraevLipatov
(BFKL) equation obtained at low x and small gluon transverse momenta outside
the saturation region. Then, we extend the input TMD gluon density to higher
mu^2 numerically using the CataniCiafoloniFioraniMarchesini (CCFM) gluon
evolution equation. A special attention is put to the phenomenological
applications of obtained TMD gluon density to some LHC processes, which are
sensitive to the gluon content of a proton.

In the framework of ktfactorization approach,the production of prompt
$\chi_c$ mesons in pp collisions at the LHC energies is studied. Our
consideration is based on the offshell amplitudes for hard partonic
subprocesses $g^*g^*\to\chi_{cJ}$ and nonrelativistic QCD formalism for bound
states. The transverse momentum dependent (unintegrated) gluon densities in a
proton were derived from CiafaloniCataniFioraniMarchesini evolution equation
or, alternatively, were chosen in accordance with KimberMartinRyskin
prescription. Taking into account both color singlet and color octet
contributions, we deduce the corresponding nonperturbative longdistance
matrix elements from the fits to the latest ATLAS data on $\chi_{c1}$ and
$\chi_{c2}$ transverse momentum distributions at $\sqrt s = 7$ TeV. We find
that these distributions at small and moderate pt are formed mainly by the
color singlet components. We successfully described the data on the
differential cross sections and relative production rates
$\sigma(\chi_{c2})/\sigma(\chi_{c1})$ presented by the ATLAS, CMS and LHCb
Collaborations. We find that the fit points to unequal wave functions of
$\chi_{c1}$ and $\chi_{c2}$ states.

In the framework of ktfactorization approach, the production and
polarization of prompt psi(2S) mesons in pp collisions at the LHC energies is
studied. Our consideration is based on the nonrelativistic QCD formalism for
bound states and offshell amplitudes for hard partonic subprocesses. The
transverse momentum dependent (unintegrated) gluon densities in a proton were
derived from CiafaloniCataniFioraniMarchesini evolution equation or,
alternatively, were chosen in accordance with KimberMartinRyskin
prescription. The nonperturbative coloroctet matrix elements were first
deduced from the fits to the latest CMS data on psi(2S) transverse momentum
distributions and then applied to describe the ATLAS and LHCb data on psi(2S)
production and polarization at 7 TeV. We perform the estimation of polarization
parameters $\lambda_\theta$, $\lambda_\phi$ and $\lambda_{\theta \phi}$ which
determine psi(2S) spin density matrix and demonstrate that taking into account
the offshellness of initial gluons in the coloroctet contributions leads to
unpolarized psi(2S) production at high transverse momenta, that is in
qualitative agreement with the LHC data.

We investigate the possibility to detect the scalar Higgs boson decay $H\to
b\bar b$ in the associated $Z$ and $b\bar b$ production at the LHC using the
$k_T$factorization QCD approach. Our consideration is based on the offshell
(i.e. depending on the transverse momenta of initial quarks and gluons)
production amplitudes of $q^* \bar q^* \to Z H \to Z q^\prime \bar q^\prime$,
$q^* \bar q^* \to Z q^\prime \bar q^\prime$ and $g^* g^* \to Z q^\prime \bar
q^\prime$ partonic subprocesses supplemented with the
CataniCiafoloniFioraniMarchesini (CCFM) dynamics of parton densities in a
proton. We argue that the $H \to b\bar b$ signal could be observed at large
transverse momenta near Higgs boson peak despite the overwhelming QCD
background, and point out an important role of angular correlations between the
produced $Z$ boson and $b$quarks.

Associated production of charged gauge bosons and charged charmed mesons at
the LHC is considered in the framework of $k_T$factorization approach.
Theoretical predictions are compared with ATLAS data, and reasonably good
agreement is found. Predictions on the samesign $W^{\pm}D^{\pm}$
configurations are presented including single parton scattering and double
parton scattering contributions. The latter are shown to dominate over the
former, thus giving evidence that the proposed process can serve as another
indicator of double parton interactions.

We consider the transverse momentum dependent (TMD) quark densities of the
proton which are very important ingredients for unpolarized DrellYan (DY)
lepton pair production. We calculate the TMD sea quark density as a convolution
of the CataniCiafaloniFioraniMarchesini (CCFM)evolved gluon distribution
and the TMD gluontoquark splitting function. Based on the ${\cal
O}(\alpha^2)$ production amplitude $q^* + \bar q^* \to Z/\gamma^* \to l^+ +
l^$, calculated by taking into account the effective $q^*\bar q^* Z/\gamma^*$
 vertex, we analyze the distributions on the dilepton invariant mass,
transverse momentum and rapidity and specific angular correlations between the
produced leptons as measured by the CMS, ATLAS and LHCb collaborations. We
briefly duscuss also the process of assiciated lepton pair and jet production
in $pp$ collisions at the LHC.

We study the associated production of real (isolated) or virtual photons
(with their subsequent leptonic decay) and hadronic jets in protonproton
collisions at the LHC using the $k_T$factorization approach of QCD. The
consideration is based on the offshell quarkgluon QCD Compton scattering
subprocesses. In the case of virtual photon production, the contributions from
Z boson exchange as well as gamma*Z interference with the full spin
correlations are included. The transverse momentum dependent (TMD) quark and
gluon densities in a proton are determined from the KimberMartinRyskin
prescription or CataniCiafoloniFioraniMarchesini (CCFM) equation. In the
latter, we restricted to the case where the gluontoquark splitting occurs at
the last evolution step and calculate the sea quark density as a convolution of
the CCFMevolved gluon distribution and the TMD gluontoquark splitting
function. Our numerical predictions are compared with the recent experimental
data taken by the ATLAS Collaboration. We discuss the theoretical uncertainties
of our calculations and argue that further studies are capable of constraining
the TMD parton densities in a proton.

We use the TMD quark densities of the proton to investigate unpolarized
DrellYan lepton pair production in protonproton collisions at the LHC
energies. We investigate the case where the gluontoquark splitting occurs at
the last evolution step and calculate the TMD sea quark density as a
convolution of the CCFMevolved gluon distribution and the TMD gluontoquark
splitting function which contains all single logarithmic smallx corrections to
the sea quark evolution for any order of perturbation theory. Based on the
O(alpha) production amplitude of quarkantiquark annihilation calculated
according to the reggeized quark approach, we analyze the distributions on the
dilepton invariant mass, transverse momentum and rapidity as well as the
specific angular correlations between the produced leptons as measured by the
CMS, ATLAS and LHCb collaborations. We argue that these measurements impose
stringent constraints on the TMD quark distributions of the proton.

We study the unintegrated, or transverse momentum dependent (TMD) gluon
distribution obtained from the best description of the LHC data on the
inclusive spectra of hadrons produced in the midrapidity region and low
transverse momenta at starting scale Q0^2 = 1 GeV^2. To extend this gluon
density at higher Q^2 we apply the CataniCiafoloniFioraniMarchesini (CCFM)
evolution equation. The influence of the initial (starting) nonperturbative
gluon distribution is studied. The application of the obtained gluon density to
the analysis of the ep deep inelastic scattering allows us to get the results
which describe reasonably well the H1 and ZEUS data on the longitudinal proton
structure function FL(x,Q^2). So, the connection between the soft processes at
LHC and small x physics at HERA has been confirmed and extended to a wide
kinematical region.

We reconsider the prompt photon photoproduction at HERA in the framework of
the $k_T$factorization QCD approach. The proposed method is based on the
${\cal O}(\alpha^2\alpha_s)$ amplitudes for $\gamma q\to\gamma gq$ and $\gamma
g^*\to\gamma q\bar q$ partonic subprocesses. Additionally, we take into account
the ${\cal O}(\alpha^2\alpha_s^2)$ box contributions $\gamma g\to\gamma g$ to
the production cross sections. The unintegrated (or transverse momentum
dependent) parton densities in the proton are determined using
KimberMartinRyskin (KMR) prescription. Our consideration covers both
inclusive and jet associated prompt photon photoproduction rates. We find that
our numerical predictions agree well with the recent data taken by H1 and ZEUS
collaborations at HERA. We demonstrate that the box contributions are sizeable
and amount up to $\sim15%$ of the calculated total cross section.

We find out the connection between the unintegrated gluon distribution at low
intrinsic transverse momenta and the inclusive spectrum of the hadrons produced
in pp collision at LHC energies in the midrapidity region and low hadron
transverse momenta. The parameters of this distribution are found from the best
description of the LHC data. Its application to the analysis of ep deep
inelastic scattering allows us to obtain the results which reasonable describe
the H1 and ZEUS data on the structure functions at low x. The connection
between the soft processes at LHC and smallx physics at HERA is found.

In the framework of the k_Tfactorization approach, the production of prompt
photons in association with a heavy (charm or beauty) quarks at high energies
is studied. The consideration is based on the O(\alpha \alpha_s^2) offshell
amplitudes of gluongluon fusion and quark(anti)quark interaction
subprocesses. The unintegrated parton densities in a proton are determined
using the KimberMartinRyskin prescription. Our numerical predictions are
compared with the D0 and CDF experimental data. Also we extend our results to
LHC energies.

In the framework of the ktfactorization approach, the isolated prompt photon
pair production in pp and p \bar p collisions at high energies is studied. The
consideration is based on the quarkantiquark annihilation, quarkgluon
scattering and gluongluon fusion subprocesses, where the nonzero transverse
momenta of incoming partons are taken into account. The unintegrated quark and
gluon densities in a proton are determined using the KimberMartinRyskin
prescription. The numerical analysis covers the total and differential
production cross sections and extends to specific angular correlations between
the produced prompt photons. Theoretical uncertainties of our evaluations are
studied and comparison with the NLO pQCD calculations is performed. The
numerical predictions are compared with the recent experimental data taken by
the D0, CDF, CMS and ATLAS collaborations at the Tevatron and LHC energies.

In the framework of the ktfactorization approach, the production of prompt
photons in association with a heavy (charm or beauty) quarks at high energies
is studied. The consideration is based on the O(\alpha \alpha_s^2) offshell
amplitudes of gluongluon fusion and quark(anti)quark interaction
subprocesses. The unintegrated parton densities in a proton are determined
using the KimberMartinRyskin prescription. The analysis covers the total and
differential cross sections and extends to specific angular correlations
between the produced prompt photons and muons originating from the semileptonic
decays of associated heavy quarks. Theoretical uncertainties of our evaluations
are studied and comparison with the results of standard NLO pQCD calculations
is performed. Our numerical predictions are compared with the recent
experimental data taken by the D0 and CDF collaborations at the Tevatron.
Finally, we extend our results to LHC energies.