osc::PMNS_NSI Class Reference

#include <PMNS_NSI.h>

Inheritance diagram for osc::PMNS_NSI:

Public Member Functions
	PMNS_NSI ()
virtual	~PMNS_NSI ()
void	SetNSI (double eps_ee, double eps_emu, double eps_etau, double eps_mumu, double eps_mutau, double eps_tautau, double delta_emu=0, double delta_etau=0, double delta_mutau=0)
Public Member Functions inherited from osc::PMNSOpt
	PMNSOpt ()
virtual	~PMNSOpt ()
virtual void	SetMix (double th12, double th23, double th13, double deltacp)
virtual void	SetDeltaMsqrs (double dm21, double dm32)
virtual void	PropMatter (double L, double E, double Ne, int anti=1)
virtual void	PropMatter (const std::list< double > &L, double E, const std::list< double > &Ne, int anti)
virtual void	PropVacuum (double L, double E, int anti=1)
virtual double	P (int flv) const
virtual void	ResetToFlavour (int flv=1)

Protected Member Functions
virtual void	SolveHam (double E, double Ne, int anti)
Protected Member Functions inherited from osc::PMNSOpt
virtual void	BuildHlv ()
virtual void	SetVacuumEigensystem (double E, int anti)

Protected Attributes
double	fEps_ee
double	fEps_mumu
double	fEps_tautau
complex	fEps_emu
complex	fEps_etau
complex	fEps_mutau
bool	fResetNSI
Protected Attributes inherited from osc::PMNSOpt
double	fDm21
double	fDm31
double	fTheta12
double	fTheta23
double	fTheta13
double	fDeltaCP
complex	fHlv [3][3]
complex	fEvec [3][3]
double	fEval [3]
complex	fNuState [3]
double	fCachedNe
double	fCachedE
int	fCachedAnti
bool	fBuiltHlv

Additional Inherited Members
Protected Types inherited from osc::PMNSOpt
typedef std::complex< double >	complex

Detailed Description

Definition at line 21 of file PMNS_NSI.h.

Constructor & Destructor Documentation

PMNS_NSI::PMNS_NSI

(

)

Definition at line 38 of file PMNS_NCI.cxx.

{
  this->SetMix(0.,0.,0.,0.);
  this->SetDeltaMsqrs(0.,0.);
  this->SetNSI(0.,0.,0.,0.,0.,0.,0.,0.,0.);
  this->ResetToFlavour(1);
  fCachedNe = 0.0;
  fCachedE =  1.0;
  fCachedAnti = 1;
}

PMNS_NSI::~PMNS_NSI ( )

virtual

Definition at line 49 of file PMNS_NCI.cxx.

                   {
}

Member Function Documentation

void PMNS_NSI::SetNSI	(	double	eps_ee,
		double	eps_emu,
		double	eps_etau,
		double	eps_mumu,
		double	eps_mutau,
		double	eps_tautau,
		double	delta_emu = 0,
		double	delta_etau = 0,
		double	delta_mutau = 0
	)

Definition at line 54 of file PMNS_NCI.cxx.

{

  fEps_ee     = eps_ee;
  fEps_mumu   = eps_mumu;
  fEps_tautau = eps_tautau;
  fEps_emu    = eps_emu   * complex(cos(delta_emu) ,   sin(delta_emu));
  fEps_etau   = eps_etau  * complex(cos(delta_etau) ,  sin(delta_etau));
  fEps_mutau  = eps_mutau * complex(cos(delta_mutau) , sin(delta_mutau));

  fResetNSI = true;

}

void PMNS_NSI::SolveHam ( double  E, double  Ne, int  anti  )

protectedvirtual

Reimplemented from osc::PMNSOpt.

Definition at line 71 of file PMNS_NCI.cxx.

{

  // Check if anything has changed before recalculating
  if(Ne!=fCachedNe || E!=fCachedE || anti!=fCachedAnti || !fBuiltHlv || fResetNSI){
    fCachedNe = Ne;
    fCachedE = E;
    fCachedAnti = anti;
    fResetNSI = false;
    this->BuildHlv();
  }
  else return;

  double lv = 2 * kGeV2eV*E / fDm31;  // Osc. length in eV^-1 
  double kr2GNe = kK2*M_SQRT2*kGf*Ne; // Matter potential in eV

  // Finish build Hamiltonian in matter with dimension of eV
  complex A[3][3];
  for(int i=0;i<3;i++){
    A[i][i] = fHlv[i][i]/lv;
    for(int j=i+1;j<3;j++){
      if(anti>0) A[i][j] = fHlv[i][j]/lv;
      else       A[i][j] = conj(fHlv[i][j])/lv;
    }
  }
  if(anti>0){
    A[0][0] += kr2GNe * (1 + fEps_ee);
    A[0][1] += kr2GNe * fEps_emu;
    A[0][2] += kr2GNe * fEps_etau;
    A[1][1] += kr2GNe * fEps_mumu;
    A[1][2] += kr2GNe * fEps_mutau;
    A[2][2] += kr2GNe * fEps_tautau;
  }
  else{
    A[0][0] -= kr2GNe * (1 + fEps_ee);
    A[0][1] -= kr2GNe * fEps_emu;
    A[0][2] -= kr2GNe * fEps_etau;
    A[1][1] -= kr2GNe * fEps_mumu;
    A[1][2] -= kr2GNe * fEps_mutau;
    A[2][2] -= kr2GNe * fEps_tautau;
  }

  // Solve Hamiltonian for eigensystem using the GLoBES method
  zheevh3(A,fEvec,fEval);

}

Member Data Documentation

double osc::PMNS_NSI::fEps_ee

protected

Definition at line 33 of file PMNS_NSI.h.

complex osc::PMNS_NSI::fEps_emu

protected

Definition at line 36 of file PMNS_NSI.h.

complex osc::PMNS_NSI::fEps_etau

protected

Definition at line 37 of file PMNS_NSI.h.

double osc::PMNS_NSI::fEps_mumu

protected

Definition at line 34 of file PMNS_NSI.h.

complex osc::PMNS_NSI::fEps_mutau

protected

Definition at line 38 of file PMNS_NSI.h.

double osc::PMNS_NSI::fEps_tautau

protected

Definition at line 35 of file PMNS_NSI.h.

bool osc::PMNS_NSI::fResetNSI

protected

Definition at line 39 of file PMNS_NSI.h.

---

The documentation for this class was generated from the following files:

• g4lbne/ProductionScripts/OscLib/PMNS_NSI.h
• g4lbne/ProductionScripts/OscLib/PMNS_NCI.cxx