In a series of papers, we have investigated the compatibility of the Kimber-Martin-Ryskin (KMR) and Martin-Ryskin-Watt (MRW) unintegrated parton distribution functions (UPDFs) as well as the description of the experimental data on the proton structure functions. The present work is a sequel to that survey, via calculation of the transverse-momentum distribution of the electroweak gauge vector bosons in the kt-factorization scheme, by the means of the KMR, the leading-order (LO) MRW, and the next-to-leadingorder (NLO) MRW UPDF, in the NLO. To this end, we have calculated and aggregated the invariant amplitudes of the corresponding involved diagrams in the NLO and counted the individual contributions in different frameworks. The preparation process for the UPDF utilizes the parton distribution functions of Martin et al., MSTW2008-LO, MSTW2008-NLO, MMHT2014-LO, and MMHT2014-NLO, as the inputs. Afterward, the results have been analyzed against each other as well as the existing experimental data, i.e., D0, CDF, ATLAS, and CMS collaborations. Our calculations show excellent agreement with the experiment data. It is, however, interesting to point out that the calculation using the KMR framework illustrates a stronger agreement with the experimental data, despite the fact that the LO MRW and the NLO MRW formalisms employ a better theoretical description of the Dokshitzer- Gribov-Lipatov-Altarelli-Parisi evolution equation. This is of course due to the use of the different implementation of the angular ordering constraint in the KMR approach, which automatically includes the resummation of lnð1=xÞ, Balitski-Fadin-Kuraev-Lipatov logarithms, in the LO Dokshitzer-Gribov- Lipatov-Altarelli-Parisi evolution equation.
