REVIEW 4 major objections 7 minor 1 cited by
Analysis Note: Directed flow $v_1$ of protons in the Xe+Cs(I) collisions at 3.8 AGeV
T0 review · 4 major / 7 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read This paper reports the first measurement of the directed flow $v_1$ of protons in 10-30% central Xe+Cs(I) collisions at 3.8 AGeV, showing that the rapidity dependence is roughly described by the JAM model with a momentum-dependent mean…
desk verdict A genuine first proton v1 point for Xe+Cs at 3.8 AGeV, using solid methods and transparent systematics, but with a self-declared preliminary centrality calibration and no numerical tables; worth refereeing conditionally. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the first-order event plane vector $Q_1$ reconstructed from the energy deposition of spectator fragments in the forward hadron calorimeter (FHCal). The method uses three (and sometimes four) sub-events defined by FHCal pseudorapidity ranges F1, F2, F3 and by charged-track groups $T^+$ and $T^-$, with the resolution correction factor $R_1^y$ computed from the $Y$-component correlations. Because the BM@N magnetic field deflects charged particles along $x$, the directed flow is extracted only from $Y$ components: $v_1 = 2\langle y_1 Y_1^a\rangle / R_1^y$. Non-uniform azimuthal acceptance is handled by recentering, twist, and rescaling corrections from the QnTools framework. The midrapidity slope is obtained by fitting $v_1(y_{cm})$ with $v_1 = a + b y_{cm} + c y_{cm}^3$.
What would settle it
Recompute the centrality boundaries after measuring the CCT2 trigger efficiency and the run-by-run FSD+GEM multiplicity drift, then re-extract $v_1(y_{cm})$ and its midrapidity slope for the revised 10-30% bin. If the slope moves outside the band of published values, or the JAM comparison degrades by more than the quoted uncertainties, the paper's central quantitative claim is falsified.
Extended reading notes
Core claim
The paper claims that the first directed-flow ($v_1$) measurement of protons from Xe+Cs(I) collisions at 3.8 AGeV has been obtained, for the 10-30% centrality bin. Using the first-order event plane reconstructed from spectator energy deposition in the FHCal, with three- and four-sub-event resolution corrections and acceptance corrections applied to both the flow vectors and the proton $u_1$ vectors, the authors measure $v_1$ as a function of rapidity and $p_T$. The resulting $v_1(y_{cm})$ is roughly reproduced by JAM in the mean-field mode with a momentum-dependent potential, while the extracted midrapidity slope $dv_1/dy_{cm}|_{y_{cm}=0}$ agrees within uncertainties with the published trend from STAR, HADES, and FOPI. The paper states these are first results and that the centrality estimate is very preliminary, so the final quantitative conclusions await the pending corrections.
Load-bearing premise
The load-bearing premise is that the 10-30% centrality selection is correct, because the measured $v_1$ is quoted for that bin and the note itself says the centrality estimate is very preliminary.
Editorial extensions
If this is right
- BM@N can produce differential flow measurements of identified protons in the $\sqrt{s_{NN}} \approx 2$--$3.5$ GeV range, adding a new data point to the energy and system-size dependence of $v_1$.
- The measured $v_1(y_{cm})$ provides a benchmark for JAM with a momentum-dependent mean field in the Xe+Cs system, extending model tests beyond Au+Au.
- The midrapidity slope of proton $v_1$ at $\sqrt{s_{NN}} = 3.26$ GeV can be added to the existing compilation of $dv_1/dy$ versus collision energy, and it is consistent with the published trend.
- The agreement among resolution factors and $v_1$ values from different FHCal sub-event planes indicates that the spectator-based event plane works reliably despite the beam-hole leakage.
- The final quantitative results require the pending centrality corrections, as the note states that the centrality estimate is very preliminary.
Reading between the lines
- One testable extension is to apply the same analysis chain to $v_2$ and to pions and kaons in the same dataset; such multi-differential data would constrain the equation of state more strongly than $v_1$ alone.
- If the preliminary centrality is later shifted by the trigger-efficiency and multiplicity corrections, the quantitative comparison with JAM could change, but the qualitative $v_1(y)$ trend is likely stable because the directed-flow slope has only weak centrality dependence except for the most central bin.
- A direct comparison of Xe+Cs with Au+Au at matched center-of-mass energy would isolate the system-size dependence of spectator shadowing, a question the paper motivates but does not quantify.
- One could also cross-check the spectator-plane resolution with a larger rapidity gap or a different sub-event grouping; if the extracted $v_1$ changed by more than the estimated non-flow contribution, the event-plane assumption would need revision.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This analysis note reports the first measurement of the directed flow v1 of protons in 10-30% central Xe+Cs(I) collisions at a beam energy of 3.8 AGeV (BM@N run8). The authors describe the run-by-run quality assurance, event and track selection, proton identification via TOF, centrality determination with MC-Glauber and Gamma-fit methods, event-plane reconstruction using the FHCal spectator calorimeter, and the scalar-product extraction of v1 as a function of rapidity and transverse momentum. The measured v1(y_cm) is compared with JAM transport model calculations with a momentum-dependent mean field, and the midrapidity slope dv1/dy is compared with published proton directed-flow slopes from other experiments. The paper is written as an internal analysis note and includes extensive figures, systematic checks, and links to analysis code and data paths.
Significance. If validated, this is a genuinely new observable: the first directed-flow measurement for the intermediate-mass, asymmetric Xe+Cs(I) system at NICA energies, and it can serve as a system-size and asymmetry benchmark for transport-model comparisons in the 2-5 GeV energy range. The note has real strengths: a very detailed run-by-run QA section, explicit event-selection statistics, a realistic GEANT4-based correction framework, public code repositories for the QA and centrality frameworks, and a systematic-uncertainty study covering tracking, PID, DCA, off-target collisions, acceptance/efficiency, and run-period stability. The main caveats are that the centrality calibration is explicitly labeled preliminary, the final v1 points and slopes are not quoted numerically, and there is an inconsistency in the resolution-correction formulas that directly affects the absolute normalization of v1.
major comments (4)
- [Section 3.4, final paragraph] The centrality estimate is explicitly stated to be 'very preliminary', and the text lists the missing ingredients: CCT2 trigger efficiency, run-dependent changes in average FSD+GEM multiplicity, and systematics from the MC-Glauber versus Gamma-fit choices. Since the entire physics claim is made for the 10-30% centrality bin, the comparison with the JAM 10-30% curve in Figure 47 and with the published slope systematics in Figure 48 is conditional on this calibration. If the true centrality of the analyzed sample is shifted, the model comparison and the energy-dependent slope comparison are not apples-to-apples. This is the single most load-bearing limitation and must be resolved before the result can be considered final.
- [Section 5, Figures 47-48] The paper never gives the numerical values of the v1(y_cm) and v1(pT) points, the fitted polynomial coefficients a, b, c from v1 = a + b y_cm + c y_cm^3, the extracted slope dv1/dy|y=0, or the corresponding statistical and systematic uncertainties. The claims of 'roughly captures' the JAM trend and 'reasonable agreement' with published slopes cannot be quantitatively checked from the text alone. A table with the binned v1 values, their uncertainties, the fit parameters, and the fit quality (e.g., chi2/ndf) should be added.
- [Section 4.1, Eq. (25) and Section 4.3, Eq. (28)] The resolution-correction formulas are inconsistent as printed. Equation (25), R^y_1{a(b,c)} = sqrt( 2<Y_b Y_c> / (2<Y_a Y_b> 2<Y_a Y_c>) ), is dimensionally inconsistent for a three-subevent resolution and does not match Eq. (28), R^y_1{a(b,c)} = sqrt( <Y_a Y_b><Y_a Y_c>/<Y_b Y_c> ). Moreover, if the latter is used together with v1 = 2<y1 Y_a*>/R^y_1, the standard scalar-product derivation for y-components gives R_y = sqrt(2<Y_aY_b><Y_aY_c>/<Y_bY_c>), i.e. a factor sqrt(2) larger than Eq. (28). As written, the formula would overestimate v1 by sqrt(2), which directly scales every v1 point and the extracted slope. The implemented formula needs to be stated unambiguously and validated, for example by reproducing the known JAM input v1 in the closed-symbol curves of Figure 37.
- [Section 3.3, QA run removal] The 3-sigma bad-run rejection (Section 3.2) removes about 18M events based on the deviation of run-averaged observables from the global mean. If bad runs cluster in specific time periods, this can bias the multiplicity distribution used for centrality, especially because the note later states that the average FSD+GEM multiplicity changed during run8. The magnitude of this effect should be quantified as part of the centrality systematics, or at least discussed in Section 3.4.
minor comments (7)
- [Abstract] The abstract says 'The systematic uncertainty study will also be presented and discussed', but the study is in fact included in Section 4.4; the future tense should be replaced with a present-tense statement.
- [Section 3.4] In the paragraph describing the application to experimental data, the sentence 'Figure shows the results' is missing the figure number; it should refer to Figure 28.
- [Throughout] There are numerous typos and grammatical errors that should be cleaned up: 'standart', 'creats', 'resgion', 'colid', 'bellow', 'calculatad', 'persented', 'previus', and 'pior' among others. These do not affect the physics but make the note harder to read.
- [Section 4.1, Eqs. (17)-(19)] The notation for the three- and four-subevent resolutions would benefit from a short derivation or a reference to the standard scalar-product formalism, because the current text moves quickly from Q-vector correlations to the component form and the definitions of a, b, c, d are implicit.
- [Section 4.4, first bullet] The statement that the momentum-reconstruction systematic uncertainty is 'bellow 2-5%' is overly broad; it would be clearer to report the uncertainty separately for the pT and rapidity ranges where it was evaluated, and to state whether the quoted range is a maximum or a typical value.
- [Section 5, Figure 47] The right panel shows v1 as a function of pT, but the text does not describe the pT dependence in words or compare it with the known behavior from HADES/STAR; a brief qualitative statement would help the reader interpret the figure.
- [References] The note relies heavily on two self-citations from the same group (Refs. [17] and [21]) for the JAM model setup and for the performance studies. This is acceptable, but the paper should state explicitly which JAM parameters and EOS variants were used in the comparison shown in Figure 47, so that the model curve is reproducible without consulting the unpublished analysis-note chain.
Circularity Check
No circularity: the v1 measurement is self-contained, and the JAM comparison and slope extraction are external benchmarks rather than fitted inputs.
full rationale
The central observable, directed flow v1 of protons, is measured directly from BM@N run8 data using standard event-plane and scalar-product techniques (Eqs. 27-29) with QnTools acceptance corrections and resolution correction factors from three- and four-sub-event methods. No parameter of the flow analysis is fitted to the v1 result, and no equation reduces v1 to an input assumption. The comparison with the JAM transport model is a benchmark, not a fit: the model is generated independently at fixed EOS options, and the note states only that JAM 'roughly captures the overall magnitude and trend' of the measured signal. The same-group papers cited for the JAM mean-field setup ([17, 21]) supply model parameters and prior validation against HADES/STAR data; they do not impose the BM@N v1 values by construction. The centrality determination is indeed labelled 'very preliminary' and depends on fits to the FSD+GEM multiplicity distribution, but this affects the interpretation of the 10-30% bin, not the derivation of v1 itself; an incorrect centrality mapping would be a systematic uncertainty, not a circular step. The midrapidity slope is extracted by a polynomial fit to measured v1 and then compared with published slopes, which is a data summary rather than a self-referential prediction. No self-definitional, fitted-input-as-prediction, or self-citation-chain circularity is present.
Assumptions & free parameters
free parameters (3)
- MC-Glauber NBD parameters f, mu, k =
not quoted
- Gamma-fit parameters theta, k0, a1, a2, a3 =
not quoted
- v1(y_cm) polynomial coefficients a, b, c =
not quoted
assumptions (6)
- standard math The azimuthal distribution of particles can be expanded as 1 + 2 sum_n v_n cos(n(phi - Psi_R)) with a well-defined reaction plane Psi_R.
- domain assumption Spectator fragments carry a positive directed flow v1 > 0 in the forward rapidity region and deposit energy in FHCal modules that tracks the first-order event plane.
- domain assumption Non-flow correlations are suppressed by requiring pseudorapidity separation between the sub-events used in the three- and four-sub-event resolution formulas.
- domain assumption The nuclear density profiles for Xe and Cs and the inelastic NN cross section sigma_inel = 27.7 mb from the PHOBOS MC-Glauber model describe the collision geometry.
- domain assumption The multiplicity fluctuations at fixed impact parameter follow the NBD parameterization (MC-Glauber) or a gamma distribution (Gamma-fit method).
- domain assumption The JAM transport model with momentum-dependent RQMD.RMF mean field gives a realistic description of the collision dynamics.
Cite this review
Pith. "Pith review of Analysis Note: Directed flow $v_1$ of protons in the Xe+Cs(I) collisions at 3.8 AGeV." pith.science (2026). https://pith.science/paper/P3TEEHQ4
@misc{pith2026241208570,
author = {Pith},
title = {Pith review of: Analysis Note: Directed flow $v_1$ of protons in the Xe+Cs(I) collisions at 3.8 AGeV},
year = {2026},
howpublished = {\url{https://pith.science/paper/P3TEEHQ4}},
note = {Machine review of arXiv:2412.08570}
}
abstract
In this note, we present the directed flow $v_1$ measurements of protons from Xe+Cs(I) collisions at 3.8 AGeV (BM@N run8). We show the datasets, event and track selection cuts, centrality definition, event plane reconstruction and resolution. The $v_1$ results are presented as function of transverse momentum ($p_T$) and rapidity ($y_{cm}$) for 10-30\% central Xe+Cs(I) collisions. The systematic uncertainty study will also be presented and discussed. The $v_1$ measurements are compared with results of JAM transport model calculations and published data from other experiments.
Figures
Figures from the paper (49 more)
Forward citations
Cited by 1 Pith paper
-
Bayesian analysis of properties of nuclear matter with the FOPI experimental data
Bayesian fits to FOPI Au+Au flow and stopping data yield m*/m0 around 0.78-0.88 and F around 0.75-0.88, while K0 remains unconstrained.
Reference graph
Works this paper leans on
- [1]
-
[2]
Poskanzer A. M., Voloshin S. A. Methods for analyzing anisotropic flow in relativistic nuclear collisions // Phys. Rev. C. — 1998. — Vol. 58. — P. 1671– 1678
work page 1998
-
[3]
Voloshin S.,Zhang Y.Flow study in relativistic nuclear collisions by Fourier ex- pansion of azimuthal particle distributions // Zeitschrift for Physik C Particles and Fields. — 1996. — May. — Vol. 70, no. 4. — P. 665–671
work page 1996
-
[4]
Voloshin S. A., Poskanzer A. M., Snellings R. Collective phenomena in non- central nuclear collisions // Landolt-Bornstein / ed. by R. Stock. — 2010. — Vol. 23. — P. 293–333. — arXiv:0809.2949 [nucl-ex]
arXiv 2010
-
[5]
Selyuzhenkov I., Voloshin S.Effects of non-uniform acceptance in anisotropic flow measurement // Phys. Rev. C. — 2008. — Vol. 77. — P. 034904. 48 Figure 43: Directed flow 𝑣1 of protons as a function of rapidity𝑦𝑐𝑚 measured for different values of the𝐷𝐶𝐴 cut and different n-𝜎 PID cuts for the proton identifi- cation: (𝑚2- ⟨︀ 𝑚2 𝑝 ⟩︀ ) < 1, 2, 3 𝜎𝑚2𝑝 cut (right)
work page 2008
- [6]
- [7]
- [8]
Show all 41 references
-
[9]
[et al.]
Reisdorf W. [et al.]. Systematics of azimuthal asymmetries in heavy ion colli- sions in the 1 A GeV regime // Nucl. Phys. A. — 2012. — Vol. 876. — P. 1– 60
2012
-
[10]
Proton, deuteron and triton flow measurements in Au+Au collisions at√𝑠NN = 2.4 GeV // Eur
Adamczewski-Musch J.[et al.]. Proton, deuteron and triton flow measurements in Au+Au collisions at√𝑠NN = 2.4 GeV // Eur. Phys. J. A. — 2023. — Vol. 59, no. 4. — P. 80. 49 Figure 44: Left: the distribution of the primary vertex in X-Y plane. Right: Directed flow 𝑣1 of protons a...
2023
-
[11]
[et al.]
Adam J. [et al.]. Flow and interferometry results from Au+Au collisions at √𝑠𝑁 𝑁= 4.5 GeV // Phys. Rev. C. — 2021. — Vol. 103, no. 3. — P. 034908
2021
-
[12]
S.[et al.]
Abdallah M. S.[et al.]. Disappearance of partonic collectivity in√𝑠NN = 3 GeV Au+Au collisions at RHIC // Phys. Lett. B. — 2022. — Vol. 827. — P. 137003
2022
-
[13]
G.Determination of the equation of state of dense matter // Science
Danielewicz P., Lacey R., Lynch W. G.Determination of the equation of state of dense matter // Science. — 2002. — Vol. 298. — P. 1592–1596
2002
-
[14]
Senger P.Heavy-Ion Collisions at FAIR-NICA Energies // Particles. — 2021. — Vol. 4, no. 2. — P. 214–226
2021
-
[15]
Constraining the nuclear matter equation of state around twice saturation density // Nucl
Le F` evre A.[et al.]. Constraining the nuclear matter equation of state around twice saturation density // Nucl. Phys. A. — 2016. — Vol. 945. — P. 112–133
2016
-
[16]
Sensitivity of Au+Au collisions to the symmetric nu- clear matter equation of state at 2–5 nuclear saturation densities // Phys
Oliinychenko D.[et al.]. Sensitivity of Au+Au collisions to the symmetric nu- clear matter equation of state at 2–5 nuclear saturation densities // Phys. Rev. C. — 2023. — Vol. 108, no. 3. — P. 034908. — arXiv:2208.11996 [nucl-th]. 50 Figure 45: Directed flow 𝑣1 of protons as ...
2023 arXiv
-
[17]
Mamaev M., Taranenko A.Toward the System Size Dependence of Anisotropic Flow in Heavy-Ion Collisions at√𝑠𝑁 𝑁= 2–5 GeV // Particles. — 2023. — Vol. 6, no. 2. — P. 622–637
2023
-
[18]
[et al.]
Afanasiev S. [et al.]. The BM@N spectrometer at the NICA accelerator com- plex. — 2023. — Dec. — arXiv:2312.17573 [hep-ex]
2023
-
[19]
Sharma S. R. First-Order Event Plane Correlated Directed and Triangular Flow from Fixed-Target Energies at RHIC-STAR // Universe. — 2024. — Vol. 10, no. 3. — P. 118
2024
-
[20]
Liu Z.Anisotropic Flow of Identified Particles in Au + Au Collisions at√𝑠𝑁 𝑁 = 3-3.9 GeV at RHIC. — 2023. — Dec. — arXiv:2312.16758 [nucl-ex]
2023 arXiv
-
[21]
Model Study of the Energy Dependence of Anisotropic Flow in Heavy-Ion Collisions at√𝑠𝑁 𝑁= 2–4.5 GeV // Particles
Parfenov P. Model Study of the Energy Dependence of Anisotropic Flow in Heavy-Ion Collisions at√𝑠𝑁 𝑁= 2–4.5 GeV // Particles. — 2022. — Vol. 5, no. 4. — P. 561–579. 51 Figure 46: Directed flow 𝑣1 of protons as a function of rapidity𝑦𝑐𝑚 measured in the different run periods (le...
2022
-
[22]
B.[et al.]
Larionov A. B.[et al.]. Squeezeout of nuclear matter in peripheral heavy ion collisions and momentum dependent effective interactions // Phys. Rev. C. —
-
[23]
Bass S. A. [et al.]. Microscopic models for ultrarelativistic heavy ion colli- sions // Prog. Part. Nucl. Phys. — 1998. — Vol. 41. — P. 255–369
1998
-
[24]
[et al.]
Aichelin J. [et al.]. Parton-hadron-quantum-molecular dynamics: A novel mi- croscopic 𝑛-bodytransportapproachforheavy-ioncollisions,dynamicalcluster formation, and hypernuclei production // Phys. Rev. C. — 2020. — Vol. 101, no. 4. — P. 044905
2020
-
[25]
[et al.]
Baznat M. [et al.]. Monte-Carlo Generator of Heavy Ion Collisions DCM- SMM // Phys. Part. Nucl. Lett. — 2020. — Vol. 17, no. 3. — P. 303–324
2020
-
[26]
JAM: an event generator for high energy nuclear collisions // EPJ Web Conf
Nara Y. JAM: an event generator for high energy nuclear collisions // EPJ Web Conf. / ed. by B. Pattison [et al.]. — 2019. — Vol. 208. — P. 11004
2019
-
[27]
Nara Y., Ohnishi A. Mean-field update in the JAM microscopic transport model: Mean-field effects on collective flow in high-energy heavy-ion collisions 52 Figure 47: Directed flow𝑣1 of protons in 10-30% central Xe+Cs(I) collisions at 3.8 A GeV as a function of rapidity𝑦𝑐𝑚 (lef...
2022
-
[28]
Nara Y., Maruyama T., Stoecker H.Momentum-dependent potential and col- lectiveflowswithintherelativisticquantummoleculardynamicsapproachbased on relativistic mean-field theory // Phys. Rev. C. — 2020. — Vol. 102, no. 2. — P. 024913
2020
-
[29]
Loizides C., Nagle J., Steinberg P.Improved version of the PHOBOS Glauber Monte Carlo // SoftwareX. — 2015. — Vol. 1/2. — P. 13–18
2015
-
[30]
Centrality determination of Pb-Pb collisions at√𝑠𝑁 𝑁= 2.76 TeV with ALICE // Phys
Abelev B.[et al.]. Centrality determination of Pb-Pb collisions at√𝑠𝑁 𝑁= 2.76 TeV with ALICE // Phys. Rev. C. — 2013. — Vol. 88, no. 4. — P. 044909
2013
-
[31]
[et al.]
Segal I. [et al.]. Using multiplicity of produced particles for centrality deter- mination in heavy-ion collisions with the CBM experiment // J. Phys. Conf. Ser. / ed. by P. Teterin. — 2020. — Vol. 1690, no. 1. — P. 012107. 53 Figure 48: The slope of 𝑣1 of protons at midrapidi...
2020
-
[32]
Centrality determination of Au + Au collisions at 1.23A GeV with HADES // Eur
Adamczewski-Musch J.[et al.]. Centrality determination of Au + Au collisions at 1.23A GeV with HADES // Eur. Phys. J. A. — 2018. — Vol. 54, no. 5. — P. 85
2018
-
[33]
A Centrality Detector Concept // Nucl
Tarafdar S., Citron Z., Milov A. A Centrality Detector Concept // Nucl. In- strum. Meth. A. — 2014. — Vol. 768. — P. 170–178
2014
-
[34]
Rogly R., Giacalone G., Ollitrault J.-Y.Reconstructing the impact parameter of proton-nucleus and nucleus-nucleus collisions // Phys. Rev. C. — 2018. — Vol. 98, no. 2. — P. 024902
2018
-
[35]
D.[et al.]
Frankland J. D.[et al.]. Model independent reconstruction of impact parameter distributions for intermediate energy heavy ion collisions // Phys. Rev. C. —
-
[36]
[et al.]
Parfenov P. [et al.]. Relating Charged Particle Multiplicity to Impact Param- eter in Heavy-Ion Collisions at NICA Energies // Particles. — 2021. — Vol. 4, no. 2. — P. 275–287
2021
-
[37]
Idrisov D., Parfenov P., Taranenko A.Centrality Selection Effect on Elliptic Flow Measurements in Relativistic Heavy-Ion Collisions at NICA Energies // Particles. — 2023. — Vol. 6, no. 2. — P. 497–514
2023
-
[38]
QnTools -A framework for multi-differential accep- tancecorrectionandanisotropicflowanalysis//https://github.com/HeavyIonAnalysis/QnTools.— 2024
Selyuzhenkov I., Kreis L. QnTools -A framework for multi-differential accep- tancecorrectionandanisotropicflowanalysis//https://github.com/HeavyIonAnalysis/QnTools.— 2024. 55
2024
-
[2000]
— Vol. 62. — P. 064611
-
[2002]
— Vol. 66. — P. 021901
-
[2021]
— Vol. 104, no. 3. — P. 034609. 54
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.