{"id":"e61bae5d-493a-4c21-aad9-50764fe99237","arxiv_id":"2412.08570","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"First directed flow v1 measurement of protons in Xe+Cs(I) at 3.8 AGeV from BM@N, compared with JAM and earlier flow data.","lead":"This note reports the first measurement of proton directed flow v1 in Xe+Cs(I) collisions at 3.8 AGeV recorded by the BM@N experiment. It adds a new system-size data point at NICA energy and compares it with JAM transport calculations and Au+Au flow data.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Centrality-dependent physics claims rest on a self-declared 'very preliminary' centrality calibration; an incorrect 10-30% bin mapping could shift v1 comparisons with JAM and published slopes.","rationale":"The paper's strongest claim is explicitly about 10-30% central Xe+Cs(I) collisions, so the physical interpretation (JAM comparison and slope vs energy) is conditional on the centrality selection being correct. The note itself concedes in Section 3.4 that the centrality is 'very preliminary' and enumerates the missing corrections: CCT2 trigger efficiency, run8 time-dependent FSD+GEM multiplicity drift, and method systematics. This is exactly the weakest assumption flagged by the reader. No internal inconsistency in the flow extraction itself was found; the v1 analysis uses standard scalar-product/event-plane methods with QnTools corrections, three/four-subevent resolution corrections, and a reasonable set of systematic checks. The concern is a calibration completeness issue, not a methodological error. The concrete test of applying a trigger-efficiency weight to the centrality fit would settle whether the 10-30% label is robust; if robust, the central physics claim stands as stated. Therefore the verdict should remain CONDITIONAL, matching the reader's assessment.","tokens_in":22715,"tokens_out":1447,"duration_ms":13789,"concrete_test":"Re-derive the centrality calibration from the run8 experimental Nch distribution with an explicit CCT2 trigger efficiency ε(Nch), applied as a weight in the fits for both the MC-Glauber and Gamma-Fit methods, and propagate the resulting 10-30% selection boundaries into <b> and the v1(y_cm) points in Figure 47. If the 10-30% bin boundaries shift by more than ~1-2% in centrality fraction, or if the recomputed v1(y_cm) curve moves outside the quoted systematic band, the paper should explicitly state the sensitivity of the JAM comparison and slope extraction to the centrality calibration.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim is the first measurement of proton v1 in 10-30% central Xe+Cs(I) collisions at 3.8 AGeV, and the key interpretive steps—comparison with JAM in the same centrality class and the midrapidity slope dv1/dy|0 versus published Au+Au data—both depend on the centrality classification. Section 3.4 explicitly states the centrality estimate is 'very preliminary' and that final results require evaluating the CCT2 trigger efficiency, accounting for run8 drift in average FSD+GEM multiplicity, and assessing MC-Glauber vs Gamma-Fit systematics. The two methods are applied to data only after validation on DCM-QGSM-SMM Monte Carlo, but the experimental Nch distribution is analyzed as if it were purely from CCT2-triggered events without a trigger efficiency correction. If the true centrality of the analyzed event sample is shifted (e.g., 5-25% rather than 10-30%), the comparison with the JAM 10-30% curve in Figure 47 and the slope comparison in Figure 48 would not be apples-to-apples, since v1 is known to vary with centrality near mid-central collisions. The reader identified this same weakest assumption. This is a condition for the physics claim, not an internal inconsistency.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":23002,"tokens_out":7115,"duration_ms":70717,"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":[{"comment":"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":"Section 3.4, final paragraph"},{"comment":"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":"Section 5, Figures 47-48"},{"comment":"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":"Section 4.1, Eq. (25) and Section 4.3, Eq. (28)"},{"comment":"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.","section":"Section 3.3, QA run removal"}],"minor_comments":[{"comment":"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":"Abstract"},{"comment":"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.","section":"Section 3.4"},{"comment":"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":"Throughout"},{"comment":"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":"Section 4.1, Eqs. (17)-(19)"},{"comment":"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":"Section 4.4, first bullet"},{"comment":"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.","section":"Section 5, Figure 47"},{"comment":"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.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"This is a preliminary analysis note rather than a finished journal article. The strongest concern for the editor is that the central physics claim is explicitly conditional on a 'very preliminary' centrality calibration, and the resolution-formula inconsistency in Sections 4.1 and 4.3 could affect the absolute normalization of v1. Both issues are fixable within the manuscript's scope, so I recommend major revision rather than rejection. The paper would also be much stronger if the numerical values of the v1 points and the extracted slope were made available in the text or in an appendix."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this is a real first measurement, not a rehash. The proton v1 for Xe+Cs(I) at 3.8 AGeV from BM@N run8 does not exist anywhere else in the literature, and the note uses standard scalar-product/event-plane methods with a careful sub-event resolution treatment. The systematic section is honest and itemized: track quality, DCA, PID, vertex position, TOF systems, run periods, efficiency weighting. The JAM-with-mean-field comparison is stated as \"roughly\" capturing the data, which is the right level of claim.\n\nThe main soft spot is exactly the one the note itself flags. Section 3.4 says the centrality estimate is \"very preliminary\" and lists the missing pieces: CCT2 trigger efficiency, run8 drift in FSD+GEM multiplicity, and MC-Glauber vs Gamma-Fit systematics. The physics comparison with JAM and the published dv1/dy energy dependence both use the 10-30% bin. Since v1 is known to vary with centrality, a miscalibrated bin could shift the comparison. Note that the v1 slope is reported to be only weakly centrality dependent except for very central collisions, so the effect might be modest, but it has to be demonstrated, not assumed.\n\nSecond soft spot: no numbers. The final v1(y) and v1(pT) points and the extracted slope dv1/dy are shown only in figures. The reader cannot fit the polynomial or compare with the literature without digitizing. For an analysis note this is a real deficiency, and it is easily fixed.\n\nThe reader's circularity concern is minor and I mostly disagree with it as a burden. The v1 signal is measured; the model comparison is downstream. There is some reliance on self-cited simulation work [17,21] from the same group, but those are public performance/model studies used as benchmarks, and the calibration methods are standard. Self-citation here is not a smoking gun.\n\nWho is this for: the BM@N/NICA flow community and people doing system-size dependence of flow at 2-5 GeV. It is a within-subfield data point, not a field reorganization. It deserves a serious referee. I would send it to review with the condition that the authors add numerical tables, state the centrality systematics as a caveat or update the calibration, and show the fit parameters for the slope.","headline":"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.","tokens_in":23510,"tokens_out":2112,"would_cite":true,"duration_ms":22510,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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…","keywords":["directed flow","proton v1","Xe+Cs collisions","BM@N","NICA","event plane","nuclear equation of state","JAM transport model"],"falsifier":"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.","tokens_in":22540,"feed_emoji":"⚛️","tokens_out":12862,"duration_ms":107569,"temperature":0.7,"pith_summary":"This analysis note reports the first measurement of the directed flow of protons ($v_1$) in Xe+Cs(I) collisions at 3.8 AGeV, for the 10-30% centrality class, using data from the BM@N run8. The authors extract $v_1$ as a function of center-of-mass rapidity and transverse momentum with the event-plane method, using spectator fragments registered in the forward hadron calorimeter (FHCal) as the symmetry plane. They find that the measured $v_1(y_{cm})$ is roughly described by the JAM transport model with a momentum-dependent mean field, and that the midrapidity slope $dv_1/dy_{cm}|_{y_{cm}=0}$ is consistent with the published energy dependence of proton directed flow from other experiments. The centrality estimate is explicitly preliminary: the note states that the CCT2 trigger efficiency, run-by-run FSD+GEM multiplicity changes, and the systematic uncertainties of the MC-Glauber and Gamma-Fit methods still need evaluation. The result extends flow measurements, a probe of the nuclear equation of state, to a new reaction system at NICA energies.","feed_headline":"First proton directed-flow measurement from Xe+Cs at 3.8 AGeV","feed_subtitle":"New NICA data match JAM and extend equation-of-state tests to a new system.","key_machinery":"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$.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the standard event-plane and scalar-product formalism for anisotropic flow.","marker":"[2]"},{"why":"Provides the flow-vector framework and the recentering, twist, and rescaling corrections used to handle the non-uniform BM@N acceptance.","marker":"[4]"},{"why":"Gives the non-uniform-acceptance correction method that the QnTools framework implements.","marker":"[5]"},{"why":"FOPI proton flow data in the 1 AGeV regime used as a comparison point for the midrapidity slope trend.","marker":"[9]"},{"why":"HADES Au+Au proton flow data at 2.4 GeV used for comparison of the rapidity dependence and slope.","marker":"[10]"},{"why":"STAR BES II proton flow results at 3 GeV used in the energy-dependence comparison of the midrapidity slope.","marker":"[12]"},{"why":"BM@N performance study for flow measurements and the JAM mean-field settings on which the present analysis and comparisons rely.","marker":"[17]"},{"why":"JAM event generator used for the simulated detector performance and for the model comparison.","marker":"[26]"},{"why":"Momentum-dependent mean-field potential in JAM that produces the directed-flow signal compared with the data.","marker":"[28]"},{"why":"MC-Glauber model used to construct the centrality classes via multiplicity fitting.","marker":"[29]"}],"fun_headline_variants":["First proton directed flow from Xe+Cs at NICA","Proton flow in Xe+Cs: first NICA measurement","NICA reveals first proton flow in Xe+Cs collisions","First proton flow from Xe+Cs at 3.8 AGeV","Proton v1 from Xe+Cs: first NICA result"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["First proton directed flow from Xe+Cs at NICA","Proton flow in Xe+Cs: first NICA measurement","NICA reveals first proton flow in Xe+Cs collisions","First proton flow from Xe+Cs at 3.8 AGeV","Proton v1 from Xe+Cs: first NICA result"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001275,"raw_usage":{"total_tokens":5178,"prompt_tokens":870,"completion_tokens":4308,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":486,"completion_tokens_details":{"reasoning_tokens":4215}},"tokens_in":486,"tokens_out":4308,"duration_ms":29916,"temperature":1.0,"reasoning_tokens":4215,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T17:43:52.121969+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"M., Voloshin S","cited_arxiv_id":null,"evidence_quote":"Supplies the standard event-plane and scalar-product formalism for anisotropic flow."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the non-uniform-acceptance correction method that the QnTools framework implements."},{"cited_title":"[et al.]","cited_arxiv_id":null,"evidence_quote":"FOPI proton flow data in the 1 AGeV regime used as a comparison point for the midrapidity slope trend."},{"cited_title":"Proton, deuteron and triton flow measurements in Au+Au collisions at√𝑠NN = 2.4 GeV // Eur","cited_arxiv_id":null,"evidence_quote":"HADES Au+Au proton flow data at 2.4 GeV used for comparison of the rapidity dependence and slope."},{"cited_title":"S.[et al.]","cited_arxiv_id":null,"evidence_quote":"STAR BES II proton flow results at 3 GeV used in the energy-dependence comparison of the midrapidity slope."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"BM@N performance study for flow measurements and the JAM mean-field settings on which the present analysis and comparisons rely."},{"cited_title":"JAM: an event generator for high energy nuclear collisions // EPJ Web Conf","cited_arxiv_id":null,"evidence_quote":"JAM event generator used for the simulated detector performance and for the model comparison."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Momentum-dependent mean-field potential in JAM that produces the directed-flow signal compared with the data."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"MC-Glauber model used to construct the centrality classes via multiplicity fitting."}],"review_version":1}