1 Introduction

Muon identification and measurement in ATLAS relies on two complementary detectors, one of which is its huge Muon Spectrometer (MS). This is based on the use of three very large air core toroidal magnets, each containing eight superconducting coils, and three measuring planes of high-precision chambers. This system is designed for efficient muon detection even in the presence of very high particle backgrounds and for excellent muon momentum resolution up to very high momenta of \(\sim {}1\) TeV. This unprecedented stand-alone performance of the ATLAS muon spectrometer is due to the large field integral (ranging between 2 and 6 Tm for most of the detector), the very low multiple scattering in the material of the air core toroids (1.3 units of radiation length over a large fraction of the acceptance in the barrel toroid), the very high precision measurements along the muon trajectory (chamber resolution \(35\,\upmu \)m) and the extreme alignment precision of the measuring planes (30 \(\upmu \)m).

The other very important component of the muon identification and measurement in ATLAS is the inner detector (ID). In ATLAS the very efficient muon detection and high momentum resolution, with nominal relative momentum resolutions of \({<}3.5\,\%\) up to transverse momenta \(p_{\mathrm {T}}\,{\sim }\,200\) GeV and \({<}10\,\%\) up to \(p_{\mathrm {T}}\sim 1\) TeV, are obtained by a combination of measurements from the ID and the MS [1, p.162]. The complementarity of these measurements can be exploited to provide measurements of the muon reconstruction efficiencies in both tracking systems. In this paper, the muon reconstruction efficiencies are measured using dimuon decays of \(J/\psi \) mesons to access the region \(p_{\mathrm {T}}\,{<}\,10\) GeV and dimuon decays of \(Z\) bosons to access the region \(20\text { GeV}\,{<}\,p_{\mathrm {T}}\,{<}\,100\text { GeV}\). The efficiency determination in the region \(10\text { GeV}\,{<}\,p_{\mathrm {T}}\,{<}\,20\text { GeV}\) is not possible due to the limited sample of muons with \(p_{\mathrm {T}}\) higher than 10 GeV in the \(J/\psi \) decays and difficulties in controlling the backgrounds in the sample of \(Z\) decays that lead to muons with \(p_{\mathrm {T}}\) smaller than 20 GeV. For these analyses, one of the decay muons is reconstructed in both detector systems and the other is reconstructed by just one of the systems in order to probe the efficiency of the other. This method (known as tag-and-probe, and described in more detail in Sect. 4) is applied to the ATLAS proton–proton (\(pp\)) collision data recorded at the Large Hadron Collider (LHC) in 2010 at a centre-of-mass energy of 7 TeV.

Muon isolation criteria are used to select muons in many physics analyses, and measurements of the isolation efficiency performed using \(Z\rightarrow \mu ^+\mu ^-\) decays are described in Sect. 9. The invariant mass distributions from these data are also used to extract the muon momentum resolutions. The analysed data sample corresponds to the full 2010 \(pp\) dataset with an integrated luminosity of 40 pb\(^{-1}\) [2] after applying beam, detector and data-quality requirements.

2 The ATLAS detector

A detailed description of the ATLAS detector can be found elsewhere [3]. Muons are independently measured in the ID and in the MS.

The ID measures tracks up to \(|\eta |=2.5\) Footnote 1 exploiting three types of detectors operated in an axial magnetic field of \(2\) T: three layers of silicon pixel detectors closest to the interaction point, four layers of semiconductor microstrip detectors (SCT) surrounding the pixel detector, and a transition radiation straw-tube tracker (TRT) covering \(|\eta |<2.0\) as the outermost part. The innermost pixel layer (known as the b-layer) has a radius of \(50.5\) mm in the barrel, whilst the outermost TRT tubes are at \(r\approx 1\) m.

The electromagnetic and hadronic calorimeters surround the ID and cover the pseudorapidity range \(|\eta |<4.9\), far beyond the range over which muons are identified. In the barrel and end-cap, in the region \(|\eta |<3.2\) the electromagnetic calorimeter consists of lead absorbers with liquid-argon (LAr) as active material. The barrel hadronic tile calorimeter is a steel/scintillating-tile detector and is extended by two end-caps with LAr as the active material and copper as absorber. The total combined thickness of 11 interaction lengths (\(\lambda \)) includes 9.7 \(\lambda \) of active calorimeter and 1.3 \(\lambda \) of outer support.

The magnetic field of the MS is produced by three large air-core superconducting toroidal magnet systems (two end-caps, where the average field integral is about 6 Tm, and one barrel, where the field integral is about 2.5 Tm). The field is continuously monitored by approximately 1800 Hall sensors distributed throughout the spectrometer volume. The deflection of the muon trajectory in this magnetic field is measured via hits in three layers of precision monitored drift tube (MDT) chambers for \(|\eta |<2.0\) and two outer layers of MDT chambers in combination with one layer of cathode strip chambers (CSCs) in the innermost end-cap wheels (\(2.0\le |\eta |<2.7\)). Three layers of resistive plate chambers (RPCs) in the barrel (\(|\eta |<1.05\)) and three layers of thin gap chambers (TGCs) in the end-caps (\(1.05<|\eta |<2.4\)) are used by the muon trigger (see below). The RPCs, TGCs and CSCs also measure the muon trajectory in the non-bending (\(\phi \)) plane of the spectrometer magnets. The following text frequently refers to chambers which make a measurement in the bending (\(\eta \)) plane as ‘precision chambers’, since these have a much better spatial resolution (important for a good momentum resolution) than the chambers used for triggering.

The chambers are monitored by an optical alignment system, designed to provide an accuracy of 30 \(\upmu \)m in the barrel and 40 \(\upmu \)m in the end-cap [4]. The absolute alignment combines information from the optical system and a track-based procedure (using cosmic data samples and special runs without the toroidal field), and then the optical alignment is used continuously to monitor relative chamber displacements with time (for example, when the toroid magnet is turned on). With the limited 2010 data sample, the muon spectrometer alignment precision was around three times larger than designed.

The ATLAS detector has a three-level trigger system: level 1 (L1), level 2 (L2), and the event filter (EF). The MS provides a L1 hardware muon trigger which is based on hit coincidences in different RPC and TGC detector layers within programmed geometrical windows which define the muon \(p_{\mathrm {T}}\). The L2 and EF muon triggers perform a software confirmation of the L1 muon trigger using refined \(p_{\mathrm {T}}\) measurements from the precision chambers.

Figure 1 shows a schematic drawing of the ATLAS MS. The barrel muon chambers are installed around the calorimeters in roughly cylindrical rings of approximately 5, 7 and 9 m radius. Large barrel chambers are mounted between the barrel toroid coil cryostats. Small barrel chambers are installed on the toroid coil cryostats. The barrel end-cap extra (BEE) chambers are mounted on the end-cap toroid cryostats. The end-cap chambers are arranged in disks with \(z\) axis positions of approximately 7, 13 and 21 m from the centre of the detector, and which are orthogonal to the proton beams.

Fig. 1
figure 1

