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

  float dd[141][3] = {
       4.7500000,      2.0000000,      1.4142135,
       5.2500000,      0.0000000,      0.0000000,
       5.7500000,      0.0000000,      0.0000000,
       6.2500000,      4.0000000,      2.0000000,
       6.7500000,      1.0000000,      1.0000000,
       7.2500000,      1.0000000,      1.0000000,
       7.7500000,      2.0000000,      1.4142135,
       8.2500000,      1.0000000,      1.0000000,
       8.7500000,      1.0000000,      1.0000000,
       9.2500000,      0.0000000,      0.0000000,
       9.7500000,      1.0000000,      1.0000000,
      10.2500000,      0.0000000,      0.0000000,
      10.7500000,      0.0000000,      0.0000000,
      11.2500000,      1.0000000,      1.0000000,
      11.7500000,      2.0000000,      1.4142135,
      12.2500000,      0.0000000,      0.0000000,
      12.7500000,      2.0000000,      1.4142135,
      13.2500000,      0.0000000,      0.0000000,
      13.7500000,      0.0000000,      0.0000000,
      14.2500000,      0.0000000,      0.0000000,
      14.7500000,      1.0000000,      1.0000000,
      15.2500000,      2.0000000,      1.4142135,
      15.7500000,      2.0000000,      1.4142135,
      16.2500000,      1.0000000,      1.0000000,
      16.7500000,      0.0000000,      0.0000000,
      17.2500000,      2.0000000,      1.4142135,
      17.7500000,      2.0000000,      1.4142135,
      18.2500000,      2.0000000,      1.4142135,
      18.7500000,      0.0000000,      0.0000000,
      19.2500000,      1.0000000,      1.0000000,
      19.7500000,      1.0000000,      1.0000000,
      20.2500000,      0.0000000,      0.0000000,
      20.7500000,      2.0000000,      1.4142135,
      21.2500000,      1.0000000,      1.0000000,
      21.7500000,      0.0000000,      0.0000000,
      22.2500000,      3.0000000,      1.7320508,
      22.7500000,      1.0000000,      1.0000000,
      23.2500000,      0.0000000,      0.0000000,
      23.7500000,      0.0000000,      0.0000000,
      24.2500000,      1.0000000,      1.0000000,
      24.7500000,      1.0000000,      1.0000000,
      25.2500000,      4.0000000,      2.0000000,
      25.7500000,      1.0000000,      1.0000000,
      26.2500000,      2.0000000,      1.4142135,
      26.7500000,      5.0000000,      2.2360680,
      27.2500000,      3.0000000,      1.7320508,
      27.7500000,      1.0000000,      1.0000000,
      28.2500000,      2.0000000,      1.4142135,
      28.7500000,      5.0000000,      2.2360680,
      29.2500000,      2.0000000,      1.4142135,
      29.7500000,      0.0000000,      0.0000000,
      30.2500000,      0.0000000,      0.0000000,
      30.7500000,      2.0000000,      1.4142135,
      31.2500000,      0.0000000,      0.0000000,
      31.7500000,      3.0000000,      1.7320508,
      32.2500000,      3.0000000,      1.7320508,
      32.7500000,      8.0000000,      2.8284271,
      33.2500000,      1.0000000,      1.0000000,
      33.7500000,      4.0000000,      2.0000000,
      34.2500000,      7.0000000,      2.6457512,
      34.7500000,      4.0000000,      2.0000000,
      35.2500000,      5.0000000,      2.2360680,
      35.7500000,      9.0000000,      3.0000000,
      36.2500000,      6.0000000,      2.4494898,
      36.7500000,      2.0000000,      1.4142135,
      37.2500000,      5.0000000,      2.2360680,
      37.7500000,      9.0000000,      3.0000000,
      38.2500000,      5.0000000,      2.2360680,
      38.7500000,      8.0000000,      2.8284271,
      39.2500000,     13.0000000,      3.6055512,
      39.7500000,      5.0000000,      2.2360680,
      40.2500000,     11.0000000,      3.3166249,
      40.7500000,     10.0000000,      3.1622777,
      41.2500000,     18.0000000,      4.2426405,
      41.7500000,     11.0000000,      3.3166249,
      42.2500000,     21.0000000,      4.5825758,
      42.7500000,     19.0000000,      4.3588991,
      43.2500000,     19.0000000,      4.3588991,
      43.7500000,     21.0000000,      4.5825758,
      44.2500000,     29.0000000,      5.3851647,
      44.7500000,     22.0000000,      4.6904159,
      45.2500000,     30.0000000,      5.4772258,
      45.7500000,     29.0000000,      5.3851647,
      46.2500000,     36.0000000,      6.0000000,
      46.7500000,     47.0000000,      6.8556547,
      47.2500000,     35.0000000,      5.9160800,
      47.7500000,     29.0000000,      5.3851647,
      48.2500000,     38.0000000,      6.1644139,
      48.7500000,     33.0000000,      5.7445626,
      49.2500000,     38.0000000,      6.1644139,
      49.7500000,     38.0000000,      6.1644139,
      50.2500000,     41.0000000,      6.4031243,
      50.7500000,     38.0000000,      6.1644139,
      51.2500000,     44.0000000,      6.6332498,
      51.7500000,     49.0000000,      7.0000000,
      52.2500000,     33.0000000,      5.7445626,
      52.7500000,     42.0000000,      6.4807405,
      53.2500000,     34.0000000,      5.8309517,
      53.7500000,     32.0000000,      5.6568542,
      54.2500000,     39.0000000,      6.2449980,
      54.7500000,     32.0000000,      5.6568542,
      55.2500000,     26.0000000,      5.0990195,
      55.7500000,     29.0000000,      5.3851647,
      56.2500000,     21.0000000,      4.5825758,
      56.7500000,     25.0000000,      5.0000000,
      57.2500000,     17.0000000,      4.1231055,
      57.7500000,     11.0000000,      3.3166249,
      58.2500000,     13.0000000,      3.6055512,
      58.7500000,     17.0000000,      4.1231055,
      59.2500000,     13.0000000,      3.6055512,
      59.7500000,     13.0000000,      3.6055512,
      60.2500000,      8.0000000,      2.8284271,
      60.7500000,      6.0000000,      2.4494898,
      61.2500000,      7.0000000,      2.6457512,
      61.7500000,      6.0000000,      2.4494898,
      62.2500000,      4.0000000,      2.0000000,
      62.7500000,      5.0000000,      2.2360680,
      63.2500000,      2.0000000,      1.4142135,
      63.7500000,      3.0000000,      1.7320508,
      64.2500000,      2.0000000,      1.4142135,
      64.7500000,      3.0000000,      1.7320508,
      65.2500000,      3.0000000,      1.7320508,
      65.7500000,      1.0000000,      1.0000000,
      66.2500000,      0.0000000,      0.0000000,
      66.7500000,      1.0000000,      1.0000000,
      67.2500000,      0.0000000,      0.0000000,
      67.7500000,      0.0000000,      0.0000000,
      68.2500000,      0.0000000,      0.0000000,
      68.7500000,      0.0000000,      0.0000000,
      69.2500000,      0.0000000,      0.0000000,
      69.7500000,      0.0000000,      0.0000000,
      70.2500000,      0.0000000,      0.0000000,
      70.7500000,      0.0000000,      0.0000000,
      71.2500000,      0.0000000,      0.0000000,
      71.7500000,      0.0000000,      0.0000000,
      72.2500000,      0.0000000,      0.0000000,
      72.7500000,      1.0000000,      1.0000000,
      73.2500000,      0.0000000,      0.0000000,
      73.7500000,      0.0000000,      0.0000000,
      74.2500000,      0.0000000,      0.0000000,
      74.7500000,      0.0000000,      0.0000000,
  };

  return dd[i][k];
}

