HNLs at DUNE, DarkQuest, and the SBN Program

Code written for the production of $Z^\prime$ and Heavy Neutral Leptons (HNL) for https://arxiv.org/abs/2410.08981, "Enhancing the Sensitivity to Seesaw Predictions in Gauged $B−L$ Scenarios", Francesco Capozzi, Bhaskar Dutta, Gajendra Gurung, Wooyoung Jang, Ian M. Shoemaker, Adrian Thompson, and Jaehoon Yu

Contact: Adrian Thompson [email protected]

Requirements:

• alplib: https://github.com/athompson-git/alplib
• git LFS (for large flux files)

Files include

• HNL branching ratios: hnl_decays.py
• Zprime branching ratios: zprime_decay_widths.py
• Flux simulation: hnl_flux.py

Fluxes generated with GEANT4, tabulated as 4-vectors (in GeV)

• Pi0 meson 4-momenta generated in the DUNE target geometry with GEANT4: data/dune_target_pi0_4vectors_1e5POT.txt
• Eta0 meson 4-momenta generated in the DUNE target geometry with GEANT4: data/dune_target_eta0_4vectors_1e5POT.txt

Constants

• darkquest_constants.py
• dune_constants.py

Data

• R ratio: data/r_ratio_by_sqrts-GeV.txt

We also acknowledge and use constraints on the HNL parameter space from https://www.hep.ucl.ac.uk/~pbolton/index.html