Pith. sign in

REVIEW 4 major objections 3 minor 2 cited by

This paper claims that a two-dimensional Bayesian centrality method, built on a product of two Gamma distributions for the joint fluctuations of track-hit count and forward spectator energy, reconstructs the impact-parameter distribution in

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · deepseek-v4-flash

2026-08-04 07:14 UTC pith:Q4MRO2KL

load-bearing objection A two-observable Gamma kernel for centrality at BM@N is a good idea, but the rotation formula is wrong and the validation is in-sample, so the paper needs major revision. the 4 major comments →

arxiv 2510.26276 v2 pith:Q4MRO2KL submitted 2025-10-30 hep-ex

A Two-dimensional Bayesian Approach to Centrality Determination in Nucleus-Nucleus Collisions

D. Idrisov , F. Guber , N. Karpushkin , P. Parfenov This is my paper
classification hep-ex PACS 25.75.-q
keywords centrality determinationimpact parameterBayesian reconstructionGamma distributionheavy-ion collisionsspectator energytrack multiplicitytwo-dimensional observable analysis
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper proposes that centrality in nucleus-nucleus collisions can be determined more reliably by using two correlated observables at once—the number of charged-particle track hits and the energy deposited by spectator nucleons in a forward calorimeter—instead of a single observable. Its central claim is that the joint fluctuation of these two quantities at a fixed impact parameter is well described by the product of two Gamma distributions, obtained after a rotation that is intended to decorrelate the observables, and that this kernel yields an accurate Bayesian reconstruction of the impact-parameter distribution. A sympathetic reader would care because at the low beam energies and small systems considered here, single-observable centrality estimators are skewed in both central and peripheral regions, and a faithful reconstruction of the impact parameter is what makes comparisons with theoretical models meaningful. The quantitative result is that, in a realistic simulation of Xe+CsI collisions at 3.8A GeV, the reconstructed mean impact parameter as a function of centrality agrees with the true generator values within 2%.

Core claim

In the paper's own terms, the discovery is an ansatz for the conditional probability of two detector signals at fixed impact parameter: after a rotation of the track-hit count and spectator energy into new variables with vanishing covariance, each is modeled by a Gamma distribution, so the joint density is the product of two Gammas. The shape and scale parameters are fixed by the impact-parameter-dependent means and variances through k = mean^2/variance and theta = variance/mean, and the whole two-dimensional density is fitted to simulated detector output with a small set of scaling parameters that absorb differences between data and simulation. The fit reproduces the two-dimensional distrib

What carries the argument

Equation (6), P(N_hits, E_spec|b) = Gamma(k1, theta1) * Gamma(k2, theta2), is the load-bearing object: a factorized kernel built from the two rotated observables. The rotation is introduced to decorrelate the two measurements, the means and variances of the rotated variables come from Eqs. (2)–(5), and the Gamma parameters are derived from those means and variances. Bayes' theorem (Eq. (12)) then turns a measured window of the two observables into a posterior probability over the impact parameter, and constrained k-means clustering divides the two-dimensional distribution into centrality classes. The paper also derives a decomposition of multiplicity fluctuations into an ancestor-number comp

Load-bearing premise

The whole construction depends on the assumption that a certain rotation of the two measured signals produces two independent, Gamma-distributed variables, so that their joint probability factorizes exactly as in Eq. (6); if the rotation does not actually eliminate the correlation, or the signals are not Gamma-shaped, the reconstructed impact-parameter distributions inherit a bias.

What would settle it

