float dat1(int i, int k) {
//translated by Rudolph Magyar from mbias_zeros_vs_eta.top 7/15/99 // 
float dd[9][4] = {
       0.5000000,      0.1375000,0.,      0.0099967,
       1.5000000,      0.1166483,0.,      0.0090812,
       2.5000000,      0.1326646,0.,      0.0093344,
       3.5000000,      0.1258936,0.,      0.0090951,
       4.5000000,      0.1641803,0.,      0.0108177,
       5.5000000,      0.1759043,0.,      0.0107929,
       6.5000000,      0.2051428,0.,      0.0122784,
       7.5000000,      0.2474457,0.,      0.0133292,
       8.5000000,      0.2918572,0.,      0.0157867,
  };

  return dd[i][k];
}

float dat2(int i, int k) {

  float dd[9][4] = {
       0.5000000,      1.90,0., .14,
       1.5000000,      1.58,0., .12,
       2.5000000,      1.74,0., .12,
       3.5000000,      1.26,0., .09,
       4.5000000,      1.56,0., .10,
       5.5000000,      0.94,0., .08,
       6.5000000,      1.43,0., .08,
       7.5000000,      1.39,0., .07,
       8.5000000,      1.53,0., .08,
  };

  return dd[i][k];
}

int mbias_zeros_vs_eta() {
   gROOT->Reset();
   c1 = new TCanvas("mbias_zeros_vs_eta","mbias_zeros_vs_eta",
		    0,0,600,400);

				// create a 2D histogram to define the range

   TH2F *hr = new TH2F("hr","\"Zeros\" vs tower Eta, P 1-2 GeV",
		       1,0.,9,1,0,0.5);

   hr->SetXTitle("Tower Eta");
   hr->SetYTitle("f_0");
   hr->Draw();

   c1->GetFrame()->SetFillColor(21);
   c1->GetFrame()->SetBorderSize(12);

   float x[1000], y[1000], ex[1000], ey[1000];

					// create first graph
   int n1 = 9;
   for (int i=0; i<n1; i++) {
     x [i] = dat1(i,0);
     y [i] = dat1(i,1);
     ex[i] = dat1(i,2);
     ey[i] = dat1(i,3);
   }

   gr1 = new TGraphErrors(n1,x,y,ex,ey);
   gr1->SetMarkerColor(kBlue);
   gr1->SetMarkerStyle(21);
   gr1->Draw("LP");


   TText* text = new TText(0.5,0.22,"Blue: f_0");
   text->Draw();  

   //This should be on a seperate graph!!!

   int n2 = 9;
   for (int i=0; i<n2; i++) {
     x [i] = dat2(i,0);
     y [i] = dat2(i,1)/4; //notice the renormalization
     ex[i] = dat2(i,2);
     ey[i] = dat2(i,3)/4;
   }

   gr2 = new TGraphErrors(n2,x,y,ex,ey);
   gr2->SetMarkerColor(kRed);
   gr2->SetMarkerStyle(20);
   gr2->Draw("LP");

   TText* text = new TText(0.5,0.32,"Red: f_0/Dead");
   text->Draw();     
}