Measurement of the t-channel single-top-quark production cross section and of the |Vtb| CKM matrix element in pp collisions at sqrt(s) = 8 TeV

Measurements are presented of the t-channel single-top-quark production cross section in proton-proton collisions at sqrt(s) = 8 TeV. The results are based on a data sample corresponding to an integrated luminosity of 19.7 inverse femtobarns recorded with the CMS detector at the LHC. The cross section is measured inclusively, as well as separately for top (t) and antitop (t-bar), in final states with a muon or an electron. The measured inclusive t-channel cross section is sigma[t-ch] = 83.6 +/- 2.3 (stat.) +/- 7.4 (syst.) pb. The single t and t-bar cross sections are measured to be sigma[t-ch,t] = 53.8 +/- 1.5 (stat.) +/- 4.4 (syst.) pb and sigma[t-ch,t-bar] = 27.6 +/- 1.3 (stat.) +/- 3.7 (syst.) pb, respectively. The measured ratio of cross sections is R[t-ch] = sigma[t-ch,t]/sigma[t-ch,t-bar] = 1.95 +/- 0.10 (stat.) +/- 0.19 (syst.), in agreement with the standard model prediction. The modulus of the Cabibbo-Kobayashi-Maskawa matrix element Vtb is extracted and, in combination with a previous CMS result at sqrt(s) = 7 TeV, a value abs(Vtb) = 0.998 +/- 0.038 (exp.) +/- 0.016 (theo.) is obtained.


Introduction
In the standard model (SM), the production of single top quarks (t or its antiparticle t) in protonproton (pp) collisions proceeds through the charged-current electroweak interaction.At leading order (LO), three different mechanisms can be distinguished, namely the t-channel, the s-channel and the associated production of a single top quark and a W boson (tW) [1][2][3].In this work, measurements are presented of t-channel production.LO diagrams contributing to t-channel single t and t production are presented in figure 1. Processes involving single top quarks provide direct probes of electroweak interactions, and thereby important tests of the SM predictions as well as excellent opportunities for searching for new physics.Since a Wtb vertex, where W and b denote the W boson and the b quark respectively, is involved in all SM single-top-quark production mechanisms, the modulus of the Cabibbo-Kobayashi-Maskawa (CKM) matrix element |V tb | can be determined from their measured cross sections.Depending on whether the b quarks are considered part of the proton or not, single-top-quark production can be studied in the 5-or 4-flavour schemes [4], respectively.In the 4-flavour (4F) scheme, the b quarks are generated in the hard scattering from gluon splitting.In the 5-flavour (5F) scheme, the b quarks are considered as constituents of the proton.An additional feature of the tand schannels, specific to pp collisions, is the difference between production cross sections of single t and t that results from a difference in parton distribution functions (PDF) of incident up and down quarks involved in the hard scattering.The ratio of t over t production cross sections in the t-channel (R t-ch. ) is therefore sensitive to the PDF of the up-and down-type quarks in the proton.The ratio R t-ch. is also directly sensitive to physics beyond the SM manifesting as anomalous couplings in the Wtb vertex [5], or to possible contributions from flavour-changing neutral current processes [6].For pp collisions at a centre-of-mass energy √ s = 8 TeV, the predicted theoretical cross section for SM t-channel single-top-quark production is σ theo.
For this calculation the top-quark mass m t is set to 173 GeV, and the factorisation and renormalisation scales are set both to m t .The uncertainty receives contributions from the PDF uncertainty and the missing higher-order corrections, estimated by varying the factorisation and renormalisation scales by a multiplicative factor of 0.5 or 2.0.The same calculations predict the following production cross sections for single t and t, separately: σ theo.t-ch.(t) = 56.4+2.1 −0.3 (scale) ± 1.1 (PDF) pb, σ theo.
The measurement performed by CMS of inclusive single-top-quark production cross section in the t-channel at √ s = 8 TeV, as well as separate measurements of single t and t production cross sections are presented.Signal events are characterised by products of top-quark decay that are accompanied by a light quark emitted at high rapidity and a soft b quark.Events are selected if a muon or electron consistent with originating from a top-quark decay chain is present in the final state.The signal yield is extracted from a maximum likelihood fit to the distribution of the absolute value of the pseudorapidity (η) of the jet (j ) originating from the light quark, |η j |.The expected distributions of |η j | are determined from data for the relevant backgrounds.Two independent fit procedures are implemented to extract the total t-channel production cross section and t and t production cross sections separately.The ratio of t-channel production cross sections at √ s = 8 TeV and 7 TeV, R 8/7 , can provide complementary information on the PDF with respect to the ratio of t over t, and can be compared to the prediction obtained using the cross sections in ref. [7], which is: R theo.