Take the same simulated events, compute the two rotated variables X1 and X2 using the paper's angle phi = arctan(2Cov/(sigma_E^2 - sigma_N^2)), and evaluate their sample covariance: if it is substantially different from zero, the product-of-Gamma kernel does not describe the true joint distribution. A complementary check is to replace phi by (1/2)arctan(2Cov/(sigma_E^2 - sigma_N^2)) and re-fit; if the reconstructed mean-impact-parameter curve shifts by more than a few percent, the 2% agreement is tied to the rotation convention rather than to the physics.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • If the kernel is right, centrality classes can be defined from two observables simultaneously, reducing the autocorrelation that arises when the same multiplicity used for classification is also used for physics analyses.
  • The method provides a parameterized bridge between simulations and data: the alpha, beta, and epsilon parameters calibrate means, variances, and trigger losses, so the approach can be applied to real detector runs after tuning.
  • The Gamma-based kernel extends reliable centrality reconstruction into peripheral and very central windows where skewness makes a normal kernel inadequate.
  • The 2% agreement of the mean impact parameter versus centrality, if it survives in real data, would permit quantitative comparisons of centrality-dependent observables across models and experiments.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The rotation angle in Section IV is not the standard diagonalizing angle: the covariance of the rotated variables is not actually zero with the stated formula, so the factorization into independent Gammas is an additional model assumption rather than a consequence of the rotation. A corrected angle or a direct numerical check of Cov(X1,X2)=0 would settle this internally.
  • If the product-of-Gamma assumption holds for this pair of observables, the same scheme could in principle be extended to three or more observables, although independence of more than two rotated Gamma variables is unlikely to survive without a more general copula.
  • The validation is done on simulations where the fit function is built from the same generator as the data; a stronger test would apply the fit to an independent event generator or to real data with the calibration parameters left free.
  • A direct, testable signature of the paper's assumption is that conditional slices of the two-dimensional distribution at fixed impact parameter should be Gamma-shaped; event generators that store the true impact parameter could verify this slice-by-slice.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

4 major / 3 minor

Summary. This paper presents a two-dimensional Bayesian centrality-determination method for BM@N Xe+CsI collisions at 3.8A GeV. The authors propose that at fixed impact parameter the joint distribution of track-hit multiplicity and FHCal spectator energy factorizes into a product of two Gamma distributions after a rotation of the observables; the means and variances of the Gamma parameters are obtained from DCM-QGSM-SMM simulation and parameterized by polynomials in centrality. Additional scaling parameters α_E, α_N, β_E, β_N, and an efficiency ε are introduced to connect experiment and simulation. The full two-dimensional distribution is fit to the DCM-QGSM-SMM data, centrality classes are defined via constrained k-means, and the resulting impact-parameter distributions and mean-impact-parameter dependence are compared with the same model, claiming agreement within 2%.

Significance. The paper addresses a relevant problem for low-energy heavy-ion fixed-target experiments, where a non-monotonic FHCal response requires more than one observable for centrality determination. Its strengths include a realistic GEANT4 setup, a complete data-analysis chain, and quantitative comparisons. However, the central mathematical step contains an incorrect rotation angle, the factorized-Gamma kernel is an unverified modeling assumption, and the validation is performed in-sample on the same simulation that provides the parameterization. As a result, the paper's main quantitative claim is not currently supported. The methodological core appears repairable, but the present draft requires substantial correction and additional validation.

