In this paper, we continue to concentrate on the quark-antiquark pair (or dijet) production in a special process that a virtual photon scattering off a hadronic target, this process typically corresponds to the quark-antiquark (or dijet) production in DIS. Due to a broad range of kinematics that can be accessed in the future EIC, this process allows to mapping out the evolutions of gluon TMDs. Especially, when the scattered energy is sufficiently large (or equivalent to small-\(x\) ), a high density regime can be probed, in this sense, gluon saturation effects will become more predominant, and hence the behaviors of gluon TMDs in the dense regime can be investigated as well. In the correlation limit, prior works have shown that this process can be consistently described in the framework of Color Glass Condensate (CGC) theory, and simultaneously, the relevant two WW gluon distributions can be expressed in terms of Wilson line correlators that making them can be calculated by using the solutions of Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner (JIMWLK) equation [1]}, [2]}, [3]}, [4]}. However, the special limit has a strong restriction on the dijets, which requires the transverse momentum imbalance \(q_\perp \)  is much smaller than the relative transverse momentum \(P_\perp \) , in this case, there is another class of effects, known as Sudakov resummations, can give sizable contribution to the cross section, and this should be taken into account when exploring the behaviors of the two gluon distributions in the considered process. Other phenomenological studies have shown that such resummations is very important in quantitatively describing experimental observables [5]}, [6]}, [7]}, [8]}, [9]}, [10]}, [11]}, [12]}, [13]}, [14]}, [15]}, [16]}, [17]}. Therefore, in the present work, we plan to conduct this direction by including the Sudakov resummations, and show how the resummations affect the azimuthal correlations of the dijets. Furthermore, as shown in Refs. [3]}, [4]}, [20]}, [21]}, thanks to the numerical implementation of JIMWLK equation, the two gluon distributions appeared in the considered process can be obtained. However, the numerical implementation is only achieved by fixing coupling constant. The running-coupling corrections, as an important QCD dynamics, has not been performed yet in their studies. In the present work, we prefer to work within the Gaussian approximation of JIMWLK evolution and use the running-coupling Balitsky-Kovchegov (rcBK) equation to take the relevant calculations, because we expect the running-coupling corrections here to be much more important than that of Gaussian approximation. Besides that, we stress the rcBK solutions are well constrained from DIS data [22]} and have been received a great successful in describing particle spectrum at hadronic collisions [23]}, [24]}, therefore, one can be firmly convinced that the gluon TMDs can be well studied by using the solutions of rcBK equation.
