projects in heidelbergplehn/includes/talks/2010/mc_ka_10.pdf · projects in heidelberg tilman plehn...

Projects...

Tilman Plehn

MadGOLEM

SHERPA

HEPTopTagger

SUSY

Higgs

Projects in Heidelberg

Tilman Plehn

Universitat Heidelberg

Terascale Meeting, Karlsruhe 12/2010

Projects...

Tilman Plehn

MadGOLEM

SHERPA

HEPTopTagger

SUSY

Higgs

MadGOLEM

Multileg ME generation modified MADGRAPH Feynman diagrams+ spinor-helicity formalism

One-loop ME generation QGRAF graph theoryReduction to scalar integrals GOLEM dim. reg. (UV & IR)

Passarino, Veltman [’79]IR subtraction modified MADGRAPH Catani, Seymour [’97]

own routinesUV renormalization modified MADGRAPH OS field renormalization

own routines MSbar αsSUSY-restoration

OS subtraction own routines PROSPINO scheme [’96]applied in MC@NLO [’08]

Projects...

Tilman Plehn

MadGOLEM

SHERPA

HEPTopTagger

SUSY

Higgs

MadGOLEM

u

u

e˜

e˜γ

u

u

e˜

e˜

Z

P1: pp → X1X2

leading order: MadGRAPH

Projects...

Tilman Plehn

MadGOLEM

SHERPA

HEPTopTagger

SUSY

Higgs

MadGOLEM

u

u

e ˜

e ˜γg

u1e˜

u1e˜

u

u

e ˜

e ˜γg

u

u

P2: pp → X1X2 + j

NLO - dipoles: MadGOLEM

Projects...

Tilman Plehn

MadGOLEM

SHERPA

HEPTopTagger

SUSY

Higgs

MadGOLEM

u

u

e ˜

e ˜γg

u1e˜

u1e˜

u

u

e ˜

e ˜γg

u

u

P2: pp → X1X2 + j

NLO - virtual corrections: QGRAF + MadGOLEM

Projects...

Tilman Plehn

MadGOLEM

SHERPA

HEPTopTagger

SUSY

Higgs

MadGOLEM

u

u

e ˜

e ˜

g

uγ

u

u

e ˜

e ˜

g

uZ

P3: pp → X1X2 + j

NLO - real emission: MadGRAPH

Projects...

Tilman Plehn

MadGOLEM

SHERPA

HEPTopTagger

SUSY

Higgs

MadGOLEM: tests

pp → eLe∗L

100 200 300 400M

squark [GeV]

47.9

48

48.1

48.2

48.3

48.4

σ [

fb]

σNLO - prospino

σNLO - MadGOLEM

MseL

= 150 GeV

Mgluino

= 200 GeV

0 200 400 600 800 1000µ

F = µ

R [GeV]

46

46.5

47

47.5

48

σ [

fb]

σNLO - prospino

σNLO - MadGOLEM

MseL

= 150 GeV

Mgluino

= 200 GeV

Msquark

= 3000 GeV

– rate as a function of the squark mass

– rate as a function of the renormalization scale

– red points from MadGOLEM, black line from Prospino:

Projects...

Tilman Plehn

MadGOLEM

SHERPA

HEPTopTagger

SUSY

Higgs

MadGOLEM: tests

pp → uχ e+e− → uu∗

1e-08

1e-06

0.0001

0.01

1

1e+02

1e+04|M|^2DipoleAbs(|M|^2 - Dipole)

Dipole Vs |M|^2, without pdf’sug -> uln1 g, m_ul = 150 GeV, m_n1 = 90 GeV

1Dipole/|M|^2

1e-12 1e-11 1e-10 1e-09 1e-08 1e-07 1e-06 1e-05 0.0001 0.001 0.01 0.1 1Eg^2 / Ecm^2

-0.04-0.020.000.020.04

|M|^2 - Dipole

0.0001 0.01 1alpha_fi

-5e-09

0

5e-09

1e-08 int dip_firealsum

ug>uln1susy dipole fi

– SUSY dipoles in pp → uχ for soft gluons

– α dependence of the SUSY dipoles in e+e− → uu∗

Projects...

Tilman Plehn

MadGOLEM

