Pith. sign in

REVIEW 3 major objections 3 minor 23 references

Reduction of the Uncertainty in the Atmospheric Neutrino Flux Prediction Below 1 GeV Using Accurately Measured Atmospheric Muon Flux

T0 review · 3 major / 3 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Accurate muon flux measurements, especially from a high mountain, can constrain the hadronic-interaction uncertainty in sub-GeV atmospheric neutrino flux predictions to below about 10 percent.

desk verdict A useful, clearly-presented framework for constraining sub-GeV atmospheric neutrino flux uncertainties with muon measurements, but the headline <10% bounds are conditional on an uncalibrated 50%-per-grid-point variation ensemble. read the letter →

arxiv 1908.08765 v4 pith:3JMZKFHK submitted 2019-08-23 astro-ph.HE

M. Honda , M. Sajjad Athar , T. Kajita , K. Kasahara , S. Midorikawa This is my paper · ORCID
classification astro-ph.HE PACS 95.85.Ry13.85.Tp13.35.Bv14.60E
keywords atmosphericneutrinofluxmuonhadronicinteractionuncertaintyintegralkernelB-splinevariationsub-GeVneutrinoscalibrationoscillation
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

Sub-GeV atmospheric neutrino fluxes are needed for precision neutrino-oscillation experiments, but hadronic-interaction model error dominates their uncertainty. This paper establishes that a precision measurement of the atmospheric muon flux can quantify and, at the right observation site, largely remove that error: if a calculated model can reconstruct the measured muon flux, the same model's neutrino flux is correspondingly close to the true one. Using millions of random deformations of the hadronic production kernel, the authors show the residual neutrino-flux uncertainty can be brought below about 10% ($\varsigma_0 < 0.1$) across $0.15 < E_\nu < 10$ GeV vertically and $E_\nu > 0.3$ GeV horizontally by muon data from a high mountain near 4500 m, whereas sea-level or balloon-altitude data leave larger irreducible errors below 1 GeV. The paper also identifies hadronic scattering-angle uncertainty as a separate source that muon flux normalization alone cannot remove, and shows it can be attacked through the zenith-angle dependence of high-altitude muon fluxes.

What carries the argument

The load-bearing object is the integral kernel $D(N_{\mathrm{proj}}, p_N, M_{\mathrm{born}}, p_M, L, p_L, x)$ of the pseudo-analytic lepton-flux formula, which bundles the meson decay probability, hadronic production probability, production cross section, air density, and projectile flux into a density distribution over the hadronic-interaction phase space $(\log_{10} p_{\mathrm{proj}}, \log_{10} p_{\mathrm{born}})$. Variations of the interaction model are generated by multiplying this kernel by $1 + \delta \sum_{ij} R_{ij} B_i(\log_{10} p_{\mathrm{proj}}) B_j(\log_{10} p_{\mathrm{born}})$, with $R_{ij}$ independent standard normal random numbers and third-order B-spline basis functions; with $\delta=1$ this gives a 50% standard deviation per grid point. The method then generates millions of random kernel variations, computes the induced neutrino and muon flux variations, and selects only those whose muon flux variation satisfies $|\Delta \Phi_\mu/\Phi_\mu| < \varepsilon$. The standard deviation of the constrained neutrino variation distribution, $\sigma_{\mathrm{shrink}}$, is fitted by $\sqrt{\varsigma_0^2 + (\varsigma_1 \varepsilon)^2}$, so $\varsigma_0$ is the neutrino-flux uncertainty irreducible by muon data. This machinery is what converts a muon reconstruction residual into a quantitative neutrino flux error.

What would settle it

Compute $\varsigma_0$ with correlated or larger random deformations of the kernel (for example, smoothing neighboring B-spline coefficients or doubling $\delta$) and check whether the high-mountain result still gives $\varsigma_0 < 0.1$ in $0.15 < E_\nu < 10$ GeV; if not, the claimed reduction is an artifact of the independence assumption. A complementary check is to calibrate two independent hadronic interaction models against the same 5%-accurate mountain muon data and see whether their sub-GeV neutrino fluxes agree to better than 10%.

Watch

Extended reading notes

Core claim