Schematic drawing of the ATLAS muon spectrometer

3 Muon reconstruction and identification in ATLAS

Muon identification in ATLAS uses independent track reconstruction in the ID and MS, which are then combined. Track reconstruction in the muon spectrometer is logically subdivided into the following stages: pre-processing of raw data to form drift-circles in the MDTs or clusters in the CSCs and the trigger chambers, pattern-finding and segment-making, segment-combining, and finally track-fitting. Track segments are defined as straight lines in a single MDT or CSC station. The search for segments is seeded by a reconstructed pattern of drift-circles and/or clusters.

Full-fledged track candidates are built from segments, typically starting from the outer and middle stations and extrapolating back through the magnetic field to the segments reconstructed in the inner stations (though other permutations are also explored). Each time a reasonable match is found, the segment is added to the track candidate. The final track-fitting procedure takes into account all relevant effects (e.g. multiple scattering, field inhomogeneities, inter-chamber misalignments, etc.). More details about the muon reconstruction can be found in Ref. [1, p. 165].

A similar approach is followed by the ID track reconstruction where the pattern recognition uses space-points formed from the pixel and SCT clusters to generate track seeds. These seeds are then extended into the TRT and drift circles are associated. Finally the tracks are refitted with the information coming from all three detectors. More details about the ID track reconstruction can be found in Ref. [1, p. 19].

The analyses presented here make use of three classes of reconstructed muons, as described below.

  • Stand-alone (SA) muon: the muon trajectory is reconstructed only in the MS. The direction of flight and the impact parameter of the muon at the interaction point are determined by extrapolating the spectrometer track back to the point of closest approach to the beam line, taking into account the energy loss of the muon in the calorimeters.

  • Combined (CB) muon: track reconstruction is performed independently in the ID and MS, and a combined track is formed from the successful combination of a SA track with an ID track.

  • Segment-tagged (ST) muon: a track in the ID is identified as a muon if the track, extrapolated to the MS, is associated with at least one segment in the precision muon chambers.

The main goals of this paper are the measurement of the reconstruction efficiencies, for combined (CB) and combined-plus-segment-tagged (CB+ST) muons, and reconstruction resolutions, for MS and ID muons. The use of the ID for CB and CB+ST muons limits their acceptance to \(|\eta |<2.5\). Stand-alone muons are employed to measure the muon reconstruction efficiency in the ID.

The CB muon candidates constitute the sample with the highest purity. The efficiency for their reconstruction is strongly affected by acceptance losses in the MS, mainly in the two following regions:

  • at \(\eta \sim 0\), the MS is only partially equipped with muon chambers in order to provide space for services of the ID and the calorimeters;

  • in the region (\(1.1<|\eta |<1.3\)) between the barrel and the end-caps, there are regions in \(\phi \) where only one layer of chambers is traversed by muons in the MS, due to the fact that some chambers were not yet installed in that region during the 2010–2012 data-taking. Here no stand-alone momentum measurement is available and the CB muon reconstruction efficiency is decreased.

The reconstruction algorithms for ST muons have higher efficiency than those for CB muons as they can recover muons which did not cross enough precision chambers to allow an independent momentum measurement in the MS. They are also needed for the reconstruction of low-\(p_\mathrm {T}\) muons which only reach the innermost layer of the muon chambers. Due to their lower purity and poorer momentum resolution, ST muons are only used in cases where no CB muon can be reconstructed.

In the early phase of the LHC operation, ATLAS used two entirely independent strategies for the reconstruction of both the CB and ST muons. These two approaches, known as chain 1 and chain 2 in the following, provide an invaluable cross-check on the performance of a very complex system, and allow ATLAS to ultimately take the best aspects of both. The chains have slightly different operating points, with chain 1 typically more robust against background, whilst chain 2 has a slightly higher efficiency.

In chain 1, the momentum of the muon is obtained from a statistical combination of the parameters of the tracks reconstructed by the ID and MS [1, p. 166]. SA muon tracks are required to have a sufficient number of hits in the precision and trigger chambers, to ensure a reliable momentum measurement. In chain 2, the combined muon momentum is the result of a simultaneous track fit to the hits in the ID and the MS. The requirements applied to the hit multiplicities in the MS are less stringent than in chain 1 because certain information, such as the trajectory in the plane transverse to the proton beams, is better provided by the ID in the simultaneous fit. In both chains, muon track segments can additionally be assigned to ID tracks to form ST muons, based on the compatibility of the segment with the extrapolated ID track.

To illustrate the high purity of the ATLAS muon identification and the size of the dimuon dataset, Fig. 2 shows the reconstructed invariant mass distribution of opposite-sign muon candidate pairs. The events are selected by an unprescaled, 15 GeV \(p_{\mathrm {T}}\) threshold single muon trigger, which is reconfirmed offline by requiring at least one muon to have \(p_{\mathrm {T}}\!>\!15\) GeV. Both muons are required to be of CB type and to pass the ID track selection criteria of Sect. 6.2. The distance of closest approach of the muon to the primary vertex is limited to 5 mm in the transverse plane and 200 mm/\(\sin \theta \) in the longitudinal direction. The \(J/\psi \), \(\Upsilon \) and \(Z\) peaks are clearly visible, and the muon reconstruction has the capability to resolve close-by resonances, such as the \(J/\psi \) and \(\psi ^\prime \) as well as the \(\Upsilon (1S)\) and \(\Upsilon (2S)\). The shoulder near \(m_{\mu \mu }\!\approx \!15\,\)GeV is caused by the kinematic selection.

Fig. 2
figure 2

Reconstructed invariant mass, \(m_{\mu \mu }\), distribution of muon candidate pairs. The number of events is normalised by the bin width. The uncertainties are statistical only

4 The tag-and-probe method

As track reconstruction is performed independently in the ID and MS, the reconstruction efficiency for CB or ST muons is the product of the muon reconstruction efficiency in the ID, the reconstruction efficiency in the MS, and the matching efficiency between the ID and MS measurements (which includes the refit efficiency in the case of chain 2). It is therefore possible to study the full reconstruction efficiency by measuring these individual contributions. A tag-and-probe method is employed, which is sensitive to either the ID efficiency or the combined MS and matching efficiency.Footnote 2 This technique is applied to samples of dimuons from the \(J/\psi \) and \(Z\) decays.

For \(Z\rightarrow \mu ^+\mu ^-\) decays, events are selected by requiring two oppositely charged isolated tracks with a dimuon invariant mass near the mass of the \(Z\) boson. One of the tracks is required to be a CB muon candidate, and to have triggered the readout of the event (see Sect. 6). This muon is called the tag. The other track, the so-called probe, is required to be a SA muon if the ID efficiency is to be measured. If the MS reconstruction and matching efficiency is to be measured the probe must be an ID track. The ID reconstruction efficiency is defined as the fraction of SA probes which can be ascribed to an inner detector track. The combined MS and matching efficiency is the fraction of ID probes which can be associated to a CB or ST muon.

