// Experimental constraints, base
//double VEV = 246.0;
//double VEV2 = pow(VEV,2);
//double mh = 125.0;
//double mh2 = pow(mh,2);
//double VEV_red = VEV/sqrt(2.0);
// --------------------------------------------------
// Masses
//  double const EXP_val_Masses_Down[3] = {0.0029, 0.055, 2.9};
//  double const EXP_unc_Masses_Down[3] = {0.0012, 0.015, 0.1};
//  double const EXP_val_Masses_Up[3] = {0.00126, 0.68, 171.5};
//  double const EXP_unc_Masses_Up[3] = {0.0005, 0.1*0.68, 3.0};
double const Masses_Down[3] = {m_qd, m_qs, m_qb};
double const Masses_Up[3] = {m_qu, m_qc, m_qt};
// ---------------
  // Mixings
  double const EXP_val_CKM[2][3] = {{0.97425, 0.2252, 0.00375}, {0.230, 1.023, 0.041}};
  double const EXP_unc_CKM[2][3] = {{0.00022, 0.0009, 0.00046}, {0.011, 0.036, 0.0013}};
  double const EXP_val_gamma = 73.2*M_PI/180.0;
  double const EXP_unc_gamma = 7.0*M_PI/180.0;
  double const EXP_val_Wdecays = 2.002;
  double const EXP_unc_Wdecays = 0.027;
  // --------------------------------------------------