The paper's central claim is that the uncertainty in the predicted atmospheric neutrino flux from hadronic-interaction model error can be read off quantitatively from how well the same model reproduces an accurately measured atmospheric muon flux: when the calculated muon flux can be made to match a precision measurement, the calculated neutrino flux is forced close to the true one. Concretely, using variations of the hadronic interaction model built from random B-spline deformations of the production kernel, the authors find that the spread of neutrino flux variations shrinks as the muon flux variation is constrained, and the irreducible part of that spread is the muon-independent uncertainty $\varsigma_0$. The value of $\varsigma_0$ depends strongly on the muon observation site: at sea level it rises steeply below 1 GeV, but for a high-mountain site at 4500 m above sea level, with muons measured down to $P_\mu \geq 0.3$ GeV/c, it stays below 0.1 for all neutrino species in $0.15 < E_\nu < 10$ GeV vertically and $E_\nu > 0.3$ GeV horizontally. The paper further shows that the hadronic scattering-angle uncertainty is a separate error source, producing roughly a 10% neutrino flux change for a 20% scattering-angle change near the horizontal, and that high-altitude muon zenith-angle data can help constrain it. The concluding claim is therefore that reconstructing an accurately measured atmospheric muon flux is sufficient to calculate the atmospheric neutrino flux accurately in the sub-GeV region.

Load-bearing premise

The quantitative bounds depend on the assumption that the true hadronic-interaction error can be represented by random, independent, 50%-per-grid-point deformations of the production kernel and that these deformations do not change the incoming cosmic-ray flux; if real model errors are correlated, larger, or shift the projectile flux, the reported error reductions do not follow.

Editorial extensions

If this is right

  • A precision muon measurement on a 3000–5000 m mountain, covering momenta down to 0.3 GeV/c, would keep the residual hadronic-interaction error of atmospheric neutrino flux below 10% for all neutrino species across $0.15 < E_\nu < 10$ GeV vertically and $E_\nu > 0.3$ GeV horizontally.
  • Sea-level muon data, even extended down to 0.3 GeV/c, leave an irreducible neutrino flux uncertainty that rises steeply below 1 GeV, so they cannot fill the same calibration role.
  • Balloon-altitude muon data produce larger irreducible neutrino uncertainty than mountain data, making balloon programs a less promising route for this calibration.
  • Reaching 5% neutrino flux accuracy requires reducing hadronic scattering-angle uncertainty to about 10%, since a 20% scattering-angle change induces roughly 10% neutrino flux change near the horizontal; zenith-angle-dependent muon data from a high mountain can contribute to that reduction.
  • The framework turns muon-flux reconstruction residuals directly into a neutrino-flux error estimate, replacing the earlier qualitative muon-calibration argument with a quantitative one.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The same random-deformation machinery could be rerun with a prior whose amplitude and correlation are tuned to accelerator data on specific production channels; the paper notes that doing so for kaons would suppress the rise in $\varsigma_0$ at higher energies, but leaves that rerun to future work.
  • The site dependence implies a design rule the authors do not state: the best muon calibration site is one where muons suffer little energy loss, so their parent-meson phase space overlaps that of sub-GeV neutrinos, which selects mountains near 3–5 km over both sea level and balloon altitude.
  • A natural next step, using the paper's own ingredients, is a joint fit of high-altitude muon flux normalization and zenith-angle distribution to separate kernel-magnitude errors (tracked by $\varsigma_0$) from scattering-angle smearing (tracked by the zenith-angle response).
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 3 minor

