Total: 1
We present a complete analytic calculation of the next-to-leading-order (NLO) QCD cross sections for the production of an electroweak vector boson ($W^\pm$, $Z/γ^*$) at finite transverse momentum, decomposed into the full set of helicity cross sections that determine the angular coefficients $A_0$-$A_7$ of the lepton-pair angular distribution in the Collins-Soper frame. The calculation includes the unpolarised cross section and all eight helicity projections for both neutral- and charged-current Drell-Yan production, the axial quark-triangle contributions in the $m_t\to\infty$ limit and with full top-quark-mass dependence, and closed-form expressions for all coefficient functions. The results are implemented in the DYTurbo program and validated at the level of a few times $10^{-4}$ against its numerical implementation derived from MCFM, as well as through the analytic $q_T\to0$ limit. The analytic implementation achieves a speed-up of more than five orders of magnitude with respect to the numerical code, enabling fast and numerically stable predictions for precision electroweak measurements at the LHC. The calculation was carried out with extensive use of large language models (Anthropic Claude) as autonomous agents; we document the methodology and the associated computational and token budgets.