Towards Precision Jet Mass Calculations Randall S. Kelley Frontiers - - PowerPoint PPT Presentation

towards precision jet mass calculations
SMART_READER_LITE
LIVE PREVIEW

Towards Precision Jet Mass Calculations Randall S. Kelley Frontiers - - PowerPoint PPT Presentation

Towards Precision Jet Mass Calculations Randall S. Kelley Frontiers in QCD (INT-11-3) Oct 5, 2011 R. S. Kelley (Harvard) Jet Mass Oct 5, 2011 1 / 40 References Resummation of jet mass with a jet veto arXiv:1102.0561v2 RK, Matthew D.


slide-1
SLIDE 1

Towards Precision Jet Mass Calculations

Randall S. Kelley Frontiers in QCD (INT-11-3) Oct 5, 2011

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 1 / 40

slide-2
SLIDE 2

References

Resummation of jet mass with a jet veto arXiv:1102.0561v2

RK, Matthew D. Schwartz, Hau Xing Zhu

The two-loop hemisphere soft function arXiv:1105.3676

RK, Robert M. Schabinger, Matthew D. Schwartz, Hau Xing Zhu

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 2 / 40

slide-3
SLIDE 3

Outline

1

Introduction

2

2-loop Hemisphere Soft function

3

Inclusive R dependent Jet Shapes

4

Exclusive Jet Masses

5

Factorization of the Soft Function

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 3 / 40

slide-4
SLIDE 4

Outline

1

Introduction

2

2-loop Hemisphere Soft function

3

Inclusive R dependent Jet Shapes

4

Exclusive Jet Masses

5

Factorization of the Soft Function

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 4 / 40

slide-5
SLIDE 5

Introduction

Very large number of jets at the LHC. Jets provide a wealth of information about QCD and exploring new physics.

  • excess in the number of jets could be a sign of new physics

Substructure may be critical in new physics searches.

  • massive boosted heavy particles can be found in jet

Jet rate distributions have been calculated to NLO, but little has been said about structure of jets (i.e. m2, R, angularity, etc.). Predictions may be spoiled by large logarithms ( logn m1

m2 , log R, etc)

Effective field theories provide a way to systematically improve calculations.

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 5 / 40

slide-6
SLIDE 6

Factorization (a preview)

τ = 1 − T ≈ m2

L + m2 R

Q2 dσ dm2

Ldm2 R

∼H(Q, µh)

  • dkLdkR

× J(m2

L − kLQ, µj)J(m2 R − kRQ, µj)S(kL, kR, µs)

( Fleming et al., Schwartz) Factorization is achieved using Soft Collinear Effective theory (SCET) Use the LO results in SCET to predict the NLO singular piece using renormalization group evolution (RGE). Compare α2

s results to EVENT2

(Catani and Seymore)

QCD SCET

0.0 0.1 0.2 0.3 0.4 0.5 5 10 15 20

Τ

Τ Σ0 dΣ dΤ

Thrust LO

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 6 / 40

slide-7
SLIDE 7

Factorization (a preview)

τ = 1 − T ≈ m2

L + m2 R

Q2 dσ dm2

Ldm2 R

∼H(Q, µh)

  • dkLdkR

× J(m2

L − kLQ, µj)J(m2 R − kRQ, µj)S(kL, kR, µs)

( Fleming et al., Schwartz) Factorization is achieved using Soft Collinear Effective theory (SCET) Use the LO results in SCET to predict the NLO singular piece using renormalization group evolution (RGE). Compare α2

s results to EVENT2

(Catani and Seymore)

QCD SCET

0.0 0.1 0.2 0.3 0.4 0.5 50 100 150 200 250 300

Τ

Τ Σ0 dΣ dΤ

Thrust NLO

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 6 / 40

slide-8
SLIDE 8

Factorization (a preview)

τ = 1 − T ≈ m2

L + m2 R

Q2 dσ dm2

Ldm2 R

∼H(Q, µh)

  • dkLdkR

× J(m2

L − kLQ, µj)J(m2 R − kRQ, µj)S(kL, kR, µs)

( Fleming et al., Schwartz) Factorization is achieved using Soft Collinear Effective theory (SCET) Use the LO results in SCET to predict the NLO singular piece using renormalization group evolution (RGE). Compare α2

s results to EVENT2

