A Search For ZH→l+l-bb
in 9.45/fb of CDF II Data

Link To Analysis Documentation Key Figures For Conference Speakers

A brief outline of the search techniques is presented below. For more details please see the analysis documentation available here.


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:

high res. image

The dijet mass for all candidate events:

high res. image

The dijet mass for all candidate events after NN correction:

high res. image

Output of the tt expert discriminant:

high res. image

Output of the lf+charm expert discriminant:

high res. image

Output of the diboson expert discriminant:

high res. image

3-stage Neural Network Signal Discrimination:

high res. image

120 GeV Optimized Network Output Before Sorting :

high res. image

120 GeV Optimized Network Output After Sorting :

high res. image

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

high res. image

95% CL Upper Limits :

high res. image