float dat2(int i, int k) {

  float dd[200][3] = {
       0.5000000,      0.0000000,      0.0000000,
       1.5000000,      0.0000000,      0.0000000,
       2.5000000,      0.0000000,      0.0000000,
       3.5000000,      4.0000000,      2.0000000,
       4.5000000,      1.0000000,      1.0000000,
       5.5000000,      3.0000000,      1.7320508,
       6.5000000,      4.0000000,      2.0000000,
       7.5000000,      3.0000000,      1.7320508,
       8.5000000,      2.0000000,      1.4142135,
       9.5000000,      2.0000000,      1.4142135,
      10.5000000,      2.0000000,      1.4142135,
      11.5000000,      3.0000000,      1.7320508,
      12.5000000,      1.0000000,      1.0000000,
      13.5000000,      3.0000000,      1.7320508,
      14.5000000,      2.0000000,      1.4142135,
      15.5000000,      4.0000000,      2.0000000,
      16.5000000,      4.0000000,      2.0000000,
      17.5000000,      1.0000000,      1.0000000,
      18.5000000,      3.0000000,      1.7320508,
      19.5000000,      3.0000000,      1.7320508,
      20.5000000,      4.0000000,      2.0000000,
      21.5000000,      3.0000000,      1.7320508,
      22.5000000,      4.0000000,      2.0000000,
      23.5000000,      5.0000000,      2.2360680,
      24.5000000,      6.0000000,      2.4494898,
      25.5000000,      7.0000000,      2.6457512,
      26.5000000,     11.0000000,      3.3166249,
      27.5000000,     11.0000000,      3.3166249,
      28.5000000,      3.0000000,      1.7320508,
      29.5000000,     12.0000000,      3.4641016,
      30.5000000,      8.0000000,      2.8284271,
      31.5000000,     12.0000000,      3.4641016,
      32.5000000,     19.0000000,      4.3588991,
      33.5000000,     27.0000000,      5.1961522,
      34.5000000,     44.0000000,      6.6332498,
      35.5000000,     40.0000000,      6.3245554,
      36.5000000,     49.0000000,      7.0000000,
      37.5000000,     66.0000000,      8.1240387,
      38.5000000,    100.0000000,     10.0000000,
      39.5000000,    112.0000000,     10.5830050,
      40.5000000,    164.0000000,     12.8062487,
      41.5000000,    165.0000000,     12.8452330,
      42.5000000,    211.0000000,     14.5258389,
      43.5000000,    234.0000000,     15.2970581,
      44.5000000,    244.0000000,     15.6204996,
      45.5000000,    263.0000000,     16.2172756,
      46.5000000,    250.0000000,     15.8113880,
      47.5000000,    227.0000000,     15.0665188,
      48.5000000,    191.0000000,     13.8202753,
      49.5000000,    187.0000000,     13.6747942,
      50.5000000,    171.0000000,     13.0766964,
      51.5000000,    149.0000000,     12.2065554,
      52.5000000,    127.0000000,     11.2694273,
      53.5000000,     99.0000000,      9.9498739,
      54.5000000,     89.0000000,      9.4339809,
      55.5000000,     61.0000000,      7.8102498,
      56.5000000,     54.0000000,      7.3484693,
      57.5000000,     45.0000000,      6.7082038,
      58.5000000,     45.0000000,      6.7082038,
      59.5000000,     21.0000000,      4.5825758,
      60.5000000,     15.0000000,      3.8729835,
      61.5000000,      6.0000000,      2.4494898,
      62.5000000,      9.0000000,      3.0000000,
      63.5000000,      4.0000000,      2.0000000,
      64.5000000,      4.0000000,      2.0000000,
      65.5000000,      1.0000000,      1.0000000,
      66.5000000,      0.0000000,      0.0000000,
      67.5000000,      0.0000000,      0.0000000,
      68.5000000,      1.0000000,      1.0000000,
      69.5000000,      0.0000000,      0.0000000,
      70.5000000,      0.0000000,      0.0000000,
      71.5000000,      0.0000000,      0.0000000,
      72.5000000,      0.0000000,      0.0000000,
      73.5000000,      0.0000000,      0.0000000,
      74.5000000,      0.0000000,      0.0000000,
      75.5000000,      0.0000000,      0.0000000,
      76.5000000,      0.0000000,      0.0000000,
      77.5000000,      0.0000000,      0.0000000,
      78.5000000,      0.0000000,      0.0000000,
      79.5000000,      0.0000000,      0.0000000,
      80.5000000,      0.0000000,      0.0000000,
      81.5000000,      0.0000000,      0.0000000,
      82.5000000,      0.0000000,      0.0000000,
      83.5000000,      0.0000000,      0.0000000,
      84.5000000,      0.0000000,      0.0000000,
      85.5000000,      0.0000000,      0.0000000,
      86.5000000,      0.0000000,      0.0000000,
      87.5000000,      0.0000000,      0.0000000,
      88.5000000,      0.0000000,      0.0000000,
      89.5000000,      0.0000000,      0.0000000,
      90.5000000,      1.0000000,      1.0000000,
      91.5000000,      0.0000000,      0.0000000,
      92.5000000,      0.0000000,      0.0000000,
      93.5000000,      0.0000000,      0.0000000,
      94.5000000,      0.0000000,      0.0000000,
      95.5000000,      0.0000000,      0.0000000,
      96.5000000,      0.0000000,      0.0000000,
      97.5000000,      0.0000000,      0.0000000,
      98.5000000,      0.0000000,      0.0000000,
      99.5000000,      0.0000000,      0.0000000,
  };

  return dd[i][k];
}

