Measurement of the ttbar Cross Section in the Lepton + Jets Channel 

Using Neural Network in 4.6 fb-1 of CDF Data

Including Ratio Over the Z Cross Section


A. Lister,
Université de Genève  

J. Conway, R. Erbacher, W. Johnson, T. Schwaz,

UC Davis

K. Lannon,
University of Notre Dame

R. Hughes, B. Winer
The Ohio State University

CDF public conference note (not yet updated since 2.8 fb-1)


NN Fit

ttbar Cross Section Results for 4.6 fb-1

σttbar = 7.63 ± 0.37(stat)± 0.35 ± 0.15 (theory) pb.

Note: All results assume Mttbar = 172.5 GeV/c2 unless otherwise stated.





Abstract

We present a measurement of the top pair production cross section  in ppbar collisions at 1.96 TeV, with an integrated luminosity of  4.6 fb-1 at the CDF experiment on the Fermilab Tevatron. We use a neural network technique to discriminate between top pair production and background processes in a sample of events with an isolated, energetic lepton, large missing transverse energy and three or more energetic jets. We then significantly reduce the dependence on the luminosity measurement and its associated large systematic uncertainty. We compute the ratio of the ttbar to Z cross section, measured using the same triggers and dataset, and then multiply this ratio by the theoretical Z cross section. In essence we replace the luminosity uncertainty with the uncertainty on the theoretical Z cross section.
We measure a top pair production cross section of 
σttbar =  7.63 ± 0.37 (stat) ± 0.35 (sys) ± 0.15 (Z theory) pb for a top mass of 172.5 GeV/c2
The total uncertainty is 7.0%, greatly surpassing the Tevatron Run II goal of 10%, and now more precise than the best theoretical calculations.

The top cross section was measured for 2 additional assumed top masses: 170 and 175 GeV/c2.
The systematics were not re-computed but were scaled from the 172.5 GeV/c2 measurement

σttbar (Mttbar = 170 GeV/c2) =  8.33 ± 0.40 (stat) ± 0.39 (sys) ± 0.17 (Z theory) pb  
σttbar (Mttbar = 175 GeV/c2) =  7.29 ± 0.35 (stat) ± 0.34 (sys) ± 0.14 (Z theory) pb.

Note that for this analysis, the central value is quoted after re-weighting to the CTEQ6.6 central PDF.




Event Selection

We use the standard CDF lepton + jets selection
  • Central Lepton (electron or muon) pT >= 20 GeV/c;
  • Missing transverse energy >= 20 GeV;
  • At least 3 jets with ET >= 20 GeV. 
 We apply additional cuts
  • Tighter missing transverse energy cut at 35 GeV;
  • Tighter leading jet ET cut at 35 GeV.
The tighter cuts were optimised to remove the most QCD background while maintaining reasonable efficiency for the ttbar signal. As the statistical uncertainty is not our dominant uncertainty, we can afford to cut quite hard. These cuts remove ~85% of the QCD, ~40% of the W+jets while maintaining ~80% efficiency for the ttbar signal.

Both central electrons and muons have been used for the final fit.
 
Note: We do not require any b-tagging.

Signal and Background Modeling

  • The ttbar signal is modeled from PYTHIA Monte Carlo with an assumed mass of 172.5 GeV/c2.
  • The W+jets background is used to model all EWK backgrounds; the previous analysis showed that the effect of including all backgrounds was very small and such a difference is included as a source of systematic uncertainty.
  • The W+jets background is modeled from ALPGEN+PYTHA MC. The W+jets sample is obtained from a combination of W+0p, W+1p, W+2p, W+3p exclusive as well as W+4p inclusive.
  • The QCD shape is modeled from a sample of di-jet events which pass our event selection criteria. These events are dominated by jets that fake electrons. The QCD contamination in the electron sample is significantly larger than in the muon sample






Neural Network Input Variables

The Neural Network uses 7 kinematic distributions as inputs, with 1 hidden node.
The NN is trained to separate W+4p from ttbar Monte Carlo. It was not retrained for this iteration of the analysis.
The input variables are
Global event variables
  • Σjets ET of all jets excluding two leading jets
  • HT of the event (sum of transverse energy of all reconstructed objects)
  • Aplanarity of the event
Variables related to the 3 leading jets
  • Σjets pz / Σjets ET
  • Minimum di-jet mass
  • Minimum di-jet separation
  • Maximum jet ET

Ratio over the Z Cross Section

The dominant systematic uncertainty of the direct measurement of the ttbar cross section is in the luminosity measurement, 5.8%. 
We can almost entirely cancel this uncertainty by considering the ratio of the ttbar to the Z cross sections.
The Z cross section is measured using the same triggers ( central electron and central muon) and the same data sample as for the ttbar cross section.
For this measurement, the signal MC for ttbar and Z is re-weighted to the CTEQ6.6 PDF (with its associated uncertainties).
The measured ttbar cross section, from the NN fit, is found to be
σttbar = 7.52) ± 0.36 (stat) ) ± 0.34 (sys) ± 0.44 (lumi) pb.
The Z cross section is measured to be
σz =  247.79 ± 0.79 (stat) ± 4.38 (sys) ± 14.59 (lumi) pb.
The ratio of the ttbar to Z cross section is computed, taking into account the correlations between the systematics
We multiply the ratio by the theoretical Z cross section
σZtheory = 251.3 ± 5.0 (sys) pb. The dominant uncertainties on the theoretical Z cross section are the scale and the PDFs.
We get a final result for the measured top pair production cross section of
σttbar = 7.63 ± 0.37 (stat)± 0.35 (sys)± 0.15 (theory) pb.
The total uncertainty is 7.0%, a significant reduction on the 8.6% uncertainty obtained using the NN kinematic fit only.

Final Fits to NN Output Distribution

Fits for the 4.6 fb-1 data sample follow bellow.

The input QCD normalisation is obtained from a fit to the missing ET distribution before the missing ET cut is applied.
The uncertainty on the constraint is set to 50% in the final fit.
The ttbar and W+jets shapes float freely in the fit.  

NNfit

NN Output distribution showing the signal and background contributions obtained from the fit to data.
The fit returns 1718 reconstructed ttbar event (after selection).

Plots of NN Input Variable

These plots are shown using the signal and background normalisations obtained from the final fit.
Note that for these plots the overflow and underflow bins are not shown, hence the slightly different number of events)
Sum Et jets 3,4,5
Σ ET of jets excluding 2 leading jets

Ht HT of the event

Sum Pz / Sum ET Σ pz / Σ ET of 3 leading jets

Aplanarity Aplanarity of the Event

minimum dijet separation Minimum dijet separation of 3 leading jets

minimum dijet mass Minimum dijet mass of 3 leading jets

Maximum jet eta Maximum jet eta of 3 leading jets

Fit to Missing Transverse Energy to get QCD input normalisation

Sum Et jets 3,4,5
This plot shows the fit to the Missing transverse momentum, before the cut on that variable is applied.
The SM cross section is assumed for ttbar, QCD and W+jets are floated freely in the fit.
The top plots shows the templates used in the fit, the bottom shows the data compared to the output of the fit.
The Missing transverse energy cut is then applied and the insert box shows the expected number of QCD events that would pass this cut.

Systematics

The main sources of systematic uncertainties in this measurement is due to the jet energy scale and the systematics associated to the modeling of the ttbar signal in Monte Carlo. In particular the difference between ALPGEN+PYTHIA vs ALPGEN+HERWIG and the difference when considering the uncertainty on the branching ratio of the W boson from the top decay (which we take as the PDG uncertainty).
systematics
Last modified on 13 August 09 by Alison Lister