We present an efficient way to calculate the effect of soft QCD radiation atone loop, which is needed for predictions at next-to-next-to-leadinglogarithmic accuracy. We use rapidity coordinates and isolate the divergencesin the integrand. By performing manipulations with cumulative variables, weavoid complications from plus distributions. We address rapidity divergences,divergences with an azimuthal dependence, complicated jet boundaries andmulti-differential measurements. The process and measurements can be easilyadjusted, as we demonstrate by reproducing many existing soft functions. Theresults for a general LHC process with multiple Wilson lines are obtained bytreating Wilson lines that are not back-to-back using a boost. We also obtain,for the first time, the N-jettiness soft function for generic jet angularities,and the collinear-soft function for the measurement of two angularities.