Pith. Sign in

REVIEW

Measurement of SUSY masses via cascade decays for SPS 1a

PaperPDF from arXiv

data files

ReviewThe objections and verdict, as JSON Paper breakdownThe authors' claims, conditions and measurements, as JSON Signed receiptProof this review has not been altered
X Bluesky LinkedIn Reddit Hacker News
B. K. Gjelsten, D. J. Miller, P. Osland

High Energy Physics - Phenomenology · Submitted 22 Oct 2004 · 50 references

Pith's verdict: A durable methods paper on inverting cascade endpoints into sparticle masses, with one honest caveat that turns out to be sharper than the authors let on.

T0 · AI review

How this review was made: An AI referee (selfhosted/deepseek-v4.1-flash) read the paper on November 11, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 3 major objections · 7 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDhep-ph/0410303 (version 2)
Pith ID2VRTSVWB, Pith's permanent record of this paper and its review
CitationB. K. Gjelsten, D. J. Miller, P. Osland. "Measurement of SUSY masses via cascade decays for SPS 1a." arXiv:hep-ph/0410303 (2004). https://pith.science/paper/2VRTSVWB
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report

How Pith checks papers → · Full audit of this paper →

The paper

The paper claims

Superpartner masses can be read off the edges of LHC decay distributions, even when the lightest particle escapes unseen.

The reading

At the LHC, if supersymmetry is real, new heavy particles will be produced in pairs and decay down through a chain of lighter partners, ending in a lightest particle that leaves no trace in the detector. The paper asks whether the masses of the unseen chain can still be recovered from the shapes of the visible decay products alone. Its answer is yes: the sharp edges and thresholds of invariant-mass distributions built from one quark and two leptons sit at positions fixed by the four masses in the chain, and those relations can be turned around to give the masses back. The authors write the inversion formulae explicitly, check them against simulated data at two points on a standard benchmark line, and show that a future electron-positron collider measurement of the lightest particle's mass would remove the remaining overall scale uncertainty. The work is a careful stress test of a method rather than a discovery claim: it maps where the method applies, and where it quietly breaks.

What carries the argument

The load-bearing device is the region-by-region inversion of the endpoint formulae. Each endpoint is a single-valued function of the masses on each patch of mass space, but the inverse is not single-valued: one set of measured endpoints can correspond to several mass sets, and the paper treats this explicitly rather than hiding it. The masses are then extracted by minimising $Σ = [E_{exp} - E_{th}(m)]^T W [E_{exp} - E_{th}(m)]$, where the weight matrix $W$ is the inverse of the error matrix and carries the correlations induced by a common jet energy scale. Because $Σ$ is a composite of nine branch functions, its minima must be checked for lying in their own region of validity; minima on region borders and 'unphysical' minima in the wrong region are catalogued and weighed. A fifth measurement, the threshold of the $m_{qll}(θ>à ¡/2)$ distribution, overconstrains the system and is used to discriminate among the competing minima.

Extended reading notes

Core claim

The paper establishes that the decay chain $ñq_L → Ç_2^0 q → é_R l q → Ç_1^0 l l q$ leaves four measurable kinematic endpoints in the invariant-mass distributions $m_{ll}$, $m_{qll}$, $m_{ql}(low)$ and $m_{ql}(high)$, and that these four numbers can be inverted into the four sparticle masses. The inversion is not a single formula but a patchwork: the endpoint expressions take different algebraic forms depending on the ratios $m_{ñq_L}/m_{Ç_2^0}$, $m_{Ç_2^0}/m_{é_R}$ and $m_{é_R}/m_{Ç_1^0}$, giving nine regions of mass space, of which six admit an inversion and three are degenerate because $(m^{max}_{qll})^2 = (m^{max}_{ql(high)})^2 + (m^{max}_{ll})^2$ there. The authors give the explicit formulae for all six usable regions, then simulate 300 fb$^{-1}$ of ATLAS data at two benchmark points to see how well the endpoints can actually be measured. At the lighter point the recovered masses track the true ones closely, with the three lightest sparticles strongly correlated and mass differences determined far better than masses. At the heavier point the method degrades: the least-squares function develops two or three competing minima, the ensemble distributions become skewed, and a small unmodelled 'foot' at the end of the $m_{ql}(high)$ edge is shown to shift the recovered masses by tens of GeV.

