We consider the quark beam function that describes collinear initial-state radiation that is constrained by a veto on reconstructed jets. As the veto is imposed on the transverse momenta of the jets, the beam function is subject to rapidity divergences, and we use the collinear-anomaly framework to extract the perturbative matching kernels to next-to-next-to-leading order (NNLO) in the strong-coupling expansion. Our calculation is based on a novel framework that automates the computation of beam functions in Mellin space and it provides the ingredients to extend jet-veto resummations for quark-initiated processes to NNLL′ accuracy.