(Catani and Seymore)

12 10 8 6 4 2 14000 12000 10000 8000 6000 4000 2000

Τ

Τ Σ0 dΣ dΤ

Thrust NLO

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 6 / 40

slide-9
SLIDE 9

What could go wrong?

The factorization theorem is valid only when mL and mR are small (i.e. the small m2 region is dominated by IR degrees of freedom) SCET does not guarantee log m2

L/m2 R are resummed by RGE

(can be calculated by brute force) Produce non-global logarithms ( Dasgupta and Salem) −CF CA αs 4π 2 16π2 3 log2 m2

L

Q2

  • Hard emissions are not included in SCET degrees of

freedom (type 1) Sharply divided phase space with separated scales mL ≪ mR (type 2) Finite jet size (R), and cutoff scales (Eout < ω) complicate the problem considerably.

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 7 / 40

slide-10
SLIDE 10

What could go wrong?

Σ(ρR) = ∞ dm1 ρRQ2 dm2

2

d2σ dm2

1dm2 2

Produce non-global logarithms ( Dasgupta and Salem) −CF CA αs 4π 2 16π2 3 log2 m2

L

Q2

  • Hard emissions are not included in SCET degrees of

freedom (type 1) Sharply divided phase space with separated scales mL ≪ mR (type 2) Finite jet size (R), and cutoff scales (Eout < ω) complicate the problem considerably.

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 7 / 40

slide-11
SLIDE 11

What could go wrong?

Σ(ρR) = ∞ dm1 ρRQ2 dm2

2

d2σ dm2

1dm2 2

Produce non-global logarithms ( Dasgupta and Salem) −CF CA αs 4π 2 16π2 3 log2 m2

L

Q2

  • Hard emissions are not included in SCET degrees of

freedom (type 1) Sharply divided phase space with separated scales mL ≪ mR (type 2) Finite jet size (R), and cutoff scales (Eout < ω) complicate the problem considerably.

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 7 / 40

slide-12
SLIDE 12

Main Points

Seek to understand non-global logarithms and how to control them. Understand how different jet shapes and jet sizes (R) affect the observables. Consider first inclusive and then exclusive observables. We perform resummation for a 2-jet observable with jets of size R. τω = m2

1 + m2 2

Q2 , E3 < ω Demonstration of factorization of the soft function: SR(k, ω, µ) = Sin

R (k, µ)Sout R (ω, µ)

and discuss limitations.

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 8 / 40

slide-13
SLIDE 13

Outline

1

Introduction

2

2-loop Hemisphere Soft function

3

Inclusive R dependent Jet Shapes

4

Exclusive Jet Masses

5

Factorization of the Soft Function

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 9 / 40

slide-14
SLIDE 14

Hemisphere Jets

try calculating: dσ dm2

Ldm2 R

∼H(Q, µh)

  • dkLdkR

× J(m2

L − kLQ, µj)J(m2 R − kRQ, µj)S(kL, kR, µs)

RG evolution only resums log m2 Q2 , but does not say anything about log m2

L

m2

R

. These logs come from kL/kR in the soft function.

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 10 / 40

slide-15
SLIDE 15

Calculation of two-loop Soft function

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 11 / 40

slide-16
SLIDE 16

Calculation of two-loop Soft function

S(kL, kR) = δ(kL)δ(kR) + αs 4π S(1)(kL, kR, µ) + α2

s

16π2 S(2)(kL, kR, µ) + · · · NLO result S(2) = C2

F SCF + CF CASCA + CF nfTF Snf

S = µ4ǫ (kLkR)1+4ǫ f kL kR , ǫ

  • +

µ2ǫ k1+2ǫ

R

δ(kL) + µ2ǫ k1+2ǫ

L

δ(kR)

  • g(ǫ)

There is a different f(r, ǫ) and g(ǫ) for each color factor, where r = kL/kR. f(r, ǫ) was calculated independently by (Hornig et al. 1105.4628)

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 12 / 40

slide-17
SLIDE 17

Cumulative the Soft Function

Terms of the form µaǫ k1+aǫ , a = 2, 4, must be thought of as distributions and integrated. R(X, Y, µ) = X dkL Y dkR S(kL, kR, µ) Result is used for integrated heavy jet mass and thrust distributions. The singular parts of the thrust and heavy jet mass distributions can be extracted (previously only known numerically ) 1 σ0 dσ dτ = δ(τ)D(τ)

