Highly-parallelized simulation of a pixelated LArTPC on a GPU

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 103 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

Identification and reconstruction of low-energy electrons in the ProtoDUNE-SP detector

Measurements of electrons from νe interactions are crucial for the Deep Underground Neutrino Experiment (DUNE) neutrino oscillation program, as well as searches for physics beyond the standard model, supernova neutrino detection, and solar neutrino measurements. This article describes the selection and reconstruction of low-energy (Michel) electrons in the ProtoDUNE-SP detector. ProtoDUNE-SP is one of the prototypes for the DUNE far detector, built and operated at CERN as a charged particle test beam experiment. A sample of low-energy electrons produced by the decay of cosmic muons is selected with a purity of 95%. This sample is used to calibrate the low-energy electron energy scale with two techniques. An electron energy calibration based on a cosmic ray muon sample uses calibration constants derived from measured and simulated cosmic ray muon events. Another calibration technique makes use of the theoretically well-understood Michel electron energy spectrum to convert reconstructed charge to electron energy. In addition, the effects of detector response to low-energy electron energy scale and its resolution including readout electronics threshold effects are quantified. Finally, the relation between the theoretical and reconstructed low-energy electron energy spectra is derived, and the energy resolution is characterized. The low-energy electron selection presented here accounts for about 75% of the total electron deposited energy. After the addition of lost energy using a Monte Carlo simulation, the energy resolution improves from about 40% to 25% at 50 MeV. These results are used to validate the expected capabilities of the DUNE far detector to reconstruct low-energy electrons

Identification and reconstruction of low-energy electrons in the ProtoDUNE-SP detector

Measurements of electrons from $\nu_e$ interactions are crucial for the Deep Underground Neutrino Experiment (DUNE) neutrino oscillation program, as well as searches for physics beyond the standard model, supernova neutrino detection, and solar neutrino measurements. This article describes the selection and reconstruction of low-energy (Michel) electrons in the ProtoDUNE-SP detector. ProtoDUNE-SP is one of the prototypes for the DUNE far detector, built and operated at CERN as a charged particle test beam experiment. A sample of low-energy electrons produced by the decay of cosmic muons is selected with a purity of 95%. This sample is used to calibrate the low-energy electron energy scale with two techniques. An electron energy calibration based on a cosmic ray muon sample uses calibration constants derived from measured and simulated cosmic ray muon events. Another calibration technique makes use of the theoretically well-understood Michel electron energy spectrum to convert reconstructed charge to electron energy. In addition, the effects of detector response to low-energy electron energy scale and its resolution including readout electronics threshold effects are quantified. Finally, the relation between the theoretical and reconstructed low-energy electron energy spectrum is derived and the energy resolution is characterized. The low-energy electron selection presented here accounts for about 75% of the total electron deposited energy. After the addition of lost energy using a Monte Carlo simulation, the energy resolution improves from about 40% to 25% at 50~MeV. These results are used to validate the expected capabilities of the DUNE far detector to reconstruct low-energy electrons.Comment: 19 pages, 10 figure

Impact of cross-section uncertainties on supernova neutrino spectral parameter fitting in the Deep Underground Neutrino Experiment

A primary goal of the upcoming Deep Underground Neutrino Experiment (DUNE) is to measure the $\mathcal{O}(10)$ MeV neutrinos produced by a Galactic core-collapse supernova if one should occur during the lifetime of the experiment. The liquid-argon-based detectors planned for DUNE are expected to be uniquely sensitive to the $\nu_e$ component of the supernova flux, enabling a wide variety of physics and astrophysics measurements. A key requirement for a correct interpretation of these measurements is a good understanding of the energy-dependent total cross section $\sigma(E_\nu)$ for charged-current $\nu_e$ absorption on argon. In the context of a simulated extraction of supernova $\nu_e$ spectral parameters from a toy analysis, we investigate the impact of $\sigma(E_\nu)$ modeling uncertainties on DUNE's supernova neutrino physics sensitivity for the first time. We find that the currently large theoretical uncertainties on $\sigma(E_\nu)$ must be substantially reduced before the $\nu_e$ flux parameters can be extracted reliably: in the absence of external constraints, a measurement of the integrated neutrino luminosity with less than 10\% bias with DUNE requires $\sigma(E_\nu)$ to be known to about 5%. The neutrino spectral shape parameters can be known to better than 10% for a 20% uncertainty on the cross-section scale, although they will be sensitive to uncertainties on the shape of $\sigma(E_\nu)$. A direct measurement of low-energy $\nu_e$-argon scattering would be invaluable for improving the theoretical precision to the needed level.Comment: 25 pages, 21 figure

Highly-parallelized simulation of a pixelated LArTPC on a GPU

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

Proof of Concept of CLIC Final Focus Quadrupoles Stabilization

International audienceThe Compact LInear Collider (CLIC) [1] luminosity requires extremely low beam emittances. Therefore, high beam position stability is needed to provide cen-tral collisions of the opposing bunches. Since ground motion (GM) amplitudes are likely to be larger than the required tolerances, an Active Vibration Control (AVC) system is required to damp quadrupole motion to the desired value of 0.2 nm RMS at 4 Hz. This paper focuses on the vertical final focus quadrupoles (QD0, QF1) stabilization and demonstrates its feasibility. An AVC system to be installed under QD0 and QF1 has been developed and successfully tested at LAPP. Based on a dedicated homemade sensor with an ex-tremely low internal noise level of 0.05 nm at 4 Hz, it damps GM in the frequency range [3;70] Hz by up to 30 dB, leading to RMS values of approximately 0.25 nm at 4 Hz. Simulations based on GM measured in the Compact Muon Solenoid (CMS) experimental hall [2] show that with such a GM level, the specifications would only be achieved with a Passive Insulation (PI) system, which would filter ground motion starting at ~ 25 H

Vibration Control Using a Dedicated Inertial Sensor

International audienceThe Compact Linear Collider (CLIC) project is right in the development phase. In this prospect, the main objective of the final focus stabilization is to succeed in damping ground motion vibration at the sub-nanometer scale to avoid any unwanted motion of the two last final focus magnets. Active vibration control using high-resolution industrial seismic sensors has already shown its limits in the desired range of attenuation [4-100] Hz. Hence, a dedicated inertial sensor has been designed, tuned to fit CLIC requirements. It doesn't contain any feedback or coil which would linearize its dynamic in its bandwidth [0.1-100] Hz, avoiding electrical noise or thermal noise, and when used for active control, enhancing its performances thanks to the internal resonance. An analytical modeling of the sensor dynamic behavior also shown experimentally gives a first internal resonance around 11 Hz. This paper is complemented by a comparative measurement with high precision industrial sensors. A noise level 2.5 times better (0.04 nm at 4 Hz) has been achieved. This sensor has also been tested in the context of vibration rejection with a sub-nanometer active vibration rejection support. Experimental results have been compared with the one obtained with the state-of-the-art industrial sensors. Vertical seismic motion attenuation results have shown unprecedented performances. A damping ratio of 8.5 has been achieved at 4 Hz, leading to an rms displacement of the support of 0.25 nm thanks to the active support prototype and the dedicated sensor