michel electron reconstruc0on in dune
play

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

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


  1. Michel Electron Reconstruc0on in DUNE Aleena Rafique , Zelimir Djurcic on behalf of the DUNE collabora7on CPAD 2019 12/10/2019

  2. Outline • Why study Michel electrons • TPC Michel Reconstruc7on • Energy Reconstruc7on and Michel Energy Spectrum • Conclusions 12/10/19 A. Rafique, ANL 2

  3. Michel electrons • 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 A. Rafique, ANL 3

  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

  5. ProtoDUNE-SP LArTPC Y= beam X= Z= 12/10/19 A. Rafique, ANL 5

  6. Candidate michel event displays in ProtoDUNE Candidate Candidate Candidate muon muon michel Candidate Other michel cosmics 12/10/19 6 A. Rafique, ANL

  7. Michel event ini0al selec0on 1. T 0 tagged (cathode crossing) tracks . beam y z x (0,0,0) 12/10/19 A. Rafique, ANL 7

  8. Michel event ini0al selec0on 1. T 0 tagged (cathode crossing) tracks 30 cm 30 cm 30 cm 2. Tracks star7ng in detector edges 12/10/19 A. Rafique, ANL 8

  9. Michel event ini0al selec0on 1. T 0 tagged (cathode 3. Tracks crossing) ending in FV tracks 2. Tracks star7ng in detector edges 12/10/19 A. Rafique, ANL 9

  10. Michel event ini0al selec0on 1. T 0 tagged (cathode 3. Tracks crossing) e e e e ending in FV tracks e e e e e e Drig direc7on 4. Tracks 2. Tracks ending not star7ng in around APA detector bounds edges 12/10/19 A. Rafique, ANL 10

  11. Candidate muon selec0on • Remove broken tracks Candidate muon θ d MC Sample purity = 27% 12/10/19 A. Rafique, ANL 11

  12. Candidate muon selec0on • Remove broken tracks True MC normalized to reco MC • Select candidate muon tracks with > 75 cm length ProtoDUNE-SP Simulation ProtoDUNE-SP Simulation, Preliminary reco events 4500 reco events 4500 true michel events true michel events 4000 4000 ProtoDUNE simula7on 3500 3500 3000 3000 # Events # Events 2500 2500 2000 2000 1500 1500 1000 1000 length > 75 cm 500 500 0 0 0 100 200 300 400 500 600 700 800 900 1000 0 100 200 300 400 500 600 700 800 900 1000 candidate muon track length (cm) candidate muon track length (cm) MC Sample purity = 27% 12/10/19 A. Rafique, ANL 12

  13. Candidate muon selec0on • Remove broken tracks • Select candidate muon tracks with > 75 cm length ProtoDUNE-SP Simulation ProtoDUNE-SP Simulation ProtoDUNE-SP Simulation, Preliminary 5 5 10 10 reco events reco events reco events • Minimum hit 7me > 200 7cks 4500 true michel events true michel events true michel events 4000 ProtoDUNE 4 4 10 10 simula7on 3500 3 3000 3 10 10 # Events # Events # Events ProtoDUNE 2500 ProtoDUNE simula7on 2 10 2 10 2000 simula7on 1500 10 10 1000 500 1 1 e e e 0 0 1000 2000 3000 4000 5000 6000 0 0 100 1000 200 300 2000 400 3000 500 600 4000 700 800 5000 900 1000 6000 e e Minimum hit peak time Minimum hit peak time candidate muon track length (cm) e MC Sample purity = 63% 12/10/19 A. Rafique, ANL 13

  14. Candidate michel selec0on • Remove broken tracks • Select candidate muon tracks with > 75 cm length ProtoDUNE-SP Simulation ProtoDUNE-SP Simulation, Preliminary reco events reco events • Minimum hit 7me > 200 7cks 4 10 4 10 true michel events true michel events ProtoDUNE • Nearby hit count > 5 & < 40 simula7on 3 10 3 10 # Events # Events ProtoDUNE 2 10 2 10 ProtoDUNE simula7on simula7on Drig dir (7me) 10 10 R = 1 1 1 0 c m 0 20 40 60 80 100 120 140 160 180 200 0 20 40 60 80 100 120 140 160 180 200 number of hits close to the track end point number of hits close to the track end point Wire (beam direc7on) MC Sample purity = 86% 12/10/19 A. Rafique, ANL 14

  15. Candidate michel selec0on • Remove broken tracks • Select candidate muon tracks with > 75 cm length ProtoDUNE-SP Simulation ProtoDUNE-SP Simulation, Preliminary • Minimum hit 7me > 200 7cks reco events reco events 7000 7000 ProtoDUNE true michel events true michel events • Nearby hit count > 5 & < 40 simula7on 6000 6000 5000 5000 • Closest reco shower distance ProtoDUNE # Events # Events ProtoDUNE ProtoDUNE simula7on 4000 4000 simula7on < 10 cm simula7on 3000 3000 2000 2000 1000 1000 R = 1 0 0 0 0 0 5 5 10 10 15 15 20 20 25 25 30 30 c m shower start distance from muon end point (cm) shower start distance from muon end point (cm) Final MC Sample purity = 89% 12/10/19 A. Rafique, ANL 15

  16. Michel energy reconstruc0on 9 , ∗; < ∗; => ∗?,/4∗@/AB_D3.7/A 𝐹 = ∑ ,-./00. 20345 6,78 ;30,E_./487∗F5./BE_D3.7/A Where Q i = charge deposited on a hit C x ,C yz = 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 ProtoDUNE-SP Simulation, Preliminary Calib_const = converts dQ/dx in ADC/cm into dE/dx in MeV/cm 2000 Recomb_factor = 0.7; to incorporate the recombina7on effects 1800 true_michel_mc 1600 Entries 12028 1400 Mean 33.58 Std Dev 14.68 # Events 1200 1000 800 600 400 200 0 0 20 40 60 80 100 120 140 160 180 200 True michel reconstructed energy (MeV) 12/10/19 A. Rafique, ANL 16

  17. Michel energy reconstruc0on True michel reco hits • Defined a cone at the end point ProtoDUNE-SP Simulation, Preliminary 50 of the parent muon 10000 40 30 8000 cone length (cm) 20 Candidate 10 6000 muon 0 - 10 4000 - 20 - 30 2000 L - 40 θ - 50 0 0 20 40 60 80 100 120 140 160 180 ProtoDUNE-SP Simulation, Preliminary cone angle 2000 Other hits of the event Missing reco MC 1800 michel true michel MC ProtoDUNE-SP Simulation, Preliminary 3 1600 ´ 10 100 Where θ=0 hits 1400 reco_mc 80 corresponds to the Entries 12028 120 # Events 1200 Mean 31.09 60 hits along the Std Dev 21.43 1000 100 40 cone length (cm) nearby wire hits 800 true_michel_mc 20 80 Entries 12028 600 0 Mean 33.58 Std Dev 14.68 60 - 400 20 - 200 40 40 - 0 60 0 20 40 60 80 100 120 140 160 180 200 cone 20 Reconstructed energy (MeV) - 80 - 100 0 0 20 40 60 80 100 120 140 160 180 cone angle 12/10/19 A. Rafique, ANL 17

  18. Recovering michel hits from parent muon Before recovery • Look at the charge deposi7on at the ProtoDUNE-SP Simulation, Preliminary 1 last 10 hits of reco muons 900 0.8 800 0.6 • Charge dep = (Q i – Q i-1 )/Q i-1 700 Charge deposition 0.4 600 0.2 • Ager removing hits beyond the 500 0 400 - 0.2 maximum truncated charge value - 300 0.4 - • Recovered 6% of the total missing 200 0.6 - 0.8 100 michel hits - 1 0 0 2 4 6 8 10 Hit number Ager recovery True muon ProtoDUNE-SP Simulation, Preliminary 1 ProtoDUNE-SP Simulation, Preliminary 1 900 0.8 900 0.8 800 0.6 800 0.6 700 Charge deposition 0.4 Charge deposition 700 0.4 600 0.2 600 0.2 500 0 500 0 - 400 0.2 - 400 0.2 - 0.4 300 - 0.4 300 - - 0.6 200 0.6 200 - - 0.8 0.8 100 100 - - 1 0 1 0 0 2 4 6 8 10 0 2 4 6 8 10 Hit number Hit number 12/10/19 A. Rafique, ANL 18 Last hit

  19. Michel energy spectrum MC ProtoDUNE-SP Simulation, Preliminary 100 120 90 Michel reco energy (MeV) 80 100 70 ProtoDUNE-SP Simulation, Preliminary 80 60 2000 reco MC 50 60 1800 40 true michel MC 1600 30 40 20 reco_mc 1400 20 10 Entries 12028 # Events 1200 Mean 36.8 0 0 0 10 20 30 40 50 60 70 80 90 100 Std Dev 22.33 True michel reco energy (MeV) 1000 true_michel_mc 800 ProtoDUNE-SP Simulation, Preliminary micheltrueparEreso Entries 12028 3000 Entries 12028 Mean 33.58 600 - Mean 3.431 Std Dev 14.68 Std Dev 23.55 2500 400 200 2000 # Events 0 0 20 40 60 80 100 120 140 160 180 200 1500 Reconstructed energy (MeV) 1000 500 0 - - - - 200 150 100 50 0 50 100 150 200 True michel reco - reco energy (MeV) 12/10/19 A. Rafique, ANL 19

  20. Michel energy data/MC comparison ProtoDUNE-SP, Preliminary ProtoDUNE-SP, Preliminary 1200 1200 MC MC Bkg MC Data 1000 1000 Data MC_copy 800 MC_copy Very good 800 Entries 13509 Entries 13509 agreement in Entries Entries Mean 38.19 Mean 38.19 Std Dev 24.3 ProtoDUNE data Std Dev 24.3 600 600 and simula7on Data Data Data Data Entries Entries 7760 7760 400 400 Entries Entries 7760 7760 Mean Mean 38.89 38.89 Mean Mean 38.89 38.89 Std Dev 24.52 Std Dev 24.52 Std Dev 24.52 Std Dev 24.52 200 200 0 0 0 20 40 60 80 100 120 140 160 180 200 0 20 40 60 80 100 120 140 160 180 200 Reconstructed energy (MeV) Reconstructed energy (MeV) 12/10/19 A. Rafique, ANL 20

Recommend


More recommend