Summary. The paper proposes a method to quantify and reduce the uncertainty in sub-GeV atmospheric neutrino flux predictions by using precision measurements of the atmospheric muon flux. The authors write a pseudo-analytic expression for the atmospheric lepton flux, identify an ``integral kernel'' in the hadronic-interaction phase space, and model model-to-real differences as random independent Gaussian B-spline deformations of that kernel with a 50% per-grid-point amplitude (δ=1, Δ=0.5 in Eqs. 8-13). They generate 3,000,000 random model variations, impose a constraint on the resulting muon-flux variation (Eq. 16), and measure how much the neutrino-flux variation shrinks. The shrinkage is characterized by σ_shrink and fitted to σ_shrink = sqrt(ς0^2 + (ς1 ε)^2) (Eq. 20), with ς0 interpreted as the muon-independent component of the neutrino-flux uncertainty. They compute ς0 for neutrino energies from 0.1 to 100 GeV at Kamioka using muon kernels at Tsukuba, Mt. Norikura, Hanle, and balloon altitude. The main quantitative result is that for muons observed at Hanle (4500 m) down to 0.3 GeV/c, ς0 < 0.1 for all neutrino kinds in 0.15-10 GeV (vertical) and above 0.3 GeV (horizontal) directions. The paper also compares with existing BESS muon data and discusses the influence of projectile-flux ratio and scattering-angle uncertainties.

Significance. If the quantitative claim is taken at face value, the paper offers a practical path to controlling hadronic-interaction uncertainties in atmospheric neutrino flux calculations below 1 GeV using high-mountain precision muon measurements, which is directly relevant to upcoming long-baseline and atmospheric neutrino experiments. The pseudo-analytic integral-kernel framework is a useful conceptual contribution, and the systematic survey of observation sites is informative. The paper is also honestly self-critical: it explicitly reports that the BESS muon data below 1 GeV cannot be simultaneously reproduced at Tsukuba and Mt. Norikura, and it discusses scattering-angle uncertainty as a residual limitation. However, the absolute scale of the reported uncertainties (ς0 values) is set by the uncalibrated assumption of independent 50%-per-grid-point variations, so the central numerical conclusions are conditional on the assumed ensemble rather than directly validated by data.

major comments (3)
  1. [Sec. III, Eqs. (8)-(13), and Figs. 7-13] The amplitude δ=1 and the independence of the normal random coefficients R_ij are assumed without calibration. Since ΔΦ_ν and ΔΦ_μ both scale linearly with δ, the fitted values of ς0 are proportional to the assumed 50% standard deviation per B-spline grid point. The paper provides no external check of this amplitude or of the independence assumption, for example by comparing the ensemble with the spread among existing hadronic interaction models or by using the BESS data to constrain δ. As a result, the absolute statement ς0<0.1 in Fig. 12 is a property of the assumed variation ensemble, not a measured atmospheric-neutrino uncertainty. The site ranking may be more robust than the absolute values, but the paper should either calibrate the ensemble or explicitly present all quantitative bounds as conditional on the assumed variation model.
  2. [Sec. V and Sec. VIII] The paper states that below 1 GeV the model fails to reproduce the BESS muon fluxes at Tsukuba and Mt. Norikura simultaneously. This is exactly the energy range most relevant to the paper's central claim. The variation study in Sec. IV conditions on being able to reconstruct the muon flux to within ε≲0.05, but the actual model cannot do that below 1 GeV. Thus the concluding statement that ``when we can reconstruct the atmospheric muon flux measured by a precision experiment, we can also calculate the atmospheric neutrino flux accurately'' is not demonstrated for the low-energy region where it matters most. The paper should either resolve the low-energy BESS discrepancy, for example by including it as a constraint on the variation ensemble, or explicitly limit the conclusions to models that already pass the low-energy muon data.
  3. [Sec. IV, Eqs. (17)-(19)] Eqs. (17)-(19) are internally inconsistent as written. Eq. 17 is not a normal distribution because the exponent lacks a minus sign and the factor 1/2. Eq. 18 then states that the normalized ratio of constrained to unconstrained distributions is exp[(1/σ_ε^2 - 1/σ_∞^2) x^2], which grows with |x| rather than decaying, so the concentration effect is lost. If corrected to a proper Gaussian ratio, the exponent becomes -(x^2/2)(1/σ_ε^2 - 1/σ_∞^2), which is what the subsequent definition of σ_shrink in Eq. 19 requires. The authors should re-derive these equations carefully and ensure the notation in Eq. 17 matches a normalized Gaussian.
