//transcribed by Rudolph Magyar from tb_57_pi_mip_peak.C om 7/16/99 //

float dat1(int i, int k) {

  float dd[100][2] = {
   0.00503, 0,
   0.0150,  0,
   0.0250,  0,
   0.0350,  0,
   0.0450,  0,
   0.0550,  0,
   0.0650,  0,
   0.0750,  0,
   0.0850,  0,
   0.0950,  0,
   0.1050,  0,
   0.1150,  0,
   0.1250,  0,
   0.1350,  0,
   0.1450,  0,
   0.1550,  0,
   0.1650,  0,
   0.1750,  0,
   0.1850,  0,
   0.1950, 5,
   0.2050, 17.000,
   0.2150, 43.000,
   0.2250, 89.000,
   0.2350, 141.00,
   0.2450, 178.00,
   0.2550, 232.00,
   0.2650, 262.00,
   0.2750, 307.00,
   0.2850, 319.00,
   0.2950, 292.00,
   0.3050, 285.00,
   0.3150, 248.00,
   0.3250, 253.00,
   0.3350, 221.00,
   0.3450, 226.00,
   0.3550, 191.00,
   0.3650, 180.00,
   0.3750, 181.00,
   0.3850, 131.00,
   0.3950, 157.00,
   0.4050, 128.00,
   0.4150, 120.00,
   0.4250, 120.00,
   0.4350, 103.00,
   0.4450, 88.000,
   0.4550, 72.000,
   0.4650, 83.000,
   0.4750, 74.000,
   0.4850, 60.000,
   0.4950, 66.000,
   0.5050, 58.000,
   0.5150, 55.000,
   0.5250, 50.000,
   0.5350, 59.000,
   0.5450, 55.000,
   0.5550, 33.000,
   0.5650, 47.000,
   0.5750, 38.000,
   0.5850, 33.000,
   0.5950, 28.000,
   0.6050, 28.000,
   0.6150, 31.000,
   0.6250, 25.000,
   0.6350, 24.000,
   0.6450, 32.000,
   0.6550, 20.000,
   0.6650, 26.000,
   0.6750, 23.000,
   0.6850, 21.000,
   0.6950, 19.000,
   0.7050, 17.000,
   0.7150, 17.000,
   0.7250, 19.000,
   0.7350, 24.000,
   0.7450, 15.000,
   0.7550, 17.000,
   0.7650, 17.000,
   0.7750, 16.000,
   0.7850, 14.000,
   0.7950, 22.000,
   0.8050, 8.00,
   0.8150, 9.00,
   0.8250, 12.0,
   0.8350, 15.0,
   0.8450, 7.00,
   0.8550, 10.0,
   0.8650, 7.00,
   0.8750, 13.0,
   0.8850, 7.00,
   0.8950, 7.00,
   0.9050, 11.0,
   0.9150, 9.00,
   0.9250, 7.00,
   0.9350, 13.0,
   0.9450, 16.0,
   0.9550, 13.0,
   0.9650, 12.0,
   0.9750, 4.00,
   0.9850, 6.00,
   0.9950, 11.0,
  };

  return dd[i][k];
}

float dat2(int i, int k) {

  float dd[100][2] = {
   0.005,    7.00,
   0.01500,  5.00,
   0.02500,  1.00,
   0.03500,  8.00,
   0.04500,  7.00,
   0.05500,  11.0,
   0.06500,  7.00,
   0.07500,  12.0,
   0.08500,  13.0,
   0.09500,  12.0,
   .105000,  10.0,
   .115000,  8.00,
   .125000,  16.0,
   .135000,  11.0,
   .145000,  15.0,
   .155000,  26.0,
   .165000,  36.0,
   .175000,  30.0,
   .185000,  41.0,
   .195000,  55.0,
   .205000,  75.0,
   .215000,  93.0,
   .225000,  121.,
   .235000,  117.,
   .245000,  135.,
   .255000,  162.,
   .265000,  174.,
   .275000,  190.,
   .285000,  216.,
   .295000,  226.,
   .305000,  200.,
   .315000,  205.,
   .325000,  193.,
   .335000,  180.,
   .345000,  181.,
   .355000,  155.,
   .365000,  142.,
   .375000,  118.,
   .385000,  117.,
   .395000,  144.,
   .405000,  120.,
   .415000,  96.0,
   .425000,  95.0,
   .435000,  96.0,
   .445000,  77.0,
   .455000,  75.0,
   .465000,  72.0,
   .475000,  56.0,
   .485000,  46.0,
   .495000,  55.0,
   .505000,  49.0,
   .515000,  47.0,
   .525000,  44.0,
   .535000,  52.0,
   .545000,  54.0,
   .555000,  45.0,
   .565000,  38.0,
   .575000,  41.0,
   .585000,  34.0,
   .595000,  32.0,
   .605000,  38.0,
   .615000,  25.0,
   .625000,  27.0,
   .635000,  26.0,
   .645000,  21.0,
   .655000,  21.0,
   .665000,  18.0,
   .675000,  19.0,
   .685000,  13.0,
   .695000,  29.0,
   .705000,  18.0,
   .715000,  23.0,
   .725000,  18.0,
   .735000,  14.0,
   .745000,  18.0,
   .755000,  14.0,
   .765000,  19.0,
   .775000,  13.0,
   .785000,  13.0,
   .795000,  14.0,
   .805000,  25.0,
   .815000,  14.0,
   .825000,  13.0,
   .835000,  13.0,
   .845000,  14.0,
   .855000,  11.0,
   .865000,  14.0,
   .875000,  8.00,
   .885000,  11.0,
   .895000,  9.00,
   .905000,  9.00,
   .915000,  10.0,
   .925000,  8.00,
   .935000,  10.0,
   .945000,  10.0,
   .955000,  13.0,
   .965000,  7.00,
   .975000,  12.0,
   .985000,  7.00,
   .995000,  6.00,

  };

  return dd[i][k];
}

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

				// create a 2D histogram to define the range

   TH2F *hr = new TH2F("hr","Min Ionizing Peak in CEM for Pions",
		       1,0.,1.0,1,0,400);

   hr->SetXTitle("Energy in GeV/c");
   hr->Draw();

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

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

					// create first 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] = 0;
   }

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

   TText* text = new TText(0.5,250,"57.1 GeV Testbeam (Red)");
   text->Draw();     

   int n2 = 100;
   for (int i=0; i<n2; i++) {
     x [i] = dat2(i,0);
     y [i] = dat2(i,1);
     ex[i] = 0;
     ey[i] = 0;
   }

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

   TText* text = new TText(0.5,200,"QFL(Blue)");
   text->Draw();     
}