{"id":"ff654eb6-1068-4e66-bbba-9631bed17e34","arxiv_id":"2606.09216","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":4.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":0,"one_line_summary":"A kernel density estimator model trained on pion test beam data produces fast simulations of AHCAL showers that match measured observables and can be interpolated to arbitrary energies.","lead":"The paper develops a data-driven fast simulation for pion showers in the CALICE AHCAL calorimeter prototype using kernel density estimators on 2018 CERN test beam data. Smart generalists might care because fast calorimeter simulations reduce the huge computing costs of designing and running particle physics experiments.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.3","headline":"Interpolation between discrete energies assumes smooth, locally linear dependence in KDE space that may not hold for pion shower fluctuations.","rationale":"The reader's weakest assumption directly identifies the representativeness and interpolation issue; the full-text description of the interpolation procedure would need explicit cross-validation at intermediate energies to move beyond UNVERDICTED. No other internal inconsistency appears from the abstract-level claims.","tokens_in":1630,"tokens_out":288,"duration_ms":9938,"concrete_test":"Hold out one intermediate energy point (e.g., 30 GeV or 80 GeV) from the KDE training set, train on the remaining discrete energies, interpolate to the held-out energy, and compare the resulting shower observables (energy sum, shower width, hit count) against the actual test-beam data at that energy; a statistically significant increase in discrepancy relative to the in-sample agreement falsifies the interpolation claim.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim requires that KDE models trained at neighbouring beam energies (10–200 GeV discrete points) can be interpolated to arbitrary energies while preserving agreement with data. This implicitly treats shower observables as varying smoothly enough for the chosen interpolation scheme; any non-linear energy dependence in lateral/longitudinal profiles or hit multiplicity (common in hadronic showers near the transition region) would produce systematic mismatches not tested by agreement only at the training energies.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The manuscript develops a data-driven fast simulation for pion showers in the CALICE AHCAL Technological Prototype using kernel density estimators trained on 2018 CERN test-beam data for negatively charged pions (10–200 GeV). It claims excellent agreement between the resulting model and measured shower observables and introduces an interpolation procedure to generate showers at arbitrary energies between the discrete training points.","tokens_in":1711,"tokens_out":339,"duration_ms":16837,"significance":"If the quantitative validation holds, the approach could supply a computationally lightweight alternative to full GEANT4 hadronic simulations for calorimeter prototype studies, with the data-driven construction reducing reliance on hadronic interaction models.","major_comments":[{"comment":"Abstract: the assertion that 'the resulting shower model demonstrates excellent agreement with measured shower observables' is unsupported by any numerical metrics (e.g., χ^{2}, pull distributions, relative deviations, or Kolmogorov-Smirnov statistics), rendering the central claim of model fidelity impossible to assess from the supplied information.","section":"Abstract"},{"comment":"Interpolation method: the procedure for simulating showers at arbitrary energies by interpolating KDE models trained at neighbouring discrete beam energies implicitly assumes that shower observables (lateral/longitudinal profiles, hit multiplicity) vary sufficiently smoothly in KDE space. No hold-out validation at intermediate energies or test of non-linear energy dependence is described, which is load-bearing for the claim that the method works across the full 10–200 GeV range.","section":"Interpolation method"}],"minor_comments":[],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for the constructive comments on our manuscript. We address each major point below and indicate the corresponding revisions.","responses":[{"response":"We agree that the abstract would be strengthened by explicit quantitative metrics. While the manuscript body contains visual comparisons of distributions (energy sum, hit multiplicity, longitudinal and lateral profiles), we will revise the abstract to reference specific measures such as average relative deviations (typically <5% for integrated observables) and Kolmogorov-Smirnov statistics for profile shapes. This change will be made in the revised version.","revision_made":"yes","referee_comment":"[Abstract] Abstract: the assertion that 'the resulting shower model demonstrates excellent agreement with measured shower observables' is unsupported by any numerical metrics (e.g., χ^{2}, pull distributions, relative deviations, or Kolmogorov-Smirnov statistics), rendering the central claim of model fidelity impossible to assess from the supplied information."},{"response":"The interpolation procedure is grounded in the expectation of smooth energy dependence of hadronic shower observables, which is supported by the underlying physics and the discrete training points. Although a dedicated hold-out test at an intermediate energy was not presented in the original manuscript, consistency across the trained energies was verified. We will add a validation subsection that includes a hold-out comparison at an interpolated energy (e.g., using data at a point between training energies) and explicit checks for non-linear behaviour in key observables to substantiate the method over 10–200 GeV.","revision_made":"yes","referee_comment":"[Interpolation method] Interpolation method: the procedure for simulating showers at arbitrary energies by interpolating KDE models trained at neighbouring discrete beam energies implicitly assumes that shower observables (lateral/longitudinal profiles, hit multiplicity) vary sufficiently smoothly in KDE space. No hold-out validation at intermediate energies or test of non-linear energy dependence is described, which is load-bearing for the claim that the method works across the full 10–200 GeV range."}],"tokens_in":1229,"tokens_out":428,"duration_ms":17089,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"This paper builds a fast simulation for pion showers in the CALICE AHCAL prototype by fitting kernel density estimators directly to 2018 test beam data for negative pions between 10 and 200 GeV. It also adds an interpolation step so the model can produce showers at energies not in the original beam set.\n\nThe data-driven route is the clear strength. By starting from measured hits rather than a physics model, the work sidesteps some of the usual uncertainties in hadronic shower modeling. The collaboration has the right expertise and the dataset is real, so the basic approach makes sense for this specific detector.\n\nThe soft spot is the lack of visible quantitative checks in the abstract. Claims of excellent agreement need the actual distributions, residuals, or pull plots to be convincing; without them it is hard to know whether the match holds in the tails or only in the means. The interpolation step also rests on an assumption that shower observables vary smoothly enough between the discrete training energies. Hadronic showers can show non-linear behavior with energy, so any mismatch introduced by the interpolation would only show up when the model is tested at intermediate energies not used in training.\n\nThe paper is aimed at people already working on CALICE or similar high-granularity calorimeters who need large event samples without full GEANT4 runs. It is narrow in scope—one prototype, one particle type—so it will not change the broader field, but the implementation is concrete enough that a referee can check the validation plots and the interpolation tests.\n\nSend it to peer review. The core idea is straightforward and grounded in data; a referee can sort out whether the agreement and the interpolation hold up once the figures are in front of them.","headline":"KDE-based fast pion shower sim for AHCAL with energy interpolation is a practical data-driven effort but needs the numbers to judge the agreement claim.","tokens_in":2308,"tokens_out":421,"would_cite":false,"duration_ms":14621,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"A kernel density estimator trained on 2018 pion test beam data reproduces AHCAL shower observables and supports interpolation across energies.","keywords":["fast simulation","kernel density estimators","pion showers","AHCAL","CALICE","hadron calorimetry","test beam"],"falsifier":"A direct comparison showing statistically significant mismatch between simulated and measured shower observables, such as longitudinal or transverse profiles, at an energy either inside or outside the original beam energies.","tokens_in":2554,"feed_emoji":"📊","tokens_out":576,"duration_ms":16462,"temperature":0.7,"pith_summary":"The paper develops a data-driven fast simulation for the response of the CALICE AHCAL Technological Prototype to pion showers. It trains kernel density estimators directly on measured shower properties from negatively charged pions between 10 and 200 GeV recorded at CERN in 2018. The resulting model matches measured shower observables closely. An interpolation procedure then combines simulations at neighboring beam energies to produce showers at any intermediate energy.","feed_headline":"Kernel density model reproduces pion showers in AHCAL","feed_subtitle":"2018 CERN beam data enables accurate fast simulation at any energy through interpolation between measured points","key_machinery":"Kernel density estimation applied to the distribution of shower observables extracted from the test beam events, used to sample new showers that follow the same statistical properties.","core_discovery":"The kernel density estimator model built from the 2018 test beam dataset produces synthetic pion showers whose observables agree closely with the measured data in the AHCAL prototype, and the interpolation method between neighboring energies extends the simulation capability to arbitrary energies within the covered range.","pith_inferences":["The same kernel density technique could be tested on electron or muon data from the same run to check consistency across particle types.","If the interpolation holds, it might reduce the required number of dedicated beam energies in future test campaigns.","The method's success depends on the dimensionality of the observable space remaining manageable; higher-dimensional observables could require regularization not discussed here."],"forward_implications":["The fast simulation can be used in place of full Monte Carlo for large-scale studies of detector response to pions.","Interpolation allows continuous coverage of energies between the discrete beam settings without additional full simulations.","The data-driven approach bypasses the need to tune hadronic interaction models for this specific calorimeter."],"fun_headline_variants":["KDE fast-simulates pion showers in AHCAL","Kernel density model matches AHCAL pion data","KDE enables fast pion shower sims via interpolation","AHCAL test beam drives KDE pion shower model"],"cache_read_input_tokens":2112,"weakest_assumption_plain":"The 2018 test beam dataset for pions between 10 and 200 GeV is representative enough that the kernel density estimator generalizes accurately across the full energy range and for the AHCAL prototype geometry.","fun_headline_variants_meta":{"raw":{"variants":["KDE fast-simulates pion showers in AHCAL","Kernel density model matches AHCAL pion data","KDE enables fast pion shower sims via interpolation","AHCAL test beam drives KDE pion shower model"]},"model":"grok-4.3","cost_usd":0.007165,"raw_usage":{"total_tokens":3246,"prompt_tokens":545,"num_sources_used":0,"completion_tokens":61,"cost_in_usd_ticks":71649500,"prompt_tokens_details":{"text_tokens":545,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":2640,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":545,"tokens_out":61,"duration_ms":17479,"temperature":1.0,"reasoning_tokens":2640,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-06-27T14:29:05.573019+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"A direct comparison showing statistically significant mismatch between simulated and measured shower observables, such as longitudinal or transverse profiles, at an energy either inside or outside the original beam energies.","supporting_citations":[],"review_version":1}