text:
# Measurement of the 𝑡𝑡̄ production cross section in 𝑝𝑝̄ collisions at √𝑠 = 1.8 TeV

## I Introduction

In the Standard Model (SM) the top quark completes the third fermion generation. The measurement of the top-quark pair production cross section $`\sigma _{t\overline{t}}`$ is of interest as a test of quantum chromodynamics (QCD) predictions. Aside from the obvious observation that a deviation from these predictions could be indicative of new physics, recent QCD calculations predict $`\sigma _{t\overline{t}}`$ with an uncertainty smaller than 15%  which motivates measurements of comparable precision.

In $`p\overline{p}`$ collisions at $`\sqrt{s}`$ = 1.8 TeV top quarks are pair produced through $`q\overline{q}`$ annihilation ($``$ 90%) or gluon fusion ($``$ 10%). In the SM framework the top quark decays into a $`\mathrm{W}`$ boson and a $`b`$ quark. When one of the $`\mathrm{W}`$ bosons decays to an electron or a muon, the final state includes a charged lepton with high transverse momentum ($`p_T`$), a large transverse energy imbalance from the undetected neutrino, referred to as $`\mathrm{}\mathrm{E}_\mathrm{T}`$, and four jets from the hadronized quarks. However, due to gluon radiation or jet merging, the number of detected jets may vary. We measure $`\sigma _{t\overline{t}}`$ using this final state, referred to in this paper as $`\mathrm{W}+`$ jets and selected from the data (105.1 pb-1) collected by CDF in the $`19921995`$ collider run.

The same data set has been used in the previous CDF measurement of $`\sigma _{t\overline{t}}`$ . This paper revises that measurement and expands on many of the analysis details. The selection of the $`\mathrm{W}+`$ jet sample follows the guide-lines used in all previous CDF measurements of the top mass and production cross section .

As done in previous analyses, we employ two techniques to enhance the relative fraction of events coming from top quark decays with respect to the background. The first method searches a jet for the presence of a secondary vertex reconstructed using the silicon vertex detector (SVX) and displaced from the primary event vertex due to the long $`b`$-quark lifetime (SECVTX tag). The second method searches a jet for the presence of a lepton, indicative of a semileptonic $`b`$-decay. Since these leptons typically have low momentum compared to the lepton from the $`W`$ decay, they are referred to as soft lepton tags (SLT).

In this analysis we use the same SECVTX and SLT algorithms as in Ref. . Differently from Refs. , we search jets and not events for soft lepton tags; this approach provides the rates of backgrounds for the top mass measurement which uses the same technique .

As a cross-check, we take advantage of a third algorithm, jet-probability, which uses the impact parameter significance of all tracks in a jet to derive a probability that the jet originates from the primary event vertex . Jets with small probability of having zero lifetime are considered jet-probability tags (JPB). The value of the jet-probability threshold is tuned to have higher tagging efficiency than SECVTX in jets originating from $`c`$ quarks and a higher rate of fake tags in jets without heavy flavor. Although similar to SECVTX in that it is based on track information from displaced vertices, it is used here to provide important cross-checks of the background determination and measured cross section.

The method used to measure $`\sigma _{t\overline{t}}`$ is outlined in Ref.  and has since been improved. As summarized in Ref. , the method relies on the calculation of all the background contributions to the tagged $`W`$+jet sample. The excess over background of the $`W+3,4`$ jet events with at least one tag is attributed to $`t\overline{t}`$ production and used to derive $`\sigma _{t\overline{t}}`$.

The major sources of background are the processes $`p\overline{p}Wg`$ with $`gb\overline{b}`$, $`c\overline{c}`$ (referred to as gluon splitting) and $`p\overline{p}Wc`$. The second largest source of background is mistags (tags in jets which do not contain heavy flavor). Smaller contributions come from other processes like non-$`W`$ production, single top production, $`WW`$, $`WZ`$, $`ZZ`$ and $`Z\tau \tau `$.

The method used to measure $`\sigma _{t\overline{t}}`$ relies on the correct calibration of the Monte Carlo generators and the detector simulation. Simulated events are produced with the Herwig  or Pythia  Monte Carlo generators. Hadrons with heavy flavor ($`b`$ and $`c`$) are decayed using the CLEO Monte Carlo (QQ) . All other particles are decayed, when appropriate, by the CDF detector simulation (QFL) which uses its own lifetime table for $`b`$ and $`c`$-hadrons. QFL simulates the interaction of all particles in the final state with the CDF detector; the detector response is based on parametrizations that are functions of the particle kinematics and have been derived using the data.

This paper describes the work done to understand and improve the calibrations used in the method to calculate the background to $`t\overline{t}`$ events using independent data samples and the corresponding simulations. This work was primarily focussed on the components with the largest influence on the determination of $`\sigma _{t\overline{t}}`$; mistags, the efficiencies of the tagging algorithms, and the fraction of $`W+`$ jet direct production which contains heavy flavor. We summarize here the relevant conclusions:

* we find that the SECVTX tagging efficiency for $`b`$-quark jets ($`b`$-jets) is ($`25\pm 13`$)% higher in data control samples than it is in the Monte Carlo simulation of the same processes. Therefore we conclude that the $`b`$-jet tagging rate in Refs.  is underestimated by this factor. This data-to-simulation discrepancy is largely due to errors in the simulation that were found a posteriori. Instead of remaking the large Monte Carlo samples used in this paper we chose to rescale the simulated $`b`$-quark tagging rate by the factor of 1.25$`\pm `$0.13.
* we find that the rate of SECVTX mistags in jets without heavy flavor is (50$`\pm `$5)% smaller than what we estimated in Refs. .
* we find the fraction of $`gb\overline{b}`$ and $`gc\overline{c}`$ in the $`\mathrm{W}+`$ jet direct production evaluated with Herwig needs to be increased by ($`39\pm 19`$)% and ($`35\pm 36`$)%, respectively.

These last two effects tend to cancel, leaving the net background to top approximately unchanged from our previous results in Refs. . In the $`W+3,4`$ jets sample we observe 29 events with one or more SECVTX tags and 25 events with one or more SLT tags. The expected backgrounds are $`8.0\pm 1.0`$ and $`13.2\pm 1.2`$ events, respectively. The excess of SECVTX tags yields the cross section $`\sigma _{t\overline{t}}`$ = $`5.08\pm 1.54`$ pb and the excess of SLT tags yields $`\sigma _{t\overline{t}}`$ = $`9.18\pm 4.26`$ pb for a top quark mass of $`175`$ GeV/$`c^2`$.

Following a brief description of the CDF detector in Section II, Section III describes the triggers and the reconstruction of leptons, jets and the missing transverse energy. The selection of $`\mathrm{W}+`$ jet events is detailed in Section IV, along with the selection of the $`Z+`$ jet sample, which will be used to check the background calculation. The selection of other data samples used to calibrate the event generators and the detector simulation is described in Section V. Sections VI and VII are dedicated to the Monte Carlo generators and the CDF detector simulation (QFL), respectively. Section VIII describes the algorithms used for the identification of the heavy flavor. The efficiency of those algorithms is calculated in Section IX, which also includes numerous checks of the result and an evaluation of its systematic uncertainty. The new method for evaluating mistags and the determination of its accuracy is described in Section X. Section XI details the calibration of the $`gb\overline{b}`$ and $`gc\overline{c}`$ cross sections predicted by the Herwig generator. Section XII describes the calculation of the backgrounds to the $`t\overline{t}`$ production. In Section XIII, we check the background calculation using the $`Z+`$ jet sample. Additional checks of the background calculation are described in Section XIV. Finally, $`\sigma _{t\overline{t}}`$ is derived in Section XV. In Section XVI, we combine the present results with previous CDF measurements of $`\sigma _{t\overline{t}}`$ that are derived using different data sets. We conclude in Section XVII.

## II  The CDF detector

CDF is a general purpose detector with azimuthal and forward-backward symmetry designed to study $`p\overline{p}`$ interactions. The CDF coordinate system has the $`z`$-axis pointing along the proton momentum and the $`x`$-axis located in the horizontal plane of the Tevatron storage ring pointing radially outward so that the $`y`$-axis points up. The coordinates $`r`$-$`\varphi `$ are the standard cylindrical coordinates. A complete description of CDF can be found in Refs. . The detector components most relevant to this analysis are summarized below.

A superconducting solenoid of length 4.8 m and radius 1.5 m generates a 1.4-T magnetic field. The solenoid contains three types of tracking chambers for detecting charged particles and measuring their momenta. A four layer silicon microstrip vertex detector (SVX) surrounds the beryllium beam pipe of radius 1.9 cm. The SVX has an active length of 51 cm; the four layers of the SVX are at distances of approximately 2.9, 4.2, 5.5 and 7.9 cm from the beamline. Axial micro-strips with 60-$`\mu `$m pitch provide accurate track reconstruction in the $`r`$-$`\varphi `$ plane transverse to the beam . Outside the SVX there is a vertex drift chamber (VTX) which provides track information up to a radius of 22 cm and for pseudo-rapidities $`|\eta |`$ 3.5. The VTX measures the $`z`$-position of the primary vertex. Both the SVX and the VTX are mounted inside the CTC, a 3.2 m long drift chamber with an outer radius of 132 cm containing 84 concentric, cylindrical layers of sense wires, which are grouped into 8 alternating axial and stereo superlayers. The solenoid is surrounded by sampling calorimeters used to measure the electromagnetic and hadronic energy of jets and electrons. The calorimeters cover the pseudo-rapidity range $`|\eta |`$ 4.2. The calorimeters are segmented in $`\eta `$-$`\varphi `$ towers pointing to the nominal interaction point. There are three separate $`\eta `$-regions of calorimeters. Each region has an electromagnetic calorimeter \[central (CEM), plug (PEM) and forward (FEM)\] and behind it a hadron calorimeter \[CHA, PHA and FHA, respectively\]. Located six radiation lengths inside the CEM calorimeter, proportional wire chambers (CES) provide shower-position measurements in the $`z`$ and $`r\varphi `$ view. Proportional chambers (CPR) located between the solenoid and the CEM detect early development of electromagnetic showers in the solenoid coil. These chambers provide $`r`$-$`\varphi `$ information only.

The calorimeter acts as a hadron absorber for the central muon detection system (CMU). The CMU consists of four layers of drift chambers located outside the CHA calorimeter. The CMU system covers the pseudo-rapidity $`|\eta |`$0.6 and can be reached by muons with $`p_T1.4`$ GeV/c. The CMU system is followed by 0.6 m of steel and four additional layers of drift chambers (CMP). The system of drift chambers CMX extends the muon detection to $`|\eta |`$ 1.0.

## III Data collection and identification of high $`p_T`$ leptons and jets

The last collider run, called Run I, lasted from August of 1992 till July of 1993 (Run 1A) and from January of 1994 till July of 1995, (Run 1B). The data collected during this run correspond to an integrated luminosity of $`105.1\pm 4.0`$ pb-1 when using the total $`p\overline{p}`$ cross section value $`80.03\pm 2.24`$ mb . We begin this section with a description of the triggers used in this analysis. This is followed by subsections on the reconstruction and identification of electrons, muons, jets and neutrinos.

### A Triggers

A three-level trigger system is used to select events originating from $`p\overline{p}`$ interactions and containing electrons, muons, jets or missing transverse energy ($`\mathrm{}\mathrm{E}_\mathrm{T}`$).

The first-level trigger (L1) accepts events based on the identification of energy clusters in the calorimeter or track segments in the muon chambers. The L1 calorimeter trigger requires a single CEM or CHA trigger tower with transverse energy greater than 8 or 12 GeV, respectively (these thresholds were set at 6 and 8 GeV during Run 1A). The L1 muon trigger infers the track momentum from the deflection of the track segment in the muon chambers due to the magnetic field; it requires a minimum transverse momentum of 6 and 10 GeV/c in the CMU and CMX chambers, respectively. A minimum energy of 300 MeV is required in the hadron calorimeter tower associated with the track segment.

The second-level trigger (L2) uses the calorimetry information with greater sophistication. The L2 trigger is a fastbus based processor  with a decision time of approximately 20 $`\mu `$s. It combines calorimetry towers forming electromagnetic and jet-like clusters. An electromagnetic cluster is constructed as a set of contiguous CEM (PEM) towers each with $`E_T`$ 7 GeV (4 GeV), including at least one seed tower with $`E_T`$ 8 GeV (6 GeV).

The L2 jet clusters are formed starting with a seed tower with $`E_T3`$ GeV and summing all contiguous towers with $`E_T1`$ GeV. A crude estimate of $`\mathrm{}\mathrm{E}_\mathrm{T}`$ is also available at this trigger level. The L2 trigger utilizes the list of $`r\varphi `$ tracks provided by the central fast tracker (CFT), a hardware processor which uses fast timing information from the CTC as input. The events used in this analysis were collected using the L2 high-$`p_T`$ electron trigger, which requires an electromagnetic cluster of transverse energy greater than 16 GeV matched by a CFT track with transverse momentum $`p_T12`$ GeV/c. A second trigger requires an electromagnetic cluster of $`E_T16`$ GeV and $`\mathrm{}\mathrm{E}_\mathrm{T}20`$ GeV and is used to recover losses due to the CFT inefficiency. The L2 high-$`p_T`$ muon trigger requires a CFT track with $`p_T12`$ GeV/c pointing within $`5^{}`$ to a L1 track segment in the muon detectors. To ensure good efficiency, additional L2 muon triggers require only a L1 track segment accompanied by at least one jet cluster with $`E_T`$ 15 GeV or $`\mathrm{}\mathrm{E}_\mathrm{T}`$ 35 GeV.

