{"id":"6ab559d6-514f-4853-b7c3-60a2f7bceacc","arxiv_id":"1909.01639","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Using a neural network trained on simulated jets, the ALICE collaboration measures jet spectra and nuclear modification factors in Pb-Pb collisions down to 40 GeV/c and up to jet radius 0.6.","lead":"ALICE tests a machine-learning correction for jet momentum in lead-lead collisions and reports jets at lower momenta and larger radii than previously possible. The new results overlap with standard methods where they can be compared, and extend measurements to a new regime.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Toy-background training with an exponential high-pT tail and a fragmentation systematic that changes only the response matrix leave the low-pT and R=0.6 claims without a direct closure test.","rationale":"The reader's weakest assumption identified the general risk that simulation training does not represent real Pb-Pb background and fragmentation. My concern agrees with that risk but makes it more specific: the toy background's exponential high-pT tail is unrealistic, and the fragmentation systematic only varies the response matrix, not the estimator itself. These are load-bearing because the novelty of the paper is the extension to lower pT and to R=0.6, where no area-based cross-check exists. The paper does have independent support: the residual-width comparisons in Sec. 3.3 and the compatibility with the area-based method at higher pT are real checks, and the method is a direct extension of a peer-reviewed method paper. But those checks do not close the low-pT/R=0.6 extrapolation. The concern is therefore not grounds for rejection, but it reinforces the need for an embedding closure test and a fragmentation systematic that includes retraining or at least re-evaluating the estimator. The conditional verdict remains appropriate; no verdict change is needed.","tokens_in":4914,"tokens_out":6867,"duration_ms":76766,"concrete_test":"Perform a closure test by embedding PYTHIA jets with known true pT into real Pb-Pb events (the Sec. 3.3 procedure) and compare, binwise in the lowest reported pT bins (for example 40-60 GeV/c in 0-10% and R=0.6), the mean ML-corrected pT minus true pT normalized by true pT. Require this residual to be compatible with zero within the quoted systematic uncertainties, and repeat using a quark-jet fragmentation sample. If the residual deviates significantly, the low-pT and R=0.6 claims are not supported by the presented validation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim depends on the ML mapping learned in Sec. 3.2 from PYTHIA jets embedded in a thermal toy background generalizing to real Pb-Pb events. The toy background in Sec. 3.1 has an exponentially falling momentum distribution above roughly 4 GeV/c, whereas real Pb-Pb events contain a power-law tail of high-pT particles from semi-hard and hard processes. Because the input features include the pT of the first eight leading jet constituents, a high-pT background or underlying-event particle in a real event can be misidentified by the estimator as part of the signal jet, producing a positive bias in the corrected pT. The validation shown in Sec. 3.3 reports only the width of residual distributions for embedded probes; it does not establish that the mean residual is unbiased in the lowest pT bins or for R=0.6. Furthermore, the fragmentation systematic described in Sec. 4 only replaces the response matrix used in unfolding with a quark-jet response. It does not retrain or re-evaluate the ML estimator itself under altered fragmentation, even though the estimator is trained on PYTHIA jets and its input-output mapping can depend on fragmentation. The advertised low-pT reach and the first R=0.6 measurement are precisely the region where no cross-check against the area-based method exists, so the compatibility seen at higher pT cannot validate the extrapolation.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This conference proceeding reports a machine-learning-based jet momentum correction for Pb-Pb collisions at sqrt(s_NN)=5.02 TeV, applied to track-based jets with R=0.2, 0.4, and 0.6. A neural network is trained on PYTHIA jets embedded in a thermal background to map raw jet features to an estimate of the true detector-level jet pT on a jet-by-jet basis. The authors present nuclear modification factors for R=0.4 (down to 40 GeV/c in 0-10% and 30 GeV/c in 30-50% central collisions) and, for the first time, for R=0.6, together with jet cross-section ratios sigma(R=0.2)/sigma(R=0.4) and sigma(R=0.2)/sigma(R=0.6). They claim compatibility with the area-based estimator at higher pT, no significant R-dependence of R_AA, and no centrality or pp deviation in the cross-section ratios.","tokens_in":5265,"tokens_out":3427,"duration_ms":33592,"significance":"If the ML correction is unbiased, the method would extend ALICE jet measurements to lower pT and larger R than the standard area-based method, and the R=0.6 result would be a first in heavy-ion collisions at the LHC. The paper contains genuine validation steps: an embedding-based comparison of residual widths against the area-based estimator (Fig. 1) and a direct comparison of R_AA at R=0.4 in two centrality classes (Fig. 2). The work is a proceedings contribution and the underlying method is documented in a separate paper (ref. [4]), which is appropriately cited. However, the new physics claims (low-pT reach and first R=0.6 measurement) rest on the untested assumption that a model trained on a simplified thermal background generalizes to real Pb-Pb events; the evidence presented in this manuscript is not sufficient to establish that. The central results are therefore plausible but not yet convincingly validated.","major_comments":[{"comment":"The training background in Sec. 3.1 uses an exponential momentum tail above roughly 4 GeV/c, whereas real Pb-Pb events contain a power-law tail of high-pT particles from semi-hard and hard processes. Because the input features include the transverse momenta of the first eight leading jet constituents, a high-pT background or underlying-event particle in a real event could be misidentified by the estimator as part of the signal jet, producing a positive bias in the corrected pT. The embedding validation in Sec. 3.3 reports only the width of residual distributions; it does not show that the mean residual is consistent with zero, especially in the lowest pT bins or for R=0.6. To support the advertised low-pT reach, the authors should provide a closure test with a realistic background model that includes a power-law component, or at least report the mean residual as a function of pT for each R and centrality.","section":"Sec. 3.1 and 3.3"},{"comment":"The fragmentation systematic variation replaces only the response matrix used in unfolding with a quark-jet response. It does not retrain or re-evaluate the ML estimator itself under altered fragmentation, even though the estimator is trained on PYTHIA jets and its input-output mapping is likely sensitive to jet fragmentation. Changing only the unfolding response cannot capture a fragmentation-dependent bias in the estimator itself. The authors should either retrain the estimator on samples with modified fragmentation and compare the final corrected spectra, or provide an argument based on the estimator's features as to why it is insensitive to fragmentation.","section":"Sec. 4, fragmentation systematic"},{"comment":"The compatibility with the area-based method shown in Fig. 2 is established only in the pT region where the area-based method is already reliable. The new claims -- jets at 40 (30) GeV/c for R=0.4 and the R=0.6 measurement -- lie precisely in the region where no such cross-check exists. The cross-section ratios in Fig. 4 and the R=0.6 R_AA in Fig. 3 therefore rest entirely on the ML estimator's extrapolation. To support these claims, the authors should provide an alternative validation in this region, for example a closure test on full heavy-ion simulations with realistic underlying events, or a comparison with an independent background subtraction technique such as constituent subtraction or soft-drop grooming.","section":"Sec. 4, Figs. 2-4"}],"minor_comments":[{"comment":"The abstract and introduction state that transverse momentum spectra will be presented, but the figures show only nuclear modification factors and cross-section ratios; the spectra themselves are not displayed. Please clarify what is shown or add the spectra.","section":"Abstract and Sec. 1"},{"comment":"The caption of Fig. 1 is incomplete: the left panel is described as 'the comparison of the different background estimators for R=0.4' but the reader cannot identify which curves correspond to which estimator or centrality. The right panel lacks axis labels and a legend. Please improve the figure captions and labels.","section":"Fig. 1"},{"comment":"The systematic uncertainties are described only qualitatively; no numeric values are given for any observable. For a measurement-oriented proceeding, at least the dominant systematic uncertainties at representative pT values should be stated.","section":"Sec. 4"},{"comment":"The regression target is defined as the reconstructed jet momentum multiplied by the momentum fraction carried by PYTHIA particles in the jet. This definition is confusing and should be clarified: it appears to approximate the true detector-level jet momentum, but the potential bias introduced by this target choice is not discussed. Please rephrase and add a sentence on the validity of this approximation.","section":"Sec. 3.2"},{"comment":"Equation (4.1) is not fully typeset: the differentials and symbols are inconsistently formatted (e.g., d2Nch jet/dpT,ch jetdηjet). Also, the definition T_AA = N_coll/sigma_tot is a rough shorthand; please use the standard definition and cite a reference.","section":"Sec. 4, Eq. (4.1)"}],"recommendation":"major_revision","confidential_remarks":"This is a conference proceeding, and the central method is documented in a separate PRC paper [4]. The main issue is that the new physics results (low-pT and R=0.6) are not validated by a direct closure test, and the fragmentation systematic does not cover the estimator's own training dependence. These are fixable within the scope of a revised proceedings by adding a closure test on a realistic background simulation and retraining-based fragmentation study, or by explicitly restricting the claims. The paper is otherwise clearly written and the comparison with the area-based method is a strength."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, the headline: this is a conference proceeding that presents the first physics results from the ML jet-background estimator introduced in [4]. Those results—R=0.6 jet R_AA in Pb-Pb, and R=0.4 down to 40/30 GeV/c—are genuinely new. The method itself isn't, and the authors don't pretend it is.\n\nWhat the paper does well: the ML-corrected R_AA is shown head-to-head with the area-based method in the overlapping pT range, and it matches. That's the most convincing evidence that the correction isn't just adding a constant. The systematic list is standard and includes a fragmentation variation, which is more than many proceedings include. The embedding validation on real Pb-Pb events is a reasonable check that the estimator transfers from the toy training background to real events.\n\nThe soft spots are the ones you'd expect. The training background has an exponential tail above 4 GeV/c; real Pb-Pb events have a power-law tail. Since the estimator uses the eight leading jet constituents, a hard background particle inside the jet cone can fool it. The validation reports the width of residual distributions, not the mean residual as a function of pT and R, so a low-pT bias wouldn't show up. And the fragmentation systematic replaces the response matrix used in unfolding but doesn't retrain or re-evaluate the ML estimator itself, even though the estimator was trained on PYTHIA jets and is sensitive to fragmentation. Those are the two places where the advertised low-pT and R=0.6 reach rely on extrapolation without a direct closure test.\n\nThat said, this is a proceedings. It's short, the plots are the substance, and the central physics claims are modest—no big R-dependence, consistent with hybrid model. I don't think these concerns sink it, but they should be addressed in a full paper: mean residual plots, an alternative background or fragmentation training sample, and numerical systematics.\n\nThe audience is the jet-quenching community. I'd take it seriously as a conference proceedings; if it grows into a journal submission, it deserves a real peer review, and a referee with an ML background would add value. I'd be more likely to cite the method paper than this one, but I'd bring this to a reading group as a concrete example of ML in heavy-ion jet reconstruction.","headline":"Solid proceedings with first Pb-Pb results from the ML jet correction; credible where checkable, but low-pT and R=0.6 claims rest on a training background and fragmentation systematic that are not fully closed.","tokens_in":5733,"tokens_out":3528,"would_cite":false,"duration_ms":35920,"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":"Machine-learning background correction extends jet measurements in Pb-Pb to record-low momenta and to R=0.6 jets.","keywords":["jet reconstruction","heavy-ion collisions","machine learning","background fluctuations","jet transverse momentum","nuclear modification factor","Pb-Pb collisions","jet radius dependence"],"falsifier":"Embed simulated events with a known true jet $p_{\\mathrm{T}}$ but with a fragmentation pattern deliberately different from the PYTHIA training sample (for example, quark-only jets or jets with strong medium-induced energy loss), run the trained ML estimator on them, and test whether the residual between corrected and true jet momentum shifts systematically with jet $p_{\\mathrm{T}}$ or radius $R$; a shift that the response-matrix and unfolding procedure cannot remove would show the estimator is not robust to the fragmentation assumption.","tokens_in":4706,"feed_emoji":"🎯","tokens_out":10511,"duration_ms":96694,"temperature":0.7,"pith_summary":"This paper argues that a machine-learning estimator, trained on simulated jets embedded in a thermal background, can correct the transverse momentum of each jet individually in Pb-Pb collisions, reducing the large background fluctuations that limit the standard area-based method. Using this jet-by-jet correction, the ALICE detector measures track-based jets at $\\sqrt{s_\\mathrm{NN}}=5.02$ TeV down to $40$ GeV/$c$ in 0-10% central collisions and $30$ GeV/$c$ in 30-50% central collisions for resolution parameter $R=0.4$, and reports the first measurement of $R=0.6$ jets in heavy-ion collisions at the LHC. The resulting nuclear modification factors are compatible with the area-based method and show no significant dependence on jet radius. The payoff of the approach is that it lowers the momentum threshold and extends the radius reach of jet measurements, opening a part of phase space that studies of jet quenching in heavy-ion collisions could not previously access.","feed_headline":"ML background correction extends jets to record-low momentum in Pb-Pb","feed_subtitle":"The new jet-by-jet estimator allows the first R=0.6 jet measurements in heavy-ion collisions at the LHC.","key_machinery":"The machinery is a supervised regression model, in this analysis a neural network with three hidden layers of 100, 100, and 50 neurons, that takes a set of per-jet observables and outputs a single corrected jet $p_{\\mathrm{T}}$. The input features are the jet $p_{\\mathrm{T}}$ corrected by the standard area-based method, the jet angularity, the number of jet constituents, and the transverse momenta of the eight hardest constituents; the training target is the detector-level true jet momentum, approximated as the reconstructed jet momentum times the PYTHIA-carried momentum fraction. The trained model is applied jet by jet, and a response matrix built from embedding vacuum jets into real Pb-Pb background handles any residual smearing and detector effects before an unfolding step produces the final spectra.","core_discovery":"The central discovery is that a neural-network regressor can learn the mapping from raw, background-contaminated jet observables to the true jet transverse momentum well enough to replace the statistical area-based background subtraction with a per-jet correction. The paper demonstrates this by training the model on PYTHIA jets reconstructed to detector level and embedded in a thermal background, with the regression target defined as the reconstructed jet momentum multiplied by the momentum fraction carried by PYTHIA particles in the jet. The trained estimator is then applied to Pb-Pb data, and the residual fluctuations are removed by unfolding through a response matrix built from embedding vacuum jets into real data backgrounds. The result is that track-based jets are measured down to $p_{\\mathrm{T}}=40$ GeV/$c$ (0-10% central) and $30$ GeV/$c$ (30-50%) for $R=0.4$, and $R=0.6$ jets are measured for the first time in heavy-ion collisions at the LHC, with spectra and nuclear modification factors compatible with the area-based method and with no significant $R$ dependence.","pith_inferences":["Because the estimator is trained on the PYTHIA momentum fraction, its jet-by-jet correction inherits the fragmentation model; training instead on a mixture of quark- and gluon-initiated jets or on jets with varied energy loss could make the correction less model-dependent.","The same regression approach could be pushed to $R=0.8$ or $0.9$ within the ALICE acceptance, exactly the regime where the area-based method's background fluctuations are most severe and where the ML method's advantage would be largest.","The jet-by-jet idea is not limited to jet $p_{\\mathrm{T}}$: similar features could estimate groomed jet momentum, subjet momentum, or the background under a jet, which would extend the method to jet-structure observables in heavy-ion collisions."],"forward_implications":["Jets with resolution parameter $R=0.6$ can be measured in Pb-Pb collisions at the LHC for the first time, extending jet quenching studies to larger angular scales.","Track-based jets at $R=0.4$ become measurable down to $40$ GeV/$c$ in 0-10% and $30$ GeV/$c$ in 30-50% central Pb-Pb collisions, a much larger low-$p_{\\mathrm{T}}$ reach than the area-based method.","The ML-corrected nuclear modification factors agree with the area-based estimator where both are available, confirming the method's consistency with an established baseline.","No significant $R$ dependence of the nuclear modification factor is observed between $R=0.2$, $0.4$, and $0.6$, and the jet cross-section ratios are consistent with the pp baseline."],"supporting_citations":[{"why":"introduces the ML-based jet momentum estimator and its toy-model training on which this analysis is built.","marker":"[4]"},{"why":"quantifies the event background momentum density and its fluctuations in Pb-Pb collisions.","marker":"[2]"},{"why":"area-based background correction and unfolding procedure used as the comparison baseline.","marker":"[3]"},{"why":"PYTHIA8 simulation used to generate the jets for training and for the response matrix.","marker":"[5]"},{"why":"GEANT3 detector simulation used to reconstruct embedded jets at detector level.","marker":"[6]"},{"why":"pp charged-jet cross-section data used as the reference for the nuclear modification factor.","marker":"[7]"},{"why":"anti-$k_{\\mathrm{T}}$ jet clustering algorithm used for jet reconstruction.","marker":"[9]"},{"why":"hybrid strong/weak-coupling model calculation to which the measured jet suppression is compared.","marker":"[12]"}],"fun_headline_variants":["ML jet correction unlocks first R=0.6 jets in heavy-ion collisions","Neural network beats background, extends jets to record-low momentum in Pb-Pb","Per-jet ML boost: first measurement of R=0.6 jets in Pb-Pb","AI-based jet momentum reconstruction reaches new low pT in ALICE Pb-Pb","Machine learning lifts jet reconstruction to new lows in Pb-Pb data"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The method's reliability rests on the assumption that jets embedded in a simplified thermal background, trained with the PYTHIA momentum fraction as the target, represent the true relationship between the measured jet features and the actual jet momentum in real Pb-Pb collisions; if heavy-ion jet fragmentation or background structure differs from this training model in ways the embedding tests do not capture, the corrected jet spectra will be biased.","fun_headline_variants_meta":{"raw":{"variants":["ML jet correction unlocks first R=0.6 jets in heavy-ion collisions","Neural network beats background, extends jets to record-low momentum in Pb-Pb","Per-jet ML boost: first measurement of R=0.6 jets in Pb-Pb","AI-based jet momentum reconstruction reaches new low pT in ALICE Pb-Pb","Machine learning lifts jet reconstruction to new lows in Pb-Pb data"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000897,"raw_usage":{"total_tokens":3872,"prompt_tokens":961,"completion_tokens":2911,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":577,"completion_tokens_details":{"reasoning_tokens":2806}},"tokens_in":577,"tokens_out":2911,"duration_ms":19616,"temperature":1.0,"reasoning_tokens":2806,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:11:21.321901+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Embed simulated events with a known true jet $p_{\\mathrm{T}}$ but with a fragmentation pattern deliberately different from the PYTHIA training sample (for example, quark-only jets or jets with strong medium-induced energy loss), run the trained ML estimator on them, and test whether the residual between corrected and true jet momentum shifts systematically with jet $p_{\\mathrm{T}}$ or radius $R$; a shift that the response-matrix and unfolding procedure cannot remove would show the estimator is not robust to the fragmentation assumption.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"quantifies the event background momentum density and its fluctuations in Pb-Pb collisions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"area-based background correction and unfolding procedure used as the comparison baseline."},{"cited_title":"Sjostrand et al.: PYTHIA 6.4 – Physics and Manual, JHEP 0605 (2006) 026","cited_arxiv_id":null,"evidence_quote":"PYTHIA8 simulation used to generate the jets for training and for the response matrix."},{"cited_title":"Brun et al.: GEANT Detector Description and Simulation Tool, CERN Program Library Long Writeup CERN-W-5013 (1994)","cited_arxiv_id":null,"evidence_quote":"GEANT3 detector simulation used to reconstruct embedded jets at detector level."},{"cited_title":"Cacciari, G.P","cited_arxiv_id":null,"evidence_quote":"anti-$k_{\\mathrm{T}}$ jet clustering algorithm used for jet reconstruction."},{"cited_title":"Casalderrey-Solana et al., A hybrid strong/weak coupling approach to jet quenching, JHEP 10 (2014) 019","cited_arxiv_id":null,"evidence_quote":"hybrid strong/weak-coupling model calculation to which the measured jet suppression is compared."}],"review_version":1}