Pith. sign in

REVIEW 3 major objections 6 minor 33 references

A charge-only maximum-likelihood method using calibration-derived nPE maps reconstructs both radius and visible energy of positron events in a ton-scale spherical liquid-scintillator detector, achieving in Monte Carlo a vertex bias below 2

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · deepseek-v4-flash

2026-08-01 14:15 UTC pith:BJMI56RA

load-bearing objection Honest but overstated application of a known likelihood template method to a small LS detector; the missing energy-resolution number and radius-dependent vertex claim need fixing before the headline results can be trusted. the 3 major comments →

arxiv 2607.18797 v1 pith:BJMI56RA submitted 2026-07-21 physics.ins-det hep-ph

A method for energy and radius reconstruction with simulated charge information in a ton-scale liquid scintillator detector like Taishan Antineutrino Observatory

classification physics.ins-det hep-ph
keywords liquid scintillator detectorSiPMcharge reconstructionnPE mapmaximum likelihoodvertex reconstructionenergy reconstructionreactor antineutrino
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

The paper tries to establish that a compact spherical liquid-scintillator detector can reconstruct an event's radius and visible energy using only the charge collected by SiPMs, without timing information. The method builds expected 'number of photoelectrons per MeV' maps from simulated calibration sources, interpolates them to arbitrary radii, and maximizes a Poisson-Gaussian likelihood over radius and energy. In Monte Carlo, positron vertices are reconstructed with bias below 2 cm and resolution below 4 cm, with energy bias growing mainly from energy leakage near the detector edge. If the method holds on real data, it would give near-detector reactor experiments a simple, calibration-driven way to measure antineutrino spectra precisely while keeping instrumentation minimal.

Core claim

On its own terms, the paper claims that the expected light pattern on the SiPM array, expressed as an nPE map per unit visible energy and parameterized only by event radius and SiPM polar angle, is a sufficient statistic for simultaneous radius and energy reconstruction. Because the detector is spherical, maps measured along one axis describe all directions, and maps at five calibration radii can be linearly interpolated. The maximum-likelihood fit over radius and energy, with angles fixed by a center-of-charge seed, yields sub-4-cm radius resolution and sub-2-cm bias for positrons in simulation, and an energy resolution whose dominant term is photoelectron statistics, with quenching the lar

What carries the argument

The central object is the nPE map, mu_hat(r, theta_SiPM): the expected number of photoelectrons per MeV of visible energy seen by each SiPM as a function of the event radius r and the angle between the event position and the SiPM direction. It feeds a likelihood that combines a Poisson distribution for the number of photoelectrons on each channel with a Gaussian single-photoelectron charge response; maximizing the negative log-likelihood over r and E gives the reconstructed values. The map's separability into radius and angle plus linear interpolation between calibration points is what carries the argument.

Load-bearing premise

The method assumes that nPE maps measured at five calibration radii can be linearly interpolated to any radius, and that maps measured along one axis describe all directions because the detector is perfectly spherically symmetric.

What would settle it

Take a real calibration source at a radius not in the template grid, say 500 mm, and reconstruct it using templates built from the other radii; if the reconstructed radius bias or resolution exceeds the quoted few-centimeter values, the linear-interpolation or spherical-symmetry assumption breaks. Separately, compare real 137Cs/68Ge calibration maps with simulated maps: a mismatch in the reflection peak or angular shape would directly invalidate the MC-based performance numbers.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • Charge-only information is sufficient for radial and visible-energy reconstruction in compact spherical liquid-scintillator detectors, so timing readout is not needed for these quantities.
  • Calibration with gamma sources at a few positions (e.g., 137Cs and 68Ge) can generate templates that work across event energies; the energy-dependent reflection region changes the overall resolution by less than 1%.
  • The same template method transfers to other detectors of similar size and spherical geometry, provided their SiPM response can be calibrated.
  • The main limit to energy resolution is photoelectron statistics; crosstalk, dark noise, and reflections each contribute roughly half a percent or less, while scintillator quenching is the largest non-statistical effect.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • A natural test beyond the paper is to build nPE maps from real calibration data at the five radii and reconstruct events from an independent source position; success would confirm that linear interpolation and one-axis maps hold in the physical detector.
  • Because the quoted performance comes from the same simulated detector used to produce the templates, a mismatch between simulation and the real optical response (e.g., unmodeled reflections or SiPM non-uniformity) would change the numbers; the method's robustness to such mismatch is not yet demonstrated.
  • The radial resolution for positrons is worse than for electrons because annihilation gammas smear the light profile; combining charge-only reconstruction with timing information could push vertex resolution below the quoted 4 cm.
  • The energy bias rising with energy near the boundary points to energy leakage; an energy-dependent template interpolation or an explicit leakage correction could extend the method's useful fiducial volume.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 6 minor

