Antiproton production in p-He collisions, and more, at LHCb LHCb on - - PowerPoint PPT Presentation

antiproton production in p he collisions and more at lhcb
SMART_READER_LITE
LIVE PREVIEW

Antiproton production in p-He collisions, and more, at LHCb LHCb on - - PowerPoint PPT Presentation

Antiproton production in p-He collisions, and more, at LHCb LHCb on a Space Mission He at 6.5 TeV rest proton antiproton Giacomo Graziani (INFN Firenze) on behalf of the LHCb Collaboration ICRC 2017, Busan, Korea July 15, 2017 The LHCb


slide-1
SLIDE 1

Antiproton production in p-He collisions, and more, at LHCb

LHCb on a Space Mission

6.5 TeV proton He at rest antiproton

Giacomo Graziani (INFN Firenze)

  • n behalf of the LHCb Collaboration

ICRC 2017, Busan, Korea July 15, 2017

slide-2
SLIDE 2

The LHCb Experiment

LHCb is the experiment devoted to heavy flavours at the LHC Focused on CP violation and rare signatures in b and c decays Exploiting LHC as the biggest b and c fac- tory on earth Detector requirements: Forward geometry optimize acceptance for bb pairs Tracking : best possible proper time and momentum resolution Particle ID : excellent capabilities to select exclusive decays Trigger : high flexibility and bandwidth (up to 15 kHz to disk) ➨ allowed to widen our physics program to include hadron spectroscopy, EW physics, kaon physics, heavy ion physics (pPb and PbPb collisions) ...

  • G. Graziani

slide 2 ICRC 2017

slide-3
SLIDE 3

SMOG: the LHCb internal gas target

LHCb is the LHC experiment with “fixed-target like” geometry very well suited for...fixed target physics!

JINST 3, (2008) S08005 Int.J.Mod.Phys.A30 (2015) 1530022

The System for Measuring Overlap with Gas (SMOG) allows to inject small amount of noble gas (He, Ne, Ar, ...) inside the LHC beam around (∼ ±20 m) the LHCb collision region Expected pressure ∼ 2 × 10−7 mbar Originally conceived for the luminosity determination with beam gas imaging

JINST 9, (2014) P12005

Became the LHCb internal gas target for a rich and var- ied fixed target physics program

  • G. Graziani

slide 3 ICRC 2017

slide-4
SLIDE 4

Soft QCD for Cosmic Rays Physics

Fixed target collisions allow to study exclusive particle production at the energy scale of ∼ 100 GeV, with access to large x in the target; can provide valuable inputs for modeling cosmic ray showers in the atmosphere and in the cosmos, in particular for antimatter production AMS02 results provide unprecedented accuracy for measurement of p/p ratio in cosmic rays at high energies

PRL 117, 091103 (2016)

hint for a possible excess, and milder en- ergy dependence than expected prediction for p/p ratio from spallation

  • f primary cosmic rays on intestellar

medium (H and He) is presently limited by uncertainties on p production cross- sections, particularly for p-He no previous measurement of p production in p-He, current predictions vary within a factor 2 the LHC energy scale and LHCb +SMOG are very well suited to perform this mea- surement

Giesen et al., JCAP 1509, 023 (2015)

  • G. Graziani

slide 4 ICRC 2017

slide-5
SLIDE 5

Detector and Acceptance

JINST 3, (2008) S08005 Int.J.Mod.Phys.A30 (2015) 1530022

p [GeV/c]

20 40 60 80 100

[GeV/c]

T

p

0.5 1 1.5 2 2.5 3 3.5 4 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 =5 η = 4 . 5 η

LHCb Preliminary

LHCb-CONF-2017-002

Total acceptance × reconstruction efficiency for antiprotons Tracking efficiency estimated from simulation, validated on (pp) data

  • G. Graziani

slide 5 ICRC 2017

slide-6
SLIDE 6

The p-He run

LHCb-CONF-2017-002

