Pith. sign in

REVIEW 3 major objections 4 minor 32 references

A novel approach for deducing the mass composition of cosmic rays from lateral densities of EAS particles

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

Pith's one-line read The paper claims that two pointwise shape parameters of electron and muon lateral density profiles—local age and segmented slope—show cosmic-ray composition changing gradually from proton- to iron-dominated across the knee.

desk verdict A useful extension of local-age methods to muon LDFs, but the exponent swap in Eq. (1) biases the central observable and must be fixed before the results can be trusted. read the letter →

arxiv 1908.04150 v1 pith:UPOSDHH6 submitted 2019-08-12 astro-ph.HE

classification astro-ph.HE PACS 96.50.sd95.75.z96.50.S
keywords cosmic-raycompositionextensiveairshowerslateraldensitydistributionlocalageparametersegmentedslopeKASCADEkneeofspectrumNKGfunction
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

This paper argues that the mass composition of cosmic rays around the knee of the energy spectrum can be read off from the radial shapes of the lateral density profiles of air-shower electrons and muons, without needing a full reconstruction of each shower. Using proton- and iron-initiated Monte Carlo showers together with the KASCADE experiment's public electron and muon density data, the authors propose two pointwise shape estimators—the local age parameter for the electron (NKG) profile and the segmented slope parameter for the muon (Greisen) profile—and show that their mean extrema (minimum local age near 44 m, maximum segmented slope near 71 m) separate light and heavy primaries. Comparing simulation to data, they report a gradual change from a proton-dominated to an iron-dominated composition as shower size and muon size increase across the knee. If this is right, the lateral shapes of particle densities provide a simple, model-light handle on cosmic-ray composition in an energy region where the transition from galactic to extragalactic sources is still debated.

What carries the argument

The central objects are the local age parameter (LAP) and segmented slope parameter (SSP), pointwise shape estimators computed from density ratios at two core distances. For any lateral distribution $f(r)$, equation (1) defines $$a_{\rm local}=\frac{\ln(F_{ij}X_{ij}^{\alpha_1}Y_{ij}^{\alpha_2})}{\ln(X_{ij}Y_{ij})}$$ with $F_{ij}=\rho(r_i)/\rho(r_j)$, $X_{ij}=r_i/r_j$, and $Y_{ij}=(r_i/R+1)/(r_j/R+1)$. For the NKG electron form $f(r)=C(r/R)^{s-2}(1+r/R)^{s-4.5}$ the authors take $\alpha_1=4.5$, $\alpha_2=2$ and call the result the LAP; for the Greisen-type muon form $f_\mu(r)=A(r/R_G)^{-\beta}(1+r/R_G)^{-2.5}$ they take $\alpha_1=0$, $\alpha_2=2.5$ and call it the SSP. These estimators are meant to convert local density shapes into the age or slope of the cascade at a given radius, and the analysis uses their extrema, the minimum LAP near 44 m and the maximum SSP near 71 m, as mass-sensitive observables.

What would settle it

Simulate showers with a known single primary mass and a known NKG age, then apply the paper's exact equations (1) and (3) to the simulated densities; if the output local age differs from the input age by the radius-dependent term $2.5/(1+2x)$, the KASCADE trends in figures 4–6 cannot be interpreted as measuring true shower age, and the composition conclusion would need to be re-derived with the corrected formula.

Watch

Extended reading notes

Core claim

The central claim is that the local age parameter and the segmented slope parameter, estimated directly from density ratios at neighboring core distances, are mass-sensitive observables and that their averages over the KASCADE data track a composition transition around the knee. The paper's Figure 4 shows that the mean minimum local ages of electrons and muons from observed events fall between the proton and iron simulation bands and move from the proton side toward the iron side with increasing shower size and muon size; Figures 5 and 6 make the same point for muon local age and for the maximum segmented slope. The authors interpret this as evidence that the cosmic-ray composition gradually shifts from lighter (proton-dominated) to heavier (iron-dominated) nuclei as primary energy increases through the knee, consistent with the earlier KASCADE result based on the electron-size versus muon-size correlation. They also report that the radial variation of the local age and segmented slope is approximately independent of shower size, which they call scaling behavior.

Load-bearing premise

The load-bearing step is equation (1)'s claim that $s_{\rm local}$ recovered from a two-point density ratio equals the true shower age for the standard electron profile; direct algebra with the paper's chosen exponents gives $s_{\rm local}=s+2.5/(1+2x)$ instead, so the reported 'minimum at 44 m' and the composition trend built on it may be measuring this offset.

Editorial extensions

If this is right

  • The mean minimum local age near 44 m and the mean maximum segmented slope near 71 m can serve as mass discriminators for experiments that measure lateral densities but not the depth of shower maximum.
  • Because the radial LAP and SSP curves are nearly size-independent, composition might be inferred from a single shape extremum per event, reducing the need for precise primary-energy reconstruction.
  • Extending the local-age analysis from electrons to muons gives a second, independent composition probe that is less affected by atmospheric attenuation.
  • If the reported composition trend is correct, the knee of the cosmic-ray spectrum is not a sudden switch but a gradual light-to-heavy transition, in line with rigidity-dependent cutoffs of galactic sources.

Reading between the lines

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

  • A direct algebra test of equation (1) with the paper's own NKG exponents gives $s_{\rm local}=s+2.5/(1+2x)$ for a pure NKG profile rather than the true age $s$, so the 'minimum at 44 m' may partly reflect this offset; a corrected estimator would need to be tested before the composition conclusion is taken quantitatively.
  • A stronger test would apply the same LAP and SSP recipe to Monte Carlo showers with known mixed compositions (helium, CNO, and silicon as well as proton and iron) and check whether the mean extrema still separate mass groups with the same ranking.
  • The 44 m and 71 m feature locations may depend on detector spacing and on the chosen Moliere and Greisen radii; repeating the analysis at a denser array would show whether the features are a property of the shower or of the sampling geometry.
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 / 4 minor

Summary. The paper introduces two local shape observables for extensive air showers: the local age parameter (LAP), computed from the NKG lateral density function, and the segmented slope parameter (SSP), computed from a Greisen-like muon lateral density function. Using CORSIKA proton and iron simulations at KASCADE energies, the authors study radial profiles of these parameters and define the mean minimum LAP and the mean maximum SSP as functions of electron shower size or truncated muon size. They then compute the same observables from public KASCADE (KCDC) data and interpret the comparison as evidence for a gradual light-to-heavy change of cosmic-ray mass composition across the knee. The paper also reports a first study of radial LAP/SSP profiles for muon lateral densities.

Significance. If the proposed shape parameters were cleanly defined, the approach could offer a useful composition-sensitive observable that uses only lateral density data, and the comparison with public KASCADE data would be a valuable external benchmark. The paper has concrete strengths: it uses a dedicated CORSIKA simulation sample, checks two high-energy hadronic interaction models, and is the first to report radial LAP/SSP profiles for muon densities. However, the central algebraic definition of the LAP is inconsistent with the NKG form as stated, so the method as written does not return the lateral shower age. Because the LAP is the main observable behind Figures 1, 3, 4, and 5 and the primary composition conclusion, this is a load-bearing technical error that must be corrected and the analysis redone before the claims can be evaluated.

major comments (3)
  1. [II, Eqs. (1)-(2)] The stated exponents α1=4.5 and α2=2 are inconsistent with the NKG lateral distribution in Eq. (3). Substituting f(r)=C x^(s-2)(1+x)^(s-4.5) into Eq. (1) gives ln(Fij Xij^{α1} Yij^{α2}) = (s-2+α1) ln Xij + (s-4.5+α2) ln Yij. For this to equal s (ln Xij + ln Yij), one needs α1=2 and α2=4.5, not the values printed in the text. With the printed values, the continuous limit in Eq. (2) yields slocal(r)=s+2.5/(1+2r/Rm), a radial-dependent offset. Quantitatively, at r=44 m with Rm=89 m the offset is about 1.26, and with Rm=320 or 420 m it is about 2.0. Thus every reported LAP value, the radial position and depth of the minimum, the Rm dependence shown in Figure 3, and the quantitative age scales in Figures 4 and 5 are biased. If the implementation actually used the swapped exponents, the text misstates the method; if it followed the text, the reported observables are not lateral shower ages. The paper must correct the definition, recompute all LAP-based quantities, and re-examine the composition conclusions.
  2. [V.A and VI.4] Section V.A first states that the KASCADE muon data 'also indicates a heavier domination with increasing muon size/energy across the knee,' but the immediately following sentence says that the idea of a gradual transition from lighter to heavier composition is 'somewhat at least not obvious in experimental results.' These statements are mutually contradictory. Given that the muon-LAP trend is acknowledged not to be obvious, the strong conclusion in Section VI.4 that Figures 4, 5, and 6 indicate a gradual light-to-heavy change is an overstatement, independent of the algebraic issue in Eq. (2). The authors should either provide a quantitative measure of the trend or temper the conclusion.
  3. [V, Figs. 4-6] The composition inference is only a qualitative bracketing of the KASCADE data between proton and iron simulation curves. No mixture fraction, goodness-of-fit statistic, or compatibility test between the data and the p/Fe templates is presented. Since the central claim is a 'gradual change' in composition, the paper should quantify how well the data points are described by mixtures of the two templates, with uncertainties, rather than relying on visual interpolation.
minor comments (4)
  1. [IV] The text contains the typo 'fist time' in the sentence introducing the first muon LAP/SSP study; it should be 'first time.'
  2. [Fig. 1 caption] The caption refers to 'Fig. c and Fig. d'; these should be labeled '(c)' and '(d)' consistently with the other panels.
  3. [Fig. 5 caption] The caption contains the typo 'QGSJeT'; it should be 'QGSJet' for consistency with the text.
  4. [II] The phrase 'each point' and the notation alocal(i,j) for distinct radii is a little confusing; the paper would benefit from an explicit statement that Eq. (2) is obtained in the limit ri -> rj and how the discrete experimental estimates are assigned to radial positions.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the composition claim rests on independent CORSIKA p/Fe templates compared against external KASCADE/KCDC data, not on self-referential fitting.

full rationale

The paper's central inference is a template comparison, not a circular reduction. The observables (LAP and SSP) are defined a priori by Eqs. (1), (3), and (6) from the NKG/Greisen LDFs, with KASCADE's external choices for Rm (89 m, 320/420 m), muon-size truncation, and LECFs. The same definition is applied to simulated CORSIKA proton/iron showers and to KASCADE/KCDC densities, and the mass-composition statement is obtained by comparing the KASCADE points to p/Fe curves in Figs. 4-6. No parameter is fitted to KASCADE data and then renamed a prediction, and no uniqueness theorem or forced ansatz is imported from the authors' prior work. Self-citations to Dey et al. [2,11] introduce the LAP concept and the scaling idea, but the scaling behavior is re-derived from the present simulations, and the composition claim is benchmarked against an independent experiment, so these self-citations are not load-bearing. Two non-circular concerns should be flagged. First, the algebra in Eqs. (1)-(3) as printed is suspect: substituting the NKG f(r) into Eq. (1) with alpha1=4.5, alpha2=2 gives alocal = s + 2.5/(1+2x), not s, so the 'local age' may be a distance-dependent transformation of the shower age rather than the age itself; this is a correctness/calibration issue, not a self-referential reduction. Second, the paper itself concedes in Sec. V.A that the muon-LAP trend 'is somewhat at least not obvious in experimental results', which weakens the strongest claim but is an evidentiary limitation, not circularity.

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

The central analysis rests on the LAP/SSP observables defined against NKG and Greisen functions. The most fragile input is the exponent assignment in the local age formula; algebra shows the stated assignment does not return the NKG age. The simulation and data also hinge on the chosen hadronic models, the scale radii, the radial cuts, and the two-component (p, Fe) interpretation of the composition.

free parameters (3)
  • Moliere radius for muon LDF = 320 m and 420 m
    Scale radius in the muon lateral distribution, adopted from KASCADE analysis conventions. The paper shows (Fig. 3) that the local age distribution changes with this radius, so the chosen value influences all LAP and SSP numbers, though relative comparisons between data and simulation share the same choice.
  • Observable positions = 44 m for mean minimum LAP, 71 m for mean maximum SSP
    The radial positions used to summarize each event were read from the simulated radial curves, so they are not fixed a priori by a physical principle; different positions would change the reported composition comparison.
  • Spectral index for MC generation = 2.7 below 3e15 eV, 3.1 above
    Used to weight the simulated event sample; not fitted to data, but it sets the energy distribution of events and therefore the statistical composition of the sample in each size bin.
assumptions (5)
  • domain assumption NKG and Greisen LDFs describe the lateral distributions of EAS electrons and muons
    Section II defines the local shape parameters against these analytical forms; if the true LDD does not follow these forms, the extracted LAP and SSP values are biased.
  • ad hoc to paper The local age formula (Eq. 1) with alpha1=4.5 and alpha2=2 returns the NKG age
    Direct algebra with the NKG form shows the stated exponents produce slocal = s + 2.5/(1+2x), not s. The correct assignment appears to be alpha1=2 and alpha2=4.5. This enters in Section II and all subsequent LAP estimates.
  • domain assumption Hadronic interaction models EPOS, QGSJet, and EPOS-LHC bracket real hadronic interactions
    Section III uses these models to generate simulations for proton and iron primaries; if the true interaction model yields different muon lateral distributions, the p/Fe templates may be biased.
  • domain assumption KASCADE muon data are reliable in the 10-120 m radial range after LECF conversion
    Section IV restricts the muon data to this range because data beyond 120 m are described as very unreliable; the analysis assumes the 10-120 m range is accurate and sufficient.
  • domain assumption Only proton and iron primaries bracket the full composition
    Section III generates only p and Fe showers; composition is inferred by comparing data to these two templates, ignoring helium, CNO, and intermediate nuclei, which may bias the inferred trend.

how reviews work

0 comments
Cite this review

Pith. "Pith review of A novel approach for deducing the mass composition of cosmic rays from lateral densities of EAS particles." pith.science (2026). https://pith.science/paper/UPOSDHH6

@misc{pith2026190804150,
  author       = {Pith},
  title        = {Pith review of: A novel approach for deducing the mass composition of cosmic rays from lateral densities of EAS particles},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/UPOSDHH6}},
  note         = {Machine review of arXiv:1908.04150}
}
read the original abstract

A Monte Carlo (MC) simulation study of cosmic ray (CR) extensive air showers (EAS) has been carried out in the energy regime of the KASCADE experiment. From the characteristics of lateral distributions of electrons and muons of simulated EAS, some important EAS observables are extracted by a novel approach, and their CR mass-sensitivity is demonstrated. The study takes into account the issue of the experimental lateral density profiles of EAS electrons and muons after introducing the notion of the local age and segmented slope parameters, aimed to extract information on CR mass composition from observed data. The estimated lateral shower age and slope from the analysis of the KASCADE data (KCDC) agrees with the idea of a gradual change of CR mass composition from light to heavy around the knee.

Figures

Figures reproduced from arXiv: 1908.04150 by the authors.

Figure 1
Figure 1. FIG. 1: Variation of the LAP (estimated from simulated elec [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 3
Figure 3. FIG. 3: Distribution of the mean minimum LAP from sim [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 2
Figure 2. FIG. 2: Variation of the SSP (estimated from simulated muon [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figures from the paper (2 more)
Figure 4
Figure 4. Figure 4: FIG. 4: Variation of the mean minimum local age (estimated [PITH_FULL_IMAGE:figures/full_fig_p006_4.png]
Figure 6
Figure 6. Figure 6: FIG. 6: Variation of the mean maximum segmented slope [PITH_FULL_IMAGE:figures/full_fig_p006_6.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

32 extracted references · 32 canonical work pages

  1. [1]

    This feature of the LAP still persists when the parameter has been estimated from the LDD of muons

    The LDD of electrons in a shower manifests some sort of scaling (shower size independent) nature in terms of the LAP. This feature of the LAP still persists when the parameter has been estimated from the LDD of muons. The radial variation of the LAP follows a config- uration where with an increasing of the core distance, the parameter decreases initially a...

  2. [2]

    Figure - 1(a) represents slocal(r) versus r variation ob- tained from the simulated LDD of electrons for different electron size intervals. To examine whether the experi- mental data also demonstrates a high-low-high kind of nature in the LAP, we have included the local ages from the LDD data of electrons of the KASCADE experiment for a particular electron...

  3. [3]

    The frequency distributions of the minimum values of slocal(r) for two different Moli´ ere radii possess different mean val- ues

    The scale distance such as the Moli´ ere/Greisen ra- dius used in the LDF functions may lead to determine the shape of the slocal(r) or βss(r) versus r curves. The frequency distributions of the minimum values of slocal(r) for two different Moli´ ere radii possess different mean val- ues. This suggests that the scale distance regulates the shape of the sloc...

  4. [4]

    The characteristic feature of SSP versus the core distance curves satisfy a complete opposite configuration compared to the radial variation of the LAP

    The LDD of muons in a shower also exhibits scaling nature in terms of the SSP. The characteristic feature of SSP versus the core distance curves satisfy a complete opposite configuration compared to the radial variation of the LAP. The SSP estimated from the LDD of muons follows a low-high-low kind of radial variation at least within the range 10 − 300 m. ...

  5. [5]

    Haungs et al., European Phys.J

    A. Haungs et al., European Phys.J. C, 78:741 (2018)

  6. [6]

    This above feature sup- ports the results obtained from the study of the Ne–Nµ variations in the KASCADE experiment

    The KASCADE experimental data in figures 4, 5 and 6 in terms of the important measured parameters (mean minimum LAP and mean maximum SSP) indicate that the CR mass composition follows a gradual change from predominantly lighter (proton-dominated) to heav- ier (iron-dominated) nuclei with the increase of shower size and muon size or energy. This above featu...

  7. [7]

    J. N. Capdevielle and J. Procureur, Proc. 18th ICRC, Bangalore, 11 307 (1983)

  8. [8]

    R. K. Dey, A. Bhadra and J. N. Capdevielle, J. Phys. G:Nucl. Part. Phys. 39 085201 (2012)

Show all 32 references
  1. [9]

    J. N. Capdevielle and J. Gawin, J. of Physics G, 8, 1317 (1982)

  2. [10]

    D. Heck, J. Knapp, J. N. Capdevielle, G. Schatz and T. Thouw, FZKA report-6019 ed. FZK The CORSIKA Air Shower Simulation Program , Karlsruhe (1998)

  3. [11]

    R. K. Dey and S. Dam, Eur. Phys.J. Plus 131: 402 (2016)

  4. [12]

    Moli‘ere,: Cosmic Radiation, W

    G. Moli‘ere,: Cosmic Radiation, W. Heisenberg, ed., Dover Publications, New York, 1st ed., (1946)

  5. [13]

    Rossi and K Greisen, Rev

    B. Rossi and K Greisen, Rev. Mod. Phys. 13, 240 (1941)

  6. [14]

    Nishimura, and K

    J. Nishimura, and K. Kamata: Progr. Theor. Phys., 5, 899 (1950)

  7. [15]

    Kamata and J

    K. Kamata and J. Nishimura, Prog. Theor. Phys. Suppl. 6, 93 (1958)

  8. [16]

    Lipari, Phys

    P. Lipari, Phys. Rev. D 79 063001 (2009)

  9. [18]

    J. N. Capdevielle, J. Cohen, J. Phys. G: Nucl. Part. Phys . 31 507 (2005)

  10. [19]

    Bourdeau, J.N

    M.F. Bourdeau, J.N. Capdevielle and J. Procureur, J. Phys.G, 6, 901, (1980)

  11. [20]

    Two smaller samples of the MC events have also been generated with combinations, QGSJet01.c [22] - UrQMD and EPOS- LHC [23]-UrQMD

    interaction model for the high-energy hadronic in- teractions, in combination with the UrQMD model [21] for the low-energy hadronic interactions. Two smaller samples of the MC events have also been generated with combinations, QGSJet01.c [22] - UrQMD and EPOS- LHC [23]-UrQMD. ...

  12. [21]

    Greisen, Prog

    K. Greisen, Prog. in Cosmic Ray Physics, North Holland, Co., Amsterdam, 3 1 (1956)

  13. [22]

    Capdevielle et al., J

    J.N. Capdevielle et al., J. Phys. G: Nucl. Part. Phys. 31 507 (2005)

  14. [23]

    Nagano et al., J

    M. Nagano et al., J. Phys. Soc. Japan 53 (1984) 1667

  15. [24]

    J. N. Capdevielle and P. Gabinski, J. Phys. G.:Nucl. Par t. Phys. 16, 769 (1990)

  16. [25]

    S. K. Gupta et al., Nucl. Instr. Meth. A, 540, 311 (2005)

  17. [26]

    Werner, et al., Phys

    K. Werner, et al., Phys. Rev. C 74, 044902 (2006)

  18. [27]

    Bleicher, et al., J

    M. Bleicher, et al., J. Phys. G: Nucl. Part. Phys. 25, 185 9 (1999)

  19. [28]

    N. N. Kalmykov, S. S. Ostapchenko and A. I. Pavlov Nucl. Phys. B (Proc. Suppl.) 52 17 (1997)

  20. [29]

    L. G. Dedenko et al., J. Phys. Conf. Series, 934 012017 (2017)

  21. [30]

    W. R. Nelson, H. Hiramaya, D. W. O. Rogers, Report 8 SLAC 265, (1985)

  22. [31]

    Apel et al., Nucl

    W.D. Apel et al., Nucl. Instrum. Methods, A 620 (23):202, (2010)

  23. [32]

    Sanyal et al., Aust

    S. Sanyal et al., Aust. J. Phys. 46 589 (1993)

  24. [33]

    Antoni et al., Astropart

    T. Antoni et al., Astropart. Phys. 14 245 (2001)

Pith tools

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