PMNS_SiderealLIV.h

Go to the documentation of this file.

 
 
#ifndef PMNS_SiderealLIV_H
#define PMNS_SiderealLIV_H
 
#include "PMNS_Fast.h"
 
namespace OscProb {
 
  class PMNS_SiderealLIV : public PMNS_Fast {
    public:
      PMNS_SiderealLIV();          
      virtual ~PMNS_SiderealLIV(); 
 
      virtual void SetA(int flvi, int flvj, int coord, double val);
      virtual void SetC(int flvi, int flvj, int coord1, int coord2, double val);
 
      virtual void SetNeutrinoDirection(double zenith, double azimuth);
 
      virtual void SetTimeHours(double hours);
 
      virtual void SetColatitude(double chi);
 
      virtual void SetColatitude(double deg, double min, double sec = 0.0);
 
      virtual double GetColatitude() const;
 
      virtual double GetA(int flvi, int flvj, int coord);
      virtual double GetC(int flvi, int flvj, int coord1, int coord2);
 
    protected:
      virtual void UpdateHam();
 
      virtual void SolveHam();
 
      virtual void FillCache() {} 
 
      static bool ValidateAndOrder(int& flvi, int& flvj, const char* name);
 
      static bool ValidateCoord(int coord, const char* coordName,
                                const char* paramName);
 
      // ------------------------------------------------------------------------
      // Sidereal-frequency constants
      // ------------------------------------------------------------------------
      static constexpr double kSiderealDayHours = 23.9344696;
 
      static constexpr double kOmegaSidereal = 2.0 * M_PI / kSiderealDayHours;
 
      // ------------------------------------------------------------------------
      // SME coefficients (real; only the upper triangle in flavour indices is
      // used)
      // ------------------------------------------------------------------------
      double fa[3][3][3];    
      double fc[3][3][3][3]; 
 
      // ------------------------------------------------------------------------
      // Geometry / time state
      // ------------------------------------------------------------------------
      double fN[3]; 
      double fTime; 
 
    private:
      double fChi; 
  };
 
} // namespace OscProb
 
#endif