Two-loop $N$-jettiness soft function for $pp \to 2j$ production

We calculate the two-loop soft function for $pp \to jj$, $e^+e^- \to 4j$ or $e p \to 3j$ when the $N$-jettiness observable is measured. The result presented here makes up the necessary piece for realizing the full next-to-next-to-leading order predictions for these processes using the $N$-jettiness subtraction scheme.