EigenPoint.cxx

Go to the documentation of this file.

 
//
// Struct to organise eigensystems for caching
//
//.................................................
//
// jcoelho\@apc.in2p3.fr
 
#include "EigenPoint.h"
 
using namespace OscProb;
 
//.............................................................................
EigenPoint::EigenPoint(int numNus, double e, NuPath p, bool n)
    : fEval(numNus, 0), fEvec(numNus, vectorC(numNus, 0))
{
  SetVars(e, p, n);
}
 
//.............................................................................
void EigenPoint::SetVars(double e, NuPath p, bool n)
{
  fEnergy = e;
  fPath   = p;
  fNubar  = n;
  SetNE();
}
 
//.............................................................................
void EigenPoint::SetNE()
{
  fNE = fEnergy * fPath.density * fPath.zoa;
  if (fNE < 1e-12) fNE = 1e-12;
  if (fNubar) fNE = -fNE;
}
 
//.............................................................................
bool EigenPoint::operator<(const EigenPoint& rhs) const
{
  if (fNE == rhs.fNE) return fPath.zoa < rhs.fPath.zoa;
  return fNE < rhs.fNE;
}
 
//.............................................................................
bool EigenPoint::operator==(const EigenPoint& rhs) const
{
  return (fNE == rhs.fNE) && (fPath.zoa == rhs.fPath.zoa);
}