Search for metastable heavy charged particles with large ionization energy loss in pp collisions at $\sqrt{s} = 13$ TeV using the ATLAS experiment

This paper presents a search for massive charged long-lived particles produced in pp collisions at $\sqrt{s}=$ 13 TeV at the LHC using the ATLAS experiment. The dataset used corresponds to an integrated luminosity of 3.2 fb$^{-1}$. Many extensions of the Standard Model predict the existence of massive charged long-lived particles, such as $R$-hadrons. These massive particles are expected to be produced with a velocity significantly below the speed of light, and therefore to have a specific ionization higher than any Standard Model particle of unit charge at high momenta. The Pixel subsystem of the ATLAS detector is used to measure the ionization energy loss of reconstructed charged particles and to search for such highly ionizing particles. The search presented here has much greater sensitivity than a similar search performed using the ATLAS detector in the $\sqrt{s}=$ 8 TeV dataset, thanks to the increase in expected signal cross-section due to the higher center-of-mass energy of collisions, to an upgraded detector with a new silicon layer close to the interaction point, and to analysis improvements. No significant deviation from Standard Model background expectations is observed, and lifetime-dependent upper limits on $R$-hadron production cross-sections and masses are set. Gluino $R$-hadrons with lifetimes above 0.4 ns and decaying to $q\bar{q}$ plus a 100 GeV neutralino are excluded at the 95% confidence level, with lower mass limit ranging between 740 GeV and 1590 GeV. In the case of stable $R$-hadrons the lower mass limit at the 95% confidence level is 1570 GeV.


Introduction
Massive, long-lived particles (LLPs) are predicted by a wide range of physics models that extend the Standard Model (SM).LLPs arise in proposed solutions to the gauge hierarchy problem [1], including supersymmetric (SUSY) models that either violate [2][3][4] or conserve [5][6][7][8][9][10][11][12] R-parity.The lifetime of these particles depends on the mass difference between the particle and the lightest stable SUSY particle, and on the size of any R-parity-violating coupling [13].
Because of their large mass, LLPs are expected to be slow (β significantly below 1) and, if charged, to have a specific ionization higher than any SM particles of unit charge at high momenta.The Pixel subsystem [14] of the ATLAS detector [15] provides measurements of ionization energy loss (dE/dx) for charged particles and can be used to distinguish such highly ionizing particles from SM particles.In this paper, this information is used to search for LLPs using data collected at √ s = 13 TeV, extending in reach a study performed on data collected at √ s = 8 TeV [16].The increase in probability for highenergy parton-parton collisions significantly enhances the production cross-section of LLPs via the strong interaction.Together with this enhancement, the search sensitivity has been considerably extended thanks to the upgrades to the Pixel detector (see Section 2) and improvements to the analysis (see Section 6).
Many strategies have been explored to find LLPs at the LHC.A search using the complete √ s = 8 TeV dataset recorded by the ATLAS detector [17] combined Pixel dE/dx information with time-of-flight measurements, made in both the calorimeters and the muon spectrometer, to search for LLPs that do not decay within the detector.Charginos that are nearly mass degenerate with the lightest stable particle have been searched for in ATLAS looking for disappearing tracks [18].A displaced vertex search [19] performed by the ATLAS Collaboration provided limits on the mass of gluinos decaying due to couplings violating both R-parity and lepton number conservation.The CMS Collaboration published an analysis searching for nondecaying LLPs based on dE/dx, β, and muon identification [20].Additional analyses at the CMS experiment searched for LLPs with short lifetimes by looking for secondary vertices [21,22] and for disappearing tracks [23].
The current study addresses many different models of new physics, in particular those that predict the production of metastable massive particles with O(ns) lifetime at LHC energies, such as mini-split SUSY [11,24,25] or anomaly-mediated supersymmetry-breaking (AMSB) models [26,27].A metastable gluino with a mass of approximately 1 TeV would be compatible with the measured Higgs boson mass according to mini-split SUSY models, which also predict squarks with masses of 10 3 -10 5 TeV, therefore making the gluino the only observable sparticle produced through the strong interaction at LHC energies.Results are presented assuming the production of R-hadrons as composite colorless states of a gluino together with SM quarks or gluons and the subsequent decay, on nanosecond timescales, of the gluino to a stable neutralino and a q q pair [28].This search could, in principle, also be sensitive to the production of LLPs through weak interactions, but these models are not considered here as the cross-sections are much lower.
In the rest of this paper, R-hadrons with lifetimes consistent with decaying inside the ATLAS detector are referred to as metastable; the search strategy focusses on decays occurring within the active detector volume, which covers lifetimes from around 1 ns to several tens of ns.R-hadrons whose lifetime disfavors decay within the ATLAS detector are termed stable.The approach described in this paper allows direct observation of charged R-hadrons which are either stable or metastable and traverse the first seven layers of the ATLAS inner tracking system.

