{"id":"a3d4619e-58b5-41d0-b95a-870137c4502d","arxiv_id":"1908.07429","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Adding the g* g* to g c cbar process to a kT-factorization improved color evaporation model makes the JH-2013 predictions for prompt J/psi transverse momentum agree with LHCb data.","lead":"Using two existing QCD frameworks together, the authors find that adding a 2-to-3 gluon process fixes the high-momentum shape of J/psi production for one family of gluon distributions. It is a short conference note that shows a working recipe, not a parameter-free prediction, for quarkonium production at the LHC.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 2→3 term's differential shape is unvalidated; the claimed JH-2013 improvement may be an artifact of p0_T and unstated scales.","rationale":"The reader's weakest_assumption already identifies the suppression factor and tree-level normalization as the structurally weak premise, and my independent reading agrees. I do not see a reason to change the CONDITIONAL verdict: the concern strengthens the need for the stated conditions (report uncertainties, specify all scales and masses) rather than invalidating the central claim outright. The paper's self-described 'somewhat arbitrary' suppression factor is an explicit limitation that should have been stress-tested in the manuscript; the proposed test would settle whether the JH-2013 agreement survives reasonable variations of p0_T and of unstated scales. The missing direct-to-prompt reference '[?]' is a normalization nuisance that is largely absorbed by the fitted P_J/ψ and is not load-bearing for the shape claim.","tokens_in":6468,"tokens_out":10107,"duration_ms":105771,"concrete_test":"Recompute Fig. 3 (top right) with p0_T = 1.0 and 2.0 GeV, re-fitting P_J/ψ to the low-pT LHCb normalization in each case, and with two scale choices (e.g., μ_F = m_T(c cbar) vs μ_F = 2 m_c) and charm masses (m_c = 1.3 vs 1.5 GeV). If the high-pT curve moves outside the LHCb uncertainties or the good description is lost, the improvement is parameter- or scale-driven. As an independent check, compute the open-charm D-meson pT distribution with the same JH-2013+2→3 setup and p0_T = 1.5 GeV and compare to 7 TeV data, where no J/ψ transition probability is fitted.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (Fig. 3, top panels) is that adding the tree-level g∗g∗→g c cbar process to g∗g∗→c cbar with the JH-2013-set2 UGDF yields a very good description of prompt J/ψ pT distributions. For this to hold, the 2→3 term must provide the correct differential shape of the c cbar pair pT, not just the correct integrated rate. The only regularization of the soft/collinear region is Fsup(pT)=pT^4/((p0_T)^2+pT^2)^2 with p0_T=1.5 GeV, 'adjusted to the exact NLO calculations of the charm production cross section at the LHC' (Sec. 2). That adjustment constrains the normalization, not the shape; in the pT≈4–14 GeV region where the improvement is visualized, Fsup is already 0.77–0.96, so the predicted slope is essentially the bare tree-level off-shell 2→3 matrix element. No comparison to an independent differential NLO calculation or to open-charm data is shown for the JH-2013+2→3 combination. Moreover the renormalization and factorization scales and the charm-quark mass are not stated, so the claimed 'very good description' cannot be reproduced or tested from the paper alone. The qualitative direction is unsurprising—an extra recoiling gluon produces a harder spectrum—but the quantitative claim is not yet robust.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript applies the improved color evaporation model (ICEM) to prompt J/psi production at LHC energies, computing the underlying c-cbar-pair cross section in kT-factorization. Two hard subprocesses are included: the standard g*g* -> c cbar channel and, as the new element, the 2 -> 3 channel g*g* -> g c cbar evaluated with the KATIE generator. The model is compared with ALICE and LHCb rapidity and pT distributions at 7 TeV using two unintegrated gluon distributions, the KMR-CT14lo/nlo UGDFs and the CCFM-based JH-2013-set2 UGDF. The main claim is that adding the 2 -> 3 process to the JH-2013-set2 UGDF changes a poor high-pT slope into a very good description of the LHCb data, whereas the KMR-based calculation already works with only the 2 -> 2 mechanism provided a kt < mu_F cut is imposed to avoid double counting.","tokens_in":6753,"tokens_out":2843,"duration_ms":30986,"significance":"If the central claim is correct, the paper offers a practical and conceptually useful way to incorporate additional hard emissions in CCFM-based UGDFs within kT-factorization, and it sharpens the interpretation of why KMR and JH-2013 UGDFs differ for correlation-type observables. The use of KATIE for the tree-level off-shell 2 -> 3 matrix elements is a concrete and reproducible technical step, and the explicit discussion of the kt > mu_F region in the KMR prescription is a genuine conceptual contribution. At the same time, the quantitative support for the headline claim is visual only, and two parameters controlling the normalization (P_J/psi and, indirectly, p0_T) are fitted to data entering the comparison, so the paper's significance rests on the shape improvement rather than on an absolute prediction.","major_comments":[{"comment":"The central claim that g*g* -> g c cbar 'completely changes the picture' and gives 'a very good description' of the pT distributions for the JH-2013-set2 UGDF is not quantitatively validated. The suppression factor Fsup(pT) = pT^4 / ((p0_T)^2 + pT^2)^2 is adjusted to the total charm production cross section, which constrains the normalization but not the differential shape; in the pT ~ 4-14 GeV region shown in Fig. 3, Fsup is already 0.77-0.96, so the predicted slope is essentially the bare tree-level off-shell 2 -> 3 matrix element. No comparison to an independent differential NLO calculation or to open-charm pT data is shown for the JH-2013 + 2 -> 3 combination, so the paper does not currently demonstrate that the improvement is robust rather than a consequence of the unvalidated 2 -> 3 shape. I recommend adding such a differential validation or, failing that, explicitly presenting the scale and mass dependence as an uncertainty band.","section":"Section 2, paragraph 3; Section 3, Fig. 3"},{"comment":"The parameter P_J/psi is not predicted but fitted to the LHCb data at small pT for each UGDF (P_J/psi = 0.0065 for JH-2013-set2 and 0.018 for KMR-CT14lo). Since the comparison in Figs. 2 and 3 is to the same data set, the absolute normalization is not a prediction. The pT-slope comparison is partly independent because P_J/psi only rescales all bins, but the statement that the model gives 'a very good description' should be qualified by this fitting procedure, and the paper should state explicitly which observables are genuinely predicted.","section":"Section 2, Eq. (2.2); Section 3, Fig. 2"},{"comment":"The paper does not specify the renormalization and factorization scales, the charm-quark mass, or the precise implementation of the kt < mu_F cut used for the KMR UGDF in the 2 -> 3 process. These choices are load-bearing for the normalization and for the claimed absence of double counting, and without them the numerical results cannot be reproduced or independently tested. The authors should state all input parameters and, ideally, show the sensitivity of the pT distributions to their variation.","section":"Section 2, paragraph 3; Section 3, Fig. 3"},{"comment":"The direct-to-prompt ratio of 0.62 is introduced with a missing reference (the text shows '[?]'). Since this factor directly multiplies the comparison with prompt-J/psi data, it is a substantive input and must be referenced or derived. If it is taken from a specific fit or measurement, the source should be cited; if it is a phenomenological choice, its uncertainty should be propagated into the comparison.","section":"Section 2, Eq. (2.2); Section 3"}],"minor_comments":[{"comment":"There are several typographical errors, including 'valueable' in the Introduction and 'respecitvely' in Section 3; these should be corrected in a revision.","section":"Throughout"},{"comment":"The theoretical curves are shown without uncertainty bands or goodness-of-fit measures such as chi2, even though the central comparison is visual. Adding a quantitative measure would substantially strengthen the paper.","section":"Figures 2 and 3"},{"comment":"The label 'NLO alpha_s + kt < mu_F' in the bottom panels is difficult to parse; the caption should explain that the cut applies to the 2 -> 3 contribution with the KMR UGDF.","section":"Section 3, Fig. 3 caption"}],"recommendation":"major_revision","confidential_remarks":"This is a short proceedings-style paper, so the technical scope is necessarily limited, but the main claim is currently supported mainly by visual inspection and by a fitted normalization. The missing scale settings and the unreferenced direct-to-prompt ratio are easy to fix, whereas the differential validation of the 2 -> 3 shape is the substantive issue that should be addressed before the paper is accepted. I do not see a fundamental error in the approach, but the manuscript needs additional quantitative support or a more cautious wording of the central claim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the paper's main point is plausible and worth knowing. In kT-factorized ICEM, the JH-2013 CCFM UGDF gives too steep a J/psi pT spectrum from the 2->2 g*g*->ccbar term; adding an explicit 2->3 g*g*->g ccbar process through KATIE flattens the slope and matches LHCb. For the KMR UGDF, the 2->2 term already works because KMR includes the kt > mu_F region, and the 2->3 contribution needs a dedicated cut to avoid double counting. That is a coherent and physically sensible story, and the authors are upfront that the ICEM-plus-kT-factorization baseline was independently done in Ref. [9]. The specific numerical implementation of the 2->3 term and the systematic UGDF comparison is the new piece.\n\nWhat the paper does well: it uses a real tree-level MC (KATIE), it frames the UGDF difference in terms of where hard emissions live, and it flags the double-counting issue explicitly rather than sweeping it under the rug. The conclusions are stated cleanly and the plots tell a visually consistent story.\n\nSoft spots: the stress test holds up. p0_T = 1.5 GeV is tuned to the total charm cross section, so it fixes normalization, not shape. In the pT = 4-14 GeV region where the JH-2013 improvement is claimed, the suppression factor is already 0.77-0.96, meaning the predicted slope is essentially the bare tree-level 2->3 matrix element. There is no independent differential check against open-charm data or an NLO calculation for that combination. Renormalization and factorization scales and the charm quark mass are not stated, so the \"very good description\" cannot be reproduced from the text. The direct-to-prompt ratio 0.62 appears as \"[?]\" with no reference. Finally, all comparisons are visual; no uncertainty bands or chi2 values. These are real limitations, but they are limitations of a proceedings-length paper, not signs of a broken calculation.\n\nWho this is for: anyone working on quarkonium production in kT-factorization, especially those using CCFM-based UGDFs. The qualitative message that 2->3 processes matter for the pair pT spectrum in that framework is useful even before the quantitative details are pinned down.\n\nRecommendation: I would send this to peer review. It deserves referee time as a legitimate phenomenological contribution, but I would ask the authors to state all scales and masses, resolve the missing reference, and show at least one independent differential validation of the 2->3 shape. If those are added, the claim would be much more robust.","headline":"Plausible qualitative fix for the JH-2013 pT slope in kT-factorized ICEM, but the quantitative claim rests on visual comparison and unreported scales.","tokens_in":7309,"tokens_out":1784,"would_cite":true,"duration_ms":19609,"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":"Adding the 2-to-3 gluon-fusion channel to the color evaporation model restores the LHC J/psi transverse-momentum spectrum.","keywords":["J/psi production","color evaporation model","kT-factorization","unintegrated gluon distributions","2 to 3 processes","transverse momentum distributions","LHC","heavy quarkonium"],"falsifier":"Recompute the JH-2013 plus $2\\to3$ prediction at $\\sqrt{s}=13$ TeV, or with $p_T^0$ varied between 1.0 and 2.0 GeV, and compare with LHCb transverse-momentum data; if the good 7 TeV agreement depends strongly on this tuned parameter, the regularization is carrying the result rather than the $2\\to3$ mechanism.","tokens_in":6172,"feed_emoji":"⚛️","tokens_out":14945,"duration_ms":127406,"temperature":0.7,"pith_summary":"The paper argues that the poor high-transverse-momentum tail of $J/\\psi$ production predicted by the $k_T$-factorization color evaporation model with CCFM-based unintegrated gluon distributions is not a failure of the model but a missing emission channel. Adding the $2\\to3$ process $g^* g^* \\to g c\\bar c$ at tree level to the standard $2\\to2$ process $g^* g^* \\to c\\bar c$ turns a poor description into a very good one for the JH-2013 gluon distribution. For the KMR gluon distribution no such addition is needed, because its construction already includes hard gluon emissions with transverse momentum above the factorization scale. The result matters because it shows how to compare two popular unintegrated gluon distributions within one quarkonium-production approach and where higher-order emissions are hidden.","feed_headline":"Adding a hard gluon emission fixes LHC J/psi momentum spectrum","feed_subtitle":"The 2 to 3 channel g* g* to g c cbar restores the high-pT slope for CCFM gluons; KMR gluons already include it.","key_machinery":"The machinery is a $k_T$-factorization calculation of $c\\bar c$-pair production followed by the improved color evaporation model mapping to $J/\\psi$. The two hard subprocesses are the $2\\to2$ fusion $g^* g^* \\to c\\bar c$, using the standard analytic off-shell matrix element, and the $2\\to3$ process $g^* g^* \\to g c\\bar c$, computed at tree level with off-shell matrix elements from a numerical Monte Carlo generator. The two unintegrated gluon distributions play different roles: KMR-type distributions include transverse momenta above the factorization scale and therefore contain hard emissions internally, while JH-2013 CCFM-type distributions contain only soft emissions and need the $2\\to3$ term added externally. A suppression factor $F_{\\rm sup}(p_T)=p_T^4/((p_T^0)^2+p_T^2)^2$ with $p_T^0=1.5$ GeV regularizes the low-$p_T$ minijet region, and a direct-to-prompt correction of 0.62 converts the prediction to the measured prompt $J/\\psi$ yield.","core_discovery":"The central claim is that the steep fall of the JH-2013 CCFM-based prediction for the $J/\\psi$ transverse-momentum distribution is caused by the omission of the $k_t>\\mu_F$ region in that unintegrated gluon distribution, and that the missing physics is restored by computing the $2\\to3$ subprocess $g^* g^* \\to g c\\bar c$ explicitly. In the $k_T$-factorization realization of the improved color evaporation model, the $2\\to2$ term alone gives a spectrum that falls much faster than the LHCb data; adding the $2\\to3$ tree-level contribution, regularized by a suppression factor with $p_T^0=1.5$ GeV, produces a very good description of both the rapidity and transverse-momentum distributions at $\\sqrt{s}=7$ TeV. The same addition leaves the KMR-based result essentially unchanged, because the KMR gluon distribution already contains the hard-emission tail, and for KMR a cut $k_t<\\mu_F$ is imposed on the $2\\to3$ term to avoid double counting. The paper therefore identifies the difference between the two unintegrated gluon distributions as the real origin of the earlier discrepancy.","pith_inferences":["The same reasoning should apply to other quarkonia such as $\\Upsilon$ and to correlation observables like $D\\bar D$ azimuthal decorrelations: wherever a CCFM-type unintegrated gluon distribution is used, an explicit hard-emission channel should be needed at high $p_T$.","A decisive test of the mechanism would be to use a CCFM-based unintegrated gluon distribution that does include the $k_t>\\mu_F$ region; the paper's logic predicts its $2\\to2$ results should match the JH-2013-plus-$2\\to3$ results without the extra channel.","The suppression parameter $p_T^0=1.5$ GeV is fitted to total charm production; refitting it directly to the $J/\\psi$ spectrum would show whether the regularization is physical or just an effective cutoff."],"forward_implications":["For JH-2013 and similar CCFM-based unintegrated gluon distributions, a $k_T$-factorization description of prompt $J/\\psi$ must include the $2\\to3$ channel $g^* g^* \\to g c\\bar c$ to reproduce the LHCb high-$p_T$ data.","For KMR-type distributions, the standard $2\\to2$ term already accounts for the relevant hard emissions, and adding the $2\\to3$ term with a cut $k_t<\\mu_F$ leaves the predictions essentially unchanged.","The fitted $c\\bar c\\to J/\\psi$ probability differs between the two unintegrated gluon distributions (0.018 for KMR-CT14lo, 0.0065 for JH-2013-set2), so the small-$p_T$ normalization depends on the gluon distribution.","Using both mechanisms and the direct-to-prompt ratio 0.62, the model describes the LHCb prompt $J/\\psi$ rapidity and transverse-momentum distributions at $\\sqrt{s}=7$ TeV."],"supporting_citations":[{"why":"Defines the improved color evaporation model used to convert $c\\bar c$ pairs into $J/\\psi$ with a fitted probability.","marker":"[5]"},{"why":"Reports the same $k_T$-factorization CEM approach and documents the poor JH-2013 $p_T$ slope that motivates the $2\\to3$ extension.","marker":"[9]"},{"why":"Supplies the analytic off-shell matrix element for $g^* g^* \\to c\\bar c$ used in the $2\\to2$ calculation.","marker":"[10]"},{"why":"Provides the Monte Carlo generator used to compute the tree-level off-shell matrix element for $g^* g^* \\to g c\\bar c$.","marker":"[11]"},{"why":"Defines the KMR unintegrated gluon distributions whose content above the factorization scale makes the $2\\to2$ term sufficient.","marker":"[12]"},{"why":"Defines the JH-2013 CCFM-based unintegrated gluon distributions that lack the hard-emission tail and need the $2\\to3$ term.","marker":"[13]"},{"why":"Supplies the ALICE rapidity data for prompt $J/\\psi$ used in the comparisons.","marker":"[15]"},{"why":"Supplies the LHCb rapidity and transverse-momentum data for prompt $J/\\psi$ used in the comparisons.","marker":"[16]"}],"fun_headline_variants":["Hard gluon emission fixes J/psi pT spectrum","Missing 2->3 process explains J/psi high-pT deficit","Explicit g*g*->g c cbar restores J/psi pT shape","CCFM gluon tail needs 2->3 channel for J/psi","New subprocess resolves J/psi momentum mismatch"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the tree-level $2\\to3$ contribution, regularized by $F_{\\rm sup}(p_T)=p_T^4/((p_T^0)^2+p_T^2)^2$ with $p_T^0=1.5$ GeV tuned to total charm production at the LHC, correctly represents the higher-order physics missing from CCFM-type UGDFs; if that suppression or the tree-level normalization is wrong, the claimed agreement is not robust.","fun_headline_variants_meta":{"raw":{"variants":["Hard gluon emission fixes J/psi pT spectrum","Missing 2->3 process explains J/psi high-pT deficit","Explicit g*g*->g c cbar restores J/psi pT shape","CCFM gluon tail needs 2->3 channel for J/psi","New subprocess resolves J/psi momentum mismatch"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000285,"raw_usage":{"total_tokens":1667,"prompt_tokens":922,"completion_tokens":745,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":538,"completion_tokens_details":{"reasoning_tokens":652}},"tokens_in":538,"tokens_out":745,"duration_ms":7221,"temperature":1.0,"reasoning_tokens":652,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:18:58.105386+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the JH-2013 plus $2\\to3$ prediction at $\\sqrt{s}=13$ TeV, or with $p_T^0$ varied between 1.0 and 2.0 GeV, and compare with LHCb transverse-momentum data; if the good 7 TeV agreement depends strongly on this tuned parameter, the regularization is carrying the result rather than the $2\\to3$ mechanism.","supporting_citations":[{"cited_title":"Ma and R","cited_arxiv_id":null,"evidence_quote":"Defines the improved color evaporation model used to convert $c\\bar c$ pairs into $J/\\psi$ with a fitted probability."},{"cited_title":"Catani, M","cited_arxiv_id":null,"evidence_quote":"Supplies the analytic off-shell matrix element for $g^* g^* \\to c\\bar c$ used in the $2\\to2$ calculation."},{"cited_title":"van Hameren, Comput","cited_arxiv_id":null,"evidence_quote":"Provides the Monte Carlo generator used to compute the tree-level off-shell matrix element for $g^* g^* \\to g c\\bar c$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the KMR unintegrated gluon distributions whose content above the factorization scale makes the $2\\to2$ term sufficient."},{"cited_title":"Hautmann and H","cited_arxiv_id":null,"evidence_quote":"Defines the JH-2013 CCFM-based unintegrated gluon distributions that lack the hard-emission tail and need the $2\\to3$ term."},{"cited_title":"Abelev et al","cited_arxiv_id":null,"evidence_quote":"Supplies the ALICE rapidity data for prompt $J/\\psi$ used in the comparisons."},{"cited_title":"Aaij et al","cited_arxiv_id":null,"evidence_quote":"Supplies the LHCb rapidity and transverse-momentum data for prompt $J/\\psi$ used in the comparisons."}],"review_version":1}