Abstract:We calculate the cross sections of associated $J/\psi + \psi^\prime$ and $J/\psi + J/\psi$ production in $pp$ collisions at $\sqrt s = 13$ TeV in the forward kinematic region. The High Energy Factorization ($k_T$-factorization) framework supplemented with the Catani-Ciafaloni-Fiorani-Marchesini evolution of gluon densities in a proton is applied. We demonstrate that latest data on $J/\psi + J/\psi$ production and first experimental data on $J/\psi + \psi^\prime$ events taken very recently by the LHCb Collaboration can be described well by the color singlet terms and contributions from the double parton scattering (DPS) with the standard choice for $\sigma_{\rm eff}$ parameter. The relative production rate $\sigma(J/\psi + \psi^\prime)/\sigma(J/\psi + J/\psi)$ is found to be sensitive to the DPS terms as well as to feeddown contributions.

What problem does this paper attempt to address?

This paper aims to solve the matching problem between the theoretical description and experimental data of the production processes of the doublets J/ψ+ψ′ and J/ψ+J/ψ in high - energy proton - proton collisions, especially in the forward kinematic region. Specifically, the authors use the High Energy Factorization (HEF) framework, combined with the Catani - Ciafaloni - Fiorani - Marchesini (CCFM) - evolved gluon density, to calculate the cross - sections of associated production in these processes. ### Main problems: 1. **Matching of theory and experimental data**: The paper attempts to verify the latest LHCb collaboration experimental data through theoretical calculations, especially the data on J/ψ+ψ′ and J/ψ+J/ψ production. These data include the distributions of multiple momentum variables, such as the rapidity difference ($\Delta y$), the azimuthal angle difference ($\Delta \phi$), the transverse momentum ($p_T$), the rapidity ($y$) and the invariant mass ($M$). 2. **Contributions of single - parton scattering (SPS) and double - parton scattering (DPS)**: The paper explores the relative contributions of the single - parton scattering and double - parton scattering mechanisms to these processes. In particular, it studies the importance of the DPS mechanism in the low - $p_T$ and large - invariant - mass regions. 3. **Comparison of color - singlet (CS) and color - octet (CO) mechanisms**: The paper analyzes the main role of the color - singlet mechanism in the forward - rapidity region and evaluates the contribution of the color - octet mechanism. ### Research methods: - **High Energy Factorization (HEF) framework**: Use the kT - factorization method, combined with the CCFM evolution equation, to calculate the transverse - momentum dependence of the gluon density. - **Numerical simulation**: Through numerical simulation, calculate the differential cross - sections under different momentum variables and compare them with the LHCb experimental data. - **Parameter selection**: Consider different TMD gluon density distributions (such as JH’2013 set 2 and LLM’2022), as well as different sources of theoretical uncertainties, such as scale uncertainties and quark - mass variations. ### Conclusions: - **Consistency between theory and experiment**: The theoretical predictions of the paper are generally consistent with the LHCb experimental data, especially in the low - invariant - mass and medium - transverse - momentum regions. - **Importance of the DPS mechanism**: In certain specific regions, such as the low - $p_T$ and large - invariant - mass regions, the contribution of the DPS mechanism is significant and has an important impact on the quality of data description. - **Dominant role of the color - singlet mechanism**: In the forward - rapidity region, the color - singlet mechanism is the main contributor, while the contribution of the color - octet mechanism is small. Through these studies, the paper not only verifies the effectiveness of the existing theoretical models but also provides an important reference for further exploring the strong interactions in high - energy collisions.

Forward $J/ψ+ J/ψ$ and $J/ψ+ ψ^\prime$ production with High Energy Factorization