The CMS detector
The CMS apparatus features a 6 m diameter superconducting solenoid that provides a magnetic field of 3.8 T and allows for the relatively compact design of the detector.The inner bore of the solenoid hosts a tracking system, composed of silicon pixel and silicon strip detectors, that allows for reconstruction of charged-particle tracks bending in the internal magnetic field.A lead tungstate crystal electromagnetic calorimeter and a brass/scintillator hadron calorimeter surround the tracker volume.Outside the solenoid, gas-ionisation detectors, i.e. resistive plate chambers, drift tubes and cathode strip chambers, are interleaved with the steel fluxreturn yoke of CMS and form the muon system.A quartz-fibre and steel absorber Cherenkov calorimeter, located outside the muon system close to the beam pipe, extends the calorimetric system angular acceptance in the region along the beam axis.A more detailed description of the CMS detector can be found in ref. [16].
The CMS experiment uses a right-handed coordinate system centred on the nominal interaction point, with the z-axis parallel to the anticlockwise-beam direction, the x axis lying in the plane of the LHC ring and pointing to its centre, and the y axis pointing upwards to the surface.The pseudorapidity η is defined as − ln [tan (θ/2)], where θ is the azimuthal angle with respect to the z axis.

Data and simulated samples
The measurement is performed on a data sample collected during 2012 at √ s = 8 TeV, selected with triggers requiring one muon (µ) or one electron (e), and corresponding to an integrated luminosity of 19.7 fb −1 .
The simulated t-channel events are generated with POWHEG 1.0 [17][18][19][20] interfaced to PYTHIA 6.4 [21] for parton shower evolution and hadronisation.Other single-top-quark processes, i.e. the s-channel and the tW, are considered as backgrounds for this measurement and simulated with the same Monte Carlo (MC) generators.Top quark pair production, single vector boson production associated with jets (W/Z+jets), and double vector boson (diboson) production are amongst the backgrounds taken into consideration and have been simulated with MAD-GRAPH 5.148 [22] interfaced to PYTHIA for parton showering.The PYTHIA generator is used to simulate QCD multijet samples enriched with isolated muons or electrons.The value of the top-quark mass used in all simulated samples is m t =172.5 GeV.All samples are generated using the CTEQ6.6M[23] PDF set.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.The production cross section used to scale the simulation of single-top-quark tW and s-channel processes is taken from ref. [7], while the tt production cross section is taken from ref. [24].