Data collected in May 2016, with proton en- ergy 6.5 TeV, √sNN = 110 GeV Most data from a single LHC fill (5 hours) Minimum bias trigger, fully efficient on can- didate events Exploit excellent particle identification (PID) capabilities in LHCb to count antiprotons in (p, pT) bins within the kinematic range 12 < p < 110 GeV /c, pT > 0.4 GeV /c Exploit excellent vertexing capabilities to separate prompt and detached components. Only the prompt component included in this preliminary result (analysis of component from hyperon decays ongoing). Residual detached component estimated to be (2.6 ± 0.6)% and subtracted Background from gas contamination measured to be 0.6 ± 0.2%

Data p 21.4 - 24.4 pt 1.2- 1.5

DLL (p -K)

  • 200
  • 100

100 200 ) π DLL (p -

  • 200
  • 150
  • 100
  • 50

50 100 150 200

LHCb Preliminary

  • 200
  • 100

100 200

  • 200
  • 150
  • 100
  • 50

50 100 150 200 20 40 60 80 100 120 140

Template for p

  • 200
  • 100

100 200

  • 200
  • 150
  • 100
  • 50

50 100 150 200 100 200 300 400 500 600 700 800 900

π Template for

  • 200
  • 100

100 200

  • 200
  • 150
  • 100
  • 50

50 100 150 200 50 100 150 200 250

Template for K

  • 200
  • 100

100 200

  • 200
  • 150
  • 100
  • 50

50 100 150 200 10 20 30 40 50 60

Template for ghost

  • G. Graziani

slide 6 ICRC 2017

slide-7
SLIDE 7

Normalization

Gas target density not precisely known, using p-e− elastic scattering Pro: LHCb sees the purely elastic regime: θ > 10mrad ➨ ϑs < 29 mrad, Q2 < 0.01 GeV2 ➨ cross-section very well known distinct signature with single low-p and very low pT electron track, and nothing else background events mostly expected from very soft collisions, where candidate comes from γ conversion or pion from central exclusive pro- duction event ➨ background expected to be charge symmetric, can use “single positrons” to model it in data

SPD hits

5 10 15 20

scattered electron candidates

1000 2000 3000 4000 5000 6000 7000 8000 9000

  • Simulation of single e

candidates

  • e

candidates

+

e

LHCb Preliminary

LHCb-CONF-2017-002

Cons: cross-section is small (order 100 µb, 3 orders of magnitude below hadronic cross section) electron has very low momentum and showers through beam pipe/detectors ➨ low acceptance and reconstruction efficiency

  • G. Graziani

slide 7 ICRC 2017

slide-8
SLIDE 8

Event display of a candidate scattered electron

  • G. Graziani

slide 8 ICRC 2017

slide-9
SLIDE 9

Electron spectra

LHCb-CONF-2017-002 p [MeV/c]

5000 10000 15000

Candidates per 260 MeV/c

1000 2000 3000 4000 5000

candidates

  • e

candidates

+

e

LHCb Preliminary

p [MeV/c]

5000 10000 15000

Candidates per 260 MeV/c

500 1000 1500 2000 2500 3000 candidates (Bkg Sub.)

  • e

Simulation (normalized) LHCb Preliminary

[MeV/c]

T

p

50 100

Candidates per 2.4 MeV/c

500 1000 1500 2000 2500 LHCb Preliminary

[MeV/c]

T

p

50 100

Candidates per 2.4 MeV/c

200 400 600 800 1000 1200 1400 1600 1800 LHCb Preliminary

Very good agreement with simulation of single scat- tered electrons Data confirm charge sym- metry of background L = 0.443 ± 0.011 ± 0.027 nb−1 Systematic from variation of selection cuts, largest dependence is on azimuthal angle equivalent gas pressure is 2.4 × 10−7 mbar, in agreement with the expected level in SMOG

  • G. Graziani

slide 9 ICRC 2017

slide-10
SLIDE 10

Total relative uncertainty per bin, in per cent

LHCb-CONF-2017-002

24.7 13.3 27.6 19.1 9.5 9.6 17.5 8.2 8.8 7.9 19.9 22.2 17.0 12.4 9.8 8.5 9.4 9.5 15.0 9.7 9.4 10.2 10.4 15.7 10.2 10.9 10.8 10.7 8.1 9.3 7.9 8.5 8.6 8.0 13.6 23.3 10.7 12.2 9.5 9.8 13.9 14.7 9.6 9.8 8.2 8.5 7.9 8.0 8.1 15.1 21.0 10.4 9.4 11.1 9.7 10.6 17.3 13.7 10.8 11.7 11.4 9.8 9.9 9.2 11.4 11.1 12.1 9.6 9.6 10.4 9.7 13.4 12.4 11.0 10.1 8.7 8.5 10.3 9.8 15.8 9.8 9.2 9.2 9.2 9.5 12.8 11.6 10.7 12.6 17.2 17.5 11.6 12.9 10.7 9.1 12.5 10.3 9.4 10.1 9.4 9.4 9.4 9.5 9.5 10.2 14.0 20.0 25.7 20.6 17.0 16.9 14.5 13.0 12.0 15.8 11.1 10.8 10.8 10.8 11.3 11.2 10.5 12.2 12.3 12.5 25.1 25.8 24.9 24.4 20.9 20.6 20.1 19.6 20.4 22.2 28.6

p [GeV/c]

2

10

[GeV/c]

T

p

1 5 10 15 20 25 = 5 η = 4 . 5 η

LHCb Preliminary

dominated by systematics largest correlated uncertainty is the 6% from normalization largest uncorrelated uncertainty from PID analysis

  • G. Graziani

slide 10 ICRC 2017

slide-11
SLIDE 11

Result for cross section, compared with EPOS LHC

LHCb-CONF-2017-002

[GeV/c]

T

p

1 2 3 4

]

