Michel analysis updates Aleena Rafique, Zelimir Djurcic ProtoDUNE - - PowerPoint PPT Presentation

▶
michel analysis updates
SMART_READER_LITE
LIVE PREVIEW

Michel analysis updates Aleena Rafique, Zelimir Djurcic ProtoDUNE - - PowerPoint PPT Presentation

Aug 29, 2022 229 likes •594 views

Michel analysis updates Aleena Rafique, Zelimir Djurcic ProtoDUNE sim/reco meeting 03/18/2020 Updates 1. Changing 2D 3D distance for additional muon hits 2. Michel energy spectrum with cone-only hits Data and MC comparison 3. Looking for

slide-1
SLIDE 1

Michel analysis updates

Aleena Rafique, Zelimir Djurcic ProtoDUNE sim/reco meeting 03/18/2020

slide-2
SLIDE 2

Updates

  • 1. Changing 2Dà3D distance for additional muon hits
  • 2. Michel energy spectrum with cone-only hits Data and MC

comparison

  • 3. Looking for kink around the muon end position
  • 4. Recombination factor Investigation
  • 5. Implementing new wire charge integral method
  • 6. Calibrate to true Michel distribution

3/18/20

  • A. Rafique, ANL

2

slide-3
SLIDE 3
  • 1. Changing 2D à 3D distance for

additional muon hits

3/18/20

  • A. Rafique, ANL

3

slide-4
SLIDE 4

Additional selection cuts to increase hit purity

  • Added fewer additional hits to Michels from

the parent muons based on new cuts on some quantities

  • hit distance in coll. plane from muon end point
  • Cone angle
  • Was suggested to change 2Dà3D hit distance

2D hit distance in

  • coll. plane

from muon end point

3/18/20

  • A. Rafique, ANL

4

https://indico.fnal.gov/event/23362/contribution/1/material/slides/0.pdf cm

slide-5
SLIDE 5

2D à 3D distance for additional muon hits

2D hit distance in

  • coll. plane

from muon end point (Hit distance <=1) Applied an additional cut: (Hit distance <=1.5) Applied an additional cut: 3D hit distance in

  • coll. plane

from muon end point

3/18/20

  • A. Rafique, ANL

5

slide-6
SLIDE 6

Updated additional hits purity

false hits True hits True hits false hits

  • Hit purity total (59%)
  • Hit purity for events with
  • E<60 MeV: 61%
  • E>=60 MeV: 40%
  • Hit purity total (62%)
  • Hit purity for events with
  • E<60 MeV: 64%
  • E>=60 MeV: 58%

3/18/20

  • A. Rafique, ANL

6

slide-7
SLIDE 7

Updated total hit purity

false hits True hits false hits True hits

  • Hit purity (77%)
  • Hit purity for events with
  • E<60 MeV: 80%
  • E>=60 MeV: 55%
  • Hit purity (76%)
  • Hit purity for events with
  • E<60 MeV: 80%
  • E>=60 MeV: 45%

3/18/20

  • A. Rafique, ANL

7

slide-8
SLIDE 8

Hit completeness

  • Hit completeness (76%)
  • Hit completeness (76%)

True hits not selected True hits selected True hits not selected True hits selected

3/18/20

  • A. Rafique, ANL

8

slide-9
SLIDE 9

Michel energy spectrum in MC

3/18/20

  • A. Rafique, ANL

9

No noticeable change in the Michel energy distribution

slide-10
SLIDE 10

Michel energy data MC comparison

Michel reconstructed energy (MeV) Michel reconstructed energy (MeV)

3/18/20

  • A. Rafique, ANL

10

Mean ~ 31 MeV Mean ~ 30 MeV

slide-11
SLIDE 11
  • 2. Michel energy spectrum with

cone-only hits Data and MC comparison

3/18/20

  • A. Rafique, ANL

11

slide-12
SLIDE 12
  • 2. Michel energy spectrum cone-only hits

Data and MC comparison

Indicates that the selection doesn’t bias the distribution for MC and data

3/18/20

  • A. Rafique, ANL

12

slide-13
SLIDE 13
  • 3. Looking for kink around the

muon end position

3/18/20

  • A. Rafique, ANL

13

slide-14
SLIDE 14
  • 3. Looking for kink around the muon end

position

Linear correlation coefficient for five hits surrounding the Bragg’s peak MC

Most of the last hits of muons are linear, Pandora probably looks for it already

Data

3/18/20

  • A. Rafique, ANL

14

slide-15
SLIDE 15
  • 4. Recombination factor

Investigation

3/18/20

  • A. Rafique, ANL

15

slide-16
SLIDE 16
  • 4. Recombination factor Investigation
  • No direct way of getting a Michel-only generated sample using Gen
  • Recombination factors of a sample of electrons with energy 20, 40, and 50 MeV
  • Took the average value = 0.6441

Aaron H.: https://indico.fnal.gov/event/23307/contribution/2/material/slides/0.pdf (slide 24) E= 20 MeV E= 40 MeV E = 50 MeV

3/18/20

  • A. Rafique, ANL

16

slide-17
SLIDE 17

After changing recombination factor

Recomb factor = 0.7 Recomb factor = 0.6441

3/18/20

  • A. Rafique, ANL

17

slide-18
SLIDE 18
  • 5. Implementing new wire charge

integral method

3/18/20

  • A. Rafique, ANL

18

slide-19
SLIDE 19
  • 5. New wire charge integral method

Wire charge

Talk by Aaron: https://indico.fnal.gov/event/23307/contrib ution/2/material/slides/0.pdf MC

3/18/20

  • A. Rafique, ANL

19

slide-20
SLIDE 20

New wire charge integral method

New Old True

3/18/20

  • A. Rafique, ANL