SHERPA

HEPTopTagger

SUSY

Higgs

Current Sherpa-1.2.3 (Dec 7th)

I POWHEG & MENLOPS [arXiv:1008.5399 & arXiv:1009.1127]I DY & W production, WW dibosonsI Higgs production in gluon fusion

I PROFESSOR tuningI hadron data from Tevatron & LHC only tuned to minbias/UE, not to DY, W, ...

I different PDF sets: CTEQ6L1 & CTEQ66 (default) LO∗ planned

I LHE output for MEsI various bugfixes

http://www.sherpa-mc.de/

Projects...

Tilman Plehn

MadGOLEM

SHERPA

HEPTopTagger

SUSY

Higgs

POWHEG vs. MEPS vs. MENLOPS

example: Drell-Yan at 14 TeV LHC

LO+PS × 1.2ME+PS (3-jet) × 1.2POWHEG

MENLOPS (3-jet)

1 10 1 10 210−6

10−5

10−4

10−3

10−2

10−1

1

10 1

pℓℓ⊥ [GeV]

dσ/dp

ℓℓ ⊥[pb/GeV

]

LO+PS × 1.2ME+PS (3-jet) × 1.2POWHEG

MENLOPS (3-jet)

10 210−9

10−8

10−7

10−6

10−5

10−4

10−3

10−2

10−1

p⊥(jet 3) [GeV]

dσ/dp⊥(

jet

3)[pb/GeV

]

MENLOPS:I inclusive cross section at NLOI multiple hard emissions through tree level MEs

Projects...

Tilman Plehn

MadGOLEM

SHERPA

HEPTopTagger

SUSY

Higgs

POWHEG vs. MEPS vs. MENLOPS

example: WW production at 14 TeV LHC

MENLOPS (3-jet)ME+PS (3-jet) × 1.04POWHEG

10−7

10−6

10−5

10−4

10−3

HT

dσ/dH

T[pb/GeV

]

10 2 10 3

0.6

0.8

1

1.2

1.4

HT [GeV]

Ratio

MENLOPS (3-jet)ME+PS (3-jet) × 1.04POWHEG

0

0.02

0.04

0.06

0.08

0.1

Azimuthal decorrelation of leading and second leading jet

dσ/d

∆φ(jet

1,jet

2)[pb]

0 0.5 1 1.5 2 2.5 3

0.6

0.8

1

1.2

1.4

∆φ(jet 1, jet 2)

Ratio

HT = ET/ +Xi=l,j

pT ,i

Projects...

Tilman Plehn

MadGOLEM

SHERPA

HEPTopTagger

SUSY

Higgs

UE at LHC 900 GeV – included in tuning

ATLAS-CONF-2010-081

b

b

b

bb

bb b b b b

b bb b

b b bb b b

b

bb

ATLAS datab

CTEQ6L1CTEQ66ATLAS preliminary

0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

Toward Nch density vs p⊥ (leading track),√

s = 900 GeV

d2〈N

ch〉/

dη

dφ

1 2 3 4 5 6 7 8 9 10

0.6

0.8

1

1.2

1.4

p⊥ (leading track) [GeV]

MC

/d

ata

bb

bb b b b b b b b b b

b b b b b bb

b

b

b

b

ATLAS datab

CTEQ6L1CTEQ66ATLAS preliminary

0

0.2

0.4

0.6

0.8

1

1.2

Transverse 〈p⊥〉 vs p⊥ (leading track),√

s = 900 GeV

〈p⊥〉

[GeV

]

1 2 3 4 5 6 7 8 9 10

0.6

0.8

1

1.2

1.4

p⊥ (leading track) [GeV]M

C/

da

ta

Projects...

Tilman Plehn

MadGOLEM

SHERPA

HEPTopTagger

SUSY

Higgs

UE at LHC 7 TeV – included in tuning

ATLAS-CONF-2010-081

[rad]wrt leading track

φ∆

­3 ­2 ­1 0 1 2 3

φ∆

dη

/dch

N2d

0.5

1

1.5

2

2.5

= 7 TeVs

> 1, 2, 3, 5 GeV, bottom to toplead

Tp

| <2.5,η> 0.5 GeV and |T