Load-bearing premise

The precision claims assume that the systematic uncertainty of fitting the edges can be driven below the statistical and energy-scale errors, even though the paper's own table lists fit systematics comparable to or larger than those errors.

Editorial extensions

If this is right

  • The LHC can measure superpartner masses without reconstructing the invisible lightest particle, using only the endpoints of one-quark-two-lepton mass distributions.
  • Mass differences within the chain ($m_{é_R} - m_{Ç_1^0}$, $m_{Ç_2^0} - m_{Ç_1^0}$, $m_{ñq_L} - m_{Ç_1^0}$) are determined much more precisely than the masses themselves, because the overall scale drops out.
  • A linear-collider measurement of the lightest neutralino mass to about 0.05 GeV fixes the scale and collapses the LHC uncertainties by roughly an order of magnitude, and removes the skewness and the false minima.
  • Over much of the mSUGRA parameter space the required mass hierarchy and a usable branching ratio are both present, so the method is not tied to one benchmark point.
  • The same inversion machinery, extended by one more quark momentum, should apply to chains headed by a gluino.

Reading between the lines

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

  • The paper's own diagnosis of the mismeasured $m_{ql}(high)$ edge suggests a general self-checking loop: take the fitted masses, regenerate the theory distributions, and look for structure the fit failed to include. That loop is a reusable recipe for any endpoint analysis, not just this chain.
  • The recurring problem is that the map from masses to endpoints is many-to-one. A natural extension is to treat the whole measured distribution, not just its endpoint, as the observable, which would break the degeneracies that four numbers cannot.
  • Because the three lightest masses come out strongly correlated, publishing the full error matrix rather than per-mass errors is what downstream work on soft-breaking parameters would actually need.
  • The border effects described here are a generic feature of any analysis whose fitting function is piecewise in the parameters, and would appear in other cascade chains with the same structure.
Original abstract

If R-parity conserving supersymmetry exists below the TeV-scale, new particles will be produced and decay in cascades at the LHC. The lightest supersymmetric particle will escape the detectors, thereby complicating the full reconstruction of the decay chains. In this paper we expand on existing methods for determining the masses of the particles in the cascade from endpoints of kinematical distributions. We perform scans in the mSUGRA parameter space to delimit the region where this method is applicable. From the examination of theoretical distributions for a wide selection of mass scenarios it is found that caution must be exerted when equating the theoretical endpoints with the experimentally obtainable ones. We provide analytic formulae for the masses in terms of the endpoints most readily available. Complications due to the composite nature of the endpoint expressions are discussed in relation to the detailed analysis of two points on the SPS~1a line. Finally we demonstrate how a Linear Collider measurement can improve dramatically on the precision of the masses obtained.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13-20260917 via OpenRouter, question bank abstract-p1-20260923, 2026-09-23: 1 request, 6,179 input tokens, $0.0003 reported.

What the abstract itself states, as probabilities: claims first 0.03 · states limitation 0.91 · human participants 0.01 · controlled comparison 0.04 · randomized assignment 0.01 · states quantitative result 0.07 · compares with prior results 0.50 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v2. Results and the main claim are underlined; .

