Pith. sign in

REVIEW 3 major objections 5 minor 37 references

In PYTHIA 8 Angantyr, d+Au activity selection produces a hard-probe bias: the most-active class is depleted and the peripheral class enhanced relative to an Ncoll-weighted reference.

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-01 05:25 UTC pith:IKTAXWT3

load-bearing objection A careful, honest generator study making a specific falsifiable prediction about activity-selected hard probes in d+Au; the central caveat is that the prediction is a property of Angantyr until the forward-activity–Ncoll correlation is validated against data. the 3 major comments →

arxiv 2607.27572 v1 pith:IKTAXWT3 submitted 2026-07-30 hep-ph nucl-exnucl-th

Photon-calibrated event-activity bias and subcollision geometry in d+Au collisions at sqrt{s_(NN)} = 200 GeV with PYTHIA 8 Angantyr

Muhammad Ajaz , Haifa I. Alrebdi , Muhammad Waqas This is my paper
classification hep-ph nucl-exnucl-th
keywords d+Au collisionsevent-activity selectionhard-bias factordirect photonsneutral pionsPYTHIA Angantyrcentrality biassubcollision geometry
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.

This paper asks whether classifying d+Au collisions by forward event activity distorts the population of hard probes such as direct photons and neutral pions. Using the PYTHIA 8 Angantyr event generator, the authors define a hard-bias factor that compares the probability for a hard-scattered event to land in an activity class with a minimum-bias reference weighted by the number of nondiffractive subcollisions. They find that the most-active class is depleted and the peripheral class is enhanced for both photons and pions, with the bias factor rising from about 0.5 in the most-active class to about 1.7 in the peripheral class. Photon-tagged events also sample a smaller mean impact parameter and more subcollisions than pion-tagged or inclusive hard events. The result matters because it quantifies, at generator level, a selection effect that experimental centrality-dependent measurements in small collision systems must understand.

Core claim

The central discovery is the pattern of hard-bias factors B_c^X: for the 0–20% most-active class, B_gamma = 0.517–0.543 and B_pi0 = 0.482–0.648, while for the 60–100% peripheral class, B_gamma = 1.652–1.735 and B_pi0 = 1.353–1.982. Relative to an N_coll^ND-weighted minimum-bias reference, the model therefore predicts a depleted most-active class and an enhanced peripheral class for both probes. A decomposition into three constructions (fully correlated events, Ncoll-reweighted reference, and a factorized diagnostic) shows that an explicitly geometric Ncoll weighting does not remove the class ordering, leaving residual hard–soft correlations of unspecified origin.

What carries the argument

The hard-bias factor B_c^X = P(c|X hard) / [PMB(c) ⟨N_coll⟩_c/⟨N_coll⟩_MB], computed with the realized nondiffractive subcollision count N_coll^ND from Angantyr event truth. The same event-level truth provides the weighted means of impact parameter and N_coll^ND that distinguish photon-, pion-, and inclusive-hard events.

Load-bearing premise

The predictions describe real d+Au collisions only if PYTHIA 8 Angantyr correctly represents the correlation between hard scattering and forward multiplicity; the pp tune was transferred without retuning, and the calculation has no detector response or nuclear parton-distribution variation.

What would settle it

Measure the activity-class probabilities for direct-photon-tagged and pion-tagged events in d+Au collisions at 200 GeV and compute the hard-bias factor of Eq. (4) using an Ncoll-weighted minimum-bias reference; observing B values near unity in all classes, or a non-monotonic ordering across classes, would falsify the central claim.

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

If this is right

  • In the model, central-to-peripheral ratios of hard-probe yields inherit a class-dependent bias even after Ncoll weighting, so comparisons of activity-selected data to minimum-bias references need this correction.
  • Direct photons sample a more central geometry than pions, making the photon-tagged sample a different activity-selection reference than the pion-tagged sample.
  • The residual A/B ratio shows a hard–soft correlation beyond geometry, so simple Ncoll reweighting does not remove the activity bias.
  • The ordering and sign of the bias are robust to the choice of forward-activity estimator, though the absolute class probabilities change.
  • The prediction is generator-level and does not by itself establish collectivity, final-state energy loss, or a detector-level centrality response.

