Pith. Sign in

REVIEW

Comparison of pair-source and correlation-function based L\'evy analyses in Ar+Sc system at SPS energies

High Energy Physics - Phenomenology · Submitted 6 Oct 2026 · 36 references

Pith's verdict: A useful but unvalidated pipeline comparison: the low-mT method differences are plausible and internally consistent, but the C(q) branch needs an analytic closure test before the offsets are treated as source properties.

T0 · AI review

How this review was made: An AI referee (deepseek/deepseek-v4-flash) read the paper on October 7, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 4 major objections · 4 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv ID2610.08141 (version 1)
Pith IDUMP2WG4X, Pith's permanent record of this paper and its review
CitationBarnab\'as P\'orfy. "Comparison of pair-source and correlation-function based L\'evy analyses in Ar+Sc system at SPS energies." arXiv:2610.08141 (2026). https://pith.science/paper/UMP2WG4X
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

Which fit pipeline you use shifts Lévy radii at low mass

The reading

The paper claims that, in UrQMD simulations of central Ar+Sc collisions, the Lévy femtoscopic parameters extracted from the same events depend on which analysis pipeline is used. At low transverse mass, fitting the FFT-reconstructed three-dimensional correlation function $C(q)$ yields systematically smaller values for the Lévy index $\alpha$, the intercept $\lambda^*$, and the longitudinal radius $R_\text{long}$, and larger values for $R_\text{out}$ and $R_\text{side}$, relative to fitting the one-dimensional projections of the pair-source distribution $D(\rho)$. Toward higher transverse mass the two routes converge for all parameters except $R_\text{long}$, and the direction-averaged radius $\bar{R}$ agrees between the two methods throughout. If this is right, experimental correlation-function results carry a method-dependent offset at low transverse mass whenever they are compared with source-based simulation extractions, and the simulated monotonic decrease of $\alpha$ with $m_T$ further signals that the cascade-only UrQMD baseline needs a hydrodynamic stage to match the measured trend.

What carries the argument

The load-bearing object is the three-dimensional Lévy stable law of Eq. (3), which describes the pair-source distribution $D(\rho)$ and, by Fourier transformation, yields the correlation function $C(q) = 1 + \lambda\, e^{-|q^T R^2 q|^{\alpha/2}}$ of Eq. (4). The paper feeds the same simulated freeze-out coordinates through two fitting pipelines defined in Section 2: one fits one-dimensional projections of $D(\rho)$ in the Bertsch–Pratt out-side-long frame, and the other reconstructs the full three-dimensional $C(q)$ with an odd-bin histogram and a fast Fourier transform, then fits it in a spherical momentum range whose $q_\text{max}$ values are listed in Table 2. All reported parameter differences are differences between these two numerical routes applied to the same events.

Extended reading notes

Core claim

The paper establishes that the two standard routes to Lévy parameters from simulated pion sources do not give identical answers in the Ar+Sc system at SPS energies. Using the same UrQMD events and kinematic cuts, the pair-source pipeline fits one-dimensional projections of the freeze-out distance distribution $D(\rho)$ with a three-dimensional Lévy stable law, while the correlation-function pipeline constructs the full three-dimensional $C(q)$ by fast Fourier transform of $D(\rho)$ and fits Eq. (4) over a spherical range $|\vec{q}| < q_\text{max}$. At low $m_T$, $C(q)$ returns systematically smaller $\alpha$, $\lambda^*$, and $R_\text{long}$ and larger $R_\text{out}$ and $R_\text{side}$ than $D(\rho)$; at high $m_T$ the sets converge apart from $R_\text{long}$, and $\bar{R}$ agrees between methods. The author presents these differences as a practical caution for comparing simulated pair-source studies with the correlation functions that experiments can actually measure, and notes that the monotonic $m_T$ dependence of $\alpha$ in UrQMD contradicts the non-monotonic trend in NA61/SHINE data.

Load-bearing premise

The load-bearing premise is that the FFT-reconstructed three-dimensional correlation function, built with the odd-bin histogram construction of Section 2.2 and the hand-chosen spherical fit ranges of Table 2, preserves the pair-source information well enough that the fitted differences between the two methods are properties of the simulated source rather than artifacts of binning, normalization, or truncation.

Editorial extensions