δ

+ αs 4π [D(1)(τ)]+ + αs 4π 2 [D(2)(τ)]+ + · · · Removes a source of theoretical uncertainty in N3LL result for heavy jet mass, improving fits to αs.

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 13 / 40

slide-18
SLIDE 18

Cumulative the Soft Function

Terms of the form µaǫ k1+aǫ , a = 2, 4, must be thought of as distributions and integrated. R(X, Y, µ) = X dkL Y dkR S(kL, kR, µ) Result is used for integrated heavy jet mass and thrust distributions. The singular parts of the thrust and heavy jet mass distributions can be extracted (previously only known numerically ) D(τ)

δ

= αs 4π 2

  • − 3π4

10 C2

F + CF CA

638ζ3 9 − 335π2 54 + 22π4 45 − 2140 81

  • + CF TF nf
  • −232ζ3

9 + 74π2 27 + 80π2 81 Removes a source of theoretical uncertainty in N3LL result for heavy jet mass, improving fits to αs.

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 13 / 40

slide-19
SLIDE 19

Asymptotic behavior and Non-Global Logarithms

Non-global logs must come from the µ-independent part of the soft function. R(X, Y, µ) = Rµ X µ , Y µ

  • + Rf

X Y

  • for z = X

Y ≫ 1,

Rz≫1

f

(z) = π4 2 C2

F +

8 3 − 16π2 9

  • | log z| + − 136

81 + 154π2 27 + 184ζ3 9

  • CF nfTF

+

  • − 4

3π2log2 z +

  • −8ζ3 − 4

3 + 44π2 9

  • | log z| − 506ζ3

9 + 8π4 5 − 871π2 54 − 2032 81

  • CF CA
  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 14 / 40

slide-20
SLIDE 20

Asymptotic behavior and Non-Global Logarithms

Non-global logs must come from the µ-independent part of the soft function. R(X, Y, µ) = Rµ X µ , Y µ

  • + Rf

X Y

  • for z = X

Y ∼ 1,

Rz∼1

f

(z) = π4 2 C2

F +

  • − 2

3 − 4π2 3 − 4 log2 2 + 44 log 2 3

  • log2 z − 32Li4

1 2

  • + 88ζ3

9 −28ζ3 log(2) − 2032 81 − 871π2 54 + 16π4 9 − 4 log4 2 3 + 4 3π2log2 z

  • CF CA

+ 4 3 − 16 log 2 3

  • log2 z + 154π2

27 − 136 81 − 32ζ3 9

  • CF nfTF + O(log3 z).

Hoang-Kluth ansatz (0806.3852) only valid at small log z.

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 14 / 40

slide-21
SLIDE 21

Asymptotic behavior and Non-Global Logarithms

Non-global logs must come from the µ-independent part of the soft function. R(X, Y, µ) = Rµ X µ , Y µ

  • + Rf

X Y

  • 10

5 5 10 40 20 20 40 60 80

ln z ln X

Y 2 ln ML MR

Rf CF nf TF

6 4 2 2 4 6 250 200 150 100 50

ln z ln X

Y 2 ln ML MR

Rf CF CA

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 14 / 40

slide-22
SLIDE 22

Outline

1

Introduction

2

2-loop Hemisphere Soft function

3

Inclusive R dependent Jet Shapes

4

Exclusive Jet Masses

5

Factorization of the Soft Function

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 15 / 40

slide-23
SLIDE 23

Jet definition

Cambridge/Aachem algorithm

  • Assign Rij = 1

2(1 − cos θij) to each pair of particles

  • For the smallest value of Rij, merge the four vectors of the pair if Rij < R.
  • Repeat until there are no pairs with Rij < R, then stop.

Order jets by energy, E1 > E2 > E3 > · · · Veto events with E3 > ω if interested in dijets.

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 16 / 40

slide-24
SLIDE 24

Inclusive R dependent Jet Shapes

τA = mpri2 + m2

pri

Q2 τA ≪ 1 forces dijets R-dependent jet shape (log R’s) Very sensitive to the choice of the primary jet, sometimes not well defined. may be useful in Hadron colliders (dynamical threshold enhancement)

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 17 / 40

