{"id":"35f84f3d-fe88-4864-8417-16d06ff4421b","arxiv_id":"2502.05050","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":27,"one_line_summary":"A CeBr3 plus plastic scintillator telescope was calibrated with proton, lithium, and carbon beams, giving Birks' law parameters, energy resolution of about 10 MeV, and time resolution between 171 and 282 ps.","lead":"This paper calibrates a detector telescope that measures the energy loss and remaining energy of particles produced by cancer treatment beams, using a plastic scintillator and a CeBr3 crystal. It reports that the telescope follows an established quenching law and can tell particle types apart, a step toward better measurement of secondary radiation in hadrontherapy.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"CeBr3 calibration fits have reduced chi-square ~10^3–10^4, so the claimed Birks' law adherence and parameter precision are not supported by the data.","rationale":"I read the paper as a detector calibration study whose central claim is that the Delta-E-E telescope is successfully calibrated, with the CeBr3 response obeying Birks' law up to 2350 MeV and an energy resolution of about 10 MeV. For this claim to hold, the Birks' law fits must describe the data. The reported reduced chi-square values show they do not. For example, in Figure 5(a), the carbon fit has chi^2/ndf about 987; the proton fit about 2593. These are not borderline cases; they signal systematic deviations from the model. The authors only propagate statistical uncertainties on the fitted parameters, producing extremely small error bars that falsely imply high precision. The physical deviations may stem from unmodeled quenching effects, PMT nonlinearity, or incorrect energy-loss values, but whatever the cause, the current analysis does not establish the quantitative calibration. The reader's weakest assumption about the computed energy scale is valid and important, but the poor fit quality is a more direct, internal inconsistency that independently invalidates the quoted parameters and resolutions. A residual analysis or model-extension test would settle whether the misfit is due to underestimated errors or a wrong functional form. If the misfit persists, the quantitative claims should be scaled back; if the misfit disappears with a systematic error term, then the paper should still report that systematic term and revise the parameter uncertainties. I therefore keep the CONDITIONAL verdict because the qualitative Delta-E-E separation and timing performance are still valuable, but the calibration constants and resolution values must be revised with a proper treatment of systematic errors and model verification.","tokens_in":15918,"tokens_out":5334,"duration_ms":52856,"concrete_test":"Using the published data points from Figure 5, re-fit the CeBr3 amplitude versus E_dep with an additional systematic error term added in quadrature to the statistical errors, and find the minimum systematic error needed to achieve chi^2/ndf near 1. If the required systematic error exceeds 10% of the amplitude, or if no systematic error can make the fit acceptable, the quoted parameter precision and the energy-resolution claim are invalid. Alternatively, test for model misspecification by adding a quadratic term to the Birks' denominator and checking whether it is significantly non-zero; a significant term would falsify the simple Birks' law parametrization.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim that the CeBr3 response follows Birks' law up to 2350 MeV (Section 3.2, Figure 5) is undermined by the fit quality. For the +350 V data, chi^2/ndf is about 987 for carbons and 2593 for protons; for +400 V it is even larger (carbons: about 1.36e4). These values are orders of magnitude above unity, indicating large systematic deviations from the three-parameter Birks' law. Despite this, Table 2 quotes parameter uncertainties at the 10^-4 level (e.g., S' = 0.1495 +/- 0.0005), which are statistical-only and meaningless when the model does not describe the data. Consequently, the extracted S', k'B, A0 values and the derived energy resolution of about 10 MeV are not reliable. The 'adherence to Birks' law' statement is therefore not demonstrated. This is more immediately fatal to the calibration claim than the unverified energy scale, since even with perfectly known deposited energies, the model fails. The two issues compound: an incorrect energy scale could also partially absorb the misfit.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports the calibration and characterization of a ΔE-E telescope consisting of a thin plastic scintillator (ΔE) and a CeBr3 crystal (E) for detecting secondary charged particles in hadron therapy. Data were taken with proton, lithium, and carbon beams at multiple facilities (Cyrcé, CAL, GSI, CNAO, GANIL). The response of both scintillators is fitted with an adapted Birks' law, and the authors claim adherence up to 50 MeV deposited energy for the plastic and up to 2350 MeV for CeBr3. They also report energy resolution on the order of 10 MeV and time resolutions of 282±1 ps for 180 MeV/u 12C and 171±1 ps for 25 MeV protons. The manuscript concludes that the calibrated telescope can identify secondary fragments by charge and will provide input for Monte Carlo benchmarking.","tokens_in":16323,"tokens_out":3962,"duration_ms":40255,"significance":"If the calibration and performance claims were validated, the detector would be a useful tool for measuring secondary charged particles produced by therapeutic ion beams, addressing a recognized gap in experimental nuclear data for hadron therapy. The paper presents data from several clinical and accelerator facilities, and the combination of a CeBr3 crystal with a plastic scintillator in a ΔE-E telescope is relatively novel. The time resolution values are competitive with existing systems. However, the central quantitative claims are currently not supported by the reported fit quality: the chi-square values of the Birks' law fits and the resolution fits are orders of magnitude larger than unity, and the quoted parameter uncertainties are unrealistically small. The energy scale used in all fits is itself computed, not independently measured, further weakening the reliability of the extracted parameters.","major_comments":[{"comment":"The CeBr3 calibration fits have reduced chi-square values of order 10^3–10^4 (for example, χ²/ndf = 1.086e4/11 for carbon at +350 V, χ²/ndf = 1.037e4/4 for protons at +350 V, and even larger values for the +400 V fits). These values indicate large systematic deviations from the three-parameter Birks' law, yet the paper quotes parameter uncertainties at the 10^-4 level (e.g., S' = 0.1495±0.0005). The claim that 'adherence to Birks' law is observed up to 2350 MeV' is therefore not supported by the reported statistics. The authors should either report fits that are statistically acceptable, including systematic uncertainties and a discussion of the deviations, or restrict the claim to the energy range where the model actually describes the data.","section":"§3.2, Fig. 5, Table 2"},{"comment":"The energy resolution fits also have extremely large chi-square per degree of freedom (e.g., χ²/ndf = 2421/10 for the plastic in Fig. 4(b), and χ²/ndf = 1.201e4/10 and 9160/31 for the CeBr3 in Fig. 6). The extracted parameters a and b therefore do not meaningfully represent the data, and the statement that the energy resolution is 'on the order of 10 MeV' is not supported. The manuscript should present the resolution data with a model that adequately fits the points and with uncertainties that reflect the scatter.","section":"§3.1 and §3.2, Figs. 4(b) and 6"},{"comment":"The deposited energies used as the x-axis in all Birks' law fits are not measured independently but are computed from nominal beam energies, degrader thicknesses, and stopping-power calculations. Any systematic error in these assumptions shifts all fitted constants (S', k'B, A0) and the extracted resolutions. The Geant4 comparison in Fig. 7 is a partial check but uses the same class of stopping-power calculations, so it does not validate the absolute energy scale. This limitation should be stated explicitly, and an independent validation (e.g., via time-of-flight or a range measurement) would strengthen the calibration considerably.","section":"§2.3, Eq. (2.2), Table 1"},{"comment":"The intrinsic energy resolution is obtained by subtracting in quadrature the Geant4 standard deviation from the measured one: σ_E = sqrt(σ_det^2 - σ_G4^2). This is valid only if the two contributions are independent and Gaussian, and if the Geant4 simulation contains exactly the energy straggling and beam-scattering contributions and nothing else. The manuscript itself shows in Fig. 7(a) that the Geant4 distribution does not reproduce the fragmentation tail, so the full distributions are not identical. The values of σ_det and σ_G4 are not reported, making it impossible to assess whether the subtraction is physical. Please show these values separately and discuss the validity of the subtraction, including any cases where σ_G4 could approach or exceed σ_det.","section":"§2.4, Eq. (2.3)"}],"minor_comments":[{"comment":"The abstract contains a spacing error: 'dosedeposition' should be 'dose deposition'.","section":"Abstract"},{"comment":"There is a typo: 'fo X-ray and gamma-ray' should be 'for X-ray and gamma-ray'.","section":"§2.1"},{"comment":"The axis label in Fig. 9 reads 'Time (ms)' while the values and text are in nanoseconds; correct the units.","section":"§3.4, Fig. 9"},{"comment":"Table 2 lists only S' and k'B but the fit model Eq. (2.2) includes a third parameter A0. Either include A0 in the table or state clearly that it was found compatible with zero for all configurations.","section":"Table 2"},{"comment":"The fit legends in Figures 3–6 use generic labels p0, p1, p2 alongside the physics parameters S', A0, kB'. Please make the notation consistent throughout.","section":"Figures 3–6"},{"comment":"The conclusion contains the typo 'the the CeBr3 crystal'; please correct it.","section":"Conclusion"},{"comment":"Reference [26] is garbled ('Cancer, Rev7.1(Realese 11.1), April July 31st, 2023'); please correct the citation to the Geant4 physics list documentation.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The manuscript addresses a relevant experimental need in hadron therapy, and the dataset from multiple facilities is valuable. However, the reported chi-square values for the central calibration and resolution fits are so large that the quantitative claims (Birks' law parameters, energy resolution) are not credible as presented. The issues are fixable in a revision: the authors need to reconsider the fitting procedure, report goodness-of-fit and systematic uncertainties, discuss the deviations, and validate or at least caveat the computed energy scale. If the fits cannot be brought to a statistically acceptable level, the claims should be restricted accordingly. I therefore recommend major revision rather than rejection, because the underlying data and the detector concept have merit, but the current analysis does not support the conclusions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is a workmanlike calibration study of a CeBr3+plastic ΔE-E telescope for hadrontherapy secondary fragments. The new bits are real: first use of this detector combination for that purpose, new measured Birks-law constants for CeBr3 under protons, lithium, and carbon, and a working Z separation in ΔE-E space. Timing resolution (171–282 ps) is competitive with similar setups. The authors also deserve credit for being explicit that the fitted constants are effective calibration parameters for their full detection chain, not intrinsic material properties, and for using Geant4 simulations that do not incorporate the fitted parameters, so there is no circularity in that comparison.\n\nThe central problem is that the quantitative claims do not survive contact with the paper's own fit statistics. The CeBr3 calibration curves in Figure 5 have reduced chi-square values around 1000 for +350 V data and even larger at +400 V (e.g., 1.086e4/11 for carbons). The resolution fits in Figure 6 are similarly catastrophic (chi2/ndf ~ 10^3–10^4). These are not mild discrepancies; they mean the three-parameter Birks form is not actually describing the data to within any reasonable statistical tolerance. Yet the paper states 'adherence to Birks’ law' as a headline result and quotes parameter uncertainties at the 10^-4 level. Those uncertainties are statistical-only and meaningless when the model is rejected by the data. At minimum, the authors need to report the chi2 values in the text, discuss what causes the deviations, and likely add systematic terms or a different functional form.\n\nThe second soft spot is the energy scale. The x-axis of every fit is computed from nominal beam energies, degrader thicknesses, and stopping-power calculations, not measured in an independent way. A systematic error in that scale shifts every fitted constant and resolution. The Geant4 comparison in Figure 7 is a partial check, but it uses the same class of stopping-power calculations, and the proton comparison shows a ~6% peak shift that is waved away. Systematic uncertainties are never propagated.\n\nWhat is solid: the qualitative behavior—monotonic response, charge-dependent bands, plausible timing—is probably right, and the paper is honest about the limits of its effective-parameter interpretation. But as a calibration paper, its quantitative output is not usable until the fits are understood and the systematics are handled. No data or code are released, which makes independent checking harder.\n\nMy recommendation: this deserves a serious referee, but only with major revision. The authors should report all reduced chi2 values, re-examine the Birks form or add systematic terms, propagate energy-scale uncertainties, and release the calibration data. If those changes are made, the paper would be a useful reference for the hadrontherapy instrumentation community. As it stands, I would not cite the fitted constants in my own work.","headline":"Useful engineering calibration, but the CeBr3 Birks-law fit is quantitatively unsupported by the paper's own chi2 values.","tokens_in":17003,"tokens_out":1616,"would_cite":false,"duration_ms":19314,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A plastic-plus-CeBr3 telescope follows Birks' law for proton, carbon, and lithium beams in hadron therapy, delivering ~10 MeV energy resolution, 171–282 ps timing, and charge separation of secondary fragments.","keywords":["hadrontherapy instrumentation","dE/dx detectors","heavy-ion detectors","CeBr3 scintillator","Birks' law","energy calibration","particle identification","time resolution"],"falsifier":"Direct a beam of a species and energy not used in the calibration, such as oxygen-16 ($^{16}$O) near 150 MeV/u or $\\alpha$ particles through the same PMMA degraders, at the telescope and check whether the pulse amplitudes fall on the extrapolated Birks' law curves within uncertainties. A second check is to reconstruct fragment energies independently by time of flight over a known flight path (or with a magnetic spectrometer) for the same events and compare with the energies the calibration assigns.","tokens_in":15643,"feed_emoji":"⚛️","tokens_out":14889,"duration_ms":130909,"temperature":0.7,"pith_summary":"This paper calibrates a $\\Delta E$–$E$ telescope—a thin plastic scintillator in front of a CeBr$_3$ crystal—for the ion beams used in hadron therapy, and argues that the calibrated device can measure and identify the secondary charged fragments those beams produce in tissue-equivalent material. The central result is that both scintillators respond according to Birks' law across the clinical range: up to 50 MeV of deposited energy for the plastic and up to 2350 MeV for the CeBr$_3$ crystal, with energy resolution on the order of 10 MeV and coincidence time resolution of 282±1 ps for 180 MeV/u carbon ions (171±1 ps for 25 MeV protons). If the calibration is right, the telescope supplies the experimental fragment-production data that treatment-planning Monte Carlo codes currently lack, which is what justifies shrinking the safety margins those codes now force. The paper is careful to treat the fitted Birks' parameters as effective calibration constants for the full detection chain, not intrinsic properties of the scintillator materials.","feed_headline":"Telescope catches therapy-beam fragments by charge","feed_subtitle":"Plastic-plus-CeBr3 pair follows Birks law to 2350 MeV with ~10 MeV resolution and 171–282 ps timing.","key_machinery":"The load-bearing object is the telescope itself: a 6×6×0.2 cm³ plastic scintillator (EJ-228) that records the energy loss $\\Delta E$, followed by a 2×2 inch CeBr$_3$ crystal that records the remaining energy $E$, both read out by photomultipliers and digitized at 3.2 GHz. The argument is carried by the adapted Birks' law $A = (S' E + A_0)/(1 + k'_B E)$, used as a phenomenological three-parameter fit of pulse amplitude to deposited energy, with $S'$ a gain factor, $k'_B$ an effective quenching constant, and $A_0$ a pedestal; the paper stresses that these are effective parameters of the whole detection chain, not intrinsic material constants. The x-axis energies are not measured directly but computed from nominal beam energies, degrader thicknesses, and stopping-power tables, and Geant4 with the INCL++ fragmentation model supplies the straggling and scattering widths subtracted in quadrature, $\\sigma_E = \\sqrt{\\sigma_{det}^2 - \\sigma_{G4}^2}$, before resolutions are quoted; the same simulation is used to validate the calibration against the measured CeBr$_3$ energy spectra.","core_discovery":"The paper's claim is that one two-layer scintillator telescope can be calibrated with protons, lithium, and carbon beams from five accelerator facilities, and that the calibration carries over to its intended task: resolving secondary charged fragments by charge in $\\Delta E$–$E$ space. For the plastic scintillator the pulse amplitude as a function of computed deposited energy follows the adapted Birks' law $A = (S' E + A_0)/(1 + k'_B E)$ up to 50 MeV, for both photomultiplier and voltage configurations tested; for the CeBr$_3$ crystal the same functional form holds up to 2350 MeV of deposited energy, with fitted quenching constants near $10^{-3}$ MeV$^{-1}$ that stay consistent across ion species at a given voltage. Energy resolution, after subtracting Monte Carlo straggling in quadrature, is on the order of 10 MeV for both detectors, and the coincidence time resolution between the two scintillators is 282±1 ps for 180 MeV/u $^{12}$C and 171±1 ps for 25 MeV protons. The authors present this as the first optimization of a CeBr$_3$-plus-plastic $\\Delta E$–$E$ telescope for hadrontherapy secondaries, and show a $\\Delta E$ versus $E$ scatter plot, from a 200 MeV/u carbon beam on a tissue-equivalent target, in which fragments with $Z = 1$ through $Z = 6$ form distinct branches.","pith_inferences":["The near-identical $k'_B$ values for protons, lithium, and carbon at a given voltage (for example 0.000857 and 0.000864 MeV$^{-1}$ at +350 V) suggest the quenching is set by the detector chain rather than the ion; if that holds, a future calibration could fix one shared $k'_B$ and fit only $S'$ per species, cutting the beam time needed to calibrate at each new facility.","The fragmentation-tail discrepancy between data and the INCL++ simulation implies the calibrated telescope could quantify underprediction of light-fragment yields; the paper notes the discrepancy but does not turn it into a measured cross section.","A testable extension the authors do not pursue is to use the CeBr$_3$ crystal's gamma-ray and neutron sensitivity in coincidence with charged fragments, connecting charged and neutral fragmentation channels in the same run."],"forward_implications":["The telescope can be deployed in clinical-like conditions to record yields, charges, and energies of secondary charged fragments ($Z = 1$ to $Z = 6$) from proton and carbon beams on tissue-equivalent targets, as the distinct $\\Delta E$–$E$ branches illustrate.","The calibrated CeBr$_3$ energy spectra become experimental benchmarks for Monte Carlo treatment-planning codes; the paper's own comparison already exposes a fragmentation tail present in the data but largely absent from the simulation.","With 171–282 ps timing, the same telescope extends to $\\Delta E$–ToF particle identification for high-energy ions, tightening fragment identification beyond charge alone.","Because Birks' law holds up to 2350 MeV of deposited energy, one calibration curve per ion species and voltage setting covers the full clinical energy range, so no per-energy recalibration is needed before each measurement campaign."],"supporting_citations":[{"why":"Supplies the original Birks' law equation on which the adapted amplitude-versus-energy fit is built.","marker":"[21]"},{"why":"Provides the expected ranges for the quenching constant and conversion factor used to justify the fit model.","marker":"[22]"},{"why":"The Geant4 toolkit used to compute deposited energies, straggling, and scattering, and to validate the calibration against measured spectra.","marker":"[6]"},{"why":"The INCL++ physics list that models nuclear fragmentation for the Monte Carlo comparisons and straggling subtraction.","marker":"[26]"},{"why":"Motivates the ion-dependent calibration functions and reduced PMT voltages chosen to avoid saturation with heavy ions.","marker":"[11]"},{"why":"Supplies the quadrature-subtraction approach separating detector resolution from beam straggling, plus CeBr3 timing context.","marker":"[27]"},{"why":"Provides the ~300 ps time-resolution benchmark from a similar high-energy ion setup against which the telescope's timing is compared.","marker":"[31]"}],"fun_headline_variants":["CeBr3-plastic ΔE-E telescope calibrated for hadrontherapy fragments","First calibration of a two-scintillator telescope for therapy secondaries","Telescope resolves fragment charge in carbon beam on tissue target","Birks law fits CeBr3 and plastic detectors up to 2350 MeV","Energy resolution ~10 MeV, timing 171-282 ps in new telescope"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calibration assumes the energy each ion deposits in the two scintillators is already known from the nominal beam energy, the degrader thickness, and stopping-power calculations rather than measured independently, so any error in those computed energies shifts every fitted constant and quoted resolution.","fun_headline_variants_meta":{"raw":{"variants":["CeBr3-plastic ΔE-E telescope calibrated for hadrontherapy fragments","First calibration of a two-scintillator telescope for therapy secondaries","Telescope resolves fragment charge in carbon beam on tissue target","Birks law fits CeBr3 and plastic detectors up to 2350 MeV","Energy resolution ~10 MeV, timing 171-282 ps in new telescope"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001014,"raw_usage":{"total_tokens":4383,"prompt_tokens":1148,"completion_tokens":3235,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":764,"completion_tokens_details":{"reasoning_tokens":3147}},"tokens_in":764,"tokens_out":3235,"duration_ms":27005,"temperature":1.0,"reasoning_tokens":3147,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T20:25:37.574115+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Direct a beam of a species and energy not used in the calibration, such as oxygen-16 ($^{16}$O) near 150 MeV/u or $\\alpha$ particles through the same PMMA degraders, at the telescope and check whether the pulse amplitudes fall on the extrapolated Birks' law curves within uncertainties. A second check is to reconstruct fragment energies independently by time of flight over a known flight path (or with a magnetic spectrometer) for the same events and compare with the energies the calibration assigns.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the original Birks' law equation on which the adapted amplitude-versus-energy fit is built."},{"cited_title":"Knoll.Radiation Detection and Measurement","cited_arxiv_id":null,"evidence_quote":"Provides the expected ranges for the quenching constant and conversion factor used to justify the fit model."},{"cited_title":"Agostinelli et al","cited_arxiv_id":null,"evidence_quote":"The Geant4 toolkit used to compute deposited energies, straggling, and scattering, and to validate the calibration against measured spectra."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The INCL++ physics list that models nuclear fragmentation for the Monte Carlo comparisons and straggling subtraction."},{"cited_title":"Investigation of the LaBr3 scintillator response to heavy ions.Radiation Measurements, 115:43–48, August 2018","cited_arxiv_id":null,"evidence_quote":"Motivates the ion-dependent calibration functions and reduced PMT voltages chosen to avoid saturation with heavy ions."},{"cited_title":"Salvador et al","cited_arxiv_id":null,"evidence_quote":"Supplies the quadrature-subtraction approach separating detector resolution from beam straggling, plus CeBr3 timing context."},{"cited_title":"Salvador et al","cited_arxiv_id":null,"evidence_quote":"Provides the ~300 ps time-resolution benchmark from a similar high-energy ion setup against which the telescope's timing is compared."}],"review_version":1}