The invariant mass spectra of \(Z\) boson tag-and-probe pairs, shown in Fig. 3, illustrate how muon isolation requirements (see Sects. 6 and 9) almost entirely remove contributions from background processes, resulting in a relatively pure sample of muon tag-and-probe pairs. Monte Carlo studies show that the contribution from other sources is below 0.1 % when MS probes are used and below 0.7 % when ID probes are used. These backgrounds arise from \(Z\rightarrow \tau ^+\tau ^-\), \(W^\pm \!\rightarrow \!\mu ^{\pm \,}{}^{\scriptscriptstyle (\!}_{}\bar{\nu }^{\scriptscriptstyle )}_{{\scriptstyle \mu }}\) \(W^\pm \!\rightarrow \!\tau ^{\pm \,}{}^{\scriptscriptstyle (\!}_{}\bar{\nu }^{\scriptscriptstyle )}_{{\scriptstyle \tau }}\) \(b\bar{b}\), \(c\bar{c}\), and \(t\bar{t}\). The presence of backgrounds in the data leads to an apparent decrease in the muon efficiency in the range \(p_{\mathrm {T}}\lesssim 30\) GeV, for both reconstruction chains. This is taken into account by comparing the measured efficiencies to efficiencies predicted using simulated samples which include these background contributions.

Fig. 3
figure 3

Invariant mass, \(m\), distribution of pairs of tag muons (chain 2) and ID track probes for different sets of muon isolation requirements for the \(Z\) boson analysis, as indicated in the legend

To investigate the reconstruction efficiency at lower transverse momenta, dimuon pairs from \(J/\psi \rightarrow \mu ^+\mu ^-\) decays are used in the same way as those from \(Z\rightarrow \mu ^+\mu ^-\) decays. Because \(J/\psi \) mesons are produced inside jets, isolation requirements cannot be used to select a pure sample. In this case, the invariant mass distribution of the tag-and-probe pairs is fitted using the sum of a quadratic background term and a Gaussian signal term [5]. This is illustrated in Fig. 4 for probe muons selected in the range \(0.1<|\eta |<1.1\) and \(3\mathrm {\,GeV}<p_{\mathrm {T}}<4\) GeV. The invariant mass spectra are shown for tag-and-probe pairs in which the probes are matched to reconstructed muons (see Sect. 6.5) and for unmatched tag-and-probe pairs. The muon reconstruction efficiency is then extracted from a simultaneous fit to the distributions obtained from the matched and unmatched tag-and-probe pairs.

Fig. 4
figure 4

Distribution of the invariant mass, \(m\), of the unmatched (upper distributions) and matched (lower distributions) tag-and-probe pairs for CB and CB+ST muons of chain 2, for the \(J/\psi \) analysis with a probe muon selection as described in the legend. Also shown are the results of the fit using a Gaussian signal and a quadratic background contribution

5 Monte Carlo samples and expectations

The measurements presented in this paper are compared with predictions of Monte Carlo (MC) simulations. For the efficiency measurements in the region \(p_{\mathrm {T}}>20\) GeV, five million \(Z\rightarrow \mu ^+\mu ^-\) events were simulated with PYTHIA 6.4 [6], passed through the full simulation of the ATLAS detector [7], based on GEANT4 [8, 9], and reconstructed with the same reconstruction programs as the experimental data.

During the 2010 data taking, the average number of \(pp\) interactions per bunch crossing was about 1.5. This “pile-up” is modelled by overlaying simulated minimum bias events on the original hard-scattering event. It is found to have a negligible impact for these measurements. The following background samples were used: \(Z\!\rightarrow \!\tau ^+\tau ^-\), \(W^\pm \!\rightarrow \!\mu ^{\pm \,}{}^{\scriptscriptstyle (\!}_{}\bar{\nu }^{\scriptscriptstyle )}_{{\scriptstyle \mu }}\) \(W^\pm \!\rightarrow \!\tau ^{\pm \,}{}^{\scriptscriptstyle (\!}_{}\bar{\nu }^{\scriptscriptstyle )}_{{\scriptstyle \tau }}\) \(b\bar{b}\), \(c\bar{c}\), and \(t\bar{t}\). More details can be found in Ref. [10].

The reconstruction efficiency at low \(p_{\mathrm {T}}\) was studied with a simulated sample of five million prompt \(J/\psi \) events generated with PYTHIA using the PYTHIA implementation of the colour-octet model. In order to increase the number of events at the higher end of the low-\(p_{\mathrm {T}}\) region, this sample was supplemented with a sample of one million \(pp\!\rightarrow \!b\bar{b}\) events also generated with PYTHIA, in which at least one \(J/\psi \) decaying into muons of \(p_{\mathrm {T}}> 2.5\) GeV was required in the \(b\)-quark decay chain.

The reconstruction efficiencies obtained from the analysis of the \(J/\psi \) Monte Carlo samples are shown in Fig. 5, as a function of \(p_{\mathrm {T}}\) and \(\eta \), for CB and CB+ST muons from chain 1. The most discernible features are the areas of lower efficiency at fixed \(\eta \) that result from the un-instrumented (‘crack’) region in the MS at \(\eta \sim 0\) and from the barrel/end-cap transition regions where the chamber configuration (\(1.1<|\eta |<1.3\)) and the magnetic field (\(1.1<|\eta |<1.7\)) are rather non-uniform. Also visible is the impact of the energy loss in the calorimeter on the efficiency, for muons with \(p_{\mathrm {T}}\) of less than 2–5 GeV (depending on the \(\eta \) region), which are absorbed in the calorimeter. For \(|\eta |<2.0\), the CB+ST muon reconstruction starts to be efficient at \(p_{\mathrm {T}}\) values lower than in the reconstruction of pure CB muons, since it includes muons reaching only the inner layer of MDT chambers. For \(|\eta |>2.0\) the CB and CB+ST efficiencies are very similar for chain 1, because cases with only one segment in the CSC chambers, corresponding to the inner layer of precision chambers in this region, are not considered for ST muons. Chain 2 does make use of these segments, and shows an improved CB+ST efficiency in this region (see Sect. 8.3). These detector features motivate the binning used for the determination of the \(p_{\mathrm {T}}\) dependence of the reconstruction efficiency at low \(p_{\mathrm {T}}\).

Fig. 5
figure 5

The chain 1 muon reconstruction efficiency from simulated \(J/\psi \) decays for CB (left) and CB+ST (right) muons as a function of \(\eta \) and \(p_{\mathrm {T}}\) for efficiency values above 50 %

For the \(J/\psi \!\rightarrow \!\mu ^+\mu ^-\) analysis the measured efficiencies are separated into five pseudorapidity intervals according to the different MS regions:

$$\begin{aligned} \left. \begin{array}{l@{\quad }l} |\eta |<0.1 &{} \text {the } \eta =0 \text { crack region};\\ 0.1<|\eta |<1.1 &{} \text {the barrel region};\\ 1.1<|\eta |<1.3 &{} \text {the transition region} \\ &{} \text {between barrel and end-cap}; \\ 1.3<|\eta |<2.0 &{} \text {the end-cap region};\\ 2.0<|\eta |<2.5 &{} \text {the forward region}.\\ \end{array} \right. \end{aligned}$$
(1)

