Pith. sign in

REVIEW 2 major objections 5 minor 10 references

On the resolution of dual readout calorimeters

T0 review · 2 major / 5 minor · reviewed 2026-08-10 · deepseek-v4-flash

Pith's one-line read One resolution formula predicts when the dual-readout correction improves hadronic calorimeters and when it hurts them.

desk verdict A clean, honest derivation of a dual-readout resolution formula, validated broadly; the covariance check the stress-test wants would be nice but the paper already earns a serious review. read the letter →

arxiv 2501.15329 v3 pith:NL5UNWIL submitted 2025-01-25 physics.ins-det hep-ex

classification physics.ins-dethep-ex PACS 29.40.Vj
keywords dualreadoutcalorimetryhadronicenergyresolutionelectromagneticfractionCherenkovscintillationformulaescapingnuclearbindingloss
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper derives a closed-form formula for the resolution of a dual-readout calorimeter's corrected energy estimate $D$. The formula combines the scintillation resolution $\sigma_S$, the Cherenkov resolution $\sigma_C$, and the shower-to-shower fluctuation $\sigma_f$ of the electromagnetic fraction, explicitly subtracting the part of the two signals that is correlated through $f$. The authors verify the formula with a toy Monte Carlo and with full simulations of five calorimeter geometries, including a poorly contained sampling calorimeter in which the dual-readout correction makes the resolution worse. The result matters because dual-readout calorimeters are leading candidates for future collider calorimetry, and the formula lets designers know in advance whether the correction will help a given geometry or be defeated by uncorrelated noise and escaping energy.

What carries the argument

The load-bearing object is the corrected estimator $D$, a linear projection of the $(C,S)$ pair along their correlation line onto the electron energy scale. The argument uses simple error propagation: because the scintillation and Cherenkov signals both depend on the same electromagnetic fraction $f$, their covariance is $(1-h_S)(1-h_C)\sigma_f^2$; inserting this covariance into the standard error-propagation formula gives Eq. (11). The secondary machinery is the decomposition of the shower into an electromagnetic part, a non-electromagnetic part, and the loss fractions $f_{NC}$ (nuclear binding) and $f_{ES}$ (escaping energy), which yields approximate formulas for $h_S$ and $\sigma_S$ and makes the meaning of the correction concrete.

What would settle it

Run a full shower simulation, compute the energy estimates $S$ and $C$ shower by shower, and measure their covariance directly; the formula requires $\mathrm{cov}(S,C) = (1-h_S)(1-h_C)\sigma_f^2$, so a disagreement larger than the statistical uncertainties would invalidate the error-propagation chain that produces Eq. (11).

Watch

Extended reading notes

Core claim

The paper's central claim is Eq. (11): $$\sigma_D = \frac{1}{h_S-h_C}\sqrt{(1-h_C)^2\$sigma_S^{2}$ + (1-h_S)^2\$sigma_C^{2}$ - 2(1-h_S)^2(1-h_C)^2\$sigma_f^{2}$},$$ which predicts the resolution of the dual-readout-corrected energy estimate $D=\frac{(1-h_C)S-(1-h_S)C}{h_S-h_C}$ for hadrons of fixed energy. Here $h_S$ and $h_C$ are the average responses of the scintillation and Cherenkov readouts to the non-electromagnetic part of the shower, and $\sigma_f$ is the rms fluctuation of the electromagnetic fraction. The paper tests this expression against simulated showers in a homogeneous crystal calorimeter, two fiber calorimeters, and two sampling-tile calorimeters, finding predicted corrected resolutions close to the simulated values for all five geometries, including the SampS case where the correction degrades resolution because escaping energy and large correlated fluctuations dominate. The same derivation shows the correction compensates not only for nuclear binding-energy loss but also for energy escaping the calorimeter and energy missed by the clustering algorithm, and it yields approximate formulas for $h_S$ and $\sigma_S$ in terms of sampling fractions and the loss fractions $f_{NC}$ and $f_{ES}$.

