Jake Calcutt Dec. 19, 2019
Jake Calcutt 2 The ProtoDUNE beamline provides a momentum measurement to be used within analyses Need to investigate resolution and bias on the measurement within MC Can give us an estimate on this within data ●
Jake Calcutt 3 Bending Magnets Scintillator Profile Monitors For Momentum/Tracking Scintillator Planes for TOF and Trigger Cerenkov Devices for Particle ID (PID)
Jake Calcutt 4 Project hit in Profiler 1 to magnet, connect hits in Profilers 2,3 ● Get θ between these trajectories ○ Get current in magnet: I B ● Momentum: p ~ f( I B ) / θ ● θ Profiler 1 Magnet Profiler 2 Profiler 3
Jake Calcutt 5 At last week’s Analysis meeting, Flavio brought up the simulated energy loss within the beamline This motivated me to look into: 1) At what point does the initial momentum in MC come from? a) How much energy is ‘already lost’ at generation time? 2) How does the true momentum compare to the reconstructed momentum (in MC) Note: using simulation files from beam group listed here
Jake Calcutt 6 Looking at dunetpc/dune/EventGenerator/ProtoDUNEBeam_module.cc The momentum comes from that as simulated at the front of the TPC
Jake Calcutt 7 Our MC events are produced at a point where energy loss has already been considered Motivates us to consider the difference between the momentum at this point and as measured in the spectrometer
Jake Calcutt 8 An additional effect to consider is the inherent resolution (+ bias) on the spectrometer measurement as compared to the true momentum at the spectrometer This effect is ‘folded in’ with the energy loss See backup for these plots alone
Jake Calcutt 9 Proton Pion μ: 6.50 μ: 2.09 σ : 23.83 σ : 21.24 Electron μ: 6.00 Source of data-MC σ : 22.54 discrepancies in Full Avg: 44.7 Aaron’s studies See full range of distribution in backup Note: p = Reconstructed - True (@ NP04 Front)
Jake Calcutt 10 Pion Proton μ: -2.62 μ: -4.58 σ : 52.67 σ : 52.26 Electron μ: 9.49 σ : 53.60 Full Avg: 299 Note: p = Reconstructed - True (@ NP04 Front)
Jake Calcutt 11 Proton Pion μ: -8.97 μ: -6.37 σ : 102.07 σ : 101.80 Electron μ: 41.19 σ : 112.44 Full Avg: 1031 Note: p = Reconstructed - True (@ NP04 Front)
Jake Calcutt 12 Have an estimate in MC for smearing from true to reco momentum within the beamline Should give more accurate comparisons between MC and data in ● various analyses (i.e. Aaron’s electron analysis) We need to consider how well the simulation describes our beamline Might need to take into account model uncertainties ● With access + understanding of the base simulation files, I’m going to implement adding the ‘ProtoDUNEBeamEvent’ class to MC within the event generator
Jake Calcutt 13
Jake Calcutt 14 Proton Pion μ: -3.40 μ: -3.06 σ : 23.66 σ : 21.23 Electron μ: -3.89 σ : 21.29
Jake Calcutt 15 Pion Proton μ: -12.03 μ: -12.75 σ : 52.29 σ : 51.63 Electron μ: -13.66 σ : 52.07
Jake Calcutt 16 Proton Pion μ: 23.75 μ: -23.89 σ : 101.80 σ : 101.86 Electron μ: -25.87 σ : 101.26
Jake Calcutt 17 Electron Electron 1 GeV 3 GeV Electron 6 GeV
Jake Calcutt Thanks for listening
Recommend
More recommend