PMNS_NUNM.h

Go to the documentation of this file.

 
 
#ifndef PMNS_NUNM_H
#define PMNS_NUNM_H
 
#include "Definitions.h"
#include "PMNS_Fast.h"
#include <Eigen/Dense>
 
namespace OscProb {
 
  class PMNS_NUNM : public PMNS_Fast {
    public:
      PMNS_NUNM(int scale = 0); 
      virtual ~PMNS_NUNM();     
 
      virtual void SetAlpha(int i, int j, double val,
                            double phase); // i, j between 0 and 2
 
      virtual complexD GetAlpha(int i, int j);
 
      void SetNUNM(double alpha_11, double alpha_21, double alpha_31,
                   double alpha_22, double alpha_32, double alpha_33);
 
      // Set the diagonal real NUNM parameters
      virtual void SetAlpha_11(double a); 
      virtual void SetAlpha_22(double a); 
      virtual void SetAlpha_33(double a); 
 
      // Set the off-diagonal complex NUNM parameters
      virtual void    SetAlpha_21(double a,
                                  double phi); 
      virtual void    SetAlpha_31(double a,
                                  double phi); 
      virtual void    SetAlpha_32(double a,
                                  double phi); 
      virtual void    SetFracVnc(double f);
      virtual matrixD ProbMatrix(int nflvi, int nflvf);
 
      virtual void SetIsOscProbAvg(bool isOscProbAvg)
      {
        fIsOscProbAvg = true;
      } 
 
    protected:
      virtual void UpdateHam();
      virtual void Propagate();
      virtual void PropagatePath(NuPath p);
      vectorC      ApplyAlphaDagger(vectorC fState);
      vectorC      ApplyAlpha(vectorC fState);
      void         InitMatrix();
 
      int     fscale;
      double  fracVnc; // set fraction of matter potential affecting NC
      vectorC fNuStateBuffer;
 
      Eigen::Matrix<std::complex<double>, 3, 3> X;
      Eigen::Matrix<std::complex<double>, 3, 3> Alpha;
      Eigen::Matrix<std::complex<double>, 3, 3> V;
      Eigen::Matrix<std::complex<double>, 3, 3> Ham;
  };
 
} // namespace OscProb
 
#endif