Muons from \(Z\!\rightarrow \!\mu ^+\mu ^-\) decays were required to have \(p_{\mathrm {T}}>20\) GeV. In contrast to the case of lower-\(p_{\mathrm {T}}\) muons from \(J/\psi \) decays, the \(\phi \) deflections of these muons by the magnetic fields in the detector are so small that one can use the muon directions of flight at the \(pp\) interaction point to associate them with specific (\(\eta ,\phi \)) regions of the MS. Ten different regions are defined, corresponding to ten different physical regions in the MS [3]. In each of these, the muon traverses a particular set of detector layers and encounters a different quality of detector alignment, a different amount of material or a different magnetic field configuration. The ten regions are described below (see also Fig. 1).

  • Barrel large: the regions containing large barrel chambers only, which are mounted between the barrel toroid coils.

  • Barrel small: the regions containing small barrel chambers only, which are mounted on the barrel toroid coils.

  • Barrel overlap: the regions where small and large barrel chambers have slight overlaps in acceptance.

  • Feet: the detector is supported by ‘feet’ on its bottom half, which results in a loss of acceptance due to missing chambers, making muon reconstruction more challenging.

  • Transition: the region \(1.1<|\eta |<1.3\), between the barrel and the end-cap wheels.

  • End-cap small: the small end-cap sectors, consisting of MDT chambers.

  • End-cap large: the large end-cap sectors, consisting of MDT chambers and which (in contrast to the Barrel large regions) contain the toroid coils.

  • BEE: the regions containing barrel end-cap extra chambers, which are mounted on the end-cap toroid cryostats.

  • CSC small: the end-cap sectors consisting of small CSC chambers.

  • CSC large: the end-cap sectors consisting of large CSC chambers.

6 Selection of tag-and-probe pairs

6.1 Event selection

The events used for the efficiency measurements were selected online with a single-muon trigger. For the studies with \(J/\psi \!\rightarrow \!\mu ^+\mu ^-\) decays, a combined muon is required, with minimum \(p_{\mathrm {T}}\) thresholds of 4, 6, 10, or 13 GeV (as it was necessary to increase the thresholds during the year, in order to keep the trigger rate within limits). For the studies with \(Z\!\rightarrow \!\mu ^+\mu ^-\) decays, events have to pass the lowest \(p_{\mathrm {T}}\) threshold muon trigger that was unprescaled. The thresholds of the selected triggers range from 10 to 13 GeV, well below the transverse momentum threshold of the tag muon in the analysis. To suppress non-collision background events, a reconstructed collision vertex with at least three associated ID tracks is required.

6.2 Inner detector track selection

Tracks in the ID are required to satisfy conditions on the number of hits in the silicon detectors in order to qualify as a muon candidate. They must have at least two pixel hits, including at least one in the b-layer, and at least six SCT hits. In order to reduce inefficiencies due to known inoperative sensors,Footnote 3 the latter are counted as hits for tracks crossing them. Within \(|\eta | < 1.9\), a good-quality extension of the muon trajectory into the TRT is enforced by requirements on the numbers of associated good TRT hits and TRT outliers.

The TRT outliers appear in two forms in the track reconstruction: as straw tubes with a signal from tracks other than the one in consideration, or as a set of TRT measurements in the extrapolation of a track which fail to form a smooth trajectory together with the pixel and SCT measurements. The latter case is typical of a hadron decay-in-flight, and can be rejected by requiring that the outlier fraction (the ratio of outliers to total TRT hits) is less than 90 %. In the region \(|\eta | <1.9\) the sum of the numbers of TRT hits and outliers is required to be greater than five, with an outlier fraction less than 90 %. At higher \(|\eta |\) the requirement on the total number of TRT hits and outliers is not applied, but tracks which do pass it are also required to pass the cut on the outlier fraction. These quality cuts suppress fake tracks and discriminate against muons from \(\pi / K \) decays.

6.3 Tag selection

For each of the two reconstruction chains, tag muons are defined as CB muons from the interaction vertex. Different selection cuts are applied for the measurements using \(J/\psi \!\rightarrow \!\mu ^+\mu ^-\) and \(Z\!\rightarrow \!\mu ^+\mu ^-\) decays to account for the different kinematics and final-state topologies. For the studies with \(J/\psi \!\rightarrow \!\mu ^+\mu ^-\) a tag muon has to pass the following requirements:

  • the tag muon triggered the readout of the event;

  • \(p_{\mathrm {T}}> 4\) GeV, \(|\eta | < 2.5\);

  • the distance of closest approach of the muon to the primary vertex, in the transverse plane, has transverse coordinate \(|d_0|<0.3\) mm, and longitudinal coordinate \(|z_0|<1.5\) mm, and significances \(|d_0|/\sigma (d_0)<3\), \(|z_0|/\sigma (z_0)<3\), respectively.

For the studies with \(Z\rightarrow \mu ^+\mu ^-\) decays an additional quantity is used, namely track isolation

$$\begin{aligned} T_{\mathrm {isol}}^{\Delta R<0.4}= \sum { p _\mathrm{{T}}(\Delta R <0.4)} / p _{\mathrm {T}}(\mathrm {tag}), \end{aligned}$$
(2)

where the sum extends over all tracks with \(p_{\mathrm {T}}>1\) GeV (excluding the track on which the tag was based), within a cone of \(\Delta R\equiv \sqrt{(\Delta \eta )^2+(\Delta \phi )^2} = 0.4\) around the tag. A tag muon must pass the following requirements:

  • the tag muon triggered the readout of the event (restricting the tag muon to the trigger acceptance, \(|\eta |<2.4\));

  • \(p_{\mathrm {T}}>20\) GeV;

  • \(T_{\mathrm {isol}}^{\Delta R<0.4}<0.2\).

6.4 Probe selection

Probes are either SA muons or ID tracks, depending on which efficiency measurement is being performed. They have to satisfy the following criteria for studies using \(J/\psi \rightarrow \mu ^+\mu ^-\) decays:

  • an ID track fulfilling the hit requirements described in Sect. 6.2 (SA muons are not used, as the ID efficiency is not measured using these decays);

  • reconstructed momentum, \(p>3\hbox { GeV}\), \(|\eta |<2.5\);

  • the tag and the probe are oppositely charged;

  • the tag and the probe must be associated with the same vertex;

  • \(\Delta R<3.5\) between the tag and probe.

  • the invariant mass of the tag-and-probe pair is within the range of \(2<m<3.6\) GeV

Different cuts are applied in case of \(Z\rightarrow \mu ^+\mu ^-\) decays:

  • an ID track fulfilling the hit requirements or a SA muon with at least one \(\phi \) measurement;

  • \(p_{\mathrm {T}}>20\) GeV, \(|\eta |<2.5\);

  • the tag and the probe are oppositely charged;

  • the tag and the probe are associated with the same vertex;

  • azimuthal separation of the tag and the probe, \(\Delta \phi >2.0\);

  • \(T_{\mathrm {isol}}^{\Delta R<0.4}<0.2\);

  • the invariant mass of the tag-and-probe pair is within 10 GeV of \(m_Z\).

