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 kt-factorization 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 c-jets 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 non-zero IC contributions and less affected by theoretical QCD scale uncertainties. One of such observables is the so-called 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 charm-quark 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 Catani-Ciafaloni-Fiorani-Marchesini (CCFM) evolution equation and adopt fragmentation functions based on the non-relativistic 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 non-prompt $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 off-shell gluon-gluon 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 quark-antiquark and quark-gluon 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 Catani-Ciafaloni-Fiorani-Marchesini (CCFM) evolution equation. We achieve reasonably good agreement with the latest data taken by CMS and ATLAS Collaborations. The comparison of our results with next-to-leading-order 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 proton-proton collisions at the CERN LHC conditions using the $k_T$-factorization approach. Our analysis is based on the dominant off-shell gluon-gluon 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 Ciafaloni-Catani-Fiorani-Marchesini equation or, alternatively, were chosen in accordance with Kimber-Martin-Ryskin prescription. We estimate the theoretical uncertainties of our calculations and compare our results with next-to-next-to-leading-order plus next-to-next-to-leading-logarithmic 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 proton-proton collisions at high energies in the framework of kt-factorization QCD approach. The attention is focused on the dominant off-shell gluon-gluon 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 Kimber-Martin-Ryskin 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 non-relativistic QCD formalism for bound states and off-shell amplitudes for hard partonic subprocesses. Both the direct production mechanism and feed-down 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 Ciafaloni-Catani-Fiorani-Marchesini evolution equation or, alternatively, were chosen in accordance with Kimber-Martin-Ryskin prescription. The non-perturbative color-octet 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 color-electric 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 non-zero 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 pseudo-rapidities $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 Drell-Yan 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)$ off-shell (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 Kimber-Martin-Ryskin prescription. We show that the sub-leading high-order ${\cal O}(\alpha \alpha_s^2)$ contributions, not covered by the non-collinear 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 same-sign $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 non-perturbative 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 Balitsky-Fadin-Kuraev-Lipatov (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 Catani-Ciafoloni-Fiorani-Marchesini (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 kt-factorization approach,the production of prompt $\chi_c$ mesons in pp collisions at the LHC energies is studied. Our consideration is based on the off-shell amplitudes for hard partonic subprocesses $g^*g^*\to\chi_{cJ}$ and non-relativistic QCD formalism for bound states. The transverse momentum dependent (unintegrated) gluon densities in a proton were derived from Ciafaloni-Catani-Fiorani-Marchesini evolution equation or, alternatively, were chosen in accordance with Kimber-Martin-Ryskin prescription. Taking into account both color singlet and color octet contributions, we deduce the corresponding non-perturbative long-distance 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 kt-factorization 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 non-relativistic QCD formalism for bound states and off-shell amplitudes for hard partonic subprocesses. The transverse momentum dependent (unintegrated) gluon densities in a proton were derived from Ciafaloni-Catani-Fiorani-Marchesini evolution equation or, alternatively, were chosen in accordance with Kimber-Martin-Ryskin prescription. The non-perturbative color-octet 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 off-shellness of initial gluons in the color-octet 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 off-shell (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 Catani-Ciafoloni-Fiorani-Marchesini (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 same-sign $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 Drell-Yan (DY) lepton pair production. We calculate the TMD sea quark density as a convolution of the Catani-Ciafaloni-Fiorani-Marchesini (CCFM)-evolved gluon distribution and the TMD gluon-to-quark 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 proton-proton collisions at the LHC using the $k_T$-factorization approach of QCD. The consideration is based on the off-shell quark-gluon 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 Kimber-Martin-Ryskin prescription or Catani-Ciafoloni-Fiorani-Marchesini (CCFM) equation. In the latter, we restricted to the case where the gluon-to-quark splitting occurs at the last evolution step and calculate the sea quark density as a convolution of the CCFM-evolved gluon distribution and the TMD gluon-to-quark 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 Drell-Yan lepton pair production in proton-proton collisions at the LHC energies. We investigate the case where the gluon-to-quark splitting occurs at the last evolution step and calculate the TMD sea quark density as a convolution of the CCFM-evolved gluon distribution and the TMD gluon-to-quark splitting function which contains all single logarithmic small-x corrections to the sea quark evolution for any order of perturbation theory. Based on the O(alpha) production amplitude of quark-antiquark 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 mid-rapidity 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 Catani-Ciafoloni-Fiorani-Marchesini (CCFM) evolution equation. The influence of the initial (starting) non-perturbative 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 Kimber-Martin-Ryskin (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 mid-rapidity 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 small-x physics at HERA is found.

In the framework of the k_T-factorization 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) off-shell amplitudes of gluon-gluon fusion and quark-(anti)quark interaction subprocesses. The unintegrated parton densities in a proton are determined using the Kimber-Martin-Ryskin 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 kt-factorization 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 quark-antiquark annihilation, quark-gluon scattering and gluon-gluon fusion subprocesses, where the non-zero transverse momenta of incoming partons are taken into account. The unintegrated quark and gluon densities in a proton are determined using the Kimber-Martin-Ryskin 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 kt-factorization 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) off-shell amplitudes of gluon-gluon fusion and quark-(anti)quark interaction subprocesses. The unintegrated parton densities in a proton are determined using the Kimber-Martin-Ryskin 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.