ATLAS Detector and Pixel dE/dx measurement
The ATLAS detector1 consists of a tracker, for measuring the trajectories of charged particles, surrounded by a solenoid magnet, followed by calorimeters for measuring the energy of particles that interact electromagnetically or hadronically.A muon spectrometer immersed in a toroidal magnetic field surrounds the calorimeters, and provides tracking for muons.The detector is hermetic and can therefore measure the missing transverse momentum associated with each event.A complete description of the ATLAS detector can be found in Ref. [15].The tracker is made of three detector systems.Moving from the solenoid magnet toward the beam collision region, there is first a transition radiation tracker with approximately 400,000 channels [29], followed by a silicon microstrip detector (SCT) with about 6 million channels [30], and at the innermost radius sits a roughly 92 million channel silicon Pixel detector, which is crucial for this measurement, and is described below in some detail.
The current Pixel detector provides at least four precision measurements for each track in the region |η| < 2.5 at radial distances of 3.4 to 13 cm from the LHC beam line.After Run 1, the ATLAS Pixel detector was upgraded with the insertion of an additional layer, the Insertable B-Layer (IBL), which was installed inside the Run 1 Pixel detector, mounted on a new beam pipe of smaller diameter.The IBL has smaller area pixels, reduced thickness, and different coding electronics, which provides charge measurements with lower resolution and dynamic range than the other Pixel layers.It therefore requires a separate treatment of the dE/dx information.At normal incidence, the average charge released by a minimum ionizing particle (MIP) in a pixel sensor is ≈ 20000 e − (≈ 16000 e − for IBL) and the charge threshold is set to 3500±40 e − (2500±40 e − for IBL).Signals are accepted if they are larger than this threshold.They are then matched to a specific beam crossing.The hit efficiency under these conditions exceeds 99%.When detector data are read out, the Time over Threshold (ToT), i.e. the time interval with the signal above the threshold, is digitized with 8 bits (4 bits for IBL).The ToT is approximately proportional to the ionization charge [31] and its dynamic range corresponds to 8.5 times (1.5 times for IBL) the average charge released by a MIP for a track normal to the silicon detectors which deposits all its ionization charge in a single pixel.If this value is exceeded, the hit charge information is either underestimated in the IBL, where the electronics signals the excess with an overflow bit, or lost in the other pixel layers.
The charge released by a track crossing the Pixel detector is rarely contained within just one pixel; neighboring pixels registering hits are joined together using a connected component analysis [32,33] to form clusters.The charge of a cluster is calculated by summing the charges of all pixels belonging to the cluster after calibration corrections; the IBL overflow information is accounted for separately.The dE/dx is estimated using the average of the individual cluster ionization measurements (charge collected in the cluster per unit track length in the sensor), for the clusters associated with a track.To reduce the effect of the tails of the Landau distribution, the average is evaluated after removing the highest dE/dx cluster, or the two highest dE/dx clusters in the rare case of more than four clusters on the track.A track is considered to have a good ionization measurement only if there are at least two remaining clusters.Including the IBL charge information reduces the tails of the dE/dx distribution, as illustrated in Figure 1, and increases the fraction of tracks for which the ionization measurement is valid from 77% to 91%.The most probable value of the measured dE/dx is 1.12 MeV g −1 cm 2 for a minimum ionizing particle.This value Figure 1: Distribution of the dE/dx computed with the Pixel algorithm in data (black dots) and Pixel without IBL (red triangles) for minimum bias event tracks selected as in Ref. [34], with the additional requirements of p > 3 GeV and at least two associated pixel clusters.The corresponding dE/dx distributions for simulation with and without the IBL are also plotted, respectively, as dashed and dotted lines.The simulation overestimates the dE/dx, which is taken into account in the systematic uncertainty evaluation.The addition of IBL clusters reduces the tails at both low and high dE/dx.The fraction of tracks with dE/dx <0.5 (>1.8)MeV g −1 cm 2 decreases with the use of the IBL charge information from 0.065% to 0.018% (from 0.89% to 0.44%).
is obtained with a gaussian fit to the data using the complete Pixel detector (see Figure 1); the variance of the gaussian fit is 0.13 MeV g −1 cm 2 .The βγ measurable with the current dE/dx method cannot be smaller than approximately 0.3 for particles with unit charge because of the ToT dynamic range.