2

/GeV

2

b c µ [

T

X)/dpdp p ( σ

2

d

19 −

10

17 −

10

15 −

10

13 −

10

11 −

10

9 −

10

7 −

10

5 −

10

3 −

10

1 −

10 10

2

10

x (12.0 < p < 14.0 GeV/c)

10 x (14.0 < p < 16.2 GeV/c)

  • 1

10 x (16.2 < p < 18.7 GeV/c)

  • 2

10 x (18.7 < p < 21.4 GeV/c)

  • 3

10 x (21.4 < p < 24.4 GeV/c)

  • 4

10 x (24.4 < p < 27.7 GeV/c)

  • 5

10 x (27.7 < p < 31.4 GeV/c)

  • 6

10 x (31.4 < p < 35.5 GeV/c)

  • 7

10 x (35.5 < p < 40.0 GeV/c)

  • 8

10 x (40.0 < p < 45.0 GeV/c)

  • 9

10 x (45.0 < p < 50.5 GeV/c)

  • 10

10 x (50.5 < p < 56.7 GeV/c)

  • 11

10 x (56.7 < p < 63.5 GeV/c)

  • 12

10 x (63.5 < p < 71.0 GeV/c)

  • 13

10 x (71.0 < p < 79.3 GeV/c)

  • 14

10 x (79.3 < p < 88.5 GeV/c)

  • 15

10 x (88.5 < p < 98.7 GeV/c)

  • 16

10 x (98.7 < p < 110.0 GeV/c)

  • 17

10

LHCb Preliminary

Result for prompt production (excluding weak decays of hy- perons) The total inelastic cross section is also measured to be σLHCb

inel

= (140 ± 10) mb The EPOS LHC prediction

[T. Pierog at al, Phys. Rev. C92 (2015), 034906]

is 118 mb, ratio is 1.19 ± 0.08.

  • G. Graziani

slide 11 ICRC 2017

slide-12
SLIDE 12

Result for cross section, ratio with models

LHCb-CONF-2017-002

DATA / PREDICTION

1 2 3 4 data/ prediction 0.5 1 1.5 2 2.5 3

12.0 < p < 14.0 GeV/c

LHCb Preliminary 1 2 3 4 0.5 1 1.5 2 2.5 3

14.0 < p < 16.2 GeV/c

LHCb Preliminary 1 2 3 4 0.5 1 1.5 2 2.5 3

16.2 < p < 18.7 GeV/c

LHCb Preliminary 1 2 3 4

data/ prediction 0.5 1 1.5 2 2.5 3

18.7 < p < 21.4 GeV/c

LHCb Preliminary 1 2 3 4 0.5 1 1.5 2 2.5 3

21.4 < p < 24.4 GeV/c

LHCb Preliminary 1 2 3 4 0.5 1 1.5 2 2.5 3

24.4 < p < 27.7 GeV/c

LHCb Preliminary 1 2 3 4

data/ prediction 0.5 1 1.5 2 2.5 3

27.7 < p < 31.4 GeV/c

LHCb Preliminary 1 2 3 4 0.5 1 1.5 2 2.5 3

31.4 < p < 35.5 GeV/c

LHCb Preliminary 1 2 3 4 0.5 1 1.5 2 2.5 3

35.5 < p < 40.0 GeV/c

LHCb Preliminary 1 2 3 4

data/ prediction 0.5 1 1.5 2 2.5 3

40.0 < p < 45.0 GeV/c

LHCb Preliminary 1 2 3 4 0.5 1 1.5 2 2.5 3

45.0 < p < 50.5 GeV/c

LHCb Preliminary 1 2 3 4 0.5 1 1.5 2 2.5 3

50.5 < p < 56.7 GeV/c

LHCb Preliminary 1 2 3 4

data/ prediction 0.5 1 1.5 2 2.5 3

56.7 < p < 63.5 GeV/c

LHCb Preliminary 1 2 3 4 0.5 1 1.5 2 2.5 3

63.5 < p < 71.0 GeV/c

LHCb Preliminary 1 2 3 4 0.5 1 1.5 2 2.5 3

71.0 < p < 79.3 GeV/c

LHCb Preliminary [GeV/c]

T

p 1 2 3 4 data/ prediction 0.5 1 1.5 2 2.5 3

79.3 < p < 88.5 GeV/c

LHCb Preliminary [GeV/c]

T

p 1 2 3 4 0.5 1 1.5 2 2.5 3

88.5 < p < 98.7 GeV/c

LHCb Preliminary [GeV/c]

T

p 1 2 3 4 0.5 1 1.5 2 2.5 3

98.7 < p < 110.0 GeV/c

LHCb Preliminary

Trasverse Momentum (GeV/c)

EPOS LHC EPOS 1.99 QGSJETII-04 HIJING 1.38

Cross section is larger by factor ∼ 1.5 wrt EPOS LHC (mostly from larger p rate per collision). Better agreement with EPOS 1.99, HIJING 1.38 and QGSJET-IIm (low energy extension

  • f QGSJET-II-04, not shown)

Many thanks to T. Pierog for his advice with EPOS/CRMC!

  • G. Graziani

slide 12 ICRC 2017

slide-13
SLIDE 13

Prospects for antimatter production

We plan to extend the study to p produced by hyperon decays (accounting for 20-30%

  • f the production). LHCb can cleanly select

decays of Λ

]

