Measurement of the top quark mass using single top quark events in proton-proton collisions at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s}= 8$$\end{document}s=8 TeV

A measurement of the top quark mass is reported in events containing a single top quark produced via the electroweak t channel. The analysis is performed using data from proton-proton collisions collected with the CMS detector at the LHC at a centre-of-mass energy of 8 TeV, corresponding to an integrated luminosity of 19.7 fb\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{-1}$$\end{document}-1. Top quark candidates are reconstructed from their decay to a \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm {W}$$\end{document}W boson and a b quark, with the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm {W}$$\end{document}W boson decaying leptonically to a muon and a neutrino. The final state signature and kinematic properties of single top quark events in the t channel are used to enhance the purity of the sample, suppressing the contribution from top quark pair production. A fit to the invariant mass distribution of reconstructed top quark candidates yields a value of the top quark mass of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$172.95 \pm 0.77\,\text {(stat)} ^{+0.97}_{-0.93}\,\text {(syst)} \,\text {GeV} $$\end{document}172.95±0.77(stat)-0.93+0.97(syst)GeV. This result is in agreement with the current world average, and represents the first measurement of the top quark mass in event topologies not dominated by top quark pair production, therefore contributing to future averages with partially uncorrelated systematic uncertainties and a largely uncorrelated statistical uncertainty.


Introduction
All previously published measurements of the top quark mass have been obtained using samples of top quark-antiquark pairs. A combination of measurements from the CDF and D0 experiments at the Tevatron and ATLAS and CMS experiments at the LHC yields a value of 173.34 ± 0.27 (stat) ± 0.71 (syst) GeV for the top quark mass m t [1]. Measuring m t in single top quark production enriches the range of available measurements, exploiting a sample which is almost statistically independent from those used by previous ones, and with systematic uncertainties partially uncorrelated from those considered in tt production. Because of the different production mechanism, the mass extraction is affected differently by the modelling of both perturbative effects, such as initiale-mail: cms-publication-committee-chair@cern.ch and final-state radiation, and nonperturbative effects, such as colour reconnection, in quantum chromodynamics (QCD). Some discussion on these topics, though mainly restricted to the case of pair production, can be found in Refs. [2,3]. In perspective, the lower level of gluon radiation and final state combinatorial arrangements with respect to tt production will make this channel a good candidate for precision measurements of m t when larger samples of events are available.
At the CERN LHC, top quarks are mainly produced as tt pairs, through gluon-gluon fusion or quark-antiquark annihilation, mediated by the strong interaction. The standard model (SM) predicts single top quark production through electroweak processes, with a rate about one third that of the tt production cross section. This has been confirmed by observations at the Tevatron [4] and LHC [5,6].
In this paper, top quark candidates are reconstructed via their decay to a W boson and a b quark, with the W boson decaying to a muon and a neutrino. The event selection is tailored, before looking at data in the signal region, to enhance the single top quark content in the final sample and so have a result as independent as possible from those obtained using tt events.
The paper is organised as follows. Section 2 describes the CMS detector, followed by information about the data sample and simulation used in the analysis in Sect. 3. The selection of events and the reconstruction of the top quark candidates is given in Sect. 4, and the description of the maximumlikelihood fit to derive the top quark mass is in Sect. 5. Section 6 describes the systematic uncertainties affecting the measurement and Sect. 7 summarises the results.

The CMS detector
The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diameter, providing a magnetic field of 3.8 T. Within the solenoid volume are a silicon pixel and strip tracker, a lead tungstate crystal electromagnetic calorimeter, and a brass and scintillator hadron calorimeter, each composed of a barrel and two endcap sections. Forward calorimeters extend the pseudorapidity coverage provided by the barrel and endcap detectors. Muons are measured in gasionisation detectors embedded in the steel flux-return yoke outside the solenoid.
A more detailed description of the CMS detector, together with a definition of the coordinate system used and the relevant kinematic variables, can be found in Ref. [7].

