Determining $\alpha_s(m_Z)$ from the Heavy Jet Mass distribution

Miguel A. Benitez

REVIEW 3 major objections 5 minor 14 references

Determining $\alpha_s(m_Z)$ from the Heavy Jet Mass distribution

T0 review · 3 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read The paper argues that heavy jet mass data, when resummed in both the dijet and three-jet shoulder regions, give $\alpha_s(m_Z)=0.1145$, matching thrust and C-parameter, and a negative trijet hadronization shift.

desk verdict A plausible alpha_s extraction from heavy jet mass that resolves the tension with thrust and C-parameter, but the negative trijet power correction is not yet robust—it flips sign with shoulder resummation and no stability test is shown. read the letter →

arxiv 2506.07723 v1 pith:RKJHFKRX submitted 2025-06-09 hep-ph

Miguel A. Benitez This is my paper

classification hep-ph

keywords alpha_sstrongcouplingheavyjetmasseventshapesresummationshoulderlogarithmspowercorrectionshadronizatione+e-annihilation

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper argues that the historically low values of $\alpha_s(m_Z)$ extracted from the heavy jet mass distribution disappear once the theory includes both dijet resummation and resummation of the shoulder logarithms around the symmetric three-jet point, together with a theory covariance matrix in the fit. Using 700 $e^+e^-$ data points from 35 to 207 GeV, the best fit gives $\alpha_s(m_Z)=0.1145$ with an uncertainty near 0.002, compatible with thrust and C-parameter results. The same fit yields a negative trijet hadronization shift $\Theta_1=-0.50\pm0.17$ GeV, but only when shoulder resummation is included; without it both nonperturbative parameters stay positive. If correct, the result removes a long-standing tension between heavy jet mass and more inclusive event shapes and provides evidence for a negative three-jet power correction.

What carries the argument

The argument rests on two factorization formulas and a matching prescription. In the dijet limit, Eq. (1) factors the cross section into hard, jet, and soft functions convoluted with a shape function $F_{1,2}^\Xi(\Omega_1^\rho)$ encoding dijet hadronization. Near the symmetric three-jet point $\rho\to 1/3$ from the left, Eq. (2) factors the perturbative cross section as $H^{sh}J_1J_2J_3\otimes S_{1,2,3}$, and nonperturbative trijet effects enter as a simple shift $\Theta_1/Q$ in Eq. (3). The matching in Eq. (4) combines dijet, fixed-order, and shoulder pieces with overlap subtractions, and profile functions interpolate the scales across regions. A theory covariance matrix built from 5000 random profile-parameter sets is added to the experimental covariance in the $\chi^2$, which is what lets the fit claim a reliable uncertainty on $\alpha_s$ and on the sign of $\Theta_1$.

What would settle it

Refit the 700 data points with the same covariance and fit parameters but exclude the shoulder contribution entirely and force $\Theta_1=0$; if the $\chi^2$ degrades by less than one unit per degree of freedom, the data do not actually demand the negative trijet shift. A sharper test is to measure the heavy jet mass distribution near $\rho=1/3$ with sub-percent precision at several center-of-mass energies and check whether the shift scales as $\Theta_1/Q$ as Eq. (3) requires.

Watch

Extended reading notes

Core claim

The central claim is that a complete resummed description of the heavy jet mass spectrum, combining a factorized dijet region, a factorized shoulder region, and fixed order through profile functions, fits the available $e^+e^-$ data with $\alpha_s(m_Z)=0.1145\pm0.0020$ (theory plus experiment) and a small fit-range uncertainty. The accompanying nonperturbative parameters are $\Omega_1^\rho=0.57\pm0.09$ GeV and $\Theta_1=-0.50\pm0.17$ GeV. The negative $\Theta_1$ appears only when shoulder resummation is included, and the fitted ratio $\Omega_1^\rho/\Theta_1$ agrees with the prediction for the three-jet power-correction function $\zeta(\rho)$ derived in Refs. [6,7]. The paper therefore claims both that heavy jet mass no longer disagrees with thrust and C-parameter determinations of the strong coupling and that the data prefer a negative hadronization shift in the three-jet shoulder region.

Load-bearing premise

The result rests on the assumption that the three-jet shoulder region is described by a factorized formula and that hadronization there is a single shift of the distribution; if either assumption fails, the fitted trijet shift and its negative sign would be biased, and the coupling would shift less but not be immune.

