#include <ostream>
void calc_all_yields(){
  //  .L /work/halld2/home/staylor/src/macros/calc_3pi_yields.C
  //.L /work/halld2/home/staylor/src/macros/calc_3g_yields.C

  double omega_3pi[45],domega_3pi[45];
  double omega_3g[45],domega_3g[45];
  double eta_3piq[45],deta_3piq[45];
  double phi_3pi[45],dphi_3pi[45]; 
  double phi_3g[45],dphi_3g[45];

  for (unsigned int i=0;i<45;i++){
    double yields[3];
    double uncs[3];
    calc_3pi_yields(i+1,yields,uncs); 
    eta_3piq[i]=yields[0];
    deta_3piq[i]=uncs[0];
    omega_3pi[i]=yields[1];
    domega_3pi[i]=uncs[1]; 
    phi_3pi[i]=yields[2];
    dphi_3pi[i]=uncs[2];
 
    calc_3g_yields(i+1,yields,uncs);
    omega_3g[i]=yields[1];
    domega_3g[i]=uncs[1]; 
    phi_3g[i]=yields[2];
    dphi_3g[i]=uncs[2];
  }
  cout << "Done!" <<endl;
  
  ofstream fout("omega_yields.dat");
  for (int j=0;j<45;j++){
    fout << 2.9+0.2*float(j) << " " << omega_3g[j] << " " << domega_3g[j] 
	 << " " << omega_3pi[j] << " " << domega_3pi[j] << endl;
  }
  fout.close();
  ofstream fout2("phi_yields.dat");
  for (int j=0;j<45;j++){
    fout2 << 2.9+0.2*float(j) << " " << phi_3g[j] << " " << dphi_3g[j] 
	  << " " << phi_3pi[j] << " " << dphi_3pi[j] << endl;
  }
  fout2.close(); 
  ofstream fout3("eta_3piq_yields.dat");
  for (int j=0;j<45;j++){
    fout3 << 2.9+0.2*float(j) << " " << eta_3piq[j] << " " << deta_3piq[j] << endl;
  }
  fout3.close(); 

  double egamma[45];
  for (int j=0;j<45;j++) egamma[j]=2.9+0.2*j;
  TGraphErrors *g1=new TGraphErrors(45,egamma,omega_3pi,0,domega_3pi);
  g1->Draw();
  TGraphErrors *g2=new TGraphErrors(45,egamma,phi_3pi,0,dphi_3pi);
  g2->Draw("same");  
  TGraphErrors *g3=new TGraphErrors(45,egamma,eta_3piq,0,deta_3piq);
  g3->Draw("same");
  
}