20

On average, around 6 MeV energy is recovered/increased

slide-21
SLIDE 21

Data/MC comparison

3/18/20

  • A. Rafique, ANL

21

Mean ~ 32 MeV Mean ~ 38 MeV

slide-22
SLIDE 22

Summary for muon-based calibration

  • 1. Changed 2Dà3D distance for additional muon hits
  • No significant difference is observed
  • 2. Cone-only hits Data and MC comparison
  • Data and MC agree which indicates that the selection doesn’t bias the

distribution much

  • 3. Looked for kink around the muon end position
  • Found that most of the last muon hits are linear
  • 4. Recombination factor Investigation
  • On average, 3 MeV energy was retrieved with new recomb factor
  • 5. Implemented new wire charge integral method
  • On average, 6 MeV energy was retrieved with new new wire charge integral

method

3/18/20

  • A. Rafique, ANL

22

slide-23
SLIDE 23
  • 6. Calibrate to true Michel

distribution

Implemented by MiniBooNE

3/18/20

  • A. Rafique, ANL

23

slide-24
SLIDE 24

Calibrate to true Michel distribution

  • We have been using stopping muon calibrations until

this point

  • Ajib: https://docs.dunescience.org/cgi-

bin/private/RetrieveFile?docid=15974&filename=prod2_calibration_constants_for_select ed_runs.pdf&version=2

  • Now switching to true Michel calibration
  • Method implemented by MiniBooNE

(https://wwwboone.fnal.gov/publications/Papers/ccox_thesis.pdf)

  • MiniBooNE calibrates charge collected by PMTs to true

Michel energy by a calibration fit procedure.

  • Charge distribution is fit to true Michel distribution

smeared with a gaussian

3/18/20

  • A. Rafique, ANL

24

MiniBooNE

slide-25
SLIDE 25

Few examples to show effects of finite energy resolution

Added a gaussian distribution to true Michel energy distribution Smeared with a Gaussian: Mean = 0 MeV

  • St. Dev = 10 MeV

Smeared with a Gaussian: Mean = 0 MeV

  • St. Dev = 5 MeV

Initial true Michel energy distribution in ProtoDUNE

3/18/20

  • A. Rafique, ANL

25

slide-26
SLIDE 26

Method implemented for ProtoDUNE-SP TPC charge

  • To show how this method could in principle work with

ProtoDUNE TPC charge, we performed the following steps:

  • Started with candidate Michel reco charge (not energy)

distribution

  • Consider a true Michel energy distribution smeared with a

gaussian with some fit parameters

  • Match peaks of the two distributions
  • Calculated χ2 in a range of the distributions
  • Best χ2 value indicatess that we found the best match
  • Fit parameters tell us the energy resolution of the detector

3/18/20

  • A. Rafique, ANL

26

Michel candidate charge distribution MC

slide-27
SLIDE 27

Energy resolution χ2

  • True Michel energy smeared with gaussian

with mean = “0” and δ(Ε) = “2√E”

  • Resolution: δ(Ε)/E = 28% resolution at 50

MeV

  • Divided reco charge to match peaks by

factor = 115.714

  • Achieved with best value of χ2 in the range

(11--64) MeV

3/18/20

  • A. Rafique, ANL

27

slide-28
SLIDE 28

Data & MC Michel reco charge distribution comparison

MC data

data divided by an additional scale factor

!"#$_&"!'( &)*+,_&"!'( !"#" !"#$_&"!'( &)*+,_&"!'( $% = 0.917

Ajib: https://docs.dunescience.org/cgi-bin/private/RetrieveFile?docid=15974&filename=prod2_calibration_constants_for_selected_runs.pdf&version=2 3/18/20

  • A. Rafique, ANL

28

slide-29
SLIDE 29

Data and MC comparison

3/18/20

  • A. Rafique, ANL

29

slide-30
SLIDE 30

Summary and next steps

  • Tried several updates/techniques throughout the analysis
  • Significant difference is not observed in purity and energy distribution
  • We illustrated how the new method for Michel energy calibration,

independent from muon dE/dx, could work

  • will add a new section in the note describing this method and results

3/18/20

  • A. Rafique, ANL

30

slide-31
SLIDE 31

Backup slides

3/18/20

  • A. Rafique, ANL

31

slide-32
SLIDE 32

Energy resolution χ2

True Michel energy smearing with gaussian = 14% = 1√E Divided reco charge to match peaks by factor = 99.7537 χ2 in the range (11--64) MeV is = 1744/18=97

3/18/20

  • A. Rafique, ANL

32

slide-33
SLIDE 33

Energy resolution χ2

True Michel energy smearing with gaussian = 21% = 1.5√E Divided reco charge to match peaks by factor = 107.143 χ2 in the range (11--64) MeV is = 447/18=25

3/18/20

  • A. Rafique, ANL

33

slide-34
SLIDE 34

Energy resolution χ2

True Michel energy smearing with gaussian = 35% = 2.5√E Divided reco charge to match peaks by factor = 115.714 χ2 in the range (11--64) MeV is = 165/18=9.1

3/18/20

  • A. Rafique, ANL

34

slide-35
SLIDE 35

Data & MC charge distribution comparison

True Michel energy smearing with gaussian = 28% = 2√E Divided reco charge to match peaks by factor = 115.714 χ2 in the range (11--64) MeV is = 121/18=6.7 Applied the same scale factor to data (115.714) MC data

3/18/20

  • A. Rafique, ANL

35

slide-36
SLIDE 36

Data & MC charge distribution comparison

After scaling data with a factor equals

!"#$_&"!'( &)*+,_&"!'( !"#" !"#$_&"!'( &)*+,_&"!'( $% = 0.917

3/18/20

  • A. Rafique, ANL

36