PMNS_Base.cxx
//......................................................................
///
/// Get the mixing angle theta_ij in radians.
///
/// Requires that i<j. Will notify you if input is wrong.
/// If i>j, will assume reverse order and swap i and j.
///
/// @param i,j - the indices of theta_ij
///
Get the mixing angle theta_ij in radians.

Requires that i<j. Will notify you if input is wrong.
If i>j, will assume reverse order and swap i and j.
Get the mixing angle theta_ij in radians . Requires that i<j . Will notify you if input is wrong . If i > j will assume reverse order and swap i and j .
double PMNS_Base::GetAngle(int i, int j)
{
  if(i>j){
    cout << "Warning: First argument should be smaller than second argument" << endl;
    cout << "         Setting reverse order (Theta" << j << i << "). " << endl;
    int temp = i;
    i = j;
    j = temp;
  }
  if(i<1 || i>fNumNus-1 || j<2 || j>fNumNus){
    cout << "ERROR: Theta" << i << j << " not valid for " << fNumNus;
    cout << " neutrinos. Returning zero." << endl;
    return 0;
  }
  return fTheta[i-1][j-1];
}
double PMNS_Base :: GetAngle ( int i , int j ) { if ( i > j ) { cout << " " << endl ; cout << " " << j << i << " " << endl ; int temp = i ; i = j ; j = temp ; } if ( i < 1 || i > fNumNus - 1 || j < 2 || j > fNumNus ) { cout << " " << i << j << " " << fNumNus ; cout << " " << endl ; return 0 ; } return fTheta [ i - 1 ] [ j - 1 ] ; }
......................................................................
Get the mixing angle theta_ij in radians.
...................................................................... Get the mixing angle theta_ij in radians .