Editorial extensions

If this is right

Heavy jet mass becomes a standard observable for $\alpha_s(m_Z)$ determinations, with a central value and uncertainty comparable to thrust and C-parameter fits.
The historical gap between heavy jet mass and inclusive event shapes is explained by missing shoulder resummation and fixed-order-only treatment rather than by new physics or data problems.
The fitted $\Theta_1<0$ supports the existence of negative three-jet power corrections in $e^+e^-$ event shapes, quantified by the ratio $\Omega_1^\rho/\Theta_1$.
Including the theory covariance matrix stabilizes the fit against the choice of fit range, so the quoted uncertainty can be interpreted as a genuine total uncertainty.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

A direct test would refit the same data with an alternative resummation scheme that does not assume the shoulder factorization of Eq. (2); if the negative $\Theta_1$ does not survive, the sign is a model-dependent artifact rather than a property of QCD.
Because $\Theta_1$ is defined as a universal trijet hadronization parameter, its negative value could be cross-checked in other observables sensitive to three-jet configurations, such as the C-parameter in the symmetric limit or thrust in the three-jet region, once analogous shoulder resummations are implemented.
The quoted central value suggests that future high-precision data at higher energies, where power corrections are suppressed, should shift $\alpha_s$ only within the quoted 0.002 band; a larger shift would indicate missing higher-order or nonperturbative effects.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Quick take: this is a conference proceedings that distills the fit machinery from the companion paper (arXiv:2502.12253) into a compact alpha_s(mZ) extraction from heavy jet mass data. The alpha_s result, 0.1145 with total uncertainty near 0.002, is plausible and consistent with thrust and C-parameter. That part holds up. What does not yet hold up is the paper's more striking claim, the negative trijet hadronization parameter Theta1 = -0.50 ± 0.17 GeV. The sign is entirely a consequence of including shoulder resummation: without it, Theta1 is +0.53 ± 0.13. That ~1 GeV shift is six times the quoted uncertainty, and the paper gives no stability test, e.g., varying the shoulder profile parameters or matching scale, to show the negative sign survives reasonable variations. I trust the stress-test note here: Theta1 is a fitted non-perturbative parameter and can absorb missing higher-order shoulder logarithms. The paper itself says data favor a negative sign "only if shoulder resummation is included," which is honest, but it does not establish that the sign is a robust physical feature rather than a degeneracy with the newly resummed logs.

The paper does several things well. The fit setup is described cleanly: a 5000-profile random scan, theory covariance matrix via Eqs. (6)-(7), and a 700-point dataset. The demonstration that the fixed-order result is linearly sensitive to the fit-range lower bound, and that dijet resummation removes that sensitivity, is a nice pedagogical point. The comparison of the fitted Omega1/Theta1 ratio with the Refs. [6,7] prediction for the 3-jet power correction is also a reasonable cross-check, though it is a check against a fitted ratio rather than an a priori prediction, so the reader should not over-weight it.

The main weakness is that the paper is not self-contained. The theoretical prediction is borrowed from Ref. [11], the fit-range prescription is in Ref. [11], and no code or data listing is provided. For a proceedings this is understandable, but it means the numerical results are not independently verified here. The 17 profile parameters are mentioned but not specified, so a reader cannot reproduce the covariance matrix without the companion paper. I would not count this as a fatal flaw, but it is a real limitation.

Bottom line: the alpha_s extraction is likely solid and worth taking seriously; the negative Theta1 is an interesting hint, not a demonstrated result. The paper is for QCD phenomenologists working on event shapes and power corrections. I'd rather referee the companion paper than this proceedings, but the proceedings is a fair summary and the underlying work deserves refereeing. I would not cite this version; I'd cite Ref. [11].

Referee Report

3 major / 5 minor

Summary. This paper presents a fit of the strong coupling constant and two nonperturbative parameters to e+e- heavy jet mass event-shape data from 35 to 207 GeV (700 points). The theoretical model combines fixed-order QCD, dijet resummation, and resummation of shoulder logarithms near rho=1/3, with nonperturbative corrections encoded by a dijet shape function (Omega_1^rho) and a trijet shift (Theta_1). Experimental correlations are modeled with the minimal overlap model (Eq. 5), and a theory covariance matrix is built from a 5000-point scan over 17 profile-function parameters (Eqs. 6-7). The main result in Table 1 is alpha_s(m_Z)=0.1145 with a combined uncertainty near 0.002 in the preferred FO+dijet+shoulder 3D fit, together with Omega_1^rho=0.57 +/- 0.09 GeV and Theta_1=-0.50 +/- 0.17 GeV. The paper explicitly states that the negative sign of Theta_1 appears only when shoulder resummation is included, and it interprets this as evidence for a negative trijet power correction consistent with Refs. [6,7].