6.5 Matching of probes to ID tracks and muons

After selecting all tag-and-probe pairs, an attempt is made to match probe tracks to the objects for which the efficiency is to be measured, i.e. SA probe tracks to ID tracks in the case of the ID efficiency, or ID tracks to CB or CB+ST muons in the case where the reconstruction efficiencies for these two classes of muons are investigated. A match between an ID probe and a reconstructed muon is considered successful if they have the same charge and are close in (\(\eta ,\phi \)) space: \(\Delta R\le 0.01\). Similarly, a match between an SA probe and an ID track is considered successful if \(\Delta R\le 0.05\).

7 Low-\(p_{\mathrm {T}}\)  reconstruction efficiency measured with \(J/\psi \rightarrow \mu ^+\mu ^-\) decays

Figures 6 and 7 show the reconstruction efficiencies for chain 1 and chain 2 with respect to ID tracks with momentum \(p>3\) GeV, as a function of the probe \(p_{\mathrm {T}}\), for the five bins in probe \(|\eta |\) described in Sect. 5. Also shown are the Monte Carlo predictions, which agree with data within the statistical and systematic uncertainties of the measurements.

Fig. 6
figure 6

Efficiency for chain 1 CB and CB+ST muons with momentum \(p>3\) GeV (from \(J/\psi \) decays), as a function of \(p_{\mathrm {T}}\), for five bins in \(|\eta |\) as described in the legend, for data and MC events. The error bars represent the statistical uncertainties while the bands around the data points represent the statistical and systematic uncertainties added in quadrature

Fig. 7
figure 7

Efficiency for chain 2 CB and CB+ST muons with momentum \(p>3\) GeV (from \(J/\psi \) decays), as a function of \(p_{\mathrm {T}}\), for five bins in \(|\eta |\) as described in the legend, for data and MC events. The error bars represent the statistical uncertainties while the band around the data points represents the statistical and systematic uncertainties added in quadrature

A number of checks were performed to study the dependence of the results on analysis details and assumptions.

  1. 1.

    Signal shapes: the means and the widths of the two (matched and unmatched) Gaussians in the fit were allowed to vary independently.

  2. 2.

    Background shape: a linear background function was used in the fit, instead of the quadratic parameterisation; in this case the fit was performed in the reduced mass range of 2.7–3.5 GeV (instead of 2.0–3.6 GeV).

  3. 3.

    Alternative fit: an independent fit to the matched and the total (matched + unmatched) distributions, rather than to matched and unmatched, was used and the efficiency estimated as the ratio of the signal normalisations in the two distributions. While this option does not provide for an easy propagation of the uncertainty from the background subtraction and does not directly account for the correlations between the two samples, it profits from a higher stability of the two simpler fits, whereas the default method needs some care in the choice of the initial conditions, in particular in cases of very high efficiency or small overall sample size.

The largest positive and negative variations obtained from any of the three checks were taken as systematic uncertainties and added in quadrature to the statistical uncertainty to obtain the total upper and lower uncertainties. The statistical uncertainties were found to be at the level of a few percent.

8 Intermediate- and high-\(p_{\mathrm {T}}\) reconstruction efficiencies measured with \(Z\rightarrow \mu ^+\mu ^-\) decays

For higher momentum muons, with \(p_{\mathrm {T}}>20\) GeV, \(Z\) decays are used to measure the reconstruction efficiencies.

8.1 Inner detector reconstruction and identification efficiency

Figure 8 shows the reconstruction and identification efficiency in the ID as a function of \(\eta \), for data and simulation, as determined using SA probes. The simulation includes all considered backgrounds. The scale factors (SF), defined as the ratio of the data efficiency to the Monte Carlo efficiency, are displayed in the lower panel (the smallness of the background correction, as described in Sect. 4, means that its effect on the SF is negligible).

Fig. 8
figure 8

Measured ID reconstruction and identification efficiency for muons (from \(Z\) decays), as a function of \(\eta \), for data and Monte Carlo simulation. The scale factors (SF), defined as the ratio of the measured efficiency to the predicted efficiency, are shown in the lower panel of the plot. The uncertainties are statistical. The systematic uncertainty is discussed in Sect. 8.4

As discussed earlier, the efficiency for the combined reconstruction varies with the detector region, and with \(p_{\mathrm {T}}\) in the range below 6 GeV. In contrast, the ID reconstruction efficiency is independent of \(\phi \) and \(p_{\mathrm {T}}\) [3], and shows only a slight dependence on \(\eta \).

The slightly lower efficiencies at \(\eta \sim 0\) and \(|\eta |\sim 1\) are caused by the ID hit requirements for muon identification described in Sect. 6.2: at \(\eta \sim 0\), ID tracks pass through an inactive region near the middle of the TRT barrel where straws produce no TRT hits; at \(|\eta |\sim 1\), there is a small region in the transition between the barrel and the end-caps of the ID in which muons cross fewer than six SCT sensors [3]. The measured ID muon reconstruction and identification efficiencies agree with the Monte Carlo predictions within 1 %, and, for the most part, within the statistical uncertainties. The average ID efficiency is 0.991 \(\pm \) 0.001 with the small loss being due to the hit requirements imposed on the ID muon tracks. These results are independent of the choice of the algorithm chain for the stand-alone muon.

8.2 Reconstruction efficiencies for CB muons

Figure 9 shows the reconstruction efficiency (relative to the ID reconstruction efficiency) for CB muons as a function of the detector region, \(p_{\mathrm {T}}\) and \(\eta \), for data and simulation (with all considered backgrounds included). The scale factors are displayed in the lower panel of each plot.

Fig. 9
figure 9

Reconstruction efficiencies (relative to the ID efficiency) and scale factors for CB muons (from \(Z\) decays) as a function of detector region, muon \(p_{\mathrm {T}}\) and muon \(\eta \) as indicated in the figure. The efficiencies for the two reconstruction chains, obtained from data (without background correction) and Monte Carlo simulation (including backgrounds) are shown in the upper part of each figure. The corresponding scale factors are shown in the lower panels. The uncertainties are statistical only. The systematic uncertainties are discussed in Sect. 8.4

The mean value of the \(\eta \)-dependent scale factor is \(0.989\pm 0.003\) for chain 1 and \(0.995\pm 0.002\) for chain 2 , where the errors are statistical. The 1 % deviation from unity in the overall efficiency scale factor of chain 1 is caused mainly by the data/MC disagreement in the transition region (\(\mathrm {SF}=0.94\)). The lower data efficiency in the transition region is attributed to the limited accuracy of the magnetic field map used in the reconstruction of the ATLAS data in this region, which leads to a small mis-measurement of the stand-alone muon momentum. This in turn may affect the combination of the MS and ID tracks, as their momenta may not be compatible. The transition region efficiency drop can be recovered, and the overall efficiency significantly increased by including ST muons, which are tagged by only one muon layer, as described in detail below.