Event selection and reconstruction
The signal events are defined by the decay of t → Wb → b ν, where = µ, e.The t → Wb → bτν contribute to the signal when a τ decays leptonically.The resulting final state includes a muon or electron, and escaping neutrinos (ν) that cause an imbalance in the momentum measured in the transverse plane.A bottom (or b) jet that stems from the hadronisation of the b quark from the top-quark decay accompanies the leptons.An additional jet originates from the light-flavoured quark recoiling against the top quark.The splitting of the gluon from the initial state produces a second b quark that recoils against the top quark, as shown in figure 1.The b jet from gluon splitting has generally a softer transverse momentum (p T ) spectrum and a broader |η| distribution compared to the one produced in top-quark decay, thus the acceptance for events with two b jets reconstructed in the final state is relatively small.In fact we can anticipate that using the selection described in this section, the number of signal events with two b jets reconstructed in the detector is one order of magnitude smaller than the number of events with just one b jet.
Events are selected online by the high-level trigger system requiring the presence of either one isolated muon with p T > 24 GeV and pseudorapidity |η| < 2.1, or one isolated electron with p T > 27 GeV and |η| < 2.5.The event is required to have at least one primary vertex reconstructed from at least four tracks, with a distance from the nominal beam-interaction point of less than 24 cm along the z axis and less than 2 cm in the transverse plane.When more than one primary vertex is found, the one with the largest ∑ p 2 T is used as leading vertex.All particles are reconstructed and identified with the CMS particle-flow (PF) algorithm [25,26].Events with exactly one good muon or electron candidate are accepted for analysis.Good muon candidates must have p T > 26 GeV and |η| < 2.1, while electron candidates must have p T > 30 GeV and |η| < 2.5, excluding the barrel-endcap transition region 1.44 < |η| < 1.57 because the reconstruction of an electron in this region is not optimal.The p T requirements on the leptons ensure that selected muons and electrons are in the plateau region of the respective trigger turn-on curves.
Muon isolation is ensured by applying requirements on the variable I rel , defined as the sum of the transverse energies deposited by stable charged hadrons, photons, and neutral hadrons in a cone of size ∆R = √ (∆η) 2 + (∆φ) 2 = 0.4, (where φ is the polar angle in radians) corrected by the average contribution of neutral particles from overlapping pp interactions (pileup), and divided by the muon p T .Muons are required to have I rel < 0.12.Electron isolation criteria are based on a variable defined analogously to the muons, with an isolation cone of size ∆R = 0.3.The isolation requirement for electrons is I rel < 0.1.Events are rejected if an additional muon (electron) candidate is present, passing looser selection requirements of p T > 10 (20) GeV, |η| < 2.5 (including the barrel-endcap transition region for electrons), and I rel < 0.2 (0.15).
The missing transverse momentum vector p / T is defined as the negative vector sum of the transverse momenta of all reconstructed particles.The missing transverse energy E T / is defined as the magnitude of p / T .The transverse mass m T for events with a muon is calculated as where p µ x and p µ y are the component of the muon momentum along the x and y axes, and p / x and p / y are the components of p / T along the x and y axes.In order to reduce the QCD multijet background, a requirement of m T > 50 GeV is applied for the muon decay channel, while a requirement of E T / > 45 GeV is applied instead for the electron channel.Control region studies, described in section 5.1, show that the procedure for the QCD multijet extraction in the electron channel yields a considerably smaller uncertainty when applying the requirement on E T / rather than on m T .
Jets are defined by clustering reconstructed particles with the anti-k T algorithm [27] with a distance parameter of 0.5.Charged particles are excluded if they have a distance with respect to any primary vertex along the z axis smaller than that with respect to the leading vertex.The average energy density in η-φ space of neutral particles not clustered into jets is used to extrapolate the energy due to pileup interactions in the jet cone.The jet energy is corrected accordingly.Further jet energy corrections are derived from the study of dijet events and pho-ton+jets events (see ref. [28]).Jets are required to be within |η| < 4.5 and to have a transverse energy E T > 40 GeV.In order to identify b-quark-induced jets, a b-tagging algorithm is used exploiting the 3D impact parameter of the tracks in the jet to define a "b-discriminator" [29].An optimised threshold is chosen on this variable with probability to misidentify jets coming from the hadronisation of light quarks (u, d, s) or gluons of 0.3% and an efficiency of selecting jets coming from b quarks of 46%, determined from simulation.Jets passing the chosen threshold are considered as "b-tagged".The majority of the background events surviving the final selection contain an actual b jet, the main exception being W+c-jet events.This algorithm is found to have good discriminating power with respect to this particular background.
Events are divided into categories according to the number of jets and b-tagged jets using the wording "n-jet m-tag ", referring to events with n jets, m of which are b-tagged.Once the event is been assigned to a category, a further selection based on the jet shape is performed to reduce the contamination due to jets coming from pileup interactions: the distance in the η-φ plane between the momenta of the particles constituting the jet and the jet axis is evaluated and its root mean square (RMS) over all the jet constituents is required to be smaller than 0.025.This requirement is applied on the jets that are not classified as b-tagged, and the event is rejected if either of those jets does not satisfy it.This requirement allows us to discriminate jets coming from u and d quarks with respect to jets coming from gluons or b quarks, which present a broader jet profile.This criteria is particularly useful in the forward region of the detector where other quality criteria making use of the tracking system cannot be applied.A top-quark candidate is reconstructed in the 2-jet 1-tag, in the 3-jet 1-tag and in the 3-jet 2-tag samples from a lepton, E T / and one b jet combination with the algorithm described in ref. [11].The b jet with the highest value of the b-discriminator is used for top-quark reconstruction in the 3-jet 2-tag sample.The mass of such a candidate "m νb " is used to define a signal region (SR) and a sideband region (SB) in each of those samples, selecting events respectively inside and outside the reconstructed top-quark mass window of 130 < m νb < 220 GeV.The variable η j in the 2-jet 0-tag sample is defined taking the pseudorapidity of each of the two jets, and two entries per event are present.In the 3-jet 1-tag sample η j is defined as the pseudorapidity of the jet with the smallest b-discriminator value.In the 2-jet 1-tag and in the 3-jet 2-tag samples it is defined as the pseudorapidity of the non-b-tagged jet.The category enriched with t-channel signal is the one with 2 jets and 1 tag.The final procedure to isolate the signal from background uses the absolute value |η j |.The pseudorapidity distribution of the outgoing jet j is typical of the t-channel processes where a light parton recoils against a much more massive particle like the top quark.Signal events populate forward regions in the |η j | spectrum that correspond to the tails of the |η j | distribution for SM processes.
The total event yields in the signal and sideband regions of the 2-jet 1-tag sample for muons and electrons are reported in table 1.The event yields in the signal region for positively and negatively charged muons and electrons separately are reported in table 2.
Table 1: Event yield for the main processes in the 2-jet 1-tag signal region (SR) and sideband region (SB), for the muon and electron decay channels.Expected yields are taken from simulation and their uncertainties are due to the finite size of the MC sample with the exception of QCD multijet yield (see section 5.1), and W/Z+jets yield (see section 5.3), whose yields and uncertainties are taken as the statistical component of the uncertainty in the estimation from data.

Background estimation and control samples
The physics processes that constitute the main backgrounds to single-top-quark production in the t-channel are tt, W+jets, and QCD multijet production.Control samples are defined for each of these contributions in order to check that the variables used in the analysis are reproduced correctly in the simulations.For the main backgrounds the most important distributions, together with constraints on their production rates, are derived from data making use of these control samples.

