The theoretical description of Drell-Yan (DY) transverse momentum spectra over wide kinematic regions in energy, mass and transverse momentum requires not only fixed-order perturbative QCD calculations but also soft-gluon QCD resummations to all orders of perturbation theory. The latter are traditionally accomplished either by Parton Showers (PS) with Monte Carlo event generators or by (different versions of) analytical procedures. In this work we focus on issues involved in the matching of the fixed-order calculation and resummation, especially in the moderate to low mass and $p_{\bot}$ region. In particular we address the DY region accessible at fixed target experiments. We present a Parton Branching (PB) formulation in which transverse momentum dependent (TMD) evolution is matched with MCatNLO calculations of NLO matrix elements. Using this formulation, we show a good theoretical description of DY data from experiments in very different kinematic ranges: NuSea, R209, Phenix, LHC $8\;\textrm{TeV}$ and $13\;\textrm{TeV}$.
Comment: Contribution to the 2021 QCD session of the 55th Rencontres de Moriond, 4 pages, 6 figures