Review of accelerator data of relevance to air shower simulations - - PowerPoint PPT Presentation

review of accelerator data of relevance to air shower
SMART_READER_LITE
LIVE PREVIEW

Review of accelerator data of relevance to air shower simulations - - PowerPoint PPT Presentation

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb2012 Review of accelerator data of relevance to air shower simulations Yoshitaka Itow STE Lab / Kobayashi-Maskawa Inst. Nagoya University UHECR 2012 Feb 13-16, 2012, CERN 1 Y.Itow,


slide-1
SLIDE 1

1

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb2012

Review of accelerator data

  • f relevance to air shower simulations

Yoshitaka Itow STE Lab / Kobayashi-Maskawa Inst. Nagoya University

“UHECR 2012” Feb 13-16, 2012, CERN

slide-2
SLIDE 2

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

2

Hadron interactions at ultra high energy Accelerator Cosmic rays

Precision improvement Hint for interactions at ultra-ultra high energy

ECM ~ ( 2 × Elab × Mp ) 1/2

√s=14TeV collision at LHC

1017eV cosmic rays

slide-3
SLIDE 3

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

3

Cosmic ray spectrum & historical colliders

LHC 14TeV Tevatron LHC 0.9TeV LHC 7 TeV SppS RHIC ISR

1010 1020 eV

>40 yrs legacy of wisdom for interactions available !

AUGER, TA TALE HEAT

slide-4
SLIDE 4

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

4

① Inelastic cross section ② Forward energy spectrum

If large k rapid development If small k deep penetrating If large σ rapid development If small σ deep penetrating

④ 2ndary interactions nucleon, π ③ Inelasticity k= 1-plead/pbeam

If softer shallow development If harder deep penetrating

Important, but irrelevant to A.S. ⑤PT ⑥multiplicity (relevant to Nµ )

slide-5
SLIDE 5

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

5

Outline of this talk

What type of interactions we concern ? Sort out the data regarding as relevance to air

showers phenomena, especially focusing on Xmax

Inelastic cross section Forward energy spectra Inelasticity low energy data

Nuclear effect is important, but …

This talk focuses just on p-p Comments on possible p-A runs at LHC before long

shutdown

slide-6
SLIDE 6

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

6

Beam 1 Beam 2

The 7 LHC experiments

IP1 : ATLAS LHCf IP5 :CMS TOTEM IP2: ALICE IP8: LHCb, MoEDAL

slide-7
SLIDE 7

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

7

A LHC detector and pseudorapidity

ln(tan ) 2 θ η = −

LHC tunnel IP

ZDC (η>~8.5)

Central detector (ATLAS) (pseudo)rapidity

slide-8
SLIDE 8

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

8

pseudorapidity and interactions

σ @7TeV

φ η

  • 10
  • 10

φ η

  • 10
  • 10

φ η

  • 10
  • 10

φ η

  • 10
  • 10

η gap

Elastic Single diffractive Double diffractive Non- diffractive

~25mb ~10mb ~10mb ~50mb

slide-9
SLIDE 9

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

9

Very forward : Majority of energy flow (√s=14TeV)

8.4 < η < ∞ Multiplicity Energy Flux All particles neutral

Most of the energy flows into very forward

( Particles of XF > 0.1 contribute 50% of shower particles )

slide-10
SLIDE 10

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

10

Energy flow for √s=14 and 7TeV

14TeV 7TeV Double diffractive π0 Single diffractive π0

ATLAS/CMS CMS HF CMS HF LHCf/ZDC LHCf/ZDC

slide-11
SLIDE 11

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

11

Inelastic cross section

  • TOTEM
  • ATLAS
  • CMS
  • ALICE
slide-12
SLIDE 12

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

12

Measurement of σinel

2

16 1 1

=

+ + =

t el inel el tot

dt dR R R π ρ σ

= t el

dt dR

inel el tot tot

R R R L + = = σ

projectile Total elastic rate: Rel Total inelastic rate : Rinel : Nproj Total rate: Rtot : Ntarg Elastic rate at 0 degree target

arg arg

2

t proj t proj rev

N N f L Σ Σ = π

( ) ( ) ( ) ( )

Im Re

el el

f f = ρ

eff

  • bs

inel inel inel

R R L R ε σ = = ,

Simple way Optical theorem VdM scan