p

PreliminaryATLAS

Transverse TowardAway Transverse Away

Data 2010PYTHIA ATLAS MC09 Tune

b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b bb b b

b bb bb bb b b b b b b b b b b b b b b b b b b b b b b b b b b b b b

b b b b bb b b b

b b b b bb b b b b b b

b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b bb bb bbb bbbbb b b b

bbbb bb b b b b b b b b b b b b b b b b b b

b b b bb b b

b b bb b b

b b bb b b b

b b b

b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b bb bb bbbbbbbbbb bbbbbbbbb b b b b b b b b b b b b b b b

b b b bb b b

b bb b b

b bb bb bb bb b b

b b

bb bb bbbbbb bbb b b b

b bb bb b b b b b b b b b b b b b b b

bb bb bbb

b

b

b

b

b

b b

b

b

b

b

b

bbb bb bb b b b b b b b b b b b

b b bb bb bb bb bb bbb bbbbbb bb bb

ATLAS datab

Sherpa-1.2.3 cteq6l1ATLAS preliminary

-3 -2 -1 0 1 2 3

0.5

1

1.5

2

2.5Nch vs ∆φ wrt leading track,

√s = 7TeV

∆φ wrt leading trackd2Nch/d

ηd

∆φ

Projects...

Tilman Plehn

MadGOLEM

SHERPA

HEPTopTagger

SUSY

Higgs

HEPTopTagger

Boosted top quarks [Kaplan, Rehermann, Schwartz, Tweedie; Princeton, Seattle...]

– hadronic tops with pT & 800 GeVisolation and b tagging challenging

– C/A algorithm with pT drop criterionall top decay jets identified3 kinematic constraints: mW ,mt , cos θhel [no b tag]

– top mass included, no sidebins

HEPTopTagger [TP, Salam, Spannowsky, Takeuchi]

– extend to lower pT & 250 GeVrealistic for t t in Standard Model

[GeV]T

p200 300 400 500 600

Num

ber

of to

ps

210

310

410

510

610

QCD

top

tagged top

W+jets

Projects...

Tilman Plehn

MadGOLEM

SHERPA

HEPTopTagger

SUSY

Higgs

HEPTopTagger

Boosted top quarks [Kaplan, Rehermann, Schwartz, Tweedie; Princeton, Seattle...]

– hadronic tops with pT & 800 GeVisolation and b tagging challenging

– C/A algorithm with pT drop criterionall top decay jets identified3 kinematic constraints: mW ,mt , cos θhel [no b tag]

– top mass included, no sidebins

HEPTopTagger [TP, Salam, Spannowsky, Takeuchi]

– extend to lower pT & 250 GeVrealistic for t t in Standard Model

1

10

210

310

[GeV]TP0 200 400 600

bjj

R∆

0

1

2

3

1

10

210

[GeV]TP0 200 400 600

bjj

R∆

0

1

2

3

Projects...

Tilman Plehn

MadGOLEM

SHERPA

HEPTopTagger

SUSY

Higgs

HEPTopTagger

Boosted top quarks [Kaplan, Rehermann, Schwartz, Tweedie; Princeton, Seattle...]

– hadronic tops with pT & 800 GeVisolation and b tagging challenging

– C/A algorithm with pT drop criterionall top decay jets identified3 kinematic constraints: mW ,mt , cos θhel [no b tag]

– top mass included, no sidebins

HEPTopTagger [TP, Salam, Spannowsky, Takeuchi]

– extend to lower pT & 250 GeVrealistic for t t in Standard Model

– top reconstruction possible

– tested and implemented by ATLAS [Kasieczka & Schatzel]

– hadronic top like tagged b

[GeV]T

p200 300 400 500 600

Eff

icie

ncie

s

0

0.2

0.4

0.6

0.8

1inside fat jet

tagged

tagged (unmatched)

Projects...

Tilman Plehn

MadGOLEM

SHERPA

HEPTopTagger

SUSY

Higgs

Stops

Stop pairs [TP, Spannowsky, Takeuchi, Zerwas] t

t– stop most important for hierarchy problemcomparison to other top partners [Meade & Reece]

– difficult semi-leptonic channel

– purely hadronic: t t∗ → tχ01 tχ0