float dat3(int i, int k) {

  float dd[104][3] = {
      4.7500000,      0.0000000,	
      6.2500000,      0.0000000,	
      6.7500000,      0.0000000,	
      7.2500000,      0.0000000,	
      7.7500000,      0.0000000,	
      8.2500000,      0.0000000,	
      8.7500000,      0.0000000,	
      9.7500000,      0.0000000,	
     11.2500000,      0.0000000,	
     11.7500000,      0.0000000,	
     12.7500000,      0.0000001,	
     14.7500000,      0.0000004,	
     15.2500000,      0.0000007,	
     15.7500000,      0.0000012,	
     16.2500000,      0.0000020,	
     17.2500000,      0.0000053,	
     17.7500000,      0.0000085,	
     18.2500000,      0.0000137,	
     19.2500000,      0.0000343,	
     19.7500000,      0.0000538,	
     20.7500000,      0.0001295,	
     21.2500000,      0.0001986,	
     22.2500000,      0.0004571,	
     22.7500000,      0.0006858,	
     24.2500000,      0.0022165,	
     24.7500000,      0.0032292,	
     25.2500000,      0.0046700,	
     25.7500000,      0.0067043,	
     26.2500000,      0.0095543,	
     26.7500000,      0.0135159,	
     27.2500000,      0.0189801,	
     27.7500000,      0.0264580,	
     28.2500000,      0.0366118,	
     28.7500000,      0.0502910,	
     29.2500000,      0.0685749,	
     30.7500000,      0.1663508,	
     31.7500000,      0.2894849,	
     32.2500000,      0.3776895,	
     32.7500000,      0.4891583,	
     33.2500000,      0.6288827,	
     33.7500000,      0.8025931,
     34.2500000,      1.0167797,
     34.7500000,      1.2786858,
     35.2500000,      1.5962706,
     35.7500000,      1.9781297,
     36.2500000,      2.4333727,
     36.7500000,      2.9714479,
     37.2500000,      3.6019130,
     37.7500000,      4.3341498,
     38.2500000,      5.1770253,
    //join dot       	       
     38.7500000,      6.1384993,
     39.2500000,      7.2251987,
     39.7500000,      8.4419537,
     40.2500000,      9.7913313,
     40.7500000,     11.2731733,
     41.2500000,     12.8841629,
     41.7500000,     14.6174574,
     42.2500000,     16.4623985,
     42.7500000,     18.4043274,
     43.2500000,     20.4245472,
     43.7500000,     22.5004158,
     44.2500000,     24.6056175,
     44.7500000,     26.7105980,
     45.2500000,     28.7831650,
     45.7500000,     30.7892494,
     46.2500000,     32.6937904,
     46.7500000,     34.4617310,
     47.2500000,     36.0590668,
     47.7500000,     37.4539413,
     48.2500000,     38.6176796,
     48.7500000,     39.5257759,
     49.2500000,     40.1587601,
     49.7500000,     40.5028687,
     50.2500000,     40.5505600,
     50.7500000,     40.3007927,
     51.2500000,     39.7590446,
     51.7500000,     38.9371223,
     52.2500000,     37.8527489,
     52.7500000,     36.5289001,
     53.2500000,     34.9930191,
     53.7500000,     33.2760544,
     54.2500000,     31.4114399,
     54.7500000,     29.4340134,
     55.2500000,     27.3789463,
     55.7500000,     25.2807293,
     56.2500000,     23.1722431,
     56.7500000,     21.0839596,
     57.2500000,     19.0432873,
     57.7500000,     17.0740776,
     58.2500000,     15.1963110,
     58.7500000,     13.4259415,
     59.2500000,     11.7748919,
     59.7500000,     10.2512007,
     60.2500000,      8.8592739,
     60.7500000,      7.6002374,
     61.2500000,      6.4723468,
    //join	       	       
     61.7500000,      5.4714451,
     62.2500000,      4.5914292,
     62.7500000,      3.8247178,
     63.2500000,      3.1626892,
     63.7500000,      2.5960867,
     64.2500000,      2.1153760,
     64.7500000,      1.7110453,
     65.2500000,      1.3738555,
     65.7500000,      1.0950305,
     66.7500000,      0.6804770,

  };

  return dd[i][k];
}

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

				// create a 2D histogram to define the range

   TH2F *hr = new TH2F("hr","Central Cal. Energy for 50 GeV Pions",
		       1,0.,100,1,0,150);

   hr->SetXTitle("E Cal, GeV/c");
   hr->Draw();

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

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

					// create graph
   int n1 = 141;
   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(20);
   gr1->Draw("P");

   TText* text = new TText(52,120,"Right CHA PMT, CEM Mips (Red)");
   text->Draw();     

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

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

   TText* text = new TText(52,130,"CEM+CHA, All Pions (Blue)");
   text->Draw();     
}