Mass Calculation
The average energy loss of massive, charged particles in matter is expected to follow the Bethe-Bloch distribution, which can be expressed as a function of the βγ of the LLPs.Their mass can be derived from a fit of the measured specific energy loss and the reconstructed momentum to a parametric Bethe-Bloch distribution in the range 0.3 < βγ < 1.5.This range overlaps the expected average βγ of LLPs produced at the LHC, which decreases with the particle mass from βγ ≈ 2.0 at 100 GeV to βγ ≈ 0.5 at 1600 GeV.The parametric function describing the relationship between the most probable value of the energy loss (dE/dx) MPV and βγ is: The p i , with i = 1 . . .5, calibration constants were measured in this data set using low-momentum pions, kaons and protons reconstructed in ATLAS as described in Ref. [35].
The distribution of the dE/dx versus momentum is shown in Figure 2 together with the fits to the pion, kaon and proton masses.The mass calculation performed through the dE/dx method is monitored by checking the stability of the proton mass measurement during the data taking.2: Two-dimensional distribution of dE/dx versus charge signed momentum (qp) for minimum-bias event tracks compatible with a primary vertex, which has the same data selection as in Ref. [34].The distributions of the most probable value for the fitted probability density functions of pions, kaons and protons are superimposed.
A mass estimate M is obtained by numerically solving the equation 1 for the unknown M. In simulated R-hadron events, the reconstructed mass is found to reproduce well the generated mass up to masses around 1800 GeV, and a residual 3% correction is applied to the reconstructed mass to improve the agreement.For particles with higher masses (beyond the current sensitivity of this analysis), the reconstructed mass is observed to slightly underestimate the generated mass.The half-width at half maximum of the reconstructed mass distribution increases with the mass value.This is due to the momentum measurement uncertainty dominating the mass resolution for masses greater than ≈200 GeV.

R-hadron Simulation
A number of samples of simulated signal events are used in this analysis to determine the expected LLP signal efficiencies and to estimate uncertainties in the efficiency.
For stable R-hadrons, pair production of gluinos with masses between 800 GeV and 1800 GeV is simulated in Pythia 6.4.27 [36] with the AUET2B [37] set of tuned parameters for the underlying event and the CTEQ6L1 [38] parton distribution function (PDF) set, incorporating dedicated hadronization routines [39] to produce final states containing R-hadrons.
The cross-sections are calculated at next-to-leading order in the strong coupling constant (NLO), including the resummation of soft-gluon emission at next-to-leading-logarithmic accuracy (NLO+NLL) [40][41][42][43][44].The nominal cross-section predictions and the uncertainty are taken from an envelope of crosssection predictions using different PDF sets and factorization and renormalization scales, as described in Ref. [45].
Selected events containing R-hadrons undergo a full detector simulation [46], where interactions of Rhadrons with matter are handled by dedicated Geant4 [47] routines based on the physics model described in Refs.[39,48].This model assumes, for each light valence quark, a cross-section of 12 mb per nucleon and neglects the heavy parton, other than as a reservoir of kinetic energy.This hadronic scattering is described through a purely phase-space-driven approach.The probability for a gluino to form a gluongluino bound state is assumed to be 10% [5].
The simulation of samples of metastable gluino-based R-hadrons is performed similarly to the stable R-hadrons.The gluinos within R-hadrons are then required to decay via the process g → qqχ 0 1 , using Pythia 6. Gluino masses (m) between 400 GeV and 3000 GeV are simulated, with the neutralino mass fixed to 100 GeV.This benchmark decay mode is chosen as it has the highest sensitivity among those studied in Ref. [16].Heavier neutralino masses would give lower signal efficiencies, mainly because of the smaller missing transverse momentum available to trigger.The gluino lifetime (τ) is varied from 0.4 ns to 50 ns.
All simulated samples include a modeling of pileup, adding the expected number of minimum bias proton-proton interactions from the same (in-time pileup) and nearby (out-of-time pileup) bunch crossings.Simulated events are reconstructed using the software used for the collision data.
Additional samples of gluinos are generated for all mass points using MG5_aMC@NLO [49], to get a more accurate description of radiative effects than Pythia 6 provides.The distribution of the transverse momentum of the gluino-gluino system simulated with Pythia 6 is then reweighted in such a way as to match that obtained in samples simulated with MG5_aMC@NLO.This is especially important for stable R-hadrons as detailed in Section 6.