The scale factors determined in bins of \(p_\mathrm {T}\) agree, within 1.5 standard deviations, with the average scale factor for the algorithm in question.

The background-corrected efficiencies for CB muons are shown in Fig. 10. The background is estimated from Monte Carlo simulation, as described in Sect. 4, and is subtracted bin by bin. The average CB muon reconstruction efficiency is \(0.928\pm 0.002\) for chain 1 and \(0.958\pm 0.001\) for chain 2 . The difference in efficiency between the two chains arises mainly from the more stringent requirements on the reconstructed MS tracks in chain 1. The ratios between data and MC efficiencies are almost identical to the SFs already discussed for Fig. 9 as a consequence of the smallness of the background correction.

Fig. 10
figure 10

Background-corrected efficiencies for CB muons (from \(Z\) decays) as a function of detector region, muon \(p_{\mathrm {T}}\) and muon \(\eta \) as indicated in the figure, obtained from data and Monte Carlo simulation for the two reconstruction chains. The uncertainties are statistical only. The systematic uncertainties are discussed in Sect. 8.4

8.3 Reconstruction efficiencies for CB+ST muons

The degree to which segment tagging can recover some muons, in particular in detector regions with only partial muon coverage, is studied by measuring the efficiency for CB+ST muons. The same tag-and-probe method is used with the only difference being that the probe is matched to a CB or ST muon. Figure 11 shows the measured CB+ST muon efficiencies as functions of the detector region, \(p_{\mathrm {T}}\) and \(\eta \), in comparison with the corresponding CB muon efficiencies. The gains in efficiency when using ST muons in addition to the CB muons are presented in the lower panels of the plots. These are largest in the ATLAS Feet (\(13\,\%\)) and Transition (\(15\,\%\)) regions of the detector for \(\mathrm{chain~1} \). For \(\mathrm{chain~2} \) the largest gain is \(3\,\%\) in the Feet and BEE regions. Figure 11 also shows that the two chains have similar overall efficiencies for CB+ST muons, 0.970 \(\pm \) 0.001 for chain 1 and 0.980 \(\pm \) 0.001 for chain 2 .

Fig. 11
figure 11

Efficiencies for CB+ST muons (from \(Z\) decays) in comparison to those for CB muons only, for the two reconstruction chains and as a function of detector region, muon \(p_{\mathrm {T}}\) and muon \(\eta \) as indicated in the figure. The relative gain is shown in the lower panel of each figure. The uncertainties are statistical only. The systematic uncertainties are discussed in Sect. 8.4

In Fig. 12, the efficiency for CB+ST muons measured from data is compared to the Monte Carlo expectations and scale factors are presented. Remarkable agreement between the measured and predicted efficiencies is achieved. The scale factors for CB+ST muons are 1.003 \(\pm \) 0.002 for chain 1 and 1.001 \(\pm \) 0.002 for chain 2 .

Fig. 12
figure 12

Efficiencies for CB+ST muons (from \(Z\) decays), for the two reconstruction chains as a function of detector region, muon \(p_{\mathrm {T}}\) and muon \(\eta \) as indicated in the figure. The efficiencies are obtained from data with background correction and from Monte Carlo simulation of the signal. The corresponding scale factors are shown in the lower panel of each plot. The uncertainties are statistical only. The systematic uncertainty is discussed in Sect. 8.4

8.4 Systematic uncertainties

Uncertainties on the background contributions and on the resolution of the detector are considered as sources of systematic uncertainties. The uncertainty due to the description of the finite detector resolution is estimated by varying the selection cuts when determining the efficiencies from MC-simulated data. For CB muons, the cuts on the mass window around \(m_Z\) and the cut on the transverse momentum of the tag are each varied within \(\pm 1\sigma \) of the \(m_{{\mu ^+ \mu ^-}}\) and \(p_{\mathrm {T}}\) resolutions. Other cuts are varied by \(\pm 10\,\%\). The resulting changes in the scale factors are quoted as systematic uncertainties. The normalisation of the background contribution inside the mass window is varied by \(\pm 10\,\%\) and the resulting differences in the scale factors are considered as additional systematic uncertainties. The individual uncertainties are considered to be uncorrelated and are added in quadrature to estimate the total systematic uncertainty. For values which result from an upwards and downwards variation, the larger value is used. The largest contribution arises from the level of background contamination, which depends primarily on the choice of the mass window and the normalisation of the backgrounds. Another important contribution is due to the variation of the probe isolation criteria. The overall systematic uncertainty on the CB muon efficiency is \(0.2\,\%\) for both chains.

As the same tag-and-probe selection is used for the measurements of the CB+ST muon efficiencies, the same systematic uncertainties are expected for the corresponding scale factors. The systematic uncertainties on the ID muon efficiency scale factors are substantially smaller, principally due to the high purity of the MS probe muons.

9 Measurement of the muon isolation efficiency

Muon isolation is a powerful tool for a high-purity event selection in many physics analyses, and is also used for rejecting muons from hadron decays in the \(Z\) decay tag-and-probe analyses presented here. It is therefore desirable to quantify the reliability of the Monte Carlo prediction of the isolation efficiency (simulated using PYTHIA).Footnote 4 This is studied using the same event selection that was used for the reconstruction efficiency measurements, up to and including the selection of the tag muon (the specific chain used is not shown, since the performance is comparable for both). In this case, the probe muon is defined as a CB muon with \(p_{\mathrm {T}}>20\) GeV that fulfils the ID hit requirements described in Sect. 6. We consider the following isolation variables:

  • track isolationFootnote 5—the summed \(p_{\mathrm {T}}\) of tracks (excluding that of the muon) in cones of size \(\Delta R=0.3\) and \(\Delta R=0.4\) around the muon, divided by the \(p_{\mathrm {T}}\) of the muon;

  • calorimeter isolation—the transverse energy (\(E_\mathrm {T}\)) deposition in the calorimeter in cones of size \(\Delta R=0.3\) and \(\Delta R=0.4\) around the muon (with the muon’s energy loss subtracted [1, p. 194]), divided by the \(p_{\mathrm {T}}\) of the muon.

The tag-and-probe selections, as described in Sect. 6, only make use of \({T_{\mathrm {isol}}^{\Delta R<0.4}<0.2}\). However, the choice of isolation criteria depends on the analysis and this section presents the comparisons of data and Monte Carlo simulations for the following combinations of isolation variables:

  • \(T_{\mathrm {isol}}^{\Delta R<0.4}<0.2\) and \(E_\mathrm {T}^{\Delta R<0.4} / p_{\mathrm {T}}(\mu )<0.2\);

  • \(T_{\mathrm {isol}}^{\Delta R<0.4}<0.1\) and \(E_\mathrm {T}^{\Delta R<0.4} / p_{\mathrm {T}}(\mu )<0.1\);

  • \(T_{\mathrm {isol}}^{\Delta R<0.3}<0.1\) and \(E_\mathrm {T}^{\Delta R<0.3} / p_{\mathrm {T}}(\mu )<0.1\).

