LANL contribution plan to the EIC Jets and Heavy Flavor working group - Xuan Li on behalf of the Los Alamos National Laboratory

LA-UR-20-21862

LANL contribution plan to the EIC Jets
 and Heavy Flavor working group

 Xuan Li on behalf of the Los Alamos
 National Laboratory

 Xuan Li (LANL) 1

LANL Electron-Ion Collider project

• We have started our
 EIC project supported
 by the Laboratory
 Directed funding by
 the Los Alamos
 National Laboratory
 LDRD office.
• PI: Ivan Vitev, Co-PI:
 Xuan Li
• Our team: 15+ staffs
 and postdocs working
 on this EIC project.
 New postdocs, visiting students,
 postdocs and scientists.
 Xuan Li (LANL) 2

LANL heavy flavor experimental studies
• The LANL team plays a leading role in the heavy flavor
 data analysis and simulation studies at PHENIX, LHCb
 and sPHENIX in the past 20 years.
 LHCb heavy flavor exotic studies.
 arXiv: 2002.01551, LHCB-CONF-2019-005
 u 

 Entries/(1 MeV/c2)
 LHCb Preliminary

 BR(c (3872)® J/y p+ p-)
 10000 LHCb Preliminary

 BR(y(2S)® J/y p p-)
 pp s = 8 TeV

 " " 
 0.12 pp s = 8 TeV

 MeV/c22)
 9000 3800
 10000 LHCb Preliminary

 Entries/(1 MeV/c
 pp s = 8 TeV

 Prompt

 +
 Entries/(1 MeV/c2 )
 9000 3800

 3600
 8000 8000
 3600

 0.1

 Entries/(1
 3400

 7000

 b decays
 3400 3200

 7000 6000
 3200
 3000

 2800
 5000 3840 3850 3860 3870 3880 3890 3900 3910

 MJ/ ψ π +π - (MeV/c2 )
 0.08
 6000 4000
 3000
 3000
 2800
 2000
 5000 3840
 1000 3850 3860 3870 3880 3890 3900 3910 0.06
 3650 3700 3750 3800 2
 MJ/ ψ π +3850 3900
 π - (MeV/c )
 3950

 4000 MJ/ ψπ+π- (MeV/c2)

 c1
 0.04
 3000

 2000 0.02 pT > 5 GeV/c

 (3872)
 sy(2S)
 1000
 3650 3700 3750 3800 3850 3900 3950 0
 0 20 40 60 80 100 120 140 160 180 200

 c1
 MJ/ ψπ+π- (MeV/c2)

 sc
 NVELO
 tracks

 sPHENIX MVTX
 bottom hadron
PHENIX FVTX quarkonia and jet studies in
and open heavy flavor p+p and A+A
data analysis in p+p, collisions.
p+A and A+A collisions.
 Xuan Li (LANL) 3

The LANL EIC LDRD project aligns well with the
 EICUG timeline

 Updated
 timeline from
 the EICUG

 Full Tracker EIC
 LANL construction operation
 EIC LDRD by DOE
 NAS EIC EIC EIC EIC
 review CD0 CD1 CD2/3 construction
 2018 2020 2021 2023 2028
 Xuan Li (LANL) 4

LANL EIC heavy flavor and jet study goals (I)
• Through measuring heavy flavor hadrons, jets which
 can be treated as surrogates of initial quarks/gluons
 and their correlations in the hadron/nuclei going
 (forward) direction at the EIC. arXiv:1610.08536
 Phys. Rev. D 96, 114005 (2017)
 • To precisely determine the
 ' + p→ ' + ( ± ) + 
 initial quark/gluon
 e-
 distribution functions in
 e- the poorly constrained
arXiv:2002.05880 kinematic region.
 
 • To precisely study the
 quark/gluon
 fragmentation/hadronizati
 on processes.
 • To provide further
 ±
 information on the gluon
 Sivers function and other
 Xuan Li (LANL) spin observables. 5

LANL EIC heavy flavor and jet study goals (II)
• Through measuring heavy flavor hadrons, jets which
 can be treated as surrogates of initial quarks/gluons
 and their correlations in the hadron/nuclei going
 (forward) direction at the EIC.
arXiv:2002.05880 ' + Au→ ' + ( ±) + • To understand the
nPDF modification
 nuclear medium effects
 e- on hadron production
 e-
 such as modification on
 nuclear PDFs, parton
 energy loss mechanisms
 Parton energy loss and hadronization
 processes through the
 comparison of
 ± measured hadron/jet
 cross section between
 e+p and e+A collisions.
 Xuan Li (LANL) 6

LANL EIC heavy flavor simulation studies (I)
• Study the heavy flavor hadron and jet reconstruction with
 the proposed Forward Silicon Tracker for the EIC.
• The full analysis framework which includes the event
 generation, detector response in fast simulation,
 background embedding, and hadron reconstruction has
 been setup.
 D-meson decay channels
 • One triggered event (trigger rate at
 500kHZ) is embedded with on
 average 0.02 p+p (~12kHZ)
 background events. Need to add
 other background such as
 synchrotron radiation later.
 • Single track efficiency set at 95%.
 • Minimum Q2 = 1GeV/c2.
 • 100% PID.
 Xuan Li (LANL) 7

Reconstructed D mesons in PYTHIA8 simulation
• Mass distributions of clusters with track transverse decay
 length matching between charged tracks. Clusters are
 required to have at least one K± tracks. The performances are
 based on K/ separation but no charge separation.
• Detector options: pixel pitch 30 m, materials per detector
 layer: 0.4%X0 and the readout rate is at 500 kHZ.
 Reconstructed cluster mass with K± Reconstructed cluster mass with K±
 5000 5000
 w/ LANL FST w/ LANL FST
 ±
 D SIG/BKG: 1.510 arXiv:2002.05880 MC
 4000 0 0 MC 4000
 D (D ) SIG/BKG: 1.360 ±
 Rec. Ds
 Total fit
 D±s (D±s) SIG/BKG: 0.670
 0 0
 3000 Rec. D (D ) 3000
 ±
 Rec. D
 Integrated luminosity: 10 fb-1
 2000 2000

 1000 1000

 1.83 1.84 1.85 1.86 1.87 1.88 1.89 1.9 1.91 1.94 1.95 1.96 1.97 1.98 1.99 2
 2
 GeV/c GeV/c2
 To extract the signal/background ratio for reconstructed D mesons.
 Xuan Li (LANL) 8

Reconstructed D mesons in PYTHIA8 simulation
 Reconstructed cluster mass with K± Reconstructed cluster mass with K±
 5000 5000
 w/ LANL FST w/ LANL FST
 D± SIG/BKG: 1.205 MC
 4000 0 0 MC 4000
 D (D ) SIG/BKG: 1.108 ±
 Rec. Ds

• Pixel size:
 Total fit
 D±s (D±s) SIG/BKG: 0.405
 0 0
 3000 Rec. D (D ) 3000

 (100 )3 Rec. D
 ±

 2000 2000

 1000 1000

arXiv:2002.05880 1.83 1.84 1.85 1.86 1.87 1.88 1.89 1.9 1.91 1.94 1.95 1.96 1.97 1.98 1.99 2
 GeV/c2 GeV/c2
 Reconstructed cluster mass with K± Reconstructed cluster mass with K±
 5000 5000
 w/ LANL FST w/ LANL FST
 D± SIG/BKG: 1.969 MC
 MC
• Material
 4000 0 0 4000
 D (D ) SIG/BKG: 1.774 ±
 Rec. Ds
 Total fit
 D±s (D±s) SIG/BKG: 0.882
 budget per 3000 Rec. D (D )
 0 0
 3000

 layer:
 ±
 Rec. D
 2000 2000

 0.3%X0
 1000 1000

 1.83 1.84 1.85 1.86 1.87 1.88 1.89 1.9 1.91 1.94 1.95 1.96 1.97 1.98 1.99 2
 GeV/c2 GeV/c2
 Xuan Li (LANL) 9

New EIC physics observables are under study
• Competing models of nuclear modification in DIS reactions
 with nuclei (e.g HERMES data). Differentiation not possible
 with light hadrons.
– Hadronization inside nuclear
 matter (dashed lines).
– Energy loss of partons,
 hadronization outside (solid lines).
• Heavy mesons have very
 different fragmentation
 functions and formation times
 – Easy to discriminate between
 larger suppression for D/B
 mesons (in-medium
 hadronization) and
 strong/intermediate z
 enhancement (E-loss).
 – Enhanced sensitivity to the
 transport properties of nuclei. arXiv:2002.05880
 Xuan Li (LANL) 10

Haitao Li

Director’s Postdoctoral Fellow at LANL

The Transverse Energy-Energy Correlations in DIS
 Define scattering plane: x-z

 XZ ET,a
 TEEC = d lp!l+a+X (cos la cos )
 AAACa3icbVFdixMxFM2MX2vVtboPivoQLAsrq2VGBX0RFpeCjxXa3UJThjuZTBs2mQzJHaWEefEn+uY/8MX/YKYtorteCJx7zr03uSd5raTDJPkRxdeu37h5a+92787de/v3+w8enjnTWC6m3ChjZzk4oWQlpihRiVltBehcifP84rTTz78I66SpJriuxULDspKl5ICByvrfmAZcWe0no9Fp+4G5RmceWspkhbSgzMmlhswrWlNm5XKFYK35StUxHM9CVWmBezrK/OQlhbb14w6pNiiFUAhMiRKPGDeOsnolwxxoX/1JtwNfZP1BMkw2Qa+CdAcGZBfjrP+dFYY3WlTIFTg3T5MaFx4sSq5E22ONEzXwC1iKeYAVaOEWfuNVSw8DU9DS2HDCihv27w4P2rm1zkNl54y7rHXk/7R5g+X7hZdV3aCo+PaislEUDe2Mp4W0gqNaBwDcyvBWylcQ/MPwPb1gQnp55avg7PUwfTNMPr8dnHzc2bFHnpLn5Iik5B05IZ/ImEwJJz+j/ehR9Dj6FR/ET+Jn29I42vUckH8iPvwNXXG6tA==
 a
 PT,l

 1 + cos
⌧ =
AAACBXicbVDLSsNAFJ3UV62vqEtdDBZBEEpSBd0IRTcuK9gHNKFMppN26GQmzEyEErJx46+4caGIW//BnX/jpM1CWw9cOJxzL/feE8SMKu0431ZpaXllda28XtnY3NresXf32kokEpMWFkzIboAUYZSTlqaakW4sCYoCRjrB+Cb3Ow9EKir4vZ7ExI/QkNOQYqSN1LcPPY0SeAW9UCKcuqceFsqLRzRL61mlb1edmjMFXCRuQaqgQLNvf3kDgZOIcI0ZUqrnOrH2UyQ1xYxkFS9RJEZ4jIakZyhHEVF+Ov0ig8dGGcBQSFNcw6n6eyJFkVKTKDCdEdIjNe/l4n9eL9HhpZ9SHieacDxbFCYMagHzSOCASoI1mxiCsKTmVohHyOShTXB5CO78y4ukXa+5ZzXn7rzauC7iKIMDcAROgAsuQAPcgiZoAQwewTN4BW/Wk/VivVsfs9aSVczsgz+wPn8ADw6Xpg==
 2 Beam Functions: initial state collinear radiation

 Jet Functions: final state collinear radiation
 Soft Functions: recoil of the collinear partilces

 It is similar to the 1-dimensional
 LO Collinear&Soft radiation
 TMD factorization
 Non-zero momentum along y-direction

The Transverse Energy-Energy Correlations in DIS
 200 80
 NLL TEEC s =318 GeV
 100
 NNLL
 60 NNNLL
ds/dln t [pb]

 d /dln (fb)
 0
 40
 -100 e+p s=141 GeV p >20 GeV
 l

 LO sing. LO Full
 dNLO Full
 20
 -200 dNLO sing.
 LO non-sing.
 dNNLO sing. dNLO non-sing.
 -300 0
 -14 -12 -10 -8 -6 -4 -2 0
 lnt
 -14 -12 -10 -8 -6 -4 -2
 ln
 ü We reproduce the full QCD result in the back-to-back region
 ü After resummation the cross section is converged

 We used anomalous dimensions up to 3 loops
 Full control in the infrared region

Name: Zelong Liu (刘泽茏)
E-mail: zelongliu@lanl.gov
Phone: (505) 665-5694

Postdoc in the group of Ivan Vitev

Research Interest:
• Collider Phenomenology (LHC and EIC)
• Jet Physics
• Precision QCD calculation

Fragmentation Functions in Medium

 He
 Kr
 Xe
 Pb

Single Inclusive Production of Hadron at EIC
Ratio of cross section of 4 production in Medium to that in vacuum
 He
 Xe
 Pb

In future, we are also interested in jet production (heavy
flavor jet, jet substructure, et al.) at EIC.

Questions
• What’s the available simulation framework for the
 EIC jets and heavy flavor studies?
• Any background evaluations to be implemented in
 the fast/full simulation?
• How to integrate or combine different approaches
 for the jets and heavy flavor studies?
• How to integrate the detector simulation such as
 the tracking detector work with the EIC jet and
 heavy flavor physics working subgroup?

 Xuan Li (LANL) 17

LANL EIC Jet and Heavy Flavor contribution plan
• Develop new heavy flavor and jet observables and
 provide theoretical predictions for the EIC Yellow
 Report preparation.
• Work on the relevant heavy flavor and jet
 reconstructions in e+p and possible e+A collisions
 in fast simulation and full simulation depending on
 the software developments.
• Work on the physics projections of the proposed
 heavy flavor and jet projections.
• Provide the detector guidance especially for the
 tracking detector based on physics requirements.

 Xuan Li (LANL) 18

Backup

Xuan Li (LANL) 19

LANL EIC program progress (I)
• Initial detector design in fast simulation:
 • Mid-rapidity silicon vertex detector: 3 barrel layers of
 Monolithic Active Pixel Sensor (MAPS) type detector.
 • Forward-rapidity silicon tracking detector (FST): 2 barrel
 layers of MAPS + other silicon detector and 5 forward planes
 of MAPS + other silicon detector.

 B = 3T
 1.0 < < 4.5

X axis (mm)
 Z axis (mm)
 Xuan Li (LANL) 20

LANL EIC program progress (II)
 • Initial track performance from the FST:
 Projected p dependent momentum resolution dp /p Projected p dependent impact parameter b resolution
 T T T 3
 T T
 10

 Impact parameter b resolution (µ m)
Momentum resolution p T

 ± in FST
 T
 dp

 = 1.5
 1
 =2
 = 2.5
 =3
 = 3.5
 =4
 1
 10 2
 10

 T
 2
 10

 ± in FST
 = 1.5
 10
 3 =2 10
 = 2.5
 =3
 = 3.5
 =4
 4
 10
 0 2 4 6 8 10 12 14 16 18 20 22 0 2 4 6 8 10 12 14 16 18 20 22
 ±
 p (GeV/c) ±
 p (GeV/c)
 T T

 • Better than 70 m resolution can be achieved by the initial
 FST design for the transverse decay length bT measurements
 for tracks with pT > 1 GeV/c over the 1.5