QCD multijet background
The vast majority of QCD multijet events are successfully rejected applying the selection described in section 4. The selected multijet events are thus found to be rare occurrences in the respective distributions, for instance populating the tails of the typical multijet lepton-p T spectra.The modelling uncertainties on the simulation have greater impact in those regions.We thus estimate the QCD multijet contribution in our signal and sideband regions directly from data, in the 2-jet 1-tag category as well as in the other control samples.The measurement is performed with a fit to the distribution of the transverse mass in the muon decay channel, and to the distribution of the missing transverse energy in the electron decay channel.A maximumlikelihood fit to the distribution either of m T in the muon case, or E T / in the electron case is performed.The data are parametrised as: The variable x is m T for the muon decay channel and E T / for the electron decay channel, while S (x) and B (x) are the expected distributions for the sum of all processes with a W boson and QCD multijet events, respectively.The distribution S (x) is derived from simulation and it includes the contribution from the signal.The distribution B (x) is obtained from a QCD multijet enriched data sample defined by taking muons and electrons with the same criteria as defined in section 4, but with reversed isolation requirements for both leptons, selecting muons or electrons with I rel > 0.2 or 0.15 respectively.The data samples defined in this way contain a fraction of events originating from QCD multijet processes of 98% in the case of the muon decay channel and of more that 99% for the electron decay channel.The residual contribution from other non-QCD multijet processes is subtracted from these samples using the expectation from simulation.The fit procedure is repeated using different QCD multijet models, obtained by either varying the isolation requirement that defines the control region or using the simulation for the QCD multijet distribution.The kinematic bias on the multijet m T (E T / ) distributions due to the extraction from the control sample is covered by the systematic uncertainty defined this way.

Top quark pair background
The tt process dominates in events with larger jet and b-tag multiplicity than the 2-jet 1-tag sample used for signal extraction.Two control samples enriched in tt are thus defined, labelled 3-jet 1-tag and 3-jet 2-tag.The distribution of |η j | in the 3-jet 1-tag and in the 3-jet 2-tag samples is shown in figure 2. Good agreement between data and simulation in the two control samples is displayed, giving confidence in the simulation of the kinematic properties of the tt background.The lepton charge in the 3-jet 1-tag and 3-jet 2-tag samples is shown in figure 3. The corresponding charge ratio in the two samples is shown in figure 4, and is close to unity as expected for tt enriched samples.To reduce the dependence of the measurements on the modelling of tt processes, the |η j | distribution (template) used for signal extraction is modified taking into account the |η j | distribution of the non-b-tagged jet in the 3-jet 2-tag sample as follows.The contribution of all SM processes except for tt in the 3-jet 2-tag is subtracted from the template |η| distribution of the non-b-tagged jet taken from data.Then the bin-by-bin ratio of the resulting template distribution and the corresponding distribution from the tt process is taken as the |η j |-dependent correction factor for the tt in the 2-jet 1-tag sample.This ratio is then applied to the simulated distribution of |η j | in the SR and SB.

The W/Z+jets background
The 2-jet 0-tag sample is enriched with W/Z+jets background and it is used to test the agreement between simulation and data on the distributions used for the signal extraction procedures.The distribution of |η j | in the 2-jet 0-tag is shown in figure 5, and good agreement between data and simulation is displayed.The lepton charge in the 2-jet 0-tag sample is shown in figure 6.The characteristic imbalance in the production of positively and negatively charged leptons in W+jets events can be seen clearly in the data, and the corresponding charge ratio is shown in figure 7. The jets in this sample mostly originate from light quarks (u, d, s) or gluons, which tend to behave differently from heavy-flavour jets (stemming from c and b quarks).For this reason, in the final fit procedure described later on in section 6 the W+jets charge ratio is extracted from data as well.The SB region in the 2-jet 1-tag sample is used in order to estimate the W/Z+jets component in a region that is expected to have a similar composition in terms of W/Z+heavy flavours with respect to the sample that is used for the cross section extraction, i.e. the 2-jet 1-tag SR.The |η j | distribution for W/Z+jets processes is taken from the sideband region by subtracting all other processes bin by bin.For this subtraction all samples except for tt and QCD multijet are derived from simulation.The latter two are estimated with the techniques described above.The scale factors between sideband region and signal region are derived from simulation.This procedure is performed for the inclusive distribution, as well as for positively and negatively charged leptons separately.The bias due to the different kinematic properties of the two regions is estimated on simulations and removed, and the uncertainty on the composition in terms of W+c-jets and W+b-jets events is taken into account as described in section 7.

