{"id":"b78b2079-4983-4a17-831a-5b0cc4c9bac3","arxiv_id":"1908.07840","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"A symmetrization procedure for radio signals of inclined air showers reconstructs electromagnetic energy with a spread below 3 percent on simulated events.","lead":"Radio signals from cosmic-ray air showers are distorted by geometry and by a second emission component. This paper introduces a two-step correction that symmetrizes the signal so a simple function can estimate the shower's energy.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed <3% energy reconstruction is demonstrated only with the true shower-max depth; real applications require an Xmax assumption, and sensitivity to Xmax errors is not quantified.","rationale":"The reader's verdict is CONDITIONAL, and my stress-test identifies the same strongest concern: the quoted 'negligible bias and <3% spread' is obtained using the true shower maximum, while a real reconstruction must supply Xmax independently. This is not an internal contradiction, but it is a scope limitation on the central claim: the symmetrization procedure is defined relative to a point source at the shower maximum, and the charge-excess parameterization depends on dmax and rho_max, so the performance metric is conditional on a quantity that is not part of the reconstructed observables. The in-sample use of the same simulations for fitting Eq. (3.3) and for calibrating Eqs. (6.1)-(6.2) is a real but secondary concern; it would strengthen the paper to split the sample or use leave-one-out validation, but even a perfect out-of-sample charge-excess parameterization would not remove the Xmax-sensitivity problem. The concrete test proposed above directly isolates whether this concern lands by perturbing Xmax and repeating the reconstruction; if the energy estimator remains within 3% for Xmax shifts of realistic size, the concern is largely mitigated. I therefore keep the reader's CONDITIONAL verdict unchanged: the paper should be accepted only with an explicit discussion and quantification of the Xmax dependence and, ideally, a demonstration that Xmax can be reconstructed or that the method is robust to Xmax errors.","tokens_in":4834,"tokens_out":4374,"duration_ms":47683,"concrete_test":"Take a subset of the CoREAS simulations spanning the energy and zenith range, and for each event deliberately perturb the assumed Xmax by ±20 g/cm^2 and ±50 g/cm^2 (shifting dmax and rho_max accordingly). Re-run the early-late correction and charge-excess subtraction with the perturbed Xmax, fit Eq. (5.1), integrate to obtain E_rad^geo, apply Eq. (6.1), and compare the reconstructed electromagnetic energy with the true value. If the bias or spread exceeds the claimed 3% for realistic Xmax uncertainties, the Sec. 7 claim should be restated as conditional on an Xmax determination of stated accuracy, or an iterative Xmax-reconstruction step must be included.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim in Sec. 7 is that area integration of the symmetrized lateral distribution yields an electromagnetic-energy estimator with negligible bias and a spread below 3%. That claim is established only under the assumption that the shower maximum Xmax is known. The early-late correction in Sec. 2 projects antenna positions along the line of sight to a point source at the shower maximum and applies an inverse-square correction using R0 and R = R0 + x; the validation shown in Figs. 2 and 3 explicitly uses the true depth of the shower maximum. The charge-excess parameterization in Eq. (3.3) also depends on dmax and rho_max, both of which are determined by Xmax. In a real reconstruction Xmax is not known a priori, and the paper gives neither an iterative scheme nor a sensitivity analysis. A biased Xmax would systematically distort the early-late correction and the charge-excess subtraction, biasing the fitted fABCD and hence the integrated radiation energy. The in-sample fit of Eq. (3.3) and the joint energy calibration of Eqs. (6.1)-(6.2), both evaluated on the same 3111 CoREAS simulations, is a further limitation on the claimed 'universal' performance, but the Xmax dependence is the more load-bearing gap because it affects the applicability of the method itself in the intended experimental setting.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents a reconstruction model for radio emission from inclined air showers (zenith angles above 60°). The method consists of three steps: a geometrical early-late correction that projects ground-plane antenna positions into the shower plane along the line of sight to a point source at the shower maximum and applies an inverse-square fluence correction; a parameterization of the charge-excess fraction a(r, dmax, ρmax) (Eq. 3.3) used to isolate the pure geomagnetic fluence (Eq. 4.1); and a rotationally symmetric lateral distribution function (Eq. 5.1) whose area integral gives the geomagnetic radiation energy, corrected for geomagnetic angle and atmospheric density (Eqs. 6.1–6.2) and calibrated to the electromagnetic shower energy (Eq. 6.2). The model is evaluated on 3111 CoREAS simulations for proton and iron primaries at the Pierre Auger site, with energies from 10^18.4 to 10^20.2 eV, zenith angles from 65° to 80°, and eight azimuth angles. The central claim is that the resulting electromagnetic-energy estimator has negligible bias and a spread below 3%.","tokens_in":5239,"tokens_out":5337,"duration_ms":50770,"significance":"If the claimed performance could be established out of sample and without assuming knowledge of Xmax, the method would be a valuable, computationally light tool for inclined-shower radio reconstruction, complementary to muon-based measurements. The manuscript's strengths are its fully analytic and explicitly documented model, the transparent definition of all fit functions, and the use of a large CoREAS simulation set spanning wide ranges of energy and geometry. However, the current evidence is in-sample: the charge-excess parameterization and the energy-calibration relation are fitted to and evaluated on the same simulations, and the early-late correction is validated with the true shower-maximum depth. These issues do not invalidate the conceptual framework, but they mean the quoted <3% spread is not yet a demonstrated property of a complete reconstruction pipeline.","major_comments":[{"comment":"The early-late correction is validated using the true depth of the shower maximum, and both the projection geometry and the inverse-square correction assume a point source located at Xmax. In an actual reconstruction Xmax is not known a priori. The paper does not quantify the sensitivity of the corrected fluences, the fitted lateral-distribution parameters, or the integrated radiation energy to errors in Xmax, nor does it propose an iterative Xmax-estimation scheme. Because Eq. (3.3) also depends on dmax and ρmax, an incorrect Xmax would bias both the symmetrization and the final energy estimator. Please provide a sensitivity scan (e.g., ±50 g/cm²) or an explicit self-consistent reconstruction loop.","section":"Sec. 2, Eq. (2.1), Figs. 2–3"},{"comment":"The charge-excess parameterization in Eq. (3.3) is fitted to the same 3111 CoREAS simulations on which the symmetrization quality and the energy-reconstruction spread are then evaluated. Similarly, the joint fit of Eqs. (6.1)–(6.2) with parameters in Table 1 uses all of these simulations. The quoted '<3% spread' is therefore a measure of the fit quality on the training set, not an unbiased estimate of reconstruction performance. An out-of-sample evaluation is needed, for example via cross-validation per energy/zenith bin or comparison with an independent CoREAS set.","section":"Secs. 3, 4, 6, Eq. (3.3), Eq. (6.2), Table 1"},{"comment":"The claimed 'universal' parameterization is established for a single site: the Pierre Auger atmosphere, magnetic-field configuration, and observer altitude. The parameterization was intentionally reformulated using dmax and ρmax to be transferable, but no test with a different atmospheric model or magnetic-field strength is shown, and the analysis excludes geomagnetic angles below 20°. Please either restrict the universality claim to the tested configuration or add concrete transferability tests.","section":"Sec. 3, Eq. (3.3); Sec. 6, Table 1"},{"comment":"The paper acknowledges that the symmetrization is not fully successful in the inner region, where the parameterization tends to overestimate the charge-excess fraction, and that small asymmetries within concentric rings remain. The effect of these residual asymmetries on the fitted parameters A–D in Eq. (5.1) and on the integrated radiation energy is not quantified. Since the energy estimator is derived from that integral, the impact of these known residuals should be estimated or bounded.","section":"Sec. 4, Figs. 5–6"}],"minor_comments":[{"comment":"The comparisons shown in these figures have no statistical uncertainties; adding error bars or confidence bands would make the claimed 2–3% agreement easier to assess.","section":"Figs. 2, 3, 6, 8"},{"comment":"The reference value ⟨ρmax⟩ = 0.4 kg/m³ is introduced without stating whether it is the mean over the simulation set or a fixed atmospheric reference value; please clarify.","section":"Eq. (3.3)"},{"comment":"The variable x in R ≡ R0 + x is not defined precisely in the text; a sentence defining x and its sign relative to the shower axis would improve reproducibility.","section":"Eq. (2.1)"},{"comment":"The phrase 'we are still investigating alternative functions' indicates that the choice of fit function is not yet settled; this is acceptable for a proceedings paper but should be flagged as ongoing optimization rather than a final model recommendation.","section":"Sec. 5"}],"recommendation":"major_revision","confidential_remarks":"This is a well-written proceedings paper whose conceptual framework is clear and potentially useful. The two load-bearing gaps are the in-sample validation of the parameterizations and the assumption of known Xmax in the early-late correction. I would encourage the editor to request an out-of-sample evaluation and an Xmax-sensitivity study rather than reject the manuscript, since the method is documented in enough detail that these additions are feasible within the paper's scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a useful methods paper for inclined radio reconstruction, but the headline sub-3% energy spread is a training-set number and assumes the true Xmax. If you read it as a proof-of-principle with an explicit missing piece, it is solid; if you read it as a universal calibration, it is not there yet.\n\nWhat is actually new is equation 3.3, a compact parameterization of the charge-excess fraction as a function of lateral distance, distance to shower maximum, and air density at the maximum. The authors combine that with their earlier early-late correction into a full pipeline: symmetrize, fit a one-dimensional lateral distribution, integrate, and convert to electromagnetic energy. The pipeline is fully analytic, closed-form, and the steps are clearly documented. The paper also honestly notes small residual asymmetries and an imperfect inner-region correction. That is real work, and the direction is sensible.\n\nThe soft spot is the one the stress test identified: the central performance claim depends on knowing the true depth of the shower maximum. The early-late correction in Sec. 2 projects antenna positions to a point source at Xmax and is validated using the true Xmax. The charge-excess parameterization also depends on dmax and rho_max, both determined by Xmax. In a real event Xmax is unknown, and the paper gives neither an iterative scheme nor a sensitivity analysis. A biased Xmax would systematically distort both the early-late correction and the charge-excess subtraction, which would bias the fitted lateral distribution and hence the integrated energy. That is not a minor issue; it is the gap between a simulation-based method and an actual reconstruction.\n\nThe second issue is that the parameterization and the energy-relation fits are evaluated on the same 3111 CoREAS simulations used to fit them. The sub-3% spread is therefore a measure of in-sample fit quality, not a predictive performance claim. The authors do not give uncertainties on the fit parameters or on the reconstructed energies. The paper would be much stronger with a train/test split, a cross-validation, or even a small independent simulation set.\n\nCitation pattern looks clean, and the authors appropriately cite their own earlier early-late paper and the relevant literature. No red flags there.\n\nWho benefits: radio air-shower experimentalists, especially those working on inclined events at the Pierre Auger Observatory or similar sparse arrays. The paper deserves a serious referee, but the referee should push for a quantified Xmax dependence and an out-of-sample validation before the sub-3% claim is taken at face value.","headline":"A useful proof-of-principle for symmetrizing inclined air-shower radio signals, but the sub-3% energy claim is in-sample and assumes the true shower maximum, so it should be treated as a solid starting point rather than a finished reconstruction.","tokens_in":5679,"tokens_out":1161,"would_cite":true,"duration_ms":13124,"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":"Inclined air-shower radio signals can be symmetrized and integrated to give cosmic-ray energies with under 3 percent spread.","keywords":["radio emission","inclined air showers","air-shower reconstruction","cosmic-ray energy","geomagnetic radiation","charge-excess radiation","early-late effect","lateral distribution function"],"falsifier":"Split the 3111 CoREAS simulations in half: fit the charge-excess parameterization and energy calibration on one half, run the full symmetrization and integration on the other, and check whether the corrected geomagnetic radiation energy still has a spread below 3 percent and no bias. If the spread grows, the claimed universal accuracy comes from in-sample tuning rather than from the model itself.","tokens_in":4671,"feed_emoji":"📡","tokens_out":11613,"duration_ms":103365,"temperature":0.7,"pith_summary":"Radio emission from inclined air showers is promising for cosmic-ray energy measurement because it is a clean probe of the electromagnetic cascade, but the signal pattern on the ground is distorted by early-late geometry and by the mixture of geomagnetic and charge-excess radiation. This paper shows that both distortions can be removed analytically, leaving a rotationally symmetric lateral distribution that can be fit in one dimension and integrated to give the radiation energy. The authors argue that the resulting estimator for the electromagnetic cascade energy is essentially unbiased, with a spread under 3 percent across 3111 simulated showers. If this holds in practice, sparse radio arrays could reconstruct inclined showers without a full two-dimensional footprint model, complementing muon-based measurements.","feed_headline":"Cosmic-ray energies from inclined air showers pinned to 3 percent","feed_subtitle":"A two-step symmetrization flattens the radio footprint into one radial curve that yields the cascade energy.","key_machinery":"The carrying object is a three-stage analytic symmetrization. First, an early-late correction moves each antenna position to the shower plane along the line of sight to a point source at the shower maximum and rescales the energy fluence by the inverse-square of the distance, so that the shower-plane footprint becomes symmetric in the radial direction. Second, a universal parameterization of the charge-excess fraction $a=\\sin^2\\alpha\\, f_{\\mathrm{ce}}/f_{\\mathrm{geo}}$, depending only on axis distance, distance to shower maximum, and air density at the maximum, is combined with the known linear polarization of geomagnetic emission to compute the pure geomagnetic fluence at each antenna. Third, the resulting rotationally symmetric distribution is fit with $f_{ABCD}(r)=A\\exp(-Br-Cr^{2}-Dr^{3})$, a one-dimensional exponential of a cubic polynomial, whose area integral defines the radiation energy. The charge-excess parameterization is the piece that converts a two-dimensional asymmetric footprint into a single radial profile.","core_discovery":"On the paper's own terms, the discovery is a complete, closed-form recipe for turning the asymmetric radio signal distribution of an inclined air shower into a rotationally symmetric one. The early-late effect is corrected by projecting antenna positions onto the shower plane along the line of sight to a point source at the shower maximum and rescaling fluences by inverse-square distance; then a universal parameterization of the charge-excess fraction (as a function of axis distance, distance to shower maximum, and air density there) is used to isolate the pure geomagnetic fluence at each antenna. The symmetrized distribution is fit with a one-dimensional exponential-of-a-cubic lateral distribution function, and its integral over area yields a corrected geomagnetic radiation energy that follows the electromagnetic cascade energy through a quadratic power law. On 3111 CoREAS simulations of proton and iron showers with energies from $10^{18.4}$ to $10^{20.2}$ eV and zenith angles from $65^\\circ$ to $80^\\circ$, this estimator is reported to be unbiased with a spread below 3%.","pith_inferences":["We infer that a fieldable reconstruction chain must supply an estimate of the shower-maximum depth for the early-late correction, because the paper's validation uses the true value; testing with realistic shower-maximum uncertainties would show how much of the 3 percent spread is reserved for ideal conditions.","We infer that the universal charge-excess parameterization is in-sample by construction, since it is fitted to the same CoREAS library used for the performance evaluation; an independent simulation set with a different magnetic-field geometry or atmosphere would be needed to establish universality.","We infer that the same symmetrization idea might extend to less inclined showers, where early-late effects are weaker but geomagnetic and charge-excess asymmetries persist; the paper only demonstrates the method for zenith angles above 65 degrees.","We infer that if the density-at-maximum factor in the charge-excess function correlates with primary composition, the energy estimator could carry a weak composition dependence that would only show up when residuals are separated by primary species; the paper reports the joint spread but does not split it that way."],"forward_implications":["Radio reconstruction of inclined air showers can be done with a one-dimensional fit, so sparse antenna arrays with a small number of stations may be sufficient for energy measurement.","The energy estimator is unbiased for both proton and iron primaries over $10^{18.4}$ to $10^{20.2}$ eV, making it a possible basis for energy assignment at large-scale radio observatories.","Because the method isolates the geomagnetic component, the derived electromagnetic energy can be combined with muon counters to probe the primary mass composition of inclined cosmic rays.","Where the signal in the $v\\times v\\times B$ polarization is strong, the charge-excess fraction can be measured directly from the data instead of using the parameterization, making the symmetrization self-calibrating in high-signal regions.","The closed-form, analytically documented nature of the model allows fast, reproducible event reconstruction suitable for real-time or offline pipelines."],"supporting_citations":[{"why":"Defines the early-late correction geometry used in Section 2 and supplies the original exponential fit that this work extends.","marker":"[5]"},{"why":"Provides the polarization decomposition (equation 3.2) used to extract the charge-excess fraction and the pure geomagnetic fluence.","marker":"[3]"},{"why":"Describes the CoREAS Monte Carlo code that generated the 3111 simulation events used to derive and test the model.","marker":"[7]"},{"why":"Introduces the atmospheric-density scaling and radiation-energy calibration underlying equations (6.1) and (6.2).","marker":"[8]"},{"why":"Contributes the canonical exponential-of-a-quadratic lateral profile that is generalized to the exponential-of-a-cubic fit function.","marker":"[4]"}],"fun_headline_variants":["Inclined shower radio maps flattened for <3% cosmic-ray energy spread","Symmetrized radio footprints deliver cosmic-ray energies within 3%","Inclined air showers: radio map symmetrized, cosmic-ray energy pinned to 3%","New symmetrization method trims cosmic-ray energy error to under 3%"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the radio emission from an inclined shower can be modeled as originating from a single point at the shower maximum, and that the charge-excess formula tuned on one simulation library remains valid when applied to real showers or other simulation conditions.","fun_headline_variants_meta":{"raw":{"variants":["Inclined shower radio maps flattened for <3% cosmic-ray energy spread","Symmetrized radio footprints deliver cosmic-ray energies within 3%","Inclined air showers: radio map symmetrized, cosmic-ray energy pinned to 3%","New symmetrization method trims cosmic-ray energy error to under 3%"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000953,"raw_usage":{"total_tokens":4084,"prompt_tokens":981,"completion_tokens":3103,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":597,"completion_tokens_details":{"reasoning_tokens":3015}},"tokens_in":597,"tokens_out":3103,"duration_ms":21501,"temperature":1.0,"reasoning_tokens":3015,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:54:49.738668+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Split the 3111 CoREAS simulations in half: fit the charge-excess parameterization and energy calibration on one half, run the full symmetrization and integration on the other, and check whether the corrected geomagnetic radiation energy still has a spread below 3 percent and no bias. If the spread grows, the claimed universal accuracy comes from in-sample tuning rather than from the model itself.","supporting_citations":[{"cited_title":"Glaser, S","cited_arxiv_id":null,"evidence_quote":"Provides the polarization decomposition (equation 3.2) used to extract the charge-excess fraction and the pure geomagnetic fluence."},{"cited_title":"Huege, M","cited_arxiv_id":null,"evidence_quote":"Describes the CoREAS Monte Carlo code that generated the 3111 simulation events used to derive and test the model."},{"cited_title":"Glaser, M","cited_arxiv_id":null,"evidence_quote":"Introduces the atmospheric-density scaling and radiation-energy calibration underlying equations (6.1) and (6.2)."},{"cited_title":"Kostunin, P.A","cited_arxiv_id":null,"evidence_quote":"Contributes the canonical exponential-of-a-quadratic lateral profile that is generalized to the exponential-of-a-cubic fit function."}],"review_version":1}