650 The ATLAS Collaboration / Physics Letters B 760 (2016) 647–665
Fig. 3. Data (black dots) and background estimates (red solid line) for m
β
(left) and m
βγ
(right) for the gluino R-hadron search (1000 GeV). The green shaded band illustrates
the statistical uncertainty of the background estimate. The blue dashed lines illustrate the expected signal (on top of background) for the given R-hadron mass hypothesis.
The black dashed vertical lines at 500 GeV show the mass selection and the last bin includes all entries/masses above. (For interpretation of the references to colour in this
figure legend, the reader is referred to the web version of this article.)
as well as βγ and β, respectively, via m = p/βγ . The minimum
mass requirements, m
min
βγ
and m
min
β
, are set to correspond to a
value about 2σ below the nominal R-hadron mass value, given
the mass resolution expected for the signal.
The
total selection efficiency depends on the sparticle mass and
varies between 9% and 15% for gluino and top-squark R-hadrons
and
6% to 8% for bottom-squark R-hadrons. The lower efficiency
for bottom squarks is expected, as R-hadrons are most likely pro-
duced
in mesonic states, where those with down-type squarks
tend to be neutral more often than those with up-type squarks,
due to light-quark production ratios of u : d : s ≈ 1 : 1 : 0.3 [12]
during
hadronisation. The expected signal yield and efficiency, es-
timated
background and observed number of events in data for
the full mass range after the final selection are summarised in
Table 3.
5. Background estimation
The background is evaluated in a data-driven manner. First,
probability distribution functions (pdf) in the momentum, and also
in the β and βγ values, are determined from data. These pdfs
are produced from candidates in data, which have passed the ini-
tial
selection mentioned earlier, but fall in sidebands of the signal
region, as described below. Background distributions in m
β
and
m
βγ
are obtained by randomly sampling the pdfs derived above
and then using the equation m = p/βγ . These mass distributions,
which are normalised to the data events outside the signal region
(i.e. not passing both mass requirements of the hypothesis in ques-
tion),
are shown in Fig. 3 along with the data and expected signal
for the 1000 GeV gluino R-hadron mass hypothesis.
Each
R-hadron mass hypothesis has a different selection, and
therefore corresponding individual background estimates are pro-
duced
accordingly. The momentum pdf is produced from events
that pass the momentum cut, but fail the β and βγ requirements
in Table 1 for the chosen R-hadron mass hypothesis, but nonethe-
less
have β<1 and βγ < 2.5. The β and βγ pdfs are produced
by selecting events which pass the respective β and βγ selec-
tion
and have momentum in the range 50 GeV < p < 200 GeV.
Since momentum is correlated with |η|, any correlation between
|η| and β (βγ ) will lead to a correlation between momentum and
Table 2
Summary
of all studied systematic uncertainties. Ranges indicate a dependency on
the R-hadron mass hypothesis (from low to high masses).
Source Relative uncertainty
[±%]
Theoretical uncertainty on signal 14–57
Uncertainty on signal efficiency 20–16
Trigger efficiency 2
QCD uncertainty (ISR, FSR) 14
Pile-up 7–1
Pixel βγ measurement 1–3
Calorimeter β measurement 10–2
Luminosity 5
Uncertainties on background estimate 30–43
β (βγ ), invalidating the background estimate. The size and impact
of such correlations are reduced by determining the three pdfs in
five equal-width bins of |η|. This procedure also ensures that dif-
ferent
detector regions are treated separately.
6. Systematic uncertainties
The systematic uncertainties are obtained from data, whenever
possible. The two major uncertainties for which this is not the case
are cross sections and ISR, the latter being folded with the trigger
efficiency curve obtained from data to produce the overall E
miss
T
trigger efficiency. The individual contributions are outlined below
and summarised in Table 2.
6.1. Theoretical cross sections
Signal cross sections are calculated to next-to-leading order in
the strong coupling constant, including the resummation of soft
gluon emission at next-to-leading-logarithmic accuracy (NLO +
NLL) [38–40]. The nominal cross section and the uncertainty are
taken from an envelope of cross-section predictions using differ-
ent
PDF sets and factorisation and renormalisation scales, as de-
scribed
in Ref. [41]. This prescription results in an uncertainty of
14% (at 600 GeV) rising to 24% (at 1600 GeV) and to 32% (at
2000 GeV) for gluino R-hadrons and marginally larger values for
squark R-hadrons.