Significance. The extraction is valuable if the theoretical framework is sound: it addresses the long-standing low alpha_s puzzle in HJM by including shoulder resummation, and it makes a concrete, falsifiable claim about a negative trijet hadronization parameter. The use of a profile-scan theory covariance matrix and the explicit comparison across fixed-order, dijet-only, and dijet+shoulder setups are strengths, and the paper is honest that the negative Theta_1 sign is conditional on shoulder resummation. If confirmed, the alpha_s value is compatible with thrust and C-parameter determinations, providing another independent input to the world average. However, the significance is tempered by the manuscript's heavy reliance on Ref. [11] for the theoretical prediction and fit-range prescription, which makes the present paper incomplete as a standalone determination. In addition, the central novel claim, the negative Theta_1, is the least stable result in Table 1, and its robustness is not yet demonstrated.

major comments (3)

[Section 4, Table 1] The central claim of a negative 3-jet power correction rests entirely on the difference between the FO+dijet 3D and FO+dijet+shoulder 3D rows of Table 1: Theta_1 changes from +0.53 +/- 0.13 GeV to -0.50 +/- 0.17 GeV, a shift of about 1.0 GeV that is roughly six times the quoted uncertainty. The paper states that data favor a negative sign only when shoulder resummation is included, but it does not provide any stability test, such as a scan over the shoulder matching scale or the shoulder profile parameters, showing that Theta_1 remains negative under reasonable variations of the shoulder treatment. Because Theta_1 is a nonperturbative parameter, it can absorb missing higher-order shoulder logarithms, and the large model-to-model shift suggests a degeneracy between the shoulder resummation prescription and the fitted Theta_1. Without such a test, the evidence for a negative Theta_1 is not robust, even though the alpha_s result is more stable.
[Section 3, Eqs. (6)-(7)] The theory covariance matrix is constructed from a 5000-point flat random scan over '17 theory parameters', but the manuscript does not specify what these parameters are, their allowed ranges, or the profile-function definitions. The central values xbar_i and uncertainties Delta_i are taken from the min/max of the same scan, and the correlation coefficients r_theo_ij are computed from the same ensemble; this is a particular estimator whose behavior for non-Gaussian, bounded profile scans is not discussed. Since the total covariance matrix enters the chi^2 in Eq. (8) and directly affects all fit results and uncertainties, the fit is not reproducible from the information given. The delegation of the prescription to Ref. [11] may be acceptable for a proceedings contribution, but for the quantitative claims made here, either the details of the scan or a comparison with an independent covariance construction should be provided.
[Section 2, Eqs. (1)-(3)] The negative Theta_1 extraction assumes that nonperturbative effects in the shoulder region are described by a single shift (Eq. 3). The paper does not assess the size of neglected higher-order power corrections or alternative modelings of trijet hadronization, so a negative fitted shift could absorb a deficiency of the shoulder resummation or of the profile matching. In addition, the agreement with the prediction of Refs. [6,7] is stated qualitatively: the ratio Omega_1^rho/Theta_1 is said to agree, but no numerical value or uncertainty for the ratio is reported. The paper should present the actual ratio with its uncertainty and, if possible, a direct comparison with the zeta(rho) prediction, rather than relying on a fitted quantity.

minor comments (5)

[Throughout] The manuscript contains many typographical and formatting artifacts, including missing spaces and broken equation renderings (for example, in Section 2 the factorization formula appears as 'dσdij =H dij ×J 1 ×J 2 ⊗S 1,2 ⊗F Ξ 1,2(Ωρ 1)' and the Introduction contains 'determinations1,2,3'). A careful copyedit is needed.
[Figure 1] The caption of Figure 1 does not identify the color coding in the right panel for Omega_1^rho and Theta_1; the text refers to green, gray, and blue curves, but without the figure the reader cannot map the colors to the two parameters.
[Table 1] The header of Table 1 is confusing: the three uncertainty entries in the alpha_s columns are not individually labeled, and the fit-range column appears merged with the chi^2/dof column. The table should be reformatted so that each uncertainty component (statistical, experimental systematic, theory, fit-range) is explicit.
[Section 3, Eq. (5)] The 'minimal overlap model' is introduced without a reference or derivation; a citation to its original definition would help the reader evaluate the treatment of experimental systematic correlations.
[Section 4, Table 1] The fit-range weighted-average prescription is mentioned only by reference to Ref. [11]; even a one-sentence summary of the averaging formula would make the central values and fit-range uncertainties in Table 1 self-contained.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the strong coupling and power-correction parameters are obtained by fitting external data; the self-cited companion paper is ancillary.

