PMNS_Base.cxx
//.....................................................................
///
/// Compute the average probability of flvi going to flvf
/// over a bin of L/E with width dLoE.
///
/// The probabilities are weighted by (L/E)^-2 so that event
/// density is flat in energy. This avoids giving too much
/// weight to low energies. Better approximations would be
/// achieved if we used an interpolated event density.
///
/// This function works best for single paths.
/// In multiple paths the accuracy may be somewhat worse.
/// If needed, average over smaller L/E ranges.
///
/// Flavours are:
/// <pre>
///   0 = nue, 1 = numu, 2 = nutau
///   3 = sterile_1, 4 = sterile_2, etc.
/// </pre>
/// @param flvi - The neutrino starting flavour.
/// @param flvf - The neutrino final flavour.
/// @param LoE  - The neutrino  L/E value in the bin center in km/GeV
/// @param dLoE   - The L/E bin width in km/GeV
///
/// @return Average neutrino oscillation probability
///
Compute the average probability of flvi going to flvf
over a bin of L/E with width dLoE.

The probabilities are weighted by (L/E)^-2 so that event
density is flat in energy. This avoids giving too much
weight to low energies. Better approximations would be
achieved if we used an interpolated event density.

This function works best for single paths.
In multiple paths the accuracy may be somewhat worse.
If needed, average over smaller L/E ranges.


Compute the average probability of flvi going to flvf over a bin of L / E with width dLoE . The probabilities are weighted by ( L / E ) ^ - 2 so that event density is flat in energy . This avoids giving too much weight to low energies . Better approximations would be achieved if we used an interpolated event density . This function works best for single paths . In multiple paths the accuracy may be somewhat worse . If needed average over smaller L / E ranges .
double PMNS_Base::AvgProbLoE(int flvi, int flvf, double LoE, double dLoE)
{
  ResetToFlavour(flvi);
  return AvgProbLoE(fNuState, flvf, LoE, dLoE);
}
double PMNS_Base :: AvgProbLoE ( int flvi , int flvf , double LoE , double dLoE ) { ResetToFlavour ( flvi ) ; return AvgProbLoE ( fNuState , flvf , LoE , dLoE ) ; }
.....................................................................
Compute the average probability of flvi going to flvf
over a bin of L/E with width dLoE.
..................................................................... Compute the average probability of flvi going to flvf over a bin of L / E with width dLoE .