slide-25
SLIDE 25

Inclusive R dependent Jet Shapes

τA = mpri2 + m2

pri

Q2 τA ≪ 1 forces dijets R-dependent jet shape (log R’s) Very sensitive to the choice of the primary jet, sometimes not well defined. may be useful in Hadron colliders (dynamical threshold enhancement)

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 17 / 40

slide-26
SLIDE 26

Try τA1

τA1 = m2

1 + m2 ¯ 1

Q2 Same as Thrust at O(αs) use Jinc(p2) Soft function depends critically on R.

Yq g Yqg Yq q

0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0

t u

QCD LO phase space

Soft function is insensitive to jet energy If soft gluon is not within R of either jet, which jet is most energetic is ambiguous.

1 σ0 dσQCD dτA1 = δ(τAq) + αs 2π CF

  • −1 + π2

3

  • δ(τA1) + αs

2π CF −4 log τA1 − 3 τA1

  • +
  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 18 / 40

slide-27
SLIDE 27

Try τA1

τA1 = m2

1 + m2 ¯ 1

Q2 Same as Thrust at O(αs) use Jinc(p2) Soft function depends critically on R.

0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0

k k

Soft phase space

Soft function is insensitive to jet energy If soft gluon is not within R of either jet, which jet is most energetic is ambiguous.

1 σ0 dσQCD dτA1 = δ(τAq) + αs 2π CF

  • −1 + π2

3

  • δ(τA1) + αs

2π CF −4 log τA1 − 3 τA1

  • +
  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 18 / 40

slide-28
SLIDE 28

Try τAq

τAq = m2

q + m2 ¯ q

Q2 define quark jet to be “primary”

Sin

R (k, µ) = δ(k)

+ αs 2π CF

  • − log2

R 1 − R + π2 6

  • δ(k)

+ αs 2π CF −8 log k

µ + 4 log R 1−R

k [k,µ]

+

Sin

1−R(k, µ) for other jet

0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0

k k

Soft phase space

1 σ0 dσSCET dτAq = δ(τAq) + αs 2π CF

  • −1 + π2

3 + log2 R 1 − R

  • δ(τAq) + αs

2π CF

  • −4 log τAq − 3

τAq

  • +

at NLO, find negative cross sections since it’s not IR safe

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 19 / 40

slide-29
SLIDE 29

Try τAq

τAq = m2

q + m2 ¯ q

Q2 define quark jet to be “primary”

Sin

R (k, µ) = δ(k)

+ αs 2π CF

  • − log2

R 1 − R + π2 6

  • δ(k)

+ αs 2π CF −8 log k

µ + 4 log R 1−R

k [k,µ]

+

Sin

1−R(k, µ) for other jet

0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0

k k

Soft phase space

1 σ0 dσSCET dτAq = δ(τAq) + αs 2π CF

  • −1 + π2

3 + log2 R 1 − R

  • δ(τAq) + αs

2π CF

  • −4 log τAq − 3

τAq

  • +

at NLO, find negative cross sections since it’s not IR safe

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 19 / 40

slide-30
SLIDE 30

Try an average

dσ dτA1 + dσ dτA2 = dσ dτAq + dσ dτA¯

q

Agrees with QCD at LO at NLO:

R 0.1 CF

2

CFCA CFn f TF

10 8 6 4 2 200 100 100 200 ln ΤA ΤAQCDSCET

R0.01 R0.1 R0.2 R0.3 x 2 x 4 x 8

10 8 6 4 2 1500 1000 500 500 1000 1500 2000 ln ΤA ΤAQCD SCET

SCET can resum all large log τA’s, for any R For small R, these may not be dominant part. Have not attempted to resum log R’s.

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 20 / 40

slide-31
SLIDE 31

Outline

1

Introduction

2

2-loop Hemisphere Soft function

3

Inclusive R dependent Jet Shapes

4

Exclusive Jet Masses

5

Factorization of the Soft Function

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 21 / 40

slide-32
SLIDE 32

Exclusive Jet Masses

Veto events with E3 > ω Trivial dependence on m2 at LO Clustered jet always has the most energy fω(R) vanishes at R → 1/2 (hemisphere case)

1 σ0

  • d2σ

dm2

1dm2 2

  • QCD

= δ(m2

1)δ(m2 2) + α