Summary. The paper proposes a charge-only maximum-likelihood reconstruction method (QMLE) for simultaneous reconstruction of the radius and visible energy of IBD positrons in a ton-scale liquid scintillator detector of the JUNO-TAO type. The method builds templates of the expected number of photoelectrons per MeV as a function of event radius and SiPM angle, using simulated calibration sources (68Ge, 137Cs), interpolates these nPE maps between calibration radii and along axes, and then minimizes a Poisson-Gaussian charge likelihood over r and E, with θ and φ fixed from the CCA algorithm. The authors report Monte Carlo results for reconstructed radial bias and resolution for e+ and e− samples and study the effect of reflections, cross-talk, dark noise, and quenching on energy resolution. The conclusion states that vertex bias below 2 cm and vertex resolution below 4 cm can be achieved for positrons, and that the method is applicable to similar spherical neutrino detectors.

Significance. If the claimed performance were fully substantiated, the method would be a useful contribution to reconstruction in compact liquid scintillator detectors: it uses only charge information, exploits spherical symmetry, and is intended to be calibrated from radioactive source data rather than from an analytical optical model. The template-based charge likelihood is a reasonable and transparent statistical framework, and the paper includes MC samples covering multiple energies and radii, and some robustness checks on cross-talk, dark noise, and reflections. However, the paper as presented does not support its headline claims: no energy resolution number is reported despite the abstract promising 'unprecedented energy resolution', and the vertex performance is qualified only in a small central region, not over the fiducial volume as the conclusion implies. These issues are central and require additional analysis or a substantial reframing of the claims.

major comments (3)
  1. [Sec. V, Sec. IV A, Fig. 6, Table II] The conclusion states 'vertex bias less than 2 cm and vertex resolution less than 4 cm can be achieved for the positrons' without radial qualification. Fig. 6(b) shows this holds only for small radii; for radii ≳500 mm the radial resolution is above 40 mm at all kinetic energies, and the bias grows to ~30–40 mm near the 650 mm fiducial boundary. Table II gives center values of 15–16 mm bias and 38–38.5 mm resolution for e+, so even at the center the resolution is only 2 mm below the claimed 4 cm. Because IBD vertices are roughly uniformly distributed in the LS volume, the volume-averaged resolution will be substantially worse than 4 cm. The paper should report the claimed performance as a function of radius and specify the fiducial radius over which the 'less than 4 cm' statement is valid.
  2. [Abstract, Sec. IV B, Figs. 9–10] The abstract's central promise is 'unprecedented energy resolution', but the paper never quotes an energy resolution value. Sec. IV B reports only relative changes ('less than 1%', 'around 0.5%', 'around 0.16%') and energy bias in Figs. 9 and 10; there is no σ(E_rec)/E or equivalent resolution metric anywhere. Without this number, the energy-reconstruction claim is unsupported. The authors should add energy resolution as a function of radius and energy for both nPE map types, and if the energy resolution is not competitive with existing methods, the abstract should be revised accordingly.
  3. [Sec. III, Figs. 3–4, Sec. IV] The performance numbers are obtained in a fully self-consistent MC setup: both the nPE-map templates and the reconstructed event samples are generated with the same Geant4 simulation and the same optical model. This is a closure test, not a validation against real calibration data. The paper emphasizes that templates will be built from calibration sources, but it does not test sensitivity to template mismatch—e.g., interpolation errors between the five calibration radii, imperfect source-position knowledge, or differences between simulated and actual reflection/absorption. A concrete test would be to perturb the optical parameters or source positions when building templates and then reconstruct the nominal samples, showing that the quoted biases and resolutions are stable at the claimed level.
