{"id":"c37c51c4-b1ad-4afa-8fb2-fc9a7cc08ba3","arxiv_id":"2507.17244","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A single hybrid station's multi-detector correlations and a transformer-based muon-direction reconstruction could give access to electromagnetic and muonic shower energy spectra, shown here only in simulation.","lead":"Two detector-analysis strategies are proposed for a single hybrid cosmic-ray station: correlating SSD, WCD, and RPC signals to sense changes in the electromagnetic shower spectrum, and using machine learning on WCD/RPC muon signals with muon production depth to probe the muon energy spectrum. The paper is a simulation-based feasibility study tied to a prototype being tested at the Pierre Auger Observatory, not yet a measured spectrum.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The muon-spectrum claim depends on a 5.8° muon-direction resolution that, at r=1000 m and typical production heights, maps to roughly 1 km of vertical uncertainty; the toy MC in Fig.","rationale":"The paper is a feasibility study, and the first method (R/θ estimators) is honestly framed as an intermediary between observables and spectral shapes rather than a direct measurement; its Tyvek-ageing robustness test is a useful, if partial, check. The second method, however, carries the stronger claim: that a single hybrid station can access the muon energy spectrum. The reader's weakest assumption correctly identifies the transfer of idealized simulated performance to real shower conditions. My stress test sharpens this into a concrete geometric problem: a 5.8° angular error, which the paper itself quotes as the current idealized resolution, corresponds to roughly a kilometer of production-height uncertainty at the stated geometry. The toy MC in Fig. 5 (right) does not explicitly show that this uncertainty is propagated into the comparison between z_time_rec and z_angle_rec, and the transformer was trained without EM contamination. Because the MPD inversion is purely geometric, this is not a matter of statistics or calibration but of whether the proposed observable can survive its own input resolution. The concern does not move the overall verdict: the paper remains a plausible feasibility proposal whose strongest claim is conditional on exactly this missing validation. I therefore keep the reader's conditional verdict unchanged.","tokens_in":5949,"tokens_out":6758,"duration_ms":80670,"concrete_test":"Re-run the toy Monte Carlo behind Fig. 5 (right) with a Gaussian smearing of σ=5.79° applied to z_angle_rec, keeping all other assumptions unchanged, and then compare the smeared z_angle vs z_time distributions for the 0.5, 1, and 3 GeV muon samples. If the energy-class separation or the ~500 MeV cutoff disappears under this smearing, the central muon-spectrum claim is not supported by the current evidence; if the separation survives, the claim gains concrete support.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim for the muonic component is that a ~5.8° muon zenith resolution from the WCD+RPC transformer, combined with the MPD, provides access to the muon energy spectrum through the kinematical delay. The load-bearing assumption is that this angular resolution can be propagated through the geometric inversion of Eq. (1) without erasing the signal. For a vertical shower and a station at r=1000 m, a muon produced at height z≈3 km arrives with tanθ≈r/z; a 5.8° error in θ shifts the inferred production height by δz≈(r/sin²θ)δθ≈1 km, and the shift grows rapidly for higher production altitudes. This is the same order as, or larger than, the separation between the muon energy classes shown in Fig. 5 (right). The paper does not state that this smearing was included in the toy MC, nor does it quantify how EM contamination and multi-muon occupancy degrade the transformer's 5° performance, which was obtained with isolated injected 1 GeV muons (as the text explicitly admits). If the smearing is not included, the apparent abrupt cutoff near 500 MeV is not a robust prediction. The abstract's claim of 'access to its energy spectrum' is therefore not yet established by the analysis shown.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes two simulation-based strategies for a single hybrid station (SSD+WCD+RPC) at a cosmic-ray array. The first strategy modifies the high-energy electromagnetic tail in CORSIKA simulations and defines two estimators, R and theta, from differences between SSD-WCD and WCD-RPC signal correlation plots; it shows that these estimators respond to artificial spectral modifications and are robust to a ~7% Tyvek reflectivity aging. The second strategy trains a Transformer network on WCD PMT time traces to reconstruct muon zenith angle, reports an angular resolution of about 5.8 degrees using injected 1 GeV muons under idealized no-EM-contamination conditions, and uses a toy Monte Carlo together with the Muon Production Depth (MPD) algorithm to suggest sensitivity to O(GeV) muon energies through the kinematical delay. The abstract claims that these approaches allow extraction of the high-energy tail of the electromagnetic spectrum, the low-energy tail of the muonic spectrum, and access to the muon energy spectrum.","tokens_in":6281,"tokens_out":5273,"duration_ms":57022,"significance":"If established, the claims would be valuable: a single hybrid station that can constrain both shower components' spectra would be a practical complement to AugerPrime and would offer a new handle on the muon puzzle. The paper's concrete strengths are the hybrid detector concept, the use of CORSIKA and Geant4 simulations, the Tyvek-aging robustness check for R and theta, and a clear machine-learning architecture for direction reconstruction. As presented, however, the evidence is at the level of a feasibility study: the first method demonstrates sensitivity to artificial spectral modifications, and the second method demonstrates a preliminary direction reconstruction in an idealized simulation. Real-data validation, propagation of the angular-resolution uncertainty, and a demonstration of the inversion from observables to a spectrum are not yet shown, so the spectral-extraction claims go beyond what the analysis supports.","major_comments":[{"comment":"The abstract and summary claim that the first approach allows 'extraction of the high-energy tail of the electromagnetic spectrum and the low-energy tail of the muonic spectrum,' but Sec. 2 only shows that R and theta respond to an artificial modification of the EM spectrum in CORSIKA simulations; no inversion, calibration, or fitting procedure that recovers a spectrum from R and theta is presented. The text itself states that R and theta 'do not directly measure the particle energy spectrum,' and the muonic low-energy-tail extraction is deferred entirely to ref. [12] with no supporting analysis in this manuscript. The claim should be either backed by an explicit extraction procedure or reformulated as a sensitivity study.","section":"Sec. 2 and Abstract"},{"comment":"The 5.8-degree resolution in Fig. 5 (left) is obtained, by the authors' own statement, with injected 1 GeV muons and 'without electromagnetic contamination.' The toy MC in Fig. 5 (right) computes z^angle_rec 'assuming the muon direction can be inferred using the neural network,' but the paper does not state whether the 5.8-degree residual distribution has been propagated into z^angle_rec. For a vertical shower with a station at r=1000 m and a typical production height of about 3 km, a 5.8-degree zenith-angle error shifts the inferred production height by roughly 1 km, which is comparable to or larger than the separation between the muon energy classes shown in Fig. 5 (right). Without this smearing and without a study of EM contamination and multi-muon occupancy, the abrupt cutoff near 500 MeV cannot be regarded as a robust prediction; the abstract's 'access to its energy spectrum' is therefore not yet established.","section":"Sec. 3, Fig. 5"},{"comment":"The inversion of Eq. (1) is the logical core of the second method, but it is not demonstrated. The kinematical delay <t_epsilon> is described as 'typically obtained from Monte Carlo shower simulations,' yet the paper does not show how z^time_rec and z^angle_rec are combined to extract <t_epsilon>, nor how uncertainties in the reconstructed direction, in the MPD input parameters, and in the Monte Carlo model of <t_epsilon> propagate into a muon energy spectrum. A closed-loop simulation that inputs a known muon spectrum and recovers it with quantified uncertainty is needed before the spectral-access claim can be supported.","section":"Sec. 3, Eq. (1)"},{"comment":"The demonstration of robustness against detector aging is limited to one mechanism (Tyvek reflectivity reduced by 7%) and one calibration procedure. The paper does not quantify how other optical-parameter changes, PMT gain drifts, or trigger-threshold variations affect R and theta, nor does it show the statistical precision with which R and theta could be measured in real events. Because the first method is explicitly framed as a way to 'identify the spectral modifications that best reproduce the observed detector responses,' a realistic uncertainty budget is required before the method can be compared with data.","section":"Sec. 2, Fig. 4"}],"minor_comments":[{"comment":"Equation (1) appears to contain a typographical or dimensional error: the expression inside the large parentheses mixes r^2/(c(t-<t_epsilon>)) with a bare c(t-<t_epsilon>) term; please correct the formula and verify all dimensions.","section":"Sec. 3, Eq. (1)"},{"comment":"The phrase 'somecm2' should read 'some cm^2'.","section":"Sec. 3, after Fig. 5"},{"comment":"The axis labels and units for the residual histogram are missing; please add them.","section":"Fig. 5 (left)"},{"comment":"The station trigger threshold '1.75 VEM' is introduced without a definition of how the VEM scale is calibrated in the simulation; please clarify.","section":"Sec. 2, trigger description"},{"comment":"The phrase 'for the first time' is used without a comparison to prior work; either substantiate it with references or remove it.","section":"Abstract"}],"recommendation":"major_revision","confidential_remarks":"This is a conference proceedings contribution that reads as a feasibility study. The main issue is the mismatch between the abstract's spectral-extraction claims and the demonstrated sensitivity-level analysis. A revision that either supplies the missing propagation/inversion studies or clearly scopes the paper as a proposal for future work would make the content appropriate for the venue."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper puts two new ideas on the table: R and theta from SSD-WCD-RPC signal correlations as handles on the electromagnetic spectral tail, and a transformer-based muon direction reconstruction feeding the MPD kinematical delay for the muon spectrum. Both are real combinations not in the cited literature, and the group is credible enough that these are worth a careful look.\n\nWhat the paper does well: the first method is honestly framed as sensitivity, not extraction. The R and theta response to artificial spectral modifications is shown, and the Tyvek aging check is a nice robustness test. The second method is also honest about its idealized training: the 5-degree resolution is quoted with no EM contamination, and the text says so. That is a good sign for scientific integrity.\n\nThe soft spot is the load-bearing step in the muon-spectrum argument. The stress test is correct: at r=1000 m and typical production heights, a ~5.8-degree zenith error translates to roughly 1 km of vertical uncertainty, which is the same order as the separation between the muon energy classes in Fig. 5. The paper does not state that this smearing is included in the toy MC, and if it isn't, the claimed cutoff near 500 MeV is not a robust prediction. EM contamination and multi-muon occupancy would only degrade the direction resolution further. So the abstract's 'access to its energy spectrum' is not yet supported by the analysis shown; 'sensitivity in a simplified simulation' is the accurate claim.\n\nThe first method has a smaller soft spot: no inversion from R and theta to the spectral shape, and no real-data validation. But it only claims sensitivity, and the cross-detector consistency argument is plausible. The companion paper [12] may carry more weight, but here it is a promise, not evidence.\n\nThis is a feasible-idea paper from an experienced group, suitable for the cosmic-ray community. It deserves a serious referee: the questions of whether the angular resolution smearing was propagated and whether the idealized direction reconstruction survives real shower conditions are exactly what peer review should probe. My recommendation: send it to review, but the authors should be pushed to either include the smearing or lower the claims to feasibility-level wording.","headline":"Two genuinely new detector-level estimators for shower spectra, but the muon-spectrum claim rides on an idealized 5-degree angular resolution that the toy MC likely does not smear.","tokens_in":6752,"tokens_out":1352,"would_cite":false,"duration_ms":16323,"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":"This paper claims that a single hybrid detector station can access the energy spectra of both the electromagnetic and muonic components of air showers, using signal-correlation estimators and muon-direction-based kinematical delay.","keywords":["extensive air showers","energy spectrum","electromagnetic component","muonic component","hybrid detector station","muon production depth","kinematical delay","transformer neural network"],"falsifier":"Take real WCD plus RPC data from the prototype station, select single-muon hits on an RPC pad, and compare the Transformer's reconstructed muon direction with the direction obtained from the arrival-time plane of the surrounding array; if the per-muon residual stays above about 5 degrees once electromagnetic contamination is included, or if the MPD-derived kinematical delay shows no correlation with muon energy in a beam-test or tagged calibration sample, the claimed access to the muon spectrum via kinematical delay does not survive.","tokens_in":5809,"feed_emoji":"💧","tokens_out":6176,"duration_ms":61836,"temperature":0.7,"pith_summary":"Ultra-high-energy cosmic rays are studied through the air showers they trigger, but ground detectors usually blur together the two shower components, the electromagnetic particles and the muons. This paper argues that a single station combining a scintillator, a water Cherenkov detector, and resistive plate chambers can separate and probe both energy spectra at once. The first method uses the different energy responses of the three detectors and two geometric estimators, R and theta, to become sensitive to modifications in the high-energy electromagnetic tail while staying robust to detector ageing. The second method reconstructs each muon's direction from PMT time traces and RPC pad hits with a Transformer neural network, then inverts the muon production depth relation to extract the kinematical delay and, with it, the muon energy distribution. If both hold, a compact station in an existing array could test hadronic-interaction models without a dense, large-area experiment.","feed_headline":"Two methods let one station read both shower energy spectra","feed_subtitle":"Scintillator, water, and plate signals separate the high-energy EM tail from low-energy muons at one site.","key_machinery":"The load-bearing object is the kinematical delay, the small extra arrival time a muon accumulates because of its finite energy and curved trajectory beyond the geometrical delay of a plane shower front. In the paper's height expression for the muon production depth, this term is the only one that depends on the muon energy distribution, so knowing the muon arrival direction lets the authors extract the delay from the measured arrival time and thereby constrain the muon spectrum. The other central objects are the estimators R and theta, defined as the distance and angle between the barycentres of the positive and negative regions in the relative-difference histograms of detector signal pairs; they translate spectral modifications into two numbers that can be compared directly between data and simulation.","core_discovery":"The paper claims that from one hybrid station one can access, for the first time, the high-energy tail of the electromagnetic spectrum and the low-energy tail of the muonic spectrum, and that combining muon direction with muon production depth gives access to the muon energy spectrum. The first claim is carried by R and theta, estimators built from the barycentres of positive and negative regions in the difference between modified and unmodified detector-signal correlation plots; these respond monotonically to changes in the high-energy electromagnetic tail and barely move under a 7% Tyvek reflectivity drop after standard vertical-equivalent-muon calibration. The second claim is carried by inverting the muon production depth height formula, whose kinematical-delay term is the only part that depends on muon energy; a Transformer trained on injected 1 GeV muons gives roughly 5-degree direction resolution with negligible bias, and a toy Monte Carlo for vertical showers at 1000 m from the core shows sensitivity to O(GeV) muons with an abrupt cutoff near 500 MeV from ionization losses. The paper presents these results as evidence of feasibility for a prototype station already under test, not as a fully validated measurement.","pith_inferences":["The authors do not test how electromagnetic contamination affects the 5-degree muon direction resolution; if contamination degrades it, the kinematical-delay inversion would still be recoverable by retraining the network on mixed showers and adding shower-distance or energy as auxiliary inputs, an extension the architecture already permits.","The R and theta method is demonstrated only at 320 m from the core for the electromagnetic tail, but the same multi-detector station could be applied to the muonic component at larger core distances, directly testing the claimed low-energy muon tail access.","One practical extension is to use the station-to-station scatter of R and theta across an array as a systematic check: genuine spectral modifications should correlate with reconstructed core distance, while detector miscalibration would not.","The MPD inversion could be cross-checked against the first method's muon spectral information whenever a station sees both components, providing an internal consistency test of the muon energy spectrum."],"forward_implications":["A single hybrid station can substitute for large-scale dedicated detectors when searching for spectral signatures of hadronic interaction models, because each station measures both the electromagnetic and muonic components locally.","The R and theta estimators can serve as simulation-to-data comparators for the high-energy electromagnetic tail, with the same station providing a built-in cross-check because genuine spectral changes must appear consistently in both detector-pair correlations.","The kinematical-delay route turns the muon production depth, normally a composition-sensitive observable, into an energy-resolving measurement for low-energy muons near the O(GeV) range.","The observed cutoff near 500 MeV, caused by ionization energy losses, can be used as an additional handle to test the low-energy muon spectrum.","The station-level trigger and analysis can be applied to existing array data without requiring full event reconstruction."],"supporting_citations":[{"why":"Supplies the air-shower Monte Carlo used to generate the reference showers and the modified spectra.","marker":"[5]"},{"why":"Provides the detector-response simulation for the three detector types in the station.","marker":"[6]"},{"why":"Documents the degradation of water-Cherenkov optical parameters over time, the ageing effect that the R and theta robustness test addresses.","marker":"[11]"},{"why":"Companion paper carrying further validation of the R and theta method for the muonic component and its dependence on shower reconstruction parameters.","marker":"[12]"},{"why":"Gives the time-to-depth reconstruction relation for muon production depth that the kinematical-delay inversion relies on.","marker":"[15]"},{"why":"Establishes that the kinematical delay materially influences the inferred muon production depth profile and depends on muon energy, motivating the inversion.","marker":"[16]"},{"why":"Demonstrates single water-Cherenkov direction reconstruction from PMT time traces with machine learning, the basis for the muon-direction network.","marker":"[17]"},{"why":"Extends the single-detector directional reconstruction setup with additional configurations, informing the WCD plus RPC approach.","marker":"[18]"}],"fun_headline_variants":["Single station splits EM and muon spectra","Two methods, one station, both spectra","Hybrid detector separates shower energy components","Machine learning aids muon energy spectrum from one station","EM and muon tails extracted at single hybrid station"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole argument assumes that the simulated detector traces, using clean injected muons with no other shower particles mixed in and reweighted spectra, represent real shower conditions closely enough that the reconstructed directions and the R and theta shifts track the true particle energy spectra rather than unrecognized contamination or calibration drift.","fun_headline_variants_meta":{"raw":{"variants":["Single station splits EM and muon spectra","Two methods, one station, both spectra","Hybrid detector separates shower energy components","Machine learning aids muon energy spectrum from one station","EM and muon tails extracted at single hybrid station"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000466,"raw_usage":{"total_tokens":2338,"prompt_tokens":972,"completion_tokens":1366,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":588,"completion_tokens_details":{"reasoning_tokens":1297}},"tokens_in":588,"tokens_out":1366,"duration_ms":10901,"temperature":1.0,"reasoning_tokens":1297,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T14:52:52.003414+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take real WCD plus RPC data from the prototype station, select single-muon hits on an RPC pad, and compare the Transformer's reconstructed muon direction with the direction obtained from the arrival-time plane of the surrounding array; if the per-muon residual stays above about 5 degrees once electromagnetic contamination is included, or if the MPD-derived kinematical delay shows no correlation with muon energy in a beam-test or tagged calibration sample, the claimed access to the muon spectrum via kinematical delay does not survive.","supporting_citations":[{"cited_title":"Heck, J.N","cited_arxiv_id":null,"evidence_quote":"Supplies the air-shower Monte Carlo used to generate the reference showers and the modified spectra."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the detector-response simulation for the three detector types in the station."},{"cited_title":"Studies on the response of a water-Cherenkov detector of the Pierre Auger Observatory to atmospheric muons using an RPC hodoscope","cited_arxiv_id":"2007.04139","evidence_quote":"Documents the degradation of water-Cherenkov optical parameters over time, the ageing effect that the R and theta robustness test addresses."},{"cited_title":"Assis, R","cited_arxiv_id":null,"evidence_quote":"Companion paper carrying further validation of the R and theta method for the muonic component and its dependence on shower reconstruction parameters."},{"cited_title":"Cazon, R","cited_arxiv_id":null,"evidence_quote":"Gives the time-to-depth reconstruction relation for muon production depth that the kinematical-delay inversion relies on."},{"cited_title":"Sensitivity of EAS measurements to the energy spectrum of muons","cited_arxiv_id":"1607.06760","evidence_quote":"Establishes that the kinematical delay materially influences the inferred muon production depth profile and depends on muon energy, motivating the inversion."},{"cited_title":"Alvarez-Muñiz, R","cited_arxiv_id":null,"evidence_quote":"Demonstrates single water-Cherenkov direction reconstruction from PMT time traces with machine learning, the basis for the muon-direction network."},{"cited_title":"Enhancing Neutrino Reconstruction in Water-Cherenkov Air Shower Arrays Using Multi-Photosensors","cited_arxiv_id":"2504.08652","evidence_quote":"Extends the single-detector directional reconstruction setup with additional configurations, informing the WCD plus RPC approach."}],"review_version":1}