4π CF δ(m2

2)

×

  • −2 + 2π2

3 −8 log R 1 − R log 2ω Q + fω(R)

  • δ(m2

1)

+    −6 + 8 log

R 1−R − 8 log m2

1

Q2

m2

1

  

+ · · ·     

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 22 / 40

slide-33
SLIDE 33

Exclusive Jet Masses

Veto events with E3 > ω Trivial dependence on m2 at LO Clustered jet always has the most energy fω(R) vanishes at R → 1/2 (hemisphere case)

Yq g Yqg Yq q

0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0

t u

QCD LO phase space

1 σ0

  • d2σ

dm2

1dm2 2

  • QCD

= δ(m2

1)δ(m2 2) + α

4π CF δ(m2

2)

×

  • −2 + 2π2

3 −8 log R 1 − R log 2ω Q + fω(R)

  • δ(m2

1)

+    −6 + 8 log

R 1−R − 8 log m2

1

Q2

m2

1

  

+ · · ·     

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 22 / 40

slide-34
SLIDE 34

Exclusive Jet Masses

dσ dm2

qdm2 ¯ q

∼H(Q, µh)

  • dkqdk¯

q

× J(m2

q − kqQ, µj)J(m2 ¯ q − k¯ qQ, µj)SR(kq, k¯ q, ω, µs)

At order αs: SR(kL, kR, ω, µ) = Sin

R (kL, µ)Sin R (kR, µ)Sout R (ω, µ)

Sin

R (k, µ) = δ(k) + αs

4π CF

  • −2 log2

R 1 − R + π2 3

  • δ(k)

+ αs 4π CF −16 log k

µ − 8 log R 1−R

k [k,µ]

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 23 / 40

slide-35
SLIDE 35

Exclusive Jet Masses

dσ dm2

qdm2 ¯ q

∼H(Q, µh)

  • dkqdk¯

q

× J(m2

q − kqQ, µj)J(m2 ¯ q − k¯ qQ, µj)SR(kq, k¯ q, ω, µs)

At order αs: SR(kL, kR, ω, µ) = Sin

R (kL, µ)Sin R (kR, µ)Sout R (ω, µ)

Sout

R (ω, µ) = 1 + αs

4π CF

  • −8 log

R 1 − R log 2ω µ + 2 log2 R 1 − R + f0(R) [k,µ]

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 23 / 40

slide-36
SLIDE 36

QCD vs SCET

1 σ0

  • d2σ

dm2

1dm2 2

  • QCD

= δ(m2

1)δ(m2 2) + α

4π CF δ(m2

2)

×

  • −2 + 2π2

3 −8 log R 1 − R log 2ω Q + fω(R)

  • δ(m2

1) +

   −6 + 8 log

R 1−R − 8 log m2

1

Q2

m2

1

  

     Combining the soft function with the hard and inclusive jet functions, we get 1 σ0 d2σ dm2

qdm2 ¯ q

= δ(m2

q)δ(m2 ¯ q) + α

4π CF

  • −2 + 2π2

3 − 8 log R 1 − R log 2ω Q + f0(R)

  • δ(m2

q)δ(m2 ¯ q)

+    −6 + 8 log

R 1−R − 8 log m2

q

Q2

2m2

q

  

δ(m2

¯ q) +

   −6 + 8 log

R 1−R − 8 log m2

¯ q

Q2

2m2

¯ q

  

δ(m2

q)

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 24 / 40

slide-37
SLIDE 37

QCD vs SCET

δ(mq)δ(m¯

q) matches δ(mq)δ(m2) with f0(R) instead of fω(R).

SCET is symmetric mq ↔ m¯

q, QCD is not

Mass of the hardest jet is not simply related to any projection of

d2σ dm2

qdm2 ¯ q

dσ dm2

  • QCD

= Q2R dm2

1

Q2R dm2

2

d2σ dm2

1dm2 2

× 1 2

  • δ(m2 − m2

1) + δ(m2 − m2 2)

  • Will have NGLs of form logn m1

m2

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 25 / 40

slide-38
SLIDE 38

Define τω

Veto events with E3 > ω Define: τω = m2

1 + m2 2

Q2 Reproduces QCD as τω → 0 Avoids NGLs of form logn m1 m2