2

c Invariant Mass [MeV/

+

π p 1100 1110 1120 1130

2

c Candidates / 0.6 MeV/ 100 200 300 ]

2

c Invariant Mass [MeV/

+

π p 1100 1110 1120 1130

2

c Candidates / 0.6 MeV/ 100 200 300

2

c 0.03 MeV/ ± = 1115.75 µ

2

c 0.33 MeV/ ± = 1.23 σ 36 ± N = 1177

LHCb

= 0.9 TeV s

Λ →pπ+ 0.25 < pT < 2.50 GeV /c 2.5 < y < 3.0

JHEP 1108 (2011) 034

Another p-He run was performed in november 2016 with a 4 TeV beam (√sNN =87 GeV) ➨ scaling violation can be constrained production of pions and kaons is also being measured ➨ positron production. Ratios of particle species, not affected by uncertainty on luminosity, can provide precise constraints to soft QCD models investigating our potential for antinuclei d, t and 3He RICH can actively identify d from p 36 GeV /c and t,3He from 54 GeV /c. dE/dx and time-of-flight information from tracking detectors at low momentum

  • G. Graziani

slide 13 ICRC 2017

slide-14
SLIDE 14

And charm, of course

Exclusive production of charm states is a specialty of LHCb

]

2

c ) [MeV/

+

π

  • m(K

1800 1820 1840 1860 1880 1900 1920 1940

)

2

c Candidates / (5 MeV/

100 200 300 400 500 600 700

= 87 GeV pHe

NN

s

LHCb preliminary

D0 → K−π+

]

2

c ) [MeV/

  • µ

+

µ m(

2950 3000 3050 3100 3150 3200

)

2

c Candidates / (10 MeV/

20 40 60 80 100

= 87 GeV pHe

NN

s

LHCb preliminary

J /ψ → µ+µ− In fixed target mode, access the large-x (target fragmenta- tion) region, where charm PDF is affected by antishadow- ing and possibly intrinsic charm (IC) effects. Important for high-energy neutrino astrophysics: back- ground for the ICECUBE experiment is dominated by charm production in atmospheric showers

  • E

2 [GeV cm-2 s-1 sr-1]

Eµ [GeV]

10-10 10-9 10-8 10-7 10-6 103 104 105 106

I n t r i n s i c C h a r m dotted grey: Conv. Atm. µ solid grey: Conv. Atm. µ + BERSS solid black: Conv. Atm. µ + Intrinsic Charm (H3A) + BERSS IceCube astrophysical flux IceCube µ

Laha and Brodsky, arXiv:1607.08240

  • G. Graziani

slide 14 ICRC 2017

slide-15
SLIDE 15

Charm in p-Ar collisions @ 110 GeV

LHCb-CONF-2017-001

x-Feynman distribution for D0 and J /ψ, compared to Pythia8 prediction

F

Feynman-x x

0.25 − 0.2 − 0.15 − 0.1 − 0.05 −

F

/dx

D

dN

0.05 0.1 0.15 0.2 0.25

6

10 ×

= 110 GeV pAr

NN

s

LHCb preliminary

sinh(y*)

NN

s

D

M = 2

F

x

F

Feynman-x x

0.45 − 0.4 − 0.35 − 0.3 − 0.25 − 0.2 − 0.15 − 0.1 − 0.05 −

F

/dx

ψ J/

dN

0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2

3

10 ×

= 110 GeV pAr

NN

s

LHCb preliminary

sinh(y*)

NN

s

ψ J/

M = 2

F

x

Obtained from the first small (few nb−1) p-Ar data sample acquired in 2015 ∼ 6500 D0 and 500 J /ψ Result limited by statistics, but demonstrates the potential for unique measurements Differential shapes can already constrain models with IC

  • G. Graziani

slide 15 ICRC 2017

slide-16
SLIDE 16

Conclusions

LHCb started its fixed target program becoming an unexpected contributor to cosmic ray physics! The p production measurement in p-He collisions is expected to narrow down significantly the uncertainty on the p/p prediction for cosmic rays Nice interplay with the cosmic rays community

thanks to O. Adriani, L. Bonechi, F. Donato and A. Tricomi for proposing this measurement

the LHCb space mission just started!

  • G. Graziani

slide 16 ICRC 2017

slide-17
SLIDE 17

Additional Material

  • G. Graziani

slide 17 ICRC 2017

slide-18
SLIDE 18

RICH Performance

  • Eur. Phys. J. C 73 (2013) 2431

Particle separation in RICH1 K/p separation vs momentum

  • G. Graziani

slide 18 ICRC 2017

slide-19
SLIDE 19

Selection of prompt component: the Vertex Detector

Current analysis limited to “prompt” component (direct production and p from strong resonance decays) Can be distinguished from p produced by weak decays of hyperons and secondary interactions using the excellent LHCb vertexing capabilities

JINST 9 (2014) P09007

  • G. Graziani

slide 19 ICRC 2017

slide-20
SLIDE 20

Background from weak hyperon decays

detached component from weak decays of hyperons is treated as a background suppressed by requir- ing small impact pa- rameter (IP)

]

  • 1

c [GeV

T

p 1/ 0.5 1 1.5 2 2.5 3 m] µ resolution [

x

IP 10 20 30 40 50 60 70 80 90 100

T

p = 11.6 + 23.4/ σ 2012 data,

T

p = 11.6 + 22.6/ σ Simulation,

LHCb VELO

JINST 9 (2014) P09007

Residual detached component estimated to be (2.6 ± 0.6)% and subtracted Systematic uncertainty estimated from data/MC comparison of IP tails

  • G. Graziani

slide 20 ICRC 2017

slide-21
SLIDE 21

Background from Residual Vacuum

Residual vacuum in LHC is not so small (∼ 10−9 mbar ) compared to SMOG pressure Can be a concern, especially for heavy contaminants (larger cross section than He), and beam-induced local outgassing Direct measurement in data: about 15% of delivered protons on target acquired before He injection (but with identical vacuum pumping configuraton)

PV Track Multiplicity

5 10 15 20 25 30 35 40

fraction of candidates

0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 p on He gas p on Residual Vacuum

LHCb Preliminary

LHCb-CONF-2017-002

Gas impurity found to be small: 0.6 ± 0.2% PV multiplicity in residual vacuum events is lower than in He events, but has longer tails ➨ confirm findings from Rest Gas Analysis that resid- ual vacuum is mostly H2, with small heavy contaminants

  • G. Graziani

slide 21 ICRC 2017

slide-22
SLIDE 22

Result for cross section: final uncertainties (relative)

LHCb-CONF-2017-002

Statistical: Yields in data/PID calibration 0.7 − 10.8% (< 3% for most bins) Normalization 2.5% Correlated Systematic: Normalization 6.0% GEC and PV cut 0.3% PV reco 0.8% Tracking 2.2% Residual Vacuum Background 0.1% Non-prompt background 0.3 − 0.7% PID 1.2 − 5.0% Uncorrelated Systematic: Tracking 3.2% IP cut efficiency 1.0% PID 0 − 26% (< 10% for most bins) MC statistics 0.8 − 15% (< 4% for pT < 2 GeV/c)

  • G. Graziani

slide 22 ICRC 2017

slide-23
SLIDE 23

Comparison with EPOS LHC, vs x-Feynman

F

x

0.2 − 0.1 −

]

2

/GeV

2

b c µ [

T

/dpdp σ

2

d

1 −

10 1 10

2

10

12<p<14 GeV/c 14<p<16.2 GeV/c 16.2<p<18.7 GeV/c 18.7<p<21.4 GeV/c 21.4<p<24.4 GeV/c 24.4<p<27.7 GeV/c 27.7<p<31.4 GeV/c 31.4<p<35.5 GeV/c 35.5<p<40 GeV/c 40<p<45 GeV/c 45<p<50.5 GeV/c 50.5<p<56.7 GeV/c 56.7<p<63.5 GeV/c 63.5<p<71 GeV/c 71<p<79.3 GeV/c 79.3<p<88.5 GeV/c 88.5<p<98.7 GeV/c 98.7<p<110 GeV/c EPOSLHC

LHCb Preliminary

  • G. Graziani

slide 23 ICRC 2017

slide-24
SLIDE 24

Comparison with models

EPOS LHC PRC92, 034906 (2015)

p [GeV/c] 20 40 60 80 100 [GeV/c]

T

p

0.5 1 1.5 2 2.5 3 3.5 4 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8

LHCb Preliminary

EPOS1.99

σ /

LHCb

σ

QGSJET-II-04 PRD83, 014018 (2011)

p [GeV/c] 20 40 60 80 100 [GeV/c]

T

p

0.5 1 1.5 2 2.5 3 3.5 4 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8

LHCb Preliminary

QGSJETII

σ /

LHCb

σ

HIJING 1.38 Comp. Phys. Comm. 83, 307 (1994)

p [GeV/c] 20 40 60 80 100 [GeV/c]

T

p

0.5 1 1.5 2 2.5 3 3.5 4 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8

LHCb Preliminary

HIJING

σ /

LHCb

σ

QGSJET-IIm Astr. J. 803:54 (2015)

p [GeV/c] 20 40 60 80 100 [GeV/c]

T

p

0.5 1 1.5 2 2.5 3 3.5 4 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8

LHCb Preliminary

QGSJETIIm

σ /

LHCb

σ

  • G. Graziani

slide 24 ICRC 2017