Figure 13 compares the distributions of the measured isolation variables for the probe muons with the Monte Carlo predictions. The experimental and simulated distributions agree well, leading to a reliable prediction as a function of \(p_{\mathrm {T}}\), of the isolation efficiency, which is defined as the fraction of probe muons passing a given set of isolation cuts.

Fig. 13
figure 13

Comparison of the measured track isolation (left) and calorimeter isolation (right) distributions of the probe muon (from \(Z\) decays) with the Monte Carlo predictions, for two different cone sizes using the isolation variables defined in the text. The upper and lower plots correspond to \(\Delta R=0.4\) and \(\Delta R=0.3\), respectively. The simulation includes the effects of pile-up, as described in the text. The uncertainties are statistical

The measured isolation efficiencies and the corresponding Monte Carlo predictions are compared for chain 1 in Fig. 14; the results for chain 2 are consistent. Experimental and simulated data agree within uncertainties. The lower efficiencies at low \(p_{\mathrm {T}}\) are mainly caused by the fact that the \(p_{\mathrm {T}}\) and \(E_\mathrm {T}\) sums, which depend only weakly on the muon \(p_{\mathrm {T}}\), are divided by this quantity, leading to isolation variables that rise with decreasing muon \(p_{\mathrm {T}}\). They are also partially due to the background, which is non-negligible in the low-\(p_{\mathrm {T}}\) region.

Fig. 14
figure 14

Isolation efficiencies for muons from \(Z\) decays as a function of \(p_{\mathrm {T}}\), for track isolation (left) and calorimeter isolation (right) requirements with different isolation cone radii, \(\Delta R\), as described in the legend. The Monte Carlo predictions include background processes as well as the \(Z\) signal. The uncertainties are statistical only

10 Measurement of the muon momentum resolution

The muon momentum resolution of the ATLAS detector depends on the \(\eta \), \(\phi \), and \(p_{\mathrm {T}}\) of the muon [3]. In the ID, the \(p_{\mathrm {T}}\) dependence of the relative momentum resolution can be parameterised to a good approximation [1] by the quadratic sum of two terms,

$$\begin{aligned} \begin{aligned} \frac{\sigma _\mathrm {ID}(p_{\mathrm {T}})}{p_{\mathrm {T}}}&= a_\mathrm {ID}(\eta )\oplus b_\mathrm {ID}(\eta )\cdot p_{\mathrm {T}}&\text {for } 0<|\eta |<2.0 ; \\ \frac{\sigma _\mathrm {ID}(p_{\mathrm {T}})}{p_{\mathrm {T}}}&= a_\mathrm {ID}(\eta )\oplus \frac{b_\mathrm {ID}(\eta )\cdot p_{\mathrm {T}}}{\mathrm {tan}^2(\theta )}&\text {for } 2.0<|\eta |<2.5. \end{aligned} \end{aligned}$$
(3)

The first term describes the multiple scattering contribution, whilst the second term describes the intrinsic resolution caused by the imperfect knowledge of the magnetic field in the ID, by the spatial resolution of the detector components, and by any residual misalignment of the detector components. For \(|\eta |>2.0\), the best parameterisation of the second term is given by \({b_\mathrm {ID}(\eta )\cdot p_{\mathrm {T}}}/\mathrm {tan}^2(\theta )\). Measurements (from data) of the material distribution in the ID [11, 12] constrain \(a_\mathrm {ID}(\eta )\) to values which agree with the Monte Carlo prediction to within \(5\) % in the barrel and \(10\) % in the end-caps. The parameter \(b_\mathrm {ID}(\eta )\) is derived from the dimuon invariant mass resolution in \(Z\!\rightarrow \!\mu ^+\mu ^-\) decays.

The stand-alone muon resolution can be parameterised as follows:

$$\begin{aligned} \frac{\sigma _\mathrm {SA}(p_{\mathrm {T}})}{p_{\mathrm {T}}} = a_\mathrm {MS}(\eta ,\phi )\oplus b_\mathrm {MS}(\eta ,\phi )\cdot p_{\mathrm {T}}\oplus \frac{c(\eta ,\phi )}{p_{\mathrm {T}}}~, \end{aligned}$$
(4)

where the first two terms parameterise the effect of the multiple scattering and the contribution of the intrinsic momentum resolution of the MS, respectively. The third term parameterises the effect of the fluctuations of the muon energy loss in the calorimeters, but this is small for the momentum range under consideration and is fixed to the value predicted by MC simulation.

A special data set, recorded in 2011, with no toroidal magnetic field in the MS, was used to simulate high-momentum (i.e. straight) tracks and estimate \(b_\mathrm {MS}(\eta ,\phi )\), yielding \(b_\mathrm {MS}(\eta ,\phi )\sim 0.2\) TeV\(^{-1}\) in the barrel and the MDT end-cap region (excluding the transition region) and \({\sim }0.4\) TeV\(^{-1}\) in the CSC end-cap region, with a relative accuracy of about 10 % in both regions. This special data set made it possible to improve the alignment of the muon chambers, leading to \(b_\mathrm {MS}(\eta ,\phi )\lesssim 0.2\) TeV\(^{-1}\) everywhere in the MS in 2011.

Figure 15 shows the dimuon invariant mass resolution of the ID in \(Z\!\rightarrow \!\mu ^+\mu ^-\) decays as a function of the pseudorapidity interval of the decay muons, where both are required to lie in the same interval. The mass resolution is the width of a Gaussian which, when convolved with the generator-level dimuon invariant mass, reproduces the dimuon invariant mass distribution observed in data. The ID dimuon invariant mass resolution is best in the barrel, where it is about 2 GeV, is better than 3 GeV for \(|\eta |<2.0\) and degrades to about 6 GeV for \(2.0<|\eta |<2.5\). The degradation of the mass resolution with increasing \(|\eta |\) is primarily caused by the fact that as \(|\eta |\) increases there is a lower field integral per track. That the dimuon invariant mass resolution measured in experimental data is worse than predicted (typically by about 30 %), is attributed to residual internal misalignments of the ID. The internal alignment of the ID was performed by minimising track residuals. This procedure has certain ambiguities which can be resolved by adding constraints such as the requirement that the energy/momentum ratio (\(E/p\)) distributions of electrons and positrons be the same. These constraints were only introduced into the alignment procedure for the 2011 data [13], in which a significantly improved dimuon invariant mass resolution is observed.

Fig. 15
figure 15

The dimuon invariant mass (\(m_{\mu \mu })\) resolutions in \(Z\!\rightarrow \!\mu ^+\mu ^-\) decays in the data and in the MC as a function of \(\eta \) region with both decay muons in the same \(\eta \) region, for the ID (left) and MS (right). The simulation assumes a perfectly aligned ATLAS detector

