REVIEW 3 major objections 6 minor 26 references
Comparative study and optimization of SDHCAL hadronic energy reconstruction methods
T0 review · 3 major / 6 minor · reviewed 2026-08-02 · deepseek-v4-flash
Pith's one-line read Split and polynomial-regression formulas give the best compromise between linearity and resolution for hadronic showers, except at high jet energies where particle-flow confusion favors the simple linear formula.
desk verdict Useful SDHCAL reconstruction study with a real but fixable flaw: the single-hadron ranking is partly in-sample, while the dijet and PerfectPFA results hold up independently. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the reconstruction formula E = αN1 + βN2 + γN3, where N1, N2, N3 are the numbers of calorimeter pads firing above each of three charge thresholds; the linear version fixes α, β, γ as constants, while the quadratic, split, and polynomial variants make the coefficients or add terms depending on total hit count N_hit, introducing saturation behavior. The coefficients are tuned by a chi-squared-like minimization on a simulated sample of contained single-hadron events. The other load-bearing mechanism is the angular correction, which multiplies hit counts by geometric factors 1/sinθ and 1/cosφ to compensate for the path-length variation of inclined tracks in the tilted barre
What would settle it
Take the same four formulas and a fixed detector geometry, then apply them to real test-beam data for 5–80 GeV single hadrons using the same thresholds and angular corrections; if the split method's low-energy resolution improvement over the linear formula does not appear, the simulation-based ranking does not transfer to the physical detector.
Extended reading notes
Core claim
The paper's central claim, stated in its conclusion, is that the split method and polynomial regression provide the best compromise across different energy regimes, combining stability with good resolution. For single neutral hadrons, linearity better than 5% is obtained above 10 GeV, with resolution from about 25% at 5 GeV to 7–8% at 80 GeV. For dijets, the optimal choice depends on the reconstruction quality: under realistic particle-flow confusion at high energies, the linear formula is favored because it is less sensitive to hit-count fluctuations, whereas with perfect truth-based reconstruction the nonlinear formulas are superior across the full energy range. The jet energy resolution r
Load-bearing premise
The entire comparison lives in a Monte Carlo simulation whose digital thresholds and hit multiplicities are tuned to reproduce test-beam data, so the formula ranking is only as reliable as that tuning, particularly at low energies where the paper reports an unresolved residual bias at 5 GeV.
Editorial extensions
If this is right
- Adopting the split or polynomial formula for single-hadron reconstruction improves resolution below 20 GeV without worsening linearity above 10 GeV.
- For jet energies above roughly 100 GeV, the resolution ceiling is set by particle-flow confusion, not by the calorimeter formula; formula choice should be made with the reconstruction algorithm in mind.
- The angular correction removes a multi-percent systematic underestimation in the barrel and should be part of any calibration of a modular hadronic calorimeter.
- If particle-flow confusion were removed (via better clustering or timing), nonlinear formulas would give better jet resolution at all energies, making the linear formula's high-energy advantage a temporary artifact of the current reconstruction.
- A jet energy resolution of 3.4–4% above 100 GeV is achievable with the semi-digital technology, within a point of the 3% benchmark of the alternative scintillator-based design.
Reading between the lines
- The same fitting machinery could take timing information (from a planned timing-capable version of the calorimeter) as an additional input to the polynomial regression, attacking confusion at its source; the paper lists this as a future direction, but the natural extension is to quantify how much confusion the polynomial's cross terms already absorb.
- The reversal of formula ranking at high jet energies is essentially a statement about the bias-variance tradeoff under noise: the linear formula is the highest-bias, lowest-variance estimator, so it wins when hit-count noise is amplified by clustering mistakes. This suggests the same qualitative ranking would appear for any high-granularity calorimeter with a noisy clustering step, not just this o
- A testable extension is to rerun the dijet comparison at the specific jet energies of planned lepton colliders (for example around 100 GeV for ZH events); the paper covers 30–280 GeV, so the formula of choice for the actual physics runs can be pinned down.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper compares four methods for reconstructing hadronic energy in the SDHCAL within the ILD_l2_v02 simulation: linear, quadratic, split, and polynomial regression. It introduces a geometrical angular correction for the barrel, calibrates coefficients on single K_L samples, and evaluates linearity and resolution on single K_L and dijet events using APRIL, with PerfectPFA and PandoraPFA cross-checks. The central claims are: (i) the angular correction is essential; (ii) all methods achieve few-percent linearity; (iii) for single neutral hadrons the split method and polynomial regression give the best compromise; and (iv) for dijets above roughly 100 GeV, PFA confusion dominates and favors the linear parametrization, while non-linear methods are superior when confusion is removed. Full parameter tables are provided.
Significance. The dijet confusion conclusion is well supported by the independent PerfectPFA and PandoraPFA checks, which is a genuine strength. The angular correction is physically motivated and addresses a real geometry effect. If the single-hadron ranking were validated on independent data, the paper would provide useful calibration guidance for the SDHCAL option in ILD. However, the single-hadron 'best compromise' claim currently rests on an in-sample evaluation; the split threshold is tuned on the same K_L sample used for Fig. 8, and no uncertainties are given for the resolution differences. The paper is therefore not yet convincing on its headline conclusion, although the methodology and the dijet analysis are sound.
major comments (3)
- [§3.1, §2.4.3] The single-K_L performance comparison is in-sample. Section 3.1 states that all methods were 'tested on the same samples of single K_L events used to fit the formulas'; Section 2.4.3 describes the split transition threshold as scanned from 200 to 600 hits on the same data, with 'optimal results' chosen. Thus the ranking in Figs. 8 and 9 is not an independent test: the split method's improvement at 5–20 GeV could be an overfit artifact. Only the polynomial regression has a within-method train/test split, which does not protect the cross-method ranking. Please re-evaluate at least the split method on a disjoint holdout sample or with k-fold cross-validation, and quote statistical uncertainties on the ΔJER values in Fig. 9. The dijet analysis is less affected because it applies pre-trained coefficients to an independent topology, but the single-hadron claim in §4 is load-bearing.
- [§2.4, Tables 1–4] Parameter uncertainties are not reported. The text says the coefficients are intended for direct use and their individual uncertainties are 'not relevant,' but the paper's central claim is a comparison of methods. Without uncertainties on the fitted coefficients or on the reconstructed-energy differences, the reader cannot judge whether the small ΔJER differences in Figs. 9 and 11 are statistically significant, even in-sample. Please provide covariance matrices or bootstrap error bars for the key energy points.
- [§2.1, §2.2] The entire study is simulation-based. The digitization thresholds are fixed to reproduce test-beam multiplicity and efficiency (Ref. [21]), but no direct validation of the simulated response against test-beam data is shown in this paper, and the 5 GeV residual in §2.2 is left unresolved. The ranking of the methods, especially at low energies, could change if the simulated multiplicity/threshold response differs from the real detector. Please state this limitation explicitly and, if possible, include a comparison with test-beam data or an alternative digitization for at least one energy point.
minor comments (6)
- [Eq. (2.1)] Define N_i^barrel and N_i^barrel,module; the current notation is not self-explanatory.
- [§2.2] The text mentions the 'standard PandoraPFA calibration procedure' but Fig. 7 uses the quadratic formula; clarify which reconstruction is used in the angular-correction study.
- [§2.4.3] The scan of transition thresholds from 200 to 600 hits is mentioned but not shown; include the scan curve or a table of the tested thresholds and corresponding resolutions.
- [Figs. 8, 9, 11] No error bars are shown; even if statistical uncertainties are small, state the procedure used to compute errors or state that they are negligible.
- [General] Typos/typesetting: 'dijet' appears as a ligature; 'K^0_L' and '𝜙_folded' are inconsistently formatted; some equations lack punctuation.
- [§3.2] The PandoraPFA confirmation is described in the text but no plot is shown; consider adding a reference to an existing figure or a supplementary plot.
Circularity Check
Single-hadron method ranking is in-sample: the calibration coefficients and the split threshold are fitted and evaluated on the same K_L samples; dijet and PerfectPFA parts are independent.
-
fitted input called prediction
[§2.3, §2.4.1-2.4.4, §3.1, Figs. 8-9]
"All the methods described in section 2.4 were first tested on the same samples of single K_L events used to fit the formulas. However, for the performance evaluation presented in this section, the containment and pure-HCAL selection criteria applied to define the calibration dataset are not imposed."
Section 2.3 constructs the calibration set by selecting contained, pure-HCAL single-K_L events and fits all formulas on it; Section 3.1 then evaluates the formulas on the same samples, only dropping the selection cuts and running APRIL. The linearity/resolution in Figs. 8-9 is therefore in-sample. The claim that the split method and polynomial regression 'offer the best compromise' for single hadrons is a description of the fitted sample, not an independent prediction; no cross-validation protects the cross-method ranking.
-
fitted input called prediction
[§2.4.3, §3.1, Fig. 9]
"Various transition thresholds, ranging from 200 to 600 hits, were tested. Ultimately, the optimal results were achieved by applying the function fitted on events with N_hit≤600 to clusters with N_hit≤400, and the function fitted on events with N_hit>400 to clusters with N_hit>400."
The 400-hit split threshold is a hyperparameter chosen by scanning 200-600 hits and keeping the 'optimal results' on the same single-K_L data used for the Fig. 8-9 evaluation. The subsequent conclusion that 'the split method yields improved performance between 5 GeV and 20 GeV' is therefore partly a restatement of that optimization, not an out-of-sample validation. This makes the headline low-energy advantage partially circular by construction.
full rationale
The single-K_L comparison is partially circular: the coefficients are fitted on the calibration subset of these events, and the paper explicitly evaluates the methods on 'the same samples ... used to fit the formulas.' The split method adds a scanned hyperparameter (200-600 hits, 400 chosen as 'optimal') on the same data, so its low-energy resolution advantage is statistically forced to look good on those curves. This does not invalidate the whole paper: the angular correction is an independent geometric correction; the dijet analysis applies pre-trained coefficients to a separate topology and is therefore a genuine transfer test; and the PerfectPFA comparison independently supports the confusion-dominance interpretation. There is no load-bearing self-citation: APRIL, CALICE test-beam data, and ILD benchmarks are external reference points. The 5 GeV residual and lack of uncertainties on ΔJER are limitations, but not circularity. Score 6 reflects a partial circularity confined to the central single-hadron ranking claim.
Assumptions & free parameters
free parameters (7)
- Linear energy factors alpha, beta, gamma =
0.03670, 0.07453, 0.36304 GeV (Table 1)
- Quadratic coefficients alpha1-alpha3, beta1-beta3, gamma1-gamma3 =
Table 2 (9 values; gamma1 bounded at 1e-10 GeV)
- Split-method coefficients (18 values) =
Table 3 (two sets of 9)
- Split transition thresholds =
N_hit = 400 (application), 600 (fit range)
- Polynomial regression coefficients (9 values) =
Table 4
- Chi-squared event weight sigma_i = sqrt(E_mc)
- Resolution-extraction window =
[mu1 +/- 1.5 sigma1] from first Gaussian fit
assumptions (5)
- domain assumption GEANT4 with QGSP_BERT physics list models hadronic showers in the SDHCAL steel/GRPC structure adequately
- domain assumption Digitization thresholds (114 fC, 6.12 pC, 16.83 pC) reproduce prototype multiplicity and efficiency
- domain assumption Geometric angular correction Eq. (2.1): hit deficit scales as 1/sin(theta) and 1/cos(phi) per module
- domain assumption Cluster-level application of event-level-calibrated formulas is valid
- domain assumption Chi-squared weighting with sigma = sqrt(E_mc) is appropriate
Cite this review
Pith. "Pith review of Comparative study and optimization of SDHCAL hadronic energy reconstruction methods." pith.science (2026). https://pith.science/paper/KZE6VMQ2
@misc{pith2026260715023,
author = {Pith},
title = {Pith review of: Comparative study and optimization of SDHCAL hadronic energy reconstruction methods},
year = {2026},
howpublished = {\url{https://pith.science/paper/KZE6VMQ2}},
note = {Machine review of arXiv:2607.15023}
}
abstract
We present a detailed study of hadronic shower energy reconstruction methods for the Semi-Digital Hadronic Calorimeter (SDHCAL) within the ILD detector concept, using the Particle Flow Algorithm (PFA) APRIL. Using samples of single $K^0_L$ and dijet ($u,d,s$) events, we compare linear, quadratic, split, and polynomial regression-based reconstruction formulas, focusing on their impact on linearity and resolution. The study also addresses angular corrections required in the barrel region due to non-perpendicular particle incidence. Results show that while all methods achieve good overall performance, the split method and the polynomial regression provide the best compromise across different energy regimes, offering improved resolution at low energies without compromising linearity at higher energies. For dijets, sensitivity to PFA confusion dominates the resolution at high energies. These findings highlight the potential of future improvements, notably the integration of precise timing information from the T-SDHCAL into APRIL, to further reduce confusion and enhance hadronic energy reconstruction for next-generation lepton colliders.
Reference graph
Works this paper leans on
-
[21]
G. Garillot,Étude des gerbes hadroniques dans un calorimètre à grande granularité et étude du canal e+e-→ HZ (Z→qq) dans les futurs collisionneurs leptoniques, theses, Université de Lyon, Feb., 2019
2019
-
[1]
T. Behnke, J.E. Brau, B. Foster, J. Fuster, M. Harrison, J.M. Paterson et al.,The International Linear Collider Technical Design Report - Volume 1: Executive Summary, 1306.6327
-
[2]
CEPC Study Group collaboration, CEPC Technical Design Report: Accelerator, Radiat. Detect. Technol. Methods8 (2024) 1 [2312.14363]
arXiv 2024
-
[3]
Abada, M
A. Abada, M. Abbrescia, S.S. AbdusSalam, I. Abdyukhanov, J. Abelleira Fernandez, A. Abramov et al.,FCC-ee: The Lepton Collider: Future Circular Collider Conceptual Design Report Volume 2, The European Physical Journal Special Topics228(2019) 261–623
2019
-
[4]
J.-C. Brient and H. Videau,The Calorimetry at the future e+ e- linear collider, eConf C010630 (2001) E3047 [hep-ex/0202004]
arXiv 2001
-
[5]
Thomson,Particle flow calorimetry and the PandoraPFA algorithm,Nucl
M. Thomson,Particle flow calorimetry and the PandoraPFA algorithm,Nucl. Instrum. Meth. A611 (2009) 25–40
2009
-
[6]
M. Ruan and H. Videau,Arbor, a new approach of the Particle Flow Algorithm, inInternational Conference on Calorimetry for the High Energy Frontier, pp. 316–324, 2013 [1403.4784]
arXiv 2013
-
[7]
B. Li, R. Été, G. Grenier and I. Laktineh,APRIL: a novel Algorithm for Particle Reconstruction at ILC,Journal of Instrumentation15(2020) C05016–C05016
2020
Show all 26 references
-
[8]
Baulieu, M
G. Baulieu, M. Bedjidian, K. Belkadhi, J. Berenguer, V. Boudry, P. Calabria et al.,Construction and commissioning of a technological prototype of a high-granularity semi-digital hadronic calorimeter, Journal of Instrumentation10(2015) P10039
2015
-
[9]
The CALICE collaboration,First results of the CALICE SDHCAL technological prototype,Journal of Instrumentation 11(2016) P04001. – 18 –
2016
-
[10]
CALICEcollaboration,Calorimetry for Lepton Collider Experiments - CALICE results and activities, 1212.5127
-
[11]
Abramowicz et al.,The International Linear Collider Technical Design Report - Volume 4: Detectors, 1306.6329
H. Abramowicz et al.,The International Linear Collider Technical Design Report - Volume 4: Detectors, 1306.6329
-
[12]
ILD Concept Group collaboration, International Large Detector: Interim Design Report, 2003.01116
2003 arXiv
-
[13]
S. Lu, F. Gaede, A. Sailer, D. Protopopescu, D. Jeans, A.T. Delgado et al.,key4hep/k4geo: v00-23, Oct., 2025. 10.5281/zenodo.17249123
2025 doi
-
[14]
Callier, J.B
S. Callier, J.B. Cizel, F. Dulucq, C.d.L. Taille, G. Martin-Chassard and N. Seguin-Moreau,ROC chips for imaging calorimetry at the International Linear Collider, Journal of Instrumentation9 (2014) C02022
2014
-
[15]
Tytgat, C
M. Tytgat, C. Combaret, C. Devanne, G. Garillot, G. Grenier, I. Laktineh et al.,Towards the T-SDHCAL hadronic calorimeter for a future Higgs factory, Nucl. Instrum. Meth. A1077 (2025) 170520
2025
-
[16]
Frank, F
M. Frank, F. Gaede, M. Petric and A. Sailer,AIDASoft/DD4hep, Oct., 2018. 10.5281/zenodo.592244
2018 doi
-
[17]
Allison, K
J. Allison, K. Amako, J. Apostolakis, P. Arce, M. Asai, T. Aso et al.,Recent developments in Geant4, Nucl. Instrum. Meth. A835 (2016) 186
2016
-
[18]
Gaede,Marlin and LCCD: Software tools for the ILC,Nucl
F. Gaede,Marlin and LCCD: Software tools for the ILC,Nucl. Instrum. Meth. A559 (2006) 177
2006
-
[19]
Breton, A
D. Breton, A. Irles, J. Jeglot, J. Maalmi, R. Pöschl and D. Zerwas,CALICE SiW ECAL—Development and performance of a highly compact digital readout system,Journal of Instrumentation15(2020) C05074
2020
-
[20]
Z. Deng, Y. Li, Y. Wang, Q. Yue, Z. Yang, D. Boumediene et al.,Resistive Plate Chamber digitization in a hadronic shower environment,Journal of Instrumentation11(2016) P06014
2016
-
[22]
Han,Induced charge signal of a glass RPC detector,Chinese Physics C38(2014) 046002
R. Han,Induced charge signal of a glass RPC detector,Chinese Physics C38(2014) 046002
2014
-
[23]
Brun and F
R. Brun and F. Rademakers,ROOT: an object oriented data analysis framework,Nucl. Instrum. Meth. A 389 (1997) 81
1997
-
[24]
James and M
F. James and M. Roos,Minuit: a system for function minimization and analysis of the parameter errors and correlations,Comput. Phys. Commun.10 (1975) 343
1975
-
[25]
Pedregosa, G
F. Pedregosa, G. Varoquaux, A. Gramfort, V. Michel, B. Thirion, O. Grisel et al.,Scikit-learn: Machine Learning in Python,Journal of Machine Learning Research12(2011) 2825
2011
-
[26]
Green,Calorimetry at a Future Linear Collider, Ph.D
S. Green,Calorimetry at a Future Linear Collider, Ph.D. thesis, Cambridge U., 2017. – 19 –
2017
Reviewed August 2, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.