major comments (4)
  1. [§IV, Eqs. (3)–(6)] The rotation angle used to diagonalize the covariance matrix is incorrect. Diagonalization requires tan(2φ)=2Cov/(σ_E^2−σ_N^2), i.e. φ=(1/2)arctan(2Cov/(σ_E^2−σ_N^2)). The paper defines φ=arctan(2Cov/(σ_E^2−σ_N^2)). With this angle, Cov(X1,X2)=0.5·sin(2φ)(σ_N^2−σ_E^2)+cos(2φ)Cov, which is not generally zero. For example, with σ_N^2=1, σ_E^2=6, and Cov=1, φ≈0.380 rad and Cov(X1,X2)≈−1. Therefore Eq. (6) does not follow from the stated construction, and this kernel feeds directly into Eq. (8) and all reconstructions in Section V. The authors must correct the angle or explicitly state the transformation actually used in the numerical implementation.
  2. [§IV, Eq. (6)] Even with the corrected half-angle, zero covariance does not imply independence of the rotated variables, nor does it imply that each rotated variable is Gamma distributed. The product-of-Gamma form in Eq. (6) is a modeling assumption, not a consequence of the covariance diagonalization. Since the entire Bayesian inversion uses this kernel, the assumption must be tested explicitly—for example by comparing the fitted two-dimensional density to the empirical distribution in narrow impact-parameter bins or by a quantitative goodness-of-fit test. The lower panels of Fig. 3 show contour comparisons but no statistical quantification of agreement.
  3. [§V, Eqs. (8)–(11), Fig. 10] The validation is circular. The polynomial coefficients describing ⟨E_spec⟩, ⟨N_hits⟩, their variances, and their covariance are obtained from the same DCM-QGSM-SMM sample later used to validate the reconstruction; the parameters α_E, α_N, β_E, β_N, and ε are also fitted to that same two-dimensional distribution. The text itself states that α_E and α_N are close to 1 and β_E, β_N close to zero because the fit function is based on the same simulation. Thus the <2% agreement in Fig. 10 demonstrates consistency of the parameterization with the sample, not predictive accuracy. To support the claim, the authors need an out-of-sample test (e.g., fit on one subset and validate on another), a closure test with injected parameters, or a comparison with a standard centrality method.
  4. [Appendix, Eqs. (14)–(17)] The derivation leading to the β_N parameter appears internally inconsistent. If one sets n=α_N n_MC, then σ_n^2=α_N^2 σ_{n_MC}^2, so the combination σ_n^2−σ_{n_MC}^2 α_N^2 in Eq. (16) vanishes identically and β_N=0. The text appears to treat σ_n^2 and σ_{n_MC}^2 as independent while simultaneously imposing a linear scaling. This undermines the parameterization in Eqs. (10)–(11). The assumptions behind the scaling should be clarified or the derivation corrected.
minor comments (3)
  1. [§V, Eq. (12)] The integral over N_hits is written with lower and upper limits both equal to N1; it should presumably be from N1 to N2. Please correct this typo.
  2. [§III and figure captions] Several typos appear: 'fitst physics run' (should be 'first'), 'dependenc'/'varianc' in Fig. 4, and 'deposed energy' in Fig. 3. These are presentation issues but should be fixed.
  3. [§IV, terminology] The text calls Eq. (6) a two-dimensional Gamma distribution, but it is a product of two independent Gamma distributions. If a true bivariate Gamma is intended, the relation should be clarified; otherwise the terminology is misleading.

Circularity Check

0 steps flagged

No significant circularity: the Bayesian reconstruction is tested as a Monte Carlo closure test, not by reducing to its fitted inputs.

full rationale

The derivation chain is not circular. Eq. (6) proposes a product-of-Gamma fluctuation kernel; the parameters in Eqs. (2)-(5) are moment-matching expressions, not outputs of the centrality reconstruction. The centrality classes and impact-parameter posteriors in Eqs. (8) and (12) follow from Bayes' theorem, with the model only supplying the conditional distribution and prior. The fact that the mean/var/cov parameterizations and alpha, beta, epsilon are fitted to the same DCM-QGSM-SMM sample is acknowledged in Sec. V ('Since the fit function was based on the distributions of the mean and variances from the same simulation...'), so the Fig. 10 agreement is an in-sample closure check of the approximate kernel, not a hidden prediction. The self-citations (refs. [8], [9]) are contextual and are not load-bearing. The strongest concern in this paper is a mathematical one, not circularity: the rotation angle quoted in Sec. IV, phi = arctan[2Cov/(sigma_E^2 - sigma_N^2)], is not the diagonalizing angle for the 2x2 covariance matrix (the correct angle is half that), so the statement that Cov(X1,X2)=0 and hence the factorization in Eq. (6) follows is not established as written. That is a correctness/derivation defect; it does not make the result equivalent to its inputs by construction.

