pith. machine review for the scientific record. sign in

arxiv: 2601.00944 · v3 · submitted 2026-01-02 · 🌌 astro-ph.HE · astro-ph.IM· hep-ex· physics.comp-ph

Recognition: 2 theorem links

· Lean Theorem

Probabilistic modeling of Cherenkov emission from particle showers

Ian Crawshaw , Tianlu Yuan , Emre Yildizci , Lu Lu , Anatoli Fedynitch
Authors on Pith no claims yet

Pith reviewed 2026-05-16 17:43 UTC · model grok-4.3

classification 🌌 astro-ph.HE astro-ph.IMhep-exphysics.comp-ph
keywords Cherenkov emissionparticle showersneutrino telescopesMonte Carlo simulationprobabilistic modelinglight yield fluctuationsicewater
0
0 comments X

The pith

Distributions of Cherenkov light yield parameters from particle showers capture event-to-event fluctuations in amplitude and shape along the shower axis.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper constructs distributions of parameters that describe the Cherenkov light yield produced by particle showers in ice or water. These distributions come from underlying Monte Carlo simulations of showers and permit sampling that reproduces observed variations in light output from one shower to the next. Full tracking of every secondary particle is too slow for large-scale work, so experiments rely on parametrized models; the new distributions add realistic fluctuations that earlier approximations omitted. A reader would care because neutrino telescopes must separate signal from background with high fidelity, and the added variability directly affects how accurately those events can be simulated. The approach exploits the broad universality of shower development across different media.

Core claim

The central claim is that parameter distributions derived from Monte Carlo simulations of particle showers can be sampled to reproduce the amplitude and longitudinal shape fluctuations of Cherenkov light emission in ice and water. This probabilistic description improves on existing parametrized approximations by incorporating event-to-event variability that is essential for accurate simulation of signal and background events in neutrino telescopes.

What carries the argument

Distributions of Cherenkov light yield parameters, constructed from Monte Carlo simulations, that are sampled to generate realistic fluctuations in amplitude and shape along the shower axis.

If this is right

  • Sampling the distributions produces shower light profiles with realistic event-to-event differences in total light and in the longitudinal development.
  • Neutrino telescope simulations gain accuracy for both signal and background events without the full computational cost of tracking every secondary particle.
  • The same distributions apply across ice and water because shower development is governed by the same underlying physics.
  • Event generators can draw from the distributions to produce varied shower realizations for large-scale detector studies.

Where Pith is reading between the lines

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

  • The approach could be extended to other media such as air if the underlying Monte Carlo shows similar parameter distributions.
  • Incorporating these sampled fluctuations into existing simulation frameworks would allow direct testing against observed event rates in operating telescopes.
  • Parameter correlations that survive the construction process might reveal simplified relations useful for analytic approximations.

Load-bearing premise

The parameter distributions derived from Monte Carlo simulations faithfully capture the true physical variability of showers without introducing biases or missing important correlations.

What would settle it

A side-by-side comparison of light-yield profiles generated by sampling the distributions against independent full Monte Carlo runs or against real data from neutrino telescopes that shows statistically significant mismatches in the amplitude or shape fluctuation statistics.

read the original abstract

Subatomic particles can interact with target nuclei in matter or decay in flight, and an individual high-energy particle can induce a particle shower composed of numerous, lower-energy secondaries. These particle showers broadly exhibit universality across diverse media, including air, water, ice, and other materials, with their development governed by the Standard Model. Full Monte Carlo simulation of particle showers, where each secondary is individually tracked and propagated, can be a computational challenge to perform at scale. Experiments thus resort to parametrized approximations when efficient simulation becomes necessary. Here, we construct distributions of parameters capable of describing the Cherenkov light yield from particle showers in ice or water. Sampling from the distributions allows for a much improved description of event-to-event fluctuations, in amplitude and shape, along the shower axis. Including these effects is essential for a more accurate simulation of signal and background events in current and next-generation neutrino telescopes.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The manuscript constructs probability distributions for parameters that describe the Cherenkov light yield from particle showers in ice or water, derived from underlying Monte Carlo simulations. Sampling from these distributions is presented as enabling a much improved description of event-to-event fluctuations in both amplitude and shape along the shower axis, which is argued to be essential for accurate and efficient simulation of signal and background events in neutrino telescopes.

Significance. If the parameter distributions can be shown to faithfully reproduce the joint variability and correlations present in full Monte Carlo showers without introducing bias, the method would offer a computationally efficient alternative to full shower tracking while preserving fluctuation statistics critical for neutrino telescope analyses. The approach builds on standard universality properties of electromagnetic and hadronic showers but currently lacks the quantitative benchmarks needed to establish its practical advantage.

