osc::PMNS_NSI Class Reference

#include "/cvmfs/nova.opensciencegrid.org/externals/osclib/v00.07/src/OscLib/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)
virtual void	SetMix (const T &th12, const T &th23, const T &th13, const T &deltacp)
virtual void	SetDeltaMsqrs (const T &dm21, const T &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 T	P (int flv) const
virtual void	ResetToFlavour (int flv=1)

Protected Types
typedef std::complex< T >	complex

Protected Member Functions
virtual void	SolveHam (double E, double Ne, int anti)
virtual void	BuildHlv ()
virtual void	SetVacuumEigensystem (double E, int anti)

Protected Attributes
double	fEps_ee
	NSI parameter ee. More...
double	fEps_mumu
	NSI parameter mumu. More...
double	fEps_tautau
	NSI parameter tautau. More...
complex	fEps_emu
	NSI parameter emu. More...
complex	fEps_etau
	NSI parameter etau. More...
complex	fEps_mutau
	NSI parameter mutau. More...
bool	fResetNSI
	True when NSI parameters are changed. More...
T	fDm21
	m^2_2 - m^2_1 in vacuum More...
T	fDm31
	m^2_3 - m^2_1 in vacuum More...
T	fTheta12
	theta12 mixing angle More...
T	fTheta23
	theta23 mixing angle More...
T	fTheta13
	theta13 mixing angle More...
T	fDeltaCP
	CP violating phase. More...
complex	fHlv [3][3]
	dimensionless matrix H*lv More...
complex	fEvec [3][3]
	Eigenvectors of the Hamiltonian. More...
T	fEval [3]
	Eigenvalues of the Hamiltonian. More...
complex	fNuState [3]
	The neutrino current state. More...
double	fCachedNe
	Cached electron density. More...
double	fCachedE
	Cached neutrino energy. More...
int	fCachedAnti
	Cached nu/nubar selector. More...
bool	fBuiltHlv
	Tag to avoid rebuilding Hlv. More...

Detailed Description

Definition at line 21 of file PMNS_NSI.h.

Member Typedef Documentation

template<typename T>

typedef std::complex<T> osc::_PMNSOpt< T >::complex

protectedinherited

Definition at line 96 of file PMNSOpt.h.

Constructor & Destructor Documentation

PMNS_NSI::PMNS_NSI

(

)

Definition at line 38 of file PMNS_NSI.cxx.

References osc::_PMNSOpt< T >::fCachedAnti, osc::_PMNSOpt< T >::fCachedE, osc::_PMNSOpt< T >::fCachedNe, osc::_PMNSOpt< T >::ResetToFlavour(), osc::_PMNSOpt< T >::SetDeltaMsqrs(), osc::_PMNSOpt< T >::SetMix(), and SetNSI().

{
  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_NSI.cxx.

                   {
}

Member Function Documentation

template<typename T >

void _PMNSOpt::BuildHlv ( )

protectedvirtualinherited

Build H*lv, where H is the Hamiltonian in vacuum on flavour basis and lv is the oscillation length

Build H*lv, where H is the Hamiltonian in vacuum on flavour basis and lv is the oscillation length

This is a dimentionless hermitian matrix, so only the upper triangular part needs to be filled

The construction of the Hamiltonian avoids computing terms that are simply zero. This has a big impact in the computation time. This construction is described in DocDB-XXXX (to be posted)

Definition at line 140 of file PMNSOpt.cxx.

References std::cos(), h11, h12, std::sin(), and T.

Referenced by SolveHam().

{

  // Check if anything changed
  if(fBuiltHlv) return;

  // Create temp variables
  T sij, cij, h00, h11, h01;
  complex expCP(0,0), h02(0,0), h12(0,0);

  // Hamiltonian in mass base. Only one entry is variable.
  h11 = fDm21 / fDm31;

  // Rotate over theta12
  sij = sin(fTheta12);
  cij = cos(fTheta12);

  // There are 3 non-zero entries after rephasing so that h22 = 0
  h00 = h11 * sij * sij - 1;
  h01 = h11 * sij * cij;
  h11 = h11 * cij * cij - 1;

  // Rotate over theta13 with deltaCP
  sij = sin(fTheta13);
  cij = cos(fTheta13);
  expCP = complex(cos(fDeltaCP), -sin(fDeltaCP));

  // There are 5 non-zero entries after rephasing so that h22 = 0
  h02 = (-h00 * sij * cij) * expCP;
  h12 = (-h01 * sij) * expCP;
  h11 -= h00 * sij * sij;
  h00 *= cij * cij  -  sij * sij;
  h01 *= cij;

  // Finally, rotate over theta23
  sij = sin(fTheta23);
  cij = cos(fTheta23);

  // Fill the Hamiltonian rephased so that h22 = -h11
  // explicit construction of complex vars to avoid problem noted at top of file
  fHlv[0][0] = complex(h00 - 0.5 * h11, 0);
  fHlv[1][1] = complex(0.5 * h11 * (cij * cij - sij * sij)  +  2 * real(h12) * cij * sij, 0);
  fHlv[2][2] = -fHlv[1][1];

  // these are all constructed from complex variables (h02 and h12) so they are ok
  fHlv[0][1] = h02 * sij  +  h01 * cij;
  fHlv[0][2] = h02 * cij  -  h01 * sij;
  fHlv[1][2] = h12 - (h11 * cij + 2 * real(h12) * sij) * sij;

  // Tag as built
  fBuiltHlv = true;

}

template<typename T >

T _PMNSOpt::P ( int  flv ) const

virtualinherited

Return the probability of final state in flavour flv

Parameters

flv	- final flavor (0,1,2) = (nue,numu,nutau)

Compute oscillation probability of flavour flv

0 = nue, 1 = numu, 2 = nutau

Definition at line 384 of file PMNSOpt.cxx.

References ana::assert(), and norm.

Referenced by osc::OscCalcPMNS_NSI::P().

{
  assert(flv>=0 && flv<3);
  return norm(fNuState[flv]);
}

template<typename T >

void _PMNSOpt::PropMatter ( double  L, double  E, double  Ne, int  anti = 1  )

virtualinherited

Propagate a neutrino through a slab of matter

Parameters

L	- length of slab (km)
E	- neutrino energy in GeV
Ne	- electron number density of matter in mole/cm^3
anti	- +1 = neutrino case, -1 = anti-neutrino case

.....................................................................

Propagate the current neutrino state over a distance L in km with an energy E in GeV through constant matter of density Ne in mole/cm^3.

Parameters

anti	- +1 = neutrino case, -1 = anti-neutrino case

Definition at line 242 of file PMNSOpt.cxx.

References plot_validation_datamc::c, std::cos(), MECModelEnuComparisons::i, calib::j, osc::kKm2eV, and std::sin().

Referenced by osc::OscCalcPMNS_NSI::P().

{

  // Solve Hamiltonian
  this->SolveHam(E, Ne, anti);

  // Store coefficients of propagation eigenstates
  complex nuComp[3];

  for(int i=0;i<3;i++){
    nuComp[i] = complex(0,0);
    for(int j=0;j<3;j++){
      nuComp[i] += fNuState[j] * conj(fEvec[j][i]);
    }
  }

  for(int i=0;i<3;i++)fNuState[i] = complex(0,0);

  // Propagate neutrino state
  for(int j=0;j<3;j++){
    auto s = sin(-fEval[j] * kKm2eV * L);
    auto c = cos(-fEval[j] * kKm2eV * L);

    complex jPart = complex(c, s) * nuComp[j];
    for(int i=0;i<3;i++){
      fNuState[i] += jPart * fEvec[i][j];
    }
  }

}

template<typename T >

void _PMNSOpt::PropMatter ( const std::list< double > &  L, double  E, const std::list< double > &  Ne, int  anti  )

virtualinherited

Do several layers in a row. L and Ne must have the same length

Definition at line 278 of file PMNSOpt.cxx.

{
  if (L.size()!=Ne.size()) abort();
  auto Li  = L.begin();
  auto Lend = L.end();
  auto Ni = Ne.begin();
  for (; Li!=Lend; ++Li, ++Ni) {
    // For very low densities, use vacumm
    static const double kRhoCutoff = 1.0E-6; // Ne in moles/cm^3
    if (*Ni<kRhoCutoff) this->PropVacuum(*Li, E, anti);
    else                this->PropMatter(*Li, E, *Ni, anti);
  }
}

template<typename T >

void _PMNSOpt::PropVacuum ( double  L, double  E, int  anti = 1  )

virtualinherited

Propagate a neutrino through vacuum

Parameters

L	- length of slab (km)
E	- neutrino energy in GeV
anti	- +1 = neutrino case, -1 = anti-neutrino case

.....................................................................

Propagate the current neutrino state over a distance L in km with an energy E in GeV through vacuum

Parameters

anti	- +1 = neutrino case, -1 = anti-neutrino case

Definition at line 337 of file PMNSOpt.cxx.

References stan::math::exp(), MECModelEnuComparisons::i, calib::j, osc::kKm2eV, CLHEP::L, and T.

{

  this->SetVacuumEigensystem(E, anti);

  complex nuComp[3];

  for(int i=0;i<3;i++){
    nuComp[i] = 0;
    for(int j=0;j<3;j++){
      nuComp[i] += fNuState[j] * conj(fEvec[j][i]);
    }
  }

  const T km2EvL = kKm2eV*L;  // needed for the templated multiplication below to work
  for(int i=0;i<3;i++){
    fNuState[i] = 0;
    for(int j=0;j<3;j++){
      complex iEval(0.0,fEval[j]);
      fNuState[i] +=  exp(-iEval * km2EvL) * nuComp[j] * fEvec[i][j];
    }
  }

}

template<typename T >

void _PMNSOpt::ResetToFlavour ( int  flv = 1 )

virtualinherited

Erase memory of neutrino propagate and reset neutrino to pure flavour flv. Preserves values of mixing and mass-splittings

Parameters

flv	- final flavor (0,1,2) = (nue,numu,nutau)

Reset the neutrino state back to a pure flavour where it starts

Definition at line 368 of file PMNSOpt.cxx.

References MECModelEnuComparisons::i.

Referenced by osc::OscCalcPMNS_NSI::P(), and PMNS_NSI().

{
  int i;
  for (i=0; i<3; ++i){
    if (i==flv) fNuState[i] = complex(1, 0);
    else        fNuState[i] = complex(0, 0);
  }
}

template<typename T >

void _PMNSOpt::SetDeltaMsqrs ( const T &  dm21, const T &  dm32  )

virtualinherited

Set the mass-splittings

Parameters

dm21	- m2^2-m1^2 in eV^2
dm32	- m3^2-m2^2 in eV^2

Set the mass-splittings. These are m_2^2-m_1^2, m_3^2-m_2^2 and m_3^2-m_1^2 in eV^2

Definition at line 115 of file PMNSOpt.cxx.

Referenced by osc::OscCalcPMNS_NSI::P(), and PMNS_NSI().

{

  fDm21 = dm21;
  fDm31 = dm32 + dm21;

  if(fDm31==0) fDm31 = 1.0e-12;

  fBuiltHlv = false;

}

template<typename T >

void _PMNSOpt::SetMix ( const T &  th12, const T &  th23, const T &  th13, const T &  deltacp  )

virtualinherited

Set the parameters of the PMNS matrix

Parameters

th12	- The angle theta_12 in radians
th23	- The angle theta_23 in radians
th13	- The angle theta_13 in radians
deltacp	- The CPV phase delta_CP in radians

Definition at line 97 of file PMNSOpt.cxx.

References th12, th13, and th23.

Referenced by osc::OscCalcPMNS_NSI::P(), and PMNS_NSI().

{

  fTheta12 = th12;
  fTheta23 = th23;
  fTheta13 = th13;
  fDeltaCP = deltacp;

  fBuiltHlv = false;

}

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_NSI.cxx.

References std::cos(), fEps_ee, fEps_emu, fEps_etau, fEps_mumu, fEps_mutau, fEps_tautau, fResetNSI, and std::sin().

Referenced by osc::OscCalcPMNS_NSI::P(), and PMNS_NSI().

{

  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;

}

template<typename T >

void _PMNSOpt::SetVacuumEigensystem ( double  E, int  anti  )

protectedvirtualinherited

Set the eigensystem to the analytic solution of the vacuum Hamiltonian

Parameters

E	- neutrino energy in GeV
anti	- +1 = neutrino case, -1 = anti-neutrino case

.....................................................................

We know the vacuum eigensystem, so just write it explicitly The eigenvalues depend on energy, so E needs to be provided in GeV

Parameters

anti	- +1 = neutrino case, -1 = anti-neutrino case

Definition at line 303 of file PMNSOpt.cxx.

References std::cos(), E, osc::kGeV2eV, std::sin(), and T.

{

  T       s12, s23, s13, c12, c23, c13;
  complex expidelta(cos(fDeltaCP), anti * sin(fDeltaCP));

  s12 = sin(fTheta12);  s23 = sin(fTheta23);  s13 = sin(fTheta13);
  c12 = cos(fTheta12);  c23 = cos(fTheta23);  c13 = cos(fTheta13);

  fEvec[0][0] =  complex(c12*c13, 0);
  fEvec[0][1] =  complex(s12*c13, 0);
  fEvec[0][2] =  s13*conj(expidelta);

  fEvec[1][0] = -s12*c23-c12*s23*s13*expidelta;
  fEvec[1][1] =  c12*c23-s12*s23*s13*expidelta;
  fEvec[1][2] =  complex(s23*c13, 0);

  fEvec[2][0] =  s12*s23-c12*c23*s13*expidelta;
  fEvec[2][1] = -c12*s23-s12*c23*s13*expidelta;
  fEvec[2][2] =  complex(c23*c13, 0);

  fEval[0] = 0;
  fEval[1] = fDm21 / (2 * kGeV2eV*E);
  fEval[2] = fDm31 / (2 * kGeV2eV*E);

}

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

protectedvirtual

Solve the full Hamiltonian for eigenvectors and eigenvalues

Parameters

E	- neutrino energy in GeV
Ne	- electron number density of matter in mole/cm^3
anti	- +1 = neutrino case, -1 = anti-neutrino case

Solve the full Hamiltonian with Non-Standard Interactions for eigenvectors and eigenvalues.

Reimplemented from osc::_PMNSOpt< T >.

Definition at line 75 of file PMNS_NSI.cxx.

References genie::units::A, osc::_PMNSOpt< T >::BuildHlv(), E, osc::_PMNSOpt< T >::fBuiltHlv, osc::_PMNSOpt< T >::fCachedAnti, osc::_PMNSOpt< T >::fCachedE, osc::_PMNSOpt< T >::fCachedNe, osc::_PMNSOpt< T >::fDm31, fEps_ee, fEps_emu, fEps_etau, fEps_mumu, fEps_mutau, fEps_tautau, osc::_PMNSOpt< T >::fEval, osc::_PMNSOpt< T >::fEvec, osc::_PMNSOpt< T >::fHlv, fResetNSI, MECModelEnuComparisons::i, calib::j, osc::kGeV2eV, osc::kGf, osc::kK2, and zheevh3().

{

  // 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 * conj(fEps_emu);
    A[0][2] -= kr2GNe * conj(fEps_etau);
    A[1][1] -= kr2GNe * fEps_mumu;
    A[1][2] -= kr2GNe * conj(fEps_mutau);
    A[2][2] -= kr2GNe * fEps_tautau;
  }

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

}

Member Data Documentation

template<typename T>

bool osc::_PMNSOpt< T >::fBuiltHlv

protectedinherited

Tag to avoid rebuilding Hlv.

Definition at line 126 of file PMNSOpt.h.

Referenced by SolveHam().

template<typename T>

int osc::_PMNSOpt< T >::fCachedAnti

protectedinherited

Cached nu/nubar selector.

Definition at line 125 of file PMNSOpt.h.

Referenced by PMNS_NSI(), and SolveHam().

template<typename T>

double osc::_PMNSOpt< T >::fCachedE

protectedinherited

Cached neutrino energy.

Definition at line 124 of file PMNSOpt.h.

Referenced by PMNS_NSI(), and SolveHam().

template<typename T>

double osc::_PMNSOpt< T >::fCachedNe

protectedinherited

Cached electron density.

Definition at line 123 of file PMNSOpt.h.

Referenced by PMNS_NSI(), and SolveHam().

template<typename T>

T osc::_PMNSOpt< T >::fDeltaCP

protectedinherited

CP violating phase.

Definition at line 118 of file PMNSOpt.h.

template<typename T>

T osc::_PMNSOpt< T >::fDm21

protectedinherited

m^2_2 - m^2_1 in vacuum

Definition at line 113 of file PMNSOpt.h.

template<typename T>

T osc::_PMNSOpt< T >::fDm31

protectedinherited

m^2_3 - m^2_1 in vacuum

Definition at line 114 of file PMNSOpt.h.

Referenced by SolveHam().

double osc::PMNS_NSI::fEps_ee

protected

NSI parameter ee.

Definition at line 37 of file PMNS_NSI.h.

Referenced by SetNSI(), and SolveHam().

complex osc::PMNS_NSI::fEps_emu

protected

NSI parameter emu.

Definition at line 40 of file PMNS_NSI.h.

Referenced by SetNSI(), and SolveHam().

complex osc::PMNS_NSI::fEps_etau

protected

NSI parameter etau.

Definition at line 41 of file PMNS_NSI.h.

Referenced by SetNSI(), and SolveHam().

double osc::PMNS_NSI::fEps_mumu

protected

NSI parameter mumu.

Definition at line 38 of file PMNS_NSI.h.

Referenced by SetNSI(), and SolveHam().

complex osc::PMNS_NSI::fEps_mutau

protected

NSI parameter mutau.

Definition at line 42 of file PMNS_NSI.h.

Referenced by SetNSI(), and SolveHam().

double osc::PMNS_NSI::fEps_tautau

protected

NSI parameter tautau.

Definition at line 39 of file PMNS_NSI.h.

Referenced by SetNSI(), and SolveHam().

template<typename T>

T osc::_PMNSOpt< T >::fEval[3]

protectedinherited

Eigenvalues of the Hamiltonian.

Definition at line 121 of file PMNSOpt.h.

Referenced by SolveHam().

template<typename T>

complex osc::_PMNSOpt< T >::fEvec[3][3]

protectedinherited

Eigenvectors of the Hamiltonian.

Definition at line 120 of file PMNSOpt.h.

Referenced by SolveHam().

template<typename T>

complex osc::_PMNSOpt< T >::fHlv[3][3]

protectedinherited

dimensionless matrix H*lv

Definition at line 119 of file PMNSOpt.h.

Referenced by SolveHam().

template<typename T>

complex osc::_PMNSOpt< T >::fNuState[3]

protectedinherited

The neutrino current state.

Definition at line 122 of file PMNSOpt.h.

bool osc::PMNS_NSI::fResetNSI

protected

True when NSI parameters are changed.

Definition at line 43 of file PMNS_NSI.h.

Referenced by SetNSI(), and SolveHam().

template<typename T>

T osc::_PMNSOpt< T >::fTheta12

protectedinherited

theta12 mixing angle

Definition at line 115 of file PMNSOpt.h.

template<typename T>

T osc::_PMNSOpt< T >::fTheta13

protectedinherited

theta13 mixing angle

Definition at line 117 of file PMNSOpt.h.

template<typename T>

T osc::_PMNSOpt< T >::fTheta23

protectedinherited

theta23 mixing angle

Definition at line 116 of file PMNSOpt.h.

---

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

• OscLib/PMNS_NSI.h
• OscLib/PMNS_NSI.cxx