Load-bearing premise

The formula assumes that the electromagnetic fraction $f$ of the shower is roughly Gaussian with moderate width and that $h_S$ and $h_C$ are constants independent of $f$ and of the noise terms; where those assumptions break down, as in the small SampS calorimeter, the quantitative agreement loosens even though the formula still describes the trend.

Editorial extensions

If this is right

  • With measured or estimated $h_S$, $h_C$, $\sigma_S$, $\sigma_C$, and $\sigma_f$, Eq. (11) predicts before full construction whether the dual-readout correction will improve a proposed calorimeter's hadronic resolution.
  • The correction is expected to hurt when the Cherenkov readout is noisy: Eq. (13) states that the dual-readout estimate is worse than the scintillation-only estimate unless $\sigma_f > \sigma_C/(\sqrt{2}(1-h_C))$.
  • Dual-readout compensation applies to scale fluctuations from escaping energy and clustering losses as well as from nuclear binding energy, so containment and algorithm efficiency enter the resolution budget explicitly.
  • For nearly compensating geometries such as Fiber2, the dual-readout correction gives only a modest improvement because the scintillation response already has little dependence on $f$.
  • The formula explains the counterintuitive SampS result from an earlier simulation, predicting the observed worsening of resolution for a small, poorly contained calorimeter.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A practical consequence not spelled out in the paper: Eq. (13) can be used as a go/no-go test during design, comparing expected Cherenkov noise to the electromagnetic-fraction fluctuation before investing in the second readout.
  • The same covariance argument transfers to any two-readout estimator that shares a common latent variable: whenever the correlated fluctuation term is small compared with the uncorrelated noises, the corrective linear combination will degrade rather than improve resolution.
  • The formula suggests a direct simulation diagnostic: measure the empirical shower-by-shower covariance of $S$ and $C$ and compare it with $(1-h_S)(1-h_C)\sigma_f^2$; a large discrepancy would identify correlations beyond the simple Gaussian-$f$ model.
  • For future collider calorimeters, the explicit appearance of $f_{ES}$ in $h_S$ turns containment into a quantitative resolution input, so optimizing length and transverse size becomes an optimization of the resolution formula itself.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 5 minor

Summary. The paper derives a closed-form expression for the energy resolution of a dual-readout calorimeter after the dual-readout correction, Eq. (11), starting from the ansatz that the scintillation and Cherenkov signals share a common dependence on the electromagnetic fraction f, plus independent noise terms. The covariance between S and C is computed in Appendix A, giving cov(S,C) = (1-hS)(1-hC)σf². The formula is tested against a toy Monte Carlo and against GEANT4 simulations of five calorimeter geometries (PbWO, fiber, and sampling calorimeters), with the predicted corrected resolution compared in Table 6. The paper also presents approximate formulae for hS and σS in terms of binding-energy loss, escaping energy, and sampling-fraction fluctuations. The stated purpose is to explain when the dual-readout correction improves or worsens the resolution.

Significance. If the formula is reliable, it provides a simple, design-oriented tool for estimating the performance of dual-readout calorimeters without full simulation, and it clarifies the counterintuitive result reported in Ref. [4] that the correction can degrade resolution. The paper is explicitly pedagogical and is likely to be useful to newcomers to the field. The full GEANT4 comparison across five geometries, including one (SampS) where the correction worsens resolution, is a valuable stress test. The toy simulation code is made publicly available on GitHub, and the derivation in Appendix A is algebraically correct. The main weaknesses are an incorrect inequality in one of the derived conditions, and a verification that tests the covariance model only indirectly through the final σD.