Where Pith is reading between the lines

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

  • If the bias pattern holds in data, measurements of centrality-dependent jet suppression in small systems would need to be reinterpreted: the activity selection alone produces peripheral enhancement without any final-state energy loss.
  • A direct testable extension is to compute the same hard-bias factors for charged hadrons or for jets with different radii, using the same event-level subcollision information; the authors note the necessary event-level trigger tags were not stored.
  • Because the residual correlation is not uniquely identified, tuning the A/B ratio to future data could distinguish among MPI, energy sharing, and forward-activity response models.
  • Comparing photon- and pion-tagged samples in the same activity class may be a cleaner diagnostic than either alone, since the two probes have different mean geometries.

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

3 major / 5 minor

Summary. The paper studies how event-activity selection biases hard-probe yields in PYTHIA 8 Angantyr d+Au collisions at 200 GeV. It defines a hard-bias factor B_c^X comparing the probability for a hard photon, pion, or jet event to fall into a fixed activity class against an Ncoll-weighted minimum-bias reference. The central results are the B_c^X values in Table III: the most-active 0–20% class is depleted (B_gamma ≈ 0.52–0.54, B_pi0 ≈ 0.48–0.65) and the peripheral 60–100% class is enhanced (B_gamma ≈ 1.65–1.74, B_pi0 ≈ 1.35–1.98). The paper also reports that photon-tagged events sample smaller impact parameters and more nondiffractive subcollisions than pion-tagged or inclusive HardQCD events, and it separates the explicit geometric contribution from residual hard–soft correlations via the A/B/C constructions. The analysis is explicitly generator-level, with no detector response and no nuclear-PDF variation, and the authors state that no unique microscopic origin for the residual is identified.

Significance. If the results are taken as a prediction of the Angantyr model, the paper provides a clean, falsifiable statement about how event-activity classes couple to hard probes in small systems. Its strengths are the large event sample (1.02e8 events), a transparent definition of the hard-bias factor, explicit propagation of MC statistical uncertainties, and a careful A/B/C closure that separates geometry weighting from residual correlations. The paper also compares minimum-bias photon and pion spectra to PHENIX data, although it does not validate the forward-activity–Ncoll relation that drives the central B_c^X pattern. The study is a useful generator-level benchmark for RHIC small-system analyses, even though its physical reach is limited by the absence of detector simulation and by the reliance on an unretuned pp tune for the d+Au activity description.

major comments (3)
  1. [Section VI, Eq. (4)] The definition of B_c^X is inconsistent with its description. The abstract and text call it the probability for a hard event to enter an activity class, but the nominal implementation weights selected photon/pi0 candidates, not events. Since the activity class is an event property, candidate weighting can bias the numerator when an event contains multiple candidates, and for pi0 samples this bias is likely correlated with activity. The event-triggered alternative is relegated to the Supplemental Material; the main text should report both or explicitly define B_c^X as a per-candidate quantity. Without this, the central numbers in Table III are not clearly connected to the stated observable.
  2. [Sections V and Fig. 4] The paper does not actually compare the activity-class dependence of the photon-normalized pion double ratio to PHENIX data, although PHENIX has centrality-dependent gamma/pi0 measurements in Ref. [8]. The text says the pp denominator is missing, but one can still compare the class dependence after normalizing each class ratio to the minimum-bias value. As written, the model's class-dependent gamma/pi0 behavior is not anchored to data, and the central B_c^X prediction lacks a direct observable check.
  3. [Sections III and IX] The only generator validation shown is the minimum-bias photon and pi0 spectra (Fig. 2), which constrain the hard process but not the forward-activity–Ncoll correlation that drives B_c^X. The paper explicitly states that there is no detector response, no nPDF variation, and that the Monash pp tune is transferred to d+Au without retuning. If the depleted-most-active/enhanced-peripheral pattern is meant to describe real d+Au collisions, the model's activity-class probabilities should be validated against measured BBC/MPC distributions, or the conclusions should be more explicitly restricted to the generator. Section X does contain this caveat, but the abstract and the main interpretation should be equally explicit.
minor comments (5)
  1. [Title/Abstract] 'ind+Au collisions' should read 'in d+Au collisions.'
  2. [Table I] The notation for Ncoll is inconsistent: the table uses 'N^ND_coll' and 'N_tot_coll' with varying superscript/subscript conventions. Please standardize the notation.
  3. [Section V, Fig. 4] The caption says 'The experimental comparison uses a minimum-bias-normalized diagnostic' but does not specify which data are overlaid. Please identify the data points and the normalization procedure.
  4. [Section IX, Table V] The 'Nch=11 convention' is listed as a separate systematic but is never defined in the text. Please define this convention or remove it.
  5. [Fig. 6] The open orange rings are not explained in the caption. Define what they indicate.

Circularity Check

0 steps flagged