minor comments (6)
  1. [Throughout] The manuscript contains many typos and grammatical errors, e.g., 'scintilling', 'the the', 'achi-eved', 'reletive', 'infromation', 'reconstrucede +', 'met-hod', and 'eutrinos'. A thorough language edit is needed.
  2. [Figs. 9–10] The y-axis label in Figs. 9 and 10 reads 'Mean(E_vis − E_rec_vis)' but the panel title says 'QMLE Vertex Bias'. This is misleading; it should be labeled 'energy bias'.
  3. [Sec. IV B] The relative changes 'less than 1%', 'around 0.5%', and 'around 0.16%' are quoted without statistical uncertainties or the baseline resolution values to which they refer. Please provide the absolute resolutions and errors for these comparisons.
  4. [Sec. II / Fig. 1] The calibration system description would benefit from a clear statement of the radius coverage of the ACU and CLS calibration points, and how the five radii in Fig. 3 (0, 200, 375, 650, 850 mm) correspond to the actual calibration path.
  5. [Eq. (4)] The notation 'unhit' and 'hit' is not formally defined; please state that unhit channels are those with zero recorded photoelectrons. Also, the Gaussian charge response in Eq. (2) uses the same S_i for all k; its validity for small k should be commented on.
  6. [Sec. III] The interpolation from five calibration radii to arbitrary radii is stated to work 'similar' for cubic and polynomial methods, but no quantitative comparison is shown. It would help to report the interpolation residuals or to state that the results are insensitive to the interpolation order.

Circularity Check

0 steps flagged

No circularity: the template likelihood is self-contained; MC-only evaluation limits generality but does not make the derivation circular.

full rationale

The claimed derivation is not circular. The nPE templates (Sec. III, Fig. 3) are generated by simulating 68Ge/137Cs sources at fixed calibration positions; Eq. (3) then writes mu_i(r,E)=E·mu_hat_i(r,theta_SiPM), and Eq. (5) minimizes the likelihood over free parameters r and E. The target radius and energy are not used to build the templates, and the test samples (Table I: e+ at 0-70 cm with 0-4 MeV) are independent of the calibration points. Thus Eq. (5) does not reduce by construction to its inputs; there is no fitted parameter renamed as a prediction. The radial degradation visible in Fig. 6 is a performance caveat, not circularity: the quoted '<2 cm bias, <4 cm resolution' holds only near the center, and the MC-only evaluation means the numbers measure self-consistency with the optical model rather than robustness to a mismatched real detector. No load-bearing self-citation appears: Ref. [19] includes an author but is an unrelated E nu ES physics citation; the template/likelihood machinery cites external prior work ([26]) and uses standard Geant4 simulation. Therefore no circular step is identified.

Axiom & Free-Parameter Ledger

4 free parameters · 4 axioms · 0 invented entities

No new entities are invented. The free parameters are the choices inside the likelihood (SPES, dark-noise subtraction, interpolation) plus the analysis cut. The main input assumption is the Poisson-Gaussian charge model and the spherical-symmetry reduction of the nPE map.

