The rest of this is old information
The code
The code is available in this directory. A "make"
should build the library and various executables.
Overview
(Note: if your Netscape/Mozilla/Galeon browser under X has trouble
reading some of the symbols in this page, seek help
here)
This page describes code useful for studying long baseline neutrino
oscillations. It allows reading in of beam spectra, applying a full
matter effect oscillation, interacting in a detector and applying of
detector effects. The various data (spectra, oscillation probs, cross
sections) can be stored in files or MySQL database for fast lookup
Currently the only beam spectra supported is that from
GNUMI
with an
emphasis on a water cerenkov detector.
Oscillations
Oscillation is done using a full 3 ν, 2 Δm2 scale
calculation including details of the Earth's density profile from:
PREM
(Phys. of the Earth and Planetary Interiors, 25 (1981) 297-356, also
see
hep-ph/0002149
for excerpted fits from that paper). This
numerical calculation is a relatively straight-forward solution of the
evolution equation:
idν(x)/dx = [(1/2E)UM2U† + V]&nu(x);
Here, ν(x) is the neutrino (weak) state, U is the mixing matrix,
M2 is the mass matrix
(diag(m12,m22,m32), and V is the
matter interaction potential
(diag(GFNe(x),0,0)).
See
hep-ph/0002149
and
hep-ph/0106139
for details and an analytical approximation. The code has been
successfully compared to Fig. 3, 4 and 5 (note: it has been confirmed
with Irina that the latter two Figs were actually generated with
Δm2sol=1.0e-4 eV2 not 0.5e-4 as
the pre-print says. Much thanks to Irina for her ongoing help and
collaboration.)
Some of the programs
Most programs will give a usage message if you run them w/out any args.
gnumi-flux
reads in GNUMI ntuples and writes out flux tables
at given baseline and off-axis angle. See load-spect.sh for example use.
spect-load
used to load MySQL, also used in load-spect.sh.
xsec-load
loads cross sections from/to file to/from MySQL.
Cross sections files are in xsec/ and are currently from Milind's
nusim program used in
E889 proposal.
prob-load
loads/generates oscillation probabilites from/to
file to/from MySQL using the method described above.
elbo
main analysis (hist generator) program.
Important Library Elements
- Oscillation and matter effects: Prob.cc nuosc.cc earth.cc
NuTransODE.cc NuoscParam.cc Wiggle.cc
- Cross sections: Xsec.cc cross-section-func.cc
- Spectrum: Spect.cc reweight.cc
- Math: matrix.cc OdeFunc.cc ode-steppers.cc InterCompoundFunc.cc
InterSplineFunc.cc
- Utils: options.cc DBInt.cc
- HBOOK Ntuple interface: Ntuple.cc
Last modified: Fri Dec 6 14:32:40 EST 2002
") ?>