Michel Electron Reconstruc0on in DUNE Aleena Rafique , Zelimir - - PowerPoint PPT Presentation

β–Ά
michel electron reconstruc0on in dune
SMART_READER_LITE
LIVE PREVIEW

Michel Electron Reconstruc0on in DUNE Aleena Rafique , Zelimir - - PowerPoint PPT Presentation

May 17, 2023 44 likes β€’305 views

Michel Electron Reconstruc0on in DUNE Aleena Rafique , Zelimir Djurcic on behalf of the DUNE collabora7on CPAD 2019 12/10/2019 Outline Why study Michel electrons TPC Michel Reconstruc7on Energy Reconstruc7on and Michel Energy

slide-1
SLIDE 1

Michel Electron Reconstruc0on in DUNE

Aleena Rafique, Zelimir Djurcic

  • n behalf of the DUNE collabora7on

CPAD 2019 12/10/2019

slide-2
SLIDE 2

Outline

  • Why study Michel electrons
  • TPC Michel Reconstruc7on
  • Energy Reconstruc7on and Michel Energy Spectrum
  • Conclusions

12/10/19

  • A. Rafique, ANL

2

slide-3
SLIDE 3

Michel electrons

  • A. Rafique, ANL

3

  • Michel electrons are electrons from the

decay of muons (0-50 MeV)

  • Common channels (in ProtoDUNE):
  • ΞΌ

+β†’ 𝑓 +πœ‰Μ…πœˆπœ‰π‘“ (80%)

  • ΞΌ

βˆ’β†’ 𝑓 βˆ’πœ‰πœˆπœ‰Μ…π‘“ π‘œΞ³ (20%)

  • Analysis makes use of low energy shower

reconstruc7on– useful for many DUNE analyses

  • Analysis goals:
  • Obtain michel electron energy spectrum
  • Correlate these events with the photon

detector data

12/10/19

slide-4
SLIDE 4

Mo0va0ons

  • Show that ProtoDUNE can use the topological / calorimetric

informa7on provided by the TPC to iden7fy a specific topology [Michel electrons].

  • Reconstruct the energy of Michel electrons using a simple,

preliminary, energy calibra7on and produce a Michel electron energy spectrum.

  • Ideal to study detector's response to electrons in the tens of MeV

energy range

  • Useful for the search of supernova events

12/10/19

  • A. Rafique, ANL

4

slide-5
SLIDE 5

ProtoDUNE-SP LArTPC

Y= X= Z=

12/10/19

  • A. Rafique, ANL

5

beam

slide-6
SLIDE 6

Candidate michel event displays in ProtoDUNE

Other cosmics Candidate muon

12/10/19

  • A. Rafique, ANL

6

Candidate michel Candidate michel Candidate muon

slide-7
SLIDE 7

Michel event ini0al selec0on

  • 1. T0 tagged

(cathode crossing) tracks beam x

.

y z

12/10/19

  • A. Rafique, ANL

7

(0,0,0)

slide-8
SLIDE 8

Michel event ini0al selec0on

  • 1. T0 tagged

(cathode crossing) tracks

12/10/19

  • A. Rafique, ANL

8

  • 2. Tracks

star7ng in detector edges

30 cm 30 cm 30 cm

slide-9
SLIDE 9

Michel event ini0al selec0on

  • 1. T0 tagged

(cathode crossing) tracks

12/10/19

  • A. Rafique, ANL

9

  • 2. Tracks

star7ng in detector edges

  • 3. Tracks

ending in FV

slide-10
SLIDE 10

Michel event ini0al selec0on

  • 1. T0 tagged

(cathode crossing) tracks

12/10/19

  • A. Rafique, ANL

10

  • 2. Tracks

star7ng in detector edges

  • 3. Tracks

ending in FV

  • 4. Tracks

ending not around APA bounds Drig direc7on

e e e e e e e e e e

slide-11
SLIDE 11

