Online (trigger) Selection
Events are recorded by an online selection algorithm which requires the presence of 1 or 2 electron candidates, a muon candidate, or missing transverse energy (>30 GeV) and jets.
Lepton Identification
We select muon and electron candidates using neural network lepton identification algorithms.
Z Candidate Selection
We require a dilepton pair with a mass between 76 and 106 GeV/c2. The dilepton mass for all candidate events is shown below.
Jet Selection
We require events to have 2 or 3 jets with a transverse energy of greater than 25 GeV and a pseudo-rapidity (η) of -2.0 < η < 2.0. The dijet mass of the two lead-ET jets is required to be greater than 25 GeV/c2. The dijet mass for all candidate events is shown below.
Jet Energy Corrections
We correct the energy of the two lead-ET jets using an artificial neural network algorithm. The dijet mass for all candidate events after NN correction is shown below.
Data Model
Major backgrounds to the ZH search are Z+2 light flavor (lf) jets, Z+cc, Z+bb, tt, and diboson processes. We model the contributions from Z+lf,Z+cc, and Z+bb using an ALPGEN simulation. Contributions from WW,WZ,ZZ, and tt are modeled using a PYTHIA simulation. We use data-derived methods to model the contribution from misidentified Z→l+l- decays.
PreTag Control Region
We validate our data model in the sample of reconstructed Z→l+l- + 2 or 3 jet events before requiring the presence of one or more identified b-jets (PreTag). In general we find reasonable agreement between our model and the observed PreTag data. The distributions of various event quantities are shown here.
b-jet Identification
We use the new CDF HOBIT algorithm to identify (tag) b-jets. We define four exclusive analysis categories based on the b-quark content of an event : events with two 'Tight' Hobit b-tags (TT), events with one Tight and one 'Loose' Hobit b-tags (TL), events with a single Tight tag (Tx), and events with two Loose tags (LL).
3-stage Neural Network Signal Discrimination
We separate ZH from background processes by using a series of neural networks. The first neural network is trained to separate tt from ZH signal and can be seen here. The second neural network is trained to separate Z+jet and Z+charm backgrounds from ZH signal and can be seen here. The third neural network separates diboson backgrounds from ZH signal and can be seen here. Based on the values returned by the three expert networks, we assign candidates to one of 4 regions in our final discriminant. The 3-stage sorting logic is summarized in the flowchart below. In each of the 4 regions of the final discriminant, we plot the output of a neural network designed to simultaneously separate ZH from all backgrounds. We deploy 13 versions of this network, each optimized for a specific Higgs boson mass. The output of this network (as optimized for a 120 GeV Higgs boson) can be seen before and after sorting. The sorted discriminant is used in the extraction of upper limits on σZH×BR(H→bb).
Systematics
The following table summarizes the systematic uncertainties affecting the shape and/or normalization of the components of our background model.
| Systematic | Effect | Affected Samples |
|---|---|---|
| Luminosity σinel(pp) | ± 3.8% | all simulated |
| Luminosity Monitor | ± 4.4% | all simulated |
| Lepton ID Efficiency | ± 1% | all simulated |
| Lepton Energy Scale | ± 1.5% | all simulated |
| Mis-reconstructed Z→l+l- Rate |
± 50% (electrons) ± 5% (muons) |
Mis-reconstructed Z Model |
| udsg-jet tight mistag rate | ± 19 to 40% | Z+lf samples |
| udsg-jet loose mistag rate | ± 10 to 20% | Z+lf samples |
| Jet Energy Scale (JES) | on the order of ± 10% includes variation in shape as seen here. |
Z+lf samples |
| tight b-tag rate | ± 3.9 to 7.8% | Z+cc,Z+bb
tt, diboson, ZH |
| loose b-tag rate | ± 3.2 to 6.3% | Z+cc, Z+bb
tt, diboson, ZH |
tt cross section | ± 10% | tt | diboson cross section | ± 6% | WW,WZ,ZZ |
| Z+cc and Z+bb cross section | ± 40% | Z+cc, Z+bb | ZH cross section | ± 5% | ZH | Amount of initial or final state radiation | ± 1 to 15% | ZH |
| Trigger Efficiency | ± 5% (muons) ± 1% (electrons) |
all simulated |
Event Totals
The predicted event totals are compared to the observed data for each of the 16 sub-channels in the following tables :
Z→e+e- + 2 jets
| Contribution (9.45/fb) | TT | TL | Tx | LL |
|---|---|---|---|---|
| tt | 20.1 ± 2.8 | 21.5 ± 2.8 | 36.1 ± 4.7 | 6.1 ± 0.8 |
| WW,WZ,ZZ | 4.7 ± 0.6 | 6.5 ± 0.9 | 19.6 ± 1.8 | 3.9 ± 0.4 |
| Z+bb | 19.1 ± 8.0 | 26.8 ± 11.3 | 81.5 ± 34.2 | 10.2 ± 4.4 |
| Z+cc | 1.5 ± 0.6 | 6.9 ± 2.9 | 39.0 ± 16.8 | 7.3 ± 3.1 |
| Z+lf jets | 0.7 ± 0.3 | 8.3 ± 2.0 | 124.9 ± 27.5 | 27.5 ± 6.6 |
| Mis-reconstructed Z | 0.1 ± 0.0 | 5.1 ± 2.6 | 7.7 ± 3.9 | 1.1 ± 0.6 |
| Total Bkg. | 46.2 ± 8.6 | 75.2 ± 12.4 | 309.2 ± 47.4 | 56.1 ± 8.6 |
| Data | 45 | 83 | 352 | 66 |
| ZH (120 GeV) | 1.1 ± 0.1 | 1.1 ± 0.1 | 1.6 ± 0.2 | 0.3 ± 0.03 |
Z→e+e- + 3 jets
| Contribution (9.45/fb) | TT | TL | Tx | LL |
|---|---|---|---|---|
| tt | 7.5 ± 1.2 | 9.3 ± 1.4 | 13.5 ± 1.9 | 2.9 ± 0.5 |
| WW,WZ,ZZ | 0.7 ± 0.1 | 1.3 ± 0.2 | 3.0 ± 0.4 | 1.0 ± 0.1 |
| Z+bb | 4.5 ± 2.0 | 6.5 ± 2.9 | 14.1 ± 6.2 | 2.5 ± 1.1 |
| Z+cc | 0.5 ± 0.2 | 1.7 ± 0.8 | 7.4 ± 3.3 | 2.4 ± 1.1 |
| Z+lf jets | 0.3 ± 0.1 | 2.8 ± 0.8 | 20.3 ± 5.5 | 8.1 ± 2.3 |
| Mis-reconstructed Z | 0.0 ± 0.0 | 2.1 ± 1.0 | 5.2 ± 2.6 | 3.0 ± 1.5 |
| Total Bkg. | 13.6 ± 2.3 | 23.6 ± 3.5 | 63.5 ± 9.5 | 19.9 ± 3.2 |
| Data | 16 | 23 | 59 | 23 |
| ZH (120 GeV) | 0.2 ± 0.04 | 0.2 ± 0.04 | 0.3 ± 0.1 | 0.1 ± 0.01 |
Z→μ +μ - + 2 jets
| Contribution (9.45/fb) | TT | TL | Tx | LL |
|---|---|---|---|---|
| tt | 20.8 ± 3.1 | 22.1 ± 3.1 | 30.4 ± 3.9 | 5.7 ± 0.8 |
| WW,WZ,ZZ | 3.8 ± 0.6 | 5.1 ± 0.7 | 15.1 ± 1.5 | 3.0 ± 0.4 |
| Z+bb | 15.0 ± 6.3 | 21.0 ± 8.8 | 64.4 ± 27.0 | 7.7 ± 3.2 |
| Z+cc | 1.0 ± 0.4 | 4.6 ± 2.0 | 30.0 ± 12.6 | 6.3 ± 2.6 |
| Z+lf jets | 0.6 ± 0.3 | 6.2 ± 1.5 | 91.7 ± 20.2 | 19.4 ± 4.5 |
| Mis-reconstructed Z | 1.0 ± 0.1 | 0.0 ± 0.0 | 10.0 ± 0.5 | 1.0 ± 0.1 |
| Total Bkg. | 42.3 ± 7.1 | 58.9 ± 9.7 | 241.5 ± 36.3 | 43.0 ± 6.2 |
| Data | 41 | 69 | 273 | 51 |
| ZH (120 GeV) | 0.9 ± 0.1 | 0.9 ± 0.1 | 1.4 ± 0.1 | 0.3 ± 0.03 |
Z→μ +μ - + 3 jets
| Contribution (9.45/fb) | TT | TL | Tx | LL |
|---|---|---|---|---|
| tt | 6.4 ± 1.2 | 7.4 ± 1.2 | 10.4 ± 1.5 | 2.4 ± 0.4 |
| WW,WZ,ZZ | 0.6 ± 0.1 | 0.9 ± 0.2 | 2.3 ± 0.3 | 0.8 ± 0.1 |
| Z+bb | 3.5 ± 1.5 | 5.2 ± 2.4 | 11.3 ± 5.0 | 2.3 ± 1.1 |
| Z+cc | 0.4 ± 0.2 | 1.5 ± 0.7 | 5.8 ± 2.5 | 1.9 ± 0.8 |
| Z+lf jets | 0.3 ± 0.1 | 2.2 ± 0.6 | 15.3 ± 4.0 | 6.3 ± 1.7 |
| Mis-reconstructed Z | 1.0 ± 0.1 | 8.0 ± 0.4 | 8.0 ± 0.4 | 5.0 ± 0.3 |
| Total Bkg. | 12.2 ± 1.9 | 25.2 ± 2.8 | 53.0 ± 7.0 | 18.8 ± 2.2 |
| Data | 15 | 24 | 46 | 25 |
| ZH (120 GeV) | 0.2 ± 0.03 | 0.2 ± 0.04 | 0.2 ± 0.05 | 0.1 ± 0.01 |
Results
We compute Bayesian 95% CL upper limits on the value of σZH×BR(H→bb), and obtain expected upper limits from the median of the upper limits resulting from a series of background-only pseudo-experiments. The observed and expected limits are listed below and are shown in the figure here.
| Higgs Mass (GeV/c2) |
Observed/SM | -2σ/SM | -1σ/SM | Median Exp./SM | +1σ/SM | +2σ/SM |
|---|---|---|---|---|---|---|
| 90 | 1.0 | 1.0 | 1.4 | 1.9 | 2.8 | 3.8 |
| 95 | 1.2 | 1.0 | 1.4 | 2.1 | 3.0 | 4.5 |
| 100 | 1.8 | 1.0 | 1.4 | 2.1 | 3.0 | 4.2 |
| 105 | 2.3 | 1.1 | 1.5 | 2.2 | 3.1 | 4.4 |
| 110 | 3.0 | 1.2 | 1.6 | 2.4 | 3.4 | 4.8 |
| 115 | 4.7 | 1.4 | 1.8 | 2.6 | 3.7 | 5.3 |
| 120 | 5.7 | 1.5 | 2.1 | 3.1 | 4.4 | 6.3 |
| 125 | 7.2 | 1.9 | 2.5 | 3.6 | 5.2 | 7.4 |
| 130 | 10.8 | 2.5 | 3.3 | 4.8 | 7.0 | 10.1 |
| 135 | 15.0 | 3.2 | 4.3 | 6.3 | 9.1 | 12.6 |
| 140 | 19.1 | 4.4 | 6.3 | 8.8 | 12.9 | 18.7 |
| 145 | 21.7 | 6.7 | 9.2 | 13.3 | 19.5 | 27.1 |
| 150 | 36.5 | 10.7 | 15.0 | 21.3 | 30.4 | 44.6 |
The dilepton mass for all candidate events:

The dijet mass for all candidate events:

The dijet mass for all candidate events after NN correction:

Output of the tt expert discriminant:

Output of the lf+charm expert discriminant:

Output of the diboson expert discriminant:

3-stage Neural Network Signal Discrimination:

120 GeV Optimized Network Output Before Sorting :

120 GeV Optimized Network Output After Sorting :

120 GeV Optimized Network Output for Total Model Under JES Variation :

95% CL Upper Limits :
