In this work we present the computation of the missing pieces to get the bHQET thrust distribution with non-vanishing secondary quark masses at NNLO: the jet and hard functions. The difference with respect to the massless case is encoded in diagrams with massive-quark bubbles. For its computation we use a Mellin Barnes representation of the dispersive integral method that permits expr... https://www.roneverhart.com/Viva-La-Fleur-by-Juicy-Couture/