The azimuthal correlation, Δ ϕ12, of high transverse momentum jets in pp collisions at s=13 TeV is studied by applying PB-TMD distributions to NLO calculations via MCatNLO together with the PB-TMD parton shower. A very good description of the cross section as a function of Δ ϕ12 is observed. In the back-to-back region of Δ ϕ12→ π, a very good agreement is observed with the PB-TMD Set 2 distributions while significant deviations are obtained with the PB-TMD Set 1 distributions. Set 1 uses the evolution scale while Set 2 uses transverse momentum as an argument in αs, and the above observation therefore confirms the importance of an appropriate soft-gluon coupling in angular ordered parton evolution. The total uncertainties of the predictions are dominated by the scale uncertainties of the matrix element, while the uncertainties coming from the PB-TMDs and the corresponding PB-TMD shower are very small. The Δ ϕ12 measurements are also compared with predictions using MCatNLO together Pythia8, illustrating the importance of details of the parton shower evolution. © 2022, The Author(s).