{"id":"00cdc207-aab3-45ed-9ae7-9ef998bb2e89","arxiv_id":"1908.01622","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"HiSPARC's school-based scintillator stations detect cosmic-ray air showers, reconstruct directions with about 6 to 8 degree uncertainty, and produce a rough energy spectrum whose slope is close to the known cosmic-ray spectrum.","lead":"HiSPARC is a network of about 140 low-cost cosmic-ray detector stations, most of them hosted by high schools, that has been taking data since 2003. This paper documents that the detectors match simulations, reconstruct shower directions to about 6 degrees, and can roughly estimate cosmic-ray energies.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Energy-scale claim rests on a proton-only simulated LDF applied to mixed-composition showers; the KASCADE cross-check was too limited to independently validate it.","rationale":"The strongest claim has two parts: direction reconstruction (~6.1 degrees) and energy reconstruction (spectrum slopes 2.85/2.86). The direction part is independently supported: KASCADE (0.3 degree truth) and the close-station comparison in Table 1 both give ~6 degree uncertainties. I therefore see no reason to question that component. The energy part is not independently validated in the same way. The paper's own text flags the necessary caveats: the LDF is fitted to proton-only CORSIKA showers (Section 7), the age parameter is constant, the core-position fit is bimodal at higher energies (Figure 27), the flux normalization is not corrected for efficiency and cuts, and the KASCADE energy comparison was too limited. These are not pipeline artifacts; they are explicit limitation statements. The agreement between stations 501 and 510 is a consistency check, not a calibration check, because both use the same simulation-derived LDF and effective area. The central claim that 'the resulting energy spectrum ... does not deviate much from the known value of 2.7' is therefore only as strong as the untested assumption that a proton-only, fixed-age LDF describes mixed-composition showers at four sampling points. A Monte Carlo closure test with mixed composition through the same pipeline would settle whether this assumption biases the slope and energy scale. If it passes, the conditional concerns are resolved; if it fails, the energy-spectrum claim should be downgraded. This supports the reader's conditional verdict rather than changing it.","tokens_in":24566,"tokens_out":7435,"duration_ms":76126,"concrete_test":"Run a CORSIKA mixed-composition (p, He, CNO, Fe) sample in the 10^14.5-10^16.5 eV range through the full HiSPARC detector simulation and the exact Section 7 reconstruction pipeline (same LDF, same cuts, same effective-area parametrisation). Compare the reconstructed log(E) and the derived spectrum slope to the input values, binning by primary mass and core distance. If the mean log(E) residual exceeds ~0.1 (i.e., ~26% in energy) or the fitted slope differs by more than ~0.1, the proton-only LDF assumption is falsified and the claimed slopes 2.85/2.86 are not secure.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central scientific-viability claim includes the energy reconstruction of Section 7, which uses the modified NKG lateral distribution (Eqs. 7-9) with constants r0=29.6, a=-0.566, b=-2.57, p=6.937, c=0.797, d=17.62. These constants are obtained only from proton-initiated, perpendicular-incidence CORSIKA showers (Fig. 26), and the age parameter is fixed. The same LDF is then applied to real, mixed-composition, inclined showers sampled at just four detector positions per station. The paper itself flags two consequences: a bimodal core-position ambiguity at higher energies (Fig. 27) and an uncorrected flux offset due to efficiency and analysis cuts. It also states that the KASCADE-HiSPARC energy comparison was too limited for a decisive analysis. Because stations 501 and 510 share the same reconstruction and the same simulation-based effective area, their mutual agreement (Fig. 28) does not validate the absolute energy scale or the fitted slopes 2.85 and 2.86. If the LDF shape or fixed age parameter is composition- or energy-dependent, the reconstructed core position, A, and hence log(E) (Eq. 9) are biased; the slope comparison to 2.7 is then not evidence of correctness.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents the HiSPARC extensive air shower experiment, a distributed network of low-cost two- and four-scintillator stations hosted largely at high schools. It describes the detector hardware, GEANT4-based single-muon and EAS response simulations, DAQ and GPS timing, trigger efficiency parametrization, direction reconstruction, and a single-station energy reconstruction based on a modified NKG lateral distribution. The main quantitative claims are that a four-detector HiSPARC station reconstructs EAS directions with an average 1-sigma uncertainty of about 6.1 degrees (validated against KASCADE), that the same station can estimate primary energies using the modified NKG function, and that the resulting energy spectrum has fitted slopes of 2.85 and 2.86 between 10^14.8 and 10^15.5 eV, close to the expected value of 2.7.","tokens_in":24897,"tokens_out":4284,"duration_ms":48705,"significance":"The paper provides a detailed, honest characterization of a low-cost air-shower detector and demonstrates meaningful scientific validation: the GEANT4 detector model is checked against table-top and single-muon measurements; the direction reconstruction is compared with KASCADE (0.3-degree accuracy) and with independent station pairs; and the trigger efficiency is studied with CORSIKA. The open software and data access are genuine strengths for an education-linked project. If the energy-reconstruction claims can be supported against composition and systematic uncertainties, the paper would show that even a single four-detector HiSPARC station yields scientifically useful energy information. The authors also clearly state several limitations, which is commendable but does not by itself resolve the load-bearing issues described below.","major_comments":[{"comment":"The authors should either add a systematic study of LDF variation (e.g., mixed-composition CORSIKA showers, several hadronic interaction models, varied age parameters), or explicitly reframe the energy spectrum as an illustrative demonstration and remove the implication that it validates the absolute energy scale.","section":"Section 7, Eqs. (7)-(9), Fig. 26"},{"comment":"As written, Fig. 28 gives the impression of a spectral-shape measurement, but the analysis is a relative demonstration whose selection efficiency has not been folded in or shown to be energy-flat.","section":"Section 7, Fig. 28"},{"comment":"This is closely related to the previous comment, but it identifies a specific missing ingredient in the chain from simulated efficiency to observed flux.","section":"Section 8.2, Eqs. (10)-(13), Figs. 32-33"}],"minor_comments":[{"comment":"The PMT response function in Eq. (1) is presented without units for the parameters a, b, c, d and for x; please specify the domain and units (volts, inferred pulse heights, etc.) so the parametrization is unambiguous.","section":"Section 2.2, Eq. (1)"},{"comment":"The y-axis label 'Measured MIP-peak value [mVns]' and the x-axis label 'Simulated MIP-peak value [# photons]' mix units; please clarify how the simulated photon count is converted to the measured mVns scale.","section":"Fig. 10"},{"comment":"The sentence 'the uncertainty obtained from the simulations slightly underestimates the real direction reconstruction performance' can be misread; since 7.7 degrees is larger than 6.1 degrees, the simulation gives a larger (more conservative) uncertainty. Please rephrase to make the direction of the comparison explicit.","section":"Section 6, Figs. 24-25"},{"comment":"The criterion 'If the best chi2 value of one of the two stations is below 5, the event is discarded' should be accompanied by the number of degrees of freedom in the fit, otherwise the cut value is not interpretable.","section":"Section 7, analysis cuts"},{"comment":"Several quantitative statements (energy-dependent detection efficiencies, mini-shower contributions, and the four-station pair analysis) rely on unpublished or internal notes [27], [53]. Please either summarize the relevant numbers in the text or add a note that these results are preliminary.","section":"References [27] and [53]"},{"comment":"For the trigger-efficiency parametrization in Eqs. (10)-(13), the fitted values of alpha, sigma, lambda, chi, and rho are listed but without uncertainties or the number of simulated showers; please provide these details so readers can assess the fit quality.","section":"Section 8.2"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid instrument-description paper and the direction-reconstruction and detector-response parts are well supported. The main risk to the central 'scientific data can be obtained' claim is Section 7: the energy spectrum is built on a proton-only LDF, uncorrected normalization, and selection cuts, and yet it is presented as being close to the known slope of 2.7. I would recommend that the authors either provide a real systematic study (mixed composition, multiple hadronic models, energy-dependent efficiency) or explicitly downgrade the energy reconstruction to a preliminary demonstration. In the latter case, the paper would still be publishable as an instrument paper, but the abstract and conclusions should not imply that the energy scale is validated. The authors might also consider citing the KASCADE-Grande energy comparison more completely, since the current sentence ('too limited for a decisive analysis') weakens the claim while the figure still shows an apparent consistency."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a genuine instrument paper, not a stunt. The genuinely new content is modest — much of the hardware, DAQ, and analysis is compiled from theses and internal notes — but the KASCADE cross-calibration is real external evidence: a four-detector HiSPARC station reconstructs shower directions with about 6.1 degrees average 1-sigma uncertainty against KASCADE, and the pair-wise station comparisons agree. That part holds up. The detector simulation also matches single-muon and EAS pulse-height data, and the public data pipeline plus the SAPPHiRE and jSparc tools are practical contributions for education and for anyone who wants to reuse the array.\n\nThe soft spot is Section 7. The energy reconstruction uses a modified NKG lateral distribution whose parameters come only from proton-induced, perpendicular CORSIKA showers, with a fixed age parameter, applied to mixed-composition, inclined showers sampled at only four detector positions. The paper itself flags the consequences: a bimodal core-position ambiguity at higher energies, an uncorrected flux offset, and analysis cuts that discard low-spread events. It also says the KASCADE energy comparison was too limited for a decisive check. Given that, the station-to-station agreement in Figure 28 validates reproducibility, not the absolute energy scale or the fitted slopes. The slopes 2.85 and 2.86 being close to 2.7 is suggestive, but it is not evidence that the energy scale is correct. I would want a systematic study of composition, LDF age, and trigger-efficiency uncertainties before using the absolute energies.\n\nMinor: the aluminum-reflectivity tuning in GEANT4 (0.88 to 0.93, to match single-muon data) is a disclosed free parameter that only scales photon counts, but it feeds the same simulation chain used for the effective area. The citation pattern leans heavily on theses and internal notes — refs 13, 42, 45, 52, 53 — which makes independent verification harder, but the paper is transparent about that rather than hiding it.\n\nWho is this for: the cosmic-ray education community, builders of low-cost arrays, and anyone using HiSPARC public data. It deserves a serious referee, not a desk rejection. If I were handling it, I would accept after revision, with the energy-spectrum claims softened to preliminary and a proper systematic-uncertainty section added or explicitly deferred to a separate paper.","headline":"A solid instrument report for a long-running school-based array: the KASCADE cross-calibration makes the direction reconstruction credible, but the energy spectrum rests on a proton-only simulated LDF and should be treated as preliminary.","tokens_in":25441,"tokens_out":1934,"would_cite":true,"duration_ms":21337,"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":"A school-based network of cheap scintillator stations can reconstruct cosmic-ray showers to about 6 degrees.","keywords":["cosmic rays","extensive air showers","scintillation detector","HiSPARC","cosmic ray energy spectrum","shower direction reconstruction","NKG lateral distribution","outreach"],"falsifier":"Take a four-detector station whose showers are also seen by a high-precision array with accurate core positions, and compare energy estimates for events with cores inside versus outside the station; if the fixed LDF is correct, the two classes agree after efficiency corrections, while a composition or shape bias would make them disagree systematically.","tokens_in":24403,"feed_emoji":"🌌","tokens_out":4186,"duration_ms":38849,"temperature":0.7,"pith_summary":"This paper makes the case that a network of small, inexpensive scintillator stations run largely by high schools can do real cosmic-ray physics, not just outreach. It argues that a single four-detector HiSPARC station reconstructs the arrival direction of an extensive air shower with an average 1-$\\sigma$ uncertainty of about 6.1 degrees, as checked against the much more precise KASCADE experiment, and can estimate the primary energy by fitting a simplified Nishimura-Kamata-Greisen lateral distribution to the particle densities at four detector positions. Using two nearby stations, the reconstructed energy spectrum between $10^{14.8}$ and $10^{15.5}$ eV has fitted slopes of 2.85 and 2.86, close to the accepted value of 2.7. The paper's broader claim is that such stations are scientifically viable detectors, not just teaching tools.","feed_headline":"School-based cosmic-ray network resolves showers to 6 degrees","feed_subtitle":"Cheap four-scintillator stations also estimate primary energies, matching the known spectrum slope near 2.85.","key_machinery":"The load-bearing object is the modified NKG lateral distribution function with all shape parameters fixed, so that fitting the four measured particle densities determines only the core position and one energy scale. Around it, the paper builds a CORSIKA-based simulation chain for shower generation, a GEANT4 detector response simulation validated against single-muon measurements, and a flat-front triangulation algorithm for directions. The fixed LDF turns a four-point measurement into a one-parameter energy fit, which is what makes single-station energy reconstruction possible.","core_discovery":"The central discovery is that a four-scintillator station, with detectors a few metres apart, samples enough of an air-shower footprint to reconstruct both direction and energy with useful accuracy. Direction reconstruction by flat-front triangulation of arrival times yields an average 1-$\\sigma$ uncertainty of $6.1^\\circ$ when validated against KASCADE, and the same algorithm applied to four closely spaced stations gives pairwise differences between $5.93^\\circ$ and $6.37^\\circ$. Energy reconstruction uses a modified NKG formula $N(r) = A (r/r_0)^a (1 + r/r_0)^b$ with fixed parameters $r_0=29.6$, $a=-0.566$, $b=-2.57$, an obliquity correction $A_\\perp = A \\exp[p(1/\\cos\\theta - 1)]$ with $p=6.937$, and $\\log E = c(\\log A_\\perp + d)$ with $c=0.797$, $d=17.62$; the resulting spectrum has slopes 2.85 and 2.86 between $10^{14.8}$ and $10^{15.5}$ eV, close to the known 2.7, with a flux offset attributed to detection efficiency and analysis cuts.","pith_inferences":["If the fixed proton-only LDF is applied to mixed-composition real showers, the energy scale will carry a composition-dependent bias; a natural extension is to include separate proton and iron templates and treat composition as a nuisance parameter.","The same four-detector geometry could be used to test shower-front curvature: the 6-degree resolution is dominated by timing jitter and the flat-front assumption, so adding curvature parameters would be a direct, testable extension.","The paper's validation method generalises: any small array can certify itself by embedding one station in a high-precision array or by cross-comparing overlapping stations, as done here with stations 501 and 510.","A classroom network with this calibration path could in principle monitor transient phenomena, such as solar-particle events or lightning-related modulation, using the existing trigger and GPS infrastructure."],"forward_implications":["A four-detector station can serve as a standalone cosmic-ray observatory, producing direction and energy information for showers above roughly $10^{14.5}$ eV.","The energy spectrum slope measured by such stations, 2.85-2.86, is consistent with the canonical value of 2.7 within the stated systematic simplifications, suggesting the method captures real spectral information.","The trigger efficiency parametrisation in eq. (10) allows a station's effective exposure to be computed, so rates can be converted into fluxes.","Clusters of stations, such as the Science Park cluster, can resolve the core-position ambiguity by comparing energy estimates from multiple stations."],"supporting_citations":[{"why":"Supplies the high-precision KASCADE reference directions (0.3 degree accuracy) used to validate HiSPARC's direction reconstruction.","marker":"[7]"},{"why":"Documents the integration of a HiSPARC four-detector station into KASCADE and the 2008 dataset used for the comparison.","marker":"[13]"},{"why":"Generates the simulated air showers that define the LDF, trigger efficiency, and effective surface area.","marker":"[31]"},{"why":"Provides the GEANT4 detector response simulation validated against single-muon measurements.","marker":"[26]"},{"why":"Supplies the accepted cosmic-ray spectrum slope of 2.7 and the NKG background used for comparison.","marker":"[1]"},{"why":"Provides the triangulation algorithm used for shower direction reconstruction.","marker":"[52]"},{"why":"Gives the pairwise station-comparison results in Table 1 that corroborate the KASCADE-based uncertainty.","marker":"[53]"},{"why":"Supplies the muon momentum and zenith-angle distributions used in the detector response simulation.","marker":"[29]"},{"why":"Sets the GPS timing accuracy of 15 ns that underpins the inter-station synchronization.","marker":"[34]"}],"fun_headline_variants":["High-school network resolves cosmic-ray directions to 6 degrees","Student-built array reconstructs shower directions within 6 degrees","Frugal scintillator stations deliver 6-degree cosmic-ray accuracy","School-based cosmic-ray array matches energy spectrum slope near 2.8","Global student network achieves 6-degree air-shower direction fix"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The energy reconstruction assumes that a lateral distribution function fitted to proton-only simulated showers, with its age parameter fixed, describes real showers of mixed composition when sampled at just four detector positions.","fun_headline_variants_meta":{"raw":{"variants":["High-school network resolves cosmic-ray directions to 6 degrees","Student-built array reconstructs shower directions within 6 degrees","Frugal scintillator stations deliver 6-degree cosmic-ray accuracy","School-based cosmic-ray array matches energy spectrum slope near 2.8","Global student network achieves 6-degree air-shower direction fix"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000288,"raw_usage":{"total_tokens":1719,"prompt_tokens":1004,"completion_tokens":715,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":620,"completion_tokens_details":{"reasoning_tokens":628}},"tokens_in":620,"tokens_out":715,"duration_ms":7739,"temperature":1.0,"reasoning_tokens":628,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:07:43.113568+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a four-detector station whose showers are also seen by a high-precision array with accurate core positions, and compare energy estimates for events with cores inside versus outside the station; if the fixed LDF is correct, the two classes agree after efficiency corrections, while a composition or shape bias would make them disagree systematically.","supporting_citations":[{"cited_title":"Antoni, W","cited_arxiv_id":null,"evidence_quote":"Supplies the high-precision KASCADE reference directions (0.3 degree accuracy) used to validate HiSPARC's direction reconstruction."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Documents the integration of a HiSPARC four-detector station into KASCADE and the 2008 dataset used for the comparison."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Generates the simulated air showers that define the LDF, trigger efficiency, and effective surface area."},{"cited_title":"Allison, K","cited_arxiv_id":null,"evidence_quote":"Provides the GEANT4 detector response simulation validated against single-muon measurements."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the accepted cosmic-ray spectrum slope of 2.7 and the NKG background used for comparison."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the triangulation algorithm used for shower direction reconstruction."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the pairwise station-comparison results in Table 1 that corroborate the KASCADE-based uncertainty."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Sets the GPS timing accuracy of 15 ns that underpins the inter-station synchronization."}],"review_version":1}