// translated by Rudolph Magyar from cdf_pi_e_vs_p.top on 7/15/99 //
float dat1(int i, int k) {
  float dd[12][4] = {
  0.5000000,      0.6360119,      0.0134859,  .5,
  1.5000000,      0.6365577,      0.0058148,  .5,
  2.5000000,      0.6626099,      0.0100115,  .5,
  3.5000000,      0.7113768,      0.0178078,  .5,
  4.5000000,      0.7296090,      0.0313966,  .5,
  5.5000000,      0.7199395,      0.0347572,  .5,
  6.5000000,      0.8224602,      0.0598570,  .5,
  7.5000000,      0.7226542,      0.0350803,  .5,
  8.5000000,      0.7603319,      0.0428771,  .5,
  9.5000000,      0.7688834,      0.0708318,  .5,
 10.50000,        0.8524919,      0.03334248, .5, 
 11.50000,        0.8559219,      0.06122317, .5,
  };
  return dd[i][k];
}
float dat2(int i, int k) {
  float dd[6][4] = {
    7.12,    0.833,0.,  0.03,
    9.59,    0.844,0.,  0.03, 
    25.30,   0.854,0.,  0.003,
    50.26,   0.870,0.,  0.002,
    151.1,   0.854,0.,  0.002,
    220.1,   0.850,0.,  0.002,
  };

  return dd[i][k];
}

