{"id":"b54e7461-8d17-4b11-bd20-8134b95d99d0","arxiv_id":"1908.05465","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For mixed-event background estimation of conserved-charge cumulants, sampling each mixed event from a global pool of all particles removes the influence of the original charge distribution even at low statistics, unlike the conventional per-event method.","lead":"This paper compares four ways to build mixed events used as backgrounds in heavy-ion collision cumulant measurements. It finds that drawing each mixed event's particles at random from one global pool is the most robust scheme.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Method-III's superiority is not established: the paper equates 'best' with erasing the original Nc distribution, but event-by-event variation of the charge fraction is itself a non-critical global/systematic effect that a mixed-event background should reproduce.","rationale":"The reader's weakest_assumption identifies the same load-bearing issue, and I agree with it. The paper's technical derivation of Method-III's independence from the original Nc distribution is credible; the problem is normative. The strongest_claim is a recommendation, and the recommendation follows only if the criterion 'least constrain' is the right one for a mixed-event background. The introduction enumerates non-critical global effects that operate through event-level charge fraction, and Method-III erases them. Therefore the paper's conclusion overreaches. A conditional verdict is appropriate: the methods paper can be accepted if it restricts its claim to 'most independent of P_o(Nc)' or supplies a physics argument that charge-fraction variation is a critical correlation. Since the reader already assigned CONDITIONAL, no adjustment is needed.","tokens_in":6667,"tokens_out":8309,"duration_ms":88038,"concrete_test":"Simulate with no intra-event correlations: for each event, draw Nch_i from P(Nch), draw f_i from a non-degenerate distribution (e.g., two values 0.15 and 0.45 with equal probability), and draw Nc_i ~ Binomial(Nch_i, f_i). This model's true non-critical cumulant kappa*sigma^2 is the expectation over f_i of the binomial cumulant. Build mixed samples from these events using Method-I and Method-III, with the same total statistics (1e6 events, with sub-samples from 1e3 to 1e6 events). Compare the mixed-sample kappa*sigma^2 to the true model value. Method-I should match within errors; a statistically significant deviation of Method-III shows that its independence from P_o(Nc) removes a non-critical global effect and that the recommendation 'best' is not universally valid.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Sections III C and IV conclude that Method-III is best because its mixed-sample cumulants are independent of P_o(Nc). This conclusion depends on an unstated criterion: that the full shape of the original Nc distribution is an unwanted correlation. The paper's own taxonomy (Sec. I) lists centrality bin width and initial size fluctuations as non-critical global/systematic effects. In realistic data, those effects vary the per-event charge fraction f_i = Nc_i/Nch_i, not just the multiplicity. Method-I/II draw each mixed-event particle from a randomly chosen original event, so the mixed-event probability of a charged particle is drawn from the distribution of f_i; the resulting Nc|Nch distribution is a mixture of binomials and retains the variance of f_i. Method-III draws from a global pool with fixed p_c = sum Nc_i / sum Nch_i, giving a pure binomial Nc|Nch; it removes the variance of f_i entirely. If f_i fluctuations are non-critical, Method-III underestimates the background and would leave those non-critical fluctuations in any subtracted cumulant. The paper does not argue that f_i variation is a critical correlation; indeed its stated goal of preserving global/systematic characters points the other way. Thus the ranking is a choice of criterion, not a derived result. The analytic baseline in Fig. 1 is also asserted without derivation, leaving the toy-model calibration unverified.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper considers four methods for constructing mixed-event samples to serve as backgrounds for higher cumulants of conserved charges in heavy-ion collisions. Method-I selects, for each mixed event, Nch original events and takes one particle from each; Method-II removes used particles from the original sample; Method-III pools all particles from all original events and draws Nch particles from the pool with replacement; Method-IV draws without replacement. The authors derive expressions for the mixed-event Nc distributions, and test the methods on toy models with fixed multiplicity and either Poisson or two-delta Nc distributions. They find that all four methods preserve the mean charged-particle number and remove the correlation between Nc and the event label. They then show that the kappa-sigma-squared cumulant of the mixed sample for Method-III is independent of the original Nc distribution even at low statistics, while Method-I only approaches the analytic expectation at high statistics. The paper concludes that Method-III, being the 'most random or least constrained', is the best method, and recommends it over the conventional Method-I.","tokens_in":6953,"tokens_out":8594,"duration_ms":75997,"significance":"The manuscript addresses a practical question of direct relevance to the STAR/RHIC cumulant program: how to construct mixed-event backgrounds without introducing statistical bias at limited event counts. The paper's explicit formulation of four methods and its transparent toy-model comparisons are useful, and the derivation of the probability distributions, despite some technical ambiguities, correctly identifies a real difference between the methods: Method-III eliminates the event-by-event variation of the charged fraction and thus converges much faster. If the recommended criterion—that the background should be independent of the original Nc distribution—is accepted, the conclusion is well-supported by the simulations. However, the paper does not justify that criterion, and it is in tension with the usual purpose of mixed-event backgrounds, which is to reproduce non-critical global/systematic effects such as centrality-bin-width and volume fluctuations. The analytic baseline in Fig. 1 is not derived, and Methods-II and -IV are not shown.","major_comments":[{"comment":"The conclusion that Method-III is 'best' rests on the assumption that the influence of the original distribution of Nc on mixed-event cumulants is a defect to be removed. This assumption is not stated or defended. In Sec. I, the paper itself identifies centrality bin width and initial size fluctuations as non-critical global/systematic effects that the mixed-event background should help subtract; those effects produce event-by-event fluctuations in the charged fraction f_i = Nc_i/Nch_i, not just in Nch. Method-I/II preserve these fluctuations, because the mixed-event Nc|Nch distribution is a mixture of binomials with event-specific probabilities f_i, whereas Method-III forces a single global probability p_c and therefore removes the variance of f_i entirely. If f_i fluctuations are non-critical, Method-III would underestimate the background, leaving those non-critical fluctuations in the subtracted cumulants. The paper provides no physics argument that the full shape of P^o(Nc) beyond its mean is a critical correlation to be erased. The ranking of methods is therefore a choice of criterion rather than a derived result.","section":"Sec. III C, Sec. IV, Abstract"},{"comment":"The blue horizontal line labeled 'expectation of analytic calculation' is introduced without derivation. Because the central comparison uses this line as the reference value to which Method-I should converge, the paper should provide the explicit analytic expression for the asymptotic kappa-sigma-squared of the mixed sample in the fixed-multiplicity toy model, and show how it depends on the original Nc distribution for Methods I and III. Without this, the claim that Method-I is 'influenced by the original distribution of Nc' at low statistics is only a simulation observation, and the reader cannot assess the size of the bias or the convergence rate claimed.","section":"Fig. 1, Sec. III C"},{"comment":"The text states that the statistics dependences for Methods II and IV are 'similar to Fig. 1(a) and (b), respectively' but does not show them. Since the paper's purpose is to compare four methods, the omission prevents the reader from verifying the conclusion that the no-replacement constraint is 'trivial and unnecessary'. A complete comparison, for example as additional panels in Fig. 1 or as a table of the four methods' kappa-sigma-squared values at representative statistics, should be added.","section":"Sec. III C, Fig. 1"},{"comment":"The probability formulas for Methods I and II are not properly defined. In the stated algorithm, an event is selected randomly from the original sample, so the probability of drawing a charged particle from the selected event is simply that event's charged fraction f_i = Nc_i/Nch_i; the factors P^o(Nc) and P^o(Nch) in Eqs. (3)-(4) do not represent the event-selection probability. As written, these formulas appear to use the product of one-dimensional marginal distributions, which is not the joint probability of drawing an event with given (Nc,Nch). The derivation of P^m(Nc) should be corrected or clarified, since it is the basis for the claim that P^o(Nc) enters the mixed-event distribution in Method-I.","section":"Sec. II, Eqs. (2)-(6)"}],"minor_comments":[{"comment":"The phrase 'Method-VI' should be 'Method-IV' in the sentence 'Obviously, Methods-II and Method-VI are more restrictive...'.","section":"Sec. II"},{"comment":"The phrase 'It is showed' should be 'It is shown'.","section":"Abstract"},{"comment":"The notation 'C(N e,N e, c, )' is confusing; please use 'C(N^e, N_c^e)'.","section":"Eq. (10)"},{"comment":"The legend text 'Poisson, δ2, )c(NoP' is garbled; it should read 'Poisson' and 'two-δ functions'.","section":"Fig. 1 caption"},{"comment":"The reported uncertainties should specify whether they are standard deviations of the distribution of the mean or standard errors of the mean.","section":"Tables I and II"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a valid and useful comparison of mixed-event construction methods, but the paper's main conclusion is over-stated relative to what is actually shown. The key issue is the unexamined evaluative criterion (independence from the original Nc distribution), and the missing analytic derivation and complete method comparison. If the authors can reframe the conclusion or justify the criterion with physics arguments, the paper could be suitable for publication. The toy-model results are reproducible in principle, but the presentation needs careful revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper lays out four ways to construct mixed events for cumulant backgrounds and compares them in a simple toy model. The construction equations are internally consistent, and the low-statistics effect they identify is real: Method-I keeps a memory of the original Nc distribution, while Method-III does not. That is a useful observation for the BES community.\n\nThe soft spot is the jump from 'independence from P_o(Nc)' to 'best.' The paper never argues that the full original Nc distribution is an unwanted influence. By its own taxonomy, non-critical global/systematic effects include centrality bin width, initial size fluctuations, and acceptance; these change the per-event charge fraction f_i = Nc/Nch. Method-I/II, by sampling particles from randomly chosen original events, retain the variance of f_i; Method-III, by drawing from a global pool with fixed p_c, removes it entirely. If f_i fluctuations are non-critical, Method-III removes part of the background and leaves those fluctuations in the subtracted cumulants. The paper's own goal of preserving global/systematic characters points the other way. The ranking is a choice of criterion, not a derived result. The stress-test concern lands exactly here.\n\nMinor issues: the analytic baseline in Fig. 1 is asserted without derivation, and Methods-II/IV are stated to behave like I/III without being shown. These are addressable, but the criterion issue is the load-bearing one.\n\nThat said, the comparison framework is transparent, the toy model is reproducible, and the low-statistics behavior of Method-I is genuinely worth knowing. The paper is not broken; it is missing an argument. With a clear statement of what the background must preserve, the criterion question could be resolved. As it stands, the conclusion is conditional on a choice the authors do not defend.\n\nWho should read this: experimentalists doing BES cumulant analyses and anyone building mixed-event estimators. It deserves serious peer review because the issue is substantive and addressable, but I would require a revised version that derives the analytic expectations, shows all four methods, and defends the criterion or modifies the recommendation. I would cite it, with care, as a cautionary comparison rather than a settled prescription.","headline":"Useful comparison of mixed-event schemes, but the 'best' claim rests on an unexamined criterion: erasing the original Nc distribution also removes non-critical global fluctuations that a background should preserve.","tokens_in":7409,"tokens_out":2996,"would_cite":true,"duration_ms":30319,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["25.75.Dw","12.38.Mh","25.75.Nq"],"model":"deepseek-v4-flash","headline":"The paper claims that the conventional mixed-event background for higher cumulants of conserved charges should be replaced by global-pool sampling, which makes the mixed sample's cumulants independent of the original charged-particle…","keywords":["mixed events","higher cumulants","conserved charges","QCD phase transition","non-critical fluctuations","heavy-ion collisions","event mixing methods","mixed-event background"],"falsifier":"Generate events with a known non-critical source of kappa-sigma-squared, for example an event-by-event fluctuation in the charged fraction tied to multiplicity, and compare the non-critical signal recovered by subtracting the conventional and global-pool backgrounds; if the global-pool background does not reproduce the known non-critical kappa-sigma-squared while the conventional one does, the paper's claim that the global-pool method is best fails.","tokens_in":6482,"feed_emoji":"🎲","tokens_out":8769,"duration_ms":79971,"temperature":0.7,"pith_summary":"This paper argues that the standard way of building mixed events, the background samples used to subtract non-critical fluctuations from higher cumulants of conserved charges in heavy-ion collisions, is not the best one. It compares four mixing schemes and claims that drawing every particle of a mixed event at random from a single pool containing all particles of all original events is superior. At low statistics, the conventional method still lets the original distribution of charged-particle numbers leak into the mixed-event cumulants, while the global-pool method does not. The authors therefore recommend the global-pool method as the least-constrained, most-random construction.","feed_headline":"Use one global particle pool for mixed-event cumulant backgrounds","feed_subtitle":"The global-pool method removes the original charged-particle distribution's influence even at low statistics.","key_machinery":"The central object is the mixed-event probability distribution for the number of charged particles. Method-III replaces per-event charged fractions with a single pool fraction, so the mixed distribution is a convolution of the original multiplicity distribution with a binomial distribution whose success probability is the global charged fraction. Because that probability is constant, the original charged-particle distribution cannot enter the mixed-event cumulants at any statistics, which is the property the paper identifies as the criterion for best.","core_discovery":"The paper's central claim is that Method-III, in which all particles of the original events are put into one pool and each mixed event is built by drawing its multiplicity of particles at random from that pool, yields mixed-event cumulants that are independent of the original distribution of charged-particle numbers at all tested statistics, including as few as ten thousand events. In the conventional Method-I, the probability of drawing a charged particle from an event depends on that event's charged fraction, so the original distribution continues to influence the mixed sample's kappa-sigma-squared until statistics are large. Because the pool probability is a single constant, Method-III erases that influence. Methods that forbid reusing original particles are found to be unnecessary constraints, and the authors conclude that the global-pool method is the best mixed-event prescription for cumulant backgrounds.","pith_inferences":["The paper's criterion of best is maximal randomness, meaning independence from the original charged-particle distribution; a complementary criterion would ask whether the background reproduces all non-critical correlations, in which case the global-pool method's erasure of the joint charged-number and multiplicity distribution would be a drawback rather than a virtue.","A directly testable extension is to simulate centrality or volume fluctuations in the charged fraction and compare how faithfully the conventional and global-pool backgrounds recover the known non-critical kappa-sigma-squared, which would show whether global-pool sampling over-subtracts.","Because the global-pool method has a closed analytic form, it could be used to derive a correction converting measured cumulants into estimates of the critical excess, a construction the paper does not pursue."],"forward_implications":["Replacing the conventional per-event selection with global-pool sampling removes the bias that the original charged-particle distribution imprints on mixed-event cumulants at low statistics.","At the lowest statistics tested, the global-pool method produces mixed-event cumulants that are indistinguishable for very different original charged-particle distributions, so background estimates no longer depend on the shape of that distribution.","The no-reuse constraints of the other two methods are unnecessary, since they add no improvement over their with-reuse counterparts.","A global-pool mixed sample preserves the multiplicity distribution and the mean charged-particle number within statistical error, retaining the global non-critical features the background is designed to estimate."],"supporting_citations":[{"why":"Defines the conventional mixed-event construction in which each particle is taken from a different original event; this is the baseline Method-I that the paper argues against.","marker":"[18]"},{"why":"One of the standard heavy-ion applications of mixed-event samples that motivates the conventional approach being re-examined.","marker":"[24]"},{"why":"Another standard mixed-event usage in heavy-ion data analysis that the conventional method inherits.","marker":"[25]"},{"why":"Supplies the sub-sample averaging procedure used to compute statistical errors on the mixed-event cumulants.","marker":"[15]"},{"why":"Provides the statistical-error method and the Poisson baseline used in the low-statistics test.","marker":"[16]"},{"why":"Establishes statistical fluctuations from finite particle number as a non-critical effect that mixed-event backgrounds are meant to estimate.","marker":"[17]"}],"fun_headline_variants":["Global pool method tops for mixed-event cumulants","One global pool erases distribution bias in cumulants","Use one random pool for all mixed-event cumulants","Simplest mixing: single global pool wins","Random global draws yield clean cumulant backgrounds"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the original distribution of event charged-particle counts is an unwanted influence that a mixed-event background should erase; if fluctuations in that distribution carry physical global effects that the background is supposed to mirror, the global-pool method would subtract real physics and calling it best would not follow.","fun_headline_variants_meta":{"raw":{"variants":["Global pool method tops for mixed-event cumulants","One global pool erases distribution bias in cumulants","Use one random pool for all mixed-event cumulants","Simplest mixing: single global pool wins","Random global draws yield clean cumulant backgrounds"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000272,"raw_usage":{"total_tokens":1534,"prompt_tokens":748,"completion_tokens":786,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":364,"completion_tokens_details":{"reasoning_tokens":715}},"tokens_in":364,"tokens_out":786,"duration_ms":8416,"temperature":1.0,"reasoning_tokens":715,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:12:43.515021+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Generate events with a known non-critical source of kappa-sigma-squared, for example an event-by-event fluctuation in the charged fraction tied to multiplicity, and compare the non-critical signal recovered by subtracting the conventional and global-pool backgrounds; if the global-pool background does not reproduce the known non-critical kappa-sigma-squared while the conventional one does, the paper's claim that the global-pool method is best fails.","supporting_citations":[{"cited_title":"(NA49 Collaboration)","cited_arxiv_id":null,"evidence_quote":"One of the standard heavy-ion applications of mixed-event samples that motivates the conventional approach being re-examined."},{"cited_title":"Phys.Rev.C, 77:064908, (2008)","cited_arxiv_id":null,"evidence_quote":"Another standard mixed-event usage in heavy-ion data analysis that the conventional method inherits."}],"review_version":1}