Search Strategy
This search is based on the signature of R-hadron events containing at least one reconstructed track with large measured dE/dx and high transverse momentum (p T ).This signature optimizes the search for metastable R-hadrons with lifetimes of O(ns), as the selection requires that the R-hadron only live long enough to traverse the seven layers of the silicon detectors, corresponding to a distance from the beam axis of 37 cm.The search is also sensitive to R-hadrons with longer lifetimes, with a small dependence on the modeling of the interaction of R-hadrons in the calorimeter and muon systems.
The expected production and interaction of R-hadrons in the detector is described in detail in Ref. [50].In the models studied in this search and described in Section 4, pair-produced gluinos hadronize into two colorless R-hadrons with the following charge states: approximately 33% events with two neutral R-hadrons, 47% events with one neutral and one charged R-hadron, and 20% events with two charged R-hadrons.The fragmentation is extremely hard and the R-hadrons are therefore isolated.As massive particles with β < 1, the charged R-hadrons are expected to deposit more ionization energy in the detector than typical particles of the same momentum.This ionization signature is used both to select signal events, as described in Section 6, and to calculate the mass of selected candidates, as described in Section 3. The estimated candidate mass is used as the final search discriminant.
If the R-hadrons decay before or in the calorimeter, the decay products of R-hadrons in models studied here include quarks, which are identified as jets in the calorimeters (see Section 6), and a 100 GeV neutralino, which carries away significant momentum undetected.If R-hadrons traverse the calorimeter without decaying, they are expected to deposit significantly less energy in the calorimeters than SM hadrons of the same momentum, as the heavy parton carries almost all the momentum of the composite particle and limited kinetic energy is available for interactions between the light SM partons and the detector.In the models studied here, over 90% of stable R-hadrons, independent of mass, deposit less than 20 GeV of energy in the calorimeters.However, the electric charge of the R-hadron can change due to the hadronic interactions, even if little energy is lost, and the charge of the R-hadron leaving the calorimeter is expected to be decoupled from the charge at production [50].As hadronic interactions are rare in the tracking detectors, the probability for an R-hadron to change its electric charge before the calorimeters is only around 3%.
Events are selected by a trigger requiring missing transverse momentum (E miss T ), calculated from energy deposits in the calorimeter [51].For events with R-hadrons that decay before or inside the calorimeter, undetected neutralinos contribute to the measured transverse momentum imbalance.As stable R-hadrons deposit little energy in the calorimeters, most stable R-hadron events are selected by the trigger only when QCD initial-state radiation (ISR) boosts the R-hadron pair system.Details of the event selection are given in Section 6.