1 [CMS TDR: leptons as spontaneous life guards]

events in 1 fb−1 t1 t∗1 t t QCD W+jets Z+jets S/B S/√

B10 fb−1

mt [ GeV] 340 390 440 490 540 640 340

pT ,j > 200 GeV, ` veto 728 447 292 187 124 46 87850 2.4 · 107 1.6 · 105 n/a 3.0 · 10−5

/ET > 150 GeV 283 234 184 133 93 35 2245 2.4 · 105 1710 2240 1.2 · 10−3

first top tag 100 91 75 57 42 15 743 7590 90 114 1.2 · 10−2

second top tag 15 12.4 11 8.4 6.3 2.3 32 129 5.7 1.4 8.3 · 10−2

b tag 8.7 7.4 6.3 5.0 3.8 1.4 19 2.6 . 0.2 . 0.05 0.40 5. 9mT 2 > 250 GeV 4.3 5.0 4.9 4.2 3.2 1.2 4.2 . 0.6 . 0.1 . 0.03 0.88 6. 1

Projects...

Tilman Plehn

MadGOLEM

SHERPA

HEPTopTagger

SUSY

Higgs

Stops

Stop pairs [TP, Spannowsky, Takeuchi, Zerwas] t

t– stop most important for hierarchy problemcomparison to other top partners [Meade & Reece]

– difficult semi-leptonic channel

– purely hadronic: t t∗ → tχ01 tχ0

1 [CMS TDR: leptons as spontaneous life guards]

– stop mass from mT2 endpoint [like sleptons or sbottoms]

– not even a hard analysis

[GeV]T2M200 300 400 500 600

[/20

GeV

]T

2dM

σdσ1

0

0.2

0.4

=340 GeVstopM

=540 GeVstopM

tt

=98 GeVLSPtestM

Projects...

Tilman Plehn

MadGOLEM

SHERPA

HEPTopTagger

SUSY

Higgs

Measuring JCP using jet substructure

0.0 0.5 1.0 1.5 2.0 2.5 3.0

Φ, full reconstruction, Herwig++

0.0000

0.0005

0.0010

0.0015

0.0020

0.0025

dσ/d

Φ [

pb/

0.1

570

5 ]

Scalar 0+

Scalar 0−

Vector 1−

Vector 1+

Tensor 2+

tt

ZZ

Z+Jet

Φ

µ−

µ+jα

jβθ�

θh

X

Z Z

p

p

ez

ez�

θ�

Φ1

1.0 0.5 0.0 0.5 1.0

cosθl, full reconstruction, Herwig++

0.0000

0.0005

0.0010

0.0015

0.0020

0.0025

dσ/d

cosθl [

pb/

0.1

]

Scalar 0+

Scalar 0−

Vector 1−

Vector 1+

Tensor 2+

tt

ZZ

Z+Jet

[Cabibbo, Maksymowicz ’65] [Dell’Aquila, Nelson ’85]

[Buszello, Fleck, Marquard, van der Bij ’04] . . .

I sub-jets plus spin analysissemi-hadronic X → ZZ → µ+µ−jj

I price for boosted signatures(many angles lost)

I discriminative powersingly-produced 0± “Higgs lookalikes”

[Englert, Hackstein, Spannowsky ’10]

Projects...

Tilman Plehn

MadGOLEM

SHERPA

HEPTopTagger

SUSY

Higgs

Outlook

MC development and LHC applications

– MadGOLEM

– SHERPA

– misc tools: SFitter, HEPTopTagger, FeynRulez, MC@NLO,...

– NLO/matching/merging for new physics

– bottom partons

– Higgs plus QCD

– SUSY Higgs predictions

– stop searches

– jets plus missing energy,...

thank you to:

Christoph Englert, Steffen Schumann, Michihisa Takeuchi, David Lopez Val,Erik Gerwick, Dorival Netto, Ioan Wigmore,Fabian Gross, Thomas Lubbert, Peter Schichtel,Kentarou Mawatari, Michael Rauch, Michael Spannowsky, Carole Weydert,...

Projects...

Tilman Plehn

MadGOLEM

SHERPA

HEPTopTagger

SUSY

Higgs