PMNS_Avg.h

Go to the documentation of this file.

 
 
#ifndef PMNS_Avg_H
#define PMNS_Avg_H
 
#include "PMNS_Fast.h"
#include "PremModel.h"
 
namespace OscProb {
 
  class PMNS_Avg : public PMNS_Fast {
    public:
      PMNS_Avg();          
      virtual ~PMNS_Avg(); 
 
      virtual void SetPremModel(OscProb::PremModel& prem);
 
      // Get probability averaged over a bin
      virtual double AvgProb(
          int flvi, int flvf, double E,
          double dE); 
 
      virtual double AvgProbLoE(
          int flvi, int flvf, double LoE,
          double dLoE); 
 
      virtual double AvgProb(
          int flvi, int flvf, double E, double cosT,
          double dcosT); 
 
      virtual double AvgProb(int flvi, int flvf, double E, double dE,
                             double cosT,
                             double dcosT); 
 
      virtual double AvgProbLoE(
          int flvi, int flvf, double LoE, double dLoE, double cosT,
          double dcosT); 
 
      // Get probability
      virtual double ExtrapolationProb(
          int flvi, int flvf, double E,
          double dE); 
 
      virtual double ExtrapolationProbLoE(
          int flvi, int flvf, double LoE,
          double dLoE); 
 
      virtual double ExtrapolationProbCosT(
          int flvi, int flvf, double cosT,
          double dcosT); 
 
    protected:
      virtual void InitializeTaylorsVectors(); 
 
      virtual void SetwidthBin(double dE,
                               double dcosT); 
 
      virtual void SetCosT(double cosT); 
 
      virtual vectorD GetSamplePoints(
          double LoE,
          double dLoE); 
 
      virtual vectorD GetSamplePoints(
          double E, double cosT,
          double dcosT); 
 
      virtual matrixC GetSamplePoints(double InvE, double dInvE, double cosT,
                                      double dcosT); 
 
      // Construction of the K matrices
      virtual void BuildKE(
          double L); 
      virtual void BuildKcosT(
          double L); 
 
      virtual double LnDerivative(); 
 
      // Rotation from mass to flavor basis
      virtual void rotateS(); 
      virtual void rotateK(); 
 
      // Multiplication rule between two pairs of (S,K)
      virtual void MultiplicationRuleS(); 
      virtual void MultiplicationRuleK(
          complexD K[3][3]); 
 
      virtual void SolveK(complexD K[3][3], vectorD& lambda,
                          matrixC& V); 
 
      // Propagating
      virtual void PropagatePathTaylor(
          NuPath p); 
      virtual void PropagateTaylor(); 
 
      virtual double AvgFormula(
          int flvi, int flvf, double dbin, vectorD flambda,
          matrixC fV); 
 
      virtual double AvgAlgo(
          int flvi, int flvf, double LoE, double dLoE,
          double L); 
 
      virtual double AvgAlgoCosT(
          int flvi, int flvf, double E, double cosT,
          double dcosT); 
 
      virtual double AvgAlgo(
          int flvi, int flvf, double InvE, double dInvE, double cosT,
          double dcosT); 
 
      virtual double AlgorithmDensityMatrix(
          int flvi, int flvf); 
 
      virtual void RotateDensityM(
          bool to_basis, matrixC V); 
 
      virtual void HadamardProduct(vectorD lambda,
                                   double  dbin); 
 
      virtual double AvgFormulaExtrapolation(
          int flvi, int flvf, double dbin, vectorD flambda,
          matrixC fV); 
 
      // Attributes
 
      matrixC fevolutionMatrixS; 
 
      matrixC fSflavor; 
      matrixC fKmass;   
      matrixC fKflavor; 
 
      double fdInvE; 
 
      complexD fKInvE[3][3]; 
      vectorD flambdaInvE;   
      matrixC fVInvE;        
 
      double fcosT;  
      double fdcosT; 
 
      complexD fKcosT[3]
                     [3];  
      vectorD flambdaCosT; 
      matrixC fVcosT;      
 
      matrixC fdensityMatrix; 
 
      // Variables for the compute of the derivation of one layer's length
      int    fLayer;
      int    fdl;
      double fDetRadius;
      double fminRsq;
 
      // Copy of the earth model used
      OscProb::PremModel fPrem;
  };
 
} // namespace OscProb
 
#endif