Signal extraction and cross section measurement
Two binned maximum-likelihood fits to the |η j | distributions of the events in the 2-jet 1-tag SR are performed.The first fit extracts the inclusive single-top-quark cross section, the second extracts the separate single t and t cross sections.
The expected number of events in each |η j | bin is modelled with the following likelihood function: In addition to the signal (indicated with subscript s), three background components are consid-      Figure 7: Charge ratio between positively and negatively charged leptons in the 2-jet 0-tag sample for muon (left) and electron (right) decay channels.The charge ratio is shown separately for each process, as well as after normalising the sum of all predictions to the data yield.Systematic uncertainty bands include all uncertainties.
The inclusive cross section is extracted from events with positively or negatively charged leptons, defining one likelihood function per lepton flavour, as in equation 5, then fitting simultaneously the two distributions for muons and electrons.The single t and t cross sections are extracted by further dividing the events by lepton charge, defining one likelihood function per lepton flavour and per charge, as in equation 5, then fitting simultaneously the four distributions.
The definition of the probability distribution functions and of the parameters included in the fit are described in the following: • Signal: P s for both fits is taken from simulation (see also section 3) as the predicted |η j | distribution.The total yield N s is fitted unconstrained in the inclusive singletop-quark cross section fit.Two parameters are introduced in the single t and t cross section fit for the positively and negatively charged lepton signal yield and fitted unconstrained.
• EW component: W/Z + jets, diboson: The P EW distribution is taken as the sum of the contribution of W/Z+jets and diboson processes.The W/Z+jets normalisation and distribution are estimated from the m νb sideband with the method described in section 5.This sideband method is applied to both muons and electrons, inclusively with respect to the lepton charge in the case of the inclusive top-quark cross section fit, and separately for positively and negatively charged leptons in the case of the single t and t cross section fit.The diboson contribution is then taken from simulation.The two contributions are summed together and the total yield N EW is derived by the fit.To take into account the prior knowledge of the normalisation obtained from the sideband a Gaussian constraint is applied to N EW in the fit, i.e. the likelihood function is further multiplied by a Gaussian function of N EW .The mean value of this function is taken from the procedure previously described in this paragraph, while the standard deviation is taken equal to the difference between the data-based yield of W/Z+jets and the expectation from simulation in the sideband region.For the single t and t cross section ratio fit, the N EW are fitted separately for positively and negatively charged leptons.
• Top quark component: tt, tW and s-channel: P t is taken from the data-based pro- cedure described in section 5, to which the single-top-quark tW and s-channel processes are added with a normalisation factor taken from simulation.This contribution is separated by lepton flavour and charge assuming charge symmetry of tt and tW events.The s-channel charge ratio is fixed to the SM prediction.The yield N t is then fitted with a Gaussian constraint, centred on the value obtained from simulation and with a variation of ±10%, which is chosen to cover both experimental and theoretical uncertainties on the tt cross section.
• QCD multijet: P MJ is taken from the QCD multijet enriched sample defined in section 5, adding an extra requirement on the angular distance of the lepton and the jets, ∆R( , j) > 0.3.The yield is fixed to the results of the m T and E T / fit.
The fit strategy driving this parametrisation is focused on constraining from data the W/Z+jets and tt backgrounds.In the particular case of the single t and t cross section fit, the event ratio of positively and negatively charged W bosons is constrained as well.The cross sections are extracted using the detector acceptance derived from the simulated signal sample.The total cross section measurement from the inclusive analysis is more precise than the one inferred from the separate-by-charge fit, due to the additional uncertainty from the W charged ratio, which is extracted from data.The |η j | distributions for the muon and electron decay channels obtained by normalising the contribution of each process to the value of the inclusive cross section and t and t cross section ratio fits are shown in figures 8 and 9, respectively.An  t-channel properties in the signal sample after normalising each process to the fit results.The reconstructed top-quark mass m νb in the region with |η j | > 2.5, after scaling each process contribution to the normalisation obtained from the fit, is shown in figure 10.This region is expected to be depleted of background events and enriched in t-channel signal events, hence displaying a characteristic peak around the top-quark mass value, which appears clearly in data for both the muon and the electron channels.

