1 GeV Beam dE/dx Study ProtoDUNE Analysis Meeting- 02/21/2019 Owen - - PowerPoint PPT Presentation

▶
1 gev beam de dx study
SMART_READER_LITE
LIVE PREVIEW

1 GeV Beam dE/dx Study ProtoDUNE Analysis Meeting- 02/21/2019 Owen - - PowerPoint PPT Presentation

Jan 30, 2023 244 likes •607 views

1 GeV Beam dE/dx Study ProtoDUNE Analysis Meeting- 02/21/2019 Owen Goodwin 1 Introduction Owen Goodwin 02/21/2019 Looking at dE/dx of reconstructed 1 GeV candidate tracks. Using all of Run 5387 (defname

slide-1
SLIDE 1

1 GeV Beam dE/dx Study

ProtoDUNE Analysis Meeting- 02/21/2019 Owen Goodwin

1

slide-2
SLIDE 2

Owen Goodwin 02/21/2019

Introduction

  • Looking at dE/dx of reconstructed 1 GeV candidate tracks.
  • Using all of Run 5387 (defname runset-5387-reco-unified-hv-180kV-beam-1GeV-v0 );

○ 1 GeV Beam Momentum ○ 180 kV drift field ○ Hadron Trigger

  • Use Justin Hugon’s beamline filter to select 1 GeV Pion events. Details here

https://indico.fnal.gov/event/19185/contribution/2/material/slides/0.pdf

  • Then use protoana::ProtoDUNEPFParticleUtils to get reconstructed TPC

beam particle

2

slide-3
SLIDE 3

Owen Goodwin 02/21/2019

Selection Cuts

3

slide-4
SLIDE 4

Owen Goodwin 02/21/2019

Calorimetry

  • Use Ajib’s calorimetry correction to

calibrate, described here. https://docs.dunescience.org/cgi-bin/pri vate/RetrieveFile?docid=12997&filenam e=Stopping_muon_calibration.pdf&versi

  • n=1
  • Uses cathode crossing T0 tagged

stopping muons. Calibration is for collection plane.

  • Done for run 5387 (this run)

4

slide-5
SLIDE 5

Owen Goodwin 02/21/2019

Residual Range vs dE/dx

5

  • Clear MIP region.
  • Most Particles don’t

stop, as expected, Pions will interact

slide-6
SLIDE 6

Owen Goodwin 02/21/2019

Selecting Beam Stopping Muons

6

slide-7
SLIDE 7

Owen Goodwin 02/21/2019

Stopping muons

Selecting tracks with track length compatible with Stopping Muon. Clear Muon like bragg peak shown, matches prediction. Beam cross check of cosmic calibration

7

slide-8
SLIDE 8

Owen Goodwin 02/21/2019

Plotting all candidates with a logarithmic Z scale show, small component consistent with stopping proton Bragg peak. Possibly protons produced by beam Pion interactions? Going to look at some event dispays

8

slide-9
SLIDE 9

Owen Goodwin 02/21/2019

Monte Carlo

  • 1 GeV Beam MC.
  • Select events where beam

particle is true pion

  • Use Ajib’s calibration for

each MCC11 sample (No SCE, SCE, SCE+FLF)

9

slide-10
SLIDE 10

Owen Goodwin 02/21/2019

Residual Range vs dE/dx (FLF MC)

10

Very low stats (~300 true pions events in each sample) 4 true pion events in stopping region

slide-11
SLIDE 11

Owen Goodwin 02/21/2019

Plans

  • Would like to run data/MC comparisons with higher stats MC.
  • Longer term goal to reconstruct energy spectrum of Michel electrons from

selected beam muons

11

slide-12
SLIDE 12

Owen Goodwin 02/21/2019

BACKUP

12

slide-13
SLIDE 13

Owen Goodwin 02/21/2019

Rough PID. Interacting pion->proton?

13

slide-14
SLIDE 14

Owen Goodwin 02/21/2019

Beam Start Position

Flipped TPC tracks with Start Z > End Z ~32 cm offset in Z as seen in other studies

14

slide-15
SLIDE 15

Owen Goodwin 02/21/2019

Beam Start Position

~32 cm offset in Z as seen in other studies

15

slide-16
SLIDE 16

Owen Goodwin 02/21/2019

Directional Match

Angle between beamline particle and TPC particle. Closest Match contains only the best matching beamline particle for each Evt

16

slide-17
SLIDE 17

Owen Goodwin 02/21/2019

Track Length

17

Distance Cuts: 6cm< x diff < 14cm and y diff < 10cm and 29cm<z diff <36 cm Angle cut: Cos(theta)>0.9 Both Cuts remove extreme high length tail.

slide-18
SLIDE 18

Owen Goodwin 02/21/2019

Track Length

Spike at 200cm cause by broken tracks at end of first APA (230cm in Z direction)

18

230cm

slide-19
SLIDE 19

Owen Goodwin 02/21/2019

Event Display (Run 5387, Ev 60537)

Track Split into two Reonstructed tracks

19

APA 3 (TPC 1) APA 2 (TPC 5)

slide-20
SLIDE 20

Owen Goodwin 02/21/2019

Event Display (Run 5387, Ev 60537)

Track Split into two Reconstructed tracks

20

slide-21
SLIDE 21

Owen Goodwin 02/21/2019

(Run 5387, Ev 49259)

21

slide-22
SLIDE 22

Owen Goodwin 02/21/2019

Track Length

Peak much beyond Pion-LAr interaction length at 1 GeV (<1m) Muon Contamination 1GeV mom Stopping Muon expected range using CSDA tables (http://pdg.lbl.gov/2012/AtomicNuclearProperties/MUON_ELOSS_TABLES/muonloss_289.pdf) is ~395cm

22

slide-23
SLIDE 23

Owen Goodwin 02/21/2019

Track Length

Calculate CSDA range of particle (muon assumption) using Beamline momentum. Divide Track Length by this. Stopping muon peak just below 1.

23

slide-24
SLIDE 24

Owen Goodwin 02/21/2019

Run 5387, Ev 8588

24

APA 3 (TPC 1) APA 2 (TPC 5)

slide-25
SLIDE 25

Owen Goodwin 02/21/2019

Run 5387, Ev 8588

Stopping Muon with reconstructed Michel electron?

25

slide-26
SLIDE 26

Owen Goodwin 02/21/2019

Summary + Possible Next Steps

Study variables separated by track end processes Matching broken tracks Muons?

26

slide-27
SLIDE 27

Owen Goodwin 02/21/2019

27

slide-28
SLIDE 28

Owen Goodwin 02/21/2019

28

slide-29
SLIDE 29

Owen Goodwin 02/21/2019

Track Reco Wrong way

7594/13122 = 57.8%

  • f Pion Candidate tracks reconstructed “wrong way” round.

(Start Z position > End Z position)

29

I flip the tracks the wrong way round. Plots shows z start position before and after flipping

slide-30
SLIDE 30

Owen Goodwin 02/21/2019

Beamline Particle info

30

Beamline particles have opposite direction to expected? Going in negative Z direction

slide-31
SLIDE 31

Owen Goodwin 02/21/2019

Beamline Particle info

31

slide-32
SLIDE 32

Owen Goodwin 02/21/2019

32

Orig is before I flip the backwards going TPC tracks Corrected after. Corrected and flipped is after I also flip direction of all beamline particles

slide-33
SLIDE 33

Owen Goodwin 02/21/2019

Pion Cand StartPos

33

slide-34
SLIDE 34

Owen Goodwin 02/21/2019

Pion Cand EndPos

34

z=230cm edge of first APA

slide-35
SLIDE 35

Owen Goodwin 02/21/2019

35