full rationale

The paper's central results are determined by minimizing the chi-squared in Eq. (8) against 700 external e+e- data points, with experimental and theory covariance matrices treated as inputs. The theoretical framework, including dijet factorization Eq. (1), shoulder factorization Eq. (2), the non-perturbative shift Eq. (3), and the matching in Eq. (4), is imported from prior published work rather than derived from the quantities being fitted. The parameters alpha_s, Omega_1^rho, and Theta_1 are extracted, not predicted, so no output is recycled as an input. The comparison with Refs. [6,7] is a post-fit check of the fitted ratio Omega_1^rho/Theta_1 against an independent external calculation, not a claim that the ratio is predicted by the present analysis. The statement that the negative Theta_1 sign appears only when shoulder resummation is included is a stability caveat about missing higher-order shoulder logarithms, not a circular step: the sign is a data-driven fit outcome. The single self-citation, Ref. [11], is used only for the prescription for averaging over fit-range variations and does not carry the central derivation. No equation reduces to itself and no fitted parameter is relabeled as a prediction.

Assumptions & free parameters 4 free parameters · 5 assumptions · 0 invented entities

The paper introduces no new fundamental entities. It fits three physical parameters (alpha_s, Omega_1^rho, Theta_1) and uses 17 profile parameters to generate a theory covariance. The main assumptions are the factorization and the non-perturbative shift model.

free parameters (4)

alpha_s(mZ) = 0.1145 (FO+dijet+shoulder 3D, Table 1)
The strong coupling constant is the main fit parameter, determined by minimizing Eq. (8) against e+e- data.
Omega_1^rho = 0.57 ± 0.09 GeV (best fit)
Non-perturbative shape-function parameter for dijet hadronization, fitted simultaneously with alpha_s.
Theta_1 = -0.50 ± 0.17 GeV (best fit)
Non-perturbative shift parameter for trijet hadronization in Eq. (3), fitted simultaneously; its sign is the paper's key finding.
17 profile-function parameters = not fitted; scanned over 5000 random sets
Profile-function parameters controlling the scale evolution in the dijet and shoulder regions. Their scan defines the theory covariance matrix (Eq. 7); the prior ranges are not specified in this paper.

assumptions (5)

domain assumption SCET factorization in the dijet and shoulder regions (Eqs. 1 and 2) holds with the stated hard, jet and soft functions.
The entire theoretical prediction rests on this factorization; deviations would propagate directly into the fitted alpha_s and Theta_1.
ad hoc to paper Non-perturbative corrections are described by a single shape function for dijet events and a single shift Theta_1 for trijet events (Eq. 3).
The shift model is a phenomenological assumption; the paper's 'negative power correction' claim depends on it being a valid description.
ad hoc to paper The minimal overlap model (Eq. 5) correctly captures experimental systematic correlations.
The fitting procedure relies on this model for the experimental covariance matrix.
ad hoc to paper The covariance matrix built from 5000 random variations of the 17 profile parameters (Eq. 7) represents the theoretical uncertainty from missing higher orders.
This is a heuristic uncertainty estimate; its faithfulness is not justified in the paper.
ad hoc to paper The weighted average over variations of the fit-range lower bound yields an unbiased central value (prescription from Ref. [11]).
The central values in Table 1 depend on this averaging prescription, which is not described in this paper.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Determining $\alpha_s(m_Z)$ from the Heavy Jet Mass distribution." pith.science (2026). https://pith.science/paper/RKJHFKRX

@misc{pith2026250607723,
  author       = {Pith},
  title        = {Pith review of: Determining $\alpha_s(m_Z)$ from the Heavy Jet Mass distribution},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/RKJHFKRX}},
  note         = {Machine review of arXiv:2506.07723}
}

abstract

