PMNS_Base.cxx
//......................................................................
///
/// Set the CP phase delta_ij in radians.
///
/// Requires that i+1<j. Will notify you if input is wrong.
/// If i>j, will assume reverse order and swap i and j.
///
/// This will check if value is changing to keep track of whether
/// the hamiltonian needs to be rebuilt.
///
/// @param i,j    - the indices of delta_ij
/// @param delta  - the value of delta_ij
///
Set the CP phase delta_ij in radians.

Requires that i+1<j. Will notify you if input is wrong.
If i>j, will assume reverse order and swap i and j.

This will check if value is changing to keep track of whether
the hamiltonian needs to be rebuilt.
Set the CP phase delta_ij in radians . Requires that i + 1<j . Will notify you if input is wrong . If i > j will assume reverse order and swap i and j . This will check if value is changing to keep track of whether the hamiltonian needs to be rebuilt .
void PMNS_Base::SetDelta(int i, int j, double delta)
{
  if(i>j){
    cout << "Warning: First argument should be smaller than second argument" << endl;
    cout << "         Setting reverse order (Delta" << j << i << "). " << endl;
    int temp = i;
    i = j;
    j = temp;
  }
  if(i<1 || i>fNumNus-1 || j<2 || j>fNumNus){
    cout << "ERROR: Delta" << i << j << " not valid for " << fNumNus;
    cout << " neutrinos. Doing nothing." << endl;
    return;
  }
  if(i+1==j){
    cout << "Warning: Rotation " << i << j << " is real. Doing nothing." << endl;
    return;
  }
  fBuiltHms *= (fDelta[i-1][j-1] == delta);
  fDelta[i-1][j-1] = delta;
}
void PMNS_Base :: SetDelta ( int i , int j , double delta ) { 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 ; } if ( i + 1 == j ) { cout << " " << i << j << " " << endl ; return ; } fBuiltHms *= ( fDelta [ i - 1 ] [ j - 1 ] == delta ) ; fDelta [ i - 1 ] [ j - 1 ] = delta ; }
......................................................................
Set the CP phase delta_ij in radians.
...................................................................... Set the CP phase delta_ij in radians .
// Check if value is actually changing