Quite recently, two sets of new experimental data from the LHCband the CMSCollaborationshave been published, concerning the production of the Z0vector boson in hadron–hadron collisions with the center-of-mass energy ECM=√s=13 TeV. On the other hand, in our recent work, we have conducted a set of semi-NLOcalculations for the production of the electro-weak gauge vector bosons, utilizing the unintegrated parton distribution functions (UPDF) in the frameworks of Kimber–Martin–Ryskin(KMR) or Martin–Ryskin–Watt(MRW) and the kt-factorization formalism, concluding that the results of the KMRscheme are arguably better in describing the existing experimental data, coming from D0, CDF, CMSand ATLASCollaborations. In the present work, we intend to follow the same semi-NLOformalism and calculate the rate of the production of the Z0vector boson, utilizing the UPDFof KMRwithin the dynamics of the recent data. It will be shown that our results are in good agreement with the new measurements of the LHCband the CMSCollaborations