Axiom & Free-Parameter Ledger

6 free parameters · 5 axioms · 0 invented entities

The method's free content is substantial: four global scaling parameters plus an efficiency and the entire set of polynomial parameterizations of b-dependent means/variances/covariance are fitted to the DCM-QGSM-SMM simulation. The central kernel itself is an unproven statistical ansatz. No new physics entities appear.

free parameters (6)
  • alpha_E = approximately 1 in simulation fit
    Global scaling of mean spectator energy vs model value, Eq. (9), fitted to the 2D distribution.
  • alpha_N = approximately 1 in simulation fit
    Global scaling of mean track-hit number vs model value, Eq. (9), fitted to the 2D distribution.
  • beta_E = approximately 0 in simulation fit
    Additive variance term for spectator energy, Eq. (10), fitted; motivated by per-spectator energy variance.
  • beta_N = approximately 0 in simulation fit
    Additive variance term for track-hit multiplicity, Eq. (11), fitted.
  • epsilon = not quoted
    Event-level efficiency/normalization parameter introduced in Section IV and fitted.
  • Polynomial coefficients f(c_b) for means, variances, covariance = not specified
    The b-dependence of mean E, mean N, variances and covariance is approximated by polynomials in true centrality c_b (Figs. 4-5); coefficients are fitted to the DCM-QGSM-SMM model but not listed.
axioms (5)
  • domain assumption DCM-QGSM-SMM + GEANT4 simulated Xe+CsI events are a faithful representation of the true BM@N response for centrality validation.
    The entire validation and centrality-dependent parameterizations come from this generator; no real data are used, so an incorrect generator response would invalidate the method test.
  • ad hoc to paper At fixed impact parameter, the observables are independent Gamma variables after a covariance-diagonalizing rotation.
    Eq. (6) assumes factorization; the paper explicitly notes there is no explicit bivariate Gamma and asserts the rotation achieves independence, but the stated angle is wrong.
  • ad hoc to paper The experiment can be related to the simulation by four global alpha,beta scaling parameters and one efficiency epsilon.
    Eqs. (9)-(11) assume a simple affine relationship between model and experimental mean/variance; no proof that this captures detector calibration and trigger effects.
  • standard math Bayes theorem and the Glauber impact-parameter distribution P(b) = 2*pi*b*P_inel(b)/sigma_inel.
    Standard setup inherited from Refs. [13-16].
  • domain assumption Ancestor-number variance decomposition of Appendix VII applies to spectator energy and multiplicity fluctuations.
    Eqs. (14)-(17) rely on a Glauber ancestor fluctuation model; used to justify the beta parameterization.

pith-pipeline@v1.3.0-alltime-deepseek · 8013 in / 14855 out tokens · 140977 ms · 2026-08-04T07:14:05.833981+00:00 · methodology

0 comments
read the original abstract

The determination of centrality in nucleus-nucleus collisions is a crucial task, as it enables the estimation of the impact parameter and thereby allows for the comparison of experimental results with predictions from theoretical models and other experiments. In this work, we present a two-dimensional approach for centrality determination based on a Bayesian framework. The observables used were the number of track hits and the deposited energy of spectators in the forward hadronic calorimeter. A distribution is proposed to describe the fluctuations of these two observables at a fixed impact parameter value. This distribution provides a more precise description of the observable distributions in both central and peripheral regions. The effectiveness of the proposed method was tested within a realistic BM@N simulation framework for Xe+CsI collisions at a beam energy of 3.8A GeV.

Figures

Figures reproduced from arXiv: 2510.26276 by D. Idrisov, F. Guber, N. Karpushkin, P. Parfenov.