Systematic uncertainties
Contributions to the total systematic uncertainty are evaluated, with the exception of the uncertainties on the background estimation described in section 5 and on the simulated samples size, with the following procedure: pseudo-experiments are constructed using for each process the distributions and the yields generated considering the altered scenario.A fit to the |η j | distribution is then performed for each pseudo-experiment with the nominal setup, and the mean shift of the fit results with respect to the value obtained for the nominal fit is taken as the corresponding uncertainty.A detailed description of each source of systematic uncertainty and of the treatment of uncertainties related to the data-based background estimation and to the size of simulated samples follows: • Jet energy scale (JES), jet energy resolution (JER), and missing transverse energy: All reconstructed jet four-momenta in simulated events are simultaneously varied according to the η-and p T -dependent uncertainties in the jet energy scale and resolution.The variation of jet momenta causes the total momentum in the transverse plane to change, thus affecting the E T / as well.The component of the missing transverse energy that is not due to particles reconstructed as leptons and photons or clustered in jets ("unclustered E T / ") is varied by ±10% [28].
• Pileup: The uncertainty in the average expected number of additional interactions per bunch crossing (±5%) is propagated as a systematic uncertainty to this measurement.
• B-tagging: B-tagging and misidentification (mis-tag) efficiencies are estimated from control samples [29].Scale factors are applied to simulated samples to reproduce efficiencies in data and the corresponding uncertainties are propagated as systematic uncertainties.
• Muon/electron trigger and reconstruction: Single-muon and single-electron trigger efficiency and reconstruction efficiency as a function of the lepton η and p T are estimated with a "tag-and-probe" method based on Drell-Yan data, as described in ref. [30].The effect of the incorrect determination of the muon charge is negligible, while for electrons the uncertainty on the determination of the charge has been measured at √ s = 7 TeV in ref. [31].• W+jets, tt, and QCD multijet estimation: The distributions and normalisations of these three main backgrounds are derived mostly from data as described in section 5.
The uncertainty related to the W+jets and tt estimation is evaluated by generating pseudo-experiments in the SB and in the 3-jet 2-tag sample.The background estimation is repeated, and then the fit to |η j | is performed and the uncertainty is taken as the RMS of the distribution of fit results.An uncertainty in the W+jets contribution is obtained from alternative |η j | shapes derived from simulation by varying the W+b-jets and the W+c-jets background fractions by ±30% independently in the SR and SB regions.An additional uncertainty in the tt estimation procedure is determined by performing the signal extraction using the tt distribution in the entire m νb range, then using two different distributions for the signal and background regions.The difference of the two results is taken as the uncertainty.The QCD multijet normalisation is varied by ±50% independently for muon and electron decay channels.This variation range is obtained by performing the multijet estimation under different conditions and assumptions as described in section 5, and taking the maximum difference with respect to the value obtained with the nominal estimation procedure.Additionally, all other systematic uncertainties are coherently propagated through • Background normalisation: An uncertainty in the tt normalisation of ±10% is considered, covering the difference between theoretical predictions in [7] and [32].For dibosons and single-top-quark tW and s-channel production the assumed uncertainty is ±30%, motivated by refs.[7,33].• Signal modelling: Renormalisation and factorisation scales used in the signal simulation are varied by a factor 2 up and down, and the corresponding variation is considered as the systematic uncertainty.The uncertainty on the simulation is obtained by comparing the results obtained with the nominal POWHEG signal samples with the ones obtained using samples generated by COMPHEP [34,35].Half of the difference is taken as systematic uncertainty.
• PDFs: The uncertainty due to the choice of the PDF set is estimated by reweighting the simulated events and repeating the signal extraction procedure.The envelope of the CT10 [36], MSTW [8], and NNPDF [37] PDF sets is taken as uncertainty, according to the PDF4LHC recommendations [38].
• Simulation sample size: The statistical uncertainty due to the limited size of simulated samples is taken into account by generating pseudo-experiments reproducing the statistical fluctuations of the model.The fit procedure is repeated for each pseudo-experiment and the uncertainty is evaluated as the RMS of the distribution of fit results.
• Luminosity: The integrated luminosity is known with a relative uncertainty of ±2.6% [39].
The contribution of each source of uncertainty to the cross section and their ratio measurements is shown in tables 3 and 4, respectively.Uncertainties due to the limited size of simulated and control samples in data for the background estimation do not cancel and thus have an impact on the ratio measurement larger than on the total cross section.Uncertainties that affect the signal efficiency in a similar way for single t and t, such as the b-tagging, or the lepton trigger and reconstruction efficiencies, tend to cancel in the cross section ratio, thus have a smaller impact on its measurement.The luminosity uncertainty cancels as well in the ratio.Uncertainties that affect the background processes that are independent from the lepton charge, like the tt or the QCD multijet, have a bigger impact on the single t cross section, for which the signal-tobackground ratio is less favourable, and for this reason they do not cancel out entirely in the ratio measurement.Since single t and t production depend on different quark PDFs, the corresponding PDF uncertainties are largely anticorrelated, and the corresponding contribution is enhanced in the charge ratio measurement.As the momentum and pseudorapidity spectra of quarks and leptons for the single t and t processes are different, the modelling uncertainties and the uncertainties from the jet energy scale and missing transverse energy do not fully cancel out in the ratio measurement.
Because of these differences, the event yields returned by the inclusive single-top-quark cross section and the single t and t cross section fits are not numerically identical.A consequence of this is that the values for the total cross section obtained in the two fits differ.In particular the uncertainty in the heavy-flavour component is anticorrelated between the two measurements, and the theoretical uncertainties tend to affect the exclusive extraction more than the inclusive one.
The choice to keep two separate procedures is motivated by the fact that the inclusive fit has a better overall performance regarding the systematic uncertainties in the inclusive cross section measurement.
The measured single t and t production cross sections in the t-channel are A comparison of the currently available measurements of the inclusive cross section with the SM expectation obtained with a QCD computation at NLO with MCR in the 5F scheme [40] and at NLO+NNLL [41] is shown in figure 11.The measurement is compared to the previous CMS t-channel cross section measurement at √ s = 7 TeV [11] and the Betatron measurements at √ s = 1.96TeV [9,15].The measurements are compared with the QCD expectations computed at NLO with MCR in the 5F scheme and at NLO+NNLL.The error band (width of the curve) is obtained by varying the top-quark mass within its current uncertainty [42], estimating the PDF uncertainty according to the HEPDATA recommendations [43], and varying the factorisation and renormalisation scales coherently by a factor two up and down.The prediction in pp collisions can be also compared with the one at pp because the inclusive single-top-quark cross section does not depend on whether the light quark originates from a proton or from an antiproton.