float dat3(int i, int k) {

  float dd[292][2] = {
  0.2000000,      0.5753636,    
  0.5000000,      0.5753636,   
   1.000000,      0.5790000,    
   1.500000,      0.5885000,    
   2.000000,      0.5980000,   
   2.500000,      0.6175000,    
   3.000000,      0.6370000,    
   3.500000,      0.6590000,    
   4.000000,      0.6810000,    
   4.500000,      0.7005000,    
   5.000000,      0.7200000,    
   5.500000,      0.7276887,    
   6.000000,      0.7357522,    
   6.500000,      0.7441835,    
   7.000000,      0.7529710,    
   7.500000,      0.7579316,    
   8.000000,      0.7632117,    
   8.500000,      0.7687868,    
   9.000000,      0.7746286,    
   9.500000,      0.7807057,    
   10.00000,      0.7869835,    
   10.50000,      0.7918916,    
   11.00000,      0.7969244,    
   11.50000,      0.8020411,    
   12.00000,      0.8072000,    
  // JOIN
  0.2000000,      0.6746364,    
  0.5000000,      0.6746364,    
   1.000000,      0.6710000,    
   1.500000,      0.6765000,    
   2.000000,      0.6820000,    
   2.500000,      0.6975000,    
   3.000000,      0.7130000,    
   3.500000,      0.7310000,    
   4.000000,      0.7490000,    
   4.500000,      0.7645000,    
   5.000000,      0.7800000,    
   5.500000,      0.7880566,    
   6.000000,      0.7957103,    
   6.500000,      0.8029402,    
   7.000000,      0.8097320,    
   7.500000,      0.8119808,    
   8.000000,      0.8137702,    
   8.500000,      0.8151023,    
   9.000000,      0.8159868,    
   9.500000,      0.8164406,    
   10.00000,      0.8164876,    
   10.50000,      0.8145048,    
   11.00000,      0.8122635,    
   11.50000,      0.8098115,    
   12.00000,      0.8072000,    
   12.50000,      0.8090000,    
  // JOIN
   5.000000,      0.7500000,    
   5.500000,      0.7541411,    
   6.000000,      0.7583244,    
   6.500000,      0.7625916,    
   7.000000,      0.7669835,    
   7.500000,      0.7673724,    
   8.000000,      0.7679620,    
   8.500000,      0.7687868,    
   9.000000,      0.7698784,    
   9.500000,      0.7712650,    
   10.00000,      0.7729709,    
   10.50000,      0.7734836,    
   11.00000,      0.7743521,    
   11.50000,      0.7755886,    
   12.00000,      0.7772000,    
   12.50000,      0.7790000,    
   13.00000,      0.7808000,    
   13.50000,      0.7826000,    
   14.00000,      0.7844000,    
   14.50000,      0.7862000,    
   15.00000,      0.7880000,    
   15.50000,      0.7898000,    
   16.00000,      0.7916000,    
   16.50000,      0.7934000,    
   17.00000,      0.7952000,    
   17.50000,      0.7970000,    
   18.00000,      0.7988000,    
   18.50000,      0.8006000,    
   19.00000,      0.8024000,    
   19.50000,      0.8042000,    
   20.00000,      0.8060001,    
   20.50000,      0.8078000,    
   21.00000,      0.8096000,    
   21.50000,      0.8114001,    
   22.00000,      0.8132000,    
   22.50000,      0.8150000,    
   23.00000,      0.8168000,    
   23.50000,      0.8186000,    
   24.00000,      0.8204000,    
   24.50000,      0.8222000,    
   25.00000,      0.8240000,    
   25.50000,      0.8245199,    
   26.00000,      0.8250400,    
   26.50000,      0.8255599,    
   27.00000,      0.8260800,    
   27.50000,      0.8266000,    
   28.00000,      0.8271200,    
   28.50000,      0.8276400,    
   29.00000,      0.8281600,    
   29.50000,      0.8286800,    
   30.00000,      0.8292000,    
   30.50000,      0.8297200,    
   31.00000,      0.8302400,    
   31.50000,      0.8307600,    
   32.00000,      0.8312800,    
   32.50000,      0.8318000,    
   33.00000,      0.8323200,    
   33.50000,      0.8328400,    
   34.00000,      0.8333600,    
   34.50000,      0.8338800,    
   35.00000,      0.8344000,    
   35.50000,      0.8349200,    
   36.00000,      0.8354400,    
   36.50000,      0.8359600,    
   37.00000,      0.8364800,    
   37.50000,      0.8370000,    
   38.00000,      0.8375200,    
   38.50000,      0.8380401,    
   39.00000,      0.8385600,    
   39.50000,      0.8390800,    
   40.00000,      0.8396000,    
   40.50000,      0.8401200,    
   41.00000,      0.8406400,    
   41.50000,      0.8411599,    
   42.00000,      0.8416800,    
   42.50000,      0.8422000,    
   43.00000,      0.8427200,    
   43.50000,      0.8432401,    
   44.00000,      0.8437601,    
   44.50000,      0.8442800,    
   45.00000,      0.8448000,    
   45.50000,      0.8453200,    
   46.00000,      0.8458400,    
   46.50000,      0.8463600,    
   47.00000,      0.8468800,    
   47.50000,      0.8474000,    
   48.00000,      0.8479200,    
   48.50000,      0.8484401,    
   49.00000,      0.8489600,    
   49.50000,      0.8494800,    
   50.00000,      0.8500000,    
   52.00000,      0.8494800,    
   54.00000,      0.8489600,    
   56.00000,      0.8484401,    
   58.00000,      0.8479200,    
   60.00000,      0.8474000,    
   62.00000,      0.8468800,    
   64.00000,      0.8463600,    
   66.00000,      0.8458400,    
   68.00000,      0.8453200,    
   70.00000,      0.8448000,    
   72.00000,      0.8442800,    
   74.00000,      0.8437600,    
   76.00000,      0.8432400,    
   78.00000,      0.8427200,    
   80.00000,      0.8422001,    
   82.00000,      0.8416800,    
   84.00000,      0.8411600,    
   86.00000,      0.8406401,    
   88.00000,      0.8401201,    
   90.00000,      0.8396000,    
   92.00000,      0.8390800,    
   94.00000,      0.8385600,    
   96.00000,      0.8380401,    
   98.00000,      0.8375200,    
   99.50000,      0.8371300,    
   100.0000,      0.8370001,    
  //   150.,  0.824,
  //   227.,  0.820,
  // JOIN dot
   5.000000,      0.7500000,    
   5.500000,      0.7609429,    
   6.000000,      0.7720094,    
   6.500000,      0.7831515,    
   7.000000,      0.7943181,    
   7.500000,      0.8013337,    
   8.000000,      0.8082809,    
   8.500000,      0.8151023,    
   9.000000,      0.8217399,    
   9.500000,      0.8281366,    
   10.00000,      0.8342368,    
   10.50000,      0.8381079,    
   11.00000,      0.8415573,    
   11.50000,      0.8445522,    
   12.00000,      0.8470666,    
   12.50000,      0.8493333,    
   13.00000,      0.8516000,    
   13.50000,      0.8538667,    
   14.00000,      0.8561333,    
   14.50000,      0.8584000,    
   15.00000,      0.8606667,    
   15.50000,      0.8629333,    
   16.00000,      0.8652000,    
   16.50000,      0.8674667,    
   17.00000,      0.8697333,    
   17.50000,      0.8720000,    
   18.00000,      0.8742667,    
   18.50000,      0.8765333,    
   19.00000,      0.8788000,    
   19.50000,      0.8810666,    
   20.00000,      0.8833333,    
   20.50000,      0.8856000,    
   21.00000,      0.8878667,    
   21.50000,      0.8901333,    
   22.00000,      0.8924000,    
   22.50000,      0.8946667,    
   23.00000,      0.8969334,    
   23.50000,      0.8992000,    
   24.00000,      0.9014667,    
   24.50000,      0.9037334,    
   25.00000,      0.9060000,    
   25.50000,      0.9064200,    
   26.00000,      0.9068400,    
   26.50000,      0.9072601,    
   27.00000,      0.9076800,    
   27.50000,      0.9081000,    
   28.00000,      0.9085200,    
   28.50000,      0.9089400,    
   29.00000,      0.9093600,    
   29.50000,      0.9097800,    
   30.00000,      0.9102000,    
   30.50000,      0.9106200,    
   31.00000,      0.9110400,    
   31.50000,      0.9114600,    
   32.00000,      0.9118800,    
   32.50000,      0.9123000,    
   33.00000,      0.9127200,    
   33.50000,      0.9131400,    
   34.00000,      0.9135600,    
   34.50000,      0.9139800,    
   35.00000,      0.9143999,    
   35.50000,      0.9148201,    
   36.00000,      0.9152400,    
   36.50000,      0.9156601,    
   37.00000,      0.9160799,    
   37.50000,      0.9165000,    
   38.00000,      0.9169200,    
   38.50000,      0.9173400,    
   39.00000,      0.9177600,    
   39.50000,      0.9181799,    
   40.00000,      0.9186000,    
   40.50000,      0.9190200,    
   41.00000,      0.9194400,    
   41.50000,      0.9198600,    
   42.00000,      0.9202800,    
   42.50000,      0.9207000,    
   43.00000,      0.9211199,    
   43.50000,      0.9215400,    
   44.00000,      0.9219599,    
   44.50000,      0.9223800,    
   45.00000,      0.9228000,    
   45.50000,      0.9232200,    
   46.00000,      0.9236401,    
   46.50000,      0.9240600,    
   47.00000,      0.9244800,    
   47.50000,      0.9249000,    
   48.00000,      0.9253200,    
   48.50000,      0.9257401,    
   49.00000,      0.9261600,    
   49.50000,      0.9265800,    
   50.00000,      0.9270000,    
   51.50000,      0.9267600,    
   53.00000,      0.9265200,    
   54.50000,      0.9262800,    
   56.00000,      0.9260400,    
   57.50000,      0.9258000,    
   59.00000,      0.9255600,    
   60.50000,      0.9253200,    
   62.00000,      0.9250800,    
   63.50000,      0.9248400,    
   65.00000,      0.9246000,    
   66.50000,      0.9243600,    
   68.00000,      0.9241200,    
   69.50000,      0.9238799,    
   71.00000,      0.9236400,    
   72.50000,      0.9234000,    
   74.00000,      0.9231600,    
   75.50000,      0.9229200,    
   77.00000,      0.9226800,    
   78.50000,      0.9224400,    
   80.00000,      0.9222000,    
   81.50000,      0.9219600,    
   83.00000,      0.9217199,    
   84.50000,      0.9214800,    
   86.00000,      0.9212400,    
   87.50000,      0.9210000,    
   89.00000,      0.9207600,    
   90.50000,      0.9205200,    
   92.00000,      0.9202800,    
   93.50000,      0.9200400,    
   95.00000,      0.9198000,    
   96.50000,      0.9195600,    
   98.00000,      0.9193200,    
   100.0000,      0.9190000,    
  // 150.,  0.911,
  // 227.,  0.907,
  // JOIN dot
  };

  return dd[i][k];
}

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

				// create a 2D histogram to define the range

   TH2F *hr = new TH2F("hr","Central Cal. Response to Single Pions",
		       1,0.25,240,1,0,1.4);

   hr->SetXTitle("POP1, GeV/c");
   hr->SetYTitle("E0P1 / P0P1");
   c1->SetLogx();
   hr->Draw();

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

   float x[1000], y[1000], ex[1000], ey[1000];
					// create first graph

   int n1 = 12;
   for (int i=0; i<n1; i++) {
     x [i] = dat1(i,0);
     y [i] = dat1(i,1);
     ex[i] = dat1(i,3);
     ey[i] = dat1(i,2);
     //Notice the order x,y,dy,dx
   }

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

   TText* text = new TText(0.01,1.2,"4 Collision Data (1991 ana?)");
   text->Draw();  
   TText* text = new TText(0.01,1.0,"P203 sys. error (Blue?)");
   text->Draw();  
   TText* text = new TText(0.01,0.2,"Testbeam sys. error (Red?)");
   text->Draw();  
   

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

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

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


   int n3 = 292;
   for (int i=0; i<n3; i++) {
     x[i] = dat3(i,0);
     y[i] = dat3(i,1);  
     ex[i] = 0;
     ey[i] = 0;
   }

   gr3 = new TGraphErrors(n3,x,y,ex,ey);
   gr3->SetMarkerColor(kGreen);
   gr3->SetMarkerStyle(20);
   gr3->SetLineColor(kGreen);
   gr3->Draw("L");

   TText* text = new TText(4.,0.22,"(Green?)");
   text->Draw();     

}