// Why doesn't paw do this?
{
  TCanvas *c1=new TCanvas("c1","Signal to Noise vs Bandwidth",200,10,700,500);
  //  c1->SetFillColor(42);
  c1->SetFillColor(0);
  c1->SetGrid();
  //  c1->GetFrame()->SetFillColor(21);
  c1->GetFrame()->SetFillColor(0);
  c1->GetFrame()->SetBorderSize(12);

  int n=6; // 6 bandwidth runs
  double bw_phi[n]={5.0, 5.0, 6.0, 7.0, 4.0, 5.0};
  double sn_phi[5][n]={11.4132, 10.4173, 11.1688, 12.9199, 10.9450, 11.0715,
		       12.1629, 12.1210, 13.0017, 14.8812, 11.9176, 12.0581,
		       12.0301, 10.5953, 13.0341, 14.0282, 11.2236, 11.9103,
		       13.9483, 13.1944, 14.1104, 15.9426, 12.1953, 13.0120,
		       12.1867, 13.1082, 12.1601, 13.2125, 11.9212, 12.2055};

  double bw_z[n]={4.0, 4.0, 5.0, 6.0, 3.0, 4.0};
  double sn_z[3][n]={8.5442,  8.3093,  8.8967,  9.0282,  8.3647,  8.9692,
		     11.2414, 11.0786, 11.9650, 11.9573, 11.2001, 11.0076,
		     11.7410, 11.1619, 11.1733, 11.7305, 11.0018, 11.1224};

  double bw_sas[n]={5.0, 5.0, 6.0, 7.0, 4.0, 5.0};
  double sn_sas[2][n]={10.1552, 10.1378, 10.3209, 11.8578,  9.3063, 10.2356,
		       11.5,    12.0,    12.7,    13.9,    11.2,    12.1};
  double sn_phi_avg[n];
  for (int i=0; i<6; i++){
    sn_phi_avg[i]=0.0;
    for (int j=0; j<5; j++){
      sn_phi_avg[i]+=sn_phi[j][i];
    }
    sn_phi_avg[i]/=5.0;
  }

  double sn_z_avg[n];
  for (int i=0; i<6; i++){
    sn_z_avg[i]=0.0;
    for (int j=0; j<3; j++){
      sn_z_avg[i]+=sn_z[j][i];
    }
    sn_z_avg[i]/=3.0;
  }

  double sn_sas_avg[n];
  for (int i=0; i<6; i++){
    sn_sas_avg[i]=0.0;
    for (int j=0; j<2; j++){
      sn_sas_avg[i]+=sn_sas[j][i];
    }
    sn_sas_avg[i]/=2.0;
  }

  double sigma_sn_phi[n]={0.22, 0.22, 0.22, 0.22, 0.22, 0.22};
  double sigma_sn_z[n]={0.17, 0.17, 0.17, 0.17, 0.17, 0.17};
  double sigma_sn_sas[n]={0.21, 0.21, 0.21, 0.21, 0.21, 0.21};
  double sigma_bw[n]={0.0, 0.0, 0.0, 0.0, 0.0, 0.0};

  TH2F *h=new TH2F("h","Signal to Noise vs Bandwidth for SVX",1,2.0,8.0,1,9.0,15.0);
  h->SetStats(0);
  h->GetXaxis()->SetTitle("Bandwidth Setting of SVX3 Chip");
  h->GetYaxis()->SetTitle("Signal/Noise");
  h->Draw();

  TGraphErrors *g_phi=new TGraphErrors(n, bw_phi, sn_phi_avg, sigma_bw, sigma_sn_phi);
  g_phi->SetMarkerSize(2.0);
  g_phi->SetMarkerColor(4);
  g_phi->SetMarkerStyle(21);
  g_phi->Draw("P");

  TGraphErrors *g_z=new TGraphErrors(n, bw_z, sn_z_avg, sigma_bw, sigma_sn_z);
  g_z->SetMarkerSize(2.0);
  g_z->SetMarkerColor(2);
  g_z->SetMarkerStyle(20);
  g_z->Draw("P");

  TGraphErrors *g_sas=new TGraphErrors(n, bw_sas, sn_sas_avg, sigma_bw, sigma_sn_sas);
  g_sas->SetMarkerSize(2.0);
  g_sas->SetMarkerColor(3);
  g_sas->SetMarkerStyle(22);
  g_sas->Draw("P");

  
  TLegend *leg=new TLegend(0.15,0.6,0.5,0.8);
  leg->AddEntry(g_phi,"Axial","p");
  leg->AddEntry(g_z,"90-degree Stereo","p");
  leg->AddEntry(g_sas,"Small Angle Stereo","p");
  //  leg->SetFillColor(21);
  leg->SetFillColor(0);
  leg->Draw();
}