major comments (2)
  1. [Abstract, §3] Abstract and §3 (parameter construction): The central claim that sampling yields a 'much improved description' of fluctuations is unsupported by any reported quantitative validation, such as Kolmogorov-Smirnov distances, RMS profile ratios, or direct comparisons of sampled vs. full Monte Carlo shower profiles. No error analysis or cross-validation against independent Monte Carlo runs is provided.
  2. [§4] §4 (results and sampling): No separation between the Monte Carlo data used to fit the parameter distributions and any test data is described, leaving the improvement claim vulnerable to circularity; the distributions may simply re-describe the same simulations without demonstrating generalization to physical variability.
minor comments (1)
  1. [§2] Notation for the Cherenkov yield parameters (e.g., amplitude and shape variables) is introduced without a consolidated table of symbols or explicit functional forms for the fitted distributions.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thoughtful comments on our manuscript. We address each major comment below and describe the revisions we plan to implement to strengthen the paper.

read point-by-point responses

  1. Referee: [Abstract, §3] Abstract and §3 (parameter construction): The central claim that sampling yields a 'much improved description' of fluctuations is unsupported by any reported quantitative validation, such as Kolmogorov-Smirnov distances, RMS profile ratios, or direct comparisons of sampled vs. full Monte Carlo shower profiles. No error analysis or cross-validation against independent Monte Carlo runs is provided.

    Authors: We agree that quantitative validation is essential to support the claim of improvement. In the revised manuscript, we will add direct comparisons of shower profiles generated by sampling from the parameter distributions against independent full Monte Carlo simulations. This will include metrics such as Kolmogorov-Smirnov distances for the distributions of key parameters, RMS differences in longitudinal profiles, and error analysis to demonstrate the fidelity of the sampled fluctuations. revision: yes

  2. Referee: [§4] §4 (results and sampling): No separation between the Monte Carlo data used to fit the parameter distributions and any test data is described, leaving the improvement claim vulnerable to circularity; the distributions may simply re-describe the same simulations without demonstrating generalization to physical variability.

    Authors: We acknowledge the need to clarify the data handling procedure. The parameter distributions were constructed from a comprehensive set of Monte Carlo simulations, and we will revise §4 to explicitly describe the use of a training set for fitting the distributions and a separate test set for validation. This will include cross-validation results to show that the sampled showers generalize beyond the fitting data and faithfully reproduce the variability seen in independent simulations. revision: yes

Circularity Check

0 steps flagged

No circularity: parametric distributions constructed from MC are used for sampling without definitional reduction

full rationale

The paper constructs parameter distributions from underlying Monte Carlo simulations to describe Cherenkov light yield and fluctuations along the shower axis. Sampling from these distributions is presented as enabling improved modeling of event-to-event variations. No equations, self-citations, or steps in the provided abstract reduce any claimed prediction or result to its inputs by construction (e.g., no fitted parameter is renamed as an independent prediction, and no uniqueness theorem or ansatz is smuggled via self-reference). The method follows standard practice for efficient parametric simulation of stochastic processes, remaining self-contained against external benchmarks like full MC runs.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The approach rests on the assumption that particle showers exhibit sufficient universality to allow a compact set of parameter distributions, with those distributions fitted from Monte Carlo simulations whose accuracy is taken as given.

free parameters (1)
  • Cherenkov yield distribution parameters
    Parameters of the probability distributions are constructed from simulation data and therefore constitute fitted quantities.
axioms (1)
  • domain assumption Particle showers exhibit universality across media including ice and water
    Invoked in the abstract as the basis for constructing broadly applicable distributions.

pith-pipeline@v0.9.0 · 5469 in / 1132 out tokens · 37544 ms · 2026-05-16T17:43:08.278969+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