Cross section ratios
The ratio of t-channel production cross sections at √ s = 8 and 7 TeV is derived with respect to the result reported in ref. [11] for the single-top-quark t-channel cross section at √ s = 7 TeV.Three measurements are combined in ref. [11]: two multivariate analyses and one, the η j analysis, making use of a strategy and a selection that are close to the ones reported in this paper.The correlations between the sources of uncertainties reported in section 7 and those in ref. [11] are determined in the following way: the uncertainties related to signal extraction and background estimation from data are treated as fully uncorrelated between 7 and 8 TeV, while for the rest of the uncertainties the 8 TeV analysis is considered fully correlated with respect to its 7 TeV η j counterpart, and the same choices for correlation as in [11] are adopted between the 8 TeV η j analysis and the two 7 TeV multivariate analyses.Taking into account the correlations as described, the measured ratio is R 8/7 = σ t-ch.(8 TeV)/σ t-ch.(7 TeV) = 1.24 ± 0.08 (stat.)± 0.12 (syst.).
The measured ratio of single t to t production cross sections at A comparison is shown in figure 12 of the measured R t-ch. to the predictions obtained with several PDF sets: MSTW2008NLO [8], HERAPDF1.5NLO [44], ABM11 [45], CT10, CT10w [36], and NNPDF [37].For MSTW2008NLO, NNPDF, ABM, and CT10w the fixed 4F scheme PDFs are used together with the POWHEG 4F scheme calculation.The POWHEG calculation in the 5F scheme is used for all other PDFs, as they are derived from a variable flavour scheme.The nominal value for the top-quark mass used is 173.0GeV.Error bars for the CMS measurement  include the statistical (light yellow) and systematic (dark green) components.Error bars for the different PDF sets include the statistical uncertainty, the uncertainty in the factorisation and renormalisation scales, derived varying both of them by a factor 1/2 and 2, and the uncertainty in the top-quark mass, derived varying the top-quark mass between 172.0 and 174.0 GeV.The different PDF sets predictions for this observable are not always compatible with each other within the respective uncertainties, thus displaying the potential for this measurement to discriminate between the different sets, should a better precision be achieved.

Extraction of |V tb |
A feature of t-channel single-top-quark production is the presence of a Wtb vertex.This allows for an interpretation of the cross section measurement in terms of the parameters regulating the strength of this coupling, most notably the CKM matrix element V tb .The presence of anomalous couplings at the Wtb vertex can produce anomalous form factors [46][47][48] which are parametrised as f Lv , where "Lv" refers to the specific left-handed vector nature of the couplings that would modify the interaction strength.In the approximation |V td |, |V ts | |V tb |, we consider the top-quark decay branching fraction into Wb, B, to be almost equal to 1, thus obtaining | f Lv V tb | = √ σ t-ch./σ theo.t-ch. .The choice of this approximation is motivated by the fact that several scenarios beyond the SM predict a deviation of the measured value of f Lv from 1, but only a mild modification of B [49].This allows to interpret a possible deviation from SM single-top-quark production cross section in terms of new physics.In the standard model case, f Lv = 1, implying that the cross section measurement yields a direct constraint on |V tb |.Thus inserting in the definition for | f Lv V tb | the measured cross section from equation 6 and the where both the experimental and the theoretical uncertainties are reported.The former comes from the uncertainties on the measurement of σ t-ch., while the latter comes from the uncertainties on σ theo.t-ch. .A similar measurement of | f Lv V tb | is performed in ref. [11].The results for | f Lv V tb | from this paper and from the three analyses in [11] are combined using the best linear unbiased estimator (BLUE) [50] method, considering the full correlation matrix amongst the four measurements and the correlations described for the R 8/7 measurement, obtaining the following result: | f Lv V tb | = 0.998 ± 0.038 (exp.) ± 0.016 (theo.)(7+8 TeV combination).
This result can be directly compared with the current world average of |V tb | from the Particle Data Group [51], which is performed without the unitarity constraints on the CKM matrix and, using the above formalism for non-SM contributions, yields | f Lv V tb | = 0.89 ± 0.07.From the result in equation 11, the confidence interval for |V tb |, assuming the constraints |V tb | ≤ 1 and f Lv = 1, is determined using the Feldman-Cousins unified approach [52], being |V tb | > 0.92 at the 95% confidence level.