major comments (2)
  1. [Sec. 2, Eq. (14)] The inequality in Eq. (14) is reversed. Starting from Eq. (11), σD > σS is equivalent to σf² < [(2−χ)σS²] / [2(1−hC)(1−hS)] + σC² / [2(1−hC)²], not the '>' condition printed. For a concrete counterexample with hS=0.9, hC=0.6, σS=σC=0.05, and σf=0.3, the RHS of Eq. (14) is 0.0625 while σf²=0.09, so the printed condition says σD > σS, but Eq. (11) gives σD≈0.039 < σS=0.05. The accompanying discussion near Eq. (13) should be re-examined as well, since the condition for the sum of the second and third terms in Eq. (12) to be positive depends on the same inequality direction.
  2. [Sec. 4.4, Table 6] The verification of Eq. (11) against the full simulation is indirect: hS, hC, and σf are all extracted from the same simulated events used to measure σD, and the comparison only checks the final σD. Since the only non-trivial input of Eq. (11) beyond the measured σS and σC is the covariance model of Eq. (10), a direct comparison between the model covariance (1−hS)(1−hC)σf² and the covariance computed from the simulated (S,C) pairs would provide a much more sensitive test of the central assumption. The statement in Sec. 4.3 that correlations between fNC and fES 'do not affect the prediction' of the corrected resolution is not demonstrated; reporting a direct covariance comparison, or at least a quantitative estimate of the neglected contributions, would substantially strengthen the verification claim.
minor comments (5)
  1. [Sec. 2, Eq. (6)] The partial derivatives in Eq. (6) are written as ∂E/∂S and ∂E/∂C, but the quantity being propagated is D; this should be ∂D/∂S and ∂D/∂C.
  2. [Sec. 4.4, Eq. (31)] Equation (31) contains a typographical error: the term (1−<S)2 is missing the closing angle bracket, and the overall layout makes the placement of the denominator ambiguous; please reformat.
  3. [Sec. 4.2, Eq. (20)] The definition fNC = (EB − EES − EI)/EB normalizes the binding-energy loss to the total beam energy, but in Eqs. (19), (22), and (24) the same symbol is used as a fractional loss of the non-EM component. Please clarify the intended normalization and reconcile it with the caption of Fig. 15, which refers to division by the non-EM energy.
  4. [Table 5] In the Fiber2 row, σf is listed as 0.12 ± 0.1; the uncertainty appears to be a factor of 10 too large and is inconsistent with the other entries; likely ±0.01 is intended.
  5. [General] There are several typographical errors, including 'simulatate' in Sec. 4.1, 'the the D resolution' in Sec. 3, 'ansantz' in Sec. 4.2, and the use of 'fN C' and 'fESC' with inconsistent spacing; a careful proofread is recommended.

Circularity Check

1 steps flagged · score 2.0 of 10

No significant circularity: the central resolution formula is independently tested against GEANT4; only the toy-MC check is a closed algebraic loop.

  1. self definitional [Section 3, 'Tests with a toy simulation']
    "We used a toy Monte Carlo simulation [5] to test Eq. 11. The toy models f as a Gaussian. The S and C signals are calculated from f using Eqs. 1 and 2. The S and C signals are then smeared using separate Gaussians to simulate sources of resolution uncorrelated between the two."

    Eq. 11 is derived in Sec. 2 by propagating errors through exactly Eqs. 1 and 2, with f Gaussian and independent noise terms x2, x3. The toy simulation generates S and C from those same equations and Gaussian smearing, so the measured toy sigma_D is the same quantity that Eq. 11 computes analytically. The agreement is therefore guaranteed by construction up to Monte Carlo statistics; the toy is a check of the algebra, not an empirical test of the physical ansatz. This loop is not load-bearing because the paper also compares Eq. 11 against independent GEANT4 full simulation results in Sec. 4.4, where the parameters are measured from simulated showers rather than generated from the formula.

full rationale