No significant circularity: the hard-bias factors are generator outputs, not fitted inputs; self-citations are contextual and non-load-bearing.

full rationale

The central observable B_c^X (Eq. 4) is defined as P(c|X hard) divided by the Ncoll-weighted minimum-bias reference PMB(c)<Ncoll>_c/<Ncoll>_MB. The numerator is the generator's class probability for hard events and the denominator is an independently formed reference; B=1 is not enforced by construction, and the reported depleted/enhanced ordering is an output of the simulation. No parameter is fitted to B_c^X, to the class probabilities, or to the photon/pion spectra used in Fig. 2; the comparison to PHENIX data is a benchmark, not a calibration of the bias. The A/B/C closure in Sec. VIII is explicitly a diagnostic: the paper says A/B "does not isolate a unique microscopic mechanism" and notes the residual is not separated into sources, so it does not disguise a definition as a result. Self-citations [31], [36], [37] support only contextual statements about tune/fragmentation/underlying-event sensitivity and are not used as a uniqueness theorem or as the proof of the main prediction; they are not load-bearing. The acknowledged limitations (Sec. IX: no detector response, no nPDF variation; Sec. X: generator-level, no unique microscopic origin) are external-validity caveats rather than circularity: they do not make any equation reduce to its inputs. Score 1 reflects only the presence of minor non-load-bearing self-citations; the derivation itself is self-contained.

Axiom & Free-Parameter Ledger

5 free parameters · 4 axioms · 0 invented entities

The paper introduces no new free parameters fitted to data; the listed items are analysis conventions chosen by hand that affect the quantitative values. The central claims depend on the domain assumption that Angantyr's MPI/subcollision model faithfully captures the hard–soft correlation in d+Au. No new entities (particles, forces, dimensions) are postulated.

free parameters (5)
  • Activity class boundaries = 0–20%, 20–40%, 40–60%, 60–100%
    Fixed multiplicity thresholds chosen by hand (Section IV). The class ordering is robust to boundary variations, but the absolute class probabilities depend on these cuts.
  • STAR-like estimator acceptance = -3.8 < η < -2.8
    Chosen to reproduce STAR-like Au-going acceptance (Section IV); affects class probabilities but not ordering.
  • PHENIX-like estimator acceptance = BBC-S acceptance proxy
    Particle-level proxy for PHENIX BBC-S (Section IV); no detector response. Class probabilities differ between estimators.
  • Jet constituent pT threshold = pT > 0.4 GeV/c
    Nominal constituent cut (Table I); variants set it to zero or remove muons, with bounded changes to jet results.
  • Nch=11 convention = Nch=11
    Separate systematic convention (Table V).
axioms (4)
  • domain assumption PYTHIA 8.316 Angantyr with the Monash tune, transferred from pp without retuning, provides a reliable model of d+Au collisions at 200 GeV.
    Load-bearing for the entire study; invoked in Section III. The paper does not quantitatively validate the model against the PHENIX/STAR data it compares with.
  • domain assumption The N_coll^ND-weighted minimum-bias sample is the correct reference for the geometric expectation.
    The hard-bias factor Eq. (4) normalizes by <N_coll>_c/<N_coll>_MB; if the true geometric baseline differs from realized nondiffractive subcollision scaling, the reported depletion/enhancement is a definitional artifact.
  • ad hoc to paper The direct-photon sample defined as HARD ME + ISR nonhadronic + FSR fragmentation is a valid hard-scattering reference.
    Photon ancestry is fixed (Table I) and not varied; fragmentation photons are included, so the 'direct photon' tag is a generator definition, not the experimental isolation definition.
  • domain assumption The anti-kT algorithm with E-scheme recombination and inclusive visible primaries is the correct jet definition for comparison with PHENIX R=0.3 jets.
    Jet observables depend on this recipe (Section IV); no detector simulation is applied.

pith-pipeline@v1.3.0-daily-deepseek · 7944 in / 14487 out tokens · 127439 ms · 2026-08-01T05:25:00.749689+00:00 · methodology

0 comments
read the original abstract

Event-activity selections in small nuclear collision systems couple collision geometry to the soft response accompanying a hard scattering. We examine this coupling in d+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV with PYTHIA 8.316 Angantyr, using direct photons, terminal pre-decay neutral pions, and anti-kT jets. A hard-bias factor compares the probability for a hard event to enter an activity class with an $N^{ND}_{coll}$ -weighted minimum-bias reference. The model predicts a depleted most-active class and an enhanced peripheral class for both photons and pions. Photon-tagged events also sample a smaller mean impact parameter and more nondiffractive subcollisions than pion-tagged or inclusive HardQCD events. Comparing fully correlated Angantyr events with an Ncoll-reweighted minimum-bias reference and a factorized diagnostic separates the explicit geometric contribution from residual hard-soft correlations. STAR-like and PHENIX-like particle-level activity proxies preserve the class ordering but give different class probabilities. The calculation is generator-level and does not identify a unique microscopic origin for the residual correlations.

