DUNE Collaboration(Abud, A. A. et al), Amedo, P., Antonova, M., Barenboim, G., Cervera-Villanueva, A., De Romeri, V., et al. (2023). Highly-parallelized simulation of a pixelated LArTPC on a GPU. J. Instrum., 18(4), P04034–35pp.
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 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.
|
DUNE Collaboration(Abud, A. A. et al), Amedo, P., Antonova, M., Barenboim, G., Cervera-Villanueva, A., De Romeri, V., et al. (2023). Identification and reconstruction of low-energy electrons in the ProtoDUNE-SP detector. Phys. Rev. D, 107(9), 092012–22pp.
Abstract: 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.
|
DUNE Collaboration(Abud, A. A. et al), Amedo, P., Antonova, M., Barenboim, G., Benitez Montiel, C., Cervera-Villanueva, A., et al. (2023). Impact of cross-section uncertainties on supernova neutrino spectral parameter fitting in the Deep Underground Neutrino Experiment. Phys. Rev. D, 107(11), 112012–25pp.
Abstract: A primary goal of the upcoming Deep Underground Neutrino Experiment (DUNE) is to measure the Oo10 thorn 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;oE & nu; thorn 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;oE & nu; thorn modeling uncertainties on DUNE's supernova neutrino physics sensitivity for the first time. We find that the currently large theoretical uncertainties on & sigma;oE & nu; thorn 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;oE & nu; thorn 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;oE & nu; thorn . A direct measurement of low-energy & nu;e-argon scattering would be invaluable for improving the theoretical precision to the needed level.
|
Renner, J., Cervera-Villanueva, A., Hernando, J. A., Izmaylov, A., Monrabal, F., Muñoz, J., et al. (2015). Improved background rejection in neutrinoless double beta decay experiments using a magnetic field in a high pressure xenon TPC. J. Instrum., 10, P12020–19pp.
Abstract: We demonstrate that the application of an external magnetic field could lead to an improved background rejection in neutrinoless double-beta (0 nu beta beta) decay experiments using a high-pressure xenon (HPXe) TPC. HPXe chambers are capable of imaging electron tracks, a feature that enhances the separation between signal events (the two electrons emitted in the 0 nu beta beta decay of Xe-136) and background events, arising chiefly from single electrons of kinetic energy compatible with the end-point of the 0 nu beta beta decay (Q(beta beta)). Applying an external magnetic field of sufficiently high intensity (in the range of 0.5-1 Tesla for operating pressures in the range of 5-15 atmospheres) causes the electrons to produce helical tracks. Assuming the tracks can be properly reconstructed, the sign of the curvature can be determined at several points along these tracks, and such information can be used to separate signal (0 nu beta beta) events containing two electrons producing a track with two different directions of curvature from background (single-electron) events producing a track that should spiral in a single direction. Due to electron multiple scattering, this strategy is not perfectly efficient on an event-by-event basis, but a statistical estimator can be constructed which can be used to reject background events by one order of magnitude at a moderate cost (about 30%) in signal efficiency. Combining this estimator with the excellent energy resolution and topological signature identification characteristic of the HPXe TPC, it is possible to reach a background rate of less than one count per ton-year of exposure. Such a low background rate is an essential feature of the next generation of 0 nu beta beta experiments, aiming to fully explore the inverse hierarchy of neutrino masses.
|
T2K Collaboration(Abe, K. et al), Antonova, M., Cervera-Villanueva, A., & Molina Bueno, L. (2021). Improved constraints on neutrino mixing from the T2K experiment with 3.13 x 10(21) protons on target. Phys. Rev. D, 103(11), 112008–59pp.
Abstract: The T2K experiment reports updated measurements of neutrino and antineutrino oscillations using both appearance and disappearance channels. This result comes from an exposure of 14.9(16.4) x 10(20) protons on target in neutrino (antineutrino) mode. Significant improvements have been made to the neutrino interaction model and far detector reconstruction. An extensive set of simulated data studies have also been performed to quantify the effect interaction model uncertainties have on the T2K oscillation parameter sensitivity. T2K performs multiple oscillation analyses that present both frequentist and Bayesian intervals for the Pontecorvo-Maki-Nakagawa-Sakata parameters. For fits including a constraint on sin(2)theta(13) from reactor data and assuming normal mass ordering T2K measures sin(2)theta(13) = 0.53(-0.04)(+0.03) and Delta m(32)(2) = (2.45 +/- 0.07) x 10(-3) eV(2) c(-4). The Bayesian analyses show a weak preference for normal mass ordering 89)% posterior probability) and the upper sin(2)theta(13) octant (80% posterior probability), with a uniform prior probability assumed in both cases. The T2K data exclude CP conservation in neutrino oscillations at the 2 sigma level.
|