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 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
1. Changing 2D à 3D distance for additional muon hits 3/18/20 A. Rafique, ANL 3
Additional selection cuts to increase hit purity • Added fewer additional hits to Michels from the parent muons based on new cuts on some quantities 2D hit distance in • hit distance in coll. plane from muon end point coll. plane • Cone angle from muon end point • Was suggested to change 2D à 3D hit distance https://indico.fnal.gov/event/23362/contribution/1/material/slides/0.pdf cm 3/18/20 A. Rafique, ANL 4
2D à 3D distance for additional muon hits 3D hit 2D hit distance in distance in coll. plane coll. plane from muon from muon end point end point Applied an additional cut: Applied an additional cut: (Hit distance <=1.5) (Hit distance <=1) 3/18/20 A. Rafique, ANL 5
Updated additional hits purity • Hit purity total (59%) • Hit purity total (62%) • Hit purity for events with • Hit purity for events with • E<60 MeV: 61% • E<60 MeV: 64% • E>=60 MeV: 40% • E>=60 MeV: 58% false True false True hits hits hits hits 3/18/20 A. Rafique, ANL 6
Updated total hit purity • Hit purity (76%) • Hit purity (77%) • Hit purity for events with • Hit purity for events with • E<60 MeV: 80% • E<60 MeV: 80% • E>=60 MeV: 45% • E>=60 MeV: 55% True false True false hits hits hits hits 3/18/20 A. Rafique, ANL 7
Hit completeness • Hit completeness (76%) • Hit completeness (76%) True hits True hits True hits True hits selected selected not not 3/18/20 selected A. Rafique, ANL 8 selected
Michel energy spectrum in MC No noticeable change in the Michel energy distribution 3/18/20 A. Rafique, ANL 9
Michel energy data MC comparison Mean ~ 31 MeV Mean ~ 30 MeV Michel reconstructed energy (MeV) Michel reconstructed energy (MeV) 3/18/20 A. Rafique, ANL 10
2. Michel energy spectrum with cone-only hits Data and MC comparison 3/18/20 A. Rafique, ANL 11
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
3. Looking for kink around the muon end position 3/18/20 A. Rafique, ANL 13
3. Looking for kink around the muon end position Linear correlation coefficient for five hits surrounding the Bragg’s peak MC Data Most of the last hits of muons are linear, Pandora probably looks for it already 3/18/20 A. Rafique, ANL 14
4. Recombination factor Investigation 3/18/20 A. Rafique, ANL 15
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 • E= 20 E = 50 E= 40 MeV MeV MeV Aaron H.: https://indico.fnal.gov/event/23307/contribution/2/material/slides/0.pdf (slide 24) 3/18/20 A. Rafique, ANL 16
After changing recombination factor Recomb factor = 0.6441 Recomb factor = 0.7 3/18/20 A. Rafique, ANL 17
5. Implementing new wire charge integral method 3/18/20 A. Rafique, ANL 18
5. New wire charge integral method Wire charge MC Talk by Aaron: https://indico.fnal.gov/event/23307/contrib ution/2/material/slides/0.pdf 3/18/20 A. Rafique, ANL 19
New wire charge integral method On average, around 6 MeV energy is recovered/increased New Old True 3/18/20 A. Rafique, ANL 20
Data/MC comparison Mean ~ 32 MeV Mean ~ 38 MeV 3/18/20 A. Rafique, ANL 21
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
6. Calibrate to true Michel distribution Implemented by MiniBooNE 3/18/20 A. Rafique, ANL 23
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 MiniBooNE • 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
Few examples to show effects of finite energy resolution Added a gaussian distribution to true Michel energy distribution Smeared with a Smeared with a Gaussian: Initial true Michel energy Gaussian: Mean = 0 MeV Mean = 0 MeV distribution in ProtoDUNE St. Dev = 10 MeV St. Dev = 5 MeV 3/18/20 A. Rafique, ANL 25
Method implemented for ProtoDUNE-SP TPC charge • To show how this method could in principle work with Michel candidate ProtoDUNE TPC charge, we performed the following charge distribution steps: • Started with candidate Michel reco charge (not energy) MC 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
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
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
Data and MC comparison 3/18/20 A. Rafique, ANL 29
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
Backup slides 3/18/20 A. Rafique, ANL 31
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
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
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
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 MC Applied the same scale factor to data (115.714) data 3/18/20 A. Rafique, ANL 35
Data & MC charge distribution comparison !"#$_&"!'( &)*+,_&"!'( !"#" After scaling data with a factor equals &)*+,_&"!'( $% = 0.917 !"#$_&"!'( 3/18/20 A. Rafique, ANL 36
Recommend
More recommend