Data Sample and Event Selection
The data sample used for this search was collected while tracking detectors, calorimeters, muon chambers, and magnets were operating normally and corresponds to an integrated luminosity of 3.2 fb −1 in 2015.The luminosity measurement was calibrated during dedicated beam-separation scans, using the same methodology as that described in Ref. [52].The uncertainty in the luminosity measurement is 5%.Candidate events are selected by an E miss T trigger, formed from energy deposits in the calorimeter [51], with a threshold of 70 GeV.The event selection, described in detail below, was optimized relative to the √ s = 8 TeV search [16] to improve the sensitivity to events with R-hadrons relative to background events from SM processes.Major improvements include the addition of selections designed to reject tracks from hadrons and electrons, as well as changes to the isolation requirements, which now reject background events from overlapping tracks by identifying clusters in the Pixel or SCT detectors which are consistent with energy deposition from multiple particles.Two signal regions are defined: one targets metastable R-hadrons with shorter lifetimes, and the other is optimized for metastable or stable R-hadrons that are expected to pass through the muon spectrometer before decaying.A primary vertex reconstructed from at least two well-reconstructed charged-particle tracks, each with p T > 400 MeV, is required in order to remove noncollision background events.If there is more than one such vertex in an event, the primary vertex is defined as the vertex with the largest scalar sum of associated track momentum.To reduce backgrounds from SM processes, events must have E miss T > 130 GeV, where E miss T is the offline missing transverse momentum.Jets are reconstructed with the anti-k t algorithm [53] with radius parameter R = 0.4 using clusters of energy depositions in the calorimeter as inputs.Events are rejected if they contain a jet with E T > 20 GeV that is consistent with noise contributions as determined based on timing and shower shape information.The offline E miss T is estimated from reconstructed electrons with |η| < 2.47, reconstructed muons with |η| < 2.5, calibrated jets with p T > 20 GeV within |η| < 4.5, and reconstructed tracks which are not associated with other objects and that pass kinematic requirements designed to reject poorly reconstructed tracks.
For signal events, the trigger efficiency varies significantly depending on the model considered.For metastable R-hadrons with a mass of 1000 GeV, decaying into a 100 GeV neutralino and quarks, the trigger efficiency increases from 65% for R-hadrons with a lifetime of 50 ns to 95% for R-hadrons with a lifetime of 0.4 ns.For an R-hadron with a lifetime of 10 ns, the trigger efficiency increases slightly with the R-hadron mass, from 91% for a mass of 1000 GeV to 96% for a mass of 1600 GeV.For stable R-hadrons, the trigger efficiency is approximately 40% for all masses studied.All efficiencies quoted in this section include both the acceptance and the selection efficiency effects.
The offline E miss T calculation includes the momentum of identified muons, which is important to reject background events from SM processes with muons, such as Z → µµ, that pass the calorimetric E miss T trigger.For events with metastable R-hadrons, the offline E miss T requirement is very efficient (around 95%) due to the neutralinos in the event.For stable R-hadrons that pass the trigger requirement, the offline E miss T cut is also efficient (around 90%), as the R-hadrons do not deposit significant energy in the calorimeter and in most events the R-hadron pair is boosted from initial-state radiation.Many R-hadrons fail to be identified as muons because they are in a neutral charge state going through either the inner detector or the muon spectrometer.Even if a charged R-hadron with β < 1 traverses the muon spectrometer, it can fail to be identified by the standard muon reconstruction algorithm as a prompt muon due to its late time of arrival.In a typical R-hadron model considered in this search, only 30% of stable R-hadrons with a reconstructed track in the inner detector are identified as prompt muons entering into the E miss T calculation.
Selected events must contain at least one candidate R-hadron track.For a track to be identified as a candidate, it must have p T > 50 GeV and meet the following selection criteria.The track must contain at least seven clusters from silicon detector layers in order to make a good measurement of the track's momentum, must have a cluster in the innermost Pixel layer if expected, and must be associated with the primary vertex.In order to ensure a good ionization measurement, the track must have at least two clusters in the Pixel detector used to measure dE/dx.These selection requirements on associated clusters are summarized as cluster requirements.
To reject background events with overlapping tracks that could produce clusters with significant measured ionization, the candidate track must not contain any cluster that is compatible with contributions from two or more tracks.This requirement is over 90% efficient for signal events with lifetimes of 3 ns or greater.For background tracks, it lessens the high ionization tail of the distribution from photon conversions and collimated particles, reducing the size of the background and leaving the remaining tracks with a dE/dx distribution largely independent of background source.In order to reduce backgrounds from particles inside hadronic jets, the scalar sum of the p T ( p T ) of other tracks, with p T > 1 GeV and consistent with the primary vertex, in a cone of size ∆R < 0.25, where ∆R ≡ (∆η) 2 + (∆φ) 2 , around the candidate track, must be less than 20 GeV.
The background from SM processes is further reduced by requiring that the candidate track momentum exceeds 150 GeV and that the relative uncertainty on the momentum measurement be less than 50%.The contribution from leptons from W boson decays is reduced by imposing a condition that the transverse mass, m T , built with the candidate track and the E miss T , be greater than 130 GeV, where: and ∆φ(E miss T , track) is the azimuthal separation between the track and the E miss T .This selection has an efficiency between 78% and 89% for R-hadrons with lifetimes of 3 ns or greater, and reduces the background by roughly a factor of three.
Background events due to tracks from electrons and hadrons are reduced by requiring that the energy of any jet associated with the candidate track must be less than the candidate track's measured momentum.Jets are associated with tracks if they are within ∆R < 0.05 of the candidate track and have p T > 20 GeV. Background events containing electrons are further reduced by rejecting candidate tracks associated with a jet with 95% or more of its energy deposited in the electromagnetic calorimeter.These selections are found to be 95% efficient for all signal models studied.
Tracks identified as well-reconstructed muons with p T > 25 GeV are rejected in the search for metastable R-hadrons with lifetimes of 30 ns or less.R-hadrons with a lifetime equal to or greater than 50 ns typically reach the muon spectrometer, and therefore no such veto is applied in the long-lived R-hadron search.To further reject hadronic background in the search for stable R-hadrons, the isolation requirement on the p T of nearby tracks is tightened from 20 GeV to 5 GeV for candidate R-hadron tracks, eliminating 30% of the remaining expected background.
The specific ionization of the candidate track, as measured by the Pixel detector, must be larger than 1.80 − 0.11|η| + 0.17η 2 − 0.05|η| 3 MeV g −1 cm −2 .This requirement corrects the slight (5%) |η| dependence of the dE/dx measurement [54] while maximizing the expected sensitivity to signal events.This selection is approximately 70-90% efficient for signal events with lifetimes of 3 ns or greater, depending on the R-hadron mass, and retains less than 1% of background events.Additionally, a time-dependent correction is applied to the ionization measurement in data to account for an observed variation in the charge measurement by the IBL electronics due to effects of radiation.The size of this correction is around 5%.
The observed numbers of events in data and the predicted numbers of simulated signal events, for a 1600 GeV R-hadron with a lifetime of 10 ns, are shown in Table 1 after each stage of the event selection.The total selection efficiency relative to generated events, which includes both the acceptance and efficiency effects, is also shown for the signal.

Selection level
Exp. signal events Observed events in 3.  1: Expected number of events at different steps of the selection for metastable 1600 GeV gluino R-hadrons decaying with a 10 ns lifetime, along with the number of events observed in data, for 3.2 fb −1 .The simulated yields are shown with statistical uncertainties only.The total efficiency × acceptance is also shown in parentheses for the signal.