The L3 trigger decision is made after the full event reconstruction. Events accepted by the L2 trigger are processed by a farm of SGI processors running the full off-line reconstruction package. The Level 3 electron trigger requires a CEM cluster with $`E_T18`$ GeV and a reconstructed track with $`p_T13`$ GeV/c pointing to it. The ratio of hadronic to electromagnetic energy in the cluster is required to be less than 0.125. The Level 3 muon trigger requires a match within 10 cm in the $`r\varphi `$ plane between a reconstructed track with $`p_T>18`$ GeV/c extrapolated to the radius of the muon detectors and a track segment in the muon chambers.

Trigger efficiencies have been measured directly using events with overlapping triggers. The electron trigger efficiency is found to be larger than 99.6% for electrons inside the detector fiducial volume. Likewise, the muon trigger efficiency is ($`70\pm 2`$)%; this includes an inefficiency due to the fact that the muon trigger does not cover the entire detector fiducial volume. The measured trigger efficiencies have been included in the detector simulation described in Section VII. A check of the muon trigger simulation was performed by comparing the rate of $`W\mu \nu `$ events in the data to that of a simulation of this process using the Herwig generator (see Section VI) normalized to the same number of $`We\nu `$ events. We observe agreement between data and simulation within 10%, and this difference is taken as the systematic error on the muon trigger simulation.

### B Electron reconstruction

The $`\mathrm{W}+`$ jet sample is selected requiring electrons reconstructed in the central pseudo-rapidity region $`|\eta |`$ 1. Stricter cuts, described in detail in Ref. , are applied to central electron candidates which passed the trigger prerequisites. The following variables are used to discriminate against charged hadrons: (1) the ratio of hadronic to electromagnetic energy of the cluster, $`E_{had}/E_{em}`$; (2) the ratio of cluster energy to track momentum, $`E/P`$; (3) a comparison of the lateral shower profile in the calorimeter cluster with that of test-beam electrons, $`L_{shr}`$; (4) the distance between the extrapolated track-position and the CES measurement in the $`r\varphi `$ and $`z`$ views, $`\mathrm{\Delta }x`$ and $`\mathrm{\Delta }z`$; (5) a $`\chi ^2`$ comparison of the CES shower profile with those of test-beam electrons, $`\chi _{strip}^2`$; (6) the interaction vertex position, $`z_{ver}`$ and the distance between the interaction vertex and the reconstructed track in the $`z`$-direction, $`z`$-vertex match; and (7) the isolation, $`I`$, defined as the ratio of the additional transverse energy in a cone of radius $`R=0.4`$ around the cluster axis to the transverse energy of the electron cluster. The electron selection criteria are listed in Table I.

Fiducial cuts on the electromagnetic shower position, as measured in the CES, are applied to insure that the electron candidate is away from the calorimeter boundaries and the energy is well measured.

Electrons from photon conversions are removed using an algorithm based on tracking information. Electron tracks close to a companion track with opposite charge are considered conversion candidates. The following variables are used to identify and remove photon conversions: (1) the difference of the polar angles, $`\delta \mathrm{cot}\theta `$; (2) the distance between the two tracks in the $`r\varphi `$ plane at the radius $`R_{conv}`$ where the tracks are parallel, $`\mathrm{\Delta }_{sep}`$; and (3) the conversion radial position, $`R_{conv}`$. If a companion track is not found, we identify conversion candidates using $`f_{VTX}`$ which is the ratio of the measured to expected number of VTX hits associated to the electron candidate. Table II summarizes the criteria used to identify and remove electrons from photon conversions. The efficiency of the conversion algorithm is measured with a sample of photon conversions selected using the CPR detector. The efficiency of the conversion removal algorithm is 90.7 $`\pm `$ 3.8%. The fraction of electrons erroneously removed is estimated using a sample of $`Ze^+e^{}`$ events to be 2.2 $`\pm `$ 0.6% and is properly accounted for by the simulation.

The total primary electron identification efficiency has been measured using a sample of $`Ze^+e^{}`$ decays and is listed in Table III.

When an electron candidate is found, the calorimeter towers belonging to the electron cluster are not used by the jet clustering algorithm.

### C Muon reconstruction

Muons are identified in the $`|\eta |1.0`$ region by extrapolating CTC tracks to the muon detectors and matching them to track segments reconstructed in the muon chambers. The following variables, described in detail in Ref. , are used to separate muon candidates from cosmic rays and from hadrons not contained by the calorimeter: (1) an energy deposition in the electromagnetic and hadronic calorimeter characteristic of minimum ionizing particles, $`E_{em}`$ and $`E_{had}`$; (2) the distance of closest approach of the reconstructed track to the beam line, $`d`$; (3) the $`z`$-vertex match; (4) the matching distance between the extrapolated track and the track segment in the muon chamber, $`\mathrm{\Delta }x=r\mathrm{\Delta }\varphi `$; and (5) the isolation $`I`$, the ratio of additional transverse energy in a cone of radius $`R=0.4`$ around the track direction to the muon transverse momentum.

The muon selection criteria are listed in Table IV. The muon identification efficiency has been measured using a sample of $`Z\mu ^+\mu ^{}`$ decays and is listed in Table III.

Leptons passing the requirements listed in Table I and IV are labeled primary leptons. As a consequence of the high luminosity of the collider, approximately 50% of the events with a primary lepton contain multiple interactions which result in more than one primary vertex in the event. The ambiguity is resolved by selecting the vertex associated with the primary lepton track to evaluate jet pseudorapidities and the missing transverse energy.

### D jet reconstruction

The CDF jet reconstruction algorithm uses a cone of fixed radius in the $`\eta \varphi `$ space. In this analysis we use a cone of radius 0.4 which has been shown to contain approximately 70% of the jet energy . A detailed description of the jet reconstruction algorithm can be found in Ref. .

The jet energy resolution can be parametrized as $`\sigma (E_T)/E_T1/\sqrt{E_T}`$, where $`E_T`$ is measured in GeV. Effects which contribute to the resolution are the lower calorimeter response at the boundaries of different towers and of different calorimeter detectors, the loss of low momentum particles inside the magnetic field, the energy deposition in towers outside the clustering cone, the contribution of the underlying-event and energy losses due to minimum ionizing particles or neutrinos present in the jet. Corrections meant to reproduce the average jet $`E_T`$ correctly (without improving the energy resolution) are often used . The average jet energy correction factor ranges from approximately 1.7 to 1.1 as the jet transverse energy increases from 15 to 100 GeV.

Checks of the jet energy corrections have been performed in Ref.  by studying the momentum balance in $`\gamma +`$ jet and $`Z+`$jet events. The energy imbalance is measured to be within 3% of the $`Z`$ or photon energy. However, the uncertainty in the modeling of the large-angle gluon emission results in a 10% systematic uncertainty of the jet energy scale.

### E $`\mathrm{}\mathrm{E}_\mathrm{T}`$ and neutrino reconstruction

The presence of neutrinos is inferred from transverse energy imbalance in the detector. The transverse missing energy is defined as

$$\mathrm{}\mathrm{E}_\mathrm{T}=|\underset{i}{}E_T^i\stackrel{}{n}_i|$$

where, $`E_T^i`$ is the magnitude of the transverse energy contained in each calorimeter tower $`i`$ in the pseudo-rapidity region $`|\eta |<3.5`$ and $`\stackrel{}{n}_i`$ is the direction of the tower in the plane transverse to the beam direction. When a muon is present in the event, $`\mathrm{}\mathrm{E}_\mathrm{T}`$ is calculated as:

$$\mathrm{}\mathrm{E}_\mathrm{T}=|\underset{i}{}E_T^i\stackrel{}{n}_i+\stackrel{}{p}_T^\mu \left(1\frac{E_T^\mu }{p_T^\mu }\right)|$$

where $`E_T^\mu `$ is the transverse energy deposited by the muon in the calorimeter and $`p_T^\mu `$ is the muon transverse momentum.

## IV Selection of the $`W`$ and $`Z+`$ jet samples

The $`\mathrm{W}+`$ jet sample, which contains the $`t\overline{t}`$ signal, is selected from the high-$`p_T`$ inclusive lepton data set by requiring at least one primary electron with $`E_T`$ 20 GeV or one primary muon with $`p_T`$ 20 $`\mathrm{GeV}/\mathrm{c}`$, $`\mathrm{}\mathrm{E}_\mathrm{T}`$ 20 GeV and at least one jet with uncorrected transverse energy $`E_T`$ 15 GeV and pseudo-rapidity $`|\eta |`$ 2. An appreciable fraction of these events is due to $`Z+`$ jet production. Some $`Z`$ events can be identifyied and removed when the second lepton from the $`Z`$ decay falls into the detector acceptance. Because $`W+`$ jet and $`Z+`$ jet events have similar production mechanisms, we will use the $`Z+`$ jet sample to check our evaluation of the backgrounds to $`t\overline{t}`$ production. It is also interesting to study this sample because events in which one of the two leptons is not identified (unidentified $`Z`$’s) are a background to $`t\overline{t}`$ production. The following subsection explains the removal of dilepton events. The events surviving dilepton removal constitute the $`W+`$ jet sample which is described in the last subsection.

### A Selection of the $`Z+`$ jet sample

$`Z`$ candidates are selected from the high-$`p_T`$ lepton data set by requiring a primary electron with $`E_T`$ 20 GeV and by searching for a second lepton of opposite charge which satisfies the criteria listed in Table V.

Searching for additional electrons we relax the isolation and $`E_{had}/E_{em}`$ cuts. We also search in the PEM and FEM detectors. Additional muons are searched for by relaxing all selection cuts defining primary muons. As shown in Table V, CTC tracks without a match to a track segment in the muon chambers but pointing to a calorimeter tower with a small energy deposition are also considered muon candidates.

Events are flagged as $`Z`$ candidates if the invariant mass of the lepton pair falls in the range $`70M_{\mathrm{}\mathrm{}}110\mathrm{GeV}/\mathrm{c}^2`$ (see Figure 1). The number of $`Z`$ candidate events as a function of the jet multiplicity is shown in Table VI.

### B Dilepton removal

All events containing a primary lepton and at least one additional lepton selected using the criteria listed in Table V are removed from the $`\mathrm{W}+`$ jet sample. These events arise from $`Z\tau ^{}\tau ^+`$, di-bosons, Drell-Yan and $`t\overline{t}`$ production. The $`\sigma _{t\overline{t}}`$ production cross section using dilepton events has been measured in Ref.  and we want to avoid obvious correlations.

We also remove events containing an isolated track with $`p_T10`$ GeV/c with charge opposite to the primary lepton . The majority of these events originates from genuine dilepton events in which one lepton is outside the region covered by the calorimeters or the muon detectors.

Finally, to remove dileptons missed due to inefficiencies of the tracking system, we remove events in which a jet with $`E_T15`$ GeV and $`|\eta |2`$ has a large electromagnetic fraction ($`E_{em}/E_{em+had}0.95`$) and less than three tracks. These types of events are mostly produced by $`Ze^+e^{}`$ decays.

The dilepton removal reduces the acceptance for $`t\overline{t}`$ events by 17.2%.

### C The $`W+`$ jet sample

The number of $`W`$ events surviving the $`Z`$ and dilepton removal is listed in Table VII as a function of the jet multiplicity. The transverse mass distribution of the $`\mathrm{W}`$ candidates is shown in Figure 2.

## V Additional data samples

In addition to the $`Z+`$ jet sample, we use a number of independent data sets for the purpose of calibrating the Monte Carlo generators and the detector simulation. The generic-jet samples are described in Subsection A. We will use these sample to derive the new parametrization of the mistag rate, to check our evaluation of the efficiency of the tagging algorithms, and to calibrate the calculation of the fraction of $`W+`$ jet events with heavy flavor. Subsection B describes the low-$`p_T`$ inclusive lepton sample which will be used to determine the efficiency of the tagging algorithms. Finally, Subsection C details the selection of the isolated photon sample. We will use this sample to check the parametrization of the mistag rate of our tagging algorithms.

### A Generic-jet samples

The samples JET 20, JET 50, JET 70 and JET 100 are data collected requiring the presence of a L2 calorimeter cluster with transverse energy $`E_T`$ 20, 50, 70 and 100 GeV, respectively.

The samples $`E_T`$ 175 and $`E_T`$ 300 are data collected requiring the scalar sum of the transverse energy of all calorimeter towers, as evaluated by the L2 trigger, to be larger than 175 and 300 GeV, respectively.

