Comparison of Methods for Extracting Three-Dimensional Lévy Source Parameters in √(sNN) = 200 GeV Au+Au Collisions with EPOS3
Abstract
Bose--Einstein femtoscopy probes the spatiotemporal geometry of the particle-emitting source in heavy-ion collisions, where experimental data favor a Lévy-stable source description. We compare four approaches for extracting 3D Lévy source parameters using π+π+ pairs from EPOS3-simulated Au+Au collisions at √(sNN) = 200 GeV: direct fits to the pair-source distribution D(ρ), and fits to the two-particle correlation function C(q) obtained via (i) the CRAB package, (ii) direct numerical integration of D(ρ), and (iii) fast Fourier transformation (FFT). On an identical momentum grid, numerical integration and FFT are analytically identical and numerically equivalent (RMSE < 10−14), while FFT reduces the computation time from hours to seconds. We extract the transverse-mass dependence of the Lévy radii, Lévy index α, and correlation strength λ in 0−−10 % central collisions. All methods reproduce the same physical trends, but CRAB results fluctuate due to event pooling, and correlation-function fits deviate systematically from direct source fits, e.g., with ∼ 7−−9 % lower α. These deviations may reflect residual non-Lévy features of the source, probed differently in coordinate and momentum space, and the comparison of event-averaged parameters with fits to the ensemble-accumulated source, a mixture of Lévy distributions with lower effective α.
Related
- Inferring Nuclear Radial Parameters from Initial-State Information in Heavy-Ion Collisions Using Deep Neural Networks
- Femtoscopic parameters of pion emission sources in EPOS4 Au+Au collisions from 7.7 to 200 GeV
- Comparison of existing methods to derive phytoplankton cell diameter from flow cytometry and ways for improvements
- Algebras of order parameters in one-dimensional spin systems
- Wide-Angle Multislice Propagation for Electron Microscopy: Comparison of Fresnel, Angular-Spectrum, and Wave-Propagation Methods
BibTeX
@misc{szilagyi2026comparison,
title = {{Comparison of Methods for Extracting Three-Dimensional Lévy Source Parameters in $\sqrts_\textNN = 200\text GeV$ Au+Au Collisions with EPOS3}},
author = {Máté Szilágyi and Máté Csanád and Dániel Kincses},
howpublished = {arXiv},
year = {2026},
month = nov,
doi = {10.48550/arxiv.2611.08069},
eprint = {2611.08069},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2611.08069},
}
RIS
TY - UNPB
TI - Comparison of Methods for Extracting Three-Dimensional Lévy Source Parameters in √(s_NN) = 200 GeV Au+Au Collisions with EPOS3
AU - Máté Szilágyi
AU - Máté Csanád
AU - Dániel Kincses
PY - 2026
DA - 2026/11/05
PB - arXiv
DO - 10.48550/arxiv.2611.08069
UR - https://arxiv.org/abs/2611.08069
AB - Bose--Einstein femtoscopy probes the spatiotemporal geometry of the particle-emitting source in heavy-ion collisions, where experimental data favor a Lévy-stable source description. We compare four approaches for extracting 3D Lévy source parameters using π^+π^+ pairs from EPOS3-simulated Au+Au collisions at √(s_NN) = 200 GeV: direct fits to the pair-source distribution D(ρ), and fits to the two-particle correlation function C(q) obtained via (i) the CRAB package, (ii) direct numerical integration of D(ρ), and (iii) fast Fourier transformation (FFT). On an identical momentum grid, numerical integration and FFT are analytically identical and numerically equivalent (RMSE < 10^−14), while FFT reduces the computation time from hours to seconds. We extract the transverse-mass dependence of the Lévy radii, Lévy index α, and correlation strength λ in 0−−10 % central collisions. All methods reproduce the same physical trends, but CRAB results fluctuate due to event pooling, and correlation-function fits deviate systematically from direct source fits, e.g., with ∼ 7−−9 % lower α. These deviations may reflect residual non-Lévy features of the source, probed differently in coordinate and momentum space, and the comparison of event-averaged parameters with fits to the ensemble-accumulated source, a mixture of Lévy distributions with lower effective α.
ER -