{"id":"40c3523c-33e6-4ff2-ae20-12739915bcb5","arxiv_id":"2502.03769","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":2,"one_line_summary":"First measurement of Upsilon meson production cross sections at sqrt(s)=500 GeV in p+p collisions, including multiplicity dependence.","lead":"STAR reports the first measurements of Upsilon(1S), Upsilon(2S), and Upsilon(3S) production in proton-proton collisions at 500 GeV, giving cross sections, ratios, and the dependence of yields on charged particle multiplicity. The results constrain competing quarkonium production models at an energy between RHIC and the LHC.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Upsilon decay electrons are included in Nch (Sec. III A), shifting the numerator multiplicity by +2 relative to the MB denominator; this autocorrelation can mimic a rising N_Y/ trend and is not tested by the 4Cx-tune systematic.","rationale":"The reader's weakest assumption was the model dependence of the Bayesian unfolding. That concern is reasonable but is partially covered by the 4Cx-tune and NBD variations and by the large systematic uncertainties assigned to the high-multiplicity bin. The concern raised here is distinct and more direct: it is a model-independent shift of the Upsilon-event multiplicity distribution relative to the minimum-bias reference, caused by the deliberate inclusion of the two signal-decay electrons in Nch. It is not listed among the sources of systematic uncertainty in Section III E and cannot be constrained by varying PYTHIA8 tunes because all tunes share the same track-inclusion definition. If the standard normalization by minimum-bias or triggered event counts is used, a flat underlying Upsilon yield would produce a rising N_Y/ roughly proportional to N_MB(Nch-2)/N_MB(Nch), which may be sizable in the tail. This does not undermine the integrated cross sections, per-state pT and y spectra, or ratios, which are the other headline results and appear carefully characterized. The multiplicity-dependent comparison to models, however, is not secure until the test above is performed. The final verdict remains CONDITIONAL, so no change from the reader's verdict is needed; the condition should include this reanalysis or a demonstrated cancellation of the signal-track autocorrelation.","tokens_in":17713,"tokens_out":17791,"duration_ms":205894,"concrete_test":"Recompute the multiplicity dependence with the two Upsilon electron tracks excluded from Nch before the Bayesian unfolding, or equivalently re-bin the Upsilon sample using Nch-2 while keeping the minimum-bias distribution unchanged, then compare the resulting N_Y/ versus Nch/ with Fig. 10. Run the PYTHIA8 STAR Heavy Flavor Tune and the CGC/Saturation model with the same exclusion. If the rising trend, or the low first bin, changes by more than the quoted bin-by-bin systematic uncertainties, the published multiplicity result is biased by signal-track autocorrelation; if the trend is unchanged, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section III A explicitly states that the electrons and positrons from Upsilon decays are included in the Nch calculation, while Section IV C uses the minimum-bias Nch distribution as the multiplicity reference. In the standard self-normalized-yield construction, the Upsilon numerator for a bin B draws events whose underlying multiplicity is roughly B-2, whereas the denominator event distribution is drawn from B. Because the minimum-bias multiplicity distribution falls steeply (mean ~ 8), the ratio N_MB(B-2)/N_MB(B) can be materially larger than 1 in the high-multiplicity tail and below 1 in the first bin, producing an apparent enhancement even if the per-event Upsilon yield were flat in underlying multiplicity. The low-Nch suppression of the first point is exactly the signature of such a shift. The systematics in Section III E (unfolding iterations, NBD shape, 4Cx tune, tracking-efficiency variations) do not test this physics-level autocorrelation: all of these variations retain the decay electrons in the response matrix and in the binning, so the 4Cx tune cannot validate the choice. The statement that the Upsilon reconstruction efficiency is flat in Nch addresses detector efficiency, not this bias. Since the multiplicity trend is one of the three headline results and is compared in Fig. 10(b) with CGC/Saturation and Percolation predictions whose multiplicity definition likely excludes the signal leptons, this is the most load-bearing concern about the central new physics claim.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports STAR measurements of Υ(1S), Υ(2S), and Υ(3S) production in p+p collisions at √s = 500 GeV using an integrated luminosity of 13 pb^-1 from the 2011 run. It presents a combined integrated cross section, per-state differential cross sections in pT and rapidity, cross-section ratios, and a measurement of the Υ yield versus charged-particle multiplicity. The results are compared with CEM, CGC+NRQCD, CSM, PYTHIA8 STAR Heavy Flavor Tune, CGC/Saturation, and String Percolation models. The main novelty is the first measurement of these observables at √s = 500 GeV and the multiplicity dependence of Υ production at RHIC.","tokens_in":17982,"tokens_out":3332,"duration_ms":35643,"significance":"If the cross-section and ratio results are correct, they fill an energy gap between fixed-target and Tevatron/LHC measurements and provide new constraints on quarkonium production models, in particular on the low-pT behavior where CGC+NRQCD overshoots the data and on the rapidity dependence where CSM undershoots. The analysis follows standard practices: Crystal Ball signal shapes fixed from full detector embedding, like-sign combinatorial background subtraction, a correlated-background model constrained by PYTHIA8, and a systematic budget that includes trigger response, polarization, tracking efficiency, and unfolding choices. The multiplicity measurement, if robust, would extend the J/ψ multiplicity studies to the Υ sector at RHIC energies and sharpen tests of saturation, percolation, and multi-parton-interaction scenarios.","major_comments":[{"comment":"Section III A explicitly states that the electrons and positrons from Υ decays are included in the Nch calculation, while Section IV C uses the minimum-bias Nch distribution as the reference. In the self-normalized ratio N_Y/ versus N_ch/, the numerator events have their measured multiplicity shifted upward by about 2 units relative to the denominator events for the same underlying multiplicity. Because the minimum-bias multiplicity distribution falls steeply (mean ≈ 8), this shift artificially suppresses the first Nch bin and enhances the high-multiplicity bins, exactly the pattern visible in Fig. 10(a). The systematic checks in Section III E (unfolding iterations, NBD shape, 4Cx tune, tracking-efficiency variations) all retain the decay leptons in the response matrix and binning, so they do not test this physics-level autocorrelation. Please either exclude the two Υ-decay leptons from the numerator multiplicity, or quantify the expected shift using the measured minimum-bias Nch distribution and apply it as a correction or an explicit systematic; without this, the multiplicity trend in Fig. 10 cannot be interpreted as a genuine enhancement of Υ production in high-multiplicity events.","section":"III A, IV C"},{"comment":"The Bayesian unfolding response matrix for the multiplicity measurement is generated with PYTHIA8 using the STAR Heavy Flavor Tune, and the same tune is later shown as a model comparison in Fig. 10(b). The 4Cx tune check in Section III E varies only one neighboring configuration within the same event-generator family; it does not validate the use of a PYTHIA8-based response matrix for the Nch dependence. This self-referentiality does not affect the cross-section results, but it weakens the abstract's statement that the multiplicity trend is 'consistent with ... PYTHIA8'. The paper should explicitly state that the model comparison shares the same simulation framework as the correction, and ideally test the unfolding with a non-PYTHIA8 response or with a larger model variation.","section":"III D and IV C, Fig. 10(b)"}],"minor_comments":[{"comment":"In the sentence describing the minimum-bias Nch distribution, 'criteia' should be 'criteria'.","section":"III A"},{"comment":"Reference [31] lists the author as 'A. Angelis and ohers'; this should read 'others'.","section":"References"},{"comment":"The colors for the CGC+NRQCD bands in the Fig. 5 captions are inconsistent with the text: Fig. 5(b) is described as a brown band in the caption but light blue in the text, and Fig. 5(c) is reversed; these should be made consistent.","section":"Fig. 5"},{"comment":"The same color inconsistency for the CGC+NRQCD bands appears in the Fig. 6 captions and should be corrected.","section":"Fig. 6"},{"comment":"The axis labels in Fig. 10 appear garbled in the manuscript (e.g., the x-axis shows a trailing '>' symbol); the rendered labels should be checked so that the reader can read 'N_ch/' and 'N_Y/' unambiguously.","section":"Fig. 10"}],"recommendation":"major_revision","confidential_remarks":"The cross-section and ratio parts of the paper are well executed and likely to be accepted after the usual polishing. The load-bearing issue is the multiplicity-dependence claim highlighted in the abstract: the inclusion of the signal leptons in the multiplicity estimator can produce the observed rising trend without any physics. Asking the authors to exclude those leptons or to quantify the shift with the measured Nch distribution is a feasible, well-scoped revision; if they cannot do so, the multiplicity part of the paper should be downgraded to a detector-level observation. The manuscript is otherwise within the journal's scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: the cross-section part of this paper is a solid, useful first measurement at sqrt(s)=500 GeV, and the systematics are handled with care. The multiplicity-dependence part has a real self-correlation issue that the paper does not address, so I would not trust that trend as published.\n\nWhat is new: first per-state Upsilon pT and y differential cross sections at 500 GeV, an integrated cross section of 199 +/- 13 +/- 33 pb, cross-section ratios versus energy and multiplicity, and N_Upsilon/ versus N_ch/. The measurement fills a genuine gap between fixed-target and Tevatron/LHC energies. The analysis is standard STAR: Crystal Ball shapes fixed from full simulation, like-sign background subtraction, embedding-based efficiencies, and a long list of correlated and uncorrelated systematics. The CEM description of Upsilon(1S) is a meaningful constraint, and the CGC+NRQCD overestimate at low pT and CSM underestimate of the rapidity dependence are worth reporting. The pT-smearing parameter being tuned to the J/psi width is a minor worry, not a flaw.\n\nSoft spots: first, the differential data are only in figures, with no tables or HEPData entry, which limits the paper's value as a reference. Second, the multiplicity unfolding relies on a PYTHIA8 response matrix checked against one alternate tune; that is a limited but honest test. The bigger issue is the autocorrelation. Section III A explicitly states that the Upsilon decay electrons are included in Nch, while the denominator uses minimum-bias multiplicity. Binning Upsilon events by measured Nch therefore shifts their effective underlying multiplicity by about +2 relative to the minimum-bias reference. Since the multiplicity distribution is steeply falling near ~ 8, a flat per-event yield would appear as a rising trend, and the low first point in Figure 10(a) is exactly the expected signature. The 4Cx-tune systematics do not test this: all variants keep the decay leptons in the response matrix and in the binning. The statement that reconstruction efficiency is flat in Nch addresses detector efficiency, not this physics-level offset. So the comparisons with CGC/Saturation, Percolation, and PYTHIA8 may be mixing different multiplicity definitions.\n\nWho this is for: the quarkonium production community and anyone benchmarking CEM, NRQCD, or CGC. The cross-section results deserve a serious referee; the multiplicity section needs a fix before publication. I would recommend unfolding with Nch excluding the signal leptons, or at least quantifying the shift and showing the trend survives.\n\nRecommendation: send it to peer review, and require the autocorrelation to be addressed and the data tables released.","headline":"The 500 GeV Upsilon cross sections are solid and worth having; the multiplicity-dependence claim has a real, unaddressed autocorrelation from including the Upsilon decay electrons in Nch.","tokens_in":18509,"tokens_out":2511,"would_cite":true,"duration_ms":29759,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper reports the first measurement of Upsilon(1S), Upsilon(2S), and Upsilon(3S) production in proton-proton collisions at 500 GeV, with a combined cross section of 199 ± 13 ± 33 pb in |y|<1 and a multiplicity trend that matches…","keywords":["Upsilon production","bottomonium","quarkonium","p+p collisions","500 GeV","dielectron channel","differential cross sections","charged-particle multiplicity"],"falsifier":"Rebuild the unfolding response matrix from a data-driven embedding or from an event generator with a very different multiparton-interaction model and recompute the multiplicity curves; a shift larger than the quoted tune systematic (up to roughly 13%) would falsify the trend. The absolute cross section could be checked independently by measuring the same 500 GeV p+p system through the dimuon decay channel.","tokens_in":17476,"feed_emoji":"⚛️","tokens_out":11845,"duration_ms":118286,"temperature":0.7,"pith_summary":"Upsilon mesons are bound states of a bottom quark and antiquark, and how they form after the hard b-bbar scattering is still unsettled. This paper fills an energy gap between fixed-target and LHC measurements by reporting the first bottomonium production data in proton-proton collisions at 500 GeV. It gives the combined Upsilon(1S+2S+3S) cross section in $|y|<1$ as $199 \\pm 13 (\\mathrm{stat}) \\pm 33 (\\mathrm{syst})$ pb, plus per-state $p_T$ and rapidity spectra and cross-section ratios. It also measures how the Upsilon yield grows when the collision produces more charged particles. The data favor the Color Evaporation Model for Upsilon(1S) over Color Singlet Model alternatives, show that CGC+NRQCD overpredicts at low $p_T$, and connect quarkonium production to bulk-event multiplicity.","feed_headline":"500 GeV Upsilon data arrive: CEM fits, CGC+NRQCD overshoots","feed_subtitle":"First p+p bottomonium measurement at this energy also shows yield growth with multiplicity that separates model classes.","key_machinery":"The measurement chain is carried by dielectron reconstruction in the time-projection chamber and barrel calorimeter, followed by a simultaneous unbinned-likelihood fit to the invariant-mass spectrum $6.6