The last generic-jet sample, $`E_T`$ 125 4CL, was collected requiring the presence of four L2 calorimeter clusters with $`E_T15`$ GeV and the scalar sum of the transverse energy of all calorimeter towers to be larger than 125 GeV.

The L2 triggers calculate the above quantities with respect to the nominal interaction point. Off-line we take as the event vertex the one with the largest $`{\displaystyle \underset{i}{}}p_T^i`$ using all tracks $`i`$ associated with the vertex. We retain the events in which the L2 requirements are also matched after the event is reconstructed using this vertex. In these events, we inspect all jets with $`E_T`$ 15 GeV and which contain at least two SVX tracks (taggable jets).

### B The low-$`p_T`$ inclusive lepton sample

The efficiency of the $`b`$-tagging algorithms needs to be measured in a sample enriched in $`b\overline{b}`$ production. The low-$`p_T`$ electron sample is collected with the L2 requirement that a CFT track with $`p_T`$ 7.5 GeV/c is matched by an electromagnetic L2 cluster with $`E_T`$ 8 GeV. The fraction of electrons coming from semileptonic $`b`$-decays is enhanced with the selection criteria listed in Table VIII. We use electrons in the CEM fiducial region and remove photon conversion candidates. We require the lepton to be in a cone of radius 0.4 around the direction of a taggable jet. We require also the presence of at least one additional taggable jet. The $`b`$-purity of this sample is approximately 50%.

We check the results obtained using the low-$`p_T`$ electron sample using a lower statistics low-$`p_T`$ muon sample collected using the inclusive muon trigger. In this case, a CFT track with $`p_T`$ 7.5 GeV/c must be matched to a reconstructed track-segment in both sets of central muon detector (CMU+CMP). Central muons which passed the trigger prerequisite are selected with the same criteria used for the high-$`p_T`$ muons listed in Table V (we require $`I`$ 0.1).

### C The isolated photon sample

The isolated photon sample was collected requiring a L2 isolated electromagnetic cluster with $`E_T`$ 16 GeV and with less than 5 GeV of additional energy in a 5$`\times `$10 grid of calorimeter towers centered on the photon direction . Photon candidates which pass the L3 trigger must be in the good fiducial region of the calorimeter and there must be less than 4 GeV in a cone of radius 0.7 around the photon direction. Table IX summarizes the off-line criteria used for the selection of the photon sample. After requiring the presence of an additional jet with $`E_T`$ 15 GeV and $`|\eta |2`$, the final sample consists of $`3000\gamma +1`$ jet events. The expected background contamination of the sample due to $`\pi ^0`$ and $`\eta `$ decays is estimated to be ($`45.0\pm 4.5`$)% .

## VI Monte Carlo generators

In this analysis we use three Monte Carlo generators, Herwig , Pythia  and Vecbos . The acceptance for $`t\overline{t}`$ events is calculated using Pythia version 5.7. The $`t\overline{t}`$ acceptance has been also evaluated using the version 5.6 of Pythia and Herwig 5.6. Herwig is used to estimate the fraction of $`W+`$ 1 jet events with heavy flavor (using the process 2100). We first calibrate the Herwig simulation using generic-jet data, as described in Section XI.

Both Herwig and Pythia generators use tree-level matrix element calculations for the parton hard scattering, convoluted with parametrizations of the parton distribution functions. The outgoing initial and final state partons are converted into a cascade of gluons and $`q\overline{q}`$ pairs with energy and angular distributions determined by the Altarelli-Parisi equations . The strength of these generators is the modeling of the parton shower which accounts for the color correlation between the initial and final state partons. The parton shower terminates when the invariant mass of the parton falls below the perturbative QCD scale. At this level the partons are turned into colorless hadrons according to phenomenological models (the process is called hadronization or fragmentation). For $`b`$ and $`c`$-quarks the fragmentation is modeled in Pythia with the Peterson parametrization . We use the fragmentation parameter $`ϵ=0.006`$ for $`b`$-quarks and $`ϵ=0.05`$ for $`c`$-quarks. Herwig uses its own hadronization model, the settings for which are listed in Ref. . Both generators include a model of the underlying event which describes the hadronization products of the beam remnants.

The Vecbos Monte Carlo is used to study the part of the phase-space in the $`W+`$ 1 jet production that is not treated correctly by parton shower Monte Carlos, specifically $`Wb\overline{b}`$ and $`Wc\overline{c}`$ events in which the two $`b`$ or $`c`$-partons produce two well separated jets. The Vecbos Monte Carlo generator provides a parton level calculation of the $`W+`$ n jet cross section based on the leading order matrix elements of the hard scattering. Infrared and collinear singularities are regulated by requiring that the final-state partons have a transverse momentum exceeding a cutoff value $`p_T^{\mathrm{min}}`$ and are separated by more than $`R_{\mathrm{min}}`$ ($`R=\sqrt{(\mathrm{\Delta }\varphi )^2+(\mathrm{\Delta }\eta )^2}`$). We use $`p_T^{\mathrm{min}}`$= 8 GeV/c and $`R_{\mathrm{min}}`$ =0.4. We use the renormalization scale $`Q^2`$=$`<p_T>^2`$, where $`<p_T>`$ is the average transverse momentum of the outgoing partons. We have verified that after our selection cuts the fraction of jets with heavy flavor calculated with Herwig matches the Vecbos prediction at the $`R_{\mathrm{min}}`$ threshold. We transform the partons produced by Vecbos into hadrons and jets using the Herwig program adapted to perform the coherent shower evolution of both initial and final-state partons .

In summary, we use Herwig to predict the fraction of $`W+`$ 1 jet events where only one jet clustered in a cone of radius 0.4 contain $`b`$ or $`c`$-hadrons while we rely on Vecbos to extend the prediction to the cases where two different jets contain both heavy-flavored hadrons.

We use the MRS D$`{}_{}{}^{}{}_{0}{}^{}`$ set of parton distribution functions  to generate $`W+`$ jet events because it has been shown to reproduce the results of the $`W`$ asymmetry measured by CDF .

The decay of hadrons with heavy flavor produced by the Monte Carlo generators is modeled using the CLEO Monte Carlo generator (QQ) . We use the QQ table of branching ratios for each decay but our own lifetime table because decay lengths are modeled inside the detector simulation.

## VII Detector Simulation

The QFL detector simulation is used to decay all generated particles and model their interactions with the various elements of the CDF detector. The detector response is based upon parametrizations and simple models which depend on the particle kinematics. The calorimeter simulation is based upon a parametrization of the calorimeter response to single particles parametrized as a function of the pseudo-rapidity and azimuthal angle (to account for cracks in the calorimetry ) and of the transverse momentum using test-beam data. After the simulation of the CDF detector, the Monte Carlo events are treated as if they were real data.

### A CTC track simulation

The CTC simulation is not a hit-level simulation. It converts each particle’s momentum vector at generator level into a reconstructed track using covariance matrices derived from the data. Not surprisingly, the track-reconstruction efficiency in the detector simulations is higher than that measured in the data. The major factor influencing the track reconstruction efficiency is the density of hits in the tracking detector. In this respect, the problem is aggravated by the fact that the Monte Carlo generators do not contain multiple interactions.

To adjust the tracking reconstruction efficiency in the simulation, CTC hits of Monte Carlo generated tracks have been embedded in generic-jet data. The efficiency is determined by counting the fraction of times the embedded track is reconstructed. The tracking efficiency is measured as a function of the hit density around the track for low luminosity runs, and then for runs of typical luminosities. Table X compares the track reconstruction efficiency in the detector simulation to the efficiency for reconstructing simulated tracks embedded in the data. The degradation of the track reconstruction efficiency is parametrized in the detector simulation as a function of the number of hits around the tracks and of the average luminosity of the data. This procedure accounts well for the dependence of the tracking efficiency on the jet transverse energy.

### B Lepton identification efficiencies

Aside from the efficiency for reconstructing a track, the primary lepton identification efficiency in the simulation depends also on how well the Monte Carlo generators reproduce the isolation distribution and how well the calorimeter response has been parametrized. In the simulation, the primary lepton identification efficiencies are measured as the ratio of the number of leptons passing the selection cuts listed in Tables I and IV to the number of leptons generated in the kinematical acceptance. The identification efficiencies in the simulation are 97 $`\pm `$ 2% for muons and 87.5 $`\pm `$ 2.0% for electrons. The identification efficiencies for primary leptons are degraded in the detector simulation to match the ones measured in the data (see Table III). Altogether, we degrade the rates of simulated primary leptons by the factor of 0.936$`\pm `$0.125 (the error includes a 10% uncertainty on the muon trigger simulation).

The efficiency for identifying soft lepton tags is a far more complicated problem because some detector responses, such as $`dE/dx`$ in the CTC and the CPR chambers, have not been parametrized in the detector simulation. The SLT simulation weights tracks corresponding to leptons from $`b`$ and $`c`$-quark decays at generator level with a parametrization of the efficiency of each selection cut measured using the data, as described in Section VIII C.

### C SVX track simulation

The detector simulation becomes unwieldy when simulating tracks that are measured by both the CTC and SVX tracking detectors as is the case for input tracks to the SECVTX and jet-probability algorithms. The SVX track reconstruction is performed by assigning hits on the silicon vertex detector to previously reconstructed CTC tracks. In the data hits are assigned if they are contained in a road around the reconstructed CTC track determined by its uncertainty (4 $`\sigma `$ in the $`r\varphi `$ plane). A CTC track with at least two associated SVX hits is defined to be a SVX track and is refitted using the SVX hits and the CTC track parameters and covariance matrix. The simulation of the SVX is a hit-level simulation in which the hit resolution is taken from the data. Simulated SVX tracks are reconstructed as in the data. However, in the data we must multiply all the elements of the covariance matrix by a factor of 2 so that the CTC-SVX matching uncertainty agrees with the measured resolution  while there is no such need in the simulation.

The efficiency for finding SVX tracks in the detector simulation also needs to be degraded, by a factor determined by measuring the efficiency for reconstructing Monte Carlo generated tracks embedded at hit-level in generic-jet data (see Table X).

Having done this, the simulation is still not a perfect reflection of the data. For example, as shown in Section VIII B, the distribution of the impact parameter significance of SVX tracks in the data and in the detector simulation are slightly different. We conclude that it is neccessary to measure the tagging efficiencies of each algorithm in the data and in the simulation and correct the detector simulation for any observed difference. This is done in Section IX.

## VIII Description of the tagging algorithms

The presence of jets originating from $`b`$ quarks is one of the characteristic signatures of $`t\overline{t}`$ events. Following previous work , we tag $`b`$-quarks using two of their distinctive properties: the relatively long lifetime and the presence of semileptonic decays. Two tagging techniques based on tracking information using the SVX detector have been developed to identify jets containing heavy flavor. The Secondary Vertex Tagging algorithm (SECVTX) is described in Subsection A. The jet-probability algorithm, used to check SECVTX results, is described in Subsection B. The Soft Lepton Tagging algorithm (SLT) is discussed in Subsection C, which also includes the evaluation of the SLT fake rate and a description of the simulation of this algorithm.

### A SECVTX algorithm

The SECVTX algorithm is described in more detail in Refs. . SECVTX is based on the determination of the primary event vertex and the reconstruction of additional secondary vertices using displaced tracks associated with jets.

The positions of the $`p\overline{p}`$ interactions (primary vertices) are distributed along the beam direction according to a Gaussian with a width of approximately 28 cm. In the plane transverse to the beam axis, these interactions follow a distribution that is a Gaussian with a width of 25 $`\mu `$m in both the $`x`$ and $`y`$ dimensions. To reconstruct the primary event vertex, we first identify its $`z`$-position using the tracks reconstructed in the VTX detector. When projected back to the beam axis, these tracks determine the longitudinal position with a precision of about 0.2 cm.

The transverse position of the primary vertex is determined for each event by a weighted fit of all SVX tracks which have a $`z`$ coordinate within 5 cm of the $`z`$-vertex position of the primary vertex associated with the trigger lepton. First, all tracks are constrained to originate from a common vertex. The position of this vertex is constrained by the transverse beam envelope described above. Tracks that have impact parameter significance $`|d|/\sigma _d`$, where $`\sigma _d`$ is the estimate of the uncertainty on the impact parameter $`d`$, larger than three with respect to this vertex are removed and the fit is repeated. This procedure is iterated until all used tracks satisfy the impact parameter requirement. At least five tracks must be used in the determination of the transverse position of the primary vertex or we use the nominal beam-line position. The primary vertex coordinates transverse to the beam direction have uncertainties in the range of 10$``$25 $`\mu `$m, depending on the number of tracks and the event topology.

