Pion scattering with the Pion scattering with the LArIAT experiment - - PowerPoint PPT Presentation

pion scattering with the pion scattering with the lariat
SMART_READER_LITE
LIVE PREVIEW

Pion scattering with the Pion scattering with the LArIAT experiment - - PowerPoint PPT Presentation

Pion scattering with the Pion scattering with the LArIAT experiment LArIAT experiment Justin Hugon (On behalf of the LArIAT experiment) Louisiana State University 1 Liquid Argon in a Test Beam (LArIAT) Liquid Argon in a Test Beam (LArIAT)


slide-1
SLIDE 1

1

Pion scattering with the Pion scattering with the LArIAT experiment LArIAT experiment

Justin Hugon (On behalf of the LArIAT experiment) Louisiana State University

slide-2
SLIDE 2

2

Liquid Argon in a Test Beam (LArIAT) Liquid Argon in a Test Beam (LArIAT)

47 cm 47 cm Drift Direction Drift Direction Beam Direction Beam Direction 90 cm 90 cm 40 cm 40 cm

LArIAT TPC LArIAT TPC

170 L 0.25 tons

  • f LAr

Reuse the ArgoNEUT TPC in a charged particle beam at Fermilab

slide-3
SLIDE 3

3

Liquid Argon in a Test Beam (LArIAT) Liquid Argon in a Test Beam (LArIAT)

Changes from ArgoNEUT:

  • New wireplanes
  • Cold front-end electronics fom MicroBooNE

C a t h

  • d

e p l a n e Wi r e / a n

  • d

e p l a n e s R e a d

  • u

t A S I C s

slide-4
SLIDE 4

4

LArIAT Goals LArIAT Goals

  • Physics Goals

– Hadron-Ar interaction cross sections

  • p+/--Ar to support ν cross-sections
  • K+/- - Ar, supporting nucleon decay
  • Geant4 validation

– e/g shower identification capabilities – Anti-proton annihilation at rest

  • Similar to BSM n-n oscillation signature

– Particle sign determination in the absence of a magnetic field, utilizing topology

  • e.g. decay vs capture
  • R&D Goals

– Ionization and scintillation light studies

  • Charge deposited vs. light collected for stopping particles of known energy

– Optimization of particle ID techniques – LArTPC event reconstruction

  • Compare 3mm, 4mm, 5mm wire pitch
slide-5
SLIDE 5

5

  • Fermilab’s main injector sends 120 GeV protons to

the Fermilab test-beam facility (FTBF)

  • FTBF creates a tunable secondary beam, 8 GeV to

80 GeV, directed toward the LArIAT experimental hall

LArIAT Testbeam LArIAT Testbeam

slide-6
SLIDE 6

6

  • Secondary target further reduces beam

energy

  • Instrumented beamline identifies and

characterizes particles

LArIAT Tertiary Beamline LArIAT Tertiary Beamline

D i p

  • l

e M a g n e t s H a l

  • v

e t

  • M

u

  • n

P u n c h t h r

  • u

g h V e t

  • M

u

  • n

R a n g e S t a c k M u l t i

  • w

i r e p r

  • p
  • r

t i

  • n

a l c h a m b e r s ( M WP C s ) C

  • l

l i m a t

  • r

s A e r

  • g

e l C e r e n k

  • v

C

  • u

n t e r s

T P C

S e c

  • n

d a r y b e a m

p’

s ( 8

  • 8

G e V )

C u t a r g e t

T i m e

  • f

F l i g h t ( T O F )

slide-7
SLIDE 7

7

Wire chambers reconstruct the position and momentum of the particles in the beamline

LArIAT Beamline: Wire Chambers LArIAT Beamline: Wire Chambers

Wire chamber reconstructed momentum compared to simulation

LArIAT Preliminary

slide-8
SLIDE 8

8

LArIAT Beamline: Time of Flight LArIAT Beamline: Time of Flight

2 scintillator counters w/ ~1ns sampling, provide the time of flight (TOF)

slide-9
SLIDE 9

9

LArIAT Beamline Detectors LArIAT Beamline Detectors

Combining the momentum and TOF allows for p/m/e, K, proton separation Additionally, using the known masses of the K and proton we can constrain the momentum scale to 1.5%

TOF vs reconstructed momentum

LArIAT Preliminary

