{"id":"424cd354-d24f-42fd-bbfc-8f5064ae3b51","arxiv_id":"2412.03724","paper_version":2,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":6,"one_line_summary":"A simple two-parameter jet energy-loss model can describe centrality-dependent jet suppression measurements, but the inferred formation time depends on the assumed path-length scaling and on how systematic errors are correlated.","lead":"This paper fits a simple two-parameter formula for how jets lose energy in hot nuclear matter to measurements of jet suppression in lead-lead collisions, and shows that the answers depend strongly on how experimental errors are assumed to be correlated. The work provides a fast, public model that the community can use to test which jet observables and error treatments matter most.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Shape-insensitivity claim is tested only for central R_AA; centrality-dependent fits use a delta distribution, leaving untested whether gamma-distribution shape alters the inferred tau_f preference.","rationale":"The paper is a clear, honest sensitivity study and the public code is a plus. The two-parameter mean energy-loss model does describe central R_AA, and the tau_f sensitivity is demonstrated with chi2 differences that are large in the fully-covariant case. My concern is narrower but load-bearing: the claim that R_AA data alone constrain only the mean energy loss is used to justify using the delta distribution for the centrality-dependent analysis, yet that claim is only tested on the central (0-10%) bin. The centrality-dependent fits never include the gamma shape parameter. Since centrality scans over a wide range of path lengths and medium densities, it could in principle provide new sensitivity to the shape of the energy-loss distribution, and if so the delta-only centrality fits could bias the inferred tau_f preference. This is an untested assumption rather than an internal inconsistency, and it can be settled by a concrete rerun of the Table IV fits with the gamma distribution. Given that the abstract and conclusions state the shape-insensitivity and tau_f results without this caveat, I recommend conditional acceptance: the authors should either perform the gamma-centrality test or explicitly restrict the shape-insensitivity conclusion to the central data.","tokens_in":13468,"tokens_out":11558,"duration_ms":121412,"concrete_test":"Rerun the Table IV Bayesian calibration using the gamma-distribution energy loss (Eq. 3) with k as a free parameter (prior [1,10]) for both linear and quadratic path weightings and for tau_f = 0.1, 0.5, 0.9 fm/c, using the same semi-diagonal and fully-covariant error matrices. If the posterior for k is broad and the tau_f ordering in chi2 is unchanged, the shape-insensitivity claim is validated; if k is constrained or the preferred tau_f flips, the central claim must be narrowed to 'central R_AA data alone' and the delta-model tau_f values treated as shape-conditioned.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper concludes that 'jet-R_AA measurements constrain only the mean value of ΔpT' (Sec. V), but the supporting evidence in Sec. III is limited to the 0-10% central bin: there, gamma- and delta-distribution energy losses fit equally well and the gamma shape parameter k is poorly constrained. The centrality-dependent extension in Sec. IV then adopts only the delta-distribution (Eq. 4) and never tests whether the full centrality-dependent R_AA data, which spans much wider path-length and medium-density ranges, can constrain k. Because the reported tau_f preference (Table IV) is obtained with this delta-only assumption, it is conditional on an untested shape hypothesis. If a gamma distribution with moderate k altered the centrality dependence of R_AA(pT), the optimal tau_f for linear versus quadratic weighting could shift, weakening the paper's central sensitivity claim. The authors honestly flag the model's simplicity, but they do not flag this particular extrapolation from central shape-insensitivity to the centrality-dependent analysis.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a two-parameter empirical model for the mean transverse-momentum loss of jets in heavy-ion collisions, with the loss distribution taken either as a gamma distribution or a delta function. Bayesian fits to the ATLAS 0-10% central Pb+Pb jet R_AA at 5.02 TeV show that the two functional forms give equivalent agreement and that the gamma shape parameter is poorly constrained, leading the authors to conclude that these central data constrain only the mean energy loss. The model is then extended to centrality-dependent R_AA by embedding a 2D TRENTo initial geometry with either linear or quadratic path-length weighting and with the formation time tau_f treated as a fixed input in a sensitivity scan. Fits to the centrality-dependent ATLAS data for several tau_f values show that longer formation times are preferred for linear weighting while shorter formation times are preferred for quadratic weighting. The paper also constructs covariance matrices for the ATLAS systematic errors under different correlation assumptions, and shows that the fitted functions and chi2_d values are sensitive to these assumptions. Appendices examine alternative covariance constructions and investigate the role of Peele's Pertinent Puzzle in the low fitted curves.","tokens_in":13649,"tokens_out":6314,"duration_ms":62443,"significance":"If the claims hold, the paper provides a fast, public, and transparent platform for sensitivity studies of jet quenching, with code available at github.com/llnl/ez-quench. The systematic comparison of error-correlation treatments is a useful methodological contribution, and the demonstration that Peele's Pertinent Puzzle is not the sole cause of the low fits is a valuable negative result. The tau_f sensitivity to the path-length weighting is a concrete, falsifiable prediction that can be tested in more complete frameworks. The main limitation is that the centrality-dependent conclusion rests on an untested shape assumption and on fixed initial-state parameters; addressing these points would considerably strengthen the paper.","major_comments":[{"comment":"The conclusion that \"R_AA measurements constrain only the mean value of ΔpT\" is supported only by the fits to the 0-10% central bin in Sec. III, where the gamma- and delta-distribution fits give equivalent agreement and the gamma shape parameter k is poorly constrained. The centrality-dependent analysis in Sec. IV and Table IV uses only the delta-distribution (Eq. 4) and never tests whether the full centrality-dependent R_AA data, which span a much wider range of path lengths and medium densities, can constrain the gamma shape parameter. Because the reported tau_f preference is obtained under this delta-only assumption, it is conditional on an untested shape hypothesis; a gamma distribution with moderate k could in principle alter the centrality dependence of R_AA(pT) and shift the optimal tau_f for linear versus quadratic weighting. Please either add a gamma-distribution sensitivity scan to the centrality-dependent fits, or soften the conclusion to state explicitly that the mean-only constraint is demonstrated for the central bin and assumed, not tested, in the centrality-dependent extension.","section":"Sec. V (and Secs. III-IV)"},{"comment":"The TRENTo initial-state parameters are fixed to values from a prior Bayesian fit [12] (Table III) without propagating their uncertainties. Since the tau_f inference is driven by the path-length and density weighting of the medium, the sensitivity of the tau_f preference to plausible variations in these geometry parameters (for example, the nucleon width w, the fluctuation parameter k, or the reduced-thickness parameter p) should be quantified, or at least the limitation should be stated explicitly, so that the reader can judge the robustness of the main sensitivity claim.","section":"Sec. IV and Table III"}],"minor_comments":[{"comment":"In the last paragraph of Sec. I, \"Centrality dependent extensions and comparisons are provided in Sec. III\" should refer to Sec. IV, not Sec. III.","section":"Introduction"},{"comment":"The delta-distribution in Eq. (4) is written as f(ΔpT) = δ(pT − µ(pT)), which is inconsistent with the argument on the left-hand side; it should presumably read δ(ΔpT − µ(pT)) or use a different symbolic convention.","section":"Eq. (4)"},{"comment":"No chi-squared values are reported for the central gamma- and delta-distribution fits shown in Fig. 2, so the claim of \"equivalent agreement\" is not quantified; adding the chi2_d values for the four fits would strengthen the comparison.","section":"Sec. III"},{"comment":"There is a typo in the sentence \"we conclude that the jet-RAA measurements are sensitive only the the mean-value of the pT-loss\" — \"only the the\" should be \"only to the\".","section":"Sec. III, text near Fig. 5"},{"comment":"The word \"scale-paramter\" in the description of the gamma distribution should be \"scale-parameter\".","section":"Sec. II"},{"comment":"The table caption does not indicate that the fully-covariant fits are performed on log(R_AA) after the logarithmic transformation described in Sec. IV; adding this note would prevent confusion when comparing alpha and beta values to the semi-diagonal fits.","section":"Table IV and Sec. IV"},{"comment":"The legend label \"ALICE 0-10\" is missing the percent sign; it should be \"ALICE 0-10%\".","section":"Fig. 10 caption"},{"comment":"The sentence \"the higher pT of 800 GeV does not permit an extrapolation of the quenched cross-section above the highest ATLAS pT bin\" is unclear; rephrasing would help, for example by stating that the pp parameterization is only fitted up to 800 GeV and the model is not extrapolated above that value.","section":"Sec. IV, final paragraph"}],"recommendation":"major_revision","confidential_remarks":"The paper is a useful sensitivity study within the scope of a heavy-ion phenomenology journal, and the public code is a clear asset. The main gap is the untested gamma-shape assumption in the centrality-dependent analysis, which directly affects the robustness of the headline tau_f sensitivity claim. The issue is fixable without new data — a gamma-shape sensitivity scan or a carefully restricted conclusion would suffice — and the rest of the paper is sound. I would not reject the manuscript, but I would require the authors to address this point before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Punchline: this is a genuinely useful sensitivity study, not a new physics claim. It ships a public, fast two-parameter benchmark (ez-quench) and shows how error-correlation choices and formation-time assumptions move centrality-dependent R_AA fits.\n\nThe paper's real contribution is the systematic comparison of covariance-matrix constructions. Appendix A's distinction between summing error contributions before or after forming the covariance is important; the Peele's Pertinent Puzzle discussion in Appendix B is the clearest treatment I've seen in this context, and their reflected-data test is a nice way to separate genuinely correlated-error bias from shape mismatch. They show the fully-correlated assumption pulls fits below the data, the log-transform doesn't fix it, and the real culprit is a subtle pT-shape mismatch between the simple model and ATLAS. That's a practical lesson for anyone doing Bayesian calibration to jet R_AA.\n\nThe central shape-insensitivity claim (gamma vs delta fit equally well for 0-10%) is supported. The centrality-dependent analysis, though, goes straight to delta-only. So the tau_f preference in Table IV is conditional on the shape of the loss distribution not mattering at other centralities. That is an untested extrapolation, and the stress-test note is right that a gamma with moderate k could shift the centrality dependence of R_AA(pT) and thus the optimal tau_f. The authors flag the model's simplicity honestly, but they don't flag this particular step.\n\nOther soft spots are minor-to-moderate: TRENTo parameters fixed from a prior JETSCAPE fit without propagating uncertainties; tau_f scanned, not fitted; fully-correlated systematic errors are recognized by the authors as an approximation. None of this is a red flag; it's a consistent, honest sensitivity study.\n\nThe math and code look fine, the citations are appropriate, and the paper never overclaims. I'd send it to a serious referee. If I were working in that space, I'd cite it for the error-matrix comparison alone.","headline":"A useful, honest sensitivity study and fast public benchmark for the jet-quenching calibration community; the formation-time preference is real but conditional on an untested delta-only shape assumption.","tokens_in":14180,"tokens_out":1988,"would_cite":true,"duration_ms":20716,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The centrality-dependent jet suppression in lead-lead collisions at 5.02 TeV is described by a two-parameter mean energy-loss model, and the data alone constrain only the mean loss, not the shape of the loss distribution.","keywords":["jet quenching","nuclear modification factor","Bayesian inference","heavy-ion collisions","quark-gluon plasma","correlated systematic errors","formation time","energy loss"],"falsifier":"Vary the TRENTo initial-state parameters within their posterior uncertainties and refit the centrality-dependent $R_{AA}$: if the optimal $\\tau_f$ no longer shifts from 0.9 fm/c under linear path weighting to 0.1 fm/c under quadratic weighting, the formation-time conclusion fails. Alternatively, a precise dijet-asymmetry measurement that cannot be described by a delta-distribution loss with the same mean would falsify the claim that $R_{AA}$ alone constrains only the mean energy loss.","tokens_in":13251,"feed_emoji":"💥","tokens_out":17062,"duration_ms":142736,"temperature":0.7,"pith_summary":"The paper asks what measured jet suppression in heavy-ion collisions can actually tell about how jets lose energy in the quark-gluon plasma. It shows that the centrality-dependent jet $R_{AA}$ for lead-lead collisions at 5.02 TeV is adequately reproduced by a model with only two parameters, both governing the mean transverse-momentum loss. Comparing fits that use a gamma distribution and a delta function for the loss shows the data cannot distinguish the shape of the loss distribution, only its mean; the gamma shape parameter is poorly constrained. Adding collision geometry from the TRENTo model makes the centrality dependence sensitive to the jet formation time, with linear path-length weighting preferring $\\tau_f = 0.9$ fm/c and quadratic weighting preferring $\\tau_f = 0.1$ fm/c. Finally, the paper shows that assumptions about how systematic errors correlate across momentum and centrality change the best-fit parameters, so such assumptions must be specified in model-to-data comparisons.","feed_headline":"Jet data fix mean energy loss, not its shape","feed_subtitle":"Centrality data fit a two-parameter model; formation-time preference depends on path-length weighting.","key_machinery":"The load-bearing object is the two-parameter mean energy-loss function $\\langle \\Delta p_T\\rangle = \\alpha (p_T)^\\beta \\log(p_T)$, which shifts the measured proton–proton jet spectrum down in transverse momentum. The same mean is inserted into two competing loss distributions, a gamma distribution (motivated as a convolution of $k$ exponential scatterings) and a delta function, to test shape sensitivity. For the centrality extension, the model takes initial geometry from the TRENTo model, with parameters fixed from an earlier Bayesian fit, and weights the medium energy density along each dijet path by an inverse-time factor starting at a formation time $\\tau_f$, using either a linear or quadratic path-length weighting. The analysis also constructs covariance error matrices from the published systematic-error components, assuming different correlation structures, and uses a log-likelihood transformation to address the known tendency of fits with strongly correlated relative errors to undershoot the data.","core_discovery":"The central claim is that measurements of jet $R_{AA}$ in lead-lead collisions, by themselves, constrain only the mean jet energy loss $\\langle \\Delta p_T \\rangle = \\alpha (p_T)^\\beta \\log(p_T)$, not the shape of the energy-loss distribution. Fitting the same mean loss with a gamma distribution and with a delta function yields equivalent agreement with the measured 0–10% central data, and the gamma shape parameter $k$ is not well constrained. The paper therefore concludes that $R_{AA}$ data fix the mean loss but that an observable such as dijet asymmetry is needed to determine the distribution's shape. Extending the model to all centralities with TRENTo geometry, the centrality dependence of $R_{AA}$ is reproduced, and the preferred formation time shifts from 0.9 fm/c for linear path-length weighting to 0.1 fm/c for quadratic weighting. The treatment of systematic errors also matters: constructing the covariance matrix by summing components versus by multiplying sums, and evaluating the likelihood on the data or its logarithm, both change the most-probable parameters and $\\chi^2_d$ values.","pith_inferences":["A natural extension is to fit the same two-parameter model to dijet asymmetry alongside $R_{AA}$, which would test whether the gamma shape parameter $k$ becomes constrained; the paper leaves exactly this as future work.","The formation-time shift with path-length weighting suggests that a single effective formation time extracted from $R_{AA}$ is model-dependent; future studies that treat jet and medium formation times separately could reinterpret the current $\\tau_f$ preference as an effective combination.","The correlated-error analysis implies that experimental papers should report per-component systematic covariance matrices; without them, model fits may systematically undershoot or overshoot the measured points, independent of the physics model."],"forward_implications":["If $R_{AA}$ alone fixes only the mean loss, then model comparisons should report $\\langle \\Delta p_T\\rangle(p_T)$ rather than the full loss distribution, and discriminating between energy-loss mechanisms will require observables such as dijet asymmetry.","Under the model, formation time and path-length weighting trade off: longer $\\tau_f$ pairs with linear weighting and shorter $\\tau_f$ with quadratic weighting, so the two cannot be separated using centrality-dependent $R_{AA}$ alone.","Fits with highly correlated systematic errors can land systematically below the data; the log-likelihood transformation mitigates but does not fully remove this, and residual shape mismatch remains the likely explanation.","The assumption about systematic-error correlations (summed components vs. full covariance) changes the most-probable parameters and $\\chi^2_d$, so any model-to-data comparison must state the correlation assumption explicitly."],"supporting_citations":[{"why":"provides the Pb+Pb jet $R_{AA}$ data (with per-component systematic errors) that are the target of all fits.","marker":"[23]"},{"why":"supplies the proton–proton inclusive jet cross-section used to set the unquenched baseline.","marker":"[27]"},{"why":"introduces the TRENTo model that supplies the initial collision geometry for the centrality-dependent extension.","marker":"[16]"},{"why":"fixes the TRENTo parameters to values from a prior Bayesian fit to LHC data.","marker":"[12]"},{"why":"introduces the original mean $\\Delta p_T$-loss parameterization that this model extends.","marker":"[20]"},{"why":"provides the Markov-chain Monte Carlo sampling used for the Bayesian calibrations.","marker":"[28]"},{"why":"documents the correlated-error undershooting bias and the log-likelihood transform used in the fits.","marker":"[29]"}],"fun_headline_variants":["Jet data fix mean loss, not loss distribution","Centrality R_AA fits simple two-parameter model","Mean jet energy loss pinned; shape unconstrained","Formation time flips with path-length assumption","Systematic error model changes jet-fit results"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The formation-time preference depends on the TRENTo geometry parameters being fixed correctly and on the assumed linear or quadratic path-length weighting; if either is wrong, the inferred $\\tau_f$ values would change.","fun_headline_variants_meta":{"raw":{"variants":["Jet data fix mean loss, not loss distribution","Centrality R_AA fits simple two-parameter model","Mean jet energy loss pinned; shape unconstrained","Formation time flips with path-length assumption","Systematic error model changes jet-fit results"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000255,"raw_usage":{"total_tokens":1589,"prompt_tokens":977,"completion_tokens":612,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":593,"completion_tokens_details":{"reasoning_tokens":540}},"tokens_in":593,"tokens_out":612,"duration_ms":6833,"temperature":1.0,"reasoning_tokens":540,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T22:08:58.903629+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Vary the TRENTo initial-state parameters within their posterior uncertainties and refit the centrality-dependent $R_{AA}$: if the optimal $\\tau_f$ no longer shifts from 0.9 fm/c under linear path weighting to 0.1 fm/c under quadratic weighting, the formation-time conclusion fails. Alternatively, a precise dijet-asymmetry measurement that cannot be described by a delta-distribution loss with the same mean would falsify the claim that $R_{AA}$ alone constrains only the mean energy loss.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"introduces the TRENTo model that supplies the initial collision geometry for the centrality-dependent extension."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"documents the correlated-error undershooting bias and the log-likelihood transform used in the fits."}],"review_version":1}