Data and simulated samples
The measurement reported here is performed using the √ s =8 TeV proton-proton collision data sample collected in 2012 with the CMS detector, corresponding to an integrated luminosity of 19.7 fb −1 .
At the lowest order in perturbation theory, single top quark production proceeds through the t-channel, s-channel, and associated tW production modes. The t channel provides the largest contribution to the single top quark cross section. The corresponding amplitude can be calculated using one of two different schemes [8][9][10]: in the 5-flavour scheme, b quarks are considered as coming from the interacting proton, and the leading-order (LO) diagram is a 2 → 2 process (Fig. 1, Fig. 1 Feynman diagrams representing the dominant single top quark production mechanisms in the t channel upper); in the 4-flavour scheme, b quarks are not present in the initial state, and the LO diagram is a 2 → 3 process (Fig. 1, lower). The predicted t-channel single top quark cross section for pp collisions at a centre-of-mass energy of 8 TeV is σ t = 54.9 +2.3 −1.9 pb for the top quark and σ t = 29.7 +1.7 −1.5 pb for the top antiquark. These values are obtained by a next-to-leading-order (NLO) calculation in quantum QCD with hathor v.2.1 [11,12], assuming a top quark mass of 172.5 GeV. The parton distribution functions (PDFs) and α S uncertainties are calculated using the PDF4LHC prescription [13,14] with the MSTW2008 68% confidence level (CL) NLO [15,16], CT10 NLO [17], and NNPDF2.3 [18] PDF sets. At 8 TeV, the predicted tt production cross section is σ (tt) = 252.9 +6.4 −8.6 (scale) ± 11.7 (PDF + α S ) pb as calculated with the Top++2.0 program to next-to-next-to-leading order in perturbative QCD, including soft-gluon resummation to next-to-next-to-leading-log order (see Ref. [19] and references therein), and assuming a top quark mass of 172.5 GeV. In this calculation, the total scale uncertainty is obtained from the independent variation of the factorisation and renormalisation scales, μ F and μ R , by a factor 2 and 1/2; the total PDF and α S uncertainties are estimated following the PDF4LHC prescription [14] with the MSTW2008 68% CL NNLO [16], CT10 NNLO [18], and NNPDF2.3 [20] FFN PDF sets.
Simulated events are used to optimise the event selection and to study the background processes and the expected performance of the analysis. The signal t-channel events are generated with the powheg generator, version 1.0 [21], in the 5-flavour scheme, interfaced with pythia [22], version 6.426, for parton showering and hadronisation. Single top quark s-channel and tW associated production are considered as backgrounds for this measurement and simulated with the same generator. Top quark pair production, single vector boson production associated with jets (referred to as W/Z+jets in the following), and double vector boson (diboson) production are amongst the background processes taken into consideration and have been simulated with the MadGraph generator, version 5.148 [23], interfaced with pythia for parton showering. The pythia generator is used to simulate QCD multijet event samples enriched with isolated muons. The value of the top quark mass used in all simulated samples is 172.5 GeV. All samples are generated using the CTEQ6.6M [24] PDF set and use the Z2* underlying event tune [25]. The factorisation and renormalisation scales are both set to m t for the single top quark samples, while a dynamic scale is used for the other samples, defined as the sum in quadrature of the transverse momentum ( p T ) and the mass of the particles produced in the central process. The passage of particles through the detector is simulated using the Geant4 toolkit [26]. The simulation includes addi-tional overlapping pp collisions (pileup) with a multiplicity that is tuned to match the one observed in data.