Candidate muon selec0on

  • Remove broken tracks
  • A. Rafique, ANL

11 12/10/19

d Candidate muon

MC Sample purity = 27%

ΞΈ

slide-12
SLIDE 12

Candidate muon selec0on

  • Remove broken tracks
  • Select candidate muon tracks

with > 75 cm length

  • A. Rafique, ANL

12 12/10/19

MC Sample purity = 27%

ProtoDUNE simula7on

True MC normalized to reco MC

length > 75 cm

100 200 300 400 500 600 700 800 900 1000 candidate muon track length (cm)

500 1000 1500 2000 2500 3000 3500 4000 4500

# Events ProtoDUNE-SP Simulation

reco events true michel events

100 200 300 400 500 600 700 800 900 1000 candidate muon track length (cm)

500 1000 1500 2000 2500 3000 3500 4000 4500

# Events ProtoDUNE-SP Simulation, Preliminary

reco events true michel events

slide-13
SLIDE 13

Candidate muon selec0on

  • Remove broken tracks
  • Select candidate muon tracks

with > 75 cm length

  • Minimum hit 7me > 200 7cks
  • A. Rafique, ANL

13 12/10/19

MC Sample purity = 63%

ProtoDUNE simula7on ProtoDUNE simula7on

e e e e e e

ProtoDUNE simula7on

100 200 300 400 500 600 700 800 900 1000 candidate muon track length (cm)

500 1000 1500 2000 2500 3000 3500 4000 4500

# Events ProtoDUNE-SP Simulation

reco events true michel events

1000 2000 3000 4000 5000 6000 Minimum hit peak time

1 10

2

10

3

10

4

10

5

10

# Events ProtoDUNE-SP Simulation

reco events true michel events

1000 2000 3000 4000 5000 6000 Minimum hit peak time

1 10

2

10

3

10

4

10

5

10

# Events ProtoDUNE-SP Simulation, Preliminary

reco events true michel events

slide-14
SLIDE 14

Candidate michel selec0on

  • Remove broken tracks
  • Select candidate muon tracks

with > 75 cm length

  • Minimum hit 7me > 200 7cks
  • Nearby hit count > 5 & < 40
  • A. Rafique, ANL

14 12/10/19

MC Sample purity = 86%

ProtoDUNE simula7on ProtoDUNE simula7on ProtoDUNE simula7on

R = 1 c m

Drig dir (7me)

20 40 60 80 100 120 140 160 180 200 number of hits close to the track end point

1 10

2

10

3

10

4

10

# Events ProtoDUNE-SP Simulation

reco events true michel events

20 40 60 80 100 120 140 160 180 200 number of hits close to the track end point

1 10

2

10

3

10

4

10

# Events ProtoDUNE-SP Simulation, Preliminary

reco events true michel events

Wire (beam direc7on)

slide-15
SLIDE 15

Candidate michel selec0on

  • Remove broken tracks
  • Select candidate muon tracks

with > 75 cm length

  • Minimum hit 7me > 200 7cks
  • Nearby hit count > 5 & < 40
  • Closest reco shower distance

< 10 cm

  • A. Rafique, ANL

15 12/10/19

Final MC Sample purity = 89%

ProtoDUNE simula7on ProtoDUNE simula7on ProtoDUNE simula7on ProtoDUNE simula7on

R = 1 c m

5 10 15 20 25 30 shower start distance from muon end point (cm)

1000 2000 3000 4000 5000 6000 7000

# Events ProtoDUNE-SP Simulation

reco events true michel events

5 10 15 20 25 30 shower start distance from muon end point (cm)

1000 2000 3000 4000 5000 6000 7000

# Events ProtoDUNE-SP Simulation, Preliminary

reco events true michel events

slide-16
SLIDE 16

Michel energy reconstruc0on

