ATLAS Heavy Flavour production Looking towards Run 2
- Heavy Flavour at the LHC
IPPP Durham 21/4/2016
ATLAS Heavy Flavour production Looking towards Run 2 Heavy - - PowerPoint PPT Presentation
ATLAS Heavy Flavour production Looking towards Run 2 Heavy Flavour at the LHC IPPP Durham 21/4/2016 Josh McFayden ATLAS Heavy Flavour production Looking towards Run 2 into Heavy Flavour at the LHC IPPP
IPPP Durham 21/4/2016
IPPP Durham 21/4/2016
Josh McFayden | HF @ LHC | 21/4/2016
Chiara already showed some of the Run 1 results and
3
Josh McFayden | MC performance | 21/05/2015 Josh McFayden | HF @ LHC | 21/4/2016
4
Josh McFayden | HF @ LHC | 21/4/2016 5
New 7 TeV ATLAS result on di-b-jet production.
One jet with pT > 270 GeV required due to trigger. b-jets with pT > 20 GeV and ∆R = 0.4. Template fit used to extract true b-b contribution.
mbb step at ~500 GeV due to “turn-on” of flavour creation. STDM-2013-03
Josh McFayden | HF @ LHC | 21/4/2016 6
New 7 TeV ATLAS result on di-b-jet production.
One jet with pT > 270 GeV required due to trigger. b-jets with pT > 20 GeV and ∆R = 0.4. Template fit used to extract true b-b contribution.
Large ΔΦ region is dominated by flavour creation and
STDM-2013-03
Josh McFayden | HF @ LHC | 21/4/2016
0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
[pb] R(b,b) ∆ d (Zbb) σ d 0.05 0.1 0.15 0.2 0.25 0.3
a Data
= 7 TeV, 4.6 fb s MCFM aMC@NLO 5FNS aMC@NLO 4FNS ALPGEN+HJ SHERPA
ATLAS 2 b-jets ≥ Z+
Data NLO
0.5 1 R(b,b) ∆
0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
Data LO multileg
0.4 0.6 0.8 1 1.2
7
These results use the same dataset…
How much do we learn about V+HF from inclusive di-b-jets?
Seems like trends might be different?
Is the large leading jet requirement good/bad?
Josh McFayden | MC performance | 21/05/2015 Josh McFayden | HF @ LHC | 21/4/2016
8
Josh McFayden | HF @ LHC | 21/4/2016 9
Generator Matrix Element Flavour number scheme Slicing/ filtering Comments Sherpa2.1 NLO@2j LO@4j 5fl pT(V) & HF
Known mismodelling
Sherpa2.2 NLO@2j LO@4j 5fl pT(V) & HF
Improvement in low pT large η jets
MG+Py8 A LO@4j 5fl N-parton
Mismodelling of jet pT (too hard)
MG+Py8 B LO@4j 5fl HT & HF
NLO PDF, difgerent shower settings. (still too hard)
MG5_aMC+Py8 FxFx NLO@2j 5fl TBD
Very promising - some N-jets mismodelling
Alpgen LO@5j 4fl N-parton & b/c/light
The new old! HFOR can be problematic.
Josh McFayden | HF @ LHC | 21/4/2016
What appears to be ~small slope at 7 TeV seems to
Important to make new measurements at 13 TeV.
10
ATLAS data, √s=7 TeV Sherpa 2.1 Sherpa 2.2 MG+Py8 A MG+Py8 B aMC@NLO FxFx 10−4 10−3 10−2 10−1 1 Transverse momentum of 1st jet dσ/dp⊥ [pb/GeV] 100 200 300 400 500 600 700 0.6 0.8 1 1.2 1.4 p⊥(1st jet) [GeV] MC/Data Sherpa 2.1 MG+Py8 A MG+Py8 B ATLAS simulation preliminary, √s = 13 TeV 10−3 10−2 10−1 1 10 1 Transverse momentum of 1st jet dσ/dp⊥ [pb/GeV] 100 200 300 400 500 600 700 0.6 0.8 1 1.2 1.4 p⊥(1st jet) [GeV] Ratio to Sherpa 2.1
7 TeV - ratio wrt data 13 TeV - ratio wrt Sherpa2.1 ATL-PHYS-PUB-2016-003
Josh McFayden | HF @ LHC | 21/4/2016 11
)
2
,jet
1
(jet φ ∆ 0.5 1 1.5 2 2.5 3 Data / Pred 0.9 1 1.1
0.5 1 1.5 2 2.5 3
rad
x10 π Events / 5000 10000 15000 20000 25000
Data 2012 =1.0) µ VH(bb) ( Diboson t t Single top Multijet W+hf W+cl W+l Z+hf Z+cl Z+l
ATLAS = 8 TeV s
L dt = 20.3 fb
[GeV]
W T
p 50 100 150 200 250 Data / Pred 0.9 1 1.1
50 100 150 200 250
Events / 5 GeV 10000 20000 30000 40000 50000 60000 70000 80000
Data 2012 =1.0) µ VH(bb) ( Diboson t t Single top Multijet W+hf W+cl W+l Z+hf Z+cl Z+l
ATLAS = 8 TeV s
L dt = 20.3 fb
7 TeV VH(→bb) analysis saw large mismodelling of Δφ(j,j)
Difference between Sherpa and Alpgen
Are in a better position with our Run 2 generators setups?
JHEP01(2015)069
Josh McFayden | HF @ LHC | 21/4/2016
Some systematic difference between Sherpa and
The data sits between the two. The modelling seems to be improved but much less data.
12
ATLAS data, √s=7 TeV Sherpa 2.1 Sherpa 2.2 MG5 aMC+Py8 CKKW-L A MG5 aMC+Py8 CKKW-L B MG5 aMC+Py8 FxFx 1 10 1 dσ/d|∆φ| [pb] 0.5 1 1.5 2 2.5 3 0.6 0.8 1 1.2 1.4 |∆φ|(1st jet, 2nd jet) MC/Data Sherpa 2.1 MG+Py8 A MG+Py8 B ATLAS simulation preliminary, √s = 13 TeV 10 1 Azimuthal distance of leading jets dσ/d|∆φ| [pb] 0.5 1 1.5 2 2.5 3 0.6 0.8 1 1.2 1.4 |∆φ|(1st jet, 2nd jet) Ratio to Sherpa 2.1
7 TeV - ratio wrt data 13 TeV - ratio wrt Sherpa2.1 ATL-PHYS-PUB-2016-003
Josh McFayden | HF @ LHC | 21/4/2016
Some systematic difference between Sherpa and
The data seems to prefer the MG5_aMC+Py8 CKKW-L shape.
13
ATLAS data, √s=7 TeV Sherpa 2.1 Sherpa 2.2 MG5 aMC+Py8 CKKW-L A MG5 aMC+Py8 CKKW-L B MG5 aMC+Py8 FxFx 10−3 10−2 10−1 dσ/dm [pb/GeV] 200 400 600 800 1000 0.6 0.8 1 1.2 1.4 m(1st jet, 2nd jet) [GeV] MC/Data Sherpa 2.1 MG+Py8 A MG+Py8 B ATLAS simulation preliminary, √s = 13 TeV 10−2 10−1 Invariant mass of leading jets dσ/dm [pb/GeV] 200 400 600 800 1000 0.6 0.8 1 1.2 1.4 m(1st jet, 2nd jet) [GeV] Ratio to Sherpa 2.1
7 TeV - ratio wrt data 13 TeV - ratio wrt Sherpa2.1 ATL-PHYS-PUB-2016-003
Josh McFayden | HF @ LHC | 21/4/2016
Main difference between Sherpa2.1 and MG5_aMC CKKW-L is in
From early analyses we see that the data seems to prefer the higher rate. Shape deviations are more important.
14
ATLAS data, √s=7 TeV Sherpa 2.1 Sherpa 2.2 MG5 aMC+Py8 CKKW-L A MG5 aMC+Py8 CKKW-L B MG5 aMC+Py8 FxFx 10−4 10−3 10−2 10−1 dσ(Zb)/dpT /Nb-jets [pb/GeV] 10 2 0.6 0.8 1 1.2 1.4 b-jet pT MC/Data Sherpa 2.1 MG+Py8 A MG+Py8 B ATLAS simulation preliminary, √s = 13 TeV 10−3 10−2 10−1 1 Z+ ≥ 1 b-jet dσ(Zb)/dpT/Nb-jets [pb/GeV] 10 2 0.6 0.8 1 1.2 1.4 b-jet pT Ratio to Sherpa 2.1
7 TeV - ratio wrt data 13 TeV - ratio wrt Sherpa2.1 ATL-PHYS-PUB-2016-003
Josh McFayden | HF @ LHC | 21/4/2016
Perhaps larger differences observed at 7 TeV than 13 TeV.
We have enough data to constrain the MC prediction here. Improvement by going to NLO.
15
ATLAS data, √s=7 TeV Sherpa 2.1 Sherpa 2.2 MG5 aMC+Py8 CKKW-L A MG5 aMC+Py8 CKKW-L B MG5 aMC+Py8 FxFx 1 10 1 dσ(Zb)/d∆φ(Z, b)/Nb-jets [pb] 0.5 1 1.5 2 2.5 3 0.6 0.8 1 1.2 1.4 ∆φ(Z, b) MC/Data Sherpa 2.1 MG+Py8 A MG+Py8 B ATLAS simulation preliminary, √s = 13 TeV 1 10 1 Z+ ≥ 1 b-jet, pT(Z) > 20 GeV dσ(Zb)/d∆φ(Z, b)/Nb-jets [pb] 0.5 1 1.5 2 2.5 3 0.6 0.8 1 1.2 1.4 ∆φ(Z, b) Ratio to Sherpa 2.1
7 TeV - ratio wrt data 13 TeV - ratio wrt Sherpa2.1 ATL-PHYS-PUB-2016-003
Josh McFayden | HF @ LHC | 21/4/2016
The familiar plot… Systematic shape differences observed
Sherpa seems to do a better job of modelling the shape of the
16
ATLAS data, √s=7 TeV Sherpa 2.1 Sherpa 2.2 MG5 aMC+Py8 CKKW-L A MG5 aMC+Py8 CKKW-L B MG5 aMC+Py8 FxFx 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 dσ(Zbb)/d∆R(b, b) [pb] 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 0.6 0.8 1 1.2 1.4 ∆R(b, b) MC/Data Sherpa 2.1 MG+Py8 A MG+Py8 B ATLAS simulation preliminary, √s = 13 TeV 0.2 0.4 0.6 0.8 1 Z+ ≥ 2 b-jets dσ(Zbb)/d∆R(b, b) [pb] 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 0.6 0.8 1 1.2 1.4 ∆R(b, b) Ratio to Sherpa 2.1
7 TeV - ratio wrt data 13 TeV - ratio wrt Sherpa2.1 ATL-PHYS-PUB-2016-003
Josh McFayden | HF @ LHC | 21/4/2016
Sherpa ~flat in m(b,b), some shape deviation from other
17
ATLAS data, √s=7 TeV Sherpa 2.1 Sherpa 2.2 MG5 aMC+Py8 CKKW-L A MG5 aMC+Py8 CKKW-L B MG5 aMC+Py8 FxFx 10−3 10−2 dσ(Zbb)/dm(b, b) [pb/GeV] 50 100 150 200 250 300 350 0.6 0.8 1 1.2 1.4 m(b, b) [GeV] MC/Data Sherpa 2.1 MG+Py8 A MG+Py8 B ATLAS simulation preliminary, √s = 13 TeV 10−3 10−2 Z+ ≥ 2 b-jets dσ(Zbb)/dm(b, b) [pb/GeV] 50 100 150 200 250 300 350 0.6 0.8 1 1.2 1.4 m(b, b) [GeV] Ratio to Sherpa 2.1
7 TeV - ratio wrt data 13 TeV - ratio wrt Sherpa2.1 ATL-PHYS-PUB-2016-003
Josh McFayden | MC performance | 21/05/2015 Josh McFayden | HF @ LHC | 21/4/2016
18
Josh McFayden | HF @ LHC | 21/4/2016
Average CPU time/event Sherpa has the highest CPU time/evt.
19
Josh McFayden | HF @ LHC | 21/4/2016
Average CPU time/event Sherpa has the highest CPU time/evt. Longer than full
20
Josh McFayden | HF @ LHC | 21/4/2016
Average CPU time/event Sherpa has the highest CPU time/evt. Longer than full
Dominated by high
NLO significantly
B- and C-hadron
21
Josh McFayden | HF @ LHC | 21/4/2016
Average CPU time/event Sherpa has the highest CPU time/evt. Longer than full
Dominated by high
NLO significantly
B- and C-hadron
22
Josh McFayden | HF @ LHC | 21/4/2016 23
Generator Matrix Element Flavour number scheme Slicing/ filtering Comments Sherpa2.1 NLO@2j LO@4j 5fl pT(V) & HF
Known mismodelling
Sherpa2.2 NLO@2j LO@4j 5fl pT(V) & HF
Improvement in low pT large η jets
MG+Py8 A LO@4j 5fl N-parton
Mismodelling of jet pT (too hard)
MG+Py8 B LO@4j 5fl HT & HF
NLO PDF, difgerent shower settings. (still too hard)
MG5_aMC+Py8 FxFx NLO@2j 5fl TBD
Very promising - some N-jets mismodelling
Alpgen LO@5j 4fl N-parton & b/c/light
The new old! HFOR can be problematic.
Josh McFayden | HF @ LHC | 21/4/2016
Generator Matrix Element Flavour number scheme Slicing/ filtering Comments Sherpa2.1 NLO@2j LO@4j 5fl pT(V) & HF
Known mismodelling
Sherpa2.2 NLO@2j LO@4j 5fl pT(V) & HF
Improvement in low pT large η jets
MG+Py8 A LO@4j 5fl N-parton
Mismodelling of jet pT (too hard)
MG+Py8 B LO@4j 5fl HT & HF
NLO PDF, difgerent shower settings. (still too hard)
MG5_aMC+Py8 FxFx NLO@2j 5fl TBD
Very promising - some N-jets mismodelling
Alpgen LO@5j 4fl N-parton & b/c/light
The new old! HFOR can be problematic.
24
Josh McFayden | MC Perf | 12/4/2016
Modelling
Significant V+jets modelling improvements since Run 1. Quite large HF variations but also large uncertainty on the data.
Conspicuous absence of 4-flavour samples due to difficulty in HFOR.
Can we do better here?
Is the variation between several 5-fl scheme MCs sufficient for
Should we always add comparison to 4-fl scheme predictions? Usually we take normalisation from a control region so only shape/
Are there any other uncertainties that are not taken into account with
25
Josh McFayden | MC performance | 21/05/2015 Josh McFayden | HF @ LHC | 21/4/2016
26
Josh McFayden | HF @ LHC | 21/4/2016
Search for VH production
Interpreted HVT model.
Selection very similar to
27
(GeV)
Z’
m 1000 1500 2000 2500 3000 ) [pb] c , c b b
× ZH)
×
1
1 10 Preliminary ATLAS
Analysis
b b
ll+J
L dt = 3.2 fb
s
Observed Expected 1 s.d. ± Expected 2 s.d. ± Expected =1
V
HVT Model A, g =3
V
HVT Model B, g
ATLAS-CONF-2015-074
Josh McFayden | HF @ LHC | 21/4/2016
Search for VH production
Interpreted HVT model.
Selection very similar to
28
(GeV)
Z’
m 1000 1500 2000 2500 3000 ) [pb] c , c b b
× ZH)
×
1
1 10 Preliminary ATLAS
Analysis
b b
ll+J
L dt = 3.2 fb
s
Observed Expected 1 s.d. ± Expected 2 s.d. ± Expected =1
V
HVT Model A, g =3
V
HVT Model B, g
High mass optimisation means a
ttbar and V+HF are dominant
ATLAS-CONF-2015-074
Josh McFayden | HF @ LHC | 21/4/2016
Search for VH production
Interpreted HVT model.
Selection very similar to
29
(GeV)
Z’
m 1000 1500 2000 2500 3000 ) [pb] c , c b b
× ZH)
×
1
1 10 Preliminary ATLAS
Analysis
b b
ll+J
L dt = 3.2 fb
s
Observed Expected 1 s.d. ± Expected 2 s.d. ± Expected =1
V
HVT Model A, g =3
V
HVT Model B, g
High mass optimisation means a
ttbar and V+HF are dominant
Fit m(T)VH distribution to look for
ATLAS-CONF-2015-074
Josh McFayden | HF @ LHC | 21/4/2016
Search for VH production
Interpreted HVT model.
Selection very similar to
30
(GeV)
Z’
m 1000 1500 2000 2500 3000 ) [pb] c , c b b
× ZH)
×
1
1 10 Preliminary ATLAS
Analysis
b b
ll+J
L dt = 3.2 fb
s
Observed Expected 1 s.d. ± Expected 2 s.d. ± Expected =1
V
HVT Model A, g =3
V
HVT Model B, g
Very low stats in 2-btag high
ATLAS-CONF-2015-074
Josh McFayden | HF @ LHC | 21/4/2016
Signal region defined
Control regions defined
Low mjet region used for
Significant
Uncertainty on V+HF
Mismodelling observed
31
Josh McFayden | MC performance | 21/05/2015 Josh McFayden | HF @ LHC | 21/4/2016
32
Josh McFayden | HF @ LHC | 21/4/2016 33
Measurements in the pipeline
8 TeV Jpsi+mu 13 TeV W/Z+HF
Flavour-labelled kT-spittings?
Josh McFayden | HF @ LHC | 21/4/2016
34
Measurements in the pipeline
8 TeV Jpsi+mu 13 TeV W/Z+HF
Flavour-labelled kT-spittings?
Josh McFayden | HF @ LHC | 21/4/2016 35
Measurements in the pipeline
8 TeV Jpsi+mu 13 TeV W/Z+HF
Flavour-labelled kT-spittings?
Josh McFayden | HF @ LHC | 21/4/2016 36
Repeat of 7 TeV measurements
Add W+HF to Z+HF measurement More statistics
dR(b,b) with smaller statistical uncertainty
More differential measurements for W+HF Add dedicated measurement in boosted regime So far plan to have the following MC comparisons
Sherpa 2.2 5fl MG5_aMC+Py8 CKKW-L 5fl MG5_aMC+Py8 FxFx 5fl Alpgen+Py6 4fl
What should we add?
MG5_aMC NLO Z+bb 4-fl Sherpa 4-fl?
Josh McFayden | HF @ LHC | 21/4/2016 37
A new measurement of B hadron pair production in
Identifying the B → J/ψ(→ μ+μ−) + X and B → μ + X decay Fit to extract differential cross
Non-prompt J/ψ Non-prompt 3rd μ
New constraints particularly in the
Sensitive to g→bb production. No jet radius restriction.
Josh McFayden | HF @ LHC | 21/4/2016 38
Signal
mass(μμ) [MeV] pseudo-proper lifetime [ps] d0 significance BDT output Extract non-prompt
Repeat in each
kinematic bin
2-D fit to J/ψ mass
Separate from
prompt J/ψ and fakes
2-D fit to 3rd μ
Separate from
prompt, decay-in- flight and pileup μ.
Josh McFayden | HF @ LHC | 21/4/2016 39
Signal cross section extracted as function of:
ΔR(J/ψ,μ) Δφ(J/ψ,μ) Δy(J/ψ,μ) pT(J/ψ,μ) m(J/ψ,μ) yboost(J/ψ,μ) m/pT pt/m
~20% xs
Josh McFayden | MC performance | 21/05/2015 Josh McFayden | HF @ LHC | 21/4/2016
40
Josh McFayden | HF @ LHC | 21/4/2016 41
Technical obstacles
HFOR HF-enhanced generation to increase efficiency?
Z+HF W+HF
More stats, more kinematic variables…
BB hadron J/ψ+μ 0.2 kT track jets?
Josh McFayden | MC performance | 21/05/2015 Josh McFayden | HF @ LHC | 21/4/2016
42
Josh McFayden | HF @ LHC | 21/4/2016
Signal region defined
Control regions defined
Low mjet region used
Significant
Uncertainty on V+HF
43
00 200 300 400 00 Events / 25 GeV 100 200 300 400 500 600 700 800
Data 2015 HVT(M=2.0 TeV) Diboson t t Single top VH W+b W+c W+l Z+b Z+c Z+l Uncertainty Pre-fit background
ATLAS Preliminary
Ldt = 3.2 fb
s 1 lep., 1 tag, 0 add. tags
[GeV]
jet
m 100 200 300 400 500 Data/Pred 0.8 1 1.2
Josh McFayden | HF @ LHC | 21/4/2016
Signal region defined
Control regions defined
Low mjet region used for
Significant
Low statistics in tail
44
00 200 300 400 00 Events / 25 GeV 100 200 300 400 500 600 700 800
Data 2015 HVT(M=2.0 TeV) Diboson t t Single top VH W+b W+c W+l Z+b Z+c Z+l Uncertainty Pre-fit background
ATLAS Preliminary
Ldt = 3.2 fb
s 1 lep., 1 tag, 0 add. tags
[GeV]
jet
m 100 200 300 400 500 Data/Pred 0.8 1 1.2
600 200 800 2400 3000 Events / 100 GeV 50 100 150 200 250
Data 2015 HVT(M=2.0 TeV) Diboson t t Single top VH W+b W+c W+l Z+b Z+c Z+l Uncertainty Pre-fit background
ATLAS Preliminary
Ldt = 3.2 fb
s 1 lep., 1 tag, 0 add. tags > 145 GeV
jet
< 65 GeV || m
jet
m
[GeV]
VH
m 500 1000 1500 2000 2500 3000 Data/Pred 0.8 1 1.2 600 200 800 2400 3000 Events / 100 GeV
1
1 10
2
10
3
10
4
10
Data 2015 HVT(M=2.0 TeV) Diboson t t Single top VH W+b W+c W+l Z+b Z+l Uncertainty Pre-fit background
ATLAS Preliminary
Ldt = 3.2 fb
s 1 lep., 1 tag, 0 add. tags < 145 GeV
jet
75 GeV < m
[GeV]
VH
m 500 1000 1500 2000 2500 3000 Data/Pred 0.5 1 1.5 2
Josh McFayden | HF @ LHC | 21/4/2016
Signal region defined by fat jet mass
Control regions defined
Low mjet region used for
Significant contribution
Very low statistics after
Will be worse with
45
00 200 300 400 00 Events / 25 GeV 100 200 300 400 500 600 700 800
Data 2015 HVT(M=2.0 TeV) Diboson t t Single top VH W+b W+c W+l Z+b Z+c Z+l Uncertainty Pre-fit background
ATLAS Preliminary
Ldt = 3.2 fb
s 1 lep., 1 tag, 0 add. tags
[GeV]
jet
m 100 200 300 400 500 Data/Pred 0.8 1 1.2
600 200 800 2400 3000 Events / 100 GeV 50 100 150 200 250
Data 2015 HVT(M=2.0 TeV) Diboson t t Single top VH W+b W+c W+l Z+b Z+c Z+l Uncertainty Pre-fit background
ATLAS Preliminary
Ldt = 3.2 fb
s 1 lep., 1 tag, 0 add. tags > 145 GeV
jet
< 65 GeV || m
jet
m
[GeV]
VH
m 500 1000 1500 2000 2500 3000 Data/Pred 0.8 1 1.2 600 200 800 2400 3000 Events / 100 GeV
1
1 10
2
10
3
10
Data 2015 HVT(M=2.0 TeV) Diboson t t Single top VH W+b W+c W+l Z+b Z+c Z+l Uncertainty Pre-fit background
ATLAS Preliminary
Ldt = 3.2 fb
s 1 lep., 2 tags, 0 add. tags < 145 GeV
jet
75 GeV < m
[GeV]
VH
m 500 1000 1500 2000 2500 3000 Data/Pred 0.5 1 1.5 2
Josh McFayden | HF @ LHC | 21/4/2016 46
ATLAS-CONF-2015-041
Josh McFayden | HF @ LHC | 21/4/2016 47
B-veto C-veto C-Filter B-veto B-Filter PtV0-70 PtV70-140 PtV140-280 PtV280-500 PtV500-700 PtV700-1000 PtV1000-2000 PtV2000+
Josh McFayden | HF @ LHC | 21/4/2016 48
χ
ATLAS
∆
Josh McFayden | HF @ LHC | 21/4/2016
LO->NLO
Sample according to the optimised integrators, but for each PS
Also the unweighting efficiency is much worse for the
High pTV
The integrators do not adapt very well to the extreme phase
Also, in these regions more multi-parton ME calculations are
49