PMNS_Base.cxx
//......................................................................
///
/// Get the neutrino mass eigenstate in vacuum
///
/// States are:
/// <pre>
///   0 = m_1, 1 = m_2, 2 = m_3, etc.
/// </pre>
/// @param mi - the mass eigenstate index
///
/// @return The mass eigenstate
///
Get the neutrino mass eigenstate in vacuum


Get the neutrino mass eigenstate in vacuum
std::vector<complexD> PMNS_Base::GetMassEigenstate(int mi){
  vector<complexD> oldState = fNuState;
  ResetToFlavour(mi);
  for(int j=0; j<fNumNus; j++){
  for(int i=0; i<j; i++){
    RotateState(i,j);
  }}
  vector<complexD> newState = fNuState;
  fNuState = oldState;
  return newState;
}
std :: vector < complexD > PMNS_Base :: GetMassEigenstate ( int mi ) { vector < complexD > oldState = fNuState ; ResetToFlavour ( mi ) ; for ( int j = 0 ; j < fNumNus ; j ++ ) { for ( int i = 0 ; i < j ; i ++ ) { RotateState ( i , j ) ; } } vector < complexD > newState = fNuState ; fNuState = oldState ; return newState ; }
......................................................................
Get the neutrino mass eigenstate in vacuum
...................................................................... Get the neutrino mass eigenstate in vacuum