The central claim of the paper is the closed-form resolution formula Eq. 11, obtained by error propagation from the ansatz Eqs. 1-2 with the covariance cov(S,C) = (1-hS)(1-hC)sigma_f^2. The derivation itself is a standard mathematical calculation, not a fit. The toy simulation in Sec. 3 is self-referential: it generates S and C from the same equations used to derive Eq. 11, so it can only verify algebraic consistency. However, the full GEANT4 simulations in Sec. 4 provide an independent test: the calorimeter response is simulated from particle physics, not from Eqs. 1-2, and the inputs hS, hC, sigma_f, sigma_S, sigma_C are extracted from the simulated distributions. Table 6 then compares the predicted and measured sigma_D for five different calorimeter geometries, including a case (SampS) where the correction worsens resolution; this is a falsifiable, non-tautological check. The paper explicitly acknowledges the conditions under which the ansatz is violated (Sec. 4.3, e.g., strong fNC-fES correlation for SampS), and those caveats concern model accuracy rather than circularity. Self-citations to Refs. [3] and [4] are not load-bearing: the relevant SampS and Fiber1 results are re-simulated and quoted directly in this paper. Overall, the derivation is self-contained and the central verification is independent, with only the toy-MC closure and a mild sharing of simulated inputs preventing a fully clean bill of health.

Assumptions & free parameters 5 free parameters · 4 assumptions · 0 invented entities

The central resolution formula itself has no fitted constants: it relates measurable quantities (σS, σC, σf, hS, hC). However, those inputs are detector-specific and in this paper are extracted from the same simulations used to measure σD, so the ledger entries record the empirical parameters and the modeling assumptions required to close the formula.

free parameters (5)
  • hS (average non-EM scintillation response) = 0.740 (PbWO), 1.192 (Fiber1), 1.075 (Fiber2), 1.29 (SampL), 0.45 (SampS)
    Extracted from simulated mean S and mean f via Eq. 29; not directly predicted from geometry except through Eq. 24, which needs fNC, fES, and g-ratio.
  • hC (average non-EM Cherenkov response) = 0.006 (PbWO), 0.21 (Fiber1), 0.14 (Fiber2), 0.08 (SampL), 0.04 (SampS)
    Extracted from simulated mean C and mean f via Eq. 30; nonzero values indicate imperfect separation of EM and non-EM shower components.
  • σf and <f> (EM fraction fluctuation and mean) = Per calorimeter: σf 0.086-0.15, <f> 0.503-0.669 (Tables 4-5)
    Measured from the simulated f distributions (Fig. 14); the resolution formula depends on σf through the negative cross-term.
  • g1, g2 (EM/non-EM sampling fraction ratios) = Per calorimeter, Table 3; e.g., Fiber1 g2 = 2.09 ± 0.06
    Fitted from linear fits to sampling fraction versus f (Fig. 11); used in the hS prediction of Eq. 24.
  • σg/g (sampling fraction fluctuation) = Per calorimeter, Table 8 (e.g., Fiber1 0.0515)
    Used in the σS prediction, Eq. 28; fitted from the spread of sampling fractions.
assumptions (4)
  • domain assumption The dual-readout response is linear in the electromagnetic fraction f: S = E(f + (1−f)hS) and C = E(f + (1−f)hC) with constant hS, hC (Eqs. 1-2).
    This is the standard dual readout model from Ref. [1], invoked in Sec. 2 as the basis of the covariance calculation.
  • domain assumption f is Gaussian-distributed with moderate width.
    Sec. 2 states the formula works best when f is somewhat Gaussian; Appendix A integrates Gaussian g1 to obtain cov(S,C) = H σf^2.
  • domain assumption Noise terms in S and C are mutually independent and independent of f.
    Eq. 9 models noise as independent Gaussians g2, g3, which makes the covariance reduce to the f-only term.
  • ad hoc to paper The ansatz S = g1 fγπ EB + g2[(1−fRπ)EB − (1−fRγ)fγπ EB](1−fNC−fES) (Eq. 19)
    This decomposition into EM and non-EM parts with separate sampling fractions and loss fractions is postulated in Sec. 4.2 to derive hS Eq. 24; the paper notes it requires uncorrelated g2, f, fNC, fES.

how reviews work

