#include <ostream>
void calc_all_yields(){
  double omega_3pi[24],domega_3pi[24];
  double omega_3g[24],domega_3g[24];
  double eta_3piq[24],deta_3piq[24];
  double phi_3pi[24],dphi_3pi[24];

  for (unsigned int i=22;i<45;i++){
    unsigned int k=i-22;
    double yields[3];
    double uncs[3];
    calc_3pi_yields(i,yields,uncs); 
    eta_3piq[k]=yields[0];
    deta_3piq[k]=uncs[0];
    omega_3pi[k]=yields[1];
    domega_3pi[k]=uncs[1]; 
    phi_3pi[k]=yields[2];
    phi_3pi[k]=uncs[2];
 
    calc_3g_yields(i,yields,uncs);
    omega_3g[k]=yields[1];
    domega_3g[k]=uncs[1];
  }
  cout << "Done!" <<endl;
  
  ofstream fout("omega_yields.dat");
  for (int j=0;j<24;j++){
    fout << omega_3g[j] << " " << domega_3g[j] << " " << omega_3pi[j] 
	<< " " << domega_3pi[j] << endl;
  }
  fout.close();
  
}