1 σ0 dσQCD dτω = δ(τω) + αs 2π CF

  • 7 − 5π2

6 + 4 log 1 − R R log 2ω Q + fω(R)

  • δ(τω)

+ αs 2π CF

  • −4 log τω − 3 − 4 log 1−R

R

τω

  • +

+ · · · 1 σ0 dσSCET dτω = δ(τω) + αs 2π CF

  • 7 − 5π2

6 + 4 log 1 − R R log 2ω Q + f0(R)

  • δ(τω)

+ αs 2π CF

  • −4 log τω − 3 − 4 log 1−R

R

τω

  • +

+ · · ·

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 26 / 40

slide-39
SLIDE 39

Define τω

Veto events with E3 > ω Define: τω = m2

1 + m2 2

Q2 Reproduces QCD as τω → 0 Avoids NGLs of form logn m1 m2

Yq g Yqg Yq q

0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0

t u

QCD LO phase space

1 σ0 dσQCD dτω = δ(τω) + αs 2π CF

  • 7 − 5π2

6 + 4 log 1 − R R log 2ω Q + fω(R)

  • δ(τω)

+ αs 2π CF

  • −4 log τω − 3 − 4 log 1−R

R

τω

  • +

+ · · · 1 σ0 dσSCET dτω = δ(τω) + αs 2π CF

  • 7 − 5π2

6 + 4 log 1 − R R log 2ω Q + f0(R)

  • δ(τω)

+ αs 2π CF

  • −4 log τω − 3 − 4 log 1−R

R

τω

  • +

+ · · ·

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 26 / 40

slide-40
SLIDE 40

Outline

1

Introduction

2

2-loop Hemisphere Soft function

3

Inclusive R dependent Jet Shapes

4

Exclusive Jet Masses

5

Factorization of the Soft Function

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 27 / 40

slide-41
SLIDE 41

Factorization of the Soft Function

SR(k, ω, µ) = Sin

R (k, µ)Sout R (ω, µ)SF R

ω k

  • All NGL’s are in SF

R

ω

k

  • When can we neglect ω/k dependence?
  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 28 / 40

slide-42
SLIDE 42

Factorization of the Soft Function

For small R, we can show SR(k, ω, µ) = Sin

R (k, µ)Sout R (ω, µ)

for ω/Q k/Q ≪ R ≪ 1 Later, we discuss log ω/k terms which violate this factorization. SCET requires ω, k to be small, but they can be far apart For small R, k is in the cone and has collinear scaling. k+ < R 1 − Rk− (k+, k−, k⊥) ∼ k R(R, 1, √ R) q is outside of either cone with Eq < ω. (q+, q−, q⊥) ∼ (ω, ω, ω)

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 28 / 40

slide-43
SLIDE 43

Factorization of the Soft Function

diagrammatic proof For small R, and ω kR

  • Mi =
  • g2 nµnν

k−q− T bT a − g2 ¯ nµnν k−q+ T bT a

  • εa

µ(q)εb ν(k)

Equivalent to the the following refactorization Y †

¯ n Yn → (Y sc ¯ n )†(Y us ¯ n )†(Y us n )(Y sc n ).

Similar to the factorization using SCET+ in Bauer et al. 1106.6047

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 29 / 40

slide-44
SLIDE 44

Soft Anomalous dimension

Γs = αs π CF Γcusp log kLkR µ2 + γout

S

+ γin

S

Extract γS from the αs calculation γout

S

= − αs 4π CF Γcusp log R 1 − R γin

S = γhemi S

+ αs 4π CF Γcusp log R 1 − R RG invariance requires the R dependence to cancel in the sum to all orders Ellis et al. 0912.062, JHEP 1011,101 (2010) Holds at two loops, suspect it holds at all orders. Refactorization gives predictive power through separating scales As R → 1

2, γin S → γhemi S

and γout

S

→ 0. At order α2

s, this form contributes terms to the expression

Γ1 log R 1 − R log τω

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 30 / 40

slide-45
SLIDE 45

Soft Anomalous dimension

Γs = αs π CF Γcusp log kLkR µ2 + γout

S

+ γin

S

Extract γS from the αs calculation γout

S

= − αs 4π CF Γcusp log R 1 − R−γR(R) γin

S = γhemi S