Event selection and reconstruction
Signal events are characterised by a single isolated muon, momentum imbalance due to the presence of a neutrino, and one central b jet from the top quark decay. In addition, events often feature the presence of a light quark jet in the forward direction, from the hard-scattering process.
The online selection requires the presence of one isolated muon candidate with p T greater than 24 GeV and absolute value of the pseudorapidity (η) below 2.1. Events are required to have at least one primary vertex reconstructed from at least four tracks, with a distance from the nominal beaminteraction point of less than 24 cm along the z axis and less than 2 cm in the transverse plane. In cases where more than one primary vertex is found, the one featuring the largest value of p 2 T is retained ("leading vertex"), where the sum runs over all the tracks assigned to that vertex.
All particles are reconstructed and identified with the CMS particle-flow algorithm [27,28]. Muon candidates are further required to have p T > 26 GeV, thus ensuring they are selected in the region of maximal trigger efficiency. Muon candidates are also required to be isolated. This is ensured by requiring I rel < 0.12, where I rel is defined as the sum of the transverse energies deposited by long-lived charged hadrons, photons, and neutral hadrons in a cone of size R = √ ( η) 2 + ( φ) 2 = 0.4 around the muon direction (φ being the azimuthal angle, in radians), divided by the muon p T itself. An offset correction is applied to remove the additional energy included in the jets that come from pileup [29]. Events are rejected if an additional muon (electron) candidate is present, passing the selection criteria p T > 10 (20) GeV, |η| < 2.5, and I rel < 0.2 (0.15).
To define jets, the reconstructed particles are clustered using the anti-k T algorithm [30] with a distance parameter of 0.5. Charged particles are excluded if they originate from a primary vertex that is not the leading vertex. The energy deposition in the jet due to neutral pileup particles is inferred and subtracted by considering charged pileup particles inside the jet cone. Additional corrections to the jet energies are derived from the study of dijet events and photon+jets events [31]. Jets are required to have |η| < 4.7 and to have a corrected transverse energy greater than 40 GeV. Jets associated with the hadronisation of b quarks ("b jets") are identified using a b tagging algorithm based on the 3D impact parameter of the tracks in the jet to define a b tagging discriminator [32]. The threshold for this variable is chosen such that the probability to misidentify jets coming from the hadronisation of light quarks (u, d, s) or gluons is small (0.1%), while ensuring an efficiency of 46% for selecting jets coming from b quarks, as determined from the simulation of events with top quark topologies. Event weights are applied to adjust the b jet yields in the simulation to account for differences in the b tagging efficiency between data and simulation.
The missing transverse momentum ( p miss T ) is calculated as the negative vector sum of the transverse momenta of all reconstructed particles. Corrections to the jet energies, as well as an offset correction accounting for pileup interactions, are propagated to p miss T . The missing transverse momentum magnitude ( p miss T ) is required to exceed 50 GeV, to suppress the QCD multijet background.
To reject jets from pileup, non b-tagged jets are rejected if the root-mean-square η-φ radius of the particles constituting the jet with respect to the jet axis is larger than 0.025. To suppress background from QCD multijet events, the transverse mass of the W boson m T (W) must be larger than 50 GeV, where m T (W) is constructed from the missing transverse momentum and muon transverse momentum vectors as (1) The same event reconstruction and selection of top quark candidates adopted by the CMS single top quark t-channel cross section measurement at 8 TeV in Ref. [5] is used. Due to the detector acceptance and jet selection requirements, most signal events are characterised by the presence of two reconstructed jets, one of which comes from the hadronisation of a b quark. Therefore, events with two reconstructed jets, exactly one of which is b tagged, constitute the "signal region" (referred to as '2J1T' in the following). Other event topologies are used to study background properties: the sample with two reconstructed jets, neither of which is b tagged ('2J0T') is dominated by W+jets events; the sample with three reconstructed jets, where two jets are b tagged ('3J2T') is dominated by tt events. For all topologies considered, the jet with the highest value of the b tagging discriminator is used for top quark reconstruction, while that with the lowest value is taken to be the light-quark jet associated with top quark production ( Fig. 1).
To enrich the sample in single top quark events, further requirements are applied to variables that exhibit good discriminating power with respect to tt events, as described below. The selection criteria have been chosen after studying their effect on the purity of the sample, while verifying that the statistical uncertainty achievable on the top quark mass would not be excessively degraded.
A feature of single top quark production in the t channel is that the top quark is accompanied by a light-quark jet (the quark labelled q in Fig. 1), which is produced in a more forward direction than jets coming from tt production or other background processes. This is reflected in the distribution of In t-channel single top quark production, top quarks are produced more frequently than top antiquarks due to the charge asymmetry of the proton-proton initial state [33], as seen in the muon charge distribution (Fig. 2, lower). To obtain as pure a sample as possible, only events with positively charged muons are retained.

Determination of the top quark mass
For each selected event, the top quark mass is reconstructed from the invariant mass m μνb calculated from the muon, the neutrino, and the b jet. The 4-momenta of the muon and the jet are measured, while, for the neutrino, the 4-momentum is determined by using the missing transverse momentum in the event and a kinematical constraint on the μν invariant mass, required to be consistent with the mass m W of the W boson [34]: where E μ is the muon energy, p μ x , p μ y and p μ z are the components of the muon momentum, p ν z is the longitudinal component of the neutrino momentum, and p miss T is used for the transverse components of the neutrino momentum. Equation 2 is quadratic in p ν z : when two real solutions are found, the one with the smallest value of | p ν z | is taken; in the case of complex solutions, the imaginary component is eliminated by modifying p miss T,x and p miss T,y independently, so as to give Figure 3 shows the m μνb distributions before and after the final event selection. According to Monte Carlo simulation, after the final selection, 73% of the reconstructed top quarks come from single top quark production, and of these about 97% come from t-channel production.
The top quark mass is measured with an extended unbinned maximum-likelihood fit to the m μνb distribution. The numbers of events for the various contributions, except for the single top quark t-channel one, are fixed to the values extracted from simulation, taking into account the different theoretical cross sections [5]. The description of the parameterisation of the signal and background components used in the fit is presented below. The free parameters of the fit are the number of single top quark signal events and the parameters of the signal shape.

