#include <ostream>

void calc_t_yields(){
  //.L /work/halld2/home/staylor/src/macros/calc_3pi_vs_t_yields.C
  double omega_3pi[8][50],domega_3pi[8][50];
  double eta_3piq[8][50],deta_3piq[8][50];
  double phi_3pi[8][50],dphi_3pi[8][50]; 
  double t[50],dt[50];
  
  for (unsigned int ebin=0;ebin<7;ebin++){
    char hname[80];
    sprintf(hname,"t-dependence for E#gamma=%f GeV",7.3+float(ebin)*0.6);
    printf("%s\n",hname);
    double max=0.;
    for (unsigned int i=1;i<=50;i++){
      double yields[3];
      double uncs[3];
      calc_3pi_vs_t_yields(ebin,i,yields,uncs); 
      unsigned int k=i-1;
      eta_3piq[ebin][k]=yields[0];
      deta_3piq[ebin][k]=uncs[0];
      omega_3pi[ebin][k]=yields[1];
      if (yields[1]>max) max=yields[1];
      domega_3pi[ebin][k]=uncs[1]; 
      phi_3pi[ebin][k]=yields[2];
      dphi_3pi[ebin][k]=uncs[2];
      t[k]=0.04*i+0.02;
      dt[k]=0.02;
    }
    cout << "Done!" <<endl;
  }
  ofstream fout("omega_t_yields.dat");
  for (int j=0;j<50;j++){
  for (int k=0;k<8;k++){
    fout << omega_3pi[k][j] << " " << domega_3pi[k][j] << " ";
  }
  fout << endl;
  }
  fout.close();
  ofstream fout2("phi_t_yields.dat");
  for (int j=0;j<50;j++){ 
  for (int k=0;k<8;k++){
    fout2 << phi_3pi[k][j] << " " << dphi_3pi[k][j] << " ";
  }
  fout2 << endl;
}
  fout2.close();

}