Signal Selection Efficiency
The overall signal selection efficiency depends on the mass and lifetime of the signal, ranging from about 19% for a 1600 GeV R-hadron with a lifetime of 10 ns to less than 1% for R-hadrons with a lifetime of 0.4 ns.The trigger and offline E miss T requirements are more efficient for R-hadrons that decay before or inside the calorimeters.However, as the R-hadron lifetime decreases, the probability of reconstructing the track from the original R-hadron in the silicon detectors decreases.For these two reasons, R-hadrons with lifetimes in the 10-30 ns range have the highest selection efficiency.
The inefficiency for selecting R-hadrons with short lifetimes can be better understood by calculating the selection efficiency relative to a simple fiducial region that separates out acceptance effects due to basic kinematic properties of the R-hadrons.An event is considered to be within the fiducial acceptance if it contains at least one R-hadron that passes the following requirements at MC generator level: charged at production, p T > 50 GeV, p > 150 GeV, |η| < 2.5, and transverse decay distance ≥ 37 cm from the origin.With this definition, the acceptance for R-hadron events with a lifetime of 1 ns ranges from 7% for a mass of 600 GeV to 4% for a mass of 1600 GeV.The selection efficiency for R-hadron events in the fiducial region has only a mild dependence on the mass; in the R-hadron mass range from 600 to 1600 GeV, the selection efficiences are 20% to 25% (τ = 1 ns), 35% to 39% (τ = 3 ns), 39% to 45% (τ = 10 ns), and 31% to 35% (τ = 30 ns), respectively.Table 2 summarizes the systematic uncertainties in the predicted signal yields.The systematic uncertainties are calculated and used independently for each mass bin and signal lifetime.In the table, only the maximum uncertainty is quoted.
The trigger and offline E miss T selections are sensitive to the modeling of ISR in the event, in particular for signal events that do not decay inside the detector volume.Half of the difference between the selection efficiency of the Pythia 6 samples and those reweighted with MG5_aMC@NLO is taken as an uncertainty due to radiative effects.
The E miss T trigger modeling is studied in a sample of Z → µµ events and an uncertainty is assigned based on differences between efficiencies measured in data and in simulation.The uncertainty in the E miss T calibration comes primarily from uncertainties on the jet energy scale.A systematic uncertainty in the measurement of the Pixel ionization is evaluated by comparing the ionization of an inclusive sample of simulated and observed low-momentum tracks.An asymmetric uncertainty is assigned due to the selection efficiency difference observed after the most probable value of the signal distribution is reweighted based on the comparison of low-momentum tracks between data and simulation.The uncertainty in the signal acceptance is asymmetric since the ionization in simulation is systematically larger than that observed in data.An uncertainty in the efficiency of the muon veto is estimated by studying the modeling of the muon spectrometer's timing measurement.
The signal selection efficiency is found to be largely independent of the number of pileup events in simulation.However, to cover any remaining difference, efficiency differences from variations in the pileup distributions are accounted for as an additional source of uncertainty.The uncertainty on the momentum scale and resolution of the tracks results in a small uncertainty in the selection efficiency.Variations of the requirements to reject tracks from electrons and hadrons are found to have a negligible effect on the signal efficiency.The theoretical uncertainty in the signal cross-section is described in Section 4.

Background Estimation
The expected shape and normalization of the mass distribution of background events from SM processes is derived from data.The distributions of key variables are extracted in two control regions in data and these templates are used to generate the expected background distribution in the signal region.The choice of control samples takes into account the measured correlations between the key variables p, dE/dx and η, and the control region selections minimize the possible signal contamination.The same regions are used in the searches for stable and metastable particles, except for the rejection of track candidates identified as muons in the latter case and for the tighter track isolation requirement in the former case.
The first control region (hereafter called CR1) is selected by applying all the selections described in Section 6, except for the requirement on high ionization, which is inverted.The tracks in this control region are kinematically similar to those in the signal region for SM processes, and the expected p and η(p) distributions of tracks from background processes are extracted from this region.The second sample (CR2 hereafter) is selected by inverting the E miss T requirement (E miss T < 130 GeV), while keeping all other selections unchanged.As more than 90% of R-hadron events that pass the trigger have E miss T > 130 GeV, while the E miss T distribution from SM background processes falls steeply in the range 70-130 GeV, CR2 is background-dominated. Tracks from this region are used to derive dE/dx templates in bins of η, as E miss T and dE/dx are uncorrelated for tracks from SM processes.The signal contamination is less than 1% in both control regions for all R-hadron masses and lifetimes considered in this search.
The shape of the background in the signal region is estimated from the control region templates with the following procedure: 1 million (p, η, dE/dx) triplets are generated, where the momentum is randomly sampled according to the template from CR1 tracks, the pseudorapidity is generated according to the η(p)-binned templates based on CR1 tracks, and the ionization is sampled according to dE/dx(η)-binned templates from CR2 tracks.This procedure maintains the measured correlations between p, η, and dE/dx for tracks from SM sources.The particle mass M is calculated given the predicted dE/dx and p values, using the technique explained in Section 3. Table 3: Summary table of the uncertainties in the background estimate.The uncertainties are evaluated integrating over the mass windows considered for the signal search (see Section 9).If the uncertainty depends on the mass, the largest value is quoted.

