We discuss peculiarities that arise in the computation of real-emission contributions to observables that contain Heaviside functions. A prominent example of such a case is the zero-jettiness soft function in SCET, whose calculation at next-to-next-to-next-to-leading order in perturbative QCD is an interesting problem. Since the zero-jettiness soft function distinguishes between emissions into different hemispheres, its definition involves θ-functions of light-cone components of emitted soft partons. This prevents a direct use of multi-loop methods, based on reverse unitarity, for computing the zero-jettiness soft function in high orders of perturbation theory. We propose a way to bypass this problem and illustrate its effectiveness by computing various non-trivial contributions to the zero-jettiness soft function at NNLO and N3LO in perturbative QCD.

We present the calculation of the next-to-next-toleading order (NNLO) zero-jettiness beam and soft functions, up to the second order in the expansion in the dimensional regularization parameter. These higher order terms are needed for the computation of the next-to-next-to-next-toleading order (N 3 LO) zero-jettiness soft and beam functions. As a byproduct, we confirm the O 0 results for NNLO beam and soft functions available in the literature by Gaunt et al.

We present an analytic calculation of the one-loop correction to the double-real emission contribution to the zero-jettiness soft function at N3LO in QCD, accounting for both gluon-gluon and quark-antiquark soft final-state partons. We explain all the relevant steps of the computation including the reduction of phase-space integrals to master integrals in the presence of Heaviside functions, and the methods we employed to compute them.

We present a calculation of all matching coefficients for N-jettiness beam functions at next-to-next-to-next-to-leading order (N3LO) in perturbative quantum chromodynamics (QCD). Our computation is performed starting from the respective collinear splitting kernels, which we integrate using the axial gauge. We use reverse unitarity to map the relevant phase-space integrals to loop integrals, which allows us to employ multi-loop techniques including integration-by-parts identities and differential equations. We find a canonical basis and use an algorithm to establish non-trivial partial fraction relations among the resulting master integrals, which allows us to reduce their number substantially. By use of regularity conditions, we express all necessary boundary constants in terms of an independent set, which we compute by direct integration of the corresponding integrals in the soft limit. In this way, we provide an entirely independent calculation of the matching coefficients which were previously computed in ref. [1].