sntools

Event generator for simulating supernova neutrino bursts in water Cherenkov or liquid scintillator detectors.

Input

Text file(s) containing information about neutrino fluxes produced by the supernova. sntools distinguishes between three flavours: nu_e, anti-nu_e and nu_x (where nu_x stands for nu_mu or nu_tau or their respective antineutrinos). The following input formats are supported; see the source files in the formats/ directory for details.

Nakazato format

Used by recent simulations by the Nakazato group. Fluxes for 13 and 20 solar mass progenitors are included as fluxes/intp1301.data and fluxes/intp2001.data. A description of the format and fluxes for more progenitors are available online. If you use these included models in your work, please cite Nakazato et al., ApJ Supp. 205 (2013) 2.

Gamma format

Text file containing time, mean energy, mean squared energy and luminosity. These parameters describe a Gamma distribution, which is a good fit to the true spectrum. See fluxes/sample-gamma.txt for an unphysical sample file.

Princeton format

Used in recent simulations by the Princeton group.

Totani format

Used in historical simulation by Totani et al., which is the baseline model in the Hyper-Kamiokande Design Report.

Interaction Channels

sntools supports the main interaction channels in water and liquid scintillator.

For water Cherenkov detectors, like Hyper-Kamiokande, these are inverse beta decay, elastic scattering on electrons and charged-current interactions of nu_e and anti-nu_e on oxygen-16 nuclei. For liquid scintillator detectors, these are inverse beta decay, elastic scattering on electrons, charged-current interactions of nu_e and anti-nu_e on carbon-12 nuclei and neutral-current interactions on carbon-12 nuclei. For details, see the files in interaction_channels/.

Output

A .kin file in the NUANCE format used by the /mygen/vecfile options in WCSim. See the format documentation for details.

Typical Usage

python genevts.py fluxes/intp2001.data --format=nakazato -o outfile.kin --hierarchy=normal --detector=HyperK

See

python genevts.py -h

for a full description of these and other options.

sntools requires Python 2.7 or 3.x, numpy 1.8 (or higher) and scipy 0.13 (or higher).