minor comments (3)
  1. [General] There are numerous typographical errors: ``neuron'' should be ``neutrino'' after Eq. (14), ``path'' should be ``pass'' below Eq. (17), ``M. Norikura'' should be ``Mt. Norikura'' in Fig. 8, and ``tat neutrino energy'' in the Fig. 3 caption should be ``at neutrino energy''.
  2. [Fig. 15 caption] The right-panel caption repeats ``atmospheric neutrino at Eν = 0.1 GeV and atmospheric muon at Pµ = 0.1 GeV/c''; based on the text in Sec. VII it should be Eν = 0.3 GeV and Pµ = 0.3 GeV/c.
  3. [Sec. IV] The criterion for choosing the muon momentum range (correlation larger than 90% of its maximum) is introduced without justification or sensitivity checks. A brief comment on how the results depend on this threshold would help.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the quantitative bounds are conditional variance decompositions under an explicitly stated random-variation ensemble, with external BESS checks and no load-bearing self-citation.

full rationale

The central quantity ς0 is not fitted to the target neutrino data and is not defined in terms of the muon measurements it is meant to constrain. It is obtained from an explicit variance decomposition of the model's own response to random B-spline deformations (Eqs. 8–13): σ_shrink is the conditional spread of ΔΦ_ν/Φ_ν given |ΔΦ_μ/Φ_μ| < ε (Eqs. 16–19), and Eq. 20 is a two-parameter fit to those generated distributions. The paper's conclusion is explicitly conditional: 'By considering our interaction model ... is a variation of the ideal interaction model which can predict the true atmospheric neutrino and muon fluxes, we may conclude that when we can reconstruct the atmospheric muon flux measured by a precision experiment, we can also calculate the atmospheric neutrino flux accurately.' This is a transparent assumption, not a concealed circular reduction. The paper also provides an external check in Sec. V with BESS data, candidly reporting that below 1 GeV the model fails to reconstruct the Tsukuba and Mt. Norikura fluxes simultaneously, and in Sec. VIII it acknowledges 'the interaction model we are using is basically constructed using the accelerator data' and 'We have assumed 50% uncertainty in the integral kernel density at each grid point.' The self-citations (Refs. [9,15,16]) are to the group's prior simulation and calibration work, but the present BESS comparison re-derives the relevant conclusion rather than importing it as an unverified uniqueness result. Thus the derivation chain is self-contained: the quantitative bounds are conditional statements about the assumed random-variation ensemble, not circular equivalents of the input data.

Assumptions & free parameters 3 free parameters · 5 assumptions · 0 invented entities

The paper's central numbers come from injecting artificial variations into its own hadronic model. The list above captures the quantities chosen by hand and the background assumptions that the quantitative conclusion inherits. Free parameters: δ, Δ, and the 50% kaon uncertainty. Axioms: factorization of the flux, first-order response, representativeness of random variations, projectile flux insensitivity, and single-nucleon superposition. No new physical entities are introduced.

free parameters (3)
  • δ (variation amplitude) = 1 (chosen)
    Sets the magnitude of the random hadronic-model variation in Eq. 8; via Eq. 13 it sets a 50% per-grid-point standard deviation. All reported ς0 and σ_shrink values scale with this choice.
  • B-spline knot separation Δ = 0.5 (chosen)
    Chosen in log10(p) for both projectile and meson momenta (Sec. III). Determines the momentum-scale correlation length of the model variations and affects how effectively muon constraints propagate to neutrino flux.
  • Kaon production uncertainty per grid point = 50% (assumed)
    Explicitly assumed in Sec. IV ('we have also assumed the uncertainty of kaon production is 50 % at the every grid point'). Inflates ς0 at high energies for νμ, where kaon contribution is largest.
assumptions (5)
  • domain assumption Atmospheric lepton flux factorizes as a convolution of the integral kernel D with the hadronic interaction probability H and projectile flux (Eq. 1).
    Standard cascade equation structure; stated in Sec. II and used to build the variation formalism.
  • domain assumption The projectile flux Φ_proj is not largely affected by variation of hadronic interactions (Sec. III).
    Stated in Sec. III and revisited in Sec. VII with a ±10% Np/Nn test. The test only partially validates the assumption.
  • ad hoc to paper True hadronic-model uncertainty is representable as independent normal-random B-spline deformations of the integral kernel with δ=1 (Eqs. 8-13).
    No empirical justification is given for the amplitude or independence structure; this assumption sets the absolute scale of all quantitative uncertainty estimates.
  • domain assumption Linear response (first order in δ) is sufficient for flux variations (Eq. 5).
    The variation of the interaction model enters linearly in Eq. 5 and Eq. 10; higher-order terms are neglected without a check.
  • standard math Nucleus-air interactions can be represented as superposition of single nucleon interactions.
    Stated in Sec. II as the standard treatment for nucleus projectiles in the simulation.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Reduction of the Uncertainty in the Atmospheric Neutrino Flux Prediction Below 1 GeV Using Accurately Measured Atmospheric Muon Flux." pith.science (2026). https://pith.science/paper/3JMZKFHK

