REVIEW 4 major objections 4 minor 17 references
LHC and DIS experimental data in the CT18(Z) global QCD analysis
T0 review · 4 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read The CT18Z global QCD analysis shows that a tuned x-dependent factorization scale, $\mu_{F,x}^2=0.82(Q^2+0.3\,\mathrm{GeV}^2/x^{0.3})$, cuts the NNLO HERA inclusive-DIS $\chi^2$ by more than 50 units and reproduces the improvement…
desk verdict A solid proceedings writeup of the CT18 PDF family, but the resummation-mimicry claim rests on a scale tuned in-sample and needs an out-of-sample check. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the x-dependent factorization scale $\mu_{F,x}^2 = 0.82(Q^2+0.3\,\mathrm{GeV}^2/x^{0.3})$, a saturation-inspired one-scale ansatz whose coefficients were chosen to minimize the HERA DIS $\chi^2$. At small $x$, the added $0.3/x^{0.3}$ term raises the scale above $Q^2$, effectively absorbing a class of enhanced power logarithms that would otherwise be treated by low-$x$ resummation; at larger $x$ it reduces to a constant factor $0.82$ times $Q^2$. The same machinery also includes fast impact-analysis programs to select eleven new LHC data sets and a parallelized fitting code with fast grid interfaces, but the scale ansatz is what carries the central claim about matching resummation.
What would settle it
Take the fitted $\mu_{F,x}$ scale into a predictor outside the fitted HERA kinematics, for example forward charm or Drell-Yan production at the LHC at small $x$, or HERA data with $Q$ below 2 GeV, and compare CT18Z predictions with data and with a resummed PDF set. If the CT18Z predictions agree with data at resummation quality, the scale is a valid proxy; if the >50-unit improvement disappears or the predictions deviate once kinematics change, the ansatz is tuned rather than physical.
Extended reading notes
Core claim
On its own terms, the paper establishes that the conventional fixed-order NNLO description of HERA data can be substantially improved without resummation by evaluating DIS cross sections with the tuned scale $\mu_{F,x}^2=0.82(Q^2+0.3\,\mathrm{GeV}^2/x^{0.3})$ instead of $\mu_F^2=Q^2$. In the region $Q>2$ GeV, $x>10^{-5}$, this reduces $\chi^2(\mathrm{HERA\ I+II})$ by more than 50 units, a quality of improvement the paper describes as comparable to that reported by resummed analyses. The resulting CT18X and CT18Z fits show reduced $u$ and $d$ (anti-)quark PDFs and increased gluon and strangeness at $x<10^{-2}$, with compensating changes at $x>0.5$ to preserve valence and momentum sum rules. In CT18Z, removal of CDHSW nuclear data and inclusion of ATLAS 7 TeV W/Z data combine to reduce the NNLO gluon-fusion Higgs cross section by about 1% relative to CT14 and CT18.
Load-bearing premise
The premise is that a scale whose two coefficients were tuned to the HERA data can genuinely reproduce the small-$x$ logarithms that resummation computes explicitly, so its benefit should persist for unmeasured kinematics rather than being an in-sample artifact.
Editorial extensions
If this is right
- At NNLO, the small-$x$ behavior of the gluon and strangeness PDFs in CT18X/Z differs from CT18; quantities such as the gluon-fusion Higgs cross section shift by about 1%, altering the central value that precision electroweak and Higgs analyses compare with.
- The reported >50-unit $\chi^2$ drop means the combined HERA data do not, by themselves, force low-$x$ resummation: a fixed-order fit with the tuned scale is an equally good description, so resummed PDFs and scale-shifted PDFs should both be used for uncertainty bands.
- Because CT18, A, X, and Z are generated under different data-selection and scale assumptions, the spread among them gives a practical systematic uncertainty on NNLO predictions for LHC processes beyond the one-sigma Hessian error.
- Including the ATLAS 7 TeV W/Z rapidity data raises the tension with the strangeness-sensitive dimuon data from CCFR and NuTeV; the elevated goodness-of-fit variables in CT18Z show that the strangeness PDF is constrained by competing data sets, not uniquely by either.
Reading between the lines
- A natural extension not drawn by the paper: the same scale ansatz could be lifted into nuclear PDF fits or future electron-ion collider projections, where the small-$x$ regime is scanned directly; if it holds there, the saturation-inspired wording gains predictive content.
- An explicit overlap test would fit with both the tuned scale and a resummed theory simultaneously and see whether the scale loses its benefit; if it does, the two mechanisms are not independent alternatives but overlapping ways of absorbing the same logarithms.
- The tuned coefficients 0.82 and 0.3 depend on the data set; one could attempt to derive them from a first-principles saturation momentum and re-fit newer datasets to see whether the same functional form recurs, rather than being a one-off adjustment.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This proceedings contribution from the CTEQ-TEA collaboration reports on the implementation of HERA I+II and LHC data in the new CT18 global QCD analysis. The paper presents four NNLO (and corresponding NLO) PDF families: CT18, CT18A, CT18X, and CT18Z, which differ in data selection and in the factorization scale used for DIS. The main physics claim is that using an x-dependent factorization scale mu_F,x^2 = 0.82(Q^2 + 0.3 GeV^2/x^0.3) in fixed-order NNLO DIS calculations reduces the HERA I+II inclusive DIS chi-square by more than 50 units, an improvement described as comparable to that obtained with low-x resummation in Refs. [6,7]. The paper also discusses tensions among HERA, Tevatron, and LHC data sets, the use of PDFSense and ePump for data selection, and the impact of the alternative scales on PDFs and on the NNLO Higgs cross section.
Significance. If the central claim is correct, CT18Z and its siblings provide a systematic family of NNLO PDFs that span a range of plausible scale and data-selection choices, and the x-dependent scale idea could offer a practical fixed-order alternative to small-x resummation for global fits. The paper is valuable for its candid discussion of data tensions, its explicit statement that the x-dependent scale coefficients are tuned to minimize the HERA chi-square, and its use of modern fast-interfacing and reweighting tools. However, the key quantitative evidence for the resummation-mimicry claim is currently in-sample: the scale is fitted to the same HERA data whose chi-square improvement is reported, and no out-of-sample test or uncertainty analysis is provided. The significance of CT18X/Z therefore rests on validation that is not presented in this manuscript.
major comments (4)
- [Combined HERA I+II DIS data and an x-dependent factorization scale] The sentence defining mu_F,x^2 states that the numerical coefficients are chosen to minimize chi-square for the HERA DIS data. The subsequent claim that this scale produces a reduction of more than 50 chi-square units and a 'comparable quality of improvement' to low-x resummation is therefore an in-sample fit statistic, not independent evidence of physical equivalence. The manuscript should provide an out-of-sample validation: for example, fix the coefficients using a training subset of HERA data or a HERA-only fit, then evaluate the held-out HERA bins and LHC observables sensitive to the small-x gluon. The improved chi-square reported for charm SIDIS and H1 FL in Fig. 1(right) does not resolve this issue because these data sets are part of the same fit.
- [Figure 1] Figure 1(left) shows ratios of CT18 PDFs obtained with the x-dependent and standard scales, but no PDF uncertainty bands or chi-square decomposition are provided. Without error bands and per-data-set chi-square values, the reader cannot determine whether the enhanced small-x gluon and altered quark distributions are statistically required or are an artifact of the tuned scale. The manuscript should show the uncertainties of CT18 and CT18Z and report the chi-square per data set for both scale choices.
- [Combined HERA I+II DIS data and an x-dependent factorization scale] The comparison with low-x resummation in Refs. [6,7] is not apples-to-apples as presented: the data selections, kinematic cuts, chi-square definitions, and perturbative orders differ between the CT18 fit and the resummed fits. The claim that the improvement is 'comparable' requires a table specifying the data sets, Q and x cuts, and chi-square definitions used in each comparison. Without this, the >50 unit improvement cannot be interpreted as equivalent to the resummation effect.
- [Selection of new LHC experiments] The abstract and body present CT18, A, X, and Z as four new PDF families, but this proceedings contribution does not provide the defining criteria for each family beyond brief descriptions: no LHAPDF names, no complete input data lists, no central chi-square values, and no error sets are given. Since these definitions are load-bearing for the paper's central claim, the manuscript should either include a compact table summarizing the data sets and scale choices for each family or explicitly refer to a publicly available companion document where these details are given.
minor comments (4)
- [Combined HERA I+II DIS data and an x-dependent factorization scale] There are typographical errors in this section: 'thee+p and e−p' should be 'the e+p and e−p', 'with e HERA' should be 'with HERA', and 'multi-prone' should be 'multi-pronged'.
- [Figure 1] The left panel's caption contains a garbled formula, 'a (x,Q)/f (2)a f', which appears to be an attempt to write the ratio of PDFs; it should be cleaned up to something like f_a^{(1)}(x,Q)/f_a^{(2)}(x,Q).
- [Combined HERA I+II DIS data and an x-dependent factorization scale] The phrase 'unconstrained region x > 0.5' is imprecise: the valence quark region is constrained by sum rules, so it would be clearer to say 'a region less tightly constrained by the fitted data' or to specify what is meant by 'unconstrained'.
- [Figure 2] The caption labels 'CT18 NNLO' and 'CT18Z NNLO' but the text refers to 'CT14HERA2' in the upper panel; the notation and the relationship between the two panels should be clarified.
Circularity Check
The x-dependent scale is tuned to the same HERA chi-square that is then quoted as evidence of resummation-level improvement.
-
fitted input called prediction
[Section 'Combined HERA I+II DIS data and an x-dependent factorization scale' (pp. 1-2), inline definition of mu_F,x and first-figure discussion]
"In our analysis, we observe that, by evaluating the DIS cross sections at NNLO with an x-dependent factorization scale, such as a tuned scale μ2F,x = 0.82(Q2 + 0.3 GeV2/x0.3), instead of the conventional choice μ2F = Q2, we achieve a comparable quality of improvement ... Namely, the χ2(HERA I+II) reduces by > 50 units ... and the numerical coefficients in μ2F,x are chosen to minimize χ2 for the HERA DIS data."
The three coefficients of μ2F,x are explicitly fitted to minimize χ2 for the HERA I+II inclusive DIS data, and the same χ2 reduction (>50 units) is then reported as evidence that this scale choice matches the improvement attributed to low-x resummation. This is an in-sample fit improvement, not an independent prediction: three free parameters tuned to a large precise data set are expected to reduce its χ2, so the comparison to fixed-scale resummation is not a controlled out-of-sample test. The paper even labels the scale 'tuned' and states that its coefficients are chosen to minimize the very χ2 that is quoted.
full rationale
The CT18, CT18A, CT18X, and CT18Z PDFs are produced by a standard global QCD fit to many independent data sets (HERA, Tevatron, LHC jets, W/Z and top production, dimuon production, etc.), so the analysis as a whole is not circular and is benchmarked against external data. The self-citations to the companion CT18 paper, CT14HERA2, ePump, and PDFSense are methodological and not load-bearing in a circular way. The one genuine circular element is the x-dependent factorization scale: its numerical constants are chosen to minimize the HERA I+II DIS chi-square, and the resulting in-sample chi-square gain is then used to claim equivalence to small-x resummation. This is a fitted input being presented as evidence for a physics conclusion without an out-of-sample check or uncertainty quantification for the resulting PDF changes. Because this affects only the specific 'mimics resummation' claim and not the central PDF-determination machinery, the overall circularity score is moderate rather than severe.
Assumptions & free parameters
free parameters (4)
- 0.82 prefactor in x-dependent scale =
0.82
- 0.3 GeV^2 offset in x-dependent scale =
0.3 GeV^2
- 0.3 exponent of x in x-dependent scale =
0.3
- PDF parametrization coefficients (Bernstein polynomial)
assumptions (4)
- domain assumption NNLO collinear factorization is valid for all included DIS and LHC observables.
- ad hoc to paper The x-dependent scale ansatz can reproduce the effect of small-x resummation.
- domain assumption The HERA I+II combined data provide unbiased constraints on quark and gluon PDFs.
- domain assumption The Hessian error PDF methodology and tolerance criteria are valid.
Cite this review
Pith. "Pith review of LHC and DIS experimental data in the CT18(Z) global QCD analysis." pith.science (2026). https://pith.science/paper/DYZH2RXS
@misc{pith2026190900001,
author = {Pith},
title = {Pith review of: LHC and DIS experimental data in the CT18(Z) global QCD analysis},
year = {2026},
howpublished = {\url{https://pith.science/paper/DYZH2RXS}},
note = {Machine review of arXiv:1909.00001}
}
read the original abstract
We discuss implementation of the LHC experimental data sets in the new CT18 global analysis of quantum chromodynamics (QCD) at the next-to-next-leading order of the QCD coupling strength. New methodological developments in the fitting methodology are discussed. Behavior of the CT18 NNLO PDFs for the conventional and "saturation-inspired" factorization scales in deep-inelastic scattering is reviewed. Four new families of (N)NLO CTEQ-TEA PDFs are presented: CT18, A, X, and Z.
Figures
Reference graph
Works this paper leans on
-
[1]
T.-J. Hou, K. Xie, J. Gao, S. Dulat, M. Guzzi, T. Hobbs et al., New CTEQ Global Analysis with High Precision Data from the LHC, in Proceedings, DIS 2019 workshop, Torino, Italy, April 2019
work page 2019
-
[2]
S. Dulat, T.-J. Hou, J. Gao, M. Guzzi, J. Huston, P. Nadolsky et al., New parton distribution functions from a global analysis of quantum chromodynamics, Phys. Rev. D93 (2016) 033006 [1506.07443]. 4 Experimental data in the CT18(Z) analysis Pavel Nadolsky
arXiv 2016
-
[3]
H1, ZEUS collaboration, Combination of measurements of inclusive deep inelastic e± p scattering cross sections and QCD analysis of HERA data, Eur. Phys. J.C75 (2015) 580 [1506.06042]
arXiv 2015
-
[4]
H1, ZEUS collaboration, Combined Measurement and QCD Analysis of the Inclusive e+- p Scattering Cross Sections at HERA, JHEP 01 (2010) 109 [0911.0884]
arXiv 2010
-
[5]
T.-J. Hou, S. Dulat, J. Gao, M. Guzzi, J. Huston, P. Nadolsky et al., CTEQ-TEA parton distribution functions and HERA Run I and II combined data, Phys. Rev. D95 (2017) 034003 [1609.07968]
arXiv 2017
-
[6]
R. D. Ball, V . Bertone, M. Bonvini, S. Marzani, J. Rojo and L. Rottoli, Parton distributions with small-x resummation: evidence for BFKL dynamics in HERA data, Eur. Phys. J.C78 (2018) 321 [1710.05935]
arXiv 2018
-
[7]
XFITTER DEVELOPERS ’ TEAM collaboration, Impact of low-x resummation on QCD analysis of HERA data, Eur. Phys. J.C78 (2018) 621 [1802.00064]
arXiv 2018
- [8]
Show all 17 references
-
[9]
H1 collaboration, Measurement of the Inclusive e± p Scattering Cross Section at High Inelasticity y and of the Structure Function FL, Eur. Phys. J.C71 (2011) 1579 [1012.4355]
2011 arXiv
-
[10]
Kovarik, P
K. Kovarik, P. M. Nadolsky and D. E. Soper, Hadron structure in high-energy collisions, 1905.06957
1905 arXiv
-
[11]
H.-L. Lai, M. Guzzi, J. Huston, Z. Li, P. M. Nadolsky, J. Pumplin et al., New parton distributions for collider physics, Phys. Rev. D82 (2010) 074024 [1007.2241]
2010 arXiv
-
[12]
ATLAS collaboration, Precision measurement and interpretation of inclusive W+ , W− and Z/γ∗ production cross sections with the ATLAS detector, Eur. Phys. J.C77 (2017) 367 [1612.03016]
2017 arXiv
-
[13]
Rojo et al., The PDF4LHC report on PDFs and LHC data: Results from Run I and preparation for Run II, J
J. Rojo et al., The PDF4LHC report on PDFs and LHC data: Results from Run I and preparation for Run II, J. Phys. G42 (2015) 103103 [1507.00556]
2015 arXiv
-
[14]
B.-T. Wang, T. J. Hobbs, S. Doyle, J. Gao, T.-J. Hou, P. M. Nadolsky et al.,Mapping the sensitivity of hadronic experiments to nucleon structure, Phys. Rev. D98 (2018) 094030 [1803.02777]
2018 arXiv
-
[15]
Schmidt, J
C. Schmidt, J. Pumplin and C.-P. Yuan, Updating and optimizing error parton distribution function sets in the Hessian approach, Phys. Rev. D98 (2018) 094005 [1806.07950]
2018 arXiv
-
[16]
Kluge, K
T. Kluge, K. Rabbertz and M. Wobisch, FastNLO: Fast pQCD calculations for PDF fits, in Proceedings, DIS 2006, Tsukuba, Japan, April 2006, p. 483, 2006, hep-ph/0609285
2006 arXiv
-
[17]
Carli, D
T. Carli, D. Clements, A. Cooper-Sarkar, C. Gwenlan, G. P. Salam, F. Siegert et al., A posteriori inclusion of parton density functions in NLO QCD final-state calculations at hadron colliders: The APPLGRID Project, Eur. Phys. J.C66 (2010) 503 [0911.2985]. 5
2010 arXiv
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.