OscProb is a small set of classes aimed at computing exact neutrino oscillation probabilities with a few different models.
OscProb contains a basic framework for computing neutrino oscillation probabilities. It is integrated into ROOT, so that each class can be used as you would any ROOT class.
Available classes are:
- PremModel: Used for determining neutrino paths through the earth
- PMNS_Fast: Standard 3-flavour oscillations
- PMNS_Iter: Standard 3-flavour oscillations (iterative approximation)
- PMNS_Sterile: Oscillations with any number of neutrinos
- PMNS_NSI: Oscillations with 3 flavours including vector Non-Standard Interactions
- PMNS_SNSI: Oscillations with 3 flavours including scalar Non-Standard Interactions
- PMNS_Deco: Oscillations with 3 flavours including a simple decoherence model
- PMNS_LIV: Oscillations with 3 flavours including Lorentz Invariance Violations
- PMNS_Decay: Oscillations with 3 flavours including neutrino decays
- PMNS_NUNM: Oscillations with 3 flavours including non-unitary neutrino mixing
- Absorption: Computes absorption probabilities for high-energy neutrinos
A few example macros on how to use OscProb are available in a tutorial directory.
OscProb is very easy to install. The only requirements is to have ROOT installed.
NEW: Thanks to Jacek Holeczek, OscProb now also builds with ROOT 6!!
Details: Eigen
In order to compile the PMNS_Decay and PMNS_Sterile classes, it is necessary to
donwload the external Eigen library. This library is added as a submodule and will be downloaded by make.
Alternatively, you can trigger them manually with:
- During cloning:
git clone --recurse-submodules https://github.com/joaoabcoelho/OscProb.git - After clonning:
git submodule update --init
Once you have ROOT setup, simply do:
cd OscProb
makeA shared library will be produced: lib/libOscProb.so
This should take a few seconds and you are all set.
Just load the shared library in your ROOT macros with:
gSystem->Load("/full/path/to/libOscProb.so");Or use the LoadOscProb.C macro (see below).
In the directory OscProb/tutorial you will find a few macros with examples using OscProb.
Two macros are particularly useful:
simpleExamples.C: Contains some short pieces of code on how to perform different tasks.MakeOscillogram.C: Runs a full example of how to plot an oscillogram with the PREM model.
Additionally, these macros contain useful tools:
LoadOscProb.C: Searches for the OscProb library in your current directory, parent directory, or library path, and then loads it.SetNiceStyle.C: Provides simple tools to make your plots look nicer. Feel free to use it anytime you're making plots, even if you're not running OscProb. This is completely independent of OscProb.
To run macros you will need to preload the OscProb library, e.g:
root -l LoadOscProb.C MakeOscillogram.Croot -l LoadOscProb.C MakeOscillogram.C+You can use OscProb in Python using PyROOT by simply loading the shared library:
import ROOT
ROOT.gSystem.Load('/path/to/OscProb/libOscProb.so')
p = ROOT.OscProb.PMNS_Fast()
print(p.Prob(1, 1, 5.0))More detailed documentation of the code can be found in a Doxygen page here: