{"id":"a9dccf36-bbd7-4404-a104-5a5e18a93560","arxiv_id":"1908.07930","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Using 22 million HAWC air showers from June 2016, the authors show that a modified NKG function fits the lateral particle distribution well and that the lateral age parameter separates proton and iron simulations with a figure of merit above 1 above about 6 TeV.","lead":"HAWC, a water Cherenkov observatory in Mexico, measured how particles from cosmic-ray air showers spread out as they hit the ground. The study finds which mathematical form best describes that spread and shows that a derived 'age' parameter could help tell light from heavy cosmic rays.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The good-description claim is not supported because every chi2/NDOF in Fig. 2 relies on the uncalibrated error model Eq. (3.1), and Sec. 7 itself admits no LDF describes the data satisfactorily.","rationale":"In good faith, the paper is a conference contribution whose central argument requires (i) a reliable metric for comparing four LDFs, and (ii) an adequate absolute fit quality for the winning LDF, before the lateral age can be promoted as a composition-sensitive observable. The reader's weakest assumption correctly targets the unvalidated error model Eq. (3.1), and I agree that this is the most load-bearing technical soft spot: it enters every chi2 value in the comparison. However, I would add that the paper's own wording creates an internal contradiction—'none of the selected LDFs describes satisfactorily' versus 'HAWC's LDF gives a good description'—and that this contradiction is not resolved by any stated goodness-of-fit threshold. Both problems are addressable rather than fatal: Eq. (3.1) can be calibrated against HAWC PMT data, and the goodness-of-fit claim can be restated as a relative, not absolute, ranking. The FOM mass-composition statement is model-dependent (QGSJET-II-03 only) and inherits the same fit-quality issues, so it should be framed as a preliminary MC-based sensitivity estimate. No evidence of fabrication or misconduct; the standard LDFs, data-driven binning, and clearly described fits are positive features. Thus the appropriate verdict remains CONDITIONAL, requiring the error-model calibration and goodness-of-fit clarification before the headline claim is accepted as stated.","tokens_in":5903,"tokens_out":6854,"duration_ms":67008,"concrete_test":"Recompute the fits and Fig. 2 (left) using two alternative error assignments: (a) the statistical error on the mean of each 2-m radial bin (with Poisson/photo-electron variance), and (b) an independently calibrated HAWC PMT charge-resolution function obtained from pedestal and single-PMT calibration runs or from the measured scatter of PMT charges within radial bins. If the rank order of Eqs. (3.2), (3.3), (3.5), (3.7), or the statement that Eq. (3.2) is among the best, changes in any energy bin, the central claim must be weakened. In addition, report a goodness-of-fit p-value or a pre-specified chi2/NDOF threshold; if Eq. (3.2) is rejected by that threshold, the phrase 'good description' in Section 7 is unsupported and should be replaced by 'least discrepant among tested functions.'","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that HAWC's LDF (Eq. 3.2) gives a good description of cosmic-ray induced EAS and that its lateral age parameter is mass-composition sensitive—depends on the chi2/NDOF comparison of Fig. 2 (left). That comparison is built on the empirical PMT-charge error model of Eq. (3.1): log10(Qerror) = 0.3 - 0.06667*log10(Qeff) for log10(Qeff) <= 3 and 0.1 otherwise. No calibration of this function against single-PMT charge resolution, photo-electron statistics, or radial-bin scatter is presented; the coefficients are simply asserted. Because the weights enter every fit, an unvalidated error model can change the ranking of Eqs. (3.2), (3.3), (3.5), and (3.7), so the conclusion that HAWC's LDF is a good description is not yet established. This is compounded by an internal tension in the paper: Section 7 states that none of the selected LDFs describes satisfactorily the measured data for all radial ranges and energies, immediately before claiming that HAWC's LDF gives a good description in 10^3.5-10^5.5 GeV. Without a pre-specified acceptance threshold, these two statements cannot both be supported by the same chi2 values. The FOM claim inherits this fragility because the lateral age s is obtained from these same fits and the separation is tested with only one hadronic interaction model (QGSJET-II-03).","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the lateral distribution function (LDF) of cosmic-ray induced air showers recorded by the HAWC observatory, using 1.3 days of 2016 data with reconstructed energies between 10^3.5 and 10^5.5 GeV and zenith angles below about 16.7 degrees. Four LDF parameterizations are fitted to the mean effective charge in radial bins via chi-square minimization: the HAWC/NKG-like function of Eq. (3.2), the KASCADE function of Eq. (3.3), the ARGO function of Eq. (3.5), and a modified scaling formalism of Eq. (3.7). The authors compare the resulting chi-square per degree of freedom as a function of energy, conclude that HAWC's LDF provides the best description of the data among the tested forms, and use the lateral age parameter from that fit to assess proton/iron separation through a figure of merit. The paper also compares the measured average lateral distributions with CORSIKA/QGSJET-II-03 predictions for proton and iron primaries.","tokens_in":6199,"tokens_out":3614,"duration_ms":35709,"significance":"If the central claims were fully supported, the paper would establish a useful result: the lateral age parameter derived from a single LDF fit could serve as a mass-composition-sensitive observable at HAWC over a multi-TeV energy range, and HAWC's standard gamma-ray LDF would be validated for cosmic-ray showers. The work has clear strengths: it compares four LDFs on the same data set, presents the energy dependence of fit quality, evaluates composition sensitivity with a quantitative FOM, and checks the data against Monte Carlo predictions. However, the main conclusions rest on chi-square comparisons that use an unvalidated empirical charge-error model, and the text contains an internal tension between the statements that no LDF describes the data satisfactorily and that HAWC's LDF gives a good description. These issues are directly load-bearing for the paper's central claims and need to be resolved before the results can be accepted.","major_comments":[{"comment":"The empirical effective-charge error function log10(Qerror) = 0.3 - 0.06667*log10(Qeff) for log10(Qeff) <= 3 and 0.1 otherwise is introduced without any calibration or validation. This function provides the uncertainties that enter every chi-square fit, so it directly controls the LDF ranking in Fig. 2 (left) and the conclusion that HAWC's LDF is the best description. The authors need to demonstrate that this error model is a faithful representation of the PMT charge uncertainties, for example by comparing it with single-PMT charge resolution measurements, repeated-event studies, or the bin-by-bin scatter of the radial distributions, and to show that the ranking is robust to plausible variations of the error model. As it stands, the central comparison is not yet established.","section":"Section 3, Eq. (3.1)"},{"comment":"There is an internal inconsistency in the goodness-of-fit claim. Section 6 concludes that for E < 10^4.5 GeV Eqs. (3.3) and (3.2) give a better description and for E > 10^4.5 GeV Eqs. (3.5) and (3.2) fit better, while Section 7 states that none of the selected LDFs describes the measured data satisfactorily for all radial ranges and energies but that HAWC's LDF gives a good description in 10^3.5-10^5.5 GeV in comparison with the other LDFs. Since all example chi2/NDOF values in Fig. 1 exceed 1, 'best among the tested functions' and 'good description in an absolute sense' are different statements. The authors need to specify a pre-defined acceptance criterion for chi2/NDOF (or equivalent), discuss the statistical and systematic errors affecting the chi-square values, and reconcile the two statements with the same numbers.","section":"Sections 6 and 7"},{"comment":"The comparison of the average lateral distributions with Monte Carlo predictions and the subsequent composition-sensitivity analysis are performed using only one high-energy hadronic interaction model, QGSJET-II-03, and only proton and iron primaries. The reconstructed energies also rely on the HAWC MC-based method described in [10], which introduces an additional model dependence into the comparison. Since the modified-scaling parameters beta, phi, and delta in Eq. (3.7) are stated to be derived from MC simulations, the comparison involving this LDF is partly circular. At minimum, the authors should state these caveats and quantify the effect of the interaction-model choice on the FOM values and on the energy dependence of the lateral age parameter.","section":"Section 5 and Fig. 2 (right)"}],"minor_comments":[{"comment":"The zenith-angle cuts are given as '< 16.7 degrees' in the abstract, 'θ < 16.71°' in Section 2, and 'θ < 16°' in the introduction; these should be made consistent.","section":"Abstract and Section 2"},{"comment":"The caption states that the error bars represent both 'the error on the mean and the sigma error' without explaining which quantity is shown; please define both and indicate which is plotted.","section":"Figure 1 caption"},{"comment":"The fitting procedure does not specify the number of radial bins, the number of degrees of freedom, the treatment of bins with low occupancy, or how the mean effective charge and its variance are computed in each bin; these details are needed for the chi-square comparison to be reproducible.","section":"Section 4"},{"comment":"The name of the hadronic interaction model appears both as 'QGSJet-II-03' and 'QGSJET-II-03'; please use a single spelling consistently.","section":"Throughout"},{"comment":"The text contains a typo ('V olcano') and would benefit from a brief explanation of why the NKG lateral age parameter, originally derived for electromagnetic showers, is applicable to hadron-induced air showers in the TeV-PeV range.","section":"Introduction"},{"comment":"The FOM definition does not include systematic uncertainties from energy reconstruction, shower-core location, or PMT saturation; given that these potentially dominate at large radii and high energies, a sentence on the impact of these systematics on the FOM would be valuable.","section":"Section 4, Eq. (4.1)"}],"recommendation":"major_revision","confidential_remarks":"This is a conference proceedings paper presenting preliminary results, and the bar for evidence may be slightly different from a full journal article. However, the central claims are not yet supported as written: the unvalidated error model in Eq. (3.1) affects all chi-square rankings, and the text explicitly undercuts its own goodness-of-fit conclusion. Both issues are fixable with additional analysis, so a major revision rather than rejection seems appropriate if the authors can provide the required calibration and reconcile the statements."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"What you should know: this is a short ICRC proceedings paper that does the first systematic comparison of four lateral distribution functions on HAWC cosmic-ray data, and then uses the lateral age parameter from the HAWC LDF to estimate proton-iron separation with a figure of merit. That is a useful and honest incremental step, not a breakthrough.\n\nThe paper earns credit for using real HAWC data with explicit quality cuts, fitting standard LDFs from KASCADE, ARGO-YBJ, and AGASA with a clear chi-square procedure, and for reporting the FOM rather than overclaiming composition resolution. It also openly states in the Discussion that none of the LDFs describes the data satisfactorily across all radial ranges and energies, which is a sign of care. Citations look appropriate, mostly to the original LDF papers and HAWC references.\n\nThe soft spots are real but addressable. The biggest one is that every chi2/NDOF value in Fig. 2 (left) is computed with the empirical PMT charge error model of Eq. (3.1), whose coefficients (0.3, -0.06667, 0.1) are simply asserted. No calibration against single-PMT charge resolution, photo-electron statistics, or bin-to-bin scatter is shown. Since the weights enter every fit, an unvalidated error model can change the ranking of the four LDFs, and therefore the claim that the HAWC LDF is the best description. That claim is further weakened by an internal tension: Section 6 says different LDFs do better at different energies, while the Conclusion first says none describes the data satisfactorily and then says the HAWC LDF gives a good description. Without a pre-specified acceptance threshold, 'good' is doing too much work. The modified scaling formalism parameters (beta, phi, delta) are taken from MC without showing the derivation, and only one hadronic interaction model (QGSJET-II-03) is tested. Also, the energy reconstruction and the MC predictions share the same HAWC simulation framework, so the agreement in Fig. 2 (right) is not an independent check of the interaction model.\n\nThat said, the central idea is plausible and the flaws are not fatal. This is a conference-length paper that should not carry the full weight of a composition analysis without further validation, but it is a legitimate first look. Who it is for: cosmic-ray and air-shower colleagues, especially those working on lateral distributions and composition with water-Cherenkov arrays. It deserves a serious referee: the data are real, the question is relevant, and the requested revisions (calibrate the error model, define a fit-acceptance criterion, show the MC parameter derivation) are reasonable. I would not desk-reject it; I would send it to review with moderate confidence that it can be strengthened.","headline":"Useful first HAWC cosmic-ray LDF comparison, but the ranking and composition claim rest on an uncalibrated error model and an internal consistency issue.","tokens_in":6787,"tokens_out":4447,"would_cite":false,"duration_ms":35925,"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":"HAWC's modified NKG lateral distribution function best describes cosmic-ray air showers from $10^{3.5}$ to $10^{5.5}$ GeV, and its fitted lateral age separates proton from iron primaries above $10^{3.8}$ GeV.","keywords":["lateral distribution function","extensive air showers","HAWC observatory","lateral age parameter","cosmic-ray mass composition","NKG function","mass composition sensitivity"],"falsifier":"Recompute the four LDF fits using charge uncertainties taken from PMT charge-resolution calibrations or from shower-to-shower fluctuations in place of Eq. (3.1), and check whether HAWC's LDF still has the lowest $\\chi^2/\\mathrm{NDOF}$; as a second check, evaluate the lateral-age proton-iron separation on independent simulated showers generated with a different hadronic interaction model.","tokens_in":5666,"feed_emoji":"🔭","tokens_out":10401,"duration_ms":88624,"temperature":0.7,"pith_summary":"This paper sets out to find which parameterization best describes the lateral distribution of particles in cosmic-ray air showers recorded by HAWC at 4100 m altitude. Using 1.3 days of 2016 data with energies between $10^{3.5}$ and $10^{5.5}$ GeV, it compares four lateral distribution functions by fitting the average deposited charge in 2 m radial bins and ranking the fits by $\\chi^2$ per degree of freedom. It concludes that HAWC's own modified NKG function, with a fixed Moliere radius $r_M = 124.21$ m and exponents $(s-3, s-4.5)$, gives the best overall description of cosmic-ray showers in that energy range. It then shows that the lateral age parameter $s$ recovered from that fit separates simulated proton and iron showers with a figure of merit above 1 for energies from $10^{3.8}$ to $10^{5.5}$ GeV. If correct, this makes a single LDF fit a usable mass-composition observable for HAWC's data.","feed_headline":"HAWC's own LDF best fits cosmic-ray showers from 3 TeV to 300 TeV","feed_subtitle":"The fitted lateral age separates proton and iron showers above 6 TeV.","key_machinery":"The load-bearing object is the modified NKG lateral distribution function $f(r)=A(r/r_M)^{s-3}(1+r/r_M)^{s-4.5}$, with $r_M=124.21$ m fixed at HAWC's altitude, where $A$ is a normalization and $s$ is the lateral age. The argument runs through three pieces: (i) a charge-error model $\\log_{10} Q_{\\mathrm{error}} = 0.3 - 0.06667\\log_{10}Q_{\\mathrm{eff}}$ for $\\log_{10}Q_{\\mathrm{eff}} \\le 3$ and $0.1$ above, which sets the weights in every fit; (ii) a $\\chi^2$ per degree of freedom comparison of four LDFs on binned mean effective charge; and (iii) the figure of merit $\\mathrm{FOM}=|s_{\\mathrm{Fe}}-s_p|/\\sqrt{\\sigma_p^2+\\sigma_{\\mathrm{Fe}}^2}$ from simulated proton and iron showers, which converts the fitted age into a composition-separation metric.","core_discovery":"The central claim is that the NKG-style function previously used by HAWC for gamma-ray showers, $f(r) = A (r/r_M)^{s-3}(1+r/r_M)^{s-4.5}$ with $r_M = 124.21$ m, also gives a good description of cosmic-ray-induced showers, and in a $\\chi^2/\\mathrm{NDOF}$ comparison over $10^{3.5}$ to $10^{5.5}$ GeV it performs at least as well as three literature LDFs across the full energy range. The fitted lateral age parameter $s$ varies with primary mass: iron showers are older (larger $s$) and protons are younger, and the separation, quantified by the figure of merit $|s_{\\mathrm{Fe}}-s_p|/\\sqrt{\\sigma_p^2+\\sigma_{\\mathrm{Fe}}^2}$, exceeds 1 for $E > 10^{3.8}$ GeV, peaks near $E = 10^{4.8}$ GeV, and remains about 1.75 at $E = 10^{5.5}$ GeV. The paper also reports that the average measured lateral distributions lie within the band predicted by the high-energy hadronic interaction model used in the simulations.","pith_inferences":["Beyond the paper: replacing the empirical charge-error function of Eq. (3.1) with PMT-level calibrated uncertainties could change the $\\chi^2$ ranking of the four LDFs, although the age-based composition result depends mainly on the HAWC LDF itself.","Beyond the paper: combining lateral age with other composition-sensitive observables, such as muon content or shower-front curvature, may extend the useful range beyond the energy where the FOM begins to decrease.","Beyond the paper: repeating the four-function comparison on showers simulated with alternative hadronic interaction models would test whether the LDF ranking is model-dependent."],"forward_implications":["The lateral age parameter from a single HAWC LDF fit can serve as a mass-composition observable between $10^{3.8}$ and $10^{5.5}$ GeV, with proton-iron separation above one standard deviation.","HAWC's existing gamma-ray LDF can be reused for cosmic-ray analyses without introducing a new parameterization.","Average lateral distributions measured by HAWC are consistent with the predictions of the QGSJET-II-03 hadronic interaction model across $10^{3.5}$ to $10^{5.5}$ GeV.","The proton-iron figure of merit peaks near $E = 10^{4.8}$ GeV and is still about 1.75 at $E = 10^{5.5}$ GeV."],"supporting_citations":[{"why":"This reference supplies the modified NKG form and the Moliere radius $r_M = 124.21$ m used for HAWC data.","marker":"[5]"},{"why":"This reference provides the alternative NKG LDF with a Gamma-function normalization factor used in the comparison.","marker":"[4]"},{"why":"This reference provides the alternative NKG-like LDF with a fixed scale radius of 30 m used in the comparison.","marker":"[13]"},{"why":"This reference provides the scaling-formalism LDF that the paper modifies for the comparison.","marker":"[14]"},{"why":"This reference is the air-shower simulation code used to generate the proton and iron showers.","marker":"[6]"},{"why":"This reference is the high-energy hadronic interaction model used for the simulated shower predictions.","marker":"[8]"},{"why":"This reference provides the energy reconstruction method and the spectrum used in the simulations.","marker":"[10]"},{"why":"This reference defines the figure of merit used to quantify proton-iron separation.","marker":"[15]"}],"fun_headline_variants":["HAWC lateral age separates proton and iron showers","Cosmic-ray mass read from HAWC's lateral age parameter","HAWC LDF fit: lateral age tracks cosmic-ray composition","One LDF fits HAWC showers, age reveals composition","HAWC finds lateral age sensitive to cosmic-ray mass"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The comparison that picks HAWC's LDF rests on the empirical charge-error function of Eq. (3.1), whose coefficients are asserted without calibration; if those errors do not represent the real PMT charge uncertainties, the $\\chi^2$ values and the LDF ranking could change.","fun_headline_variants_meta":{"raw":{"variants":["HAWC lateral age separates proton and iron showers","Cosmic-ray mass read from HAWC's lateral age parameter","HAWC LDF fit: lateral age tracks cosmic-ray composition","One LDF fits HAWC showers, age reveals composition","HAWC finds lateral age sensitive to cosmic-ray mass"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000271,"raw_usage":{"total_tokens":1665,"prompt_tokens":1015,"completion_tokens":650,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":631,"completion_tokens_details":{"reasoning_tokens":570}},"tokens_in":631,"tokens_out":650,"duration_ms":7087,"temperature":1.0,"reasoning_tokens":570,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:53:16.406977+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the four LDF fits using charge uncertainties taken from PMT charge-resolution calibrations or from shower-to-shower fluctuations in place of Eq. (3.1), and check whether HAWC's LDF still has the lowest $\\chi^2/\\mathrm{NDOF}$; as a second check, evaluate the lateral-age proton-iron separation on independent simulated showers generated with a different hadronic interaction model.","supporting_citations":[{"cited_title":"Electron, Muon, and Hadron Lateral Distributions Measured in Air-Showers by the KASCADE Experiment","cited_arxiv_id":"astro-ph/0004233","evidence_quote":"This reference provides the alternative NKG LDF with a Gamma-function normalization factor used in the comparison."},{"cited_title":"EAS age determination from the study of the lateral distribution of charged particles near the shower axis with the ARGO-YBJ experiment","cited_arxiv_id":"1707.01412","evidence_quote":"This reference provides the alternative NKG-like LDF with a fixed scale radius of 30 m used in the comparison."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"This reference provides the scaling-formalism LDF that the paper modifies for the comparison."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"This reference defines the figure of merit used to quantify proton-iron separation."}],"review_version":1}