Abstract:We introduce DarkNews , a lightweight Python -based Monte-Carlo generator for beyond-the-Standard-Model neutrino-nucleus scattering. The generator handles the production and decay of heavy neutral leptons via additional vector or scalar mediators, as well as through transition magnetic moments. DarkNews samples pre-computed neutrino-nucleus upscattering cross sections and heavy neutrino decay rates to produce dilepton and single-photon events in accelerator neutrino experiments. We present two case studies with differential distributions for models that can explain the MiniBooNE excess. The aim of this code is to aid the neutrino theory and experimental communities in performing searches and sensitivity studies for new particles produced in neutrino upscattering. Program summary Program Title: DarkNews CPC Library link to program files: https://doi.org/10.17632/k9nh5fn4pw.1 Developer's repository link:: https://github.com/mhostert/DarkNews-generator Licensing provisions: MIT Programming language: Python Nature of problem: In many theories beyond the Standard Model of particle physics, new light particles are introduced with masses below the 100 GeV scale. These can be produced by neutrino-nucleus interactions inside neutrino detectors, leaving a visible signature through their decay. To search for these signatures, experiments have to simulate the underlying physical processes, which requires detailed knowledge of the cross sections and decay rates in the model, as well as the ability to generate Monte Carlo events. Solution method: DarkNews provides the necessary tools to simulate the production of heavy neutral leptons in accelerator neutrino experiments, focusing on neutrino-nucleus scattering processes. The program's scope includes models with new mediators and three heavy neutral leptons, parameterizing a large class of theories through generic neutrino, quark, and lepton interactions. It generates weighted Monte-Carlo events using vegas and its pre-calculated differential rates.

Atmospheric Lepton Fluxes via Two-Dimensional Matrix Cascade Equations

Calculation of the Flux and Charge Ratio of High-Energy Atmospheric Muons at Sea Level*

Neutrinos from muon-rich ultra high energy electromagnetic cascades: The MUNHECA code

Neutrino initiated cascades at mid and high altitudes in the atmosphere

Efficient Modeling of Heavy Cosmic Rays Propagation in Evolving Astrophysical Environments

A Non-Parametric Approach to Infer the Energy Spectrum and the Mass Composition of Cosmic Rays

Novel Simulation Framework for Analyzing Cosmic Ray Particle Distributions at a Global Scale

Lepton-Nucleus Interactions within Microscopic Approaches

Neutrino Transport with the Monte Carlo Method. II. Quantum Kinetic Equations

QED radiative corrections for accelerator neutrinos

Inclusive fluxes of secondary air-shower particles

Reconstruction of air shower muon lateral distribution functions using integrator and binary modes of underground muon detectors

Constraints on non-unitary neutrino mixing in light of atmospheric and reactor neutrino data

Neutrino Oscillations with Two $Δm^2$ Scales

Mathematical methods for neutrino cross-section extraction

DarkNews: a Python-based event generator for heavy neutral lepton production in neutrino-nucleus scattering

AM$^3$: An Open-Source Tool for Time-Dependent Lepto-Hadronic Modeling of Astrophysical Sources

Dark fluxes from electromagnetic cascades

From ray to spray: augmenting amplitudes and taming fast oscillations in fully numerical neutrino codes