@misc{pith2026190808765,
  author       = {Pith},
  title        = {Pith review of: Reduction of the Uncertainty in the Atmospheric Neutrino Flux Prediction Below 1 GeV Using Accurately Measured Atmospheric Muon Flux},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/3JMZKFHK}},
  note         = {Machine review of arXiv:1908.08765}
}
abstract

We examine the uncertainty of the calculation of the atmospheric neutrino flux and present a way to reduce it using accurately measured atmospheric muon flux. Considering the difference of the hadronic interaction model and the real one as a variation of hadronic interaction, we find a quantitative estimation method for the error of the atmospheric neutrino flux calculation from the residual of the reconstruction of the atmospheric muon flux observed in a precision experiment, by the study of atmospheric neutrino and muon fluxes response to the variation of hadronic interaction. However, the efficiencty of this method is largely dependent on the observation site of the atmospheric muon flux, as the relation of the error of the atmospheric neutrino flux calculation and the residual of the reconstruction of the atmospheric muon flux is also largely dependent on the muon observation site, especially for the low energy neutrinos. We calculate several observation sites, near Kamioka at sea level, same but 2770m a.s.l.., Hanle India (4500m a.s.l.), and at Balloon altitude ($\sim$ 32km). Then we estimate how stringently can the atmospheric muon reduce the error in the calculation of the atmospheric neutrino flux. We also discuss on the source of error which is difficult to reduce by only the observation of atmospheric muon.

Figures

Figures reproduced from arXiv: 1908.08765 by the authors.

