{"id":"f7fbb655-822a-4c4a-ab97-acb8a3bf7daf","arxiv_id":"2411.10539","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Working order-by-order in the strong coupling, the paper constructs the likelihood ratio for boosted H->gg jets versus QCD jets and validates a new color-flow observable, d2/z^2, that improves signal over background by roughly a factor of 150 in Pythia.","lead":"This paper derives the optimal discrimination observable, order by order in the strong coupling, for separating boosted Higgs-boson decays to gluons from ordinary QCD jets. It shows that a mass-and-subjet-energy cut gives a fixed signal-to-background gain, and introduces a new observable that improves the ratio by roughly two orders of magnitude in simulation.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Perfect-discriminant proof covers only g→gg background; inclusive QCD claim (including quark jets) is not established at NLO or asymptotically.","rationale":"The reader's verdict is CONDITIONAL and its weakest-assumption field already identifies the g→gg restriction and the reliance on generator simulation. My stress-test sharpens this into a concrete coverage gap: the analytic proof of the perfect discriminant in Sec. 5.1 is explicitly restricted to massive jets from g→gg fragmentation, whereas the abstract and conclusions state the result for QCD jets broadly, and the simulated sample used to validate the several-hundred-fold improvement contains quark-initiated jets as well. Because the simulation stops at p_T > 2 TeV, it cannot verify the asymptotic p_T → ∞ behavior for the inclusive background; if quark jets have a less strongly p_T-suppressed L distribution, the background at very high p_T would be quark-dominated and the perfect-discriminant limit would not hold for inclusive QCD. This does not invalidate the finite-p_T simulation results or the LO analysis, which are separately well supported, so I do not change the reader's CONDITIONAL verdict; I recommend maintaining that condition until the quark-jet scaling is checked. The secondary point about the x→0 contribution in Eqs. (49)-(51) is also worth noting: the paper asserts that this limit vanishes without presenting the full integral, but even if that assertion is wrong, the fixed-background-efficiency scaling L ~ p_T^4 would likely survive because the argument of the double log still involves p_T^2/√L; hence I do not make this the primary concern.","tokens_in":18458,"tokens_out":18959,"duration_ms":192224,"concrete_test":"Use the MadGraph/Pythia samples to label background jets by the flavor of the initiating hard parton before showering. In each p_T bin (500, 1000, 1500, 2000 GeV), compute the L_cut value at a fixed background efficiency (e.g., 1%) separately for quark-only and gluon-only subsets, and test whether each threshold grows like p_T^4/m_H^4. In parallel, recompute the quark-jet analog of Eq. (51) using the q→qg dipole eikonal to check whether its leading-log scaling has the same p_T dependence; if quark-only thresholds grow more slowly, the inclusive 'perfect discriminant' claim fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing gap is between the NLO analytic derivation and the inclusive-QCD claim. The 'perfect discriminant' conclusion in Sec. 5.1 is derived only for the g→gg background: Eq. (33) uses the g→ggg soft matrix element, Eq. (47) uses the g→gg splitting function, and the scaling L ~ R^4 p_T^4/m_H^4 in Eq. (53) follows from that channel's leading double log. Quark-initiated background jets are never analyzed at NLO; the simulation validation in Fig. 4 mixes quark and gluon jets and stops at p_T > 2 TeV, so it cannot establish the p_T → ∞ claim for inclusive QCD jets. If quark jets have a different soft-emission color structure (e.g., q→qg dipoles), their L distribution may not sharpen as p_T^4 at fixed background efficiency, and the inclusive background at very high p_T would be dominated by the least-suppressed jet flavor, invalidating 'perfect discrimination' for QCD jets generally. In addition, the leading-double-log estimate in Eqs. (49)-(51) keeps only the x→∞ limit and asserts the x→0 contribution vanishes; the full O(α_s) integral in Eq. (48) is not evaluated, so the coefficient of log^2 and the p_T^4 scaling are not verified analytically.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a systematic, matrix-element-based approach to discriminating boosted H→gg decays from massive QCD jets. At leading order, using a jet-mass window and the softer-subjet momentum fraction z, it proves that the signal-to-background improvement is independent of the jet kinematics at high boosts and is approximately 1/αs(m_H). At next-to-leading order, it constructs IRC-safe observables from the soft-emission matrix elements, in particular d2 (the inverse of the H→ggg matrix element) and the ratio L=d2/z^2, and argues that L becomes a perfect discriminant as p⊥→∞. The theoretical results are compared with MadGraph+Pythia simulations for p⊥ > 500, 1000, 1500, and 2000 GeV.","tokens_in":18689,"tokens_out":11709,"duration_ms":106406,"significance":"The leading-order analysis is clean, and the comparison with simulation in Fig. 2 is convincing. The construction of an IRC-safe observable from the inverted signal matrix element is a useful and well-motivated idea that connects anomaly detection to analytic QCD, and the paper explicitly derives the functional form of the observable rather than fitting it to data. However, the headline asymptotic claim—perfect discrimination for inclusive QCD jets—rests on a restricted g→gg background analysis and on a leading-double-log estimate that contains an unjustified zero for the small-angle (x→0) contribution. If the x→0 term is included, the scaling L ~ p⊥^4/m_H^4 may be weakened, which would eliminate the perfect-discriminant conclusion. The paper's finite-p⊥ results are still valuable, but the central asymptotic claim needs substantial revision.","major_comments":[{"comment":"The claim in Eq. (50) that the x→0 limit of the integral in Eq. (48) vanishes is incorrect. For x→0, the integrand behaves as (1/x) log[x^2/(z^2(1-z)L)] with a theta function that sets a lower bound x > z√((1-z)L) (up to O(1) factors), while the jet-radius constraint in Eq. (47) supplies an upper bound x < √(z(1-z)) R p⊥/m_H. The resulting integral is (1/2) log^2[ R^2 p⊥^2/(m_H^2 z L) ] plus terms of order log^2 L, and is not zero; it is a leading double logarithm of the same parametric order as the x→∞ contribution kept in Eq. (49). Consequently, the approximation leading to Eq. (51) is not a controlled leading-double-log estimate, and the scaling L ~ R^4 p⊥^4/m_H^4 in Eq. (53) is not established by the given calculation. The full integral in Eq. (48) should be evaluated, or at least estimated with both end-point limits retained, before drawing the perfect-discriminant conclusion.","section":"§5.1, Eqs. (48)–(51)"},{"comment":"The perfect-discriminant conclusion is derived only for the g→gg background channel: Eq. (33) uses the g→ggg soft matrix element, and Eq. (47) uses the g→gg splitting function. Quark-initiated jets are never analyzed at NLO, despite the Introduction's statement that the high-boost problem reduces to discrimination from 'color triplet or color octet jets of QCD' and the abstract's unqualified 'perfect discriminant' claim. Since quark jets can dominate the inclusive background at large transverse momentum, the asymptotic claim for inclusive QCD jets is unsupported. The simulation in Fig. 4 includes both quark and gluon jets but stops at p⊥ > 2 TeV, so it cannot establish the p⊥→∞ behavior. The paper should either explicitly restrict the asymptotic claim to the g→gg subchannel or extend the color analysis to quark-initiated backgrounds, for example by treating the q→qg dipole.","section":"§5.1 and Introduction"},{"comment":"The conclusion from Eq. (53) that L=d2/z^2 has 'arbitrarily good signal efficiency with arbitrarily low background efficiency' as p⊥→∞ is drawn from a fixed-order, leading-double-log approximation of the background cumulative distribution. The author notes that this expansion is only meaningful where the O(αs) correction is below unity, but the asymptotic statement requires the cumulative distribution to be valid at arbitrarily small background efficiencies, where αs log^2(...) is no longer small. A resummation or a higher-order argument is needed to justify the perfect-discriminant limit, and the finite-p⊥ simulation cannot validate it. The paper should either provide such an argument or explicitly weaken the claim to a scaling prediction for the background tail at fixed, moderate background efficiencies.","section":"§5.1, Eq. (53)"}],"minor_comments":[{"comment":"The sentence 'p⊥/mH≪ 1' should read 'm_H/p_\\(\\perp\\) ≪ 1', since the highly-boosted regime is defined by the jet mass being much smaller than its transverse momentum.","section":"§2"},{"comment":"There is a typo: 'signficant' should be 'significant'.","section":"§2"},{"comment":"The word 'simualted' should be 'simulated'.","section":"§4.1"},{"comment":"In the paragraph before Eq. (42), 'an observable formed from their combination that performs between than either individually' should read 'performs better than either individually'.","section":"§5.1"},{"comment":"The left-panel legend entry for the observable (1+O_NLO)/z is garbled in the rendered text; please ensure the math displays correctly.","section":"Fig. 3"},{"comment":"The phrase 'hadronc top decay' should be 'hadronic top decay'.","section":"§6"},{"comment":"Reference [48] is listed as 'E∞ Scheme, unpublished'; if this is the intended source for Winner-Take-All recombination, please provide a full reference or remove it.","section":"Ref. [48]"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within scope for SciPost Physics, and the leading-order analysis is a solid contribution with a convincing simulation comparison. The main concern is that the NLO asymptotic claim is not currently supported: the small-angle integral in Eq. (50) is evaluated incorrectly, and the analysis does not cover quark-initiated backgrounds. These issues are fixable, but they change the strength of the conclusions, so I recommend major revision rather than acceptance. The author should also be encouraged to distinguish more carefully between the g→gg subchannel result, which is well-supported, and the inclusive-QCD claim, which is only validated at finite p⊥."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a solid, mostly honest analytic benchmark paper. The d2/z^2 observable is new, the construction is well explained, and the simulation comparison is genuine. The main caveat is the 'perfect discriminant' claim: it is derived only for g→gg background using leading double logarithms, and the paper is more confident in that limit than the derivation supports.\n\nWhat's actually new: a first NLO likelihood construction for H→gg versus QCD jets, including the observable d2, and clean analytic ROC expressions at LO. The inverse-matrix-element/anomaly-detection idea is a nice trick and is made concrete. The paper extends the author's earlier program, but the H→gg case has its own wrinkles, and the simulation validation with MadGraph+Pythia is independent enough that the central result—that d2/z2 is a strong discriminant—is credible.\n\nThe soft spots: First, the NLO analysis restricts the background to g→gg fragmentation. That's stated in Sec. 5.1, but the intro and the 'perfect discriminant' language in Sec. 5.1 go beyond what that calculation shows. Quark-initiated jets are never analyzed at NLO. Simulation at pT>2TeV includes quark jets and still shows strong discrimination, so the empirical result is fine, but the asymptotic claim for inclusive QCD is not established analytically. Second, the pT^4 scaling is derived by keeping only leading double logarithms and discarding the x→0 region without evaluating the full integral. That is a parametric estimate, not a proof. The paper is fairly transparent about the approximation, but the phrase 'perfect discriminant' overstates it. Minor: no statistical uncertainties on the ROC/S/B curves, no public code or data, and the abstract's 'several hundred' is a bit strong—the text gives ~150 with the mass cut and ~30 from d2/z2 alone at reasonable signal efficiency.\n\nFor whom: anyone working on jet substructure or on ML tagger benchmarks. It's not going to change experimental searches, and the paper says so. The LO part is a nice pedagogical example of likelihood-ratio construction; the NLO part is a useful baseline for what a simple IRC-safe observable can do.\n\nI'd send it to review. It deserves referee time; the gaps I named can be fixed with a more careful statement of what is proven and what is estimated, and maybe a footnote on quark jets. I would not desk reject.","headline":"A clean analytic benchmark for H→gg tagging, but the perfect-discriminant claim is a leading-log estimate for g→gg jets, not an inclusive-QCD proof.","tokens_in":19263,"tokens_out":6468,"would_cite":true,"duration_ms":60867,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A simple infrared-safe jet observable constructed from the inverse of the H→ggg matrix element becomes an essentially perfect discriminator between boosted Higgs decays and QCD jets in the high-energy limit, and yields a…","keywords":["boosted Higgs decays","H→gg tagging","jet substructure","likelihood ratio","infrared and collinear safety","anomaly detection","color flow","soft gluon emission"],"falsifier":"Measure the cumulative distribution of $L = d_2/z^2$ on QCD background jets at several jet transverse momenta (for example $p_\\perp = 1$, $2$, and $3$ TeV) after a Higgs mass window. The paper predicts that the value of $L$ at fixed background efficiency scales as $p_\\perp^4/m_H^4$; if that scaling is absent, or if the simulated signal-to-background improvement of several hundred does not appear in data, the central claim is falsified.","tokens_in":18218,"feed_emoji":"🎯","tokens_out":10219,"duration_ms":89934,"temperature":0.7,"pith_summary":"This paper asks how well a boosted Higgs boson decaying to two gluons can be separated from ordinary QCD jets that fake a Higgs, and answers with a systematic expansion of the likelihood ratio in the strong coupling. At leading order, it proves that a jet mass window plus the energy fraction of the softer of the two leading subjets improves signal over background by roughly $1/\\alpha_s(m_H) \\sim 10$, independent of jet kinematics at high boost. At next-to-leading order, it constructs the observable $d_2/z^2$ by inverting the $H\\to ggg$ matrix element with soft gluon emission; the inversion turns the matrix element's soft and collinear divergences into zeros, making the observable infrared and collinear safe. The paper argues analytically that this observable approaches the perfect discriminant as the jet transverse momentum grows, and shows in simulated data from standard event generators and parton showering that it improves the signal-to-background ratio by a factor of several hundred at $p_T > 2$ TeV. If correct, a simple three-emission measurement approaches the likelihood-ratio bound for $H\\to gg$ tagging and provides a concrete physics benchmark for machine-learning jet taggers.","feed_headline":"One jet measurement nears perfect Higgs tagging at high boost","feed_subtitle":"A mass cut plus the ratio d2/z2 improves signal over background by several hundred at jet pT above 2 TeV.","key_machinery":"The load-bearing object is the inverse of the next-to-leading-order $H\\to ggg$ matrix element, packaged as an IRC-safe jet observable $d_2 = \\frac{p_\\perp^2}{m_H^2}\\sum_k z_k \\frac{\\theta_{1k}^2\\theta_{2k}^2}{\\theta_{12}^2}$; the full discriminant is $L = d_2/z^2$, where $z$ is the softer-subjet momentum fraction. Its power comes from color flow: the Higgs is a color singlet, so its soft gluon emission concentrates between the two hard gluons, whereas a QCD gluon jet is a color octet and emits at wide angles. Inverting the signal matrix element turns the soft and collinear poles into zeros, which is what makes the observable infrared and collinear safe and therefore calculable order-by-order in $\\alpha_s$. The $z^2$ denominator compensates for the background's strong phase-space enhancement at small $z$, and the overall normalization makes the observable boost invariant along the jet.","core_discovery":"The central discovery is that the likelihood ratio for boosted $H\\to gg$ jets versus QCD jets can be constructed order-by-order in $\\alpha_s$ and, at next-to-leading order, reduces to a simple IRC-safe observable. Working in the collinear, high-boost limit, the paper computes the soft-gluon matrix elements for a color-singlet $H\\to gg$ signal and for color-octet $g\\to gg$ background. The difference of these matrix elements is proportional to $\\theta_{1k}^2 + \\theta_{2k}^2 - \\theta_{12}^2$, which is positive where a QCD gluon jet emits and negative where the color-singlet Higgs emits. Even simpler, the inverse of the signal matrix element defines $d_2 = \\frac{p_\\perp^2}{m_H^2}\\sum_k z_k \\frac{\\theta_{1k}^2\\theta_{2k}^2}{\\theta_{12}^2}$, which is IRC safe because soft and collinear divergences of the matrix element become zeros. The paper proves that, after a fixed jet mass window, the ratio $L = d_2/z^2$ has a background cumulative distribution whose fixed-efficiency contour scales as $L \\propto p_\\perp^4/m_H^4$, so as $p_\\perp\\to\\infty$ signal efficiency can grow while background efficiency stays fixed, meaning the observable becomes the perfect discriminant. In simulated events with $p_\\perp > 2$ TeV, cuts on this observable improve signal over background by a factor of several hundred.","pith_inferences":["The author leaves implicit that the same color-flow asymmetry should tag any color-singlet resonance decaying to two partons at high boost, not just the Higgs decaying to gluons.","The simulation validation uses default showering and hadronization settings; varying those models or adding underlying-event and pileup conditions would test whether the several-hundred-fold improvement survives in a more detector-like environment, a check the paper does not perform.","A natural extension of the paper's moment analysis would be to construct the full two-dimensional likelihood on the $(d_2, z)$ space rather than the ratio $d_2/z^2$; the paper's moments suggest the ratio is near-optimal, but the full likelihood could be computed from the same simulated samples."],"forward_implications":["After a Higgs mass window, a cut on the softer subjet energy fraction alone reduces background by about $1/\\alpha_s(m_H) \\sim 10$, independent of jet kinematics at high boost.","On simulated jets with $p_\\perp > 2$ TeV, the observable $d_2/z^2$ improves the signal-to-background ratio by a factor of several hundred over inclusive jet selection.","The inverse-matrix-element construction is IRC safe in general, so the same anomaly-detection recipe can produce calculable discriminants for other signal processes whose next-to-leading-order matrix elements are known.","The analytic ROC curves give a ground reference that machine-learning taggers for $H\\to gg$ should reproduce at minimum, which the paper frames as a step toward interpretability of jet taggers."],"supporting_citations":[{"why":"Supplies the Neyman-Pearson lemma, the principle that the likelihood ratio is the optimal discrimination observable that organizes the whole construction.","marker":"[24]"},{"why":"Provide the collinear splitting functions used to compute the leading-order signal and background energy-fraction distributions.","marker":"[25–29]"},{"why":"Provide the eikonal soft-gluon amplitudes from which the next-to-leading-order signal and background matrix elements are built.","marker":"[50–53]"},{"why":"Establishes the method of expanding the likelihood ratio through next-to-leading order that this paper applies and extends.","marker":"[49]"},{"why":"Supplies the energy-correlation power-counting construction and the D2 observable after which d2 is named.","marker":"[58]"},{"why":"Constructs an optimal color-singlet identification observable, the direct precursor to the color-flow-sensitive d2 construction.","marker":"[56]"},{"why":"Generates the simulated signal and background event samples used to validate the analytic predictions.","marker":"[35]"},{"why":"Showers and hadronizes the generated events, producing the jet samples on which the observables are tested.","marker":"[36]"},{"why":"Clusters final-state particles into the anti-kT jets analyzed in the validation.","marker":"[37]"},{"why":"Supply the grooming energy-fraction distributions and the Sudakov-safety framework underlying the leading-order likelihood.","marker":"[39,40]"}],"fun_headline_variants":["Mass cut plus d2/z2 boosts Higgs signal by hundreds at high pT","One simple ratio separates Higgs jets from QCD at high boost","Analytic proof: d2/z2 gives near-perfect Higgs tagging","Higgs tagger: one ratio, several hundred-fold boost","At high boost, a simple observable nears perfect Higgs tagging"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The analytic prediction that $d_2/z^2$ becomes a perfect discriminant assumes that background QCD jets at high transverse momentum are dominated by collinear gluon splitting with soft wide-angle gluon emission, and that hadronization does not erase the color-flow difference; the paper tests this only in simulation, not in collider data.","fun_headline_variants_meta":{"raw":{"variants":["Mass cut plus d2/z2 boosts Higgs signal by hundreds at high pT","One simple ratio separates Higgs jets from QCD at high boost","Analytic proof: d2/z2 gives near-perfect Higgs tagging","Higgs tagger: one ratio, several hundred-fold boost","At high boost, a simple observable nears perfect Higgs tagging"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00082,"raw_usage":{"total_tokens":3657,"prompt_tokens":1083,"completion_tokens":2574,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":699,"completion_tokens_details":{"reasoning_tokens":2482}},"tokens_in":699,"tokens_out":2574,"duration_ms":20089,"temperature":1.0,"reasoning_tokens":2482,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T19:35:41.405515+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the cumulative distribution of $L = d_2/z^2$ on QCD background jets at several jet transverse momenta (for example $p_\\perp = 1$, $2$, and $3$ TeV) after a Higgs mass window. The paper predicts that the value of $L$ at fixed background efficiency scales as $p_\\perp^4/m_H^4$; if that scaling is absent, or if the simulated signal-to-background improvement of several hundred does not appear in data, the central claim is falsified.","supporting_citations":[{"cited_title":"Binary Discrimination Through Next-to-Leading Order","cited_arxiv_id":"2309.14417","evidence_quote":"Establishes the method of expanding the likelihood ratio through next-to-leading order that this paper applies and extends."}],"review_version":1}