//transcribed by Rudolph Magyar from tb_50_pi_in_crack.C on 7/16/99 //

float dat1(int i, int k) {
  float dd[100][3] = {
   1.1250,  4789.0,  69.2026,
   3.3750,  938.00,  30.6268,
   5.6250,  542.00,  23.2809,
   7.8750,  859.00,  29.3087,
   10.125,  990.00,  31.4643,
   12.375,  1006.0,  31.7175,
   14.625,  1173.0,  34.2491,
   16.875,  1358.0,  36.8511,
   19.125,  1740.0,  41.7133,
   21.375,  2238.0,  47.3075,
   23.625,  2868.0,  53.5537,
   25.875,  3550.0,  59.5819,
   28.125,  3970.0,  63.0079,
   30.375,  3982.0,  63.1031,
   32.625,  3667.0,  60.5558,
   34.875,  3024.0,  54.9909,
   37.125,  2256.0,  47.4974,
   39.375,  1358.0,  36.8511,
   41.625,  752.00,  27.4226,
   43.875,  355.00,  18.8414,
   46.125,  142.00,  11.9164,
   48.375,  71.000,  8.42615,
   50.625,  31.000,  5.56776,
   52.875,  13.000,  3.60555,
   55.125,  11.000,  3.31662,
   57.375,  5.0000,  2.23607,
   59.625,  8.0000,  2.82843,
   61.875,  4.0000,  2.00000,
   64.125,  4.0000,  2.00000,
   66.375,  3.0000,  1.73205,
   68.625,  5.0000,  2.23607,
   70.875,  4.0000,  2.00000,
   73.125,  1.0000,  1.00000,
   75.375,  4.0000,  2.00000,
   77.625,  3.0000,  1.73205,
   79.875,  0.0000,  0.00000,
   82.125,  0.0000,  0.00000,
   84.375,  0.0000,  0.00000,
   86.625,  0.0000,  0.00000,
   88.875,  0.0000,  0.00000,
   91.125,  2.0000,  1.41421,
   93.375,  0.0000,  0.00000,
   95.625,  0.0000,  0.00000,
   97.875,  0.0000,  0.00000,
   100.12,  0.0000,  0.00000,
   102.37,  0.0000,  0.00000,
   104.62,  0.0000,  0.00000,
   106.87,  0.0000,  0.00000,
   109.12,  0.0000,  0.00000,
   111.37,  0.0000,  0.00000,
   113.62,  0.0000,  0.00000,
   115.87,  0.0000,  0.00000,
   118.12,  0.0000,  0.00000,
   120.37,  0.0000,  0.00000,
   122.62,  0.0000,  0.00000,
   124.87,  0.0000,  0.00000,
   127.12,  0.0000,  0.00000,
   129.37,  0.0000,  0.00000,
   131.62,  0.0000,  0.00000,
   133.87,  0.0000,  0.00000,
   136.12,  0.0000,  0.00000,
   138.37,  0.0000,  0.00000,
   140.62,  0.0000,  0.00000,
   142.87,  0.0000,  0.00000,
   145.12,  0.0000,  0.00000,
   147.37,  0.0000,  0.00000,
   149.62,  0.0000,  0.00000,
   151.87,  0.0000,  0.00000,
   154.12,  0.0000,  0.00000,
   156.37,  0.0000,  0.00000,
   158.62,  0.0000,  0.00000,
   160.87,  0.0000,  0.00000,
   163.12,  0.0000,  0.00000,
   165.37,  0.0000,  0.00000,
   167.62,  0.0000,  0.00000,
   169.87,  0.0000,  0.00000,
   172.12,  0.0000,  0.00000,
   174.37,  0.0000,  0.00000,
   176.62,  0.0000,  0.00000,
   178.87,  0.0000,  0.00000,
   181.12,  0.0000,  0.00000,
   183.37,  0.0000,  0.00000,
   185.62,  0.0000,  0.00000,
   187.87,  0.0000,  0.00000,
   190.12,  0.0000,  0.00000,
   192.37,  0.0000,  0.00000,
   194.62,  0.0000,  0.00000,
   196.87,  0.0000,  0.00000,
   199.12,  0.0000,  0.00000,
   201.37,  0.0000,  0.00000,
   203.62,  0.0000,  0.00000,
   205.87,  0.0000,  0.00000,
   208.12,  0.0000,  0.00000,
   210.37,  0.0000,  0.00000,
   212.62,  0.0000,  0.00000,
   214.87,  0.0000,  0.00000,
   217.12,  0.0000,  0.00000,
   219.37,  0.0000,  0.00000,
   221.62,  0.0000,  0.00000,
   223.87,  0.0000,  0.00000,
  };	 
	 
  return dd[i][k];
}

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

				// create a 2D histogram to define the range

   TH2F *hr = new TH2F("hr","Gap and Skin",
		       1,0.,60,1,0,5000);

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

   c1->GetFrame()->SetFillColor(20);
   c1->GetFrame()->SetBorderSize(10);

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

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

   //   TText* text = new TText(40,4000,"");
   // text->Draw();     
  
}