If this is right

  • Experimental femtoscopic results at low transverse mass should not be compared directly against source-based simulation extractions without accounting for a method-dependent offset of the size documented here for $\alpha$, $\lambda^*$, and $R_\text{long}$.
  • The direction-averaged radius $\bar{R}$ is robust to the choice of pipeline, so collision-energy and system-size trends in $\bar{R}$ can be quoted with less caveat.
  • The most pipeline-sensitive parameters at low $m_T$ are the Lévy index $\alpha$ and the halo parameter $\lambda^*$, so physical interpretations of their low-mass trends in experimental data will need a systematic cross-check between extraction methods.
  • A cascade-only UrQMD baseline is insufficient for reproducing the measured $m_T$ dependence of $\alpha$, and a hydrodynamic phase should be included before using such baselines to constrain the onset of the quark–gluon plasma.
  • At higher transverse mass the two methods agree except for $R_\text{long}$, meaning that high-$m_T$ comparisons between measured correlation functions and simulated sources can rely on the extracted radii with much less method bias.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A decisive check the paper does not perform is to feed an analytic Lévy source of known $\alpha$ and $R^2$ through the same odd-bin histogram, FFT reconstruction, spherical $q_\text{max}$ cuts, and Minuit2 fits; if the low-$m_T$ offsets persist on that exactly known input, they are numerical artifacts of the two pipelines rather than properties of the UrQMD source.
  • The convergence of the two methods at high $m_T$ hints that the offset tracks the size of $q_\text{max}$ relative to the inverse source size, which is a testable prediction: re-running the comparison with matched dimensionless fit ranges should move the two parameter sets together at all $m_T$.
  • Because $\bar{R}$ is pipeline-independent while the individual radii are not, a practical editorial recommendation would be to quote $\bar{R}$ in low-$m_T$ model-data comparisons and to treat the directional split between $R_\text{out}$, $R_\text{side}$, and $R_\text{long}$ at those masses as method-sensitive.
  • The monotonic-versus-non-monotonic $\alpha(m_T)$ discrepancy is likely tied to the absence of a hydrodynamic phase in UrQMD, and a natural follow-up is to re-run both pipelines on the same collision system with a hydrodynamically coupled generator to see whether the two methods also converge differently there.
Original abstract

Non-Gaussian features in femtoscopic observables have been observed in high-energy collisions, yet their origin remains not fully understood. This study focuses on the phenomenological approach to femtoscopy, going beyond the Gaussian approximation by utilizing L\'evy-stable source distributions. Using UrQMD simulations, three-dimensional analyses of pion sources are done in central Ar+Sc collisions, and the extracted parameters are studied as functions of transverse mass and collision energy. This analysis contributes to a hadronic baseline motivated by past and future research work in femtoscopy at the NA61/SHINE experiment. A direct methodological comparison is presented between parameters extracted from the projections of three-dimensional pair-source distribution and those obtained from the three-dimensional correlation function, reconstructed via Fast Fourier Transform.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-10-07: 1 request, 4,302 input tokens.

What the abstract itself states, as probabilities: claims first 0.03 · states limitation 0.19 · human participants 0.01 · controlled comparison 0.26 · randomized assignment 0.01 · states quantitative result 0.05 · compares with prior results 0.92 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

Non-Gaussian features in femtoscopic observables have been observed in high-energy collisions, yet their origin remains not fully understood. This study focuses on the phenomenological approach to femtoscopy, going beyond the Gaussian approximation by utilizing L\'evy-stable source distributions. Using UrQMD simulations, three-dimensional analyses of pion sources are done in central Ar+Sc collisions, and the extracted parameters are studied as functions of transverse mass and collision energy. This analysis contributes to a hadronic baseline motivated by past and future research work in femtoscopy at the NA61/SHINE experiment. A direct methodological comparison is presented between parameters extracted from the projections of three-dimensional pair-source distribution and those obtained from the three-dimensional correlation function, reconstructed via Fast Fourier Transform.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 Non-Gaussian features in femtoscopic observables have been observed in high-energy collisions, yet their origin remains not fully understood. background 0.65 0.07 0.03 asserted 0.99
s1 This study focuses on the phenomenological approach to femtoscopy, going beyond the Gaussian approximation by utilizing L\'evy-stable source distributions. contribution 0.72 0.56 0.04 asserted 0.97
s2 Using UrQMD simulations, three-dimensional analyses of pion sources are done in central Ar+Sc collisions, and the extracted parameters are studied as functions of transverse mass and collision energy. method 0.91 0.94 0.06 asserted 0.99
s3 This analysis contributes to a hadronic baseline motivated by past and future research work in femtoscopy at the NA61/SHINE experiment. contribution 0.63 0.29 0.04 asserted 0.98
s4 A direct methodological comparison is presented between parameters extracted from the projections of three-dimensional pair-source distribution and those obtained from the three-dimensional correlation function, reconstructed via Fast Fourier Transform. contribution 0.50 0.92 0.10 asserted 1.00

