{"id":"7dd74342-5098-4199-9369-3c52dba5f144","arxiv_id":"1908.04447","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The center-of-gravity position estimator's discretization error is exactly expressed in terms of the signal's Fourier transform, and a triangular crosstalk response defines an ideal detector with zero error.","lead":"This paper derives exact formulas for the systematic error in center-of-gravity position measurements from segmented detectors, and shows that triangular crosstalk can eliminate that error for any signal shape. It matters because it turns a previously empirical calibration problem into a calculable one.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Ideal-detector definition omits the uniform-crosstalk condition; the zero-first-moment condition alone does not imply xg=epsilon for extended signals.","rationale":"The central Eq. (6) and its Poisson-summation derivation are correct under the stated infinite periodic, shift-invariant assumptions; I verified the algebra including the k=0 term that produces epsilon and the k != 0 terms that produce the Fourier series. The ideal-detector theorem, as proved in the body, is also correct because the proof uses both sum_n p=1 and sum_n (x-n tau)p=0. The problem is that the formal definition of an ideal detector, and the abstract's summary, state only the second condition. This is not a stylistic omission: without M0(x)=sum_n p(x-n tau) being constant, the same derivation yields xg(epsilon) as the center of gravity of M0(y)phi(y-epsilon), which generally differs from epsilon for extended signals. The explicit counterexample above demonstrates the failure. The fix is straightforward: define an ideal detector as a uniform-crosstalk function (sum_n p=1) satisfying sum_n (x-n tau)p=0, and state that qualification in the abstract and in the definition line. Since the underlying mathematics is sound and the correction is textual, I would not reject the paper, but acceptance should be conditional on making the uniformity condition part of the definition. This concern is distinct from the reader's infinite-array caveat, which I also agree with; hence partial agreement.","tokens_in":25291,"tokens_out":32466,"duration_ms":324908,"concrete_test":"Set tau=1 and define p(x)=(1-x)(1-0.5 cos(2 pi x)) for 0<=x<=1, p(-x)=p(x), zero for |x|>1. Then sum_n (x-n)p(x-n)=0 but sum_n p(x-n)=1-0.5 cos(2 pi x), which is not constant. Take phi(x)=2 on [-1/4,1/4], epsilon=1/4, so that phi is symmetric and supported away from zero. Evaluating xg(epsilon)=integral_0^{1/2} y(1-0.5 cos(2 pi y)) dy / integral_0^{1/2} (1-0.5 cos(2 pi y)) dy gives approximately 0.3006, not 0.25. This settles that the bare moment condition does not imply the 'any signal distribution' claim; adding sum_n p(x-n tau)=1 removes the counterexample.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"Section 6.6 proves xg(epsilon)=epsilon for crosstalk functions that are both uniform (sum_n p(x-n tau, tau1)=1) and have vanishing first lattice moment (sum_n (x-n tau)p(x-n tau, tau1)=0). The derivation explicitly uses the uniformity to collapse the term involving sum_n p into the identity, and the triangular example satisfies both conditions. However, the 2019-added standalone definition of an ideal detector and the abstract state only the vanishing-moment condition, and the reader's strongest claim repeats it without uniformity. Read literally, the condition is insufficient: for a symmetric p supported on [-tau,tau], one can choose p(x)=(tau-x)r(x) on [0,tau], p(-x)=p(x), with r symmetric about tau/2 but nonconstant. Then sum_n (x-n tau)p(x-n tau)=0 identically, while M0(x)=sum_n p(x-n tau)=r(x) is not constant. Repeating the paper's interchange gives xg(epsilon)=integral y M0(y) phi(y-epsilon) dy / integral M0(y) phi(y-epsilon) dy, which is not epsilon for an extended symmetric phi. Thus the ideal-detector characterization is under-specified unless the uniform-crosstalk condition is made part of the definition.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes the center-of-gravity (COG) algorithm for one-dimensional position measurements in periodic detector arrays. The central result of Section 3.3 is an explicit formula, Eq. (6), for the systematic discretization error of the COG as a Fourier series whose coefficients are values of the Fourier transform of the average signal distribution at multiples of the sampling frequency. The paper then studies conditions under which this error vanishes, derives a resummed form of the series for finite-support signals in Section 5.3, generalizes the analysis to crosstalk and non-perfect spatial integrators in Section 6, treats finite sensor sets, noise, and fluctuations, and concludes with a claimed construction of an 'ideal detector' whose crosstalk eliminates the discretization error for any signal distribution.","tokens_in":25501,"tokens_out":4260,"duration_ms":47601,"significance":"If the results hold, the paper provides a useful analytical framework for correcting COG systematic errors in position-sensitive detectors, with direct applications to calorimeters and silicon strip trackers. The derivation of Eq. (6) is clean and is benchmarked against known special cases, including the Euler sine series in the delta-function limit and the Lau-Pyrlik conjecture. The paper also gives explicit simulation formulas. However, the 'ideal detector' claim, which is one of the advertised headline results, is under-specified as stated: the vanishing-moment condition alone is insufficient, and the paper's own derivation in Section 6.6 relies on an additional uniformity condition that is not included in the abstract or in the standalone definition. This is fixable, but it must be corrected before the central claim is accepted as stated.","major_comments":[{"comment":"The definition of an ideal detector as a crosstalk function satisfying sum_n (x - n tau) p(x - n tau, tau_1) = 0 is insufficient for the claimed conclusion x_g(epsilon) = epsilon for any signal distribution. The derivation on the page with Eq. (14) explicitly uses the uniformity condition sum_n p(x - n tau, tau_1) = 1 to cancel the term involving x' after interchanging the sum and integral. Without uniformity, the zero-first-moment condition does not imply x_g = epsilon: for symmetric p supported on [-tau, tau] of the form p(x) = (tau - x) r(x) on [0, tau] with r symmetric about tau/2 but not constant, one has sum_n (x - n tau) p(x - n tau) = 0 identically, but M_0(x) = sum_n p(x - n tau) = tau r(x) is not constant, and the resulting COG is a weighted average of y M_0(y) rather than epsilon. The uniformity condition should be added to the definition of ideal detector and to the abstract.","section":"Section 6.6 and abstract"},{"comment":"The summation leading to Eq. (13) is compressed, and the sawtooth function is defined inconsistently: the text states Theta(xi) = xi - floor(xi + 1/2), while Eq. (13) states Theta(xi) = xi - floor(xi - 1/2). These differ on half-integers and affect numerical implementation. The authors should specify the exact convention used in Eq. (13) and, if both forms are acceptable up to integer shifts, state that explicitly and verify that the residues are unchanged.","section":"Section 5.3, Eq. (13)"},{"comment":"The reference to 'Eq. (24)' in the ideal-detector derivation is incorrect: Eq. (24) is the probability-density transformation formula, not the expression for x_g(epsilon). The derivation uses the expression for x_g following from Eq. (14), so the equation numbering should be corrected and the derivation should be made self-contained.","section":"Section 6.6"}],"minor_comments":[{"comment":"The condition 'phi(x) = 0 for |x| <= D/2' should read 'phi(x) = 0 for |x| > D/2', or equivalently 'phi(x) != 0 for |x| < D/2', to match the intended finite-support statement.","section":"Section 5.1"},{"comment":"The section numbering in the introductory text does not match the actual structure of the paper; for example, the text refers to 'Section 2' and 'Section 3' in ways that do not correspond to the displayed section headings.","section":"Introduction"},{"comment":"The phrase 'the average distribution omega (x)' appears to be a typo for phi(x), and 'Molier radius' should be 'Moliere radius' throughout the figure captions and text.","section":"Section 7.2"}],"recommendation":"major_revision","confidential_remarks":"This manuscript is a 2019 arXiv posting of a paper published in 2002 in Nucl. Instrum. Methods A. The editor may wish to verify that the 2019 additions, including the 'ideal detector' definition, constitute sufficiently new material relative to the published version to warrant this submission, and that the self-citations to the author's later work are presented appropriately. The ideal-detector issue identified in the main report is a genuine gap in the current text, but it is easily repaired by adding the uniformity condition to the definition and abstract."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things. First, the central derivation is solid: Eq. (6), the Fourier-series form of the COG systematic error, follows cleanly from Poisson summation and the stated shift-invariant infinite-array assumptions, and it genuinely generalizes the Lau–Pyrlik conjecture to arbitrary symmetric signal distributions. Second, the paper's headline claim about an 'ideal detector' is under-specified as written. The 2019-added definition and the abstract state only the vanishing first moment condition, sum_n (x - n tau) p(x - n tau) = 0, but the proof in Section 6.6 also uses the uniform-crosstalk condition sum_n p(x - n tau) = 1 to derive xg = epsilon. The stress-test counterexample is correct: a symmetric p with zero first lattice moment but nonconstant sum_n p(x - n tau) does not give xg = epsilon for extended signals. The triangular crosstalk example satisfies both conditions, so the result survives, but the definition needs repairing.\n\nWhat is genuinely new and good: the derivation of Eq. (6) is nontrivial and the paper carefully benchmarks it against known special cases. The observation that certain crosstalk shapes, like triangular crosstalk of range 2 tau, can eliminate the discretization error for any signal distribution is a real and useful result, and the paper is honest about the assumptions (infinite periodic array, shift-invariant response) under which it holds. The treatment of finite sensor sets and the two-/three-strip discontinuities is also valuable and goes beyond what I have seen in the usual COG literature.\n\nSoft spots, in proportion: the ideal-detector gap is the most important, but it is a one-line fix. The sawtooth definition in Eq. (13) is inconsistent with the text (floor argument sign), and the residue-summation step leading to Eq. (13) is compressed to the point of being hard to check. No code or data accompany the paper, so independent reproduction is on the reader. The 2019 introduction is heavy on self-citation and editorializing, but the core mathematics does not depend on those later papers.\n\nWho this is for: anyone working on position reconstruction in segmented detectors, calorimetry, or silicon trackers. It deserves a serious referee, and I would send it to one, with the expectation that the author tightens the ideal-detector definition and fixes the typo. If the referee does that work, the paper is a useful contribution to the literature.","headline":"A correct and useful derivation of the COG discretization error, with one real definitional gap in the 'ideal detector' section that should be fixed before this is relied upon.","tokens_in":26042,"tokens_out":1470,"would_cite":true,"duration_ms":16991,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["07.05.Kf","06.30.Bp","42.30.Sy"],"model":"deepseek-v4-flash","headline":"Centroiding error is one Fourier series—and crosstalk can kill it","keywords":["center of gravity","centroiding","position measurements","discretization error","crosstalk","ideal detector","Fourier transform","sampling theorem"],"falsifier":"Scan a known, sharp source across one pitch of a strip or pixel detector, measure both the mean signal shape and the centroid $x_g(\\varepsilon)$, and compare the residual $x_g-\\varepsilon$ with the series in Eq. (6) using the independently measured Fourier transform $\\Phi$. Agreement within noise supports the model; systematic position-dependent deviations that do not follow the predicted coefficients would falsify it. A separate direct test: a detector with triangular crosstalk of range $2\\tau$ should show $x_g=\\varepsilon$ to within noise for every signal distribution; any periodic wiggle in the residual rules out the ideal-detector claim.","tokens_in":25058,"feed_emoji":"🎯","tokens_out":12948,"duration_ms":111050,"temperature":0.7,"pith_summary":"The center-of-gravity (COG) algorithm estimates a particle's impact point as the signal-weighted average of sensor positions. This paper derives, from first principles, exactly how that estimate deviates from the true position when the signal is collected by an infinite periodic row of integrating detectors. The deviation is a Fourier sine series in the impact position $\\varepsilon$, with coefficients given by the Fourier transform $\\Phi(2\\pi k/\\tau)$ of the average signal distribution at the sampling harmonics. That closed form turns the usual empirical centroid correction into a calculable, invertible relation, and it isolates the special cases where the error vanishes: band-limited signals, particular rectangular or triangular signal shapes, and sensors with a triangular crosstalk of range twice the pitch. The same Fourier machinery gives explicit formulas for finite clusters of sensors, crosstalk, noise, and fluctuations, so the math can be used both to correct measurements and to simulate realistic detectors.","feed_headline":"Centroiding error is one Fourier series—and crosstalk can kill it","feed_subtitle":"Exact corrections for strip and crystal detectors follow from the signal's spectrum; one crosstalk shape removes the bias entirely.","key_machinery":"The load-bearing object is the sampled signal train $s_\\varepsilon(x)=f(x-\\varepsilon)\\sum_n\\delta(x-n\\tau)$, whose Fourier transform is evaluated by the Poisson identity $S_\\varepsilon(\\omega)=\\tau^{-1}\\sum_k F_\\varepsilon(\\omega-2\\pi k/\\tau)$. Because the integrating detector converts the signal to $f(x-\\varepsilon)$ through a convolution with a rectangular window, $F_\\varepsilon(\\omega)=[2\\sin(\\omega\\tau/2)/\\omega]\\Phi(\\omega)e^{-i\\varepsilon\\omega}$, and the COG formula $x_g=i S_\\varepsilon'(0)/S_\\varepsilon(0)$ turns the derivative into the sine series of Eq. (6). The same machinery is reused in two further forms: the dual Shannon/WKS series expresses $\\Phi$ through the form factors $\\Phi(2\\pi n/D)$, and a contour-residue summation converts the resulting double series into the practical expression Eq. (13). For crosstalk, the response function $p$ enters through its Fourier transform $P$ and its derivative, and the condition for an ideal detector is stated directly as a vanishing weighted sum of shifted response copies.","core_discovery":"On the paper's own terms, the central result is Eq. (6): for an infinite periodic array of ideal integrating sensors with pitch $\\tau$, the COG estimate is $x_g = \\varepsilon + (\\tau/\\pi)\\sum_{k\\ge1}(-1)^k k^{-1}\\sin(2\\pi k\\varepsilon/\\tau)\\Phi(2\\pi k/\\tau)$, where $\\varepsilon$ is the true impact position and $\\Phi$ is the Fourier transform of the mean signal distribution. The discretization error is therefore not a statistical nuisance but a deterministic function of the signal spectrum at the sampling frequencies. The paper further proves that the error vanishes when $\\Phi(2\\pi k/\\tau)=0$ for all $k>0$, a condition broader than the Whittaker-Kotelnikov-Shannon band limit, and exhibits finite-support shapes such as rectangular signals with integer-multiple widths and triangular signals with even-multiple widths that satisfy it. It then extends the analysis to crosstalk: if the detector response $p(x,\\tau_1)$ satisfies $\\sum_n (x-n\\tau)p(x-n\\tau,\\tau_1)=0$, in particular triangular crosstalk of range $2\\tau$, the COG equals $\\varepsilon$ for any signal distribution. Finite clusters of sensors are treated too, and they produce discontinuities in $x_g(\\varepsilon)$ at the points where the set of sensors changes.","pith_inferences":["The ideal-detector condition suggests a concrete design target for silicon strip and pixel detectors: tune inter-strip capacitive coupling toward a triangular profile with range exactly twice the pitch, so centroids need no correction at all; this is an engineering consequence the paper states mathematically but does not develop as a fabrication recipe.","Because only the values $\\Phi(2\\pi k/\\tau)$ enter the correction, the scheme should transfer directly to imaging and star-tracker centroiding wherever the point-spread function is known, not just to particle detectors.","The predicted infinite peak in the $x_g$ probability density when the signal support $D\\le\\tau$ is a testable signature: a histogram of centroids from a source narrower than one pixel should show a sharp spike near zero, distinguishing this model from generic smoothing.","The boundary-discontinuity analysis implies that any empirical bias curve subtracted from cluster centroids is only valid for a fixed sensor count; switching cluster sizes in the same analysis introduces uncorrectable systematic jumps."],"forward_implications":["For any detector that matches the infinite periodic model, the systematic position error can be removed by inverting Eq. (6) or using Eq. (13), without per-detector Monte Carlo tuning.","A detector engineered to have triangular crosstalk of range $2\\tau$, or any response satisfying $\\sum_n (x-n\\tau)p(x-n\\tau,\\tau_1)=0$, would be unbiased for every signal distribution.","Cluster-size cuts create discontinuities in $x_g(\\varepsilon)$ at sensor boundaries; subtracting a fixed bias mixes incompatible two-sensor and three-sensor algorithms and cannot be fixed by a smooth polynomial correction.","Given a large sample of equivalent hits with a uniform distribution of true positions, the histogram of $x_g$ yields $\\varepsilon(x_g)$ by Eq. (25), and when the signal support fits within one pitch it also reconstructs the signal shape.","The explicit formulas for crosstalk, noise, inter-strip calibration errors, cracks, and off-axis crystals make the same equations a compact simulation tool for real strip and crystal calorimeters."],"supporting_citations":[{"why":"It supplies the Fourier-transform and convolution properties used to obtain the COG derivative formula and the transform of the integrating window.","marker":"[8]"},{"why":"It provides the Poisson identity that turns the sampled signal train into an alias sum, and the unitary periodic function notion used for crosstalk.","marker":"[9]"},{"why":"It states the Shannon/WKS sampling theorem and its dual, used to connect the form factors with the Fourier transform values at sampling harmonics.","marker":"[11]"},{"why":"It gives the earlier conjectured periodic form of the COG error and reports discontinuities for finite sensor sets that Eq. (6) reproduces and explains.","marker":"[10]"},{"why":"It introduces the two-strip COG and the eta-algorithm correction that the paper generalizes to arbitrary signal widths and to probability-distribution inversion.","marker":"[16]"},{"why":"It supplies an example of extracting form factors from experimental data, which the paper's histogram relations formalize.","marker":"[12]"},{"why":"It supplies the contour-integration and residue techniques used to sum the series in Eq. (12) into the practical form Eq. (13).","marker":"[14]"},{"why":"It provides a silicon-strip detector with unread strips and capacitive crosstalk as a physical example of the uniform-crosstalk class.","marker":"[15]"}],"fun_headline_variants":["Centroid error is a Fourier series; one crosstalk shape cancels it","Exact centroid corrections from signal spectrum, crosstalk fix","Crosstalk can eliminate centroid discretization error for any signal","Centroiding bias solved by Fourier analysis and ideal crosstalk"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole Fourier-series correction rests on treating the detector as an infinite periodic row of identical sensors whose collected signal is exactly a shift-invariant convolution of the mean signal with a fixed response function; real detectors with finite size, varying strip widths or gains, calibration errors, or unread gaps only approximate this model.","fun_headline_variants_meta":{"raw":{"variants":["Centroid error is a Fourier series; one crosstalk shape cancels it","Exact centroid corrections from signal spectrum, crosstalk fix","Crosstalk can eliminate centroid discretization error for any signal","Centroiding bias solved by Fourier analysis and ideal crosstalk"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000201,"raw_usage":{"total_tokens":1396,"prompt_tokens":979,"completion_tokens":417,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":595,"completion_tokens_details":{"reasoning_tokens":341}},"tokens_in":595,"tokens_out":417,"duration_ms":4774,"temperature":1.0,"reasoning_tokens":341,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:42:28.432070+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Scan a known, sharp source across one pitch of a strip or pixel detector, measure both the mean signal shape and the centroid $x_g(\\varepsilon)$, and compare the residual $x_g-\\varepsilon$ with the series in Eq. (6) using the independently measured Fourier transform $\\Phi$. Agreement within noise supports the model; systematic position-dependent deviations that do not follow the predicted coefficients would falsify it. A separate direct test: a detector with triangular crosstalk of range $2\\tau$ should show $x_g=\\varepsilon$ to within noise for every signal distribution; any periodic wiggle in the residual rules out the ideal-detector claim.","supporting_citations":[{"cited_title":"Bracewell, The Fourier Transform and Its Applicati on (McGraw-Hill, New Y ork, NY , 1986)","cited_arxiv_id":null,"evidence_quote":"It supplies the Fourier-transform and convolution properties used to obtain the COG derivative formula and the transform of the integrating window."},{"cited_title":"Champeney, A Handbook of Fourier Theorems (Cambrid ge University Press, Cambridge, UK, 1987)","cited_arxiv_id":null,"evidence_quote":"It provides the Poisson identity that turns the sampled signal train into an alias sum, and the unitary periodic function notion used for crosstalk."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It states the Shannon/WKS sampling theorem and its dual, used to connect the form factors with the Fourier transform values at sampling harmonics."},{"cited_title":"Lau and J","cited_arxiv_id":null,"evidence_quote":"It gives the earlier conjectured periodic form of the COG error and reports discontinuities for finite sensor sets that Eq. (6) reproduces and explains."},{"cited_title":"Belau et al., Nucl","cited_arxiv_id":null,"evidence_quote":"It introduces the two-strip COG and the eta-algorithm correction that the paper generalizes to arbitrary signal widths and to probability-distribution inversion."},{"cited_title":"Lednev, Nucl","cited_arxiv_id":null,"evidence_quote":"It supplies an example of extracting form factors from experimental data, which the paper's histogram relations formalize."},{"cited_title":"Henrici, Applied Computational Complex Analysis (J ohn Wiley, New Y ork, NY , 1974)","cited_arxiv_id":null,"evidence_quote":"It supplies the contour-integration and residue techniques used to sum the series in Eq. (12) into the practical form Eq. (13)."},{"cited_title":"Alpat et al., Nucl","cited_arxiv_id":null,"evidence_quote":"It provides a silicon-strip detector with unread strips and capacitive crosstalk as a physical example of the uniform-crosstalk class."}],"review_version":1}