0 comments
Cite this review

Pith. "Pith review of On the resolution of dual readout calorimeters." pith.science (2026). https://pith.science/paper/NL5UNWIL

@misc{pith2026250115329,
  author       = {Pith},
  title        = {Pith review of: On the resolution of dual readout calorimeters},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/NL5UNWIL}},
  note         = {Machine review of arXiv:2501.15329}
}
read the original abstract

Dual readout calorimeters allow state-of-the-art resolutions for hadronic energy measurements. Their various incarnations are leading candidates for the calorimeter systems for future colliders. In this paper, we present a simple formula for the resolution of a dual readout calorimeter, which we verify with a toy simulation and with full simulation results. This formula can help those new to dual readout calorimetry understand its strengths and limitations. The paper also highlights that the dual readout correction works not just to compensate for binding energy loss, but also for energies escaping the calorimeter or clustering algorithm. Formulae are also presented for approximate resolutions and energy scales in terms of different sources of response.

Figures

Figures reproduced from arXiv: 2501.15329 by the authors.

Figure 1
Figure 1. [left] Inspired by Ref. [1], an illustration of how fluctuations in the subcom￾ponents of a shower can dominate the resolutions of hadronic calorimeters. The two Gaussians represent the typical measured signal for a fixed incident particle energy for two different values of the electromagnetic fraction. [right] An illustration of how two measurements, one with sensitivity to all components of the shower (“S”) and an… view at source ↗
Figure 2
Figure 2. [left] Cherenkov versus scintillation signals from the toy MC using the parameters [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Energy resolution from the dual-readout estimate divided by that from the [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figures from the paper (15 more)
Figure 4
Figure 4. Figure 4: Shower from a 20 GeV charged pion for crystal calorimeter PbWO. The colors for [PITH_FULL_IMAGE:figures/full_fig_p009_4.png]
Figure 5
Figure 5. Figure 5: Shower from a 20 GeV charged pion for Fiber calorimeter Fiber1. The colors for [PITH_FULL_IMAGE:figures/full_fig_p010_5.png]
Figure 6
Figure 6. Figure 6: Shower from a 20 GeV charged pion for Tile calorimeter SampS. The colors for [PITH_FULL_IMAGE:figures/full_fig_p011_6.png]
Figure 7
Figure 7. Figure 7: Fraction of deposited ionizing energy in 20 GeV electron showers deposited by [PITH_FULL_IMAGE:figures/full_fig_p012_7.png]
Figure 8
Figure 8. Figure 8: Cherenkov versus scintillation signal for calorimeters [top] PbWO, [middle] [PITH_FULL_IMAGE:figures/full_fig_p016_8.png]
Figure 9
Figure 9. Figure 9: Cherenkov (green) and scintillation (red) signal versus number of nuclear inter [PITH_FULL_IMAGE:figures/full_fig_p017_9.png]
Figure 10
Figure 10. Figure 10: For calorimeter Fiber 1, [left] the distribution of the times of energy deposits for [PITH_FULL_IMAGE:figures/full_fig_p017_10.png]
Figure 11
Figure 11. Figure 11: Sampling fraction versus f for the scintillating (red) and Cherenkov (green) active media for calorimeters [top, left] Fiber1, [top, right] Fiber2, [bottom, left] SampL, [bottom, right] SampS. . 18 [PITH_FULL_IMAGE:figures/full_fig_p018_11.png]
Figure 12
Figure 12. Figure 12: Scintillation (red) and Cherenkov (green) signals versus the EM object fraction [PITH_FULL_IMAGE:figures/full_fig_p020_12.png]
Figure 13
Figure 13. Figure 13: Scintillation, Cherenkov, and dual-readout corrected energy distributions for [PITH_FULL_IMAGE:figures/full_fig_p021_13.png]
Figure 14
Figure 14. Figure 14: Shower EM object fraction for pion showers for the various calorimeters. [PITH_FULL_IMAGE:figures/full_fig_p022_14.png]
Figure 15
Figure 15. Figure 15: [left] Total energy loss to binding energies etc. for calorimeters divided by the [PITH_FULL_IMAGE:figures/full_fig_p022_15.png]
Figure 16
Figure 16. Figure 16: The fraction of non-EM energy utilized for overcoming nuclear binding energies [PITH_FULL_IMAGE:figures/full_fig_p023_16.png]
Figure 17
Figure 17. Figure 17: The fraction of non-EM energy escaping the calorimeter [PITH_FULL_IMAGE:figures/full_fig_p024_17.png]
Figure 18
Figure 18. Figure 18: [left] The fraction of non-EM energy going to overcoming nuclear binding [PITH_FULL_IMAGE:figures/full_fig_p025_18.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

10 extracted references · 7 canonical work pages

  1. [4]

    Chekanov, S

    S. Chekanov, S. Eno, S. Magill, C. Palmer, L. Wu, M. Aamir, Geant4 simulations of sampling and homogeneous hadronic calorimeters with dual readout for future colliders, ”Nucl. Instrum. Meth. A” 1072 (2025) 170200. doi:https://doi.org/10.1016/j..2025.170200

  2. [1]

    S. Lee, M. Livan, R. Wigmans, Dual-readout calorimetry, Rev. Mod. Phys. 90 (2) (2018) 025002. arXiv:1712.05494, doi:10.1103/ RevModPhys.90.025002

  3. [2]

    Aly, et al., First test-beam results obtained with IDEA, a detector concept designed for future lepton colliders, Nucl

    R. Aly, et al., First test-beam results obtained with IDEA, a detector concept designed for future lepton colliders, Nucl. Instrum. Meth. A 958 (2020) 162088, proceedings of the Vienna Conference on Instrumenta- tion 2019. doi:https://doi.org/10.1016/j.nima.2019.04.042

  4. [3]

    M. T. Lucchini, C. G. Tully, W. Chung, S. C. Eno, Y. Lai, L. Lucchini, M.-T. Nguyen, New perspectives on segmented crystal calorimeters for future colliders, JINST 15 (11) (2020) P11005.arXiv:2008.00338, doi: 10.1088/1748-0221/15/11/P11005

  5. [5]

    URL https://github.com/saraheno/DualReadoutToy

    github repository for toy (2024). URL https://github.com/saraheno/DualReadoutToy

  6. [6]

    Agostinelli, et al., GEANT4 – a simulation toolkit, Nucl

    S. Agostinelli, et al., GEANT4 – a simulation toolkit, Nucl. Instrum. Meth. A 506 (2003) 250–303. doi:10.1016/S0168-9002(03)01368-8

  7. [7]

    URL https://dd4hep.web.cern.ch/dd4hep/

    Detector description toolkit for high energy physics (2024). URL https://dd4hep.web.cern.ch/dd4hep/

  8. [8]

    Abdullin, et al., Design, performance, and calibration of CMS hadron- barrel calorimeter wedges, Eur

    S. Abdullin, et al., Design, performance, and calibration of CMS hadron- barrel calorimeter wedges, Eur. Phys. J. C 55 (2008) 159. doi:10.1140/ epjc/s10052-008-0573-y

Show all 10 references
  1. [9]

    Bedeschi, A detector concept proposal for a circular e+e− collider, PoS ICHEP2020 (2021) 819

    F. Bedeschi, A detector concept proposal for a circular e+e− collider, PoS ICHEP2020 (2021) 819. doi:10.22323/1.390.0819

  2. [10]

    J. B. Birks, The theory and practice of scintillation counting, Interna- tional Series of Monographs on Electronics and Instrumentation, Vol- ume 27, Pergamon Press, The Macmillan Company, New York, 1964. doi:10.1016/C2013-0-01791-4 . 30 Appendix A. Derivation of the Covarianc...

Pith tools

Reviewed August 10, 2026 · model on record in the stance chip above.