Figure 2
Figure 2. Figure 2: Correlation between the deposited spectator energy in FHCal [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 1
Figure 1. Figure 1: Schematic view of the BM@N setup in the 2023 Xe run [ [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: Fit results (red curves) of the track hit number distribution (upper panels, a-c) and the deposed energy in the FHCal (middle panels, [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figure 5
Figure 5. Figure 5: The left plot shows distribution of the mean values of the [PITH_FULL_IMAGE:figures/full_fig_p004_5.png] view at source ↗
Figure 4
Figure 4. Figure 4: The left plot shows dependenc of the varianc of the energy of [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗
Figure 7
Figure 7. Figure 7: The left plot shows the two-dimensional fit results projected [PITH_FULL_IMAGE:figures/full_fig_p005_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Centrality classes for the two-dimensional distribution. [PITH_FULL_IMAGE:figures/full_fig_p005_8.png] view at source ↗
Figure 6
Figure 6. Figure 6: Results of fitting the two-dimensional distribution of the [PITH_FULL_IMAGE:figures/full_fig_p005_6.png] view at source ↗
Figure 9
Figure 9. Figure 9: Distributions of the impact parameter for three centrality [PITH_FULL_IMAGE:figures/full_fig_p006_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Dependence of the mean impact parameter on centrality. [PITH_FULL_IMAGE:figures/full_fig_p006_10.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Extracting the speed of sound of QCD from transverse momentum fluctuations

    hep-ph 2026-03 conditional novelty 7.0

    From ATLAS data on transverse-momentum fluctuations in ultra-central Pb+Pb collisions, the QGP speed of sound is extracted as 0.496 ± 0.008 at 221 ± 13 MeV, matching lattice QCD.

  2. Methods for Centrality Determination Using Forward Detectors in the BM@N Experiment

    hep-ex 2026-05 unverdicted novelty 5.0

    Modified Bayesian and forward-detector methods for centrality in BM@N Xe+CsI collisions at 3.8 A GeV agree with Glauber model within 5%.

Reference graph

Works this paper leans on

22 extracted references · 1 canonical work pages · cited by 2 Pith papers · 1 internal anchor

  1. [1]

    M. Kapishin,Heavy Ion BM@N and MPD Experiments at NICA, inProceedings of 13th International Conference on Nucleus-Nucleus Collisions, (Saitama, Japan), Journal of the Physical Society of Japan, July, 2020, DOI

    2020

  2. [2]

    P. Senger and BM@N collaboration,The heavy-ion program at the upgraded Baryonic Matter@Nuclotron Experiment at NICA, inProceedings of The International conference on Critical Point and Onset of Deconfinement — PoS(CPOD2021), (zoom), p. 033, Sissa Medialab, Aug., 2022, DOI

    2022

  3. [3]

    Broniowski and W

    W. Broniowski and W. Florkowski,Geometric relation between centrality and the impact parameter in relativistic heavy-ion collisions,Physical Review C65(2002) 024905

    2002

  4. [4]

    Abelev, J

    B. Abelev, J. Adam, D. Adamov ´a, A.M. Adare, M.M. Aggarwal, G. Aglieri Rinella et al.,Centrality determination of Pb-Pb collisions at s N N = 2.76 TeV with ALICE,Physical Review C88(2013) 044909

    2013

  5. [5]

    Lacey, A

    R.A. Lacey, A. Taranenko, J. Jia, D. Reynolds, N.N. Ajitanand, J.M. Alexander et al.,Beam energy dependence of the viscous damping of anisotropic flow,

  6. [6]

    Kagamaster, R

    S. Kagamaster, R. Reed and M. Lisa,Centrality determination with a forward detector in the RHIC Beam Energy Scan, Physical Review C103(2021) 044902

    2021

  7. [7]

    Loizides, J

    C. Loizides, J. Nagle and P. Steinberg,Improved version of the PHOBOS Glauber Monte Carlo,SoftwareX1-2(2015) 13

    2015

  8. [8]

    Parfenov, I

    P. Parfenov, I. Segal, D. Idrisov, V . Luong, A. Taranenko, A. Demanov et al.,Centrality Determination in Heavy-ion Collisions with MPD Detector at NICA,Acta Physica Polonica B Proceedings Supplement14(2021) 503

    2021

  9. [9]

    Segal,Centrality Determination in Heavy-Ion Collisions Based on Monte-Carlo Sampling of Spectator Fragments, Particles6(2023) 568

    I. Segal,Centrality Determination in Heavy-Ion Collisions Based on Monte-Carlo Sampling of Spectator Fragments, Particles6(2023) 568

    2023

  10. [10]

    collaboration, J

    H. collaboration, J. Adamczewski-Musch, O. Arnold, C. Behnke, A. Belounnas, A. Belyaev et al.,Centrality determination of Au+Au collisions at 1.23A GeV with HADES,The European Physical Journal A54(2018) 85

    2018

  11. [11]

    STAR COLLABORATIONcollaboration,Measurements of proton high-order cumulants in √sN N= 3 GeV Au + Au collisions and implications for the qcd critical point,Phys. Rev. Lett.128(2022) 202303

    2022

  12. [12]

    Miller, K

    M.L. Miller, K. Reygers, S.J. Sanders and P. Steinberg, Glauber modeling in high-energy nuclear collisions,Annual Review of Nuclear and Particle Science57(2007) 205

    2007

  13. [13]

    S.J. Das, G. Giacalone, P.-A. Monard and J.-Y . Ollitrault, Relating centrality to impact parameter in nucleus-nucleus collisions,Physical Review C97(2018) 014905

    2018

  14. [14]

    Rogly, G

    R. Rogly, G. Giacalone and J.-Y . Ollitrault,Reconstructing the impact parameter of proton-nucleus and nucleus-nucleus collisions,Physical Review C98(2018) 024902

    2018

  15. [15]

    Frankland, D

    INDRA Collaboration, J.D. Frankland, D. Gruyer, E. Bonnet, B. Borderie, R. Bougault et al.,Model independent reconstruction of impact parameter distributions for intermediate energy heavy ion collisions, 2020. 10.48550/ARXIV .2011.04496. 8

    doi:10.48550/arxiv.2011.04496 2020

  16. [16]

    X. Chen, L. Li, Y . Cui, J. Yang, Z. Li and Y . Zhang,Bayesian reconstruction of impact parameter distributions from two observables for intermediate energy heavy ion collisions, Sept., 2023. 10.48550/arXiv.2201.12586

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2201.12586 2023

  17. [17]

    with the impact parameter. The observed difference between these correlations is at- tributed to the geometry of the BM@N experiment and the positioning of the FHCal in the fitst physics run Xe+CsI at 3.8 AGeV . Due to this configuration, small parts of Xe ion beam and the nuclear fragments from peripheral heavy-ion collisions partially strike the rare pa...

  18. [18]

    Afanasiev et al.,The BM@N spectrometer at the NICA accelerator complex,Nucl

    S. Afanasiev et al.,The BM@N spectrometer at the NICA accelerator complex,Nucl. Instrum. Meth. A1065(2024) 169532 [2312.17573]

    arXiv 2024

  19. [19]

    Tashkandy, M

    Y . Tashkandy, M. Omair and A. Alzaid,Bivariate and bilateral gamma distributions,International Journal of Statistics and Probability7(2018) 66

    2018

  20. [20]

    Baznat, A

    M. Baznat, A. Botvina, G. Musulmanbekov, V . Toneev and V . Zhezher,Monte-Carlo Generator of Heavy Ion Collisions DCM-SMM,Physics of Particles and Nuclei Letters17(2020) 303

    2020

  21. [21]

    Agostinelli, J

    S. Agostinelli, J. Allison, K. Amako, J. Apostolakis, H. Araujo, P. Arce et al.,Geant4—a simulation toolkit, Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment506(2003) 250

    2003

  22. [22]

    Bradley, K

    P. Bradley, K. Bennett and A. Demiriz,Constrained k-means clustering,