The search for a secondary vertex in a jet is a two stage process. In both stages, tracks in the jet are selected based on the significance of their impact parameter with respect to the primary vertex. The first stage (see Table XI) requires at least three candidate tracks for the reconstruction of the secondary vertex. Tracks consistent with coming from the decay $`K_s\pi ^+\pi ^{}`$ or $`\mathrm{\Lambda }\pi ^{}p`$ are not used as candidate tracks. Two candidate tracks are constrained to pass through the same space point to form a seed vertex. If at least one additional candidate track is consistent with intersecting this seed vertex, then the seed vertex is used as the secondary vertex. If the first stage is not successful in finding a secondary vertex, a second pass is attempted. More stringent track requirements (on $`|d|/\sigma _d`$ and $`p_T`$, for example) are imposed on the candidate tracks. All candidate tracks satisfying these stricter criteria are constrained to pass through the same space point to form a seed vertex. This vertex has an associated $`\chi ^2`$. Candidate tracks that contribute too much to the $`\chi ^2`$ are removed and a new seed vertex is formed. This procedure is iterated until a seed vertex remains that has at least two associated tracks and an acceptable value of $`\chi ^2`$. Table XI lists the selection criteria used for the determination of the secondary vertex candidates.

The decay length of the secondary vertex $`L_{xy}`$ is the projection of the two-dimensional vector pointing from the primary vertex to the secondary vertex on the jet axis; if the cosine of the angle between these two vectors is positive (negative), then $`L_{xy}`$ is positive (negative). Most of secondary vertices from the decay of $`b`$ and $`c`$-hadrons are expected to have positive $`L_{xy}`$. Secondary vertices from random combination of mismeasured tracks are expected to have a symmetric distribution around $`L_{xy}`$=0 . To reduce the background from false secondary vertices (mistags), a jet is considered tagged by SECVTX if it contains a secondary vertex with $`\frac{L_{xy}}{\sigma _{L_{xy}}}3.0`$, where $`\sigma _{L_{xy}}`$ is the estimated uncertainty on $`L_{xy}`$ ($``$ 130 $`\mu `$m). The mistag contribution to positive SECVTX tags is evaluated starting from the rate of negative SECVTX tags and detailed in Section X.

### B Jet-probability algorithm

The jet-probability $`b`$-tagging algorithm  is used to cross-check the SECVTX results. The jet-probability algorithm compares track impact parameters to measured resolution functions in order to calculate for each jet a probability that there are no long lived particles in the jet cone. This probability is uniformly distributed between 0 and 1 for light quark or gluon jets, but is very small for jets containing displaced vertices from heavy flavor decays. We briefly describe the transformation from the track impact parameters to the jet-probability measure.

The track impact parameter significance $`S`$ is defined as the value of the impact parameter $`d`$ divided by its uncertainty $`\sigma _d`$. Tracks used in the calculation of jet-probability are required to satisfy the quality criteria listed in Table XII. The sign of the impact parameter significance is defined to be positive if the point of closest approach to the primary vertex lies in the same hemisphere as the jet direction, and negative otherwise. Figure 3 shows the distribution of the impact parameter significance of tracks in the JET 50 sample. This distribution is fitted with the resolution function $`(S)`$.

The negative side of the resolution function $`(S)`$ derived using JET 50 data is used to determine the probability $`P(S_0)`$ that the impact parameter significance $`S_0`$ of a given track is due to the detector resolution:

$$P(S_0)=\frac{_{\mathrm{}}^{|S_0|}(S)𝑑S}{_{\mathrm{}}^0(S)𝑑S}$$

Figure 4 shows that the impact parameter significance distribution of tracks in the JET 50 data and in the corresponding simulation are slightly different. The resolution functions $`(S)`$ are therefore defined separately for the data and the simulation in order to account for the differences in the resolution between the true and the simulated detector performance.

The probability that a jet is consistent with a zero lifetime hypothesis is calculated as:

$$\text{JPB}=\underset{k=0}{\overset{N1}{}}\frac{(\mathrm{ln})^k}{k!}$$

where $``$ is the product of the individual probabilities of the $`N`$ SVX tracks in a jet which satisfy the criteria listed in Table XII. Jet-probability is defined using tracks with positive impact parameter and requiring $`N`$ 2. We also define a negative jet-probability in which we select only tracks with negative impact parameter in the calculation. This is used as a control sample and a check of our method.

Figure 5 shows the positive and negative jet-probability distributions in a sample of JET 50 and JET 140 data. The positive jet-probability distribution shows jets containing hadrons with heavy flavor as a large excess at jet-probabilities smaller than 0.05 over a flat distribution. A jet has a positive JPB tag if the jet-probability value is smaller than 0.05.

The negative jet-probability distribution is quite flat, as expected, since the resolution files were constructed using tracks with negative impact parameter. The small excess at negative jet-probability smaller then 0.05 (negative JPB tags) is due to the increase of the fraction of jets from heavy flavor in the JET 140 with respect to the JET 50 data. This excess largely disappears, as shown in Figure 6, when plotting the negative probability of jets which have a large positive jet-probability (0.1$``$1). Since tracks with negative signed impact parameter in JET 50 data are used to define the resolution function, the small contribution to negative tags from jets with heavy flavor is uncorrectly attributed to the detector resolution by this procedure. It will be accounted for in the evaluation of the JPB mistags in Section X.

Ideally JPB tags corresponding to jet-probability values smaller than 0.05 should contain a 5% mistag rate. This expectation is tested in Figure 6 fitting a first order polynomial function to the jet-probability distribution in the interval 0.1$``$1.0. The extrapolation of the fitted function predicts 4441$`\pm `$34 negative JPB tags while 4455 are observed; this corresponds to 4.94% of the total number (101050) of jets in the sample.

### C SLT algorithm

The SLT algorithm tags $`b`$ quarks by searching for an electron or muon from their decay. Low momentum or soft leptons can also result from $`b`$-hadron decay through sequential $`c`$-decays, or $`\tau `$ and $`J/\psi `$ cascade decays. This analysis follows the guidelines for the identification of soft electrons or soft muons documented in Refs. . While previous measurements of the $`t\overline{t}`$ cross section used rates of events with SLT tags, in this analysis we search for soft lepton candidates only in a cone of radius 0.4 around the axis of a jet with $`E_T15`$ GeV and $`|\eta |2`$.

To search for soft electrons, every CTC track with $`p_T2`$ GeV/c, which is associated to a jet, is extrapolated into the fiducial region of the calorimeter and is matched to a CES cluster. The matched CES cluster is required to be consistent in shape and position to the expectations for an electron. In addition, we require $`0.7E/P1.5`$ and $`E_{had}/E_{em}0.1`$. The energy deposited by the track in the preradiator (CPR) is required to be consistent with an electron shower. The track ionization rate ($`dE/dx`$), derived from the charge deposition of the CTC hits associated with the track, is also required to be consistent with the electron hypothesis. Electrons from photon conversions are removed. Photon conversions are identified as combinations of the electron candidate and an additional track with opposite charge passing the criteria listed in Table IV with the additional requirement that the invariant mass be smaller than 500 MeV/c2. The selection criteria used to define the soft electron are described in more detail in Ref. . The efficiency of each criteria used to select soft electron candidates has been measured using a sample of electrons produced by photon conversions  (the efficiency of the $`E/P`$ and $`E_{had}/E_{em}`$ cuts is calculated using the simulation).

Soft muons are identified by matching CTC tracks with $`p_T`$ 2 $`\mathrm{GeV}/\mathrm{c}`$ to track segments in the CMU, CMP and CMX muon chambers. Muons candidate tracks with $`p_T`$ 3 $`\mathrm{GeV}/\mathrm{c}`$ that are extrapolated to the fiducial volume of both the CMU and CMP system are required to be matched to track segments in both muon detectors. To maintain high efficiency for non-isolated muons, we do not impose minimum-ionization requirements on the calorimeter deposition. However, in order to reduce hadronic punch-through in the region not covered by the CMP system, we check that the energy, $`E_{had}`$, in the tower traversed by muon candidates with $`p_T`$ 6 $`\mathrm{GeV}/\mathrm{c}`$ is consistent with the muon hypothesis; we require $`E_{had}6+p`$, where $`p`$ is the scalar sum of the momenta of all tracks contained in a cone of radius 0.2 around the muon direction. The efficiency of each selection cut has been measured using a sample of $`J/\psi \mu ^+\mu ^{}`$ and $`Z\mu ^+\mu ^{}`$ decays .

Figures 7 and 8 show distributions of the invariant mass between primary and soft leptons in all $`W+`$1 jet events. As shown in Figure 8, there is a handful of events where the soft muon is consistent with being the second leg of a $`Z`$ boson decay embedded in a jet. Soft muons which, when combined with the primary muons of opposite charge, yield an invariant mass $`70M_{\mu \mu }110`$ GeV/c2 are not considered tags.

#### 1  Fake soft lepton tags

This background includes hadrons which pass the lepton selection cuts (such as pions which fake an electron or a muon) as well as electrons from conversions or muons from pions or kaons which decay in the detector. Because of the complication of accurately modeling these processes in the QFL detector simulation, this background is estimated using the data.

The SLT fake rate is measured starting from the ratio of the number of tracks passing the soft lepton selection criteria to the total number of tracks which satisfy the soft lepton fiducial requirements in generic-jet data . In the JET 20, JET 50 and JET 70 samples the probability $`P`$ that a track produces a SLT tag is computed separately for electrons and for different types of muon detectors (CMU,CMP and CMX). This probability is parametrized as a function of the track $`p_T`$ and isolation . Since in this analysis we search a jet for SLT candidates in a cone of radius of 0.4 around its axis, we define a SLT probability per jet $`P_{SLT}^{jet}(N)=_{i=1}^NP_{SLT}^{jet}(i1)+(1P_{SLT}^{jet}(i1))\times P_i`$ where $`N`$ is the number of tracks contained in a cone of radius 0.4 around the jet axis.

In Table XIII the observed rates of SLT tags in various generic-jet samples are compared to the rates predicted by the probability $`P_{SLT}^{jet}`$ described above. Since in generic-jet data the trigger jet is biased toward a lower yield of soft muons (a jet containing a muon has a lower energy deposition in the calorimeter and therefore is less likely to be the trigger jet) the comparison is performed with and without the trigger jet. However, when more than one jet is above the trigger threshold, all jets are considered. Excluding trigger jets from the comparison one observes agreement between the observed and predicted rates of tagged jets. The last seven samples shown in Table XIII were not used to determine the SLT probability per track. Predicted and observed yields of SLT tags in all samples agree within 15%. As the amount and type of heavy flavor changes appreciably in different QCD samples (see Section X) the apparent agreement suggests that the rate of SLT tags in generic-jet data is dominated by fakes.

The SLT fake probability is obtained by removing the contribution of SLT tags due to heavy flavor decays in the generic-jet data used to construct the SLT probability per track. For this purpose, we use the signed impact parameter significance distribution of the soft lepton tracks. The distribution observed in the data is fitted with the shape expected for leptons coming from the decay of $`b`$ and $`c`$-hadrons, derived in simulated events, in addition to the shape of fake SLT tags. The shape for fake SLT tags is derived using all tracks taggable by the SLT algorithm in events which do not contain any SECVTX, JPB or SLT tags.