Source of uncertainty
The normalization of the generated background is obtained by scaling the background to the data in the shoulder region of the mass distribution (i.e.M < 160 GeV), where a possible signal has already been excluded [16].The normalization is performed on the samples before the ionization requirement.The expected mass distributions of the background and the observed data, before the ionization requirement, are shown in Figure 3 for the metastable R-hadrons selection.
Entries / 50 GeV  Systematic uncertainties in the background estimate are evaluated by varying the relative fraction of muons in the control regions, using an analytical description to vary the shape of the high ionization tail in the dE/dx template, and by changing the IBL ionization correction by ±1σ.The statistical uncertainty in the template shape, which is estimated by Poisson fluctuations of the templates, dominates.There is also a small statistical uncertainty from the size of the normalization region.The breakdown of uncertainties in the background estimate is shown in Table 3.
The complete procedure is tested on signal-depleted validation regions.The validation regions are selected with the same requirements as for the nominal signal and control regions, but requiring the track momentum to be in the range 50 < p < 150 GeV.Applied to these validation samples, the procedure described above yields a calculated background of 20.0 ± 3.0 (stat.)± 1.2 (syst.)events in the validation region for metastable R-hadrons, while 14 events are observed in the data.For the control region for stable R-hadrons, 28.1 ± 4.2 (stat.)± 1.9 (syst.)events are estimated from this procedure and 20 are observed.While the observed numbers of events in the validation regions are lower than predicted, they are consistent with statistical fluctuations at the level of 1.7σ.

Results
In data, 16 events are observed for the stable R-hadron selection and 11 are observed for the metastable selection.Table 4 summarizes the background estimates with total statistical and systematic uncertainty as well as the observed events for the metastable and stable R-hadron selection.The observed mass distribution in data is shown in Figure 4 for both selection regions, along with the background expectation.No evidence of a signal above the background is observed.Expected and observed upper limits on Rhadron production cross-sections are evaluated at 95% CL for a discrete set of R-hadron mass values by counting simulated signal, background and observed data events that survive the selection in a ±1.4 σ wide mass window around the position of each simulated signal mass peak for the corresponding mass hypothesis being probed.The peak position and width are estimated by a Gaussian fit to each individual signal mass distribution (σ ≈280 [≈500] GeV at an R-hadron mass of 1200 [1600] GeV).When the lifetime of the R-hadron is long enough to reach the muon spectrometer, the stable selection becomes more efficient than the metastable one and improves the expected sensitivity.This transition happens around 50 ns.The choice of whether to apply the metastable or stable R-hadron selection to each lifetime is based on the best expected limit.Based on this procedure, the stable selection is applied to both the 50 ns and stable lifetimes.
The cross-section upper limits are extracted, separately for the stable and the metastable searches (in the latter case, for each lifetime value), with the CL S method [55], using the profile likelihood ratio as a test statistic.In the procedure, the systematic uncertainties in the signal and background yields, as evaluated in Section 7 and Section 8, are treated as Gaussian-distributed nuisance parameters.The statistical uncertainty in the background distribution also takes into account the uncertainty due to the normalization.Lower limits on the R-hadron mass are then derived, for each lifetime, by comparing the measured cross-section upper limits to the theoretically predicted production cross-section.The resulting lower limits set on the mass of the stable and metastable particles are summarized in Table 5.
The strongest mass limits are obtained for lifetimes of 10 ns or more.The cross-section upper limits for 10 ns lifetime R-hadrons decaying to q q plus a light neutralino of mass m( χ0 )=100 GeV are shown in Figure 5.For this signal, the expected and observed lower mass limits are 1600 GeV and 1570 GeV, respectively.The 8 TeV result excluded R-hadrons with a 10 ns lifetime and masses up to 1185 GeV, using the lower edge of the ±1σ band around the theoretical cross-section.Table 5: Observed and expected 95% CL lower limit on mass is listed for gluino R-hadrons that are either stable or decaying to q q plus a light neutralino of mass m( χ0 ) = 100 GeV.The stable selection is more efficient for lifetimes of 50 ns or more.

Mass [GeV]
The excluded regions on the lifetime-mass plane for R-hadrons decaying to q q plus a light neutralino of mass m( χ0 )=100 GeV are shown in Figure 6.

ATLAS
Figure 5: cross-section as a function of mass for gluino R-hadrons with lifetime τ=10 ns, decaying to q q plus a light neutralino of mass m( χ0 ) = 100 GeV.Theoretical values for the cross-section are shown with their error.The expected upper limits (UL) in the case of background only is shown by the solid black line, with its ±1σ and ±2σ bands, green and yellow respectively.The observed 95% UL is shown as a solid red line.The cross-section for the same process and the corresponding 95% UL measured at √ s = 8 TeV [16] are shown as dashed lines.

