Higgs Boson Searches at the Tevatron
Harald Fox Department of Physics h.fox@lancaster.ac.uk
Higgs Boson Searches at the Tevatron Harald Fox Department of - - PowerPoint PPT Presentation
Higgs Boson Searches at the Tevatron Harald Fox Department of Physics h.fox@lancaster.ac.uk Contents evatron, D and CDF T Higgs production Search for heavy Higgs H WW ll b jet Search for light Higgs WH l bb ZH
Harald Fox Department of Physics h.fox@lancaster.ac.uk
T evatron, DØ and CDF Higgs production Search for heavy Higgs H → WW → llνν Search for light Higgs WH → lνbb̅ ZH → ννbb̅ ZH → llbb̅ Outlook Conclusion
W + W − W +
Kane, Scientific American, June 2003
vacuum expectation value
field and acquire an effective mass
V(Ф)=μ2|Ф|+λ(|Ф|2)2
γ, W and Z bosons is determined
ratios are determined. v =
2λ = 246GeV
mγ = mW = 1 2vg mZ = 1 2vg 1 cos θW mH = √ 2λv2 mf = 1 √ 2gfv
1 2 3 4 5 6 100 30 300
mH [GeV] ∆χ2
Excluded
Preliminary
∆αhad = ∆α(5)
0.02758±0.00035 0.02749±0.00012
Theory uncertainty
mLimit = 144 GeV
mH < 144 GeV 95%C.L. LEP EWWG
Kolda, Murayama: JHEP 0007 (2000) 035
100 200 300 400 500 600 1 10 10
2
Higgs mass (GeV) (TeV)
Vacuum Stability Triviality Electroweak 10% 1%
Fine Tuning ΔFT < 10%/1%
Run IIa Run IIb
Delivered Recorded Run IIa 1.6 fb-1 1.3 fb-1 Run IIb (so far) 1.9 fb-1 1.7 fb-1 Total 3.5 fb-1 3.0 fb-1
2006 shutdown:
April 02 Jan 08 Passed 3fb-1 milestone in recorded luminosity on 16 January 2008
Two General Purpose Detectors: CDF DØ Electron acceptance |η|<2.0 |η|<3.0 Muon acceptance |η|<1.5 |η|<2.0 Silicon Precision tracking |η|<2.0 |η|<3.0 Hermetic Calorimeter |η|<3.6 |η|<4.2 Powerful trigger systems (2.5MHz →50Hz) Dilepton triggers with pT>4GeV
protons antiprotons 3 Layer Muon System Tracker Solenoid Magnet
The Higgs cross section is 10-11 orders of magnitudes lower than the total inelastic cross section.
Evidence of single top production is an important milestone towards the Higgs boson. Light quarks are ubiquitous. Plenty of W and Z bosons → calibration. Total inelastic cross section.
cross section (pb)
ℓ ℓ ℓ ℓ
W+ e+ W- e- ν n
Angular correlation of leptons due to V−A as H is a spin 0 particle:
(e,e) φ ∆
0.5 1 1.5 2 2.5 3 3.5 4
entries
10 1 10
210
310
410
510
(e,e) φ ∆
0.5 1 1.5 2 2.5 3 3.5 4
entries
10 1 10
210
310
410
510
data 10 × WW → 160 H e e → Z Diboson γ W+jets/ QCD ttbare+e−
L=1.2fb−1
DØ Run II Preliminary
ℓ ℓ ℓ ℓ
MET / GeV 10 20 30 40 50 60 70 80 90 100
10
10 1 10
210
310
410 MET / GeV 10 20 30 40 50 60 70 80 90 100
10
10 1 10
210
310
410
=160
H
After Preselection, M
WZ t t QCD WW
ZZ
µ µ
Data H160
S/B ≈ 15/300k
}(250-500)/ GeV
µ µM 20 40 60 80 100 120 140 160 180 200 220 240
10
10 1 10
210 / GeV
µ µM 20 40 60 80 100 120 140 160 180 200 220 240
10
10 1 10
210
Before Cut
QCD t t ZZ WZ WW
µ µ
Data H120
S/B ≈ 5/50
50 100 150 200 250 300 350 400 450 500
10 1 10 HT / GeV 50 100 150 200 250 300 350 400 450 500
10 1 10
Before Cut
QCD t t ZZ WZ µ µ
WW Data H160
µ ,
1µ (
1 1.5 2 2.5 3
10
10 1 10 )
2µ ,
1µ (
1 1.5 2 2.5 3
10
10 1 10
=160
H
After HT Cut, M
QCD t t ZZ WZ µ µ
WW Data H160
Selection criterion mH = 120 mH = 140 mH = 160 mH = 180 mH = 200 Cut 1 Preselection Trigger, ID, leptons with opposite charge, zV T X < 60 cm, Mµµ > 17 GeV pT > 20/10GeV 20/15 25/15 25/15 25/15 Cut 2 Missing trans- verse energy E / T 25 < E / T < 70 25 < E / T < 80 30 < E / T < 90 35 < E / T < 100 35 < E / T < 110 Cut 3 Sig(E / T ) Sig(E / T ) > 5 (for NJet > 0) Cut 4 M T
min (l, E
/ T ) M T
min > 30
M T
min > 30
M T
min > 40
M T
min > 45
M T
min > 45
Cut 5 Invariant mass Mµµ 17 < Mµµ < 60 17 < Mµµ < 70 17 < Mµµ < 75 17 < Mµµ < 85 17 < Mµµ < 95 Cut 6 ΣpT = pl
T + pl T +
E / T 60 < ΣpT < 135 70 < ΣpT < 160 80 < ΣpT < 170 90 < ΣpT < 180 90 < ΣpT < 200 Cut 7 HT (scalar sum of pJet
T )
HT < 60 HT < 60 HT < 60 HT < 60 HT < 50 Info Neural Net NN > 0.5
e+e−
L=1.2fb−1
DØ Run II Preliminary
[GeV] T miss E 20 40 60 80 100 120 140 160 180 200 entriese+e−
L=1.2fb−1
DØ Run II Preliminary
[GeV] T min M 20 40 60 80 100 120 entriese+e−
L=1.2fb−1
DØ Run II Preliminary
(e,e) φ ∆ 0.5 1 1.5 2 2.5 3 3.5 4 entriese+e−
L=1.2fb−1
DØ Run II Preliminary
[GeV] T p∑
50 100 150 200 250 300 entries∑
50 100 150 200 250 300 entriese+e−
L=1.2fb−1
DØ Run II Preliminary
NN
NN
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
entries
10 1 10
210
310
NN
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
entries
10 1 10
210
310 e+e−
L=1.2fb−1
DØ Run II Preliminary
NN Output ΣpT Mll E̷T MTmin (l,E̷T) Δϕll ≈30% improvement from NN
mH=160 x 10
data 10 × WW →
160
H
e e → Z Diboson γ W+jets/ QCD ttbar
e+e−
L=1.2fb−1
DØ Run II Preliminary
[GeV] T miss E 20 40 60 80 100 120 140 160 180 200 entriese+e−
L=1.2fb−1
DØ Run II Preliminary
[GeV] T min M 20 40 60 80 100 120 entriese+e−
L=1.2fb−1
DØ Run II Preliminary
(e,e) φ ∆ 0.5 1 1.5 2 2.5 3 3.5 4 entriese+e−
L=1.2fb−1
DØ Run II Preliminary
[GeV] T p∑
50 100 150 200 250 300 entries∑
50 100 150 200 250 300 entriese+e−
L=1.2fb−1
DØ Run II Preliminary
NN
NN
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
entries
10 1 10
210
310
NN
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
entries
10 1 10
210
310 e+e−
L=1.2fb−1
DØ Run II Preliminary
NN Output ΣpT Mll E̷T MTmin (l,E̷T) Δϕll
pt1 pt2 met dphiMetL1 dphiMetL2 M minMt dphiL1L2 type
rain against more backgrounds
nnout
0.2 0.4 0.6 0.8 1 1.2 1.4
10
10 1 10
210
Data (551) =160) h Signal (M Signal+Bkgd. Total Bkgd. Z+jets Bkgd. Diboson Bkgd. W+jets Bkgd. QCD Bkgd. Top Bkgd.µ+µ−
L=1.2fb−1
DØ Run II Preliminary
NN Output
Improved NN
:
ME dscr.
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
entries
10 1 10
210
310
ME dscr.
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
entries
10 1 10
210
310 e+e−
L=1.2fb−1
DØ Run II Preliminary
ME discr.
Matrix Element Additional input to NN
Run IIa Run IIb
NN Output
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1Events
NN Output
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1Events
NN Output NN Output
0.2 0.4 0.6 0.8 1 1.2Events
NN Output
Events
NN Output
0.2 0.4 0.6 0.8 1 1.2Events
NN Output
Events
ee eμ μμ Combination of
Uncertainties
LEP CLs Method
NN Output NN Output NN Output
Events
DØ Preliminary, L=1.20 fb , RunIIb ν µ , ν e →
W → H
(d)
NN Output1-CLb CLsb LLR = −2 ln Q CLs = 1 − CLs+b CLb Q( s, b, d) =
NChan
Nbins
(s + b) dij
ij
e (s + b)ij dij! / b d
ij e b
dij! Systematics taken into account via Gaussian marginalisation Correlations taken into account
0.02 0.04 0.06 0.08 0.1 0.12 0.14
5 10 15
Probability density
Observed Expected for background Expected for signal (mH=115.6 GeV/c2) + background
LEP
15
250 500 750 1000(a)
(b)
(c) Figure 1. Left: The pdfs of the combined Higgs search at LEP for the background (right) and signal + background hypotheses (left) for mH = 115.6 GeV/c2. The light grey region to the left
Right: Illustration of the evolution of the pdfs with falling search sensitivity from (a) to (c) as the Higgs mass hypothesis is increased and the production cross-section falls.
200
)
2
(GeV/c
H
m
120 130 140 150 160 170 180 190 200
LLR
1 2 3 4
σ 2-
b
LLR σ 1-
b
LLR
b
LLR
s+b
LLR
LLR
DØ Preliminary, L=2.3 fb
+
W → H
LLR
)
2
(GeV/c m
120 130 140 150 160 170 180 190 200
)
+
W → BR(H × H) → p (p σ Limit /
1 10
DØ Preliminary, L=2.3 fb
+
W → H
Observed Limit Expected Limit σ 1- ± Expected σ 2- ± Expected
Standard Model = 1.0
LLR
exp: 2.4
Limits relative to SM expectation!
)
2
(GeV/c
H
m 110 120 130 140 150 160 170 180 190 200 (SM) ! 95% C.L. Limit / 1 10
2
10
: 1.1 fb b llb " ZH
: 2.1 fb b b # # " ZH
: 1.7 fb b b # l " WH
: 1.1 fb
+
W W " WH
: 2.3 fb
+
W " H
: 2.3 fb $ $ " H DZero Combination
Standard Model = 1.0
DØ Preliminary, L=1.1-2.3 fb (SM) ! 95% C.L. Expected Limit /
signal ee eµ µµ etrk µtrk
signal
electrons muons
signal NN output
signal
s Higgs at Tevatron 1
exp: 2.4
ZH → l+l- bb 2 b jets ~ 1/2 MH each 2 leptons ~ 45 GeV each Z mass constraint Cleanest signal WH → lνbb 2 b jets ~ 1/2 MH each 1 lepton ~ 50 GeV each Missing ET ~ 50 GeV Highest production X-sec ZH → νν bb 2 b jets ~ 1/2 MH each 0 leptons Missing ET ~ 100 GeV Largest expected signal
lepton
frag
PV jet
B hadron
K
frag
B h a d r
K π D SV b b
Lxy
Fake Rate (%) 2 4 6 8 10 12 b-Jet Efficiency (%) 20 30 40 50 60 70 80
46 % 28 %Tagger NN JLIP
p
Z bb
T
P
20 40 60 80 100 120 140 160 180
Events
1000 2000 3000 4000
T
P
20 40 60 80 100 120 140 160 180
Events
1000 2000 3000 4000
DØ Preliminary
L = 1.7 fb
W + 2 jets
Data W + jets QCD SM bkgd
T
P
20 40 60 80 100 120 140 160 180
Events
1000 2000 3000 4000 (a)
2 20 30 40 50 60
46 % 28 %
Fake Rate % b-Jet Efficiency %
4 different analyses:
pT 2nd jet Δϕ jets pT di-jet pT (l,ETmiss)
Dijet Mass (GeV)
50 100 150 200 250 300
Events
20 40 60
Dijet Mass (GeV)
50 100 150 200 250 300
Events
20 40 60
DØ Preliminary
L = 1.7 fb
W + 2 jets / 2 b-tags
Data W + jets QCD t t b Wb
WH
115 GeV (x10)Dijet Mass (GeV)
50 100 150 200 250 300
Events
20 40 60
(d)R
1 1.5 2 2.5 3 3.5 4 4.5 5
Events
10 20 30 40
R
1 1.5 2 2.5 3 3.5 4 4.5 5
Events
10 20 30 40
DØ Preliminary
L = 1.7 fb
W + 2 jets / 2 b-tags
Data W + jets QCD t t b Wb
WH
115 GeV (x10)R
1 1.5 2 2.5 3 3.5 4 4.5 5
Events
10 20 30 40
(b)T
P
50 100 150 200 250
Events
50 100 150
T
P
50 100 150 200 250
Events
50 100 150
DØ Preliminary
L = 1.7 fb
W + 2 jets / 2 b-tags
Data W + jets QCD t t b Wb
WH
115 GeV (x10)T
P
50 100 150 200 250
Events
50 100 150
(a)NN output - 2 tags
0.2 0.4 0.6 0.8 1 1.2 1.4
Events
1 10
210
310
NN output - 2 tags
0.2 0.4 0.6 0.8 1 1.2 1.4
Events
1 10
210
310
NN output - 2 tags
0.2 0.4 0.6 0.8 1 1.2 1.4
Events
1 10
210
310
DØ Preliminary
L = 1.7 fb
W + 2 jets / 2 b-tags
Data W + jets QCD t t b Wb
WH
115 GeV (x10)NN output - 2 tags
0.2 0.4 0.6 0.8 1 1.2 1.4
Events
1 10
210
310
WH115 x10
NN
pT leading jet ΔR jets m(jet1,jet2)
single b-tag double b-tag
Dijet Mass (GeV) 50 100 150 200 250 300 Events 1 10 2 10 3 10 Dijet Mass (GeV) 50 100 150 200 250 300 Events 1 10 2 10 3 10 Dijet Mass (GeV) 50 100 150 200 250 300 Events 1 10 2 10 3 10 DØ Preliminaryexp: 7.3
CDF
)
2(GeV/c
Hm 105 110 115 120 125 130 135 140 145 ) b b
(p
10 20 30 40 50 60 70
DØ Preliminary, L=1.7 fb b b
Observed Limit Expected Limit
(a)
exp: 9.05
Limits relative to SM expectation
b jet b jet
ETmiss + 2 jets b jet b jet MET 1 tight b-tag + 1 loose b-tag S/B ~ 3.7/443
DiJet Invariant Mass (GeV)
50 100 150 200 250 300Events / 12.00 GeV
0.5 1 1.5 2 2.5 3 3.5 4 4.5 3 10)
DØ preliminary (2.1 fb DiJet Invariant Mass (GeV)
50 100 150 200 250 300Events / 12.00 GeV
10 20 30 40 50 60 50 100 150 200 250 300 10 20 30 40 50 60 Data Top Z+b/c-jets Z+jets(l.f.) W+b/c-jets W+jets(l.f.) Diboson Multijet Hx10 (115 GeV))
DØ preliminary (2.1 fb DT discriminant
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1Events / 0.05
10 20 30 40 50DT discriminant
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1Events / 0.05
10 20 30 40 50 Data Top Z+b/c-jets Z+jets(l.f.) W+b/c-jets W+jets(l.f.) Diboson Multijet VHx25 (115 GeV))
DØ preliminary (1.2 fb
Decision Tree (Run IIb)
Boosting: as single top (adaptive boosting, AdaBoost) give mis-classified events a higher weight before re- training to make the tree work harder . Asymmetric b-tagging: ϵ=73%/48%; f=5%/0.5% (@pT>30, η<0.8)
/ ET
pT jets
scalar ET pT jet1 / HT pT jet2 HT ηjet1 A(/ ET , / HT ) . = (/ ET − / HT )/(/ ET + / HT ) ηjet2 / HT /HT dijet invariant mass min ∆φ(/ ET , jets) dijet transverse mass ∆φ(/ ET , jet1) (/ ET − pT tracks)/(/ ET + pT tracks)) ∆φ(/ ET , jet2) max(∆φ(/ ET , jets)) − min(∆φ(/ ET , jets)) ∆φ(/ ET , pT jet1 + pT jet2) max(∆φ(/ ET , jets)) + min(∆φ(/ ET , jets)) ∆φ(jet1, jet2)
pT tracks
∆R(jet1, jet2) pT tracks from dijets ( pT tracks − dijets pT tracks)/ pT tracks
Table 15: Variables used as input to the Decision Tree
H H H H H H H H HT
T T T T T T T T>212
>212 >212 >212 >212 >212 >212 >212 >212 P F P F p p p p p p p p pt
t t t t t t t t<31.6
<31.6 <31.6 <31.6 <31.6 <31.6 <31.6 <31.6 <31.6 P F M M M M M M M M Mt
t t t t t t t t<352
<352 <352 <352 <352 <352 <352 <352 <352 purity purity purity purity purity purity purity purity purity
hes, XLIII Recontres de Moriond (QCD)
ℓ ed from data between
2 b-tagging requirements: both jets with secondary vertex tag 1 jet with SVT , 1 jet with low probability that all tracks originate from the primary vertex 2 separate NN:
jets and un-correlated to tracks; track based quantities enter the NN
, m(jj), E̷T(cal), met-dot-product: E̷T(cal) · E̷T(trk), dR(jj)
)
2= 115 GeV/c
HNeural Network Output (m
0.2 0.4 0.6 0.8 1 5 10 15 20 25 30 35)
2= 115 GeV/c
HNeural Network Output (m
0.2 0.4 0.6 0.8 1 5 10 15 20 25 30 35Vertex + Probability Tag (Signal Region)
))
2Neural Network Output (115 GeV/c
0.2 0.4 0.6 0.8 1 5 10 15 20 25)
2Neural Network Output (115 GeV/c
0.2 0.4 0.6 0.8 1 5 10 15 20 25Double Vertex Tag (Signal Region)
))
2= 115 GeV/c
HNeural Network Output (m
0.2 0.4 0.6 0.8 1 5 10 15 20 25 30 35
)
2= 115 GeV/c
HNeural Network Output (m
0.2 0.4 0.6 0.8 1 5 10 15 20 25 30 35
Vertex + Probability Tag (Signal Region)
)
CDF Run II Preliminary (1.7 fb
Ttbar W+h.f. Z+h.f Single Top WZ/WW ZZ Mistags QCD Bckgnd Err VH*8 (115 GeV) Data)
2
(GeV/c
H
M
110 115 120 125 130 135 140 145 150 155
95% CL Upper Limit/SM
1 10 )
CDF Run II Preliminary (1.7 fb Observed Limit σ 1 ± Expected Limit
Met+Jets Search for ZH/WH
exp: 8.3
(GeV)
H
m 105 110 115 120 125 130 135 140 145 ) b b
(p
5 10 15 20 25 30 35 40
)
DØ Preliminary (2.1 fb DT, VH Signal b b
Observed Limit Expected Limit
exp: 8.4
Upcoming improvements:
7
)
2
(GeV/c
H
m 105 110 115 120 125 130 135 140 145 ) b b → BR(H × ZH) → p (p σ Limit / 10 20 30 40 50 60 70 80 90 100
DØ Preliminary, L=1.1 fb b ll b → ZH
Observed Limit Expected Limit
exp: 20
Mee (GeV)
50 100 150 200 250
Events / 2 GeV
10
10 1 10
210
310 50 100 150 200 250
10
10 1 10
210
310
Data QCD Z+jets Z+2b tt ZZ WZ ZH 115
)
DØ Preliminary (920 pb
Leading-Pt di-jet mass (GeV)
40 60 80 100 120 140 160Number of Events / 10 GeV
DØ Preliminary Leading-Pt di-jet mass (GeV)
40 60 80 100 120 140 160Number of Events / 10 GeV
1 10 2 10 3 10 40 60 80 100 120 140 160 1 10 2 10 3 10 Data QCD Z+jets Z+bb(cc) tt WZ ZZ ZH 115DØ Preliminary
)
Neural Network output
Number of Events / 0.05
Data QCD Z+jets Z+bb(cc) tt WZ ZZ ZH 115 DØ Preliminary
No b-tag Double b-tag Neural Net S/B ~ 0.53/74 before NN
exp: 16
Two independent neural nets are trained to separate ZH from
Z+Jets Vs. ZH 0.2 0.4 0.6 0.8 1 t Z H V s . t 0.2 0.4 0.6 0.8 1 Number of Events 5 10 15 20
Ldt = 1 fb
∫
CDF II Preliminary NN Output - Data (Single Tag)
Ldt = 1 fb
∫
CDF II Preliminary NN Output - Data (Single Tag)
Z+Jets
t t̅
ZH
Use τlepτhad mode.
3 Neural Nets are trained: Signal vs Z-> ττ + jets Signal vs ttbar Signal vs QCD Select Minimum of 3 NN to fit data.
(GeV)
γ γ
M
50 100 150 200 250 300 350 400 450 500
Events/5 GeV
1 10
210
310 50 100 150 200 250 300 350 400 450 500 1 10
210
310
data γ γ QCD j γ jj *->ee γ Z/ signal(M=130GeV)
SM Higgs mass (GeV)
100 110 120 130 140 150
SM value ) γ γ BR( × σ 95% CL
20 40 60 80 100 120 140 160
Observed Limit Expected Limit σ 1 ± Expected Limit σ 2 ± Expected Limit preliminary
DØ, 2.27 fb
2 isolated em clusters QCD and γj background estimated from data
1 10 10 2 110 120 130 140 150 160 170 180 190 200 1 10 10 2 mH(GeV/c2) 95% CL Limit/SM
Tevatron Run II Preliminary, L=1.0-2.4 fb-1 D∅ Exp CDF Exp Tevatron Expected Tevatron Observed ±1σ ±2σ
LEP Limit SM
March 2, 2008
exp: 3.3
exp: 1.6
37
Improvements assumed in projections
✦ b-tagging
✦
Acceptance
✦
Analysis techniques
✦ scaling of systematic uncertainties as a function of luminosity
Additional improvements not yet included in projection
inclusion of tau channels charm rejection in single b-tag analyses
…
=160 GeV
Summer 2005 Channels Summer 2006 Channels Summer 2007 Channels Winter 2008 Channels With Improvements
2 4 6 8 10 12 14 16 1 2 3 4 5 6 7 Integrated luminosity/Experiment (fb-1) Expected Limit/SM
Summer 2005 Channels Summer 2006 Channels arXiv:0712.2383 (2007) With CDF Update Winter 2008 With Improvements
!"#$%&%'%&()*+,-.'"+"#&$ #$%&'(
115 GeV CDF & DØ combined Rob Roser, P5 Meeting, 01/02/08
2009 2010
Sensitivity factors Minimum = x1.5 Further = x2.25
CDF+D0 combined
95% CL