18 extracted references · 18 canonical work pages · 4 internal anchors

  1. [1]

    Frank and I.E

    I.M. Frank and I.E. Tamm,Coherent visible radiation of fast electrons passing through matter,Compt. Rend. Acad. Sci. URSS14(1937) 109

    work page 1937

  2. [2]

    Calculation of the Cherenkov light yield from electromagnetic cascades in ice with Geant4

    L. Radel and C. Wiebusch,Calculation of the Cherenkov light yield from electromagnetic cascades in ice with Geant4,Astropart. Phys.44(2013) 102 [1210.5140]

    work page internal anchor Pith review Pith/arXiv arXiv 2013

  3. [3]

    L. Rädel,Simulation studies of the cherenkov light yield from relativistic particles in high-energy neutrino telescopes with Geant4, master’s thesis, RWTH Aachen University, Aachen, Germany, 2012

    work page 2012

  4. [4]

    Ferrari, P.R

    A. Ferrari, P.R. Sala, A. Fasso and J. Ranft,FLUKA: A multi-particle transport code, CERN Yellow Reports: Monographs, CERN, Geneva (10, 2005), 10.5170/CERN-2005-010

    work page doi:10.5170/cern-2005-010 2005

  5. [5]

    Measurement of South Pole ice transparency with the IceCube LED calibration system

    F. Ballarini et al.,The FLUKA code: Overview and new developments,EPJ Nuclear Sci. Technol.10(2024) 16. [6]IceCubecollaboration,Measurement of South Pole ice transparency with the IceCube LED calibration system,Nucl. Instrum. Meth. A711(2013) 73 [1301.5361]. [7]IceCubecollaboration,Muon Production of Hadronic Particle Showers in Ice and Water, in Proceedi...

    work page internal anchor Pith review Pith/arXiv arXiv 2024

  6. [6]

    Cerutti, A

    F. Cerutti, A. Empl, A. Fedynitch, A. Ferrari, R. Garcia Alia, P.R. Sala et al.,Nuclear model developments in FLUKA for present and future applications,EPJ Web Conf.146(2017) 12005

    work page 2017

  7. [7]

    Eilers and B.D

    P.H.C. Eilers and B.D. Marx,Flexible smoothing with B-splines and penalties,Statistical Science11(1996) 89

    work page 1996

  8. [8]

    Currie, M

    I.D. Currie, M. Durban and P.H.C. Eilers,Generalized linear array models with applications to multidimensional smoothing,Journal of the Royal Statistical Society: Series B (Statistical Methodology)68(2006) 259 [https://rss.onlinelibrary.wiley.com/doi/pdf/10.1111/j.1467-9868.2006.00543.x]

    work page doi:10.1111/j.1467-9868.2006.00543.x 2006

  9. [9]

    Azzalini and A

    A. Azzalini and A. Capitanio,Statistical applications of the multivariate skew normal distribution,Journal of the Royal Statistical Society. Series B (Statistical Methodology)61 (1999) 579

    work page 1999

  10. [10]

    O.E. Barndorff-Nielsen,Normal inverse gaussian distributions and stochastic volatility modelling,Scandinavian Journal of Statistics24(1997) 1 [https://onlinelibrary.wiley.com/doi/pdf/10.1111/1467-9469.00045]

    work page doi:10.1111/1467-9469.00045 1997

  11. [11]

    Virtanen, R

    P. Virtanen, R. Gommers, T.E. Oliphant, M. Haberland, T. Reddy, D. Cournapeau et al., SciPy 1.0: Fundamental Algorithms for Scientific Computing in Python,Nature Methods17 (2020) 261

    work page 2020

  12. [12]

    Peyré and M

    G. Peyré and M. Cuturi,Computational optimal transport,Found. Trends Mach. Learn.11 (2019) 355–607

    work page 2019

  13. [13]

    An Introduction to PYTHIA 8.2

    T. Sjöstrand, S. Ask, J.R. Christiansen, R. Corke, N. Desai, P. Ilten et al.,An introduction to PYTHIA 8.2,Comput. Phys. Commun.191(2015) 159 [1410.3012]

    work page internal anchor Pith review Pith/arXiv arXiv 2015

  14. [14]

    Koehne, K

    J.H. Koehne, K. Frantzen, M. Schmitz, T. Fuchs, W. Rhode, D. Chirkin et al.,PROPOSAL: A tool for propagation of charged leptons,Comput. Phys. Commun.184(2013) 2070. – 21 –

    work page 2013

  15. [15]

    Hallen,On the Measurement of High-Energy Tau Neutrinos with IceCube, Ph.D

    P. Hallen,On the Measurement of High-Energy Tau Neutrinos with IceCube, Ph.D. thesis, Master’s thesis, RWTH Aachen, 2013. [18]IceCubecollaboration,Improved modeling of in-ice particle showers for IceCube event reconstruction,JINST19(2024) P06026 [2403.02470]

    work page arXiv 2013

  16. [16]

    ACEM: A Python package for modeling Cherenkov light yields from particle showers

    I. Crawshaw and T. Yuan, “ACEM: A Python package for modeling Cherenkov light yields from particle showers.”https://pypi.org/project/acem/, 2025

    work page 2025

  17. [17]

    Harris, K.J

    C.R. Harris, K.J. Millman, S.J. van der Walt, R. Gommers, P. Virtanen, D. Cournapeau et al.,Array programming with NumPy,Nature585(2020) 357

    work page 2020

  18. [18]

    Penalized Splines for Smooth Representation of High-dimensional Monte Carlo Datasets

    N. Whitehorn, J. van Santen and S. Lafebre,Penalized Splines for Smooth Representation of High-dimensional Monte Carlo Datasets,Comput. Phys. Commun.184(2013) 2214 [1301.2184]. – 22 –

    work page internal anchor Pith review Pith/arXiv arXiv 2013