Parameterisation of top quark components
The shapes of the m μνb distributions for samples where a top quark is present are studied using simulated events.
The tt component exhibits a wider peak, with a larger high-mass tail, compared to the single top quark t-channel component. The simulation shows that the number of muon and b jet pairs correctly assigned to the parent top quark is around 55% for tt events, while this fraction exceeds 90% for signal events. Both contributions can be fitted by Crystal Ball functions [36], with independent parameters μ and σ representing the Gaussian core, and α and n describing where the power-law tail begins and the exponent of the tail, respectively. The distributions obtained from the simulated samples before the final selection are shown in Fig. 4. The difference between the values of the μ parameter of the Crystal Ball function obtained from the fits is The parameter μ of the Crystal Ball function describing the single top quark t-channel component is used to estimate the top quark mass. The mass is obtained by shifting the value of μ resulting from the fit by an amount m depending on μ itself. In order to calibrate the magnitude of the shift, the fit has been repeated on a set of simulated samples including all signal and background processes, where the t-channel single top quark and tt events were generated with different values of the top quark mass, all other events remaining unchanged. Figure 5 shows the resulting values of μ as a function of the generated top quark mass (upper) and the mass calibration curve from a fit to these values (lower). The m shift to be applied to the fitted value of μ is expressed as a linear function of μ itself. The shaded grey area represents the uncertainty associated with m, obtained from the statistical uncertainties of the fits.

Parameterisation of the non-top-quark background
The W+jets events are expected to provide the largest contribution to the residual background. The '2J0T' sample is mostly populated by such events and contains a large number of events, making it in principle a suitable control region to study the expected contribution of W+jets events to the background in the signal region. However, the simulation shows that the reconstructed invariant mass distribution for W+jets events in the '2J0T' sample differs from that of the '2J1T' sample. Thus, simulation has been used for the characterisation of the W+jets component, as well as for all other nontop-quark background contributions. The shape of the invariant mass distribution for the sum of all non-top-quark background sources is well reproduced by a Novosibirsk function [37], with parameters μ and σ representing the Gaussian core, and τ describing the skewness of the distribution. The option to use the full simulated sample before the final selection, as is done for events containing top quarks, has not been chosen, as the parameters of the fitted function vary significantly with the requirement on |η j |, as shown in Fig. 6. Therefore, the sample obtained after applying the final selection is used to determine the shape parameters in the final fit.

Determination of the top quark mass from the fit
The invariant mass distribution of the selected top quark candidates is fitted with three components corresponding to signal, tt and non-top-quark processes, using the probability density functions described above. The mass is obtained from the resulting value of the mean of the Gaussian core of the Crystal Ball function fitting the single top quark contribu- tion, applying the calibration procedure described above. All parameters of the single top quark component are left free in the fit. The difference between the peak position of the tchannel and tt components is kept fixed to the value measured in simulation, to reduce the statistical fluctuations due to the small number of residual tt events. All remaining parameters (including normalisations) are fixed to the values extracted from simulation, after applying the final event selection.
The results of the fits to the simulated sample and to the collision data sample are shown in Fig. 7. The number of tchannel events returned by the fit is N fit t-ch = 2188 ± 72, in agreement with the number expected from simulation, N MC t-ch = 2216 +94 −78 . A value of m t = 172.95 ± 0.77 (stat) GeV is obtained after applying the mass calibration (Fig. 5). A systematic uncertainty of 0.39 GeV is associated to the mass calibration procedure.

