//transcribed by Rudolph Magyar from tb_7_pi.top on 7/15/99 //

float dat1(int i, int k) {

  float dd[40][3] = {
       0.0244125,      0.0000000,      0.0000000,
       0.0732375,      1.0000000,      1.0000000,
       0.1220625,      0.0000000,      0.0000000,
       0.1708875,      1.0000000,      1.0000000,
       0.2197125,      1.0000000,      1.0000000,
       0.2685375,      9.0000000,      3.0000000,
       0.3173625,     10.0000000,      3.1622777,
       0.3661875,      4.0000000,      2.0000000,
       0.4150125,      7.0000000,      2.6457512,
       0.4638375,      8.0000000,      2.8284271,
       0.5126625,      8.0000000,      2.8284271,
       0.5614875,     17.0000000,      4.1231055,
       0.6103125,     28.0000000,      5.2915025,
       0.6591375,     34.0000000,      5.8309517,
       0.7079625,     52.0000000,      7.2111025,
       0.7567875,     51.0000000,      7.1414285,
       0.8056125,     55.0000000,      7.4161983,
       0.8544375,     59.0000000,      7.6811457,
       0.9032626,     64.0000000,      8.0000000,
       0.9520875,     66.0000000,      8.1240387,
       1.0009125,     57.0000000,      7.5498343,
       1.0497375,     42.0000000,      6.4807405,
       1.0985626,     33.0000000,      5.7445626,
       1.1473875,     39.0000000,      6.2449980,
       1.1962125,     17.0000000,      4.1231055,
       1.2450374,     13.0000000,      3.6055512,
       1.2938625,     11.0000000,      3.3166249,
       1.3426876,      6.0000000,      2.4494898,
       1.3915125,      3.0000000,      1.7320508,
       1.4403375,      4.0000000,      2.0000000,
       1.4891624,      3.0000000,      1.7320508,
       1.5379875,      2.0000000,      1.4142135,
       1.5868126,      2.0000000,      1.4142135,
       1.6356375,      0.0000000,      0.0000000,
       1.6844625,      1.0000000,      1.0000000,
       1.7332875,      2.0000000,      1.4142135,
       1.7821125,      1.0000000,      1.0000000,
       1.8309376,      0.0000000,      0.0000000,
       1.8797625,      0.0000000,      0.0000000,
       1.9285876,      0.0000000,      0.0000000,
  };

  return dd[i][k];
}

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

				// create a 2D histogram to define the range

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

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

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

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

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

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

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