Where Qi = charge deposited on a hit Cx,Cyz = correc7on factors, remove non-uniformity in dQ/dx values Wion = 23.6e-6 MeV Norm_factor = normalizes the dQ/dx values to the dQ/dx at anode Calib_const = converts dQ/dx in ADC/cm into dE/dx in MeV/cm Recomb_factor = 0.7; to incorporate the recombina7on effects

12/10/19

  • A. Rafique, ANL

16

𝐹 = βˆ‘,-./00. 20345 6,78

9, βˆ—;< βˆ—;=> βˆ—?,/4βˆ—@/AB_D3.7/A ;30,E_./487βˆ—F5./BE_D3.7/A

20 40 60 80 100 120 140 160 180 200 True michel reconstructed energy (MeV)

200 400 600 800 1000 1200 1400 1600 1800 2000

# Events

true_michel_mc Entries 12028 Mean 33.58 Std Dev 14.68

ProtoDUNE-SP Simulation, Preliminary

slide-17
SLIDE 17

True michel reco hits Other hits of the event

Michel energy reconstruc0on

cone

12/10/19

  • A. Rafique, ANL

17

  • Defined a cone at the end point
  • f the parent muon

ΞΈ L

Candidate muon

20 40 60 80 100 120 140 160 180 cone angle 100

  • 80
  • 60
  • 40
  • 20
  • 20

40 60 80 100 cone length (cm)

20 40 60 80 100 120

3

10 Β΄

ProtoDUNE-SP Simulation, Preliminary 20 40 60 80 100 120 140 160 180 cone angle 50

  • 40
  • 30
  • 20
  • 10
  • 10

20 30 40 50 cone length (cm)

2000 4000 6000 8000 10000

ProtoDUNE-SP Simulation, Preliminary 20 40 60 80 100 120 140 160 180 200 Reconstructed energy (MeV)

200 400 600 800 1000 1200 1400 1600 1800 2000

# Events

true_michel_mc Entries 12028 Mean 33.58 Std Dev 14.68 reco_mc Entries 12028 Mean 31.09 Std Dev 21.43

ProtoDUNE-SP Simulation, Preliminary

reco MC true michel MC

Missing michel hits

Where ΞΈ=0 corresponds to the hits along the nearby wire hits

slide-18
SLIDE 18

Recovering michel hits from parent muon

  • Look at the charge deposi7on at the

last 10 hits of reco muons

  • Charge dep = (Qi – Qi-1)/Qi-1
  • Ager removing hits beyond the

maximum truncated charge value

  • Recovered 6% of the total missing

michel hits

12/10/19

  • A. Rafique, ANL

18 2 4 6 8 10 Hit number 1

  • 0.8
  • 0.6
  • 0.4
  • 0.2
  • 0.2

0.4 0.6 0.8 1 Charge deposition

100 200 300 400 500 600 700 800 900

ProtoDUNE-SP Simulation, Preliminary

2 4 6 8 10 Hit number 1

  • 0.8
  • 0.6
  • 0.4
  • 0.2
  • 0.2

0.4 0.6 0.8 1 Charge deposition

100 200 300 400 500 600 700 800 900

ProtoDUNE-SP Simulation, Preliminary

2 4 6 8 10 Hit number 1

  • 0.8
  • 0.6
  • 0.4
  • 0.2
  • 0.2

0.4 0.6 0.8 1 Charge deposition

100 200 300 400 500 600 700 800 900

ProtoDUNE-SP Simulation, Preliminary

True muon Before recovery Ager recovery

Last hit

slide-19
SLIDE 19

Michel energy spectrum MC

12/10/19

  • A. Rafique, ANL

19

20 40 60 80 100 120 140 160 180 200 Reconstructed energy (MeV)

200 400 600 800 1000 1200 1400 1600 1800 2000

# Events

true_michel_mc

Entries 12028 Mean 33.58 Std Dev 14.68

reco_mc

Entries 12028 Mean 36.8 Std Dev 22.33

ProtoDUNE-SP Simulation, Preliminary