If R-parity conserving supersymmetry exists below the TeV-scale, new particles will be produced and decay in cascades at the LHC. The lightest supersymmetric particle will escape the detectors, thereby complicating the full reconstruction of the decay chains. In this paper we expand on existing methods for determining the masses of the particles in the cascade from endpoints of kinematical distributions. We perform scans in the mSUGRA parameter space to delimit the region where this method is applicable. From the examination of theoretical distributions for a wide selection of mass scenarios it is found that caution must be exerted when equating the theoretical endpoints with the experimentally obtainable ones. We provide analytic formulae for the masses in terms of the endpoints most readily available. Complications due to the composite nature of the endpoint expressions are discussed in relation to the detailed analysis of two points on the SPS~1a line. Finally we demonstrate how a Linear Collider measurement can improve dramatically on the precision of the masses obtained.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 If R-parity conserving supersymmetry exists below the TeV-scale, new particles will be produced and decay in cascades at the LHC. background 1.00 0.05 0.03 hedged possibility 0.67
s1 The lightest supersymmetric particle will escape the detectors, thereby complicating the full reconstruction of the decay chains. background 0.99 0.05 0.02 asserted 0.99
s2 In this paper we expand on existing methods for determining the masses of the particles in the cascade from endpoints of kinematical distributions. contribution 0.98 0.60 0.10 asserted 0.96
s3 We perform scans in the mSUGRA parameter space to delimit the region where this method is applicable. method 0.79 0.94 0.04 asserted 0.97
s4 From the examination of theoretical distributions for a wide selection of mass scenarios it is found that caution must be exerted when equating the theoretical endpoints with the experimentally obtainable ones. result 0.68 0.90 0.24 asserted 0.99
s5 We provide analytic formulae for the masses in terms of the endpoints most readily available. contribution 0.91 0.94 0.14 asserted 1.00
s6 Complications due to the composite nature of the endpoint expressions are discussed in relation to the detailed analysis of two points on the SPS~1a line. contribution 0.65 0.76 0.06 asserted 0.91
s7 Finally we demonstrate how a Linear Collider measurement can improve dramatically on the precision of the masses obtained. result 0.51 0.84 0.35 asserted 1.00

What kind of paper the classifier reports this is: method or tool 0.67 · evidence from simulation or computation 0.79.

What was examined, and what was not

