Collider physics II: Jets
1
Tuesday, September 11, 12
Collider physics II: Jets 1 Tuesday, September 11, 12 Two aspects - - PowerPoint PPT Presentation
Collider physics II: Jets 1 Tuesday, September 11, 12 Two aspects of new developments Better QCD jet. - Smarter jet algorithm. - Noise suppression with jet grooming. Jet substructure. - Boosted top. - Higgs. Boston Jet Workshop:
1
Tuesday, September 11, 12
Boston Jet Workshop: http://jets.physics.harvard.edu/workshop/Main.html Northwest Terascale workshop http://www.physics.uoregon.edu/~soper/Jets2011/talks.html Boost 2011, May, 23-27, Princeton. http://boost2011.org Boost 2012, July, 23-27, Valencia, Spain. http://ific.uv.es/boost2012/
2
Tuesday, September 11, 12
Madgraph, Alpgen, ...
Pythia, Herwig, Sherpa, ...
PGS, Delphes, “by hand”.
Fastjet. SpartyJet
http://www.lpthe.jussieu.fr/~salam/fastjet/ http://projects.hepforge.org/spartyjet/
3
Tuesday, September 11, 12
Many of them only have hadronic channels.
˜ q ˜ q∗ ¯ q q ... ...
˜ q ˜ q∗ ¯ q q .... .... ¯ q q ˜ g ˜ g
jet jet jet jet jet jet
4
Tuesday, September 11, 12
Jets with substructure.
q q top q ... q... W ± top b b
q, ... q, ... W ± q ... q... W ± ... q, ... ¯ q, ... Z q, ... ¯ q, ... Z
5
Tuesday, September 11, 12
LHC jet: 50 GeV - several TeV Tevatron jet: 50 - 100s GeV
6
Tuesday, September 11, 12
Drawing: F. Krauss
7
Tuesday, September 11, 12
Soft or collinear radiation should not be able to induce “large” changes in the observable. Otherwise, we cannot compare calculation with
8
Tuesday, September 11, 12
9
Tuesday, September 11, 12
9
Hard interaction time (distance) scale Q-1
Tuesday, September 11, 12
9
Hard interaction time (distance) scale Q-1 “talking” to the rest of the proton time(distance) scale mproton-1
Tuesday, September 11, 12
9
Hard interaction time (distance) scale Q-1 “talking” to the rest of the proton time(distance) scale mproton-1
If Q-1 ≪ mproton-1 (hard interaction) Two processes should not affect each other → Factorization!
Tuesday, September 11, 12
9
Hard interaction time (distance) scale Q-1 “talking” to the rest of the proton time(distance) scale mproton-1
If Q-1 ≪ mproton-1 (hard interaction) Two processes should not affect each other → Factorization!
A similar story for final states fragmentation, q,g⇒ hadrons (pion, K...)
Tuesday, September 11, 12
2
Parton densities Threshold matrix elements
×phase-space
Partonic cross section
10
Tuesday, September 11, 12
11
Hard interaction, gg⇒ g h t tbar ⇒h t tbar decay PDF PDF
Tuesday, September 11, 12
12
PDF PDF Clusters of hadronic energy final state object: jet pjet = Σ p of constituents Inclusive: independent of final states, just energy
Tuesday, September 11, 12
12
PDF PDF Clusters of hadronic energy final state object: jet pjet = Σ p of constituents Inclusive: independent of final states, just energy Very important: need pjet ≈ pparton Can use parton level calculation to predict jet properties
Tuesday, September 11, 12
13
PDF PDF Initial state radiation soft, long distance interactions Fragmentation (q,g ⇒ hadrons) ...
Tuesday, September 11, 12
B1 B2 J1 J2
14
Tuesday, September 11, 12
2 3 4 5 6 [pb] ! 10
2
10
3
10
4
10
5
10
6
10 2 3 4 5 6 [pb] ! 10
2
10
3
10
4
10
5
10
6
10 ATLAS Preliminary
L dt=2.43 pb
"
R=0.4, =7 TeV)+syst. s Data ( 1.11 × ALPGEN+HERWIG AUET1 0.65 × PYTHIA AMBT1 1.22 × ALPGEN+PYTHIA MC09’
Inclusive Jet Multiplicity 2 3 4 5 6 MC/Data 0.5 1 1.5 Inclusive Jet Multiplicity 2 3 4 5 6 MC/Data 0.5 1 1.5
100 200 300 400 500 600 700 800[pb/GeV]
T
/d p ! d
10 1 10
2
10
3
10
4
10
5
10
100 200 300 400 500 600 700 800[pb/GeV]
T
/d p ! d
10 1 10
2
10
3
10
4
10
5
10 ATLAS Preliminary
L dt=2.43 pb
"
R=0.4, =7 TeV)+syst. s Data ( 1.11 × ALPGEN+HERWIG AUET1 0.65 × PYTHIA AMBT1 1.22 × ALPGEN+PYTHIA MC09’ 2 #
jetsN
(leading jet) [GeV]
T
p 100 200 300 400 500 600 700 800 MC/Data 0.5 1 1.5 (leading jet) [GeV]
T
p 100 200 300 400 500 600 700 800 MC/Data 0.5 1 1.5
ATLAS-CONF-2011-043
ATLAS-CONF-2011-043, 7 TeV, 2.43 pb-1
15
Tuesday, September 11, 12
Tuesday, September 11, 12
Tuesday, September 11, 12
Tuesday, September 11, 12
jet jet jet jet
Tuesday, September 11, 12
Multiple interaction, underlying events, pile-up
jet jet jet jet “beam”
Tuesday, September 11, 12
Multiple interaction, underlying events, pile-up
jet jet jet jet “beam”
Tuesday, September 11, 12
ISR (beam) clustered
Multiple interaction, underlying events, pile-up
jet jet jet jet “beam”
Tuesday, September 11, 12
Part of the beam? ISR (beam) clustered
Multiple interaction, underlying events, pile-up
jet jet jet jet “beam”
Tuesday, September 11, 12
Proper choice of cone size? Part of the beam? ISR (beam) clustered
Multiple interaction, underlying events, pile-up
jet jet jet jet “beam”
Tuesday, September 11, 12
Proper choice of cone size? Part of the beam? ISR (beam) clustered
Multiple interaction, underlying events, pile-up
jet jet jet jet “beam”
Tuesday, September 11, 12
17
Tuesday, September 11, 12
θ pM pA pB
M = (pA + pB)2
Relevant kinematical variables
Tuesday, September 11, 12
1 2 n − 2 n − 1 M 1 2 n − 2 n − 1 M A B A B A B |Mn+1|2 = |M(p1, ...pA, pB)|2 |Mn(p1, ...pM)|2 |M(pM → pApB)|2 ×
collinear limit
|Mn+1|2dΠn+1 ' |Mn|2dΠn dt t αS 2π P(z)dzdφ
Tuesday, September 11, 12
1 2 n − 2 n − 1 M 1 2 n − 2 n − 1 M A B A B A B |Mn+1|2 = |M(p1, ...pA, pB)|2 |Mn(p1, ...pM)|2 |M(pM → pApB)|2 ×
collinear limit
|Mn+1|2dΠn+1 ' |Mn|2dΠn dt t αS 2π P(z)dzdφ
collinear singularity: t⇒0
Tuesday, September 11, 12
1 2 n − 2 n − 1 M 1 2 n − 2 n − 1 M A B A B A B |Mn+1|2 = |M(p1, ...pA, pB)|2 |Mn(p1, ...pM)|2 |M(pM → pApB)|2 ×
collinear limit
|Mn+1|2dΠn+1 ' |Mn|2dΠn dt t αS 2π P(z)dzdφ
collinear singularity: t⇒0 Splitting function IR singularity: z⇒0, 1
P(z) ∝ |M(pM → pApB)|2
Tuesday, September 11, 12
Pq→qg(z) = CF 1 + z2 1 − z , Pg→gg(z) = CA 1 − z z + z 1 − z + z(1 − z)
Pg→q¯
q(z) = TR
,
|Mn+1|2dΠn+1 ' |Mn|2dΠn dt t αS 2π P(z)dzdφ
Tuesday, September 11, 12
Prefers collinear radiation
P ~ (z)-1 prefers soft radiation
21
QM2 = t
Tuesday, September 11, 12
jets.
2 4 6 1 2 3 4 5 6 0 5 10 15 20 25 SISCone, R=1, f=0.75 y [GeV]
t
p φ
22
Tuesday, September 11, 12
23
soft radiation collinear splitting
Tuesday, September 11, 12
Soft or collinear radiation should not be able to induce “large” changes in the observable. Otherwise we cannot compute and compare with experiments.
23
soft radiation collinear splitting
Tuesday, September 11, 12
Soft or collinear radiation should not be able to induce “large” changes in the observable. Otherwise we cannot compute and compare with experiments.
23
soft radiation collinear splitting
Tuesday, September 11, 12
24
Tuesday, September 11, 12
100 200 300 400 500 −1 1
pT (GeV/c)
stable cones from seeds 100 200 300 400 500 −1 1
pT (GeV/c)
add soft particle
100 200 300 400 500 −1 1
pT (GeV/c)
resolve overlaps
25 No other radiation with the radius of cones centered on the seeds “stable cone”, clustering stops. 2 jets.
an event with 2 jets becomes an event with one jet because of a soft radiation
Tuesday, September 11, 12
26
Tuesday, September 11, 12
27
Tuesday, September 11, 12
28
Tuesday, September 11, 12
29
Tuesday, September 11, 12
30
Tuesday, September 11, 12
Done!
31
Tuesday, September 11, 12
η = − ln ⇤ cot θ 2 ⇥⌅
32
Tuesday, September 11, 12
dij = min(p2
T i, p2 T j)
∆R R0 ⇥2 , diB = p2
T i
dij = min(p−2
T i , p−2 T j )
∆R R0 ⇥2 , diB = p−2
T i
dij = ∆R R0 ⇥2 , diB = 1
A B
anti−kT
B A
C/A
B A
kT
33
Tuesday, September 11, 12
2 4 6 1 2 3 4 5 6 0 5 10 15 20 25
, R=1
t
k
y [GeV]
t
p φ
2 4 6 1 2 3 4 5 6 0 5 10 15 20 25
Cam/Aachen, R=1
y [GeV]
t
p φ
2 4 6 1 2 3 4 5 6 0 5 10 15 20 25 SISCone, R=1, f=0.75 y [GeV]
t
p φ
2 4 6 1 2 3 4 5 6 0 5 10 15 20 25
, R=1
t
anti-k
y [GeV]
t
p φ
34
Tuesday, September 11, 12
! Example: here’s an event with 500 GeV dijets (left), and
the same event with fifty pileup events (right).
! We’ll encounter this level of pileup next year, ! Somehow we’re going to have to find new physics in
this mess!
35
Some “clean up” procedure, filtering, pruning, trimming.
Tuesday, September 11, 12
!"pt#2
pert + !"pt#2 h + !"pt#2 UE [GeV2]
R Tevatron quark jets pt = 50 GeV 1 2 3 4 5 6 7 8 9 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 !"pt#2
pert
!"pt#2
h
!"pt#2
UE
0.4 0.5 0.6 0.7 0.8 0.9 1 50 500 100 1000 best R pt [GeV] Tevatron, gluon jets Tevatron, quark jets LHC, gluon jets LHC, quark jets
“best” R
Tuesday, September 11, 12
!"pt#2
pert + !"pt#2 h + !"pt#2 UE [GeV2]
R Tevatron quark jets pt = 50 GeV 1 2 3 4 5 6 7 8 9 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 !"pt#2
pert
!"pt#2
h
!"pt#2
UE
0.4 0.5 0.6 0.7 0.8 0.9 1 50 500 100 1000 best R pt [GeV] Tevatron, gluon jets Tevatron, quark jets LHC, gluon jets LHC, quark jets
“best” R
radiation out of the cone
Tuesday, September 11, 12
!"pt#2
pert + !"pt#2 h + !"pt#2 UE [GeV2]
R Tevatron quark jets pt = 50 GeV 1 2 3 4 5 6 7 8 9 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 !"pt#2
pert
!"pt#2
h
!"pt#2
UE
0.4 0.5 0.6 0.7 0.8 0.9 1 50 500 100 1000 best R pt [GeV] Tevatron, gluon jets Tevatron, quark jets LHC, gluon jets LHC, quark jets
“best” R
radiation out of the cone Universal noise ∝ R2
Tuesday, September 11, 12
Mass [GeV]
900 920 940 960 980 1000 1020 1040 1060 1080 1100
Cross Section [A.U.]
0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4
R=0.9 R=1.1 R=1.3 R=1.5
Mass [GeV]
900 920 940 960 980 1000 1020 1040 1060 1080 1100
Cross Section [A.U.]
0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16
R=0.9 R=1.1 R=1.3 R=1.5
FSR only Including ISR, MI, UE, pile-up
Tuesday, September 11, 12
pT i < fcut · Λhard
ISR argument.
38
Tuesday, September 11, 12
39
Tuesday, September 11, 12
! "
0.5 1 1.5 # "
0.5 1 1.5 ! "
0.5 1 1.5 # "
0.5 1 1.5 ! "
0.5 1 1.5 # "
0.5 1 1.5 ! "
0.5 1 1.5 # "
0.5 1 1.5
Start Cluster into subjets Discard soft subjets Reassemble 1 2 3 4
40
Tuesday, September 11, 12
Improvement fcut, Ncut Rsub R0, ρ Γ [GeV] M [GeV] anti-kT
71 522 anti-kT (N) 40% 5∗ 0.2∗ 1.5∗ 62 499 anti-kT (f, pT ) 59% 3 × 10−2∗ 0.2 1.5 52 475 anti-kT (f, H) 61% 1 × 10−2∗ 0.2 1.5 50 478 VR 30%
62 511 VR (N) 53% 5 0.2 275∗ GeV 53 498 VR (f, pT ) 68% 3 × 10−2 0.2 300∗ GeV 49 475 VR (f, H) 73% 1 × 10−2 0.2 300∗ GeV 47 478 Filtering 27% 2 R0/2 1.3∗ 61 515
Mass [GeV]
400 420 440 460 480 500 520 540 560 580 600
Cross Section [A.U.]
0.05 0.1 0.15 0.2 0.25
Tanti-k trimmed
Tanti-k
Mass [GeV]
400 420 440 460 480 500 520 540 560 580 600
Cross Section [A.U.]
0.05 0.1 0.15 0.2 0.25
VR VR trimmed
http://jthaler.net/jets/VR_Jets.html http://jthaler.net/jets/Jet_Trimming.html
41
Tuesday, September 11, 12
More faithful (smaller) jet mass for the background.
Jet Mass [GeV]
20 40 60 80 100 120 140 160 180 200
Cross Section [A.U.]
100 200 300 400 500
FSR only
Tanti-k ISR/MI/pileup
Tanti-k Trimming FSR only Trimming ISR/MI/pileup
Without contamination With “trimming” With contamination
Tuesday, September 11, 12
43
Tuesday, September 11, 12
X ¯ t t
e.g. Lower combinatorics, SM background boost differently. W, Z b b
For example, boost tops Brooijmans; Lillie, Randall, LTW; Thaler, LTW;
Virzi
... Butterworth, Davidson, Bubin, Salam
For a summary of recent developments: C. Vermilion,1001.1335 e.g. h
Tuesday, September 11, 12
Prefers collinear radiation
P ~ (z)-1 prefers soft radiation
45
Tuesday, September 11, 12
W + b t q E1 E2
Zooming in near the first splitting Soft radiation: Top. First splitting Jet mass: QCD. Decay: Jet mass:
46
Tuesday, September 11, 12
47
Tuesday, September 11, 12
∼
48
Tuesday, September 11, 12
Thaler and LTW, arXiv:0806.0023. Almeida, Lee, Perez, Sterman, Sung, Virzi, arXiv:0807.0234
Ikl
w =
⇤
i
wi pi,k wi pi,l wi
Pf = 4λ1λ2 (λ1 + λ2)2
49
Tuesday, September 11, 12
Planar Flow
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Cross Section [A.U.]
100 200 300 400 500
FSR only
Tanti-k ISR/MI/pileup
Tanti-k Trimming FSR only Trimming ISR/MI/pileup
With “trimming” With no contamination With contamination
50
Tuesday, September 11, 12
CMS Coll. CMS PAS JME-09-001
Virzi, arXiv:0807.0243
Barger, Huang, 1102.3183
efficiency 0.1 0.2 0.3 0.4 0.5 0.6 0.7 mistag rate
10
10
10 1
Hopkins CMS Pruning ATLAS Thaler/Wang
Boost 2010 proceeding, 1012.5421
51
Tuesday, September 11, 12
Mass (GeV)
20 40 60 80 100 120 140 160 180 200
Events / 8GeV / 30fb
20 40 60 80 100 120 140
Mass (GeV)
20 40 60 80 100 120 140 160 180 200
Events / 8GeV / 30fb
20 40 60 80 100 120 140
q q V+jets VV V+Higgs
= 4.5 B S/ in 112-128GeV
(d)
Butterworth, Davison, Rubin, Salam, 0802.2470
]
2
Higgs mass [GeV/c 20 40 60 80 100 120 140 160 180 200
Events / 8GeV / 30fb 10 20 30 40 50 ]
2
Higgs mass [GeV/c 20 40 60 80 100 120 140 160 180 200
Events / 8GeV / 30fb 10 20 30 40 50
tt V+jets VV Higgs
Total S = 16.3 B = 104.2 Range 104-136GeV
ATLAS preliminary
(simulation) (c)
ATL-PHYS-PUB-2009-088
52
Tuesday, September 11, 12
TN = ⇤
k
| pkT | min
⇥ ≡ T a
N + T b N + T 1 N + · · · + T N N Jet 2 Jet b Jet a Soft Jet 3 Jet 1
b a 1 3 2
p p − +
(c) pp → leptons plus jets.
Stewart, Tackmann, Waalewijn, 1004.2489 N-jet like event
Tuesday, September 11, 12
pp->H->bb pp->g->bb
Gallicchio, Schwartz, 1001.5027 Relative enhanced radiation consentration e.g., ttbar at Dzero, Haas Boston Jet Workshop
Tuesday, September 11, 12
F
e x a m p l e : B . B e l l a z z i n i , C . C s a k i , A . F a l k
s k i , A . W e i l e r , a r X i v : 9 1 . 3 2 1 , a r X i v : 9 6 . 3 2 6 For example: P . Graham, A. Pierce, J. Wacker, hep-ph/0605162
h
a a
Soft gluon jets, considered impossible.
55
Tuesday, September 11, 12
F
e x a m p l e : B . B e l l a z z i n i , C . C s a k i , A . F a l k
s k i , A . W e i l e r , a r X i v : 9 1 . 3 2 1 , a r X i v : 9 6 . 3 2 6 For example: P . Graham, A. Pierce, J. Wacker, hep-ph/0605162
Two “equal ”clusters
h
Less radiation
Boosting the Higgs.
56
Tuesday, September 11, 12
0.5 1 1.5 2 2.5 60 70 80 90 100 110 120 130 140
Cross Section [fb/10-GeV] Mass [GeV]
Signal Background 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 60 70 80 90 100 110 120 130 140
Cross Section [fb/10-GeV] Mass [GeV]
Signal Background
W/Z+h Chen, Nojiri, Sreethawong, 1006.1151 Falkowski, Krohn, Shelton, Thalapillil, and LTW, 1006.1650 ttbar+h
57
Tuesday, September 11, 12
Kribs, Martin, Roy, Spannowsky , 0912.4731, 1006.1656 Kribs, Martin, and Roy, 1012.2886 Cui, Han, Schwartz, 1012.2077 Katz, Son, Tweedie, 1010.5253 Butterworth, Ellis, Raklev, Salam, 0906.0728 Fan, Krohn, Mosteiro, Thalapillil, 1102.0302
Tuesday, September 11, 12