slide-10
SLIDE 10

10

Matching Beamline to the TPC Matching Beamline to the TPC

  • We can take this track

reconstructed in the beamline and extrapolate it to the LArTPC and look for a match

– We match in both position

(+/- 5cm about the mean) and angle (< 10o)

10

slide-11
SLIDE 11

11

Calibrating our sample Calibrating our sample

  • Using the first few centimeters of a matched track we

can characterize the dE/dX response as a function of the track’s initial momentum in both data and simulation

  • Calibrate detector response to follow Bethe-Block

formula by selecting events with different particle type and momentum

  • Calibrate using pions; check on kaons/protons

LArIAT Preliminary

slide-12
SLIDE 12

12

  • Cosmic-ray muons provide

another calibration sample

– Width of dE/dx distribution

can be compared between data and simulation

Smearing Simulated Data Smearing Simulated Data

  • Additional 20% smearing

makes the simulation match the data

LArIAT Preliminary Events / (0.1 MeV/cm) LArIAT Preliminary

slide-13
SLIDE 13

13

  • Our MC allows us to

estimate what our fractional beam composition and our selection efficiencies are for the various particle species

Pion Event Selection Pion Event Selection

LArIAT Preliminary

slide-14
SLIDE 14

14

A b s

  • r

p t i

  • n

C a n d i d a t e ( π

  • >

3 p )

LArIAT Data

C h a r g e E x c h a n g e C a n d i d a t e

LArIAT Data LArIAT Data

  • The total p––Argon Cross-Section includes

σTotal=σelastic+σinelastic+σch-exch+σabsorp.+σp -production

π P r

  • d

u c t i

  • n

C a n d i d a t e

LArIAT Data

I n e l a s t i c S c a t t e r i n g C a n d i d a t e

LArIAT Data

+ + +

Pion Cross-Section Pion Cross-Section

+

LArIAT Data

E l a s t i c S c a t t e r i n g C a n d i d a t e

LArIAT Data

slide-15
SLIDE 15

15

LArIAT Data

  • Backgrounds are:

π D e c a y C a n d i d a t e

LArIAT Data

π C a p t u r e C a n d i d a t e

LArIAT Data

Mu

  • n

B a c k g r

  • u

n d

LArIAT Data

Note: Pion decay backgrounds are small component which remain in our result. Capture dominates the lowest energy bin and is thus excluded

Pion Cross-Section Pion Cross-Section

LArIAT Simulation LArIAT Simulation Pion Interaction Type per Kinetic Energy Pion Interaction Fraction per Kinetic Energy

slide-16
SLIDE 16

16

  • The probability of the pion interacting is thus:

where we measure the probability of interacting for that thin slab as the ratio of the number of interacting pions to the number of incident pions:

PSurvival=e

−σn z

PInteracting=1−PSurvival

Ninteracting NIncident =PInteracting=1−e

−σ nz

Thin Slice TPC Method Thin Slice TPC Method

  • Generally, the survival

probability of a pion traveling through a thin slab of argon is given by: Where σTOT is the cross-section per nucleon, z is the depth of the slab, and n is the density

slide-17
SLIDE 17

17

  • Thus you can extract the pion cross-section as a

function of energy as

Where n = rNA / A

  • Using the granularity of the LArTPC, we can treat the

wire-to-wire spacing as a series of “thin-slab” targets if we know the energy of the pion incident to that target PInteracting=1−(1−σ nδ z+...) σ(E)≈ 1 nz PInteracting=1 nz N interacting NIncident

Thin Slice TPC Method Thin Slice TPC Method

slide-18
SLIDE 18

18 Analyze the reconstructed tracks

  • Now we have a matched WC track and TPC track
  • We calculate the

p-candidate's initial kinetic energy as we take into account energy loss due to material upstream of the TPC (argon, steel, beamline detectors, etc)

  • We then follow p-candidate track treating each

point as a “thin slice” of argon which the pion is incident to at a known energy

KEi=√ p

2+mp 2−mp−EFlat

Interacting Incident

Kinetic Energy (MeV) Kinetic Energy (MeV)

KE Interaction=KE i−∑

i=0 nSpts

dE/dX i×Pitchi

Pion Cross-Section Pion Cross-Section

slide-19
SLIDE 19

19

Interacting Incident

