//transcribed by Rudolph Magyar from tb_10_pi.C on 7/15/99 //
float dat1(int i, int k) {

  float dd[40][4] = {
       0.0247025,      2.0000000,0.,      1.4142135,
       0.0741075,      1.0000000,0.,      1.0000000,
       0.1235125,      0.0000000,0.,      0.0000000,
       0.1729175,      2.0000000,0.,      1.4142135,
       0.2223225,      1.0000000,0.,      1.0000000,
       0.2717275,      1.0000000,0.,      1.0000000,
       0.3211325,      0.0000000,0.,      0.0000000,
       0.3705375,      1.0000000,0.,      1.0000000,
       0.4199425,      3.0000000,0.,      1.7320508,
       0.4693475,      2.0000000,0.,      1.4142135,
       0.5187525,      3.0000000,0.,      1.7320508,
       0.5681575,     12.0000000,0.,      3.4641016,
       0.6175625,     12.0000000,0.,      3.4641016,
       0.6669675,     18.0000000,0.,      4.2426405,
       0.7163725,     39.0000000,0.,      6.2449980,
       0.7657775,     35.0000000,0.,      5.9160800,
       0.8151825,     48.0000000,0.,      6.9282031,
       0.8645875,     48.0000000,0.,      6.9282031,
       0.9139925,     48.0000000,0.,      6.9282031,
       0.9633975,     37.0000000,0.,      6.0827627,
       1.0128025,     31.0000000,0.,      5.5677643,
       1.0622076,     28.0000000,0.,      5.2915025,
       1.1116126,     26.0000000,0.,      5.0990195,
       1.1610175,     18.0000000,0.,      4.2426405,
       1.2104225,      9.0000000,0.,      3.0000000,
       1.2598275,      9.0000000,0.,      3.0000000,
       1.3092325,      4.0000000,0.,      2.0000000,
       1.3586376,      8.0000000,0.,      2.8284271,
       1.4080426,      4.0000000,0.,      2.0000000,
       1.4574475,      0.0000000,0.,      0.0000000,
       1.5068525,      2.0000000,0.,      1.4142135,
       1.5562575,      2.0000000,0.,      1.4142135,
       1.6056626,      0.0000000,0.,      0.0000000,
       1.6550676,      0.0000000,0.,      0.0000000,
       1.7044725,      0.0000000,0.,      0.0000000,
       1.7538775,      0.0000000,0.,      0.0000000,
       1.8032825,      0.0000000,0.,      0.0000000,
       1.8526876,      0.0000000,0.,      0.0000000,
       1.9020926,      0.0000000,0.,      0.0000000,
       1.9514976,      0.0000000,0.,      0.0000000,
  };

  return dd[i][k];
}

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

				// create a 2D histogram to define the range

   TH2F *hr = new TH2F("hr","E/P for all 10 GeV pions, Test Beam",
		       1,0,2,1,0,60);

   //   hr->SetXTitle("");
   hr->Draw();

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

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

					// create graph
   int n1 = 40;
   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(kGreen);
   gr1->SetMarkerStyle(21);
   gr1->Draw("LP");

   //   TText* text = new TText(4.,0.25,"");
   //   text->Draw();     

}