2

) (

  • ut

in

p p t − =

slide-13
SLIDE 13

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

13

ATLAS σinela measurement

Nature Commun. 2 (2011) 463

sel trig BG

  • bs

inel

f Ldt N N ε ε ξ σ

ξ

6

10 5 6

1 ) 10 5 (

× < −

− × × − = × >

Minimum Bias Trigger Scintillator(MBTS)

mb

inel

.) (exp 10 . 2 33 . 60 ) 10 5 (

6

± = × >

ξ σ mb extr s m p

inel

) ( 9 . 6 .) (exp 4 . 2 1 . 69 ) / (

2

± ± = > ξ σ

  • Use MB events with 20.3 +- 0.7µb-1
  • Cut diffractive events ( ξ < 5E-6)
  • extrapolate entire ξ
slide-14
SLIDE 14

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

14

TOTEM “Roman Pod” measurement

RP (147 m) RP (220m)

14 m

Roman Pot stations in the LHC tunnel

( F.Ferro, Diffraction 2010)

slide-15
SLIDE 15

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

15

Roman pod measurement

= t el

dt dσ

= 503.7 +- 1.5 +- 26.7 mb/GeV2

2 2 2

1 ) ( 16

=

+ =

t el tot

dt dR c ρ π σ

  • =

t el

dt dσ

ρ =0.14 +0.01-0.08 (COMPETE collaboration)

mb

tot 8 . 2 7 . 2

2 . 3 . 98

+ −

± = σ

mb dt dt d

el el

2 . 1 2 . 8 . 24 ± ± = = ∫ σ σ

Integrated over entire “t” region

mb

el tot inel 8 . 1 3 . 1

6 . 5 . 73

+ −

± = − = σ σ σ

EPL, 95 (2011) 41001

TOTEM σinel

slide-16
SLIDE 16

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

16

σinel result @ 7TeV

TOTEM 73.5+-0.6+1.8-1.3 mb dσ/dt(t=0) ATLAS 69.4+-2.4+-6,9 mb MBTS sample CMS 68.0+-2.0+-2.4+-4 mb Ntrk sample ALICE 72.7+-1.1+-5.1 mb VZERO sample

Tevatron UA4 ISR LHC

slide-17
SLIDE 17

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

17

ALICE σdiffraction

  • J. Phys. G: Nucl. Part. Phys. 38 (2011) 124044
slide-18
SLIDE 18

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

18

Forward energy spectra / Inelastiscity / PT

LHCf UA7 CMS FCAL RHIC BRAHMS Forward neutron spectra

slide-19
SLIDE 19

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

19

LHC zero degree experimental site

ATLAS

140m

LHCf/ZDC

TAN absorber 140m

96mm

Charged particles Charged particles (+) (+) Neutral particles Neutral particles Beam pipe Beam pipe Protons Protons Charged particles Charged particles ( (-

  • )

)

slide-20
SLIDE 20

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

20

LHCf: location and detector layout

INTERACTION POINT INTERACTION POINT IP1 (ATLAS) IP1 (ATLAS)

Detector II Detector II Tungsten Tungsten Scintillator Scintillator Silicon Silicon µ µstrips strips Detector I Detector I Tungsten Tungsten Scintillator Scintillator Scintillating Scintillating fibers fibers

Arm#1 Detector 20mmx20mm+40mmx40mm 4 SciFi tracking layers

44X0, 1.6 λint

Arm#2 Detector 25mmx25mm+32mmx32mm 4 Silicon strip tracking layers

140 m 140 m n π0 γ γ 8 cm 6 cm Front Counter Front Counter

slide-21
SLIDE 21

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

21

Rapidity vs Forward energy spectra

η η = 8 . 4 η=8.77 η=7.60 η η = 6 . 9 1 η=5.99 η=7.60 η η = 6 . 9 1 η=5.99 η η = 8 . 4 η=8.77

θ [μrad] 310

Projected edge

  • f beam pipe

Viewed from IP1 (red:Arm1, blue:Arm2)

slide-22
SLIDE 22

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

22

LHCf single γ spectra at 7TeV

DPMJET 3.04 QGSJETII-03 SIBYLL 2.1 EPOS 1.99 PYTHIA 8.145

Blue hatch: Statistics errors of MC Gray hatch : Systematic Errors

Phys.Lett. B703 (2011) 128-134

slide-23
SLIDE 23

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

23

New LHCf single γ spectra at 900 GeV

( to be submitted PLB)

slide-24
SLIDE 24

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

24

LHCf 900GeV single γ spectra: Data/MC

900GeV

slide-25
SLIDE 25

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

25

LHCf 900GeV single γ spectra: Data/MC

High η low η 900GeV 7TeV

slide-26
SLIDE 26

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

26

UA7 π0 PT at 630GeV

UA7 630GeV p p γ

Phys.Lett. B242 (1990) 531-535 6.5 4.5 2.5 1.5 Y

dσ/dY (mb) PT ( MeV/c)

slide-27
SLIDE 27

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

27

New LHCf π0 PT at 7TeV ( Preliminary )

slide-28
SLIDE 28

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

28

CMS HF : forward energy flow

ATLAS/CMS CMS HF LHCf/ZDC

slide-29
SLIDE 29

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

29

CMS HF: Forward energy flow

0.9TeV 7TeV

CERN-PH-EP/2011-086, arXiv/0329842

3.15< η <4.19

slide-30
SLIDE 30

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

30

RHIC BRAHMS : charged spectra at √s=200GeV

PRL 98, 252001 (2007)

slide-31
SLIDE 31

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

31

Inelasticity~ 0 degree neutron spectra

Important for Xmax and also Nµ Measurement of inelasticity at LHC energy

Neutral hadrons at 14 TeV (LHCf acceptance, no resolution) Neutral hadrons at 14 TeV (LHCf acceptance, 30% resolution)

slide-32
SLIDE 32

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

32

Future data for forward energy flow

Forward neutrons by LHC ZDC’s

So far working well for centrality in HI runs Potentially they can work nicely (PID ?) Combined LHCf+ATLAS ZDC may benefit

Other LHC forward detectors

CMS CASTOR :Only coverage for η ~ 6 TOTEM T1, T2, LHCb VELO( 1.6<η< 4.9 ? )

New LHC detectors ?

CMS Forward Shower Calorimeter (FSC) ? Roman Pod type calorimeter (a la UA7) ?

RHIC 0dgree measurement ?

√s= 500GeV with larger PT acceptance Possible π0 measurement

RHIC IP

slide-33
SLIDE 33

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

33

Summary : forward spectra coverage

10 100 1000 10000

√s η

4 2 6 8 10

LHCf 7TeV γ LHCf 0.9TeV γ LHCf 7TeV π0 CMS HF ATLAS CMS ALICE dn/dη BRAHMS PT UA7 π0 CDF dn/dη UA5 dn/dη ISR dn/dη θ =<PT>/ Pbeam <PT>~0.4GeV/c

slide-34
SLIDE 34

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

34

Low energy data

Exp beam target NA49 CERN SPS 128GeV/c p C, p NA61 CERN SPS 31GeV/c p, etc π+- C, Be… HARP CERN PS 3,5,8,9,12GeV/c p π+- C,Be,p.. MIPP FNAL-MI 58,120GeV/c p C,Be… NA61 pbeam = 31GeV/c HARP p+Cπ +X pbeam = 12GeV/c Relevant E (1015eV shower) Plab=10~103GeV

slide-35
SLIDE 35

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

35

Available accelerator data summary

slide-36
SLIDE 36

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

36

Available accelerator data summary

HARP NA61 NA49 MIPP

slide-37
SLIDE 37

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

37

Summary

Synergy btw UHECR and LHC is so important to solve

UHECR problems and to explore ultra high E interactions at beyond-LHC.

Key parameters for understanding air showers, σinel, forward

spectra, inelasticity, multiplicity, PT should be measured in various energy ranges.

First TOTEM σinel LHCf forward spectra Recent progress in various energy range (i.e. NA61, HARP, etc.)

Striking impacts on UHECR analysis has been given by

recent LHC data as well as legacy data by various accelerator experiments. Stay tuned.

Nuclear effects (QGP, shadowing, etc.) not address here are

also important. LHC A-A, p-A run will be able to address it.

slide-38
SLIDE 38

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

38

p+C π X data

NA49 ( pbeam=158GeV/c ) NA61 ( pbeam=31GeV/c ) NA61(SHINE) pbeam data sets NA49 pbeam = 158GeV/c NA61 pbeam = 31GeV/c

slide-39
SLIDE 39

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

39

HARP / MIPP

MIPP detector

slide-40
SLIDE 40

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

40

Available accelerator data summary

slide-41
SLIDE 41

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

41

slide-42
SLIDE 42

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

42

Parent particles relevant for LHCf

  • bservations

Sybill at 7TeV

slide-43
SLIDE 43

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

43

slide-44
SLIDE 44

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

44

Inelasticity

slide-45
SLIDE 45

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

45

LHCf forward spectra: Data/MC

High η low η

slide-46
SLIDE 46

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

46

LHCf forward spectra: Data/MC

High η low η

slide-47
SLIDE 47

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

47

Elasticity ( XF of leading baryon)

slide-48
SLIDE 48

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

48

p+p p+X at Pbeam= 205GeV (FNAL bubble chamber)

FNAL bubble chamber p + p p + X (at 205GeV/c) Whitmore et al, PRD11(1975)3124

slide-49
SLIDE 49

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

49

UA7 vs. PYTHIA8

PYTHIA8 does not reproduce UA7 Can we confirm/update/improve UA7?

Pare et al. PYTHIA8 (histos) vs UA7 fit

49

slide-50
SLIDE 50

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

50

LHCf forward spectra: Data/MC

High η low η

slide-51
SLIDE 51

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

51

Setup in IP1-TAN (side view)

LHCf Front Counter LHCf Calorimeter BRAN-IC ZDC type1 ZDC type2

Beam pipe

TAN

Neutral particles

Side view

BRAN-Sci

IP1

slide-52
SLIDE 52

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

52

Forward E spectra forseen at 14TeV (MC for ~0.1nb-1) π0

single γ

Neutrons (w/ 30% resolution) Neutrons (true energy)

slide-53
SLIDE 53

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

53

The single photon energy spectra at 0 degree

DATA

15 May 2010 17:45‐21:23, at Low Luminosity 6x1028cm‐2s‐1, no beam

crossing angle

0.68 nb‐1 for Arm1, 0.53nb‐1 for Arm2

MC

DPMJET3.04

, QGSJETII03, SYBILL2.1, EPOS1.99 PYTHIA 8.145 with the default parameters.

107 inelastic p-p collisions by each model.

Analysis

Two pseudo-rapidity, η>10.94 and 8.81<η<8.9. No correction for geometrical acceptance. Combine spectra between Arm1 and Arm2.

Normalized by number of inelastic collisions

with assumption as σ inela = 71.5mb. (c.f. 73.5±0.6. mb by TOTEM )

+1.8

  • 1.3

(O.Adriani et al., PLB703 (2011) 128-134 Arm1 Arm2

slide-54
SLIDE 54

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

54

Particle Identification

Event selection and correction

– Select events <L90% threshold and multiply P/ε

ε (photon detection efficiency) and P (photon purity)

– By normalizing MC template L90% to data,

ε and P for certain L90% threshold are determined.

dE Integral of dE

Photon Hadron

Calorimeter layers Calorimeter layers

Elemag: 44r.l. Hadronic: 1.7λ Calorimeter Depth L90% Distribution

slide-55
SLIDE 55

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

55

Very forward – connection to low-x physics

Very forward region : collision of a low-x parton with a large-x

parton

Small-x gluon become dominating in higher energy collision by

self interaction.

But they may be saturated (Color Glass Condenstation)

Low-x high-x Very forward

Naively CGC-like suppression may

  • ccur in very forward at high energy

However situation is more complex (not simple hard parton collsions, but including soft + semi-hard ) soft semi- hard hard

slide-56
SLIDE 56

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

56

What PT range LHCf sees ? γ π0

pp 7TeV, EPOS

slide-57
SLIDE 57

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

57

RHIC:π0 at √s =500GeV

ZDC space at PHENIX (by Goto‐san): 10cm radius beam pipe aperture at 18m => η>5.9

57

slide-58
SLIDE 58

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

58

Pi0 in 500GeV p-p collisions by PYTHIA8

Eta vs. Energy Eta vs. number flux Eta vs. energy flux Vertical lines at |η|=6

Multiplicity Energy Flux All particles neutral

slide-59
SLIDE 59

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

59

Comparison between the two detector

Pseudo-rapidity selection, η>10.94 and 8.81<η<8.9 Normalized by number of inelastic collisions

with assumption as σ inela = 71.5mb ( <->73.5±0.6. mb by TOTEM )

Spectra in the two detectors are

consistent within errors.

  • Combined between spectra of Arm1 and Arm2

by weighted average according to errors Arm1 detector Arm2 detector

+1.8

  • 1.3
slide-60
SLIDE 60

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

60

Forward production spectra vs Shower curve

XF = E/Etot

Half of shower particles comes from large XF γ

Measurement at very forward region is needed

slide-61
SLIDE 61

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

61

Transition curve for 10 17eV proton

π0 reconstructed spectrum

O N C R

⇔ΔXmax (DPMJETⅢ-QGSJETⅡ) = 36g/ cm2

Xmax = 696.6g/ cm2 ± 2.3 +15.9, -7.8 g/ cm2

statistical systematic

P H Y S

OSCAR ADRIANI LHCC MEETING, CERN 23 MARCH 2011

slide-62
SLIDE 62

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

62

slide-63
SLIDE 63

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

63

Contribution from very forward production

5×1019 eV proton showers ( 60 deg zenith) # of electrons

No cut low XF γ origin ( xF < 0.05 ) π,Κ origin ( xF < 0.1 ) Half of shower particles comes from large XF γ

Measurement at very forward region is needed

slide-64
SLIDE 64

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

64

Forward energy spectra

η η = 8 . 4 η=8.77 η=7.60 η η = 6 . 9 1 η=5.99

slide-65
SLIDE 65

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

65

Impact of parameters of interactions

R.Ulrich et al., PRD83(2011)054026

Xmax RMS Xmax Xmax RMS Xmax p p Fe Fe

Cross section multiplicity elasticity

slide-66
SLIDE 66

Y.Itow, Review of Accelerator data UHECR2012@ 14Feb 2012

66

Big LHC detectors