reco MC true michel MC 10 20 30 40 50 60 70 80 90 100 True michel reco energy (MeV) 10 20 30 40 50 60 70 80 90 100 Michel reco energy (MeV)

20 40 60 80 100 120

ProtoDUNE-SP Simulation, Preliminary

200

  • 150
  • 100
  • 50
  • 50

100 150 200 True michel reco - reco energy (MeV) 500 1000 1500 2000 2500 3000 # Events

micheltrueparEreso Entries 12028 Mean 3.431

  • Std Dev

23.55

ProtoDUNE-SP Simulation, Preliminary

slide-20
SLIDE 20

Michel energy data/MC comparison

12/10/19

  • A. Rafique, ANL

20

20 40 60 80 100 120 140 160 180 200

Reconstructed energy (MeV)

200 400 600 800 1000 1200

Entries

MC_copy

Entries 13509 Mean 38.19 Std Dev 24.3

Data Entries 7760 Mean 38.89 Std Dev 24.52 Data Entries 7760 Mean 38.89 Std Dev 24.52

ProtoDUNE-SP, Preliminary

MC Data

20 40 60 80 100 120 140 160 180 200

Reconstructed energy (MeV)

200 400 600 800 1000 1200

Entries

MC_copy

Entries 13509 Mean 38.19 Std Dev 24.3

Data Entries 7760 Mean 38.89 Std Dev 24.52 Data Entries 7760 Mean 38.89 Std Dev 24.52

ProtoDUNE-SP, Preliminary

MC Bkg MC Data

Very good agreement in ProtoDUNE data and simula7on

slide-21
SLIDE 21

Summary

  • Michel energy spectrum is obtained with a high purity

michel sample

  • Good agreement in simula7on and ProtoDUNE data is
  • bserved for michel energy spectrum
  • In the near future, will study how photon detector sees

these michel events

  • Note on the analysis is coming soon. Stay tuned!

12/10/19

  • A. Rafique, ANL

21

slide-22
SLIDE 22

Thank you for your aFen0on

12/10/19

  • A. Rafique, ANL

22

slide-23
SLIDE 23

Backup slides

12/10/19

  • A. Rafique, ANL

23

slide-24
SLIDE 24

Deep Underground Neutrino Experiment (DUNE)

DUNE will be a world-class neutrino observatory

  • 1300 km baseline
  • Consists of a large LArTPC far (40 kTon) and near detectors
  • A broad and rich physics program: Neutrino oscilla7ons, CP viola7on

searches in the neutrino sector, neutrino mass hierarchy, supernova neutrinos, baryon number viola7on searches

  • The world’s most intense neutrino beam from Fermilab
  • A deep underground site, massive liquid argon detectors and a precision

near detector

12/03/19

  • A. Rafique, ANL

24

slide-25
SLIDE 25

ProtoDUNE Single Phase

  • ~7Γ—6Γ—7 m3 (770 tons of LAr) in charged test

beam at CERN

  • ProtoDUNE-SP opera7ng since September

2018

  • Accumula7ng test-beam data to

understand/calibrate response of detector to different par7cle species

  • A crucial part of the DUNE effort towards

the construc7on of the first DUNE

  • Prototyping produc7on and installa7on

procedures for DUNE Far Detector Design

  • Valida7ng design from perspec7ve of basic

detector performance

  • Demonstra7ng long term opera7onal

stability of the detector ProtoDUNE-SP at CERN neutrino plaworm

12/03/19

  • A. Rafique, ANL

25

slide-26
SLIDE 26

Principle of LArTPC

LArTPCs make 3D reconstrucOon possible.

  • Wire planes give 2D posi7on informa7on
  • The third dimension is obtained by combining 7ming

informa7on (t0) with drig velocity (vd) β†’ hence, a β€œTime projecOon chamber”

26 12/03/19

  • A. Rafique, ANL

ProtoDUNE TPC Y= X= Z=