{"id":"7ec43a33-019e-4a2f-8e9a-9ac68506d196","arxiv_id":"2607.14365","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A first analytical solution for dual-mirror Cherenkov-telescope muon-ring calibration is derived, including secondary-mirror shadowing and mirror-curvature corrections.","lead":"This paper derives the first complete analytical formulas for calibrating dual-mirror Cherenkov telescopes with muon rings, including the shadow cast by the secondary mirror and its baffles. The new formulas matter because previous dual-mirror calibrations had to ignore or approximate this shadowing, and the corrections reach up to 40% for inclined muons.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The dual-mirror shadowing predictions rely on an unvalidated flat-disk approximation of the curved M2/baffles; without ray tracing against actual SCT/SST optics, the ±40% claim and 'directly applicable to CTAO' are not established.","rationale":"The paper's central claim is that a complete analytical description of dual-mirror muon-ring shadowing has been derived, with deviations up to ±40% relative to previous hole-based algorithms, and that it is directly applicable to CTAO's SCT/SST telescopes. For this to hold, the geometric model of the secondary mirror, baffles, and primary mirror must be a faithful representation of the actual CTAO optics. The weakest point is that §4 explicitly approximates the entire M2 assembly as two flat circular disks at fixed heights, ignoring the curvature of the aspheric secondary (sagitta ~0.5 m for the SCT) and the real structure of the baffles/support. The paper does not quantify the error introduced by this simplification, and the 40% figure is a comparison between two analytic approximations, not a validation against the real telescope. The internal checks (recovering Vacanti's flat-mirror solution, reproducing Seidel coma) demonstrate algebraic consistency but do not constrain the dual-mirror shadowing geometry. Therefore the central quantitative claim is not yet established for actual CTAO hardware. The concrete test is a ray-tracing comparison with the true optical surfaces; this would settle whether the simplified geometry is adequate. Consequently the verdict should remain CONDITIONAL until that check is performed.","tokens_in":21244,"tokens_out":8123,"duration_ms":83185,"concrete_test":"Run a ray-tracing simulation (e.g., sim_telarray/ROBAST) with the actual SCT and SST optical prescriptions, including the curved aspheric secondary and the real baffle/support geometry. Generate muons with the same parameters as Figs. 9/10 (ν = 0°, 3°, ψ = 0°, 180°, ρ/R1 in [0.1, 0.95]) and record, for each photon emission angle φ, the maximum and minimum emission heights l for which light reaches the camera. Compare these measured L2,max/L2,min and shadow fractions with Eqs. 27/34/35/40 and the last two rows of Figs. 9/10. If the boundary heights differ by more than 5% or the shadow fraction by more than 2 percentage points, the flat-disk approximation is insufficient for the claimed CTAO applicability.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"Section 4 reduces the entire secondary-mirror assembly to 'two flat circular disks of radius R2' at heights Des and Dpb. The real SCT/SST secondary is a strongly curved asphere (for the SCT, Fs = 3.35 m gives sagitta ≈ R2²/(2·2Fs) ≈ 0.55 m), and the baffles/support structure are not flat disks; the paper itself marks Des as 'approximate average of a structured surface' and Rsb as 'best guess' (Table 1). Since the shadow boundary is determined by the height at which the incoming photon crosses the obstruction, a 0.5 m height error translates into an emission-height error Δl ≈ 0.5 m/θc ≈ 17 m for θc = 1.3°, i.e., several percent of the total track length Lmax ≈ 400 m. The paper's headline 40% deviation is computed relative to a different simplified model (M2 as a hole in the primary), not against the actual CTAO optics, and the flat-disk simplification is never quantified. Because the central contribution of the dual-mirror section is precisely these shadowing predictions, the 'complete analytical description ... directly applicable to CTAO' claim inherits the unvalidated geometry. Independent support exists (Vacanti limit, Seidel coma), but it does not constrain the dual-mirror shadowing step.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives closed-form analytical expressions for the Cherenkov light produced by atmospheric muons and imaged by dual-mirror IACTs, using a vector-geometry formalism with SageMath-based symbolic manipulation and Taylor expansion. The authors recover the Vacanti flat-mirror solution as a limiting case, derive a first-order correction to the maximum emission height due to primary-mirror curvature, characterize third-order coma aberration, and provide shadowing formulas for the secondary mirror, baffles, and central hole. They report that for inclined muons the new dual-mirror shadowing treatment differs by up to ±40% from earlier algorithms that treated M2 as a hole in the primary. The paper claims this is the first complete analytical description of dual-mirror muon-ring calibration and states the results are directly applicable to CTAO.","tokens_in":21555,"tokens_out":3567,"duration_ms":38459,"significance":"If the central result is correct, this is a valuable contribution to CTAO calibration: it replaces numerical or approximate treatments of dual-mirror shadowing with explicit formulas, provides a physically transparent derivation, and ships reproducible Python/Jupyter code (GitHub link in §5). The recovery of Vacanti's formula (Eq. 18) and the consistency of the coma term with third-order Seidel expectations (Eq. 14) are strong internal checks. However, the headline application relies on a simplified geometric model of the secondary-mirror assembly that has not yet been validated against the real SCT/SST optics; this gap controls the recommendation.","major_comments":[{"comment":"The dual-mirror shadowing calculation replaces the secondary mirror and its baffles with 'two flat circular disks of radius R2' at heights Des and Dpb. Table 1 itself describes Des as an 'approximate average of a structured surface' and Rsb as a 'best guess'. For the SCT, the secondary is a strongly curved asphere; using Fs=3.35 m gives a sagitta ≈0.55 m at the disk radius. Since the shadow boundary is set by the height at which a photon crosses the obstruction, a 0.5 m height error maps to an emission-height error of ≈0.5 m/θc ≈ 17 m for θc=1.3°, i.e. several percent of the total track length (~400 m). This unquantified simplification is load-bearing for the central claim: the shadow fractions, L2,max, L2,min, and the ±40% deviations in Figs. 9–10 are all predictions for two flat disks, not for the real CTAO secondary/baffle geometry. Please add a quantitative comparison against ray tra","section":"§4, Eqs. (25)–(44) and Table 1"},{"comment":"The abstract and §5 state that this is 'a complete analytical description' that is 'directly applicable to muon-based calibration of CTAO.' This overstates what has been demonstrated. The derivation is exact only for the idealized geometry defined in §4 (flat circular obstruction disks, parabolic primary approximation z≈cρ², and approximate Des/Rsb values). The paper presents no validation of the dual-mirror shadowing step against ray tracing, Geant4/CODE simulations, or real SCT/SST optical data. The Vacanti-limit recovery and the Seidel-coma check constrain only the single-mirror part, not the dual-mirror shadow fractions. The ±40% figure is a deviation between two simplified analytical models, not a measured or ray-traced discrepancy. Please add a simulation/data comparison, or rephrase the claims to describe the result as an analytically solved simplified model that requires a dedica","section":"§5 (Conclusions)"},{"comment":"The abstract says the description accounts for 'the secondary mirror and the camera.' The camera shadow, however, is treated only in Appendix A for a square camera and is not integrated into the main SCT/SST flowcharts (Figs. 8–10) or the headline shadow fractions. The main text uses camera-body parameters from Table 1 (Rcam, Dcam) but these are not used in the central derivation. This is a mismatch between the claim of completeness and the actual scope of the calculation. Either integrate camera shadowing into the main results or clarify that the camera is treated separately in an appendix.","section":"§4 and Appendix A"}],"minor_comments":[{"comment":"The notation 'h.c.' is used repeatedly (e.g., Eqs. 5b, 17e) but never defined. In this context it appears to mean 'higher-order corrections', not 'Hermitian conjugate'; please define it explicitly at first use.","section":"Throughout"},{"comment":"The two branches of Eq. (40) for Dc rely on the sign of ρ−Rsb. The second branch (ρ > Rsb) uses a negative sign before the square root; please state explicitly which root is physical and confirm the sign convention matches the geometry. A brief explanatory sentence would help.","section":"§4, Eq. (40)"},{"comment":"Equation (16) contains ratios γy/γx and γx/γy. These are undefined when the respective component vanishes (e.g., along the azimuthal axes). Please add a note on the limiting behavior or express the result in a form that avoids division by zero.","section":"§3, Eq. (16)"},{"comment":"The figures are informative but the color-coding makes several curves degenerate (as the caption admits). Consider labeling curves directly or using line styles rather than color alone, especially since some panels have overlapping families.","section":"§3, Fig. 2 and Eq. (23)"},{"comment":"Typos: 'the the' in the paragraph after Eq. (23) and 'opposite opposite' in the same area; 'expression were expanded' in §1. These are minor but should be corrected in a revised version.","section":"§1"}],"recommendation":"major_revision","confidential_remarks":"The paper is technically interesting and the symbolic derivation is well structured, with meaningful internal consistency checks. The main concern is not mathematical correctness of the single-mirror part but the external validity of the dual-mirror shadowing step. The flat-disk approximation is acknowledged in the text but its severity is not quantified, and the claim of direct applicability to CTAO is not supported without a ray-tracing or simulation validation. I would encourage the authors to add such a validation or to reframe the contribution as an analytical model of a simplified dual-mirror geometry. This is a fixable issue within the paper's scope, hence major_revision rather than reject."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThe short version: this is real new work, not a repackaging. The authors derive the first analytic dual-mirror muon-ring formalism, recover Vacanti's flat-mirror solution as a limit, and match third-order Seidel coma. The algebra is internally consistent and the code is public. What should make you hesitate is the dual-mirror shadowing model: the secondary and baffles are treated as two flat disks at heights that the authors themselves call approximate or 'best guess,' and the headline ±40% is relative to the previous hole-based approximation, not to ray tracing. So the derivation is likely correct for the simplified geometry, but the claim to be 'directly applicable to CTAO' is not yet established.\n\nThe genuinely new pieces: closed-form shadow conditions for L2,max and L2,min, baffle-crossing cases, central-hole losses, a first-order curvature correction to Lmax, and a coma-bias estimate for ring radii. The Vacanti recovery and Seidel match are strong checks. The paper has no free fitting parameters and the release of SageMath/Python code is a plus.\n\nThe soft spot is proportionate. Flat disks are a reasonable first-order model, and the paper is transparent about it. But the real SCT/SST secondaries are curved aspherics with struts and ring supports, and the table values for Des and Rsb carry sizeable uncertainty. Because a small height error maps to tens of meters in emission height at θc ~ 1.3°, the shadow fractions could shift by several percent. A comparison against full ray tracing—sim_telarray or a good optical model—would settle this. Without it, the dual-mirror predictions are exactly that: predictions for a toy geometry, with unknown error for the actual telescopes.\n\nThe single-mirror results (curvature correction, coma bias) are the most immediately usable part. For the dual-mirror claims, treat as promising but provisional.\n\nWho gets value: the CTAO muon calibration working group, and anyone fitting muon-ring light distributions in dual-mirror IACTs.\n\nRecommendation: send it to peer review. A good referee should verify the vector algebra and ask for the ray-tracing comparison and a sensitivity study on Des and Rsb. With those additions the paper could be a reference for the field.","headline":"Genuinely new analytic dual-mirror muon-ring formalism with strong internal checks; the dual-mirror shadowing claims still need ray-tracing validation before 'directly applicable to CTAO' is justified.","tokens_in":22077,"tokens_out":2547,"would_cite":true,"duration_ms":26456,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A complete analytical description of muon-ring Cherenkov light for dual-mirror telescopes is now derived, accounting for secondary-mirror and baffle shadows and correcting prior flat-mirror approximations by up to 40% for inclined muons.","keywords":["muon rings","IACT calibration","dual-mirror telescopes","Cherenkov light","shadowing","vector geometry","mirror curvature correction","Schwarzschild-Couder"],"falsifier":"A ray-tracing Monte Carlo that uses the detailed mechanical model of an actual dual-mirror telescope—curved secondary mirror, struts, and real baffle geometry—and simulates inclined muons over the same parameter ranges would settle the claim: if the simulated per-emission-angle shadowed light fraction differs from the analytical expressions by more than a few percent, the flat-disk shadow model is falsified; agreement at the claimed sub-percent precision would confirm it.","tokens_in":21128,"feed_emoji":"🔭","tokens_out":5250,"duration_ms":54577,"temperature":0.7,"pith_summary":"This paper tries to establish that the light distribution of muon rings in dual-mirror gamma-ray telescopes can be predicted analytically, not just for the flat single-mirror case solved decades ago but for the full dual-mirror geometry with shadowing by the secondary mirror, its baffles, and the central hole. The authors build the solution with vector geometry and expand it in the small parameters—Cherenkov angle, muon inclination, and mirror curvature—to better than 1 percent. Their central result is that previous algorithms, which treated the secondary mirror as a large hole in the primary, are off by up to 40 percent for inclined muons. This matters because muon rings are the main continuous calibration source for monitoring the optical throughput of these telescopes, and a 5 percent optical calibration accuracy is required for the observatory's energy-scale goal.","feed_headline":"Dual-mirror muon calibration is now fully analytic, shadows included","feed_subtitle":"A vector-geometry solution replaces the old hole approximation, which misses inclined muons by up to 40 percent.","key_machinery":"The central machinery is a vector-geometry construction of the muon track and the emitted Cherenkov photon direction, built by rotating the photon direction around the muon velocity vector with a standard rotation formula, then intersecting the photon ray with the primary mirror surface and with two flat circular disks representing the secondary-mirror housing and baffles at two fixed heights. Analytic expansions in the small angles and curvature yield explicit shadow conditions and closed-form expressions for the shadow boundaries, L2,max and L2,min, which directly determine how much ring light is lost for each photon emission angle.","core_discovery":"The paper claims the first complete analytical description of Cherenkov light emitted by atmospheric muons and detected by a dual-mirror imaging Cherenkov telescope. Using vector geometry plus symbolic algebra, it derives closed-form expressions for the photon impact points, the shadow conditions produced by the secondary mirror and its baffles, and the resulting maximum and minimum emission heights that contribute to the recorded ring. The formalism is validated by reproducing the known flat-mirror, non-inclined solution in the appropriate limit and by recovering the expected third-order coma aberration for a parabolic mirror. The largest new numerical effect is that inclined muons experien","pith_inferences":["A natural extension the paper leaves implicit: the same derivation could be repeated for more realistic aspherical primary surfaces beyond the parabolic approximation, yielding higher-order corrections that matter for the actual Schwarzschild-Couder mirrors.","The flat-disk idealization of the secondary assembly is the most likely place where real hardware will differ; a direct comparison of the analytical shadow fractions with ray traced the actual curved secondary, struts, and baffle boxes would test how much of the 40 percent effect survives real geometry.","For telescopes without baffles, such as the current small-telescope design, the shadowing model simplifies considerably, so measurements from existing dual-mirror telescopes could validate the inclination dependence with relatively clean data.","The curvature correction to the maximum emission height may also improve muon calibration of single-mirror Davies-Cotton telescopes, not just the dual-mirror designs emphasized here."],"forward_implications":["Muon-ring calibration for dual-mirror telescopes can be done analytically, without Monte Carlo ray tracing or the common shortcut of treating the secondary mirror as a simple central hole.","Inclined muons, previously neglected, change the predicted ring light by up to 40 percent; including them is necessary for reaching sub-percent-level calibration in the new telescope designs.","The first-order mirror-curvature correction to the maximum emission height, though about 1 percent, enters directly into throughput calibration and should be included in existing reconstruction algorithms.","Coma aberration creates a small bias in the reconstructed muon ring radius in single-mirror parabolic telescopes, averaging near 1 percent and rising to several percent for extreme impact parameters unless plate-scale corrections are applied.","The framework gives a unified way to compute the shadowed and unshadowed contributions to the ring as a function of emission angle, which is needed for interpreting measured muon images in the next-generation dual-mirror instruments."],"fun_headline_variants":["Analytic dual-mirror muon calibration, shadows included","Vector geometry cracks muon ring calibration for dual-mirror","First analytic muon rings with shadows for dual-mirror IACTs","Muon calibration gets analytic shadows: 40% deviations captured"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The entire dual-mirror shadow calculation rests on representing the secondary-mirror assembly as two flat circular disks at fixed heights, so if the real curved secondary mirror, support struts, or finite baffles cast a different shadow, the predicted shadow fractions and the 40 percent deviations describe an idealized telescope rather than the built one.","fun_headline_variants_meta":{"raw":{"variants":["Analytic dual-mirror muon calibration, shadows included","Vector geometry cracks muon ring calibration for dual-mirror","First analytic muon rings with shadows for dual-mirror IACTs","Muon calibration gets analytic shadows: 40% deviations captured"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000262,"raw_usage":{"total_tokens":1482,"prompt_tokens":840,"completion_tokens":642,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":584,"completion_tokens_details":{"reasoning_tokens":568}},"tokens_in":584,"tokens_out":642,"duration_ms":6474,"temperature":1.0,"reasoning_tokens":568,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T02:17:00.950161+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A ray-tracing Monte Carlo that uses the detailed mechanical model of an actual dual-mirror telescope—curved secondary mirror, struts, and real baffle geometry—and simulates inclined muons over the same parameter ranges would settle the claim: if the simulated per-emission-angle shadowed light fraction differs from the analytical expressions by more than a few percent, the flat-disk shadow model is falsified; agreement at the claimed sub-percent precision would confirm it.","supporting_citations":[],"review_version":1}