Kinetic Energy (MeV) Kinetic Energy (MeV)

KE Interaction=KE i−∑

i=0 nSpts

dE/dX i×Pitchi

Pion Cross-Section Pion Cross-Section

Interacting Incident

Kinetic Energy (MeV) Kinetic Energy (MeV)

We ignore other tracks in the event not matched to the Wire Chamber Track

When you encounter the interaction point you now fill the interacting and incident histogram for the energy the pion has at that point You now repeat this process for your entire sample

slide-20
SLIDE 20

20

  • You now take the ratio of

these two histograms to extract the cross-section

Pion Cross-Section Pion Cross-Section

Interacting Incident

Where n = rNA / A

σ(E)≈ 1 nz PInteracting=1 nz N interacting NIncident

slide-21
SLIDE 21

21

Pion Cross-Section Pion Cross-Section

Systematics Considered Here dE/dX Calibration: 3% (previously was 5%) (previously was 5%) Energy Loss Prior to entering the TPC: 3.5% Through Going Muon Contamination: 3% Wire Chamber Momentum Uncertainty: 3%

slide-22
SLIDE 22

22

Toward Exclusive Pion Channels Toward Exclusive Pion Channels

  • Working on absorption + charge

exchange: p + Ar → 0π± + X

– Useful for modeling contamination of ν CC

QE from CC RES where a π is absorbed in the interaction nucleus

– Need to identify outgoing pions v. protons

A b s

  • r

p t i

  • n

C a n d i d a t e ( π

  • >

3 p )

LArIAT Data

C h a r g e E x c h a n g e C a n d i d a t e

LArIAT Data LArIAT Data

Signal Events: 0 Secondary π± Background Events: Contain Secondary π±

slide-23
SLIDE 23

23

Likelihood-Based Particle ID Likelihood-Based Particle ID

  • Likelihood of dE/dx versus

residual range of each track hit

– Constructed from simulated tracks – Evaluate using likelihood-ratio of

all hits on a track

Proton Likelihood Pion Likelihood

slide-24
SLIDE 24

24

Anti-proton Annihilation Anti-proton Annihilation

p-candidate

g g

g g g g g g

halo muon

multi-nucleon

  • LArIAT has identified

O(20) anti-proton annihilation at rest candidates

– O(70) annihilation in flight

  • Similar to BSM n-n
  • scillation signature

– DUNE planning search

LArIAT Data LArIAT Data p-candidate Work ongoing to reconstruct these final state topologies

slide-25
SLIDE 25

25

Kaon Cross-sections Kaon Cross-sections

  • Inclusive K+ cross-

section has O(2000) Elastic/Inelastic interactions identified

– Inclusive cross-section

coming soon

  • First time measured on

argon

– Work ongoing to

reconstruct these final state topologies

  • DUNE plans search for proton decay:

p → K+ ν

  • Cross-section information will help ensure

signal efficiency is modeled properly

slide-26
SLIDE 26

26

  • Inclusive p--Argon Cross-section

– New result has x100 the initial statistics

  • Inclusion of Run-II data
  • Tuning of reconstruction cuts and improvement in dE/dX

calibration

– Paper in preparation

Conclusions Conclusions

  • LArIAT recently completed its 3rd

physics run

– Run-I / Run-II: 4mm wire pitch

  • Hadronic cross-sections
  • Scintillation Light R&D

– Run-III: 3mm / 5mm wire pitch

comparison

  • LArTPC particle ID R&D
  • New mesh cathode
  • New ARAPUCA Light Detection System
slide-27
SLIDE 27

27

  • Many other physics results following close

behind this result

– Inclusive K+ Cross-section – Inclusive p+ Cross-section

  • Absorption and charge exchange exclusive channels coming

along too

– Anti-proton annihilation at rest – e/g shower characterization

  • Inclusion of 3mm/5mm wire pitch comparisons

Future Plans Future Plans

slide-28
SLIDE 28

28

Thank you!

slide-29
SLIDE 29

29

Backup Slides

slide-30
SLIDE 30

30 Neutrino interaction in LAr produces ionization and scintillation light Drift the ionization charge in a uniform electric field Read out charge and light produced using precision wires and PMT's

Electric Field Electric Field Electric Field

g g g g g g

nm

g

Induction Induction Plane Plane Collection Collection Plane Plane