What kind of paper the classifier reports this is: simulation study 0.99 · evidence from simulation or computation 1.00.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. One-dimensional projections (“out”, “side”, “long”) of the three
Figure 2. The L´evy index α for 0–10% central collisions in Ar+Sc as a function of mT. The dashed lines represent values extracted from the simulated D(ρ) distributions, while the solid lines represent the C(q) results. The colored bands denote systematic uncertainties. with increasing mT. Finally, as a way of comparing three-dimensional results to one-dimensional experimental data of NA61/SHINE [35], the average of the three directions as R¯ = q (R2 out + R2 side + R2 long)/3 is also checked. As shown in
Figure 3. The L´evy scale parameter Rout for 0–10% central collisions in Ar+Sc as a function of mT. The dashed lines represent values extracted from the simulated D(ρ) distributions, while the solid lines represent the C(q) results. The colored bands denote systematic uncertainties. 0.2 0.4 0.6 (GeV) mT 2 4 6 8 10 (fm) side R 150A GeV/c 0.2 0.4 0.6 (GeV) mT 2 4 6 8 10 (fm) side R 30A GeV/c C(q) D(ρ) PLB877(2026)140446 0.2 0.4 0.6 (GeV) mT 2 4 6 8 10 (fm) side R 75A GeV/c 0.2 0.4 0.6 (GeV) mT 2 4 6 8 10 (fm) side R 19A GeV/c 0.2 0.4 0.6 (GeV) mT 2 4 6 8 10 (fm) side R 40A GeV/c 0.2 0.4 0.6…
Figure 4. The L´evy scale parameter Rside for 0–10% central collisions in Ar+Sc as a function of mT. The dashed lines represent values extracted from the simulated D(ρ) distributions, while the solid lines represent the C(q) results. The colored bands denote systematic uncertainties. 7
Figure 5. The L´evy scale parameter Rlong for 0–10% central collisions in Ar+Sc as a function of mT. The dashed lines represent values extracted from the simulated D(ρ) distributions, while the solid lines represent the C(q) results. The colored bands denote systematic uncertainties. 0.2 0.4 0.6 (GeV) mT 2 4 6 8 10 (fm) R 150A GeV/c 0.2 0.4 0.6 (GeV) mT 2 4 6 8 10 (fm) R 30A GeV/c C(q) D(ρ) PLB877(2026)140446 0.2 0.4 0.6 (GeV) mT 2 4 6 8 10 (fm) R 75A GeV/c 0.2 0.4 0.6 (GeV) mT 2 4 6 8 10 (fm) R 19A GeV/c 0.2 0.4 0.6 (GeV) mT 2 4 6 8 10 (fm) R 40A GeV/c 0.2 0.4 0.6 (GeV) mT 2 4 6 8 10 (fm) R…
Figure 6. The L´evy scale parameter R¯ = q (R2 out + R2 side + R2 long)/3 for 0–10% central collisions in Ar+Sc as a function of mT. The dashed lines represent values extracted from the simulated D(ρ) distributions, while the solid lines represent the C(q) results. The colored bands denote systematic uncertainties. 8
Figure 7. The intercept λ ∗ for 0–10% central collisions in Ar+Sc as a function of mT. The dashed lines represent values extracted from the simulated D(ρ) distributions, while the solid lines represent the C(q) results. The colored bands denote systematic uncertainties. simulated events. The key goal of this work was to illustrate the difference between using the pair-source distribution or correlation function to estimate the L´evy fit parameters. Parameters extracted directly from D(ρ) fits are compared with those ob￾tained from FFT-reconstructed correlation functions, C(q), in the Ar+Sc sys

Figures reproduced from arXiv:2610.08141 by the authors.

References

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.