Cross-checks
The consistency and stability of the fit are assessed using pseudo-experiments. Ensembles of experiments are simu- In each pseudo-experiment, the same fit described above is repeated and the top quark mass and the signal yield are derived. The resulting distributions of the top quark mass and its root-mean-square show that the fit does not have any significant bias, with the difference between fitted and generated top quark masses, normalised to the fitted mass uncertainty ("pull"), distributed as expected. Additionally, a test has been made where both the single top quark contribution and the tt components are fitted with a single Crystal Ball function. The results do not change appreciably within the present uncertainties with respect to the nominal fit. The mass measurement for the single top quark contribution is derived after having removed the single top antiquark events. As a check, the analysis has been repeated and the top quark mass has been measured using single top antiquark events. The difference between the two measurements is 0.8 ± 1.2 GeV, with a difference of −0.6 ± 1.5 GeV expected from simulation. Furthermore, the fit has been performed by simultaneously fitting single top quark and single top antiquark candidates: the fitted mass does not statistically differ with respect to the result obtained with the nominal fit. These studies confirm that the selection of only the top quark candidates does not introduce any bias in the measured top quark mass. In the following, the sources of uncertainties identified as relevant for the measurement are described, as well as the procedure adopted to evaluate their impact. All the uncertainties are then combined in quadrature to derive the total systematic uncertainty.
Jet energy scale (JES): JES factors are applied to the jet energy response in simulation to match that observed in data. The JES uncertainties are p T -and η-dependent, and are taken into account by scaling the energies of all jets up and down according to their individual uncertainties, as determined in dedicated studies [31]. The scaling is then propagated to the calculation of p miss T , and all other quantities dependent on the jet energies. The mass fit is repeated on the 'scaled' simulated sample and the shift with respect to the nominal fit is taken as a measure of the uncertainty. The uncertainties in the JES are subdivided into independent sources and grouped into different categories following the prescription defined in Ref. [39], aimed at simplifying the combination of measurements reported by the different LHC experiments. A total of five categories are identified referring to the effect of uncertainties related to the absolute scale determination using Drell-Yan events ("in-situ correlation group"), relative (η-dependent) calibration, and high-and lowp T extrapolation ("inter-calibration group"), flavour-specific corrections ("flavour-correlation group"), pileup corrections using an offset dependence on the jet p T ("pileup p T uncertainty"), and remaining sources, uncorrelated between ATLAS and CMS ("uncorrelated group"). b quark hadronisation model: This is the term that accounts for the flavour-dependent uncertainties arising from the simulation of the parton fragmentation. The total uncertainty can be decomposed into two separate contributions: the b quark fragmentation uncertainty and the uncertainty from b hadron decays. The b quark fragmentation uncertainty has been derived in the same way as in the top quark mass measurement using semileptonic tt events [38]. The Bowler-Lund fragmentation function for b hadrons is retuned to agree with the x B data measured by the ALEPH [40] and DEL-PHI [41] Collaborations, where x B represents the fraction of the b quark energy retained by a b hadron. A weight is attributed to each event, according to the x B value, and the difference with respect to the nominal setup is taken as the systematic uncertainty. The systematic uncertainty from the semileptonic branching ratio of b hadrons is taken from Ref.
[38], in which the branching fractions were varied by −0.45 and +0.77% to give the possible range of values and the associated uncertainty. Jet energy resolution (JER): After correcting for the mismatch between the data and simulation for the energy resolution, the uncertainty is determined by varying the corrected JER within its η-dependent ±1 standard deviation uncertainties. These changes are then propagated to the calculation of p miss T . Muon momentum scale: This contribution is determined by varying the reconstructed muon momenta by their uncertainties. These are determined as a function of the muon η and p T with a "tag-and-probe" method based on Drell-Yan data, as described in Ref. [42]. Unclustered missing transverse momentum: The uncertainty arising from the component of the missing transverse momentum that is not due to particles reconstructed as leptons and photons or clustered in jets ("unclustered p miss T ") is determined by varying it by ±10%.
Pileup: This is the uncertainty coming from the modelling of the pileup in data. It is taken as the sum of the effect of the uncertainty in the pileup rate (evaluated with pseudo-experiments in which the effective inelastic pp cross section has been varied by ±5%) and the pileup term extracted from the JES 'uncorrelated' group (see above). b tagging efficiency: To calculate the uncertainties from the b tagging efficiency and the misidentification rate, the p T -and η-dependent b tagging and misidentification scale factors are varied within their uncertainties for heavy-and light-flavour jets, as estimated from control samples [32]. The resulting changes are propagated to the event weights applied to the simulated events to obtain the uncertainties. Fit calibration: The mass is derived from the value of μ returned by the fit according to the mass calibration procedure described before: the same procedure provides an associated systematic uncertainty (Fig. 5, lower).
Background estimate: This is the uncertainty resulting from the use of simulated backgrounds in the mass determination. One contribution to the systematic uncertainty is determined by varying the background normalisations by ±1 times their standard deviation uncertainties. In addition, in the fit, the shape parameters of both the tt and the W+jets components are fixed: these parameters are varied by ±1 times their standard deviation uncertainties. As there are theoretical uncertainties on the inputs to the simulation which may alter the background shapes used in the mass fit, an additional 'radiation and matrix-element to parton-shower matching' uncertainty is included, as described below.
Generator model: Depending on whether the b quarks are considered part of the proton or not, the production of single top quarks can be studied in the 5-or 4-flavour schemes [10], respectively. The signal sample used in this work is produced with the 5-flavour scheme, where the b quarks are considered as constituents of the proton. To estimate the systematic uncertainty due to treating the b quarks like the lighter quarks, a comparison between a 5and a 4-flavour-scheme (2 → 2 and 2 → 3, respectively) samples has been performed: in the latter, the b quarks are generated in the hard scattering from gluon splitting.  The resulting difference is in agreement with what is obtained for the JES uncertainty. Radiation and matrix-element to parton-shower matching: This is the category which covers the QCD factorisation and renormalisation scale (with the nominal values of μ R = μ F = Q 2 ) and initial-and final-state radiation uncertainties.
For the renormalisation and factorisation scale uncertainty determination, dedicated samples with μ R and μ F scales shifted up or down by a factor of 2 are used. The uncertainty is determined by comparing the central result with the shifted ones.
For the matrix-element to parton-shower matching thresholds, tt and W+jets samples in which the thresholds have been shifted up or down by a factor of 2 are used, with the systematic uncertainty evaluated in the same way as for the scale uncertainty. This is not relevant for the signal data set, which does not have a matrix-element to parton-shower matching. This procedure covers initialand final-state radiation uncertainties. All variations are applied independently of each other and the corresponding uncertainties are treated as uncorrelated.
Underlying event: This term represents the uncertainty coming from the modelling of the underlying event (UE), the particles from the interaction that do not enter into the hard parton-parton interaction. It is evaluated by comparing the results from two different tunes of pythia, the default Z2* tune and the Perugia tune [46]. The differences in the value of the fitted mass are within the statistical uncertainty determined by the size of the simulated samples. In fact, the two opposite variations result in mass shifts with the same sign. For this reason, the uncertainty from this source is estimated as the maximum statistical uncertainty of the changes. Colour reconnection: This uncertainty is evaluated by comparing two different UE tunes in which one has the nominal colour-reconnection effects and the other has these turned off. Parton distribution functions: The PDF4LHC [14] prescriptions are followed to calculate the uncertainty due to the choice of the PDFs. The variation of the fitted top quark mass is estimated by using alternative sets of PDFs with respect to the nominal one, namely the MSTW2008CP [17], CT10 [24], and NNPDF2.3 [15] sets.
The systematic uncertainties that affect the top quark mass measurement are summarised in Table 1. Sources of systematic uncertainties that are totally or partially uncorrelated with the top quark mass measurements using tt events are the fit calibration, the background estimate, the generator model and theoretical parameters for the simulation of signal events, and the colour-reconnection effects.

Summary
The top quark mass is measured in a sample enriched in events with a single top quark for the first time. Top quarks are reconstructed from decays to a W boson and a b quark, with the W boson decaying to a muon and a neutrino. In the final sample, events with a top quark from single production in the t-channel account for 73% of the total number of events with a top quark. The measurement is obtained from a fit to the distribution of the reconstructed mass of top quark candidates, where the t-channel single top quark component is modelled separately from the contribution of other top quark production channels. The measured value is m t = 172.95 ± 0.77 (stat) +0.97 −0.93 (syst) GeV. This is in agreement with the current combination of Tevatron and LHC results, 173.34±0.27 (stat)±0.71 (syst) GeV, which is based on measurements using tt events. Because many of the systematic uncertainties affecting the measurement of m t using single top quark t-channel events are totally or partially uncorrelated with the measurements using tt events, and in addition the data sample analysed is largely statistically independent of the samples previously used, the result presented in this paper will be useful in the determination of world averages of the top quark mass. PhysRevD