⊕ ⊕ Drift

Drift Time Time = ✔ 3D imaging with mm 3D imaging with mm space resolution space resolution ✔ Calorimetry information Calorimetry information ✔ PID capabilities PID capabilities

Wire Number Time Tick Wire Number Time Tick

LArIAT Data LArIAT Data K+ → µ+ → e+ Candidate K+ → µ+ → e+ Candidate

Bragg peak

Liquid Argon Time Projection Chamber Liquid Argon Time Projection Chamber

slide-31
SLIDE 31

31

T i me

  • f

F l i g h t ( T O F )

H a l

  • v

e t

  • M

u

  • n

R a n g e S t a c k

g e l k

  • v

t e r s

T P C

Energy Corrections Energy Corrections

  • Adding up all the energy which a pion loses in

the region before it enters the TPC (TOF, Halo, Cryostat Cryostat, Argon Argon) gives us the “energy loss” by the pion in the upstream region

Gaussian Fit from 35 – 55 MeV

KEi=√ p

2+mp 2−mp−EFlat

slide-32
SLIDE 32

32

n=1.11 Aerogel n=1.057 Aerogel

200-300 MeV/c µ π µ π 300-400 MeV/c µ π µ π

A G

T O F

W C 1 W C 2 W C 3 W C 4

T O F

LArTPC

Muon Range Stack

Magnet Magnet

α α

LArIAT Beamline Detectors LArIAT Beamline Detectors

✔ Allows to perform p/m separation

  • ver a range of momentum

✔ Currently under investigation

π π μ μ

✔ Four layers of XY planes sandwiched between (pink) steel slabs ✔ Each plane is composed by 4 scintillating bars connected to a PMT ✔ Allows to discriminate p/m exiting the cryostat ✔ Currently under investigation

slide-33
SLIDE 33

33

Inside the cryostat: TPC and light collection system Inside the cryostat: TPC and light collection system

1 2 3 4

Cathode Cathode WirePlanes WirePlanes Pulse shaping Pulse shaping and amplifying and amplifying cold ASICs cold ASICs

Light Light Collection Collection System port System port

  • 1. PMT: Hamamatsu R-11065 (3” diameter)
  • 2. PMT: ETL D757KFL (2” diameter)
  • 3. SiPM: SensL MicroFB-60035 w/preamp
  • 4. SiPM: Hmm. S11828-3344M 4x4 array (Run I)

SiPM: Hmm. VUV-sensitive (Run II)

slide-34
SLIDE 34

34

Conversion-on-PMTs only LArIAT solution

Light Collection System Light Collection System

TPB Field Cage Wall Reflector

Reflector Foils Reflector Foils

Credit: W. Foreman

✔ Wavelength shifting (evaporated) reflected foils on the four field cage walls

✔ Technique borrowed from dark matter experiments

✔ Provides greater (~ 40 pe/MeV at zero field) and more uniform light yield respect to “conversion-on-PMTs-only” light systems ✔ R&D for future neutrino experiments as a way to improve calorimetry and triggering

Beam direction

PMTs

x

  • y

z

Beam direction Beam direction

34

slide-35
SLIDE 35

35

New ARAPUCA Light Collection System New ARAPUCA Light Collection System

  • Dichoric filter +

wavelength shifter

– Trap light inside device

  • Inner walls made of

Teflon

– Trapped light reflected

until detected by SiPM

slide-36
SLIDE 36

36

New ARAPUCA Light Collection System New ARAPUCA Light Collection System

  • ARAPUCA mounted near

existing PMTs

– Compare ARAPUCA performance

to PMTs

2x Ganged SiPM

slide-37
SLIDE 37

37

  • We begin by looking at the bin content of the cross-section from MC

– Here we show events / 50 MeV bin to mimic the binning used in the data – Plot the true kinetic energy

  • Pion captrure-at-rest dominate in the lowest energy bin (0 MeV < KE < 50 MeV)

– Constitutes ~80% of the interactions in that bin – This is not a process we want to include in the cross-section measurement

Cross-Section Cross-Section

slide-38
SLIDE 38

38

  • About 1% of the time the pion actually stops

before reaching the TPC

– The remaining portion there is actually an interaction

What happens in the upstream What happens in the upstream

slide-39
SLIDE 39

39

Validation Plots Validation Plots