Due to the toroidal magnetic field, the relative momentum resolution of SA muons (and hence the corresponding dimuon invariant mass resolution—as shown in Fig. 15) is expected to be independent of the \(\eta \) of the decay muons, except in the magnet transition region (\(1.05<|\eta |<1.7\)) where the magnetic field in the MS is highly non-uniform, with a field integral approaching 0 in certain (\(\eta ,\phi \)) regions [3]. Furthermore, some chambers in the region \(1.05<|\eta |<1.3\) were not yet installed,Footnote 6 which means that the momentum measurement relies on only two layers of chambers, causing a significant degradation in the momentum resolution.

Figure 15 also shows that the MS dimuon invariant mass resolution is consistently worse in data than in simulation (typically between 30 and 50 % worse, depending on \(\eta \) region). Two sources for this effect were identified.

  1. 1.

    Asymmetry of the magnetic field: in the MC simulation, a perfectly aligned detector is assumed. In reality, the two end-cap toroid systems are not symmetric with respect to the plane orthogonal to the major axis of the ID, and situated at the centre of the detector. This small asymmetry translates into an asymmetry of the magnetic field integrals, in particular in the transition regions. The reconstruction of the 2010 data with a corrected field map improves the dimuon invariant mass resolution in the transition region by 0.4 GeV.

  2. 2.

    Residual misalignment of the muon chambers: even after the MS alignment procedures are applied, residual misalignments remain, which limit the attainable momentum resolution. The analysis of a special set of 2011 data with no magnetic field in the MS was used to produce a Monte Carlo simulation of \(Z\!\rightarrow \!\mu ^+\mu ^-\) events with the addition of a realistic residual misalignment of the MS. The results of this simulation are in agreement with the experimentally determined invariant mass resolutions.

The dimuon invariant mass resolution obtained with CB muons profits from the complementary momentum measurements of the ID and MS. As shown in Fig. 16, a dimuon invariant mass resolution between 1.4 and 2.5 GeV is achieved, with little dependence on \(\eta \).

Fig. 16
figure 16

Dimuon invariant mass (\(m_{\mu \mu })\) resolution for combined muons in \(Z\!\rightarrow \!\mu ^+\mu ^-\) decays in the data and in the MC as a function of \(\eta \) region with both decay muons in the same \(\eta \) region. The simulation assumes a perfectly aligned ATLAS detector

The measured dimuon invariant mass resolutions can be translated into muon momentum resolutions. This was done by smearing the generated muon momenta, according to Eqs. (3) and (4), by the amounts necessary to reproduce the measured dimuon invariant mass resolutions. Only the parameters \(b_\mathrm {ID}(\eta )\) and \(a_\mathrm {MS}(\eta ,\phi )\) were varied during this procedure. The parameter \(a_\mathrm {ID}(\eta )\) was set to the Monte Carlo prediction and varied within its uncertainty (see above) to evaluate the impact on the result for \(b_\mathrm {ID}(\eta )\). The parameter \(b_\mathrm {MS}(\eta ,\phi )\) was set to the value derived from the special straight-track data set while \(c(\eta ,\phi )\) was set to its predicted value. In order to gain additional sensitivity to the momentum resolutions of the ID and MS, in addition to the dimuon mass spectrum of \(Z\) boson decays, the distributions of the differences between the ID and SA momenta of muons from \(W\!\rightarrow \!\mu \nu _\mu \) decays were compared between the experimental and smeared MC data. The \(W\) boson selection and the MC samples for the analysis are the same as in Ref. [10]. As the use of \(W\) boson decays correlates the SA and ID resolutions, those are extracted simultaneously in the fit. The results are displayed for the different detector regions in Figs. 17 and 18, with the uncertainty of the curves computed from the uncertainties of the parameters in the resolution functions [Eqs. (3) and (4)]. Also shown is the expected resolution beyond the region in \(p_{\mathrm {T}}\) probed by the \(Z\)-boson decays. As discussed earlier, the momentum resolution in experimental data is worse than in the Monte Carlo simulation, which is attributed, in part, to the residual misalignments of the ID and MS.

Fig. 17
figure 17

Muon momentum resolution as a function of \(p_{\mathrm {T}}\) for different barrel and transition \(|\eta |\) regions as denoted in the legend. The dot-dash line is from a simulation which assumes perfect alignment of the ATLAS detector, whilst the solid/dotted line shows simulation smeared to reproduce the invariant mass resolution measured in data. The solid section of the line shows the \(p_{\mathrm {T}}\) range measured by \(Z\) and \(W\) decays, and the dotted section the ‘extrapolation’ regions. The shaded bands show the uncertainty of the curves, computed from the uncertainties of the parameters derived in the resolution functions shown in Eqs. (3) and (4)

Fig. 18
figure 18

Muon momentum resolutions as a function of \(p_{\mathrm {T}}\) for different end-cap \(|\eta |\) regions as denoted in the legend. For the ID region \(|\eta |\,{>}\,2.0\) (bottom-left), the best parameterisation of the resolution depends on \(p_{\mathrm {T}}/\mathrm {tan}^2(\theta )\) instead of \(p_{\mathrm {T}}\). The dot-dash line is simulation which assumes perfect alignment of the ATLAS detector, whilst the solid/dotted line shows simulation smeared to reproduce the invariant mass resolution measured in data. The solid section of the line shows the \(p_{\mathrm {T}}\) range measured by \(Z\) decays, and the dotted section the ‘extrapolation’ regions. The shaded bands show the uncertainty of the curves, which are computed from the uncertainties of the parameters in the resolution functions shown in Eqs. (3) and (4)

11 Summary

The ATLAS muon reconstruction efficiencies were studied with \(J/\psi \!\rightarrow \!\mu ^+\mu ^-\) and \(Z\!\rightarrow \!\mu ^+\mu ^-\) decays using 40 pb\(^{-1}\) of \(\sqrt{s}=7\) TeV \(pp\) LHC collision data recorded in 2010.

Samples of \(J/\psi \) and \(Z\) decays were used to access the transverse momentum regions of \(p_{\mathrm {T}}<10\) GeV and \(20\) GeV\({}<p_{\mathrm {T}}<100\) GeV respectively. The muon reconstruction efficiency is found to be \({>}96\,\%\) and agrees with the MC prediction to better than 1 %. The reconstructed quantities used to ensure muon isolation are shown to be well modelled in Monte Carlo simulations, and the corresponding muon isolation efficiencies are in excellent agreement with the MC predictions.

The muon momentum resolutions for \(p_{\mathrm {T}}>20\) GeV are derived from the dimuon mass resolutions in \(Z\!\rightarrow \!\mu ^+\mu ^-\) decays and from the differences between the ID and SA momenta of muons from \(W\!\rightarrow \!\mu \nu _\mu \) decays. The resolutions are worse in data than in simulation for the entire momentum range considered. For instance, at \(p_{\mathrm {T}}\approx 30\) GeV and \(1.7<|\eta |<2.0\) the resolutions in experimental data are found to be about 30 % worse than predicted by the simulation. These differences are attributed to mis-modelling of the magnetic field and residual misalignments of the inner detector and muon spectrometer. An improved magnetic field map was used from 2011 onwards, and there have since been several iterations of the alignment.