Summary
The total cross sections for production in the t-channel of single top quarks and individual single t and t have been measured in proton-proton collisions at the LHC at √ s = 8 TeV.The inclusive single-top-quark t-channel cross section has been measured to be σ t-ch.= 83.6 ± 2.3 (stat.)± 7.4 (syst.)pb.The single t and t cross sections have been measured to be σ t-ch.(t) = 53.8± 1.5 (stat.)± 4.4 (syst.)pb and σ t-ch.(t) = 27.6 ± 1.3 (stat.)± 3.7 (syst.)pb, respectively.Their ratio has been found to be R t-ch.= 1.95 ± 0.10 (stat.)± 0.19 (syst.).The ratio of t-channel single-top-quark production cross sections at √ s = 8 and 7 TeV has been measured to be R 8/7 = 1.24 ± 0.08 (stat.)± 0.12 (syst.).These measurements are in agreement with the standard model predictions.From the measured single-top-quark production cross section, the modulus of the CKM matrix element V tb has been determined.This result has been combined with the previous CMS measurement at 7 TeV, yielding the most precise measurement of its kind up to date: |V tb | = 0.998 ± 0.038 (exp.) ± 0.016 (theo.).Assuming |V tb | ≤ 1, the 95% confidence level limit has been found to be |V tb | > 0.92.

Figure 1 :
Figure 1: Leading-order Feynman diagrams for (left) single t and (right) t production in the t-channel.

Figure 2 :
Figure 2: Distribution of |η j | in the 3-jet 1-tag (upper left, upper right), and 3-jet 2-tag (lower left, lower right) samples for muon and electron decay channels.The yield of the simulated processes is normalised to the results of the fit described in section 6. Systematic uncertainty bands include all uncertainties.

Figure 3 :
Figure 3: Charge of the lepton in the 3-jet 1-tag (upper left, upper right), 3-jet 2-tag (lower left, lower right) samples for muon and electron decay channels.The sum of all predictions is normalised to the data yield.Systematic uncertainty bands include all uncertainties on the charge ratio.

Figure 4 :
Figure 4: Charge ratio between positively and negatively charged leptons in the 3-jet 1-tag (upper left, upper right), 3-jet 2-tag (lower left, lower right) samples for muon and electron decay channels.The charge ratio is shown separately for each process, as well as after normalising the sum of all predictions to the data yield.Systematic uncertainty bands include all uncertainties.

Figure 5 :
Figure 5: Distribution of |η j | in the 2-jet 0-tag sample for muon (left) and electron (right) decay channels.The QCD multijet contribution is derived from the fit to m T and E T / .Systematic uncertainty bands include pre-fit uncertainties, both on the normalisation and on the shape of the distributions.ered: the electroweak component (with subscript EW composed of W/Z+jets and dibosons), the top quark component (with subscript t composed of tt and single-top-quark tW and schannel processes), and the QCD multijet component (with subscript MJ).In equation 5, N s , N EW , N t and N MJ are the yields of the signal and of the three background components; P s , P b (b=EW, t, MJ) are the binned probability distribution functions for the signal and for the different background components.

Figure 6 :
Figure 6: Charge of the lepton in the 2-jet 0-tag sample for muon (left) and electron (right) decay channels.The sum of all predictions is normalised to the data yield.Systematic uncertainty bands include all uncertainties on the charge ratio.

Figure 8 :
Figure 8: Fitted |η j | distributions for muon (left) and electron (right) decay channels, normalised to the yields obtained from the combined total cross section fit.Systematic uncertainty bands include the shape uncertainties on the distributions.indication of the validity of the fit extraction procedure comes from the study of characteristic

Figure 9 :
Figure 9: Fitted |η j | distributions for muon (upper left, lower left) and electron (upper right, lower right) decay channels, normalised to the yields obtained from the combined single t and t cross section ratio fit.Systematic uncertainty bands include the shape uncertainties on the distributions.

Figure 10 :
Figure10: Distribution of reconstructed top-quark mass m νb for muon (left) and electron (right) decay channels, in the region with |η j | > 2.5, the contribution of each process is scaled to the cross section derived from the fit.Systematic uncertainty bands include the shape uncertainties on the distributions and uncertainties on the normalisation in the |η j | > 2.5 region.

Figure 11 :
Figure 11: Single-top-quark production cross section in the t-channel versus collider centre-ofmass energy.

Table 2 :
Event yield for the main processes in the 2-jet 1-tag signal region, for events with positively and negatively charged muons and electrons.Expected yields are taken from simulation and their uncertainties are due to the finite size of the MC sample with the exception of QCD multijet yield (see section 5.1), and W/Z+jets yield (see section 5.3), whose yields and uncertainties are taken as the statistical component of the uncertainty in the estimation from data.

Table 3 :
Relative impact of systematic uncertainties for the combined muon and electron decay channels.