Reconstruction of Low Energy Neutrino Events with GPUs at IceCube

IceCube is a cubic kilometer neutrino observatory located at the South Pole that produces massive amounts of data by measuring individual Cherenkov photons from neutrino interaction events in the energy range from few GeV to several PeV. The actual reconstruction of neutrino events in the GeV range is computationally challenging due to the scarcity of data produced by single events. This can lead to run times of several weeks for the state-of-the-art reconstruction method – Pegleg – on CPUs for typical workloads of many ten-thousand events. We propose a GPU version of Pegleg that probes the likelihood space with several hypotheses in parallel while adapting the amount of parallel sampled hypotheses dynamically in order to reduce computation time significantly. Our results show an average speedup of 14 (with a maximum of over 200) for 5262 reconstructed neutrino events of different flavors on a Titan V GPU compared to the multithreaded CPU version, which enables quicker and broader analysis of IceCube events.

[1]  P. O. Hulth,et al.  Measurement of atmospheric neutrino oscillations with IceCube. , 2013, Physical review letters.

[2]  T. B. Watson,et al.  Development of an analysis to probe the neutrino mass ordering with atmospheric neutrinos using three years of IceCube DeepCore data , 2019, The European Physical Journal C.

[3]  M. P. Hobson,et al.  Importance Nested Sampling and the MultiNest Algorithm , 2013, The Open Journal of Astrophysics.

[4]  Marco Maggioni,et al.  Dissecting the NVIDIA Volta GPU Architecture via Microbenchmarking , 2018, ArXiv.

[5]  Tjalling J. Ypma,et al.  Historical Development of the Newton-Raphson Method , 1995, SIAM Rev..

[6]  A. C. Kaboth,et al.  Accelerated event-by-event neutrino oscillation reweighting with matter effects on a GPU , 2013, 1311.7579.

[7]  Richard J. Cook,et al.  The Statistical Analysis of Recurrent Events , 2007 .

[8]  A. Schukraft,et al.  The IceCube Neutrino Observatory: Instrumentation and Online Systems , 2016, 1612.05093.

[9]  C. Giunti,et al.  Fundamentals of Neutrino Physics and Astrophysics , 2007 .

[10]  P. O. Hulth,et al.  Energy reconstruction methods in the IceCube neutrino telescope , 2013, 1311.4767.

[11]  T. B. Watson,et al.  Measurement of Atmospheric tau Neutrino Appearance With Icecube/deepcore , 2019, 1901.05366.

[12]  M. Leuermann,et al.  Testing the neutrino mass ordering with IceCube DeepCore , 2019 .

[13]  Carl de Boor,et al.  On “best” interpolation☆ , 1976 .

[14]  Bertil Schmidt,et al.  Massively parallel computation of atmospheric neutrino oscillations on CUDA-enabled accelerators , 2019, Comput. Phys. Commun..

[15]  G. H. Collin,et al.  Computational techniques for the analysis of small signals in high-statistics neutrino oscillation experiments , 2018 .

[16]  F. Feroz,et al.  MultiNest: an efficient and robust Bayesian inference tool for cosmology and particle physics , 2008, 0809.3437.

[17]  F. Feroz,et al.  Multimodal nested sampling: an efficient and robust alternative to Markov Chain Monte Carlo methods for astronomical data analyses , 2007, 0704.3704.