Figure 9 shows the signed impact significance distribution of SLT tags in JET 50 data along with the fit result. The composition of the SLT tags determined from these fits is ($`74.0\pm 3.2)\%`$ fakes, ($`10.5\pm 2.3)\%`$ $`b`$’s and ($`14.5\pm 4.3)\%`$ $`c`$’s for all three generic-jet samples used to evaluate the SLT tagging probability. Due to the systematic associated with using simulated shapes for the impact parameter distribution of lepton coming from $`b`$ and $`c`$-decays, we estimate a 10% error on the fraction of fake SLT tags determined by the fits. Based on this result, the SLT mistag probability per jet is obtained by rescaling the SLT tagging probability in generic-jet data by ($`74.0\pm 7.4`$)%.

#### 2 Simulation of the SLT algorithm

A main feature of the soft lepton tagger is that it was developed from studying real leptons from photon conversions and $`J/\psi `$ mesons. The efficiency of each selection criteria is well measured in the data. Therefore, the simulation of the soft lepton tagger does not need to rely on the QFL modeling of the detector response in order to estimate the tagging efficiency. The SLT simulation matches tracks produced by QFL to electrons and muons at generator level. The electrons or muons are required to come from $`b`$ or $`c`$ decay or any of their cascade decays. Electron tracks are extrapolated to the CPR and CES detectors, and required to pass fiducial cuts. Electron candidates are eliminated if they are consistent with arising from photon conversions. Muon tracks are extrapolated to the muon detectors, required to pass the fiducial cuts and classified according to the muon detector type. Finally tracks are weighted with the measured efficiencies of the selection criteria, which are functions of the track transverse momentum . This procedure ensures that the simulation accurately models the soft lepton tagging efficiency.

In Section XI we compare rates of SLT tags in generic-jet data to the corresponding simulation to verify that the procedure has been implemented correctly. By construction, the SLT simulation does not produce mistags. However, when applying the SLT mistag probability to tracks in simulated generic-jet data, we predict the same rate of fake SLT tags per jet that is observed in the data. This implies that tracks in simulated jets have the same multiplicity and transverse momentum distributions as in the data.

## IX Efficiency of the SECVTX and JPB taggers

We first describe the calibration of the efficiency of the tagging algorithms in the simulation. For this purpose, we use the low-$`p_T`$ inclusive electron sample described in Section V B and the corresponding simulation. A large fraction of the events in this sample is expected to originate from $`b\overline{b}`$ production in which a jet containing an electron from a semileptonic $`b`$-decay, called e-jet, recoils against a jet from the other $`b`$, called the away-jet or a-jet. The tagging efficiency in the simulation, $`\epsilon _b^{\mathrm{MC}}`$, is adjusted to the value $`\epsilon _b`$ of the tagging efficiency in the data using the scale factor

$`SF={\displaystyle \frac{\epsilon _b}{\epsilon _b^{\mathrm{MC}}}}.`$

Following the derivation of the scale factor, Subsections A-I discuss the various sources of systematic uncertainty and also present cross-checks. In Subsection J we provide an explanation for the deviation of the scale factor from unity.

The data sample consists of 55248 events. The simulated sample is generated with Herwig  (process 1500, $`22`$ hard scattering with transverse momentum threshold $`p_T^{\mathrm{min}}13`$ GeV/c). Using the generic hard parton scattering, $`b\overline{b}`$ and $`c\overline{c}`$ pairs are produced through processes of order $`\alpha _s^2`$ as $`ggb\overline{b}`$ (direct production). Processes of order $`\alpha _s^3`$ are implemented in the generator through flavor excitation processes such as $`gbgb`$ or gluon splitting, where the process $`gggg`$ is followed by $`gb\overline{b}`$. We use the MRS(G) set of parton distribution functions . Apart from the parton distribution functions, the simulation package is the same as that used to generate $`W+`$ jet events. The generated hard scattering sample corresponds to an integrated luminosity of 83.5 pb-1. In this sample we select events with an e-jet containing hadrons with heavy flavor. After applying the same selection used for the data, the simulated sample low-$`p_T`$ electron sample contains 16547 events.

Table XIV shows the composition of the simulated inclusive electron sample. One notices that 80% of the e-jet are due to $`b\overline{b}`$ production and that only 33% of the away-jets contain heavy flavor.

In the simulated events where the away jet is tagged by SECVTX ($`ST_{ajet}`$), 94% of the electron-jets are due to $`b\overline{b}`$ production. It is therefore convenient to measure the $`b`$-tagging efficiency as the fraction of these events in which the electron-jet is tagged by SECVTX or JPB

$`\epsilon _b^{MC}={\displaystyle \frac{SDT}{ST_{ajet}}}`$ (1)
where $`SDT`$ is the number of events where both the electron and away-jet contain heavy flavor and are tagged. The $`b\overline{b}`$ production accounts for 99% of the simulated events with a double tag.

Table XV lists rates of tags in the data and in the simulation. In the simulation there are very few mistags and they are easily identified because the jet does not contain $`b`$ or $`c`$-hadrons in a cone of radius 0.4 around its axis. In the data, the rate of mistags is evaluated using the parametrization described in Section X.

We use the simulation to describe $`F_{hf}`$, the fraction of data in which electron-jets contain hadrons with heavy flavor. The data contain also a relevant number of e-jets in which the electron is not associated with the production of hadrons with heavy flavor (mostly from photon conversions in jets due to light quarks or gluons). In these events, the electron-jet contributes only mistags. To describe the remaining fraction ($`1F_{hf}`$) of the data, we make the additional assumption that away-jets in these events contain the same fraction of heavy flavor as generic-jets. The parametrization of the probability of tagging jets with heavy flavor in generic-jet data is derived in Section X. The 10% uncertainty associated with this parametrization is discussed in Section IX A.

We use the following procedure to derive the tagging efficiency scale factor separately for SECVTX and jet-probability, together with the heavy flavor purity $`F_{hf}`$ of the data. The data and the simulation are normalized to the same number of tagged electron-jets that contain heavy flavor, $`T_{ejet}`$ and $`ST_{ejet}`$, through the coefficient

$`\alpha ={\displaystyle \frac{T_{ejet}}{ST_{ejet}}}.`$

Before tagging, the heavy flavor purity of the data is therefore given by

$`F_{hf}={\displaystyle \frac{\alpha SN_{ejet}}{SFN_{ejet}}}`$ (2)
where $`N_{ejet}`$ and $`SN_{ejet}`$ are the number of e-jets in the data and the simulation and $`SF`$ is the tagging efficiency scale factor. We start assuming $`SF=1`$. It follows that in the data the number of events in which a tagged away-jet with heavy flavor is associated to an electron-jet without heavy flavor is

$`T_{ajet}^{QCD}=(1F_{hf})N_{ejet}P_{QCD}`$

and that the number of events in which a tagged away-jet containing heavy flavor is associated with an electron-jet also containing heavy flavor is

$`T_{ajet}^{SEC}T_{ajet}^{QCD}.`$

For the data the $`b`$-tagging efficiency, analogy of equation (1), is then

$`\epsilon _b={\displaystyle \frac{DT}{T_{ajet}^{SEC}T_{ajet}^{QCD}}}.`$

where, as before, $`DT`$ is the number of events in which the a-jet is tagged by SECVTX and the e-jet has a SECVTX or JPB tag.

The ratio of the tagging efficiencies in the data and simulation yields the scale factor

$`SF={\displaystyle \frac{\epsilon _b}{\epsilon _b^{MC}}}.`$

The value of the scale factor is inserted again in the equation (2) and we iterate until the scale factor value is stable to within 1%.

Using the numbers of electron and away-jets listed in Table XV, we derive $`SF`$=1.23$`\pm `$0.07 for SECVTX and 0.96$`\pm `$0.05 for jet-probability. The error accounts for the sample statistics (with the largest contribution coming from the simulation) and for 10% uncertainties in the evaluation of the mistag rates and in the prediction of the rate of tags in generic-jets with heavy flavor.

The $`b`$-purity of the e-jets before tagging, $`F_{hf}`$ = (43.5 $`\pm `$ 2.9)%, is in agreement with the measurement in Ref. , 37$`\pm `$8%, using the fraction of tagged electron-jets that also contain a muon of opposite charge.

The average SECVTX tagging efficiency is 36.7$`\pm `$1.9% in the data and (29.8$`\pm `$1.1)% in the simulation. The corresponding numbers for jet-probability are (39.2$`\pm `$2.1)% and (40.7$`\pm `$1.1)%, respectively.

Since the tagging efficiencies depend on the jet energy, it is important to verify the agreement between the tagging efficiencies as a function of the transverse energy in the data and the simulation as is done in Figures 10 and 11. The distributions of the lifetime and invariant mass of the SECVTX tags is given in Figure 12 and supports our determination of the $`b`$-purity of the sample. The lifetime of a SECVTX tag is defined as

$`\mathrm{pseudo}\tau ={\displaystyle \frac{L_{xy}M^{SVX}}{cp_T^{SVX}}}`$

where $`M^{SVX}`$ and $`p_T^{SVX}`$ are the invariant mass and the transverse momentum of all tracks forming the SECVTX tag.

### A Check of the background parametrization using a photon conversion sample

In events where the e-jet does not contain heavy flavor, we predict the rate of tagged away-jets containing heavy flavor using the probability of tagging jets with heavy flavor as measured in generic-jet data. We test this method in a sample of data where the electrons in the e-jet are due to photon conversions. The criteria used to identify photon conversions are listed in Table II. In this case we require that an electron is matched by a second track consistent with a photon conversion and that it is not matched by a track segment in the VTX detector. Otherwise, we select this sample as the inclusive electron sample where in contrast conversions were removed.

Following the procedure used in the previous section, we determine the fraction of events with heavy flavor to be $`F_{hf}`$ = 8.7$`\pm `$0.9% from the number of e-jets with a SECVTX or JPB tag. Tagging rates in events due to heavy flavor production are described using the Herwig simulation as used in the previous section. In the remaining 91.3% of the events, we describe the rates of tagged away-jets using the parametrization derived from generic jets.

Table XVII shows that this procedure correctly predicts the rates of tags observed in the data. We take the 10% statistical error of this comparison as the systematic uncertainty of the method.

### B Sensitivity of the scale factor to the modeling of $`c`$-jets

In the simulation the tagging efficiency is defined as the ratio of events with double tags to all events where the away-jet is tagged by SECVTX. As shown in Table XIV, the Herwig simulation predicts that 94% of the a-jets with a SECVTX tag are due to $`b\overline{b}`$ production. The remaining 6% of the a-jets are due to $`c\overline{c}`$ production and are accounted for by the simulation but in principle this could be improperly modeled. In events where a-jets have a JPB tag, the fraction of $`c\overline{c}`$ production increases to 11% (see Table XIV). If SECVTX and JPB scale factors are determined using a-jets tagged by JPB instead of a-jets tagged by SECVTX, the scale factors change each by less than 1%. Therefore, we conclude that the modeling of $`c`$-jets is satisfactory for the determination of the $`b`$-tagging efficiency scale factor.

### C Dependence of the scale factor on the gluon splitting cross section

As shown in Table XIV, a fraction of the events in the inclusive electron sample is due to gluon splitting to heavy flavor quarks. The calibration of the Herwig simulation using generic-jet data in Section XI shows that direct production and the heavy flavor excitation as implemented in Herwig provide a fair description of the data, but the gluon splitting cross section requires a (40 $`\pm `$ 20)% correction. We repeat the calculation of the scale factor using this larger gluon splitting cross section. We find that the SECVTX scale factor increases from 1.23 to 1.25. The final scale factor we use will be this latter value.

### D  $`E_T`$ dependence of the scale factor

Jets produced directly in association with a $`\mathrm{W}`$ boson have transverse energies comparable to the jets in the low-$`p_T`$ inclusive electron sample. However, $`b`$-jets produced by top decay have substantially higher transverse energies. In this section, we investigate a possible $`E_T`$ dependence of the scale factor using two methods.

First, we derive the value of the SECVTX scale factor in four different bins of the electron-jet transverse energy. In each bin, we calculate the average e-jet transverse energy $`<E_T>`$ and the scale factor using the iterative procedure previously described. The result of the study is shown in Figure 13. A fit of the scale factor as a function of the transverse energy with a first order polynomial yields a $`\chi ^2`$ of 0.3 for 2 d.o.f. and

$`SF(E_T)=(1.23\pm 0.17)(0.1\pm 4.0)\times 10^3\times <E_T>(\mathrm{GeV})`$

with a correlation $`\rho =0.95`$ between the two fit parameters. The result of this fit is therefore consistent with a constant scale factor.

In the second method, we compare the fraction of jets with heavy flavor tagged by SECVTX in JET 50 and in JET 100 data and in the corresponding Herwig simulation tuned as in Section XI. The $`b`$-tagging efficiency in the detector simulation is increased by the factor 1.25 independently of the jet transverse energy. The ratio $`RSF`$ of the fractions of tagged jets in the data and in the simulation is sensitive to any residual $`E_T`$ dependence of the scale factor. The result of this method is shown in Figure 14. We fit the ratio $`RSF`$ of the tagging efficiencies in the data to the simulation as a function of the jet transverse energy with a first order polynomial. The fit yields a $`\chi ^2`$ of 51 for 49 d.o.f and

$`RSF(E_T)=(1.01\pm 0.05)+(1.3\pm 4.6)\times 10^4\times E_T(\mathrm{GeV})`$

with a correlation $`\rho =0.92`$ between the two fit parameters. The data are consistent with a constant scale factor.

### E Uncertainty of the scale factor

The SECVTX $`b`$-tagging efficiency scale factor measurement using the inclusive electron sample has a 5.6% uncertainty. The uncertainty in the gluon splitting cross section predicted by Herwig results in an additional 0.8% uncertainty of the scale factor. By folding the $`E_T`$ spectrum of $`b`$-jets from top decays with the $`E_T`$ parametrization of the scale factor from the fit shown in Figure 13 (a variation of the fit parameters by $`\pm 1\sigma `$ yields a $`\pm `$4.2% change in the efficiency for tagging $`b`$-jets and $`\pm `$3.9% change in the efficiency to tag top events), we estimate a 4% uncertainty from any residual $`E_T`$ dependence. These errors are mostly systematic and in general highly correlated. Altogether, we assign a 10% error to the determination of the scale factor after combining linearly the above contributions. Our final estimate of the $`b`$-tagging efficiency scale factor for the SECVTX algorithm is $`SF=1.25\pm 0.13`$ and for the jet-probability algorithm is $`SF=0.96\pm 0.10`$. The latter is consistent with unity.

### F Check of the scale factor using a low-$`p_T`$ inclusive muon sample

The low-$`p_T`$ inclusive muon sample is analogous to the electron sample in that a muon with $`p_T`$ 10 GeV/c is required in place of an electron with $`E_T`$ 10 GeV (see Section V B). It provides an independent sample for checking the tagging efficiency scale factor. The low-$`p_T`$ muon sample consists of 10393 events. In these events muon-jets without heavy flavor are due to fake muons arising from non-interacting hadrons or in-flight decays of $`K`$ and $`\pi `$ mesons. We compare to a simulated sample also generated using the option 1500 of Herwig which consists of 4280 events. The same procedure described above yields a SECVTX tagging efficiency scale factor of 1.24$`\pm `$0.10, in agreement with the value 1.23 $`\pm `$0.07 derived in the inclusive electron sample (before correcting the gluon splitting cross section). At the same time the heavy flavor purity of the low-$`p_T`$ muon sample is measured to be $`F_{hf}`$=59.7$`\pm `$3.6%.

### G Check of the scale factor in jets containing inclusive $`b`$-decays

In this section we investigate whether the scale factor is different in jets containing semileptonic $`b`$-decays and inclusive $`b`$-decays. We use the low-$`p_T`$ inclusive electron samples and normalize the data and the simulation to the same number of electron-jets with a SECVTX tag after mistag removal. In the simulation, the rate of gluon splitting to $`b\overline{b}`$ and $`c\overline{c}`$ pairs is corrected as in Section XI. We compare rates of away-jets which are taggable and which are tagged by SECVTX. We find that the simulation predicts correctly the amount of taggable away-jets but it underestimates by a factor 1.23 $`\pm `$ 0.08 the rate of SECVTX tags with respect to the data.

### H Check of the scale factor using rates of double tags in generic-jet data

The studies of the $`E_T`$ dependence of the SECVTX scale factor performed in Section IX D depends upon the assumption that Herwig models correctly the fractional yield of jets with heavy flavor as a function of their transverse energy. We use the JET 50 and JET 100 data and simulation for a test independent of this assumption. We select events with only two jets: one taggable jet with transverse energy larger than the trigger threshold and one taggable jet with $`E_T`$ 15 GeV in the opposite hemisphere. We compare the number of events with double JPB tags and double SECVTX tags in the data and in the Herwig simulation after mistag removal. In the simulation, 92% of these double tags are due to $`b\overline{b}`$ production. The ratio of double SECVTX to double JPB tags in the data and the simulation is:

$`R_{\mathrm{data}}=0.92\pm 0.18\text{ and }R_{\mathrm{sim}}=0.61\pm 0.05`$

This ratio does not depend on the absolute cross section for producing jets with heavy flavor. From the equivalence

$`{\displaystyle \frac{R_{data}}{R_{sim}}}=\left({\displaystyle \frac{SF^{SEC}}{SF^{JPB}}}\right)^2,`$

we measure $`\frac{SF^{SEC}}{SF^{JPB}}`$=1.24$`\pm `$0.13 using generic jets with high transverse energy, in agreement with the value $`\frac{SF^{SEC}}{SF^{JPB}}`$=1.28$`\pm `$0.10 measured in the low-$`p_T`$ inclusive electron sample.

### I SECVTX efficiency for tagging $`c`$-jets

Since we need to apply a large correction to the simulated SECVTX efficiency for tagging $`b`$-jets, it is worth investigating differences between data and simulation for tagging $`c`$-jets. For this purpose, we compare rates of tags in the JET 50 and JET 100 data to the corresponding Herwig simulation, described in Section XI, normalized to the same number of events.

We define $`R`$ as the ratio of the number of SECVTX to JPB tags after mistag removal. In the data $`R=0.77\pm 0.07`$. Under the assumption that the heavy flavor composition of the data is modeled correctly by Herwig, the SECVTX scale factor for $`c`$-jets, $`SF_c^{SEC}`$, can then be derived solving the equivalence

$`R={\displaystyle \frac{T_b^{SEC}SF_b^{SEC}+T_c^{SEC}SF_c^{SEC}}{T^{JPB}SF^{JPB}}}`$

where $`T_b^{SEC}`$=5354 and $`T_c^{SEC}`$=2477 are the number of simulated $`b`$ and $`c`$-jets tagged by SECVTX, and $`T^{JPB}`$=11958 is the number of JPB tags. Using $`SF_b^{SEC}`$=1.23 $`\pm `$ 0.07 and $`SF^{JPB}`$=0.96 $`\pm `$ 0.05, we derive that the SECVTX scale factor for tagging $`c`$-jets is $`SF_c`$ = 0.92 $`\pm `$ 0.28. The error is determined by the uncertainty of the heavy flavor composition (see Section XI) and by the errors of the scale factors $`SF_b^{SEC}`$ and $`SF^{JPB}`$.

### J Understanding of the scale factor

In an effort to explain the 25% difference of the SECVTX tagging efficiency in the data and the simulation we uncovered three oversights in the simulation package used in this and in some previous CDF analyses . A significant fraction of the difference is due to the use of an outdated version of the CLEO decay tables and to outdated $`B`$-lifetimes in the CDF particle database. The above two inaccuracies account for $`40`$% of the difference of the SECVTX scale factor from unity. Small inconsistencies in the implementation of the SVX geometry in the simulation contribute an additional 16% to this difference. If we corrected for these effects, the new determination of the SECVTX scale factor would be 1.09$`\pm `$0.11; the uncertainty includes the error on the $`b`$-lifetime ($``$3%) and the uncertainty of the track degradation procedure described in Section VII ($``$8%). The efficiency of jet-probability is not affected by these changes in the QFL simulation.

## X SECVTX and JPB mistags

In this section we estimate the SECVTX and JPB mistag rate in a variety of control samples before applying it to $`W+`$ jet and $`Z+`$ jet events in Sections XII and XIII. Tags in jets without heavy flavor, which we call mistags, are caused by detector resolution effects. SECVTX mistags are poorly reproduced by our detector simulation and traditionally CDF removed this background from the data using a parametrization of the probability of finding negative SECVTX tags in JET 50 data . We derive a new parametrization of the mistag rate using the JET 20, JET 50, JET 70, JET 100 and $`E_T`$ 300 data described in Section V A. Even if JPB mistags are well reproduced by the detector simulation, we derive a mistag parametrization also for JPB tags because jet-probability has a higher rate of mistags than SECVTX and provides a better check of the method.

The method to evaluate the mistag probability starts with the measurement of the number of positive and negative tags in generic-jet data and their parametrization as a function of the jet $`E_T`$ and the jet track multiplicity, $`N_{TRK}^{SVX}`$. The tagging probability is derived as a ratio of the number of tags to the number of taggable jets in bins of transverse energy and track multiplicity. We use only jets which are far away from calorimeter cracks and correct the jet energy for the detector response and out-of-cone losses (see Section III D).

Negative tags are also produced in jets containing heavy flavor. In particular, the probability of producing negative tags is different for jets initiated by a heavy-quark or by gluon splitting to a pair of heavy quarks. Since this contribution to negative tags must be accounted for and subtracted in order to obtain the mistag rate due to jets without heavy flavor, it is important to parametrize the rate of negative tags in a sample in which the composition of quark and gluon jets is well understood and is not subject to additional uncertainty due to the simulation. For this reason, in each generic-jet sample, we use only jets with transverse energy above the trigger threshold (leading jets): jets with corrected $`E_T`$ 30, 70, 90, 120 and 160 GeV in the JET 20, JET 50, JET 70, JET 100 and $`E_T`$ 300 data, respectively. In the generic-jet simulation, 95% of the leading jets with a tag contain just one heavy-flavored hadron (a large fraction of these leading jets is produced by heavy quarks from flavor excitation or direct production). The $`E_T`$ region below 30 GeV is mapped selecting events containing two leading jets, but using only the additional jets in the event ; in the simulation, 96% of these non-leading jets which are tagged contain two hadrons with heavy flavor produced by a gluon splitting process.

Transverse energy distributions of the jets used to measure the tagging probability are shown in Figures 15 and 16. Projections of the tagging probability matrices are shown in Figures 17 and 18.

Figure 19 shows that the tagging probability parametrization derived using jets with well measured energies works well for all jets.

Since the heavy flavor contribution to negative tags is expected to be small, the number of tags due to heavy flavor in a given $`E_T`$ bin of the tagging probability matrix is estimated fairly well by the difference, $`PN`$, between the numbers of positive (P) and negative (N) tags in this bin. In simulated jets with heavy flavor, we measure the ratio, $`R=\frac{N}{PN}`$, as a function of the jet transverse energy. We measure this ratio separately for jets which contain only one hadron with heavy flavor ($`R_1`$) and for jets which contain two hadrons with heavy flavor ($`R_2`$). The following empirical parametrization provides a good description of $`R`$ for jets containing $`b`$ as well as $`c`$-hadrons:

$$R_1(E_T)=\{\begin{array}{cc}0.0088+0.000158E_T\hfill & \text{for SECVTX}\hfill \\ 0.039+0.00117E_T\hfill & \text{for JPB}\hfill \end{array}$$

$$R_2(E_T)=\{\begin{array}{cc}0.075+0.000158E_T\hfill & \text{for SECVTX}\hfill \\ 0.14+0.00117E_T\hfill & \text{for JPB}\hfill \end{array}$$

With this parameterization we construct the mistag probability matrix by correcting each bin of the negative tagging probability matrix by the factors:

$`\mathrm{N}(\mathrm{P}\mathrm{N})R_1(E_T)`$ $`\text{ for jets with }E_T30\mathrm{GeV}`$
$`\mathrm{N}(\mathrm{P}\mathrm{N})R_2(E_T)`$ $`\text{ for jets with }E_T30\mathrm{GeV}`$
The fraction of negative tags contributed by heavy flavors is shown in Table XVIII.

In the generic-jet samples used to derive the mistag matrices, approximately 70% of the events contain additional interactions. The rate of multiple interactions is different in other samples, e.g. $`W+`$ multi-jet events where we require an isolated primary lepton. The negative tagging rate in the generic-jet data depends on the number of additional interactions.

Figure 20 shows the relative negative tagging probability, normalized to the average, as a function of the sum of the transverse momenta of all tracks associated with additional vertices displaced by more than 5 cm from the primary vertex $`p_T^V`$. Accordingly, the mistag rate is parametrized with the additional empirical function for both SECVTX and JPB:

$$F(p_T^V)=\{\begin{array}{cc}0.8+0.0128p_T^V\hfill & \text{for }p_T^V<60\text{ GeV/c}\hfill \\ 1.57\hfill & \text{for }p_T^V60\text{ GeV/c}\hfill \end{array}$$

### A Check of the SECVTX mistag parametrization

In this section, we test the capability of our model to predict the rate of negative tags in all available generic-jet samples.

Figure 21 serves to illustrate the procedure followed to predict the rates of negative tags. They are evaluated as the sum of the mistags plus the heavy flavor contribution using the $`R_1`$ and $`R_2`$ parametrizations derived in the previous section. This procedure requires the knowledge of the fraction of quark and gluon jets as a function of jet-$`E_T`$ in each data sample (literally, we need to know the fraction of jets containing one or two hadrons with heavy flavor). In the JET 20, JET 50, JET 70, JET 100, $`E_T`$ 175 and $`E_T`$ 300 samples we make the assumption, corroborated by the corresponding simulations, that all jets below trigger threshold are gluon jets and all jets above trigger threshold are quark jets.

Figure 21(a) shows the rate of observed positive tags and predicted mistags as a function of the jet $`E_T`$. Figure 21(b) compares rates of negative tags to the predicted mistags. The mistag rate does not include any heavy flavor contribution and is lower than the observed rate of negative tags. Figure 21(c) compares the rate of mistags and the heavy flavor contribution to the negative tags obtained by multiplying the difference between positive tags and predicted mistags in Figure 21(a) by $`R_1`$ ($`R_2`$) if the jet-$`E_T`$ is above (below) the trigger threshold. Figure 21(d) compares the observed and predicted yield of negative tags. The predicted yield is derived by adding the two distributions shown in Figure 21(c).

Following the same procedure, comparisons between the corrected jet $`E_T`$ distributions of observed and predicted negative SECVTX tags are shown in Figures 22 to 24. In the case of the $`E_T`$ 125 4CL sample, the ratio of quark to gluon jets (1/6, independent of $`E_T`$) is evaluated using the corresponding Herwig simulation. In the inclusive photon sample, we use only the $`R_1`$ parametrization as the simulation shows that tagged jets are mostly contributed by the $`\gamma c`$ Compton production.

The inclusive low-$`p_T`$ electron sample, used to measure the tagging efficiency scale factor, is also a good sample to test the validity of the $`R_1`$ and $`R_2`$ parametrizations because it is enriched in heavy flavor content. We compare rates of observed and predicted negative tags both in the data and the corresponding Herwig simulation. The fraction of gluon jets in the simulation is taken from Table XIV. In the data, this fraction is increased by 40% according to the calibration of the Herwig simulation performed in Section XI. Comparisons between observed and predicted rates of negative tags are shown in Figure 25 for the data and Figure 26 for the simulation.

Table XIX summarizes the rates of observed and predicted negative SECVTX tags in all generic-jet samples. Based on the observed agreement a 10% systematic error is assigned to the estimate of the SECVTX mistag probability.

### B Check of the JPB mistag parametrization

We follow the same procedure of the previous section to test the parametrization of the mistag rate of jet-probability. Figures 27 trough 29 compare $`E_T`$ distributions of observed and predicted jets with negative JPB tags for all generic-jet samples. Rates of JPB tags are summarized in Table XX. As before, by comparing the observed and predicted number of negative tags, we assign a 10% systematic error to the parametrization of the JPB mistags.

## XI Calibration of the fraction of $`\mathrm{W}+jet`$ events with heavy flavor

$`Wb\overline{b}`$ and $`Wc\overline{c}`$ events are produced through the so-called gluon splitting process, where a final state gluon branches into a heavy quark pair. In this analysis the fraction of $`W+`$ jet events containing heavy flavor is estimated using the Herwig generator. The uncertainty in the rate of gluons splitting into heavy quarks based on the parton shower approach is estimated to be approximately 40% in Ref.  and approximately 25% in Ref. . Because of this large uncertainty we calibrate the gluon splitting cross section modeled with Herwig using generic-jet data. Heavy flavor in generic-jet data stems from three primary sources: (1) direct production (e.g. $`ggb\overline{b}`$); (2) flavor excitation (e.g. $`gbgb`$); and (3) gluon splitting. The calibration of the simulation package is performed by tuning the various cross sections calculated by Herwig to reproduce the tagging rate observed in the JET 50 and JET 100 data. In these samples, the gluon splitting contribution is comparable to the other production mechanisms. In the JET 20 simulation, the gluon splitting contribution is negligible; we compare observed and predicted rates of tags in this sample using the tuned simulation as a check that we disentangled correctly the different heavy flavor production mechanisms.

In each generic-jet sample we count the number of SECVTX tags in taggable jets. Mistags are evaluated using the mistag probability evaluated in Section X.

The simulated samples (corresponding to the JET 20, JET 50 and JET 100 data) are generated using the option 1500 of Herwig and requiring hard scattering partons with $`|\eta |`$ 4.5 and $`p_T^{\mathrm{min}}`$ 10, 40 and 80 GeV/c respectively . We use the MRS(G) set of structure functions . Generated events are simulated with the standard package discussed in Section VII. As in the data, we select events containing at least one jet above the trigger threshold.

In the simulation a jet is classified as a $`b`$ or a $`c`$-jet if it contains a $`b`$ or a $`c`$-hadron in a cone of radius 0.4 around its axis. Hadrons with heavy flavor resulting from the fragmentation of one of the hard scattering partons are indicative of direct production or flavor excitation (if one of the incoming partons of the hard scattering has heavy flavor we attribute the process to flavor excitation; in this case a second hadron of the same flavor is produced by the backward-evolution of the structure functions). All pairs of hadrons with heavy flavor of the same type which do not come from the hadronization of the hard scattering partons are attributed to gluon splitting. Table XXI lists the rate of jets containing heavy flavor per event in the simulated JET 50 and JET 100 samples.

In the data, in order to distinguish gluon splitting from the rest of the heavy flavor production we use intuitive kinematical differences. Jets from heavy flavor direct production are expected to be produced back-to-back and are more likely to produce double tags. In events produced by heavy flavor excitation, jets produced by the backward-evolution of the structure functions tend to be at large pseudo-rapidities and out of the SVX acceptance. On the other hand, gluon splitting produces pairs of jets with heavy flavor at small separation $`\mathrm{\Delta }R=\sqrt{(\mathrm{\Delta }\varphi )^2+(\mathrm{\Delta }\eta )^2}`$. Most of the time the two hadrons with heavy flavor produced by gluon splitting reside in the same jet. Figure 30 shows distributions of the distance between two $`b`$-jets for the different production mechanisms in the simulated sample. In addition, the relative gluon splitting contribution increases with the jet multiplicity.

This motivates us to compare data and simulation in the following classes of SECVTX tags:

* (1): number of tagged jets per event with at least one taggable jet.
* (2): number of tagged jets per event with at least one taggable jet and with three or more jets with $`E_T`$ 15 GeV and $`|\eta |`$ 2.
* (3): twice the number of events with two tagged jets per event with two or more taggable jets.

We also compare the data to the simulation for:

* (4): the fraction of (1) in which the tagged jet has a companion jet with $`E_T10`$ GeV in a cone of radius 1.2 around its axis.
* (5): the fraction of events with double tags where the two tagged jets are at a distance $`\mathrm{\Delta }R`$1.2.

Table XXII lists the yields of tags in the data and in the simulation.

In the simulation, one notes that after tagging with SECVTX the contribution of $`c`$-jets is reduced by more than a factor of four and becomes negligible in events with double tags. However, the ratio of double to single SECVTX tags does not discriminate between $`b\overline{b}`$ and $`c\overline{c}`$ production for this ratio is similarly small for $`b\overline{b}`$ production through flavor excitation and gluon splitting.

We discriminate the flavor type with the additional comparison of rates of JPB tags (JPB has about the same tagging efficiency of SECVTX for $`b`$-jets and is more than twice as efficient for tagging $`c`$-jets). Since we use JPB tags only to disentangle between $`b`$ and $`c`$-production, we compare data and simulations in only two classes of JPB tags:

* (6): number of tagged jets per event with at least one taggable jet.
* (7): twice the number of events with two tagged jets per event with two or more taggable jets.

Table XXIII lists the yields of JPB tags in the data and in the simulation.

We proceed to fit the data with the simulation in order to extract the rate of $`gb\overline{b}`$ and $`gc\overline{c}`$. When fitting the simulation to the data, the yield of simulated SECVTX and JPB tags is corrected for the tagging efficiency scale factors measured in Section IX. The 10% uncertainty in the scale factor determination is included in the error of the simulated rates of tags. In the fit, we also compare five distributions in each generic-jet sample and in the corresponding simulation:

1. the yield of the fraction of SECVTX tags per taggable jet as a function of the jet-$`E_T`$.
2. The distributions of the distance $`\mathrm{\Delta }R`$ between a jet tagged by SECVTX and a companion jet as defined above.
3. The distributions of the distance $`\mathrm{\Delta }R`$ between a jet tagged by JPB and a companion jet as defined above.
4. The distributions of the distance $`\mathrm{\Delta }R`$ between two jets tagged by SECVTX.
5. The distributions of the distance $`\mathrm{\Delta }R`$ between two jets tagged by JPB.

In the comparison, the area of each distribution is normalized to unity. For each distribution we compute a reduced $`\chi ^2`$

$`\chi _D^2={\displaystyle \frac{1}{N}}{\displaystyle \underset{i=1}{\overset{N}{}}}{\displaystyle \frac{(d(i)sd(i))^2}{ed(i)^2+esd(i)^2}}`$

where $`N`$ is the number of bins in each distribution, $`d(i)`$ and $`sd(i)`$ are the contents of the bin $`i`$ of the distribution in the data and the simulation respectively, with $`ed(i)`$ and $`esd(i)`$ their errors. The simulated jet-$`E_T`$ distributions have a systematic uncertainty due to the trigger simulation which is cumbersome to account for in the fit. Simulated distribution of distances between tagged jets have unaccounted systematic uncertainties due to how well the parton shower generator models gluon splitting at distances $`\mathrm{\Delta }R`$ 1.2. We use the reduced $`\chi _D`$ to diminish the importance of these comparisons with respect to the classes of absolute tagging rates. The data are fitted to the simulation using a minimum $`\chi ^2`$ method. We minimize the function

$`\chi ^2={\displaystyle \underset{JET50}{\overset{JET100}{}}}({\displaystyle \underset{j=1}{\overset{7}{}}}{\displaystyle \frac{(F(j)SF(j))^2}{EF(j)^2+ESF(j)^2}}+{\displaystyle \underset{k=1}{\overset{5}{}}}\chi _D^k)`$

where the index $`k`$ runs over the 5 kinematic distributions described in the previous paragraph, $`F(j)`$ are the yields of tags observed in the data for the classes listed in Tables XXII and XXIII, and

$`SF(j)={\displaystyle \underset{n=1}{\overset{6}{}}}P(n)C_H(j,n)SF^\alpha `$

is the corresponding yield of simulated tags. The contributions $`C_H(j,n)`$ of different flavor types and production mechanisms, as listed in Tables XXII and XXIII, are weighted with the fit parameters $`P(n)`$. $`SF`$ is the tagging efficiency scale factor and $`\alpha `$=0 for $`c`$-jet, 1 for events with one tagged $`b`$-jet and 2 for events with two tagged $`b`$-jets.

In the fit, the $`b`$-to-$`c`$ ratio for direct production is constrained to the default value with a 14% Gaussian error. Option 1500 of Herwig evaluates the direct production cross section of massless quarks. The 14% uncertainty accounts for having neglected the quark masses (estimated using option 1700 of Herwig) and for the uncertainty in the fragmentation process (estimated using the Pythia generator).

The $`b`$-to-$`c`$ ratio for flavor excitation is also constrained to the default value with a 28% Gaussian error. This uncertainty accounts for the largest variation of this ratio observed using a wide range of structure functions in the PDF library .

The ratio of the $`gb\overline{b}`$ to $`gc\overline{c}`$ is also constrained to the default value with a 28% Gaussian error. The uncertainty accounts for a $`\pm `$0.5 GeV change of the $`b`$ and $`c`$-quark masses around the default value.

The fit has 21 degrees of freedom and yields a $`\chi ^2`$ of 22. The fit results are shown in Table XXIV. The weights of the gluon splitting cross sections will be used to rescale the fraction of $`W+`$ jet events with heavy flavor predicted by Herwig.

Figure 31 compares the $`E_T`$ distributions of tagged jets in the data and in the fitted simulation. Similarly, Figures 32 and 33 compare distributions of distances between tagged jets.

Table XXV compares rates of tags in generic-jet data and in the Herwig simulation calibrated according to Table XXIV. The JET 20 sample was not used to calibrate the simulation package. Similarly, the SLT algorithm was not used in the Herwig calibration. The comparison of the number of SLT tags in the data and simulation serves to check independently the calibration of the Herwig production cross section and the SLT tagging efficiency in the simulation.

## XII Composition of the $`W+`$1 jet sample

The background to the $`t\overline{t}`$ production is determined using the data or the simulation calibrated as described in Sections IX and XI. The $`t\overline{t}`$ production cross-section is determined by attributing the excess of tagged $`W+`$ 3 jet events to $`t\overline{t}`$ production. $`\mathrm{W}+1`$ and $`\mathrm{W}+2`$ jet events provide a check of the background calculation. The evaluation of the backgrounds to $`t\overline{t}`$ events is detailed in the subsections A-H. The results of these background determinations are listed in the tables in Section XV where the cross section is calculated. The following two Sections XIII and XIV provide checks of the background estimates.

### A non-$`W`$ background

As in previous analyses , the background from non-$`W`$ sources, including $`b\overline{b}`$ production, is determined directly from the data by studying the isolation of primary lepton candidates in the low ($`\mathrm{}\mathrm{E}_\mathrm{T}`$ 10 GeV) and in the high ($`\mathrm{}\mathrm{E}_\mathrm{T}`$ 20 GeV) $`\mathrm{}\mathrm{E}_\mathrm{T}`$ region. The number of non-$`W`$ events in each jet-bin is evaluated as

$`N_{nonW}=N_C\times {\displaystyle \frac{N_A}{N_B}}`$

where $`N_A`$, $`N_B`$ and $`N_C`$ are the number of events in regions A, B and C of Figure 34. The corresponding number of tagged events is

$`N_{\mathrm{non}\mathrm{W}}^{\mathrm{tag}}=N_{\mathrm{non}\mathrm{W}}\times P^{\mathrm{tag}}`$

where $`P^{\mathrm{tag}}`$ is the tagging probability measured in region A. The yield of $`P^{\mathrm{tag}}`$ as a function of the lepton isolation is shown in Figure 35.

### B  $`Z\tau ^+\tau ^{}`$ events

The $`Z\tau ^+\tau ^{}`$ contribution is estimated using the Pythia generator (option MSEL=13). The simulation is normalized to the same number of $`Z\mu ^+\mu ^{},e^+e^{}`$ events observed in the data.

### C Single top production

The single top contribution via the $`W`$gluon fusion channel is estimated with Herwig using the process 2000. The single top production for the annihilation process $`q\overline{q}W^{}t\overline{b}`$ is estimated using the Pythia generator (option MSEL=12). We use the cross sections $`\sigma _{Wg}`$=1.5$`\pm `$0.4 pb for $`Wg`$ fusion derived using the NLO calculation in Ref.  and $`\sigma _{W^{}t\overline{b}}`$=0.74$`\pm `$0.05 pb for the annihilation process .

### D Diboson production

The contribution of the $`ZZ`$, $`WZ`$ and $`WW`$ production is estimated using the Pythia generator (options MSEL=15 and ISUB=22, 23 and 25, respectively). We use the cross sections $`\sigma (WW)`$=9.5$`\pm `$0.7 pb, $`\sigma (WZ)`$=2.6$`\pm `$0.34 pb and $`\sigma (ZZ)`$=1.0$`\pm `$0.2 pb .

### E Mistags

The SECVTX and SLT mistags are calculated weighting each jet in the $`W`$ sample with the mistag probability matrices derived in Section X and Section VIII C 1, respectively. The reevaluation of the SECVTX mistag matrix has resulted in a reduced estimate of this background in the signal region by (50 $`\pm `$ 5)% compared with the previous estimates of Ref. .

For the jet-probability algorithm, each simulated background also includes the contribution of mistags. The number of JPB mistags is evaluated only for the fraction of $`W+`$ jet events which is not simulated, i.e. $`W+`$ jet direct production without heavy flavor.

### F  The $`W+b\overline{b}`$ and $`W+c\overline{c}`$ contribution

We use the Herwig generator (process 2100 with hard scattering $`p_T^{\mathrm{min}}`$ 10 GeV/c) to estimate the fraction of $`W+`$ n jet events, $`F_1^\alpha `$, in which only one jet contains hadrons with heavy flavor resulting from gluon splitting ($`\alpha `$ refers to the flavor type). The fraction of $`W+`$ 2 jet events, $`F_2^\alpha `$, in which two different jets contain hadrons with heavy flavor is calculated using the Vecbos generator (see Section VI). The fractions of $`W+`$ jet events with heavy flavor content are listed in Table XXVI. We use the Herwig and Vecbos simulations also to determine the efficiency for finding events with one or two tagged jets, as listed in Table XXVII.

It follows that the number of tagged $`Wb\overline{b}`$ and $`Wc\overline{c}`$ events is

$`N_{1\mathrm{t}\mathrm{a}\mathrm{g}}^\alpha `$ $`=`$ $`N_W(F_1^\alpha \epsilon _{1\mathrm{t}\mathrm{a}\mathrm{g}}^{1\alpha }+F_2^\alpha \epsilon _{1\mathrm{t}\mathrm{a}\mathrm{g}}^{2\alpha })`$
$`N_{2\mathrm{t}\mathrm{a}\mathrm{g}}^\alpha `$ $`=`$ $`N_WF_2^\alpha \epsilon _{2\mathrm{t}\mathrm{a}\mathrm{g}}^{2\alpha }`$
where $`N_W`$ is the number of $`W`$ events in the data after removing the predicted number of non-$`W`$, di-bosons, single top, unidentified $`Z`$ and $`t\overline{t}`$ events.

### G The $`Wc`$ contribution

The fraction $`F_{Wc}`$ of $`gsWc`$ and $`gdWc`$ events is evaluated using the Herwig $`W+`$ 1 jet simulation and is shown in Table XXVIII. The estimated uncertainty on $`F_{Wc}`$ is dominated by the uncertainty in the strange sea content in the proton, computed by examining a wide range of different structure functions in Ref. . The tagging efficiencies for this process are listed in Table XXIX.

### H Direct production of $`Z+`$ jet with heavy flavor

We use the Pythia generator (option MSEL=13) to estimate the fraction of unidentified $`Z+`$ jet events passing our selection. The simulation is normalized to the number of $`Z\mathrm{}\mathrm{}`$ observed in the data. The fraction of $`Z+`$ jet events containing heavy flavor is determined by rescaling the Herwig prediction for the fraction of $`W+`$ jet events with heavy flavor. The rescaling factor is derived by comparing $`W+`$ 1 jet and $`Z+`$ 1 jet events generated with Pythia (option MSEL=14 and 13, respectively). For the same number of generated events, we find $`0.3\pm 0.15`$ $`Zc`$ events for every $`Wc`$ event, $`1\pm 0.3`$ $`Zc\overline{c}`$ events for every $`Wc\overline{c}`$ event, and $`2.0\pm 0.5`$ $`Zb\overline{b}`$ events for every $`Wb\overline{b}`$ event.

## XIII Check of the background calculation using the $`Z+1`$ jet sample

The production mechanisms of $`W`$ and $`Z`$ bosons in association with jets are very similar. The contribution of $`t\overline{t}`$ events to the $`Z+`$ jet events is negligible. This sample provides a good benchmark for our background calculation. The selection of the $`Z+`$1 jet event sample is described in Section IV A. Table XXX shows the predicted composition of the $`Z+`$1 jet sample before tagging. In this Table the number of $`Z+`$ jet events is derived from the data by subtracting the $`WW`$, $`WZ`$, $`ZZ`$, $`t\overline{t}`$, and single top contribution. The measured and predicted rates of events with SECVTX, JPB and SLT tags are shown in Tables XXXI to XXXIII. The product of the four Poisson probabilities, one per jet bin, for the observed numbers given the prediction is $`P_0=1.2\times 10^3`$ for Table XXXI, $`P_0=2.1\times 10^4`$ for Table XXXII, and $`P_0=1.0\times 10^3`$ for Table XXXIII. With a Monte Carlo simulation, in which we fluctuate the predicted rates by their uncertainty according to a gaussian distribution, we estimate that the likelihood of observing a probability no larger than $`P_0`$ is 33.8% for events with SECVTX tags, 17.9% for events with JPB tags and 41.1% for events with SLT tags. In $`Z+`$ jet events the background prediction agrees with the data.

## XIV Rates of negative tags in the $`W+`$ 1 jet sample

As shown in Section X, the mistag rates plus the estimated heavy flavor contribution to the negative tags account for the observed rates of negative tags in all generic-jet data. A similar test in the $`W+`$1 jet sample offers an additional check of the mistag rate predictions and a complementary test of the method used to estimate the background contribution to the $`t\overline{t}`$ signal.

The rate of negative tags for each process is calculated from the corresponding simulation or using the data as we do for positive tags. We use the sample composition before tagging listed in Tables XXXVI and XXXVIII for SECVTX and JPB, respectively. Table XXXIV compares numbers of observed and predicted negative SECVTX tags as a function of the jet multiplicity. The analogous comparison for negative JPB tags is shown in Table XXXV. Data and predictions agree within the estimated uncertainties.

## XV Measurement of the $`t\overline{t}`$ production cross section

The $`t\overline{t}`$ production cross section is derived using

$`\sigma _{t\overline{t}}={\displaystyle \frac{N_{\mathrm{tag}}^{\mathrm{obs}}N_{\mathrm{tag}}^{\mathrm{bkg}}}{𝒜_{t\overline{t}}ϵ_{\mathrm{tag}}𝑑t}}`$

where $`N_{\mathrm{tag}}^{\mathrm{obs}}`$ is the number of tagged $`W+`$ 3 jet events, $`N_{\mathrm{tag}}^{\mathrm{bkg}}`$ is the background prediction, $`𝒜_{t\overline{t}}`$ is the detector acceptance for $`t\overline{t}`$ events, $`ϵ_{\mathrm{tag}}`$ is the efficiency for tagging top events and $`𝑑t`$ = 105.1 $`\pm `$ 4.0 pb-1 is the total integrated luminosity.

The acceptance for $`t\overline{t}`$ events is evaluated with a simulation which uses the Pythia generator and is (7.8$`\pm `$1.3)% for a top mass of 175 $`\mathrm{GeV}/\mathrm{c}^2`$. The 17% systematic error accounts for all uncertainties in the simulation which come from the following: lepton identification and trigger simulation ($`\pm `$ 14%), jet energy scale ($`\pm `$ 5%), modeling of initial state gluon radiation ($`\pm `$ 2%), final state gluon radiation ($`\pm `$ 5%), Monte Carlo modeling of the $`t\overline{t}`$ production ($`\pm `$ 5%), detector resolution effects ($`\pm `$ 2%), and instantaneous luminosity dependence ($`\pm `$ 2%).

The tagging efficiencies are evaluated using the same simulation and are $`0.505\pm 0.051`$, $`0.455\pm 0.046`$ and $`0.157\pm 0.016`$ for SECVTX, JPB and SLT respectively.

In the background calculation the rate of $`W+`$ jet events with heavy flavor is estimated from the number of events due to $`W+`$ jet direct production using the fraction of heavy flavor determined in Section XII F. Therefore the contribution of $`t\overline{t}`$ events must be removed from the data. This is done by iterating. The $`t\overline{t}`$ cross section is estimated from the excess of tagged $`W+`$ 3 jet events over the background calculated assuming $`\sigma _{t\overline{t}}=0`$. The resulting $`\sigma _{t\overline{t}}`$ is used to evaluate the number of $`t\overline{t}`$ events before tagging; this contribution is subtracted from the data to obtain the contribution of the $`W+`$ jet direct production. The amount of $`W+`$ jet with heavy flavor is recalculated and $`\sigma _{t\overline{t}}`$ is updated. The procedure is repeated until $`\sigma _{t\overline{t}}`$ is stable to within 1%.

In the sample of 252 $`W+`$ 3 jet events, there are 29 events with at least one jet tagged by the SECVTX algorithm. Using the procedure described above, the background estimate is 8.0 $`\pm `$ 1.0 events. Assuming that all the excess is due to $`t\overline{t}`$ production, the resulting $`t\overline{t}`$ cross section is 5.08 $`\pm `$ 1.54 pb (the statistical error is $`\pm `$ 1.30 pb and the systematic $`\pm `$ 0.82 pb). The estimated break-down of the $`W+`$ 1 jet sample before and after tagging is shown in Tables XXXVI and XXXVII.

As a cross-check, we calculate $`\sigma _{t\overline{t}}`$ using rates of JPB tags. There are 41 events with at least one jet tagged by the JPB algorithm with a background of 11.1 $`\pm `$ 1.3 events. The observed excess of events yields $`\sigma _{t\overline{t}}`$ = 8.02 $`\pm `$ 2.16 pb. The estimated break-down of the $`W+`$ 1 jet sample before and after tagging is shown in Tables XXXVIII and XXXIX.

There are 25 events with at least one jet tagged by the SLT algorithm with a background of 13.2 $`\pm `$ 1.2 events. The observed excess of events yields $`\sigma _{t\overline{t}}`$ = 9.18 $`\pm `$ 4.26 pb (the statistical error is $`\pm `$ 3.89 pb and the systematic $`\pm `$ 1.72 pb). The estimated break-down of the $`W+`$ 1 jet sample before and after tagging is shown in Tables XL and XLI.

Because of a small dependence of the acceptance and the tagging efficiencies on the top mass, for a $`\pm `$5 $`\mathrm{GeV}/\mathrm{c}^2`$ change in the top mass the cross sections calculated using SECVTX and JPB tags change by $`\pm `$1.8% and the cross section calculated using SLT tags changes by $`\pm `$2.3%.

## XVI Combined $`t\overline{t}`$ production cross section

The best measurement of the $`t\overline{t}`$ cross section comes from combining the results of this analysis with the dilepton and all-hadronic analyses . The revised $`b`$-tagging efficiency reported in this paper effects the cross section measured in the all-hadronic channel. The details of this analysis have not changed from those reported in Ref. . The cross section measurement from the dilepton channel  does not require $`b`$-tagging information and so is unchanged by the revisions reported here. It is affected slightly by the revised determination of the total integrated luminosity as were all measurements. A comparison between the previously published results and the revised cross sections used for the new combined result is shown in Table XLII.

We combine the measurements from the SVX and SLT tagged lepton + jets, all-hadronic, and dilepton channels, using a maximum likelihood technique similar to that described in Refs. . This procedure properly accounts for correlated systematic uncertainties, such as the uncertainties in the $`b`$-tagging efficiency, the luminosity, the kinematic acceptance, and some of the calculated backgrounds. In all channels the acceptances are calculated with $`M_{top}=175`$ $`\mathrm{GeV}/\mathrm{c}^2`$. The resulting combined $`t\overline{t}`$ production cross section is

$`\sigma _{t\overline{t}}=6.5_{1.4}^{+1.7}\mathrm{pb}`$

where the quoted uncertainties include both statistical and systematic effects, which are approximately equal in magnitude.

## XVII Conclusions

Having improved the method for determining the $`b`$-tagging efficiency, and the method for calculating the backgrounds to $`t\overline{t}`$ production, we revise our previous measurements of $`\sigma _{t\overline{t}}`$ in the lepton + jets channel . We find $`\sigma _{t\overline{t}}=5.08\pm 1.54\mathrm{pb}`$ and $`\sigma _{t\overline{t}}=9.18\pm 4.26\mathrm{pb}`$ using events with SECVTX and SLT tags, respectively. We have used the jet-probability algorithm as a cross-check and find that it gives a result consistent with these measurements. The measurement of the $`t\overline{t}`$ cross section, obtained by combining the results of this analysis with the dilepton and all-hadronic analyses, is $`\sigma _{t\overline{t}}=6.5_{1.4}^{+1.7}\mathrm{pb}`$, in agreement with the SM predictions  and the measurement performed by the DØ collaboration .

## ACKNOWLEDGMENTS

We thank the Fermilab staff and the technical staff of the participating Institutions for their contributions. This work was supported by the U.S. Department of Energy and National Science Foundation; the Istituto Nazionale di Fisica Nucleare; the Ministry of Education and Culture of Japan; the Natural Science and Engineering Research Council of Canada; the National Science Council of the republic of China; the Swiss National Science Foundation; the A.P. Sloan Foundation; the Bundesministerium fuer Bildung und Forschung; and the Korea Science and Engineering Foundation.