Conclusion
This paper presents the results of a search for massive, long-lived particles with lifetimes of O(ns) or longer, produced in 3.2 fb −1 of pp collision data at √ s = 13 TeV at the LHC, and identified by anomalous specific ionization energy loss in the ATLAS Pixel detector.This analysis takes advantage of the addition of a fourth silicon layer to the Pixel detector, the increased parton luminosity, and analysis improvements  to significantly extend a search performed at √ s = 8 TeV [16].No excess of events is observed over the background estimate.Gluino R-hadrons with lifetimes above 0.4 ns and decaying to q q plus a 100 GeV neutralino are excluded at the 95% confidence level with lower mass limit range between 740 GeV and 1590 GeV.In the case of stable R-hadrons the lower mass limit at the 95% confidence level is 1570 GeV.The observed lower limit on R-hadron masses increases by up to approximately 400 GeV relative to the equivalent analysis at √ s = 8 TeV.The crucial computing support from all WLCG partners is acknowledged gratefully, in particular from CERN and the ATLAS Tier-1 facilities at TRIUMF (Canada), NDGF (Denmark, Norway, Sweden), CC-IN2P3 (France), KIT/GridKA (Germany), INFN-CNAF (Italy), NL-T1 (Netherlands), PIC (Spain), ASGC (Taiwan), RAL (UK) and BNL (USA) and in the Tier-2 facilities worldwide.
[16] ATLAS Collaboration, Search for metastable heavy charged particles with large ionisation energy loss in pp collisions at √ s = 8 TeV using the ATLAS experiment, Eur.Phys.J. C 75 (2015) 407.

Figure
Figure2: Two-dimensional distribution of dE/dx versus charge signed momentum (qp) for minimum-bias event tracks compatible with a primary vertex, which has the same data selection as in Ref.[34].The distributions of the most probable value for the fitted probability density functions of pions, kaons and protons are superimposed.

Figure 3 :
Figure3: Expected mass distribution from background events in the signal region for the metastable R-hadrons, before the ionization requirement is imposed.Data are superimposed.The yellow band around the background estimation includes both the statistical and systematic uncertainties.The ratio of data to expected background is shown in the lower panel with statistical uncertainties only.

Figure 4 :
Figure 4: Mass distribution for data and background for stable (top) and metastable (bottom) particle searches.The yellow band around the background estimation includes both the statistical and systematic uncertainties.Also shown are two examples for signals as expected for gluino R-hadrons in the explored mass range.
e V cr o ss -s e c. p re d ic tio n 8 T e V c ro s s -s e c .p re d ic ti o n 8 TeV cross-sec.95% CL UL

Figure 6 :
Figure 6: Excluded range of lifetimes as a function of gluino R-hadron mass.The expected lower limit (LL), with its experimental ±1σ band, is given with respect to the nominal theoretical cross-section.The observed 95% LL obtained at √ s = 8 TeV [16] is also shown for comparison.

[ 18 ]= 8 = 8
ATLAS Collaboration, Search for charginos nearly mass degenerate with the lightest neutralino based on a disappearing-track signature in pp collisions at √ s TeV with the ATLAS detector, Phys.Rev. D 88 (2013) 112006.[19] ATLAS Collaboration, Search for massive, long-lived particles using multitrack displaced vertices or displaced lepton pairs in pp collisions at √ s TeV with the ATLAS detector, Phys.Rev. D 92 (2015) 072004.[20] CMS Collaboration, Search for heavy long-lived charged particles in pp collisions at √ s = 7 TeV, Phys.Lett.B 713 (2012) 408.[21] CMS Collaboration, Search for long-lived neutral particles decaying to quark-antiquark pairs in proton-proton collisions at √ s = 8 TeV, Phys.Rev. D 91 (2015) 012007.[22] CMS Collaboration, Search for long-lived particles that decay into final states containing two electrons or two muons in proton-proton collisions at √ s = 8 TeV, Phys.Rev. D 91 (2015) 052012.[23] CMS Collaboration, Search for disappearing tracks in proton-proton collisions at √ s = 8 TeV, JHEP 1501 (2015) 096.ak Also at Institute for Nuclear Research and Nuclear Energy (INRNE) of the Bulgarian Academy of Sciences, Sofia, Bulgaria al Also at Faculty of Physics, M.V.Lomonosov Moscow State University, Moscow, Russia am Also at Institute of Physics, Academia Sinica, Taipei, Taiwan an Also at National Research Nuclear University MEPhI, Moscow, Russia ao Also at Department of Physics, Stanford University, Stanford CA, United States of America ap Also at Institute for Particle and Nuclear Physics, Wigner Research Centre for Physics, Budapest, Hungary aq Also at Flensburg University of Applied Sciences, Flensburg, Germany ar Also at University of Malaya, Department of Physics, Kuala Lumpur, Malaysia as Also at CPPM, Aix-Marseille Université and CNRS/IN2P3, Marseille, France * Deceased

Table 2 :
Systematic uncertainties on the predicted signal yields.The uncertainty depends on the mass and on the lifetime; only the maximum negative and positive values across all signal models (see Section 4) are reported here.

Table 4 :
Estimated number of background events and the number of observed events in data in the final selection regions.The background predictions show both the statistical and systematic uncertainties.