+ αs 4π CF Γcusp log R 1 − R+γR(R) γR(R) should approach a constant in small R The structure of γR(R) is not known beyond 1 loop

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 31 / 40

slide-46
SLIDE 46

Predictions from Refactorization

Now consider ω ≁ τωQ SCET agrees with QCD up to powers in ω/Q and τω (brute force if necessary) Neglecting powers of ω/τωQ is consistent with numerics. Could be important log

ω τωQ terms

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 32 / 40

slide-47
SLIDE 47

Predictions from Refactorization

When R is not small, “in” jet radiation is not small and there is no obvious factorization. R → 1

2 (hemisphere case), the ω dependence vanishes

Factorization captures the log R log τωQ

2ω , but not the terms constant in R wrong.

The factorization holds at small and large R and is a good approximation for moderate R. Much of the R dependence of full QCD is captured by the small R limit.

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 33 / 40

slide-48
SLIDE 48

Thrust-like jets

Find thrust axis Cluster particles within R of thrust axis Same as Cambridge/Aachem at αs, similar at α2

s

NGL’s structure is different than CA (Hornig et al. 1110.0004)

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 34 / 40

slide-49
SLIDE 49

Numerical Check of Ansatz

The α2

s predictions from SCET where compared to EVENT2

(Catani and Seymore) Checked both Cambridge/Aachem jets and Thrust-like jets We expect SCET to agree with EVENT2 up to powers in τω and ω/Q. Highly non-trivial check of the factorization theorem Holds independently various color factors C2

F , CACF and CF nfTF .

Checked for a large range of R values

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 35 / 40

slide-50
SLIDE 50

Comparison with EVENT2 QCD single scale

0.00 0.02 0.04 0.06 0.08 100 50 50 100 150 200 250

ΤΩ

ΤΩ Σ0 dΣ dΤΩ

NLO

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 36 / 40

slide-51
SLIDE 51

Comparison with EVENT2 QCD single scale two scales

0.00 0.02 0.04 0.06 0.08 100 50 50 100 150 200 250

ΤΩ

ΤΩ Σ0 dΣ dΤΩ

NLO

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 36 / 40

slide-52
SLIDE 52

Comparison with EVENT2 QCD single scale two scales with 1

0.00 0.02 0.04 0.06 0.08 100 50 50 100 150 200 250

ΤΩ

ΤΩ Σ0 dΣ dΤΩ

NLO

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 36 / 40

slide-53
SLIDE 53

Comparison with EVENT2

R 0.1, Ω 104Q

With refactorization

CF

2

CFCA CFn f TF

10 8 6 4 2 500 500 ln ΤΩ ΤΩQCDSCET

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 36 / 40

slide-54
SLIDE 54

Comparison with EVENT2

R 0.1, Ω 104Q

With refactorization and cusp ansatz

CF

2

CFCA CFn f TF

10 8 6 4 2 500 500 ln ΤΩ ΤΩQCDSCET

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 36 / 40

slide-55
SLIDE 55

Comparison with EVENT2: Thrust-Axis clustering

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 37 / 40

slide-56
SLIDE 56

Comparison with EVENT2: Thrust-Axis clustering

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 37 / 40

slide-57
SLIDE 57

Comparison with EVENT2: Thrust-Axis clustering

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 37 / 40

slide-58
SLIDE 58

Comparison with EVENT2: Thrust-Axis clustering

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 37 / 40

slide-59
SLIDE 59

Comparison with EVENT2: Thrust-Axis clustering

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 37 / 40

slide-60
SLIDE 60

Comparison with EVENT2: Thrust-Axis clustering

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 37 / 40

slide-61
SLIDE 61

Comparison with EVENT2: Thrust-Axis clustering

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 37 / 40

slide-62
SLIDE 62

Comparison with EVENT2: Thrust-Axis clustering

too far!

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 37 / 40

slide-63
SLIDE 63

Comparison with EVENT2 Cambridge/Aachem clustering

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 38 / 40

slide-64
SLIDE 64

Comparison with EVENT2 Cambridge/Aachem clustering

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 38 / 40

slide-65
SLIDE 65

Comparison with EVENT2 Cambridge/Aachem clustering

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 38 / 40

slide-66
SLIDE 66

Comparison with EVENT2 Cambridge/Aachem clustering

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 38 / 40

slide-67
SLIDE 67

