This is the development page for the Python tool nuCraft.

NuCraft calculates oscillation probabilities for atmospheric neutrinos, taking into account matter effects and the Earth's atmosphere, and supports an arbitrary number of sterile neutrino flavors with easily configurable continuous Earth models.
Continous modelling of the Earth instead of the often-used approximation of four layers with constant density and consideration of the smearing of baseline lenghts due to the variable neutrino production heights in Earth's atmosphere each lead to deviations of 10% or more for conventional neutrinos between 1 and 10 GeV.

NuCraft heavily relies on NumPy and SciPy. It has been tested with several versions of NumPy, beginning at 1.5.1, and SciPy, beginning at 0.8.0. It is possible that it will work well with earlier versions of both packages. NuCraft is compatible with both Python 2 (2.6+) and Python 3 (tested with 3.3 and 3.7).

The source code is available through SVN, further information can be found in the accompanying README file, in class and in-line documentation, and in the IceCube internal report icecube/201303001, which was made available to the public here.
A more concise, dedicated publication has been published in Computer Physics Communication 197 (2015) 185-189, DOI: 10.1016/j.cpc.2015.07.010, with a preprint available at arXiv.

Note: A recent update of HEPForge complicated anonymous SVN checkouts. You can still use the repository viewer to download raw files manually, or get this tarball of the most recent revision: nuCraft-r22.tar.gz

If you have any questions, please contact the author Marius Wallraff, mwallraff()physik.rwth-aachen.de.