{"id":"a733b612-74d0-4a77-9f13-1e917da84fee","arxiv_id":"2505.14542","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"The authors introduce matrix-element-based cluster fission and decay models and show which angularity and energy-correlation observables are most sensitive to each piece.","lead":"This paper proposes new physics-based building blocks for the cluster hadronization model used in the Herwig event generator, covering how clusters split into lighter clusters and decay into hadrons. It also identifies jet shape variables and energy correlations that can isolate these two stages, which may help future precision tuning at the LHC and future colliders.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The fission matrix element's assumed momentum-conserving continuation of the soft limit is unvalidated, and the observed Belle improvement is not yet shown to be robust to corrections or to other untuned model parameters.","rationale":"The reader identified the momentum-conservation/soft-limit assumption in Sec. 2.4 as the central weakness, and my independent re-reading of the paper confirms that this is the most load-bearing point. The strongest claim (improved Belle description via new fission matrix element) depends entirely on this untested continuation of a soft-factorized formula into full phase space. I also add a concrete, falsifiable test: compare the factorized formula to the exact matrix element of the same diagrams, and check stability of the Belle comparison under a small retune of the three default Herwig parameters. No ad hominem is intended; the authors are explicit about the assumption, but the paper currently presents it without a validation. I agree with the conditional verdict: the model-development and observable-design contribution is real and should be published, but the central claim should be re-scoped or supported by the checks above.","tokens_in":21224,"tokens_out":1501,"duration_ms":15422,"concrete_test":"Implement two independent checks in Herwig: (1) Evaluate the fission matrix element of Eqs. (9)-(12) against the fully momentum-conserving evaluation of the same Feynman diagrams (compute |M|^2 from the two diagrams in Fig. 6 with all constituent masses included, using the 2->4 phase space of Eq. (7)) and quantify the ratio pointwise in (M1,M2,cos(theta)); if the factorized soft formula deviates by large factors in regions that dominate the BELLE z~0.9-0.95 bin, the improved data description is not robustly tied to the claimed building block. (2) Retune or minimally vary Herwig's default parameters (e.g. Clmax, Clpow, PSplit) with the new fission matrix element and re-run the Fig. 11 observables; if the plateau can be reintroduced or the Belle improvement washes out under a small retune, the conclusion that the matrix element fixes the plateau is premature.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper's central claim is that the proposed matrix elements are theoretically motivated building blocks; the key load-bearing step is the factorization of the cluster-fission matrix element in Sec. 2.4: |M_{2->4}|^2 = A0 * (|M_{2->2}|_t^2 / t_denom) * S(q1,q2,q,qbar), built from the Catani-Grazzini soft approximation, but then applied over the full phase space with momentum-conserving on-shell kinematics. The authors explicitly acknowledge (Sec. 2.4) that the soft limit 'would normally not provide us with such an expression' and that they assume energy-momentum conservation can be implemented. This is precisely the assumption on which the removal of the unphysical plateau and the Belle data improvement (Figs. 11-12) rest; the paper provides no check of whether subleading-power (1/Q) corrections, the recoil scheme, or finite-mass corrections beyond the q.barq -> q.barq + m^2 replacement would preserve the qualitative improvement. Because the new building blocks are then fed into an untuned model with the t-channel cluster decay, the comparison to Belle data is not a standalone validation of the fission matrix element.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes new building blocks for the cluster hadronization model in Herwig, focusing on cluster fission and cluster decay. For fission, the authors factorize the phase space and introduce a matrix element built from a t-channel gluon exchange and a soft q-qbar function (Eqs. 9-12), evaluated with on-shell constituent masses and exact momentum conservation. For decay, they introduce a t-channel-like matrix element with a pseudo-mass parameter M_S. The model is implemented in Herwig and compared with DELPHI and Belle data; the new fission model removes an unphysical plateau in the cluster mass distribution and improves the Belle di-hadron invariant-mass spectra (Figs. 11-12). The paper also proposes angularity and energy-correlation observables intended to discriminate between fission and decay dynamics, but these are shown only as model-variant comparisons, not compared to data.","tokens_in":21562,"tokens_out":10303,"duration_ms":96340,"significance":"If the proposed fission matrix element is a valid continuation of the infrared-factorized soft limit, the work is a meaningful step toward a more principled cluster hadronization model, and the external Belle comparison provides a valuable anchor for the improvement claim. The clean phase-space factorization in Sec. 2.3 is a useful contribution in itself, and the proposed observables in Sec. 4 could be genuinely discriminating for future tuning. However, the central 'phenomenological constraints' advertised in the title are not yet delivered: the new observables are validated only against the model itself, and the key fission matrix element rests on an explicitly acknowledged assumption that is not stress-tested. The improvement over the default model is real but is not yet isolated to the matrix element as opposed to the corrected phase space.","major_comments":[{"comment":"The load-bearing step of the paper is the continuation of the infrared-factorized soft matrix element to full phase space with momentum-conserving on-shell kinematics. The authors state explicitly that the soft limit 'would normally not provide us with such an expression' and that they assume energy-momentum conservation can be implemented. This assumption underlies the removal of the plateau in Fig. 12 and the Belle improvement in Fig. 11, yet no estimate or validation of subleading-power corrections is provided. The only robustness check mentioned, varying the gluon constituent mass, is not performed. I ask for at least one concrete test, such as a comparison with the full tree-level 2-to-4 matrix element for the diagrams in Fig. 6, or a scan over the regulator parameter epsilon and m_g, to show that the qualitative improvement is not accidental to the chosen continuation.","section":"Sec. 2.4, Eqs. (9)-(12)"},{"comment":"The title and abstract promise 'phenomenological constraints' and observables with 'constraining power on the individual building blocks,' but none of the observables studied in Sec. 4 is compared to experimental data. All conclusions there are drawn from comparisons among Herwig variants with different fission/decay models, which demonstrates sensitivity but not constraint. To support the advertised claims, the authors should either overlay data from DELPHI or Belle for the angularities and energy correlators (at least where they exist), or explicitly reframe Sec. 4 as a model-diagnostic study for future tuning and adjust the title and abstract accordingly.","section":"Sec. 4, Figs. 14-21"},{"comment":"The improvement over the Herwig default is presented for an otherwise untuned model, as the authors acknowledge ('still untuned' in Sec. 3.2). Moreover, the 'CF phase space' variant (blue) also removes the plateau, so the specific role of the new matrix element, as opposed to the corrected phase space alone, is not cleanly isolated. A quantitative comparison of the CF phase-space and CF matrix-element variants (for example, chi-square per bin for the Belle distributions in Figs. 11 and 12) is needed to support the claim that the matrix element itself is the relevant building block.","section":"Sec. 3.2, Figs. 11-12"}],"minor_comments":[{"comment":"The text says the masses M1 and M2 are sampled uniformly in the 'triangular phase space of Fig. 5,' but the allowed region in Fig. 5 is not triangular due to the Källén-function thresholds; please rephrase.","section":"Sec. 2.5, step 2"},{"comment":"The fitted proposal-distribution parameters A, beta1, beta2 and the overestimate lambda_OE used for the rejection sampling are not reported; including their values, or a link to the implementation, would make the algorithm reproducible.","section":"Sec. 2.5, Eq. (15)"},{"comment":"The pseudo-mass choice M_S = max{(m1-mh1),(m2-mh2)} is introduced without derivation; since the authors label it a model assumption and the LEP/Belle comparisons show no new tensions, this is acceptable but should be flagged more visibly as a phenomenological input.","section":"Sec. 2.6, cluster decay matrix element"},{"comment":"The claim that the STA correlation is approximately independent of sqrt(s) is supported only by visual inspection; a quantitative ratio with uncertainties would be more convincing.","section":"Sec. 4.2, Fig. 21"},{"comment":"There are several typographical and formatting issues, including 'infrafred' in Sec. 5, 'ap1p2-dipole' in Sec. 2.4, and stray '/Bullet' artifacts in Fig. 3 and elsewhere; these should be cleaned up.","section":"Whole paper"}],"recommendation":"major_revision","confidential_remarks":"The main gap is between the title's promise of phenomenological constraints and the actual content: Sec. 4 contains no data comparisons, and the central fission matrix element is an explicitly assumed continuation of the soft limit without a robustness check. This is fixable within the manuscript's scope by adding sensitivity studies and either adding data comparisons or retitling. The paper is otherwise a solid contribution to the hadronization-modeling literature and the external Belle comparison is a genuine strength."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The one thing to know: this paper is a real step forward for the cluster hadronization model, but it does not deliver what the title promises. There are no phenomenological constraints in the fitting sense—the new observables in Sec. 4 are only shown as model variations, never compared to data—and the central fission matrix element rests on an assumption the authors themselves flag in Sec. 2.4. The soft-limit factorization is used with full momentum conservation and on-shell constituent masses, which the soft limit \"would normally not provide.\" No check of subleading-power corrections, recoil scheme, or finite-mass effects is given, so the Belle improvement could in principle be fragile.\n\nWhat is genuinely new and good: the paper formulates cluster fission and decay as Lorentz-invariant matrix elements rather than ad hoc longitudinal or isotropic kinematics. The t-channel ansatz and the Catani–Grazzini-based soft function are concrete and clearly specified, with a detailed rejection-sampling algorithm. The Belle comparison (Figs. 11–12) shows the new fission model removes the unphysical plateau in the cluster mass distribution and describes the di-hadron invariant mass data qualitatively better than the Herwig default. That is a real, externally anchored result, not just a re-description of the model. The WTA/STA energy-correlation observables in Sec. 4 are also a useful new toolkit for dissecting hadronization stages; the fact that STA correlations are insensitive to cluster decay while WTA correlations are sensitive to both is a genuinely informative diagnostic.\n\nWhere it is soft: the stress-test concern holds up on reading. The momentum-conservation assumption is load-bearing, and the paper does not test its robustness. The \"phenomenological constraints\" language is overclaimed: no fit is performed, the new observables are not compared to data, and the model is left untuned, so the mid/high-mass discrepancies in Fig. 11 are attributed to tuning but never checked. The pseudo-mass M_S in the cluster decay is ad hoc, though the authors are upfront about it. No code is released, which limits reproducibility, but the algorithm description is detailed enough to reimplement.\n\nWho should read it: anyone working on hadronization models, event generators, or precision fragmentation at e+e− and LHC. It deserves serious peer review. I would send it out, asking the authors to either test the soft-limit assumption's robustness (e.g., vary the gluon-mass regulator and recoil scheme) or re-scope the claims to a model-development and observable-design study. The central idea is sound enough to warrant that investment.","headline":"Genuine progress on cluster hadronization building blocks with a real Belle improvement, but the title overpromises and the fission matrix element's key momentum-conservation assumption goes untested.","tokens_in":22043,"tokens_out":2892,"would_cite":true,"duration_ms":38160,"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":"The paper proposes perturbatively motivated matrix elements for cluster fission and decay that remove an unphysical plateau in the cluster mass distribution and improve low-energy e+e− di-hadron data.","keywords":["cluster hadronization","cluster fission","cluster decay","soft factorization","matrix element building blocks","di-hadron invariant mass","jet angularities","energy correlators"],"falsifier":"A precise measurement of the di-hadron invariant mass distribution in the highest fractional-energy bin at 10.58 GeV would settle it: if the matrix-element fission, after a complete tune with colour reconnection included, does not remove the low-mass plateau or does not match the turn-on shape, the central improvement claim fails.","tokens_in":20992,"feed_emoji":"⚛️","tokens_out":9627,"duration_ms":90678,"temperature":0.7,"pith_summary":"Hadronization in event generators is a major source of systematic uncertainty for precision measurements, yet the cluster models behind it remain largely phenomenological. This paper argues that the two central steps of cluster hadronization—cluster fission and cluster decay—can be built from theoretically motivated matrix elements rather than ad hoc power laws and fixed angular prescriptions. The proposed fission matrix element, a soft quark-antiquark emission with a t-channel gluon exchange, removes an unphysical plateau in the cluster mass distribution and improves the description of low-energy di-hadron invariant mass data relative to the default model. The paper also constructs angularity and energy-correlator observables that separate the effects of fission from decay, giving a way to constrain each building block individually.","feed_headline":"Cluster fission matrix elements erase a hadronization plateau","feed_subtitle":"These matrix elements fix the mass plateau and split observables into fission- and decay-sensitive probes.","key_machinery":"The central object is the factorized $2\\to4$ phase space $f_{\\rm PS}(M_1,M_2)$ built from Källén functions, combined with the soft-emission matrix element of Eqs. (9)–(12). The fission matrix element carries a t-channel gluon exchange regulated by a gluon constituent mass and a soft function $S(q_1,q_2,q,\\bar q)$ constructed from eikonal factors $I_{ij}$; it simultaneously drives the cluster mass distribution toward lower masses and produces collinear angular distributions that smoothly match the parton shower. The cluster decay uses the t-channel-like matrix element $1/[(p_1-h_1)^2 - M_S^2]^2$, replacing isotropic decay with a kinematics that continues the fission picture. These objects define the new building blocks that the paper tests against data.","core_discovery":"At its core, the paper claims that cluster fission should be viewed as a low-scale, perturbatively motivated continuation of the parton shower rather than a purely longitudinal splitting with power-law-distributed cluster masses. It writes the fission rate as a factorized phase space times a tree-level soft quark-antiquark emission matrix element, $|\\mathcal{M}_{2\\to4}|^2 = A_0 |\\mathcal{M}_{2\\to2}|^2_t\\,S(q_1,q_2,q,\\bar q)$, where the t-channel gluon denominator is regulated by a gluon constituent mass and the soft function is built from eikonal factors. This matrix element shifts the cluster mass distribution away from the default triangular behaviour, eliminating the flat plateau that is not seen in the measured di-hadron spectra. For cluster decay the paper proposes a t-channel-like hadron matrix element proportional to $1/[(p_1-h_1)^2 - M_S^2]^2$ with $M_S = \\max\\{(m_1-m_{h_1}),(m_2-m_{h_2})\\}$, which smoothly interpolates between the fission and hadron kinematics. Together these building blocks, with default parameters otherwise untouched, improve the description of the di-hadron invariant mass distribution at B-factory energies and introduce no new tensions in high-energy event-shape observables.","pith_inferences":["If these matrix elements survive a global tune, the hadronization start could be defined at a factorization scale rather than at a shower cutoff, making predictions less dependent on where the shower stops.","The observed energy independence of the soft-takes-all correlation could be turned into a direct test of hadronization universality: any measured energy dependence of that correlator would signal energy-dependent hadronization.","The decay matrix element's pseudo-mass parameter $M_S$ is set by a kinematical maximum; a natural extension would be to promote it to a fitted form factor and constrain it with the decay-sensitive winner-takes-all correlations.","Since colour reconnection is not yet included, the extracted sensitivities of the discriminating observables may shift once it is added, so the observables should be re-run in the unified model before finalizing the building blocks."],"forward_implications":["The default cluster mass distribution plateau near the constituent-mass threshold disappears when the soft matrix element drives fission, so di-hadron invariant mass spectra at low energies match data without new tuning.","Cluster fission can be treated as a smooth, perturbative continuation of the parton shower rather than a longitudinal split with power-law masses, reducing the number of tunable parameters in the hadronization model.","Infrared-dangerous soft-takes-all energy correlations are nearly insensitive to cluster decay, while winner-takes-all correlations are sensitive to both fission and decay, so the two probes together isolate the two building blocks.","The t-channel-inspired cluster decay leaves established high-energy event-shape distributions without new tensions while improving some regions.","The same building-block logic extends to colour reconnection and cluster propagators in a future unified hadronization model."],"supporting_citations":[{"why":"Defines the default cluster hadronization model and its parameters that the new building blocks replace and are compared against.","marker":"[4]"},{"why":"Supplies the infrared-safety and factorization argument that hadronization should continue the parton shower smoothly, motivating the fission matrix element.","marker":"[18]"},{"why":"Previous matching study that identified cluster fission as a continuation of the shower and exhibits the triangular cluster-mass distribution and plateau problem.","marker":"[19]"},{"why":"Provides the infrared factorization of tree-level QCD amplitudes from which the eikonal soft function for the fission matrix element is taken.","marker":"[28]"},{"why":"The e+e- di-hadron invariant-mass measurement used as the main low-energy comparison that the new fission model improves.","marker":"[33]"},{"why":"Defines the 10.58 GeV centre-of-mass energy used for the phase-space plots and the low-energy simulations.","marker":"[27]"},{"why":"Describes baryon production from cluster hadronisation and colour rearrangement, providing the comparison point for baryon correlation observables.","marker":"[24]"},{"why":"Provides the soft-gluon-evolution basis for the planned colour-reconnection extension and for subtracting soft-gluon effects from reconnection.","marker":"[20]"},{"why":"The high-energy e+e- event-shape and momentum data used to check that the decay-model variation introduces no new tensions.","marker":"[31]"}],"fun_headline_variants":["Fission matrix elements kill the hadronization plateau","Cluster fission: a soft-emission fix for hadronization","New building blocks sharpen hadronization predictions","Fission replaces longitudinal splitting in hadronization"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the simplified soft-emission formula for cluster fission continues to be correct when applied to real clusters with full momentum conservation and massive constituents, even though that formula is normally derived only for very low-energy emissions.","fun_headline_variants_meta":{"raw":{"variants":["Fission matrix elements kill the hadronization plateau","Cluster fission: a soft-emission fix for hadronization","New building blocks sharpen hadronization predictions","Fission replaces longitudinal splitting in hadronization"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000217,"raw_usage":{"total_tokens":1407,"prompt_tokens":891,"completion_tokens":516,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":507,"completion_tokens_details":{"reasoning_tokens":454}},"tokens_in":507,"tokens_out":516,"duration_ms":5813,"temperature":1.0,"reasoning_tokens":454,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T15:31:54.918373+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A precise measurement of the di-hadron invariant mass distribution in the highest fractional-energy bin at 10.58 GeV would settle it: if the matrix-element fission, after a complete tune with colour reconnection included, does not remove the low-mass plateau or does not match the turn-on shape, the central improvement claim fails.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the 10.58 GeV centre-of-mass energy used for the phase-space plots and the low-energy simulations."},{"cited_title":"Abreu et al., Tuning and test of fragmentation models based on identified particles and precision event shape data , Z","cited_arxiv_id":null,"evidence_quote":"The high-energy e+e- event-shape and momentum data used to check that the decay-model variation introduces no new tensions."}],"review_version":1}