sources read
  • arXiv listing abstract v2
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. Classification of different hierarchies, labeled (i)–(vii), for four combinations of tan β and A0, such that the four panels contain respectively the SPS 1a line, the SPS 1b point, the SPS 3 line, and the SPS 5 point. The regions marked ‘TF’ are theoretically forbidden. (See text for details.) (ii) ˜g > max( ˜dL, u˜L, ˜b1,t˜1) and ˜lR > χ˜ 0 2 > τ˜1 (iii) ˜g > max( ˜dL, u˜L, ˜b1,t˜1) and min(˜lR, τ˜1) > χ˜ 0 2 (iv) ˜dL > g >˜ max(˜uL, ˜b1) and min(˜lR, τ˜1) > χ˜ 0 2 (v) min( ˜dL, u˜L) > g >˜ ˜b1 and min(˜lR, τ˜1) > χ˜ 0 2 (vi) min( ˜dL, u˜L, ˜b1) > g >˜ t˜1 and min(˜lR, τ˜
Figure 2. Decay channels of ˜g with A0 = −m0, tan β = 10 (left) and A0 = 0, tan β = 30 (right). In the left panel the SPS 1a line is shown together with the two points (α) and (β) marked with triangles. In the right panel the triangle marks the SPS 1b point. The branching ratio of a given decay channel in a small neighbourhood of the m1/2–m0 plane is equal to the fraction which the corresponding colour occupies in that neighbourhood. The region where ˜χ 0 1 is not the LSP is denoted ‘Charged LSP’ and is discarded. Some regions are also forbidden theoretically, in that e.g. it is not possible…
Figure 3. Decay channels of ˜uL (top) and ˜uR (bottom) with A0 = −m0, tan β = 10 (left) and A0 = 0, tan β = 30 (right). (See the caption of
Figure 4. Decay channels of ˜b1 (upper panels) and ˜b2 (lower panels) with A0 = −m0, tan β = 10 (left) and A0 = 0, tan β = 30 (right). (See the caption of
Figure 5. Decay channels of ˜χ 0 2 with A0 = −m0, tan β = 10 (left) and A0 = 0, tan β = 30 (right). (See the caption of
Figure 6. Masses (left) and relative widths (right) of relevant sparticles as m1/2, m0 and A0 are varied along the SPS 1a slope, defined by Eq. (3.1). The vertical dotted lines represent SPS 1a points (α) and (β). In addition, we define two points (α) and (β) on the SPS 1a line according to: (α) : m0= 100 GeV, m1/2 = 250 GeV, (β) : m0= 160 GeV, m1/2 = 400 GeV. (3.2) The first point (α) is the ‘basic’ SPS 1a point of Ref. [21] and studied in Ref. [27, 28], while the second (β) is a new, less optimistic scenario with a reduced cross-section for the decay chain. The masses of particles relevant…
Figure 7. Cross-sections as m1/2, m0 and A0 are varied along the SPS 1a slope, defined by Eq. (3.1). The vertical dotted lines represent SPS 1a points (α) and (β). pair production cross-sections, as m1/2 is varied along the SPS 1a line. Notice that these cross-sections fall very rapidly as m1/2 is increased, which will cause repercussions in the analysis of SPS 1a (β). The cross-sections for gluino–gluino, gluino–squark and squark–squark pair produc￾tions are detailed in
Figure 8. Branching ratios of ˜χ 0 2 as m1/2, m0 and A0 are varied along the SPS 1a slope. The vertical dotted lines represent SPS 1a points (α) and (β). ( )  = 32:8 pb q~L 31:4%  = 7:7 pb ~b1 35:5%  = 4:3 pb ~b2 18:0% 9 >= >; ~0 2 12:1% ~lR 100% ~0 1 1245 fb 329 fb 94 fb 1669 fb ( )  = 3:21 pb q~L 31:9%  = 0:50 pb ~b1 23:6%  = 0:31 pb ~b2 8:8% 9 >= >; ~0 2 5:3% ~lR 100% ~0 1 54:3 fb 6:3 fb 1:5 fb 62:0 fb q ln lf
Figure 9. The SPS 1a cascade decay chain. 3.3 The cascade The cross-sections and branching ratios of our chosen decay chains q˜L → χ˜ 0 2 q → ˜l ∓ R l ± n q → χ˜ 0 1 l ∓ f l ± n q (3.3) ˜b1 → χ˜ 0 2 b → ˜l ∓ R l ± n b → χ˜ 0 1 l ∓ f l ± n b (3.4) ˜b2 → χ˜ 0 2 b → ˜l ∓ R l ± n b → χ˜ 0 1 l ∓ f l ± n b (3.5) are summarized in
Figure 10. Theoretical mass distributions for SPS 1a (α) and (β), as well as for three other mass scenarios, denoted (i), (ii) and (iii). Kinematic endpoints are given in units of mmax qll . (More details will be given in [42].) 4.1 Theory curves of invariant mass distributions In
Figure 11. Different-flavour subtraction for mll and mqll at (α) (left) and (β) (right). The solid/dashed red curves are the same-flavour (‘SF’)/different-flavour (‘DF’) distributions. In black, their difference, the different-flavour-subtracted distribution (‘SF-DF’), is shown with error bars. The blue curve shows the part of the subtracted distribution which contains a signal chain (‘SC’). The solid/dashed green curves (‘SM’) give the Standard Model part (completely dominated by tt¯) of the same-flavour/different-flavour distributions. They are statistically equal and will cancel each…
Figure 12. Quark selection cuts for mqll (left) and mql(low) (right). Left: The two lower black curves in solid and dashed show the mqlowll (‘qlow’) and mqhighll (‘qhigh’) distributions. The highest black curve shows the sum of the two, mqbothll (‘qboth’). Notation is defined in the text. Blue: theory distribution (‘TH’). Red: inconsistency cuts (‘incons. cut’) have been applied on mql(low) and mql(high) (the resulting distribution is scaled up). Right: An inconsistency cut is applied on mqll. The jet chosen for mqlowll is the one used in mqlowl(low) etc. The inconsistency cut is seen to work…
Figure 13. Invariant mass distributions for SPS 1a (α). See the text for details. other distributions, it does not do so well. Nevertheless, there are many ways to construct the mixed sample. Further study of the method seems worthwhile. – 33 –
Figure 14. Invariant mass distributions for SPS 1a (β). See the text for details. mqll: The mqlowll distribution is shown in
Figure 15. Sparticle masses and mass differences at SPS 1a (α) for solutions with ∆Σ ≤ 1. The unfilled distributions in black show from left to right mχ˜ 0 1 , m˜lR , mχ˜ 0 2 , m˜b1 and mq˜L for solutions in the nominal region (1,1). We will have such a solution in 94% of the experiments, see
Figure 16. Sparticle masses and mass differences at SPS 1a (β). All masses of solutions with ∆Σ ≤ 1 which lie in regions (1,2), (1,3) and on their common border are shown. From left to right the unfilled distributions show mχ˜ 0 1 , m˜lR , mχ˜ 0 2 and mq˜L . The filled distributions show the narrower mass differences m˜lR − mχ˜ 0 1 , mχ˜ 0 2 − mχ˜ 0 1 and mq˜L − mχ˜ 0 1 . Skewness of mass distributions is visible. small, but for (β) the effect is large. The reason why we naively would expect a symmetric distribution around the nominal masses in the first place, is that the endpoint…
Figure 17. Border effect between region (1,2) and (1,3); the mass of ˜χ 0 1 as a function of the border parameter b, see Eq. (6.4), demonstrating border effects. Filled boxes represent physical solutions, empty boxes represent unphysical solutions. when the unphysical (1,2) solutions are replaced by the physical (1,3) solutions which lie at higher masses, the ensemble mean increases. Here, this effect brings the average value for mχ˜ 0 1 to 183 GeV, an additional increase of nearly 20 GeV. When the threshold measurement is added, the border effect is reduced, giving an average value of 173 G

Figures reproduced from arXiv:hep-ph/0410303 by the authors.

References
How later papers used it

49 later papers cite it

Evidence: substantial. 48 of the citing papers share no author with it, and 13 of those carry a confident label from Pith's citation classifier.

Pith read 26 of the 115 citing sentences, choosing only those whose citation markers point at this paper alone, and wrote the paraphrases from them. The summary, the counts, the evidence level, the links and the quotes come from code; the summary counts the papers listed below. A later paper is listed under a label when Pith's citation classifier gives one of its sentences that label at probability 0.8 or more (0.95 for support), Pith's reading finds that the citing paper itself did it, and a second check confirms the paraphrase against the sentence. This section lists uses only. It does not list criticism or disagreement, so their absence here says nothing about how the paper has been received. Papers that share an author with this one are listed apart, under "By the same authors". A quote is the citing authors' own sentence, shown when their paper's licence is CC BY, CC BY-SA or CC0. JSON

Uses its method or tool 5 of 12 papers with this label
Uses its data or values 1 paper
How the literature uses it

46 arXiv papers cite this in 107 sentences

Each citing sentence is taken from the citing paper's LaTeX source and typed by Pith's classifier. A label appears at probability 0.8 or more (0.95 for support), a contest at 0.7 or more; the quoted words are the citing authors'. How each label was checked · JSON

Papers citing it · 63 sentences without a confident label

Uses its method or tool 10 papers
Uses its data or values 1 paper
  • [ht]|c||r|r|r||r|r||r|r| 3c||SPS1a($α$) [ref] 2c|SU1 [ref] 2c|SU3 [ref] [1mm] 2-8 Variable Nominal True False True False True False [1mm] $m_χ̃^0_1$ 96.1 96.3 85.3 137.0 122.1 118.0 346.8 [1mm] $m_l̃_R$ 143.0 143.2 130.4 254.0 127.5 155.0 411.1 [1mm] $m_χ̃^0_2$ 176.8 177.0 165.5 264.0 245.9 219.0 451.6 [1mm] $m_q̃_L$ 537.2 537.5 523.5 760.0 743.6 631.0 899.9 [1mm] Region (1,1) (1,1) (1,2) (1,1) (1,3) (1,3) (1,1) [1mm] $m_ll^max$ 77.0 77.0 77.1 2c||61 2c|100 [1mm] $m_jl(lo)^max$ 298.3 298.3 299.6 2c||194 2c|322 [1mm] $m_jl(hi)^max$ 375.6 375.6 375.7 2c||600 2c|418 [1mm] $m_jll^max$ 425.8 425.8 425.6 2c||609 2c|499 [1mm] $m_jll(θ>π/2)^min$ 200.6 200.6 205.1 143 148 247 214 [1mm]

    Using kinematic boundary lines for particle mass measurements and disambiguation in SUSY-like events with missing energy · Introduction · 1.00

Compares with it 2 papers
Cites it as support 1 paper
  • As stressed in Ref. [ref] , even in favorable cases applying analytic functions too readily may lead to inaccurate measurements and underestimated errors, since endpoints can often exhibit tails or get smeared.

    Pinning down the Invisible Sneutrino · Sneutrino mass determination from on-shell kinematics · 1.00

Cites it as background 22 papers

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.