Figure 1
Figure 1. FIG. 1: The scatter plot of the projectile and meson momenta, [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2: The scatter plot of the projectile and meson momenta, [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3: Correlation coefficient for each combination of neutr [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figures from the paper (12 more)
Figure 4
Figure 4. Figure 4: FIG. 4: The muon momenta at which the muon show the maximum cor [PITH_FULL_IMAGE:figures/full_fig_p007_4.png]
Figure 5
Figure 5. Figure 5: FIG. 5: The ∆Φ [PITH_FULL_IMAGE:figures/full_fig_p007_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6: The [PITH_FULL_IMAGE:figures/full_fig_p008_6.png]
Figure 7
Figure 7. Figure 7: FIG. 7: The [PITH_FULL_IMAGE:figures/full_fig_p009_7.png]
Figure 8
Figure 8. Figure 8: FIG. 8: Left panel: vertically downward atmospheric muon flu [PITH_FULL_IMAGE:figures/full_fig_p010_8.png]
Figure 9
Figure 9. Figure 9: FIG. 9: The [PITH_FULL_IMAGE:figures/full_fig_p010_9.png]
Figure 10
Figure 10. Figure 10: FIG. 10: The [PITH_FULL_IMAGE:figures/full_fig_p011_10.png]
Figure 11
Figure 11. Figure 11: FIG. 11: The [PITH_FULL_IMAGE:figures/full_fig_p011_11.png]
Figure 12
Figure 12. Figure 12: FIG. 12: The [PITH_FULL_IMAGE:figures/full_fig_p012_12.png]
Figure 13
Figure 13. Figure 13: FIG. 13: The [PITH_FULL_IMAGE:figures/full_fig_p012_13.png]
Figure 14
Figure 14. Figure 14: FIG. 14: Left panel; relative ratio of projectile particle r [PITH_FULL_IMAGE:figures/full_fig_p013_14.png]
Figure 15
Figure 15. Figure 15: FIG. 15: The flux ratios of atmospheric neutrino and muon calc [PITH_FULL_IMAGE:figures/full_fig_p014_15.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

23 extracted references · 16 canonical work pages

  1. [1]

    Ahmed et al

    S. Ahmed et al. (ICAL), Pramana 88, 79 (2017), 1505.07380

    arXiv 2017

  2. [2]

    M. G. Aartsen et al. (IceCube), J. Phys. G44, 054006 (2017), 1607.02671

    work page Pith review arXiv 2017

  3. [3]

    Abe et al

    K. Abe et al. (Hyper-Kamiokande) (2018), 1805.04163

    arXiv 2018

  4. [4]

    Acciarri et al

    R. Acciarri et al. (DUNE) (2016), 1601.05471

    arXiv 2016

  5. [5]

    Aguilar et al

    M. Aguilar et al. (AMS Collaboration), Phys. Rev. Lett. 114, 171103 (2015), URL http://link.aps.org/doi/10.1103/PhysRevLett.114.171103

    doi:10.1103/physrevlett.114.171103 2015

  6. [6]

    Aguilar et al

    M. Aguilar et al. (AMS), Phys. Rev. Lett. 115, 211101 (2015)

    work page 2015

  7. [7]

    Adriani et al., The Astrophysical Journal 765, 91 (2013), URL http://stacks.iop.org/0004-637X/765/i=2/a=91

    O. Adriani et al., The Astrophysical Journal 765, 91 (2013), URL http://stacks.iop.org/0004-637X/765/i=2/a=91

    work page 2013

  8. [8]

    Abe et al., Adv

    K. Abe et al., Adv. Space Res. 60, 806 (2017)

    work page 2017

Show all 23 references
  1. [9]

    Sanuki, M

    T. Sanuki, M. Honda, T. Kajita, K. Kasahara, and S. Midori kawa, Phys. Rev. D75, 043005 (2007), astro-ph/0611201

    2007 arXiv

  2. [10]

    D. H. Perkins, Astropart. Phys. 2, 249 (1994)

    1994

  3. [11]

    D. H. Perkins, Nucl. Phys. B399, 3 (1993)

    1993

  4. [12]

    Fedynitch, J

    A. Fedynitch, J. Becker Tjus, and P. Desiati, Phys. Rev. D86, 114024 (2012), 1206.6710

    2012 arXiv

  5. [13]

    Y´ a˜ nez, A

    J.-P. Y´ a˜ nez, A. Fedynitch, and T. Montgomery, PoSICRC2019, 881 (2019), 1909.08365

    2019 arXiv

  6. [14]

    See https://www.iiap.res.in/centers/iao

  7. [15]

    Honda, T

    M. Honda, T. Kajita, K. Kasahara, and S. Midorikawa, Phy s. Rev. D 83, 123001 (2011), astro-ph/1102.2688, URL http://link.aps.org/doi/10.1103/PhysRevD.83.123001

    2011 arXiv

  8. [16]

    Honda, M

    M. Honda, M. Sajjad Athar, T. Kajita, K. Kasahara, and S. Midorikawa, Phys. Rev. D92, 023004 (2015), 1502.03916

    2015 arXiv

  9. [17]

    Haino et al

    S. Haino et al. (BESS), Phys. Lett. B594, 35 (2004)

    2004

  10. [18]

    Sanuki et al., Phys

    T. Sanuki et al., Phys. Lett. B541, 234 (2002), see also erratum [23]. 16

    2002

  11. [19]

    Achard et al

    P. Achard et al. (L3), Phys. Lett. B598, 15 (2004)

    2004

  12. [20]

    Abe et al

    K. Abe et al. (BESS), Phys. Lett. B564, 8 (2003)

    2003

  13. [21]

    Honda, T

    M. Honda, T. Kajita, K. Kasahara, and S. Midorikawa, Phy s. Rev. D52, 4985 (1995), hep-ph/9503439

    1995 arXiv

  14. [22]

    Battistoni, A

    G. Battistoni, A. Ferrari, T. Montaruli, and P. R. Sala, Astropart. Phys. 19, 269 (2003)

    2003

  15. [23]

    Sanuki et al., Phys

    T. Sanuki et al., Phys. Lett. B581, 272 (2004)

    2004

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.