Two loop soft function for secondary massive quarks

We present the calculation of the O(α2sCFTF) massive quark corrections to the soft function for the double hemisphere jet mass distribution in e+e− collisions, a necessary ingredient for the calculation of several event shape distributions at N3LL order. The use of the mass as an infrared regulator allows us to derive the momentum space results for the massless quark structures at O(α2sCFTFnf) and the gluonic structures at O(α2sCACF), which have not been given so far in the literature. Furthermore, we compute the corresponding corrections in the soft function for thrust, the most prominent projection of the double hemisphere mass distribution. Finally we give expressions for the corresponding renormalon subtractions in the gap scheme.