PMNS_Base.cxx
//......................................................................
///
/// Set the mass-splitting dm_j1 = (m_j^2 - m_1^2) in eV^2
///
/// Requires that j>1. Will notify you if input is wrong.
///
/// This will check if value is changing to keep track of whether
/// the hamiltonian needs to be rebuilt.
///
/// @param j    - the index of dm_j1
/// @param dm   - the value of dm_j1
///


Requires that j>1. Will notify you if input is wrong.

This will check if value is changing to keep track of whether
the hamiltonian needs to be rebuilt.
Requires that j > 1 . Will notify you if input is wrong . This will check if value is changing to keep track of whether the hamiltonian needs to be rebuilt .
void PMNS_Base::SetDm(int j, double dm)
{
  if(j<2 || j>fNumNus){
    cout << "ERROR: Dm" << j << "1 not valid for " << fNumNus;
    cout << " neutrinos. Doing nothing." << endl;
    return;
  }
  fBuiltHms *= (fDm[j-1] == dm);
  fDm[j-1] = dm;
}
void PMNS_Base :: SetDm ( int j , double dm ) { if ( j < 2 || j > fNumNus ) { cout << " " << j << " " << fNumNus ; cout << " " << endl ; return ; } fBuiltHms *= ( fDm [ j - 1 ] == dm ) ; fDm [ j - 1 ] = dm ; }
......................................................................
Set the mass-splitting dm_j1 = (m_j^2 - m_1^2) in eV^2
...................................................................... Set the mass - splitting dm_j1 = ( m_j^2 - m_1^2 ) in eV^2
// Check if value is actually changing