free parameters (4)
  • nPE map interpolation order = linear
    The nPE map at arbitrary radii is obtained by linear interpolation between calibration points; cubic and polynomial were tried and "give similar results" (Sec. III). The choice affects the reconstruction but is not optimized or justified with a number.
  • theta_cut for reflection region = 160 deg
    When studying the effect of SiPM reflections, the reconstruction excludes SiPMs with theta > 160 deg (Sec. IV.B). This is an analysis choice that changes the energy resolution by <1%, but the final quoted method appears to keep the reflection region included.
  • Single photoelectron spectrum Qi, Si = not stated
    The likelihood requires the mean and sigma of the single-photoelectron spectrum (SPES) of each SiPM (Eq. 2). The paper does not state the values used in the simulation. These are detector parameters that enter the fit and affect the charge likelihood.
  • Dark noise rate = 20 Hz/mm^2
    Dark noise is "considered as 20Hz/mm^2 in simulation" and is subtracted when building nPE maps (Sec. III). This is a detector parameter from the simulation, not a free parameter fitted to data, but it is an input assumption affecting the templates.
axioms (4)
  • domain assumption The number of photoelectrons k on each SiPM follows a Poisson distribution with mean mu_i(r,E).
    This is the central statistical model in Eq. (1). It ignores any overdispersion from correlated noise or non-Poisson photon statistics; the paper does not test this assumption against simulation output.
  • domain assumption The single-photoelectron charge response for k photoelectrons is Gaussian with mean k*Qi and variance k*Si^2.
    Eq. (2) assumes Gaussian SPES and that the charge variance scales linearly with k. The paper does not justify this against the simulated SPES.
  • domain assumption The expected nPE map mu_hat(r, theta_SiPM) is independent of the event direction and azimuthal angle.
    The method relies on spherical symmetry to reduce the map to (r, theta_SiPM) and to use maps generated along a single axis (Fig. 4). The paper verifies this in simulation for 137Cs at 400 mm, but real detectors may have asymmetric support structures and SiPM gaps.
  • domain assumption The energy dependence of the nPE map is negligible over the fit range.
    The likelihood uses mu_hat(r, theta) per MeV and multiplies by E (Eq. 3). Fig. 5 shows the nPE map depends on energy in the reflection region. The paper argues the effect is <1% after excluding theta>160, but the nominal method (with reflections included) uses an energy-independent template per MeV.

pith-pipeline@v1.3.0-alltime-deepseek · 11707 in / 7524 out tokens · 53977 ms · 2026-08-01T14:15:37.897994+00:00 · methodology

0 comments
read the original abstract

Small neutrino detectors are a ton level detectors which can be placed very close to the core of Nuclear Power Plant. In some detectors liquid scintilling (LS) material is used as the detecting material. The antineutrinos from the reactor core fall on the liquid scintillator of the detector where they deposit energy via Inverse Beta Decay process (IBD). The energy absorbed by liquid scintillator is re-emitted in the form of scintillation. These photons then travel through the scintillating material and hit the Silicon Photo Multipliers (SiPMs) which are installed on the inner surface of detector's spherical copper shell. These SiPMs absorb the photons to give a charge output signal. The energy and radius reconstruction is done using the information of charge collected by the SiPMs. Due to factors like large photo-coverage with large photon detection efficiency, small spherical detector size and low temperature operation, small size LS detectors can achieve an unprecedented energy resolution. In this paper, we have used a template-dependent method exploiting simulated data to reconstruct the event radius and energy. This method uses response functions generated using radioactive source calibration data and the charge information to reconstruct the energy and the radius of the event by constructing maximizing a likelihood function. This methodology is applicable to all similar size spherical neutrino detector experiments.

Figures

Figures reproduced from arXiv: 2607.18797 by Randhir Singh, Xiaochuan Xie, Yichen Li.