We present a state-of-the-art analysis on determining the strong coupling constant $\alpha_s(m_Z)$ from available $e^+e^-$ data for the Heavy Jet Mass (HJM) distribution. Dijet resummation is supplemented with additional resummation of shoulder logarithms appearing around the symmetric 3-jet configuration of the HJM spectrum. In addition to $\alpha_s(m_Z)$ we further obtain numerical results for two non-perturbative parameters, $\Omega_1^{\rho}$ and $\Theta_1$, encompassing effects from dijet and trijet hadronization effects, respectively.

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

14 extracted references · 12 canonical work pages

[11]

M. A. Benitez, A. Bhattacharya, A. H. Hoang, V. Mateu, M. D. Schwartz, I. W. Stewart and X. Zhang, [arXiv:2502.12253 [hep-ph]]

arXiv
[1]

G. P. Salam and D. Wicke, JHEP05(2001), 061

work page 2001
[2]

Dissertori, A

G. Dissertori, A. Gehrmann-De Ridder, T. Gehrmann, E. W. N. Glover, G. Heinrich, G. Luisoni and H. Stenzel, JHEP08(2009), 036

work page 2009
[3]

Y. T. Chien and M. D. Schwartz, JHEP08(2010), 058

work page 2010
[4]

Bhattacharya, M

A. Bhattacharya, M. D. Schwartz and X. Zhang, Phys. Rev. D106(2022) no.7, 074011

work page 2022
[5]

Bhattacharya, J

A. Bhattacharya, J. K. L. Michel, M. D. Schwartz, I. W. Stewart and X. Zhang, JHEP11 (2023), 080

work page 2023
[6]

Caola, S

F. Caola, S. Ferrario Ravasio, G. Limatola, K. Melnikov and P. Nason, JHEP01(2022), 093

work page 2022
[7]

Caola, S

F. Caola, S. Ferrario Ravasio, G. Limatola, K. Melnikov, P. Nason and M. A. Ozcelik, JHEP12(2022), 062

work page 2022

Show all 14 references

[8]

Catani, G

S. Catani, G. Turnock and B. R. Webber, Phys. Lett. B272(1991), 368-372

1991
[9]

G. P. Korchemsky and G. F. Sterman, Nucl. Phys. B555(1999), 335-351

1999
[10]

M. D. Schwartz, Phys. Rev. D77(2008), 014026

2008
[12]

Abbate, M

R. Abbate, M. Fickinger, A. H. Hoang, V. Mateu and I. W. Stewart, Phys. Rev. D83 (2011), 074021

2011
[13]

M. A. Benitez, A. H. Hoang, V. Mateu, I. W. Stewart and G. Vita, [arXiv:2412.15164 [hep-ph]]

arXiv
[14]

A. H. Hoang, D. W. Kolodrubetz, V. Mateu and I. W. Stewart, Phys. Rev. D91(2015) no.9, 094018

2015

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.

[1] [11]

M. A. Benitez, A. Bhattacharya, A. H. Hoang, V. Mateu, M. D. Schwartz, I. W. Stewart and X. Zhang, [arXiv:2502.12253 [hep-ph]]

arXiv

[2] [1]

G. P. Salam and D. Wicke, JHEP05(2001), 061

work page 2001

[3] [2]

Dissertori, A

G. Dissertori, A. Gehrmann-De Ridder, T. Gehrmann, E. W. N. Glover, G. Heinrich, G. Luisoni and H. Stenzel, JHEP08(2009), 036

work page 2009

[4] [3]

Y. T. Chien and M. D. Schwartz, JHEP08(2010), 058

work page 2010

[5] [4]

Bhattacharya, M

A. Bhattacharya, M. D. Schwartz and X. Zhang, Phys. Rev. D106(2022) no.7, 074011

work page 2022

[6] [5]

Bhattacharya, J

A. Bhattacharya, J. K. L. Michel, M. D. Schwartz, I. W. Stewart and X. Zhang, JHEP11 (2023), 080

work page 2023

[7] [6]

Caola, S

F. Caola, S. Ferrario Ravasio, G. Limatola, K. Melnikov and P. Nason, JHEP01(2022), 093

work page 2022

[8] [7]

Caola, S

F. Caola, S. Ferrario Ravasio, G. Limatola, K. Melnikov, P. Nason and M. A. Ozcelik, JHEP12(2022), 062

work page 2022