Comparison with EVENT2 Cambridge/Aachem clustering

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 38 / 40

slide-68
SLIDE 68

Comparison with EVENT2 Cambridge/Aachem clustering

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 38 / 40

slide-69
SLIDE 69

Comparison with EVENT2 Cambridge/Aachem clustering

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 38 / 40

slide-70
SLIDE 70

Two loop τω Soft Function preliminary

General form of two loop soft function SR(k, ω, µ) = Sin

R

k µ

  • Sout

R

ω µ

  • SF

R

k ω

  • SF

R = SNGL R

+ finite consider cumulative distribution instead (τω → σ) CF nfTF channel R → 0 then ΣQ 2ω → 0 (recall Σ < R)

  • −32π2

9 + 16 3

  • log ΣQ

1010 106 0.01 100 106 1010 600 400 200 200 Q 2 Ω SNGL R0.1 R0.01 R0.001 R0.0001

NGL in R → 0 is twice hemisphere case

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 39 / 40

slide-71
SLIDE 71

Two loop τω Soft Function preliminary

General form of two loop soft function SR(k, ω, µ) = Sin

R

k µ

  • Sout

R

ω µ

  • SF

R

k ω

  • SF

R = SNGL R

+ finite consider cumulative distribution instead (τω → σ) CF nfTF channel ΣQ 2ω → 0 then R → 0

  • −32π2

9 + 16 3

  • log ΣQ

2ωR2

1010 106 0.01 100 106 1010 600 400 200 200 Q 2 Ω SNGL R0.1 R0.01 R0.001 R0.0001

NGL in R → 0 is twice hemisphere case

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 39 / 40

slide-72
SLIDE 72

Two loop τω Soft Function preliminary

General form of two loop soft function SR(k, ω, µ) = Sin

R

k µ

  • Sout

R

ω µ

  • SF

R

k ω

  • SF

R = SNGL R

+ finite consider cumulative distribution instead (τω → σ) CF nfTF channel ΣQ 2ω → ∞, no R dependence. −

  • −32π2

9 + 16 3

  • log ΣQ

2ω NGL in R → 0 is twice hemisphere case

1010 106 0.01 100 106 1010 600 400 200 200 Q 2 Ω SNGL R0.1 R0.01 R0.001 R0.0001

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 39 / 40

slide-73
SLIDE 73

Two loop τω Soft Function preliminary

General form of two loop soft function SR(k, ω, µ) = Sin

R

k µ

  • Sout

R

ω µ

  • SF

R

k ω

  • SF

R = SNGL R

+ finite consider cumulative distribution instead (τω → σ) CF CA channel NGL agrees with Hornig et al. 1110.0004

  • −8π2

3 + 16Li2

  • R2

(1 − R)2

  • log2 ΣQ

2ω +

  • −16ζ3 − 8

3 + 88π2 9 + · · ·

  • log ΣQ

0.0 0.1 0.2 0.3 0.4 0.5 25 20 15 10 5

R fNGLR

NGL in R → 0 is twice hemisphere case

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 39 / 40

slide-74
SLIDE 74

Two loop τω Soft Function preliminary

General form of two loop soft function SR(k, ω, µ) = Sin

R

k µ

  • Sout

R

ω µ

  • SF

R

k ω

  • SF

R = SNGL R

+ finite consider cumulative distribution instead (τω → σ) CF CA channel NGL agrees with Hornig et al. 1110.0004

  • −8π2

3 + 16Li2

  • R2

(1 − R)2

  • log2

ΣQ 2ωR2 +

  • −16ζ3 − 8

3 + 88π2 9 + · · ·

  • log ΣQ

2ωR2

0.0 0.1 0.2 0.3 0.4 0.5 25 20 15 10 5

R fNGLR

Possible log R dependence in leading NGL missed.

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 39 / 40

slide-75
SLIDE 75

Conclusions

Inclusive observables (e.g. τA) seemed amenable to resummation Soft function factorization held in limit ω/Q τω ≪ R ≪ 1 but was not a bad approximation elsewhere. Non-global structures are present, but numerically small for a large choice of parameters The results extrapolated away from R → 0 limit provides good agreement with QCD. Calculation of τω soft function almost finished.

  • R. S. Kelley (Harvard)

Jet Mass Oct 5, 2011 40 / 40