Figure 1
Figure 1. Figure 1: FIG. 1: The detector setup, which consists of the central detector, the calibration system, consisting of an [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Event parameters in CD [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Simulated nPE maps generated by [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: nPE maps generated by [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: nPE maps generated at (400, 0, 0) mm by [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: Bias(left panel) and Resolution (right panel) of the reconstructed radius for [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: Bias(left panel) and Resolution (right panel) of the reconstructed radius for [PITH_FULL_IMAGE:figures/full_fig_p008_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8: Comparison of the nPE maps generated with and without considering cross-talk effect. [PITH_FULL_IMAGE:figures/full_fig_p009_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9: Energy Bias for the reconstructed [PITH_FULL_IMAGE:figures/full_fig_p011_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10: Energy Bias for the reconstructed [PITH_FULL_IMAGE:figures/full_fig_p011_10.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

33 extracted references · 3 canonical work pages

  1. [1]

    First results from KamLAND: Evidence for reactor anti-neutrino disappearance,

    K. Eguchiet al.[KamLAND], “First results from KamLAND: Evidence for reactor anti-neutrino disappearance,” Phys. Rev. Lett.90, 021802 (2003) doi:10.1103/PhysRevLett.90.021802 [arXiv:hep-ex/0212021 [hep-ex]]

  2. [2]

    Looking inside the Sun with the Borexino experiment: detection of solar eutrinos from the proton-proton chain and the CNO cycle,

    M. Redchuk, “Looking inside the Sun with the Borexino experiment: detection of solar eutrinos from the proton-proton chain and the CNO cycle,” doi:10.18154/R WTH-2020-12244

  3. [3]

    Indication of Reactor ¯ν e Disappearance in the Double Chooz xperiment,

    Y. Abeet al.[Double Chooz], “Indication of Reactor ¯ν e Disappearance in the Double Chooz xperiment,” Phys. Rev. Lett. 108, 131801 (2012) doi:10.1103/PhysRevLett.108.131801 [arXiv:1112.6353 [hep-ex]]

  4. [4]

    Observation of electron-antineutrino disappearance at Daya Bay,

    F. P. Anet al.[Daya Bay], “Observation of electron-antineutrino disappearance at Daya Bay,” Phys. Rev. Lett.108, 171803 (2012) doi:10.1103/PhysRevLett.108.171803 [arXiv:1203.1669 [hep-ex]]

  5. [5]

    Observation of Reactor Electron Antineutrino Disappearance in the RENO Experiment,

    J. K. Ahnet al.[RENO], “Observation of Reactor Electron Antineutrino Disappearance in the RENO Experiment,” Phys. Rev. Lett.108, 191802 (2012) doi:10.1103/PhysRevLett.108.191802 [arXiv:1204.0626 [hep-ex]]

  6. [6]

    JUNO physics and detector,

    A. Abuslemeet al.[JUNO], “JUNO physics and detector,” Prog. Part. Nucl. Phys.123, 103927 (2022) doi:10.1016/j.ppnp.2021.103927 [arXiv:2104.02565 [hep-ex]]

  7. [7]

    First measurement of reactor neutrino oscillations at JUNO,

    A. Abuslemeet al.[JUNO], “First measurement of reactor neutrino oscillations at JUNO,” [arXiv:2511.14593 [hep-ex]]

  8. [8]

    G. h. Huang, W. Jiang, L. j. Wen, Y. f. Wang and W. M. Luo, Nucl. Sci. Tech.34(2023) no.6, 83 doi:10.1007/s41365- 023-01240-0 [arXiv:2211.16768 [physics.ins-det]]

  9. [9]

    Prediction of Energy Resolution in the JUNO Experiment,

    A. Abuslemeet al.[JUNO], “Prediction of Energy Resolution in the JUNO Experiment,” Chin. Phys. C49(2025) no.1, 013003 doi:10.1088/1674-1137/ad83aa [arXiv:2405.17860 [hep-ex]]

  10. [10]

    Constraints on non-standard neutrino interactions from Borexino extended data-set,

    V. Antonelliet al.[BOREXINO], “Constraints on non-standard neutrino interactions from Borexino extended data-set,” [arXiv:2602.08685 [hep-ex]]

  11. [11]

    Determination of the Neutrino Mass Hierarchy at an Intermediate Baseline,

    L. Zhan, Y. Wang, J. Cao and L. Wen, “Determination of the Neutrino Mass Hierarchy at an Intermediate Baseline,” Phys. Rev. D78, 111103 (2008) doi:10.1103/PhysRevD.78.111103 [arXiv:0807.3203 [hep-ex]]

  12. [12]

    Neutrinoless Double-Beta Decay: Status and Prospects,

    M. J. Dolinski, A. W. P. Poon and W. Rodejohann, “Neutrinoless Double-Beta Decay: Status and Prospects,” Ann. Rev. Nucl. Part. Sci.69, 219-251 (2019) doi:10.1146/annurev-nucl-101918-023407 [arXiv:1902.04097 [nucl-ex]]

  13. [13]

    Light sterile neutrinos,

    S. Gariazzo, C. Giunti, M. Laveder, Y. F. Li and E. M. Zavanin, “Light sterile neutrinos,” J. Phys. G43, 033001 (2016) doi:10.1088/0954-3899/43/3/033001 [arXiv:1507.08204 [hep-ph]]

  14. [14]

    Future searches for light sterile neutrinos at nuclear reactors,

    J. M. Berryman, L. A. Delgadillo and P. Huber, “Future searches for light sterile neutrinos at nuclear reactors,” Phys. Rev. D105, no.3, 035002 (2022) doi:10.1103/PhysRevD.105.035002 [arXiv:2104.00005 [hep-ph]]

  15. [15]

    JUNO sensitivity to the annihilation of MeV dark matter in the galactic halo,

    A. Abuslemeet al.[JUNO], “JUNO sensitivity to the annihilation of MeV dark matter in the galactic halo,” JCAP09 (2023), 001 doi:10.1088/1475-7516/2023/09/001 [arXiv:2306.09567 [hep-ex]]

  16. [16]

    q i”generated by thei th SiPM after being hit by “k

    and also has some potential through elastic neutrino-electron scattering (EνES) [19]. The reconstruction of the individual events is one of the most important task in order to obtain neutrino information such as their energies. In general, the event reconstruction should provide us the information about the energy, vertex and the interaction time (relativ...

  17. [17]

    TAO Conceptual Design Report: A Precision Measurement of the Reactor Antineutrino Spectrum with Sub-percent Energy Resolution,

    A. Abuslemeet al.[JUNO], “TAO Conceptual Design Report: A Precision Measurement of the Reactor Antineutrino Spectrum with Sub-percent Energy Resolution,” [arXiv:2005.08745 ]

  18. [18]

    On the determination of anti-neutrino spectra from nuclear reactors,

    P. Huber, “On the determination of anti-neutrino spectra from nuclear reactors,” Phys. Rev. C84, 024617 (2011) [erratum: Phys. Rev. C85, 029901 (2012)] doi:10.1103/PhysRevC.85.029901 [arXiv:1106.0687 [hep-ph]]

  19. [19]

    Improved Predictions of Reactor Antineutrino Spectra,

    T. A. Mueller,et al.“Improved Predictions of Reactor Antineutrino Spectra,” Phys. Rev. C83, 054615 (2011) doi:10.1103/PhysRevC.83.054615 [arXiv:1101.2663 [hep-ex]]

  20. [20]

    Elastic neutrino-electron scattering perspectives at nuclear reactors,

    L. A. Delgadillo, Q. Liu and R. Singh, “Elastic neutrino-electron scattering perspectives at nuclear reactors,” Phys. Rev. D113, no.11, 113007 (2026) doi:10.1103/m6yl-nzd9 [arXiv:2602.21523 [hep-ph]]

  21. [21]

    Maximum likelihood reconstruction of a detector with reflective panels,

    L. J. Wen, J. Cao, Y. Wang, Z. Wang and C. Yang, “Maximum likelihood reconstruction of a detector with reflective panels,” Nucl. Instrum. Meth. A629, 296-302 (2011) doi:10.1016/j.nima.2010.11.030

  22. [22]

    Vertex and energy reconstruction in JUNO with machine learning methods,

    Z. Qian, V. Belavin, V. Bokov, R. Brugnera, A. Compagnucci, A. Gavrikov, A. Garfagnini, M. Gonchar, L. Khatbullina and Z. Li,et al.“Vertex and energy reconstruction in JUNO with machine learning methods,” Nucl. Instrum. Meth. A 1010(2021), 165527 doi:10.1016/j.nima.2021.165527 [arXiv:2101.04839 [physics.ins-det]]

  23. [23]

    Vertex reconstruction in the TAO experiment,

    H. Shi, J. Wang, G. Cao, W. Wang and Y. Wei, “Vertex reconstruction in the TAO experiment,” [arXiv:2508.06293 [physics.ins-det]]

  24. [24]

    First-principle event reconstruction by time-charge readouts for TAO,

    X. Liu, W. Dou, B. Xu, H. Wang and G. Cao, “First-principle event reconstruction by time-charge readouts for TAO,” Eur. Phys. J. C85(2025) no.4, 438 [erratum: Eur. Phys. J. C85(2025) no.6, 653] doi:10.1140/epjc/s10052-025-14161-6 [arXiv:2403.01239 [hep-ex]]

  25. [25]

    Measurement of electron antineutrino oscillation based on 1230 days of operation of the Daya Bay experiment,

    F. P. Anet al.[Daya Bay], “Measurement of electron antineutrino oscillation based on 1230 days of operation of the Daya Bay experiment,” Phys. Rev. D95, no.7, 072006 (2017) doi:10.1103/PhysRevD.95.072006 [arXiv:1610.04802 [hep-ex]]

  26. [26]

    Calibration strategy of the JUNO-TAO experiment,

    H. Xu,et al.“Calibration strategy of the JUNO-TAO experiment,” Eur. Phys. J. C82, no.12, 1112 (2022) doi:10.1140/epjc/s10052-022-11069-3 [arXiv:2204.03256 ]

  27. [27]

    A new method of energy reconstruction for large spherical liquid scintillator detectors,

    W. Wu, M. He, X. Zhou and H. Qiao, “A new method of energy reconstruction for large spherical liquid scintillator detectors,” JINST14, no.03, P03009 (2019) doi:10.1088/1748-0221/14/03/P03009 [arXiv:1812.01799 ]. 13

  28. [28]

    A side-by-side comparison of Daya Bay antineutrino detectors,

    F. P. Anet al.[Daya Bay], “A side-by-side comparison of Daya Bay antineutrino detectors,” Nucl. Instrum. Meth. A685, 78-97 (2012) doi:10.1016/j.nima.2012.05.030 [arXiv:1202.6181 ]

  29. [29]

    GEANT4 - A Simulation Toolkit,

    S. Agostinelliet al.[GEANT4], “GEANT4 - A Simulation Toolkit,” Nucl. Instrum. Meth. A506, 250-303 (2003) doi:10.1016/S0168-9002(03)01368-8

  30. [30]

    Optimization of the JUNO liquid scintillator composition using a Daya Bay antineutrino detector,

    A. Abuslemeet al.[JUNO and Daya Bay], “Optimization of the JUNO liquid scintillator composition using a Daya Bay antineutrino detector,” Nucl. Instrum. Meth. A988, 164823 (2021) doi:10.1016/j.nima.2020.164823 [arXiv:2007.00314 ]

  31. [31]

    A complete optical model for liquid-scintillator detectors,

    Y. Zhang, Z. Y. Yu, X. Y. Li, Z. Y. Deng and L. J. Wen, “A complete optical model for liquid-scintillator detectors,” Nucl. Instrum. Meth. A967, 163860 (2020) doi:10.1016/j.nima.2020.163860 [arXiv:2003.12212 ]

  32. [32]

    Scintillations from Organic Crystals: Specific Fluorescence and Relative Response to Different Radiations,

    J. B. Birks, “Scintillations from Organic Crystals: Specific Fluorescence and Relative Response to Different Radiations,” Proc. Phys. Soc. A64, 874-877 (1951) doi:10.1088/0370-1298/64/10/303

  33. [33]

    The Nature of the Saturation Effect of Fluorescent Scintillators,

    C. N. Chou, “The Nature of the Saturation Effect of Fluorescent Scintillators,” Phys. Rev.87, no.5, 904-905 (1952) doi:10.1103/PhysRev.87.904