Figures

Figures reproduced from arXiv: 2607.27572 by Haifa I. Alrebdi, Muhammad Ajaz, Muhammad Waqas.

Figure 1
Figure 1. Figure 1: FIG. 1. Weighted hard-event geometry from Angantyr event [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: shows the minimum-bias spectra. Their MC statistical precision ranges from 0.58% to 1.17% for direct photons and from 2.59% to 6.36% for pions. Model and experimental uncertainties are displayed separately. With the normalization in Eq. (1), the model predicts Q γ dAu = 1.586–1.715 and Qπ 0 dAu = 1.843–2.180 ( [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Photon-normalized pion double ratio from Eq. (2). All [PITH_FULL_IMAGE:figures/full_fig_p003_4.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Jet [PITH_FULL_IMAGE:figures/full_fig_p004_6.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Hard-bias factors for direct photons and terminal [PITH_FULL_IMAGE:figures/full_fig_p004_5.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Ratio of the fully correlated Angantyr construction [PITH_FULL_IMAGE:figures/full_fig_p005_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Direct-photon class probabilities for the STAR-like and [PITH_FULL_IMAGE:figures/full_fig_p005_9.png] view at source ↗

discussion (0)

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

Reference graph

Works this paper leans on

37 extracted references · 17 linked inside Pith

  1. [1]

    M. L. Miller, K. Reygers, S. J. Sanders, and P. Steinberg, Annual Review of Nuclear and Particle Science57, 205 (2007), arXiv:nucl-ex/0701025

    Pith/arXiv arXiv 2007

  2. [2]

    PHENIX Collaboration, Physical Review C90, 034902 (2014)

    2014

  3. [3]

    Kordell and A

    M. Kordell and A. Majumder, Physical Review C97, 054904 (2018)

    2018

  4. [4]

    Bzdak, V

    A. Bzdak, V. Skokov, and S. Bathe, Physical Review C 93, 044901 (2016)

    2016

  5. [5]

    Alvioli and M

    M. Alvioli and M. Strikman, Physics Letters B722, 347 (2013), arXiv:1301.0728

    Pith/arXiv arXiv 2013

  6. [6]

    McGlinchey, J

    D. McGlinchey, J. L. Nagle, and D. V. Perepelitsa, Physi- cal Review C94, 024915 (2016)

    2016

  7. [7]

    PHOBOS Collaboration, Physical Review Letters93, 082301 (2004), arXiv:nucl-ex/0311009

    Pith/arXiv arXiv 2004

  8. [8]

    PHENIX Collaboration, Physical Review Letters134, 022302 (2025), arXiv:2303.12899

    arXiv 2025

  9. [9]

    PHENIX Collaboration, Physical Review Letters116, 122301 (2016), arXiv:1509.04657

    arXiv 2016

  10. [10]

    STAR Collaboration, Physical Review C76, 054903 (2007), arXiv:nucl-ex/0609021

    Pith/arXiv arXiv 2007

  11. [11]

    PHENIX Collaboration, Physical Review C88, 024906 (2013), arXiv:1304.3410

    arXiv 2013

  12. [12]

    PHENIX Collaboration, Physical Review C105, 064902 (2022)

    2022

  13. [13]

    PHENIX Collaboration, Physical Review D86, 072008 (2012)

    2012

  14. [14]

    PHENIX Collaboration, Physical Review C87, 054907 (2013)

    2013

  15. [15]

    Arleo, K

    F. Arleo, K. J. Eskola, H. Paukkunen, and C. A. Salgado, Journal of High Energy Physics04, 055 (2011)

    2011

  16. [16]

    Sj¨ ostrand, S

    T. Sj¨ ostrand, S. Ask, J. R. Christiansen, R. Corke, N. De- sai, P. Ilten, S. Mrenna, S. Prestel, C. O. Rasmussen, and P. Z. Skands, Computer Physics Communications191, 159 (2015), arXiv:1410.3012

    Pith/arXiv arXiv 2015

  17. [17]

    Bierlich, S

    C. Bierlich, S. Chakraborty, N. Desai, L. Gellersen, I. Helenius, P. Ilten, L. L¨ onnblad, S. Mrenna, S. Pres- tel, C. T. Preuss, T. Sj¨ ostrand, P. Skands, M. Utheim, and R. Verheyen, SciPost Physics Codebases , 8 (2022), arXiv:2203.11601

    Pith/arXiv arXiv 2022

  18. [18]

    Bierlich, G

    C. Bierlich, G. Gustafson, L. L¨ onnblad, and H. Shah, Journal of High Energy Physics10, 134 (2018), arXiv:1806.10820

    Pith/arXiv arXiv 2018

  19. [19]

    Sj¨ ostrand and M

    T. Sj¨ ostrand and M. van Zijl, Physical Review D36, 2019 (1987)

    2019

  20. [20]

    Andersson, G

    B. Andersson, G. Gustafson, G. Ingelman, and T. Sj¨ ostrand, Physics Reports97, 31 (1983)

    1983

  21. [21]

    S. R. Nayak, A. Das, and B. K. Singh, European Physical Journal C86, 845 (2026), arXiv:2503.23019

    Pith/arXiv arXiv 2026

  22. [22]

    PHENIX Collaboration, Disentangling centrality bias and final-state effects in the production of high- pt π0 using direct γ in d+au collisions at √sN N= 200 gev (2025), version 1

    2025

  23. [23]

    PHENIX Collaboration, Centrality-dependent modifica- tion of jet-production rates in deuteron-gold collisions at√sN N= 200 gev (2025), version 1

    2025

  24. [24]

    STAR Collaboration, Rapidity and species dependence of particle production at large transverse momentum for d+au collisions at √sN N= 200 gev (2021), version 1

    2021

  25. [25]

    PHENIX Collaboration, Spectra and ratios of identified particles in au+au and d+au collisions at √sN N= 200 gev (2020), version 1. 7

    2020

  26. [26]

    Skands, S

    P. Skands, S. Carrazza, and J. Rojo, European Physical Journal C74, 3024 (2014), arXiv:1404.5630

    Pith/arXiv arXiv 2014

  27. [27]

    Buckley, H

    A. Buckley, H. Hoeth, H. Lacker, H. Schulz, and J. E. von Seggern, European Physical Journal C65, 331 (2010), arXiv:0907.2973

    Pith/arXiv arXiv 2010

  28. [28]

    J. R. Christiansen and P. Z. Skands, Journal of High Energy Physics08, 003 (2015), arXiv:1505.01681

    Pith/arXiv arXiv 2015

  29. [29]

    L¨ onnblad and H

    L. L¨ onnblad and H. Shah, European Physical Journal C 83, 575 (2023), arXiv:2303.11747

    Pith/arXiv arXiv 2023

  30. [30]

    L¨ onnblad and H

    L. L¨ onnblad and H. Shah, European Physical Journal C 83, 639 (2023)

    2023

  31. [31]

    Waqar, H

    M. Waqar, H. I. Alrebdi, M. Waqas, K. S. Al-Mugren, and M. Ajaz, Chinese Physics C48, 093109 (2024), arXiv:2406.17384

    Pith/arXiv arXiv 2024

  32. [32]

    Bierlichet al., SciPost Physics8, 026 (2020), arXiv:1912.05451

    C. Bierlichet al., SciPost Physics8, 026 (2020), arXiv:1912.05451

    Pith/arXiv arXiv 2020

  33. [33]

    Bierlichet al., European Physical Journal C80, 485 (2020), arXiv:2001.10737

    C. Bierlichet al., European Physical Journal C80, 485 (2020), arXiv:2001.10737

    Pith/arXiv arXiv 2020

  34. [34]

    Cacciari, G

    M. Cacciari, G. P. Salam, and G. Soyez, Journal of High Energy Physics04, 063 (2008), arXiv:0802.1189

    Pith/arXiv arXiv 2008

  35. [35]

    Cacciari, G

    M. Cacciari, G. P. Salam, and G. Soyez, European Physi- cal Journal C72, 1896 (2012), arXiv:1111.6097

    Pith/arXiv arXiv 2012

  36. [36]

    H. I. Alrebdi, M. Ajaz, M. Waqas, M. A. Ahmad, M. Waqar, A. M. Quraishi, J. H. Baker, S. Jagnandan, and A. Jagnandan, European Physical Journal Plus140, 371 (2025)

    2025

  37. [37]

    Waqar, H

    M. Waqar, H. I. Alrebdi, M. Waqas, M. A. Ahmad, and M. Ajaz, European Physical Journal Plus140, 523 (2025)

    2025