Highly-parallelized simulation of a pixelated LArTPC on a GPU

Abstract

The rapid development of general-purpose computing on graphics processing units (GPGPU) is allowing the implementation of highly-parallelized Monte Carlo simulation chains for particle physics experiments. This technique is particularly suitable for the simulation of a pixelated charge readout for time projection chambers, given the large number of channels that this technology employs. Here we present the first implementation of a full microphysical simulator of a liquid argon time projection chamber (LArTPC) equipped with light readout and pixelated charge readout, developed for the DUNE Near Detector. The software is implemented with an end-to-end set of GPU-optimized algorithms. The algorithms have been written in Python and translated into CUDA kernels using Numba, a just-in-time compiler for a subset of Python and NumPy instructions. The GPU implementation achieves a speed up of four orders of magnitude compared with the equivalent CPU version. The simulation of the current induced on $10^3$ pixels takes around 1 ms on the GPU, compared with approximately 10 s on the CPU. The results of the simulation are compared against data from a pixel-readout LArTPC prototype.

arXiv comment26 pages, 15 figures

References in corpus (10)

• Array Programming with NumPy
• NEST: A Comprehensive Model for Scintillation Yield in Liquid Xenon
• A study of electron recombination using highly ionizing particles in the ArgoNeuT Liquid Argon TPC
• Oxygen contamination in liquid Argon: combined effects on ionization electron charge and scintillation light
• Measurement of Longitudinal Electron Diffusion in Liquid Argon
• LArPix: Demonstration of low-power 3D pixelated charge readout for liquid argon time projection chambers
• Effects of Nitrogen contamination in liquid Argon
• Properties of Liquid Argon Scintillation Light Emission
• Enhancing Neutrino Event Reconstruction with Pixel-Based 3D Readout for Liquid Argon Time Projection Chambers
• Optimizing the Hit Finding Algorithm for Liquid Argon TPC Neutrino Detectors Using Parallel Architectures

Cited by in corpus (2)

• Performance of an Optical TPC Geant4 Simulation with Opticks GPU-Accelerated Photon Propagation
• Contrastive Learning for Robust Representations of Neutrino Data