NLO Massive Event-Shape Differential and Cumulative Distributions

TL;DR

This paper introduces a universal, high-precision method for calculating event-shape distributions involving massive quarks at electron-positron colliders, with an efficient algorithm that avoids Monte Carlo methods.

Contribution

It provides a novel analytic and computational framework for NLO event-shape distributions with massive quarks, including universal coefficients and an efficient, non-Monte Carlo algorithm.

Findings

Derived universal coefficient for plus distribution in event shapes

Presented an analytic expression for the Dirac delta coefficient

Developed an algorithm achieving arbitrary precision in distribution calculations

Abstract

We provide a general method to effectively compute differential and cumulative event-shape distributions to $\mathcal{O}(\alpha_s)$ precision for massive quarks produced primarily at an $e^+e^-$ collider. In particular, we show that at this order, due to the screening of collinear singularities by the quark mass, for all event shapes linearly sensitive to soft dynamics, there appear only two distributions at threshold: a Dirac delta function and a plus distribution. Furthermore, we show that the coefficient of the latter is universal for any infra-red and collinear safe event shape, and provide an analytic expression for it. Likewise, we compute a general formula for the coefficient of the Dirac delta function, which depends only on the event-shape measurement function in the soft limit. Finally, we present an efficient algorithm to compute the differential and cumulative distributions, which does not rely on Monte Carlo methods, therefore achieving a priory arbitrary precision even in the extreme dijet region. We implement this algorithm in a numeric code and show that it agrees with analytic results on the distribution for 2-jettiness, heavy jet mass and a massive generalization of C-parameter.