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

float dat1(int i, int k) {
  //EtotAllPions
  float dd[100][2] = {
   0.50000,  0.0000, //# of data pairs 100
   1.50000,  0.0000,
   2.50000,  0.0000,
   3.50000,  0.0000,
   4.50000,  0.0000,
   5.50000,  6.0000,
   6.50000,  8.0000,
   7.50000,  8.0000,
   8.50000,  7.0000,
   9.50000,  9.0000,
   10.5000,  11.000,
   11.5000,  6.0000,
   12.5000,  11.000,
   13.5000,  5.0000,
   14.5000,  7.0000,
   15.5000,  4.0000,
   16.5000,  9.0000,
   17.5000,  8.0000,
   18.5000,  9.0000,
   19.5000,  9.0000,
   20.5000,  6.0000,
   21.5000,  11.000,
   22.5000,  14.000,
   23.5000,  8.0000,
   24.5000,  13.000,
   25.5000,  14.000,
   26.5000,  17.000,
   27.5000,  17.000,
   28.5000,  17.000,
   29.5000,  22.000,
   30.5000,  24.000,
   31.5000,  27.000,
   32.5000,  30.000,
   33.5000,  34.000,
   34.5000,  42.000,
   35.5000,  39.000,
   36.5000,  44.000,
   37.5000,  58.000,
   38.5000,  66.000,
   39.5000,  77.000,
   40.5000,  108.00,
   41.5000,  109.00,
   42.5000,  164.00,
   43.5000,  247.00,
   44.5000,  290.00,
   45.5000,  371.00,
   46.5000,  469.00,
   47.5000,  581.00,
   48.5000,  619.00,
   49.5000,  747.00,
   50.5000,  797.00,
   51.5000,  835.00,
   52.5000,  885.00,
   53.5000,  917.00,
   54.5000,  830.00,
   55.5000,  798.00,
   56.5000,  742.00,
   57.5000,  647.00,
   58.5000,  569.00,
   59.5000,  523.00,
   60.5000,  431.00,
   61.5000,  353.00,
   62.5000,  349.00,
   63.5000,  227.00,
   64.5000,  224.00,
   65.5000,  176.00,
   66.5000,  116.00,
   67.5000,  110.00,
   68.5000,  99.000,
   69.5000,  52.000,
   70.5000,  40.000,
   71.5000,  29.000,
   72.5000,  27.000,
   73.5000,  9.0000,
   74.5000,  12.000,
   75.5000,  4.0000,
   76.5000,  3.0000,
   77.5000,  0.0000,
   78.5000,  0.0000,
   79.5000,  0.0000,
   80.5000,  1.0000,
   81.5000,  0.0000,
   82.5000,  0.0000,
   83.5000,  0.0000,
   84.5000,  0.0000,
   85.5000,  0.0000,
   86.5000,  0.0000,
   87.5000,  0.0000,
   88.5000,  0.0000,
   89.5000,  0.0000,
   90.5000,  0.0000,
   91.5000,  0.0000,
   92.5000,  0.0000,
   93.5000,  0.0000,
   94.5000,  0.0000,
   95.5000,  0.0000,
   96.5000,  0.0000,
   97.5000,  0.0000,
   98.5000,  0.0000,
   99.5000,  0.0000,
  };

  return dd[i][k];
}

float dat2(int i, int k) {

  float dd[100][2] = {
   1.00000,  0.0000, //# of data pairs 100
   3.00000,  0.0000,
   5.00000,  11.000,
   7.00000,  9.0000,
   9.00000,  10.000,
   11.0000,  15.000,
   13.0000,  10.000,
   15.0000,  14.000,
   17.0000,  18.000,
   19.0000,  15.000,
   21.0000,  14.000,
   23.0000,  8.0000,
   25.0000,  11.000,
   27.0000,  20.000,
   29.0000,  25.000,
   31.0000,  22.000,
   33.0000,  35.000,
   35.0000,  43.000,
   37.0000,  52.000,
   39.0000,  69.000,
   41.0000,  164.00,
   43.0000,  306.00,
   45.0000,  604.00,
   47.0000,  986.00,
   49.0000,  1331.0,
   51.0000,  1666.0,
   53.0000,  1866.0,
   55.0000,  1689.0,
   57.0000,  1515.0,
   59.0000,  1156.0,
   61.0000,  920.00,
   63.0000,  637.00,
   65.0000,  412.00,
   67.0000,  288.00,
   69.0000,  146.00,
   71.0000,  82.000,
   73.0000,  33.000,
   75.0000,  11.000,
   77.0000,  3.0000,
   79.0000,  0.0000,
   81.0000,  0.0000,
   83.0000,  0.0000,
   85.0000,  0.0000,
   87.0000,  0.0000,
   89.0000,  0.0000,
   91.0000,  0.0000,
   93.0000,  0.0000,
   95.0000,  0.0000,
   97.0000,  0.0000,
   99.0000,  0.0000,
   101.000,  0.0000,
   103.000,  0.0000,
   105.000,  0.0000,
   107.000,  0.0000,
   109.000,  0.0000,
   111.000,  0.0000,
   113.000,  0.0000,
   115.000,  0.0000,
   117.000,  0.0000,
   119.000,  0.0000,
   121.000,  0.0000,
   123.000,  0.0000,
   125.000,  0.0000,
   127.000,  0.0000,
   129.000,  0.0000,
   131.000,  0.0000,
   133.000,  0.0000,
   135.000,  0.0000,
   137.000,  0.0000,
   139.000,  0.0000,
   141.000,  0.0000,
   143.000,  0.0000,
   145.000,  0.0000,
   147.000,  0.0000,
   149.000,  0.0000,
   151.000,  0.0000,
   153.000,  0.0000,
   155.000,  0.0000,
   157.000,  0.0000,
   159.000,  0.0000,
   161.000,  0.0000,
   163.000,  0.0000,
   165.000,  0.0000,
   167.000,  0.0000,
   169.000,  0.0000,
   171.000,  0.0000,
   173.000,  0.0000,
   175.000,  0.0000,
   177.000,  0.0000,
   179.000,  0.0000,
   181.000,  0.0000,
   183.000,  0.0000,
   185.000,  0.0000,
   187.000,  0.0000,
   189.000,  0.0000,
   191.000,  0.0000,
   193.000,  0.0000,
   195.000,  0.0000,
   197.000,  0.0000,
   199.000,  0.0000,
    };

  return dd[i][k];
}
//end of data for the first plot

//EEM: All Pions

float dat3(int i, int k) {

  float dd[100][2] = {
   1.50000,  470.00, //# of data pairs 100
   2.50000,  297.00,
   3.50000,  240.00,
   4.50000,  219.00,
   5.50000,  222.00,
   6.50000,  243.00,
   7.50000,  230.00,
   8.50000,  207.00,
   9.50000,  230.00,
   10.5000,  212.00,
   11.5000,  240.00,
   12.5000,  226.00,
   13.5000,  211.00,
   14.5000,  212.00,
   15.5000,  221.00,
   16.5000,  226.00,
   17.5000,  245.00,
   18.5000,  230.00,
   19.5000,  236.00,
   20.5000,  234.00,
   21.5000,  243.00,
   22.5000,  263.00,
   23.5000,  218.00,
   24.5000,  231.00,
   25.5000,  237.00,
   26.5000,  223.00,
   27.5000,  204.00,
   28.5000,  210.00,
   29.5000,  193.00,
   30.5000,  168.00,
   31.5000,  172.00,
   32.5000,  147.00,
   33.5000,  142.00,
   34.5000,  140.00,
   35.5000,  116.00,
   36.5000,  111.00,
   37.5000,  100.00,
   38.5000,  93.000,
   39.5000,  78.000,
   40.5000,  67.000,
   41.5000,  64.000,
   42.5000,  63.000,
   43.5000,  60.000,
   44.5000,  42.000,
   45.5000,  35.000,
   46.5000,  25.000,
   47.5000,  32.000,
   48.5000,  22.000,
   49.5000,  25.000,
   50.5000,  24.000,
   51.5000,  11.000,
   52.5000,  10.000,
   53.5000,  8.0000,
   54.5000,  10.000,
   55.5000,  3.0000,
   56.5000,  3.0000,
   57.5000,  2.0000,
   58.5000,  1.0000,
   59.5000,  0.0000,
   60.5000,  0.0000,
   61.5000,  0.0000,
   62.5000,  0.0000,
   63.5000,  0.0000,
   64.5000,  0.0000,
   65.5000,  0.0000,
   66.5000,  0.0000,
   67.5000,  0.0000,
   68.5000,  0.0000,
   69.5000,  0.0000,
   70.5000,  0.0000,
   71.5000,  0.0000,
   72.5000,  0.0000,
   73.5000,  0.0000,
   74.5000,  0.0000,
   75.5000,  0.0000,
   76.5000,  0.0000,
   77.5000,  0.0000,
   78.5000,  0.0000,
   79.5000,  0.0000,
   80.5000,  0.0000,
   81.5000,  0.0000,
   82.5000,  0.0000,
   83.5000,  0.0000,
   84.5000,  0.0000,
   85.5000,  0.0000,
   86.5000,  0.0000,
   87.5000,  0.0000,
   88.5000,  0.0000,
   89.5000,  0.0000,
   90.5000,  0.0000,
   91.5000,  0.0000,
   92.5000,  0.0000,
   93.5000,  0.0000,
   94.5000,  0.0000,
   95.5000,  0.0000,
   96.5000,  0.0000,
   97.5000,  0.0000,
   98.5000,  0.0000,
   99.5000,  0.0000,
   100.500,  0.0000,
  };

  return dd[i][k];
}

float dat4(int i, int k) {

  float dd[100][2] = {   
   1.00000,  6300.0, //# of data pairs 100
   3.00000,  471.00,
   5.00000,  388.00,
   7.00000,  452.00,
   9.00000,  442.00,
   11.0000,  512.00,
   13.0000,  525.00,
   15.0000,  521.00,
   17.0000,  473.00,
   19.0000,  439.00,
   21.0000,  427.00,
   23.0000,  397.00,
   25.0000,  396.00,
   27.0000,  425.00,
   29.0000,  370.00,
   31.0000,  332.00,
   33.0000,  300.00,
   35.0000,  238.00,
   37.0000,  203.00,
   39.0000,  169.00,
   41.0000,  132.00,
   43.0000,  87.000,
   45.0000,  83.000,
   47.0000,  45.000,
   49.0000,  45.000,
   51.0000,  17.000,
   53.0000,  24.000,
   55.0000,  3.0000,
   57.0000,  0.0000,
   59.0000,  0.0000,
   61.0000,  0.0000,
   63.0000,  0.0000,
   65.0000,  0.0000,
   67.0000,  0.0000,
   69.0000,  0.0000,
   71.0000,  0.0000,
   73.0000,  0.0000,
   75.0000,  0.0000,
   77.0000,  0.0000,
   79.0000,  0.0000,
   81.0000,  0.0000,
   83.0000,  0.0000,
   85.0000,  0.0000,
   87.0000,  0.0000,
   89.0000,  0.0000,
   91.0000,  0.0000,
   93.0000,  0.0000,
   95.0000,  0.0000,
   97.0000,  0.0000,
   99.0000,  0.0000,
   101.000,  0.0000,
   103.000,  0.0000,
   105.000,  0.0000,
   107.000,  0.0000,
   109.000,  0.0000,
   111.000,  0.0000,
   113.000,  0.0000,
   115.000,  0.0000,
   117.000,  0.0000,
   119.000,  0.0000,
   121.000,  0.0000,
   123.000,  0.0000,
   125.000,  0.0000,
   127.000,  0.0000,
   129.000,  0.0000,
   131.000,  0.0000,
   133.000,  0.0000,
   135.000,  0.0000,
   137.000,  0.0000,
   139.000,  0.0000,
   141.000,  0.0000,
   143.000,  0.0000,
   145.000,  0.0000,
   147.000,  0.0000,
   149.000,  0.0000,
   151.000,  0.0000,
   153.000,  0.0000,
   155.000,  0.0000,
   157.000,  0.0000,
   159.000,  0.0000,
   161.000,  0.0000,
   163.000,  0.0000,
   165.000,  0.0000,
   167.000,  0.0000,
   169.000,  0.0000,
   171.000,  0.0000,
   173.000,  0.0000,
   175.000,  0.0000,
   177.000,  0.0000,
   179.000,  0.0000,
   181.000,  0.0000,
   183.000,  0.0000,
   185.000,  0.0000,
   187.000,  0.0000,
   189.000,  0.0000,
   191.000,  0.0000,
   193.000,  0.0000,
   195.000,  0.0000,
   197.000,  0.0000,
   199.000,  0.0000,
  };

  return dd[i][k];
}
//end of second graph data

//EHad: All Pions
float dat5(int i, int k) {

  float dd[100][2] = {   
   0.5000 , 6.00000, //# of data pairs 100
   1.5000 , 8.00000,
   2.5000 , 14.0000,
   3.5000 , 39.0000,
   4.5000 , 36.0000,
   5.5000 , 41.0000,
   6.5000 , 59.0000,
   7.5000 , 71.0000,
   8.5000 , 90.0000,
   9.5000 , 95.0000,
   10.500 , 95.0000,
   11.500 , 110.000,
   12.500 , 140.000,
   13.500 , 145.000,
   14.500 , 135.000,
   15.500 , 169.000,
   16.500 , 168.000,
   17.500 , 179.000,
   18.500 , 200.000,
   19.500 , 198.000,
   20.500 , 201.000,
   21.500 , 204.000,
   22.500 , 204.000,
   23.500 , 237.000,
   24.500 , 234.000,
   25.500 , 236.000,
   26.500 , 239.000,
   27.500 , 241.000,
   28.500 , 222.000,
   29.500 , 247.000,
   30.500 , 249.000,
   31.500 , 217.000,
   32.500 , 242.000,
   33.500 , 232.000,
   34.500 , 216.000,
   35.500 , 201.000,
   36.500 , 225.000,
   37.500 , 215.000,
   38.500 , 183.000,
   39.500 , 226.000,
   40.500 , 205.000,
   41.500 , 238.000,
   42.500 , 223.000,
   43.500 , 233.000,
   44.500 , 241.000,
   45.500 , 216.000,
   46.500 , 211.000,
   47.500 , 246.000,
   48.500 , 275.000,
   49.500 , 266.000,
   50.500 , 313.000,
   51.500 , 319.000,
   52.500 , 373.000,
   53.500 , 369.000,
   54.500 , 378.000,
   55.500 , 376.000,
   56.500 , 376.000,
   57.500 , 357.000,
   58.500 , 361.000,
   59.500 , 326.000,
   60.500 , 337.000,
   61.500 , 278.000,
   62.500 , 260.000,
   63.500 , 193.000,
   64.500 , 193.000,
   65.500 , 124.000,
   66.500 , 113.000,
   67.500 , 102.000,
   68.500 , 82.0000,
   69.500 , 45.0000,
   70.500 , 34.0000,
   71.500 , 32.0000,
   72.500 , 16.0000,
   73.500 , 12.0000,
   74.500 , 8.00000,
   75.500 , 3.00000,
   76.500 , 2.00000,
   77.500 , 0.00000,
   78.500 , 0.00000,
   79.500 , 1.00000,
   80.500 , 0.00000,
   81.500 , 0.00000,
   82.500 , 0.00000,
   83.500 , 0.00000,
   84.500 , 0.00000,
   85.500 , 0.00000,
   86.500 , 0.00000,
   87.500 , 0.00000,
   88.500 , 0.00000,
   89.500 , 0.00000,
   90.500 , 0.00000,
   91.500 , 0.00000,
   92.500 , 0.00000,
   93.500 , 0.00000,
   94.500 , 0.00000,
   95.500 , 0.00000,
   96.500 , 0.00000,
   97.500 , 0.00000,
   98.500 , 0.00000,
   99.500 , 0.00000,
  };

  return dd[i][k];
}

float dat6(int i, int k) {

  float dd[100][2] = {   
   1.0000 , 10.0000, //# of data pairs 100
   3.0000 , 30.0000,
   5.0000 , 54.0000,
   7.0000 , 70.0000,
   9.0000 , 112.000,
   11.000 , 159.000,
   13.000 , 188.000,
   15.000 , 252.000,
   17.000 , 298.000,
   19.000 , 360.000,
   21.000 , 341.000,
   23.000 , 434.000,
   25.000 , 454.000,
   27.000 , 507.000,
   29.000 , 484.000,
   31.000 , 500.000,
   33.000 , 451.000,
   35.000 , 450.000,
   37.000 , 462.000,
   39.000 , 466.000,
   41.000 , 453.000,
   43.000 , 442.000,
   45.000 , 451.000,
   47.000 , 525.000,
   49.000 , 575.000,
   51.000 , 646.000,
   53.000 , 733.000,
   55.000 , 829.000,
   57.000 , 840.000,
   59.000 , 769.000,
   61.000 , 645.000,
   63.000 , 468.000,
   65.000 , 332.000,
   67.000 , 204.000,
   69.000 , 126.000,
   71.000 , 65.0000,
   73.000 , 24.0000,
   75.000 , 5.00000,
   77.000 , 2.00000,
   79.000 , 0.00000,
   81.000,  0.00000,
   83.000,  0.00000,
   85.000,  0.00000,
   87.000,  0.00000,
   89.000,  0.00000,
   91.000,  0.00000,
   93.000,  0.00000,
   95.000,  0.00000,
   97.000,  0.00000,
   99.000,  0.00000,
   101.00,  0.00000,
   103.00,  0.00000,
   105.00,  0.00000,
   107.00,  0.00000,
   109.00,  0.00000,
   111.00,  0.00000,
   113.00,  0.00000,
   115.00,  0.00000,
   117.00,  0.00000,
   119.00,  0.00000,
   121.00,  0.00000,
   123.00,  0.00000,
   125.00,  0.00000,
   127.00,  0.00000,
   129.00,  0.00000,
   131.00,  0.00000,
   133.00,  0.00000,
   135.00,  0.00000,
   137.00,  0.00000,
   139.00,  0.00000,
   141.00,  0.00000,
   143.00,  0.00000,
   145.00,  0.00000,
   147.00,  0.00000,
   149.00,  0.00000,
   151.00,  0.00000,
   153.00,  0.00000,
   155.00,  0.00000,
   157.00,  0.00000,
   159.00,  0.00000,
   161.00,  0.00000,
   163.00,  0.00000,
   165.00,  0.00000,
   167.00,  0.00000,
   169.00,  0.00000,
   171.00,  0.00000,
   173.00,  0.00000,
   175.00,  0.00000,
   177.00,  0.00000,
   179.00,  0.00000,
   181.00,  0.00000,
   183.00,  0.00000,
   185.00,  0.00000,
   187.00,  0.00000,
   189.00,  0.00000,
   191.00,  0.00000,
   193.00,  0.00000,
   195.00,  0.00000,
   197.00,  0.00000,
   199.00,  0.00000, 
  };

  return dd[i][k];
}
//end of third graph data

//Emip

float dat7(int i, int k) {

  float dd[100][2] = {   
   0.5000,  0.00000, //# of data pairs 100
   1.5000,  0.00000,
   2.5000,  0.00000,
   3.5000,  0.00000,
   4.5000,  0.00000,
   5.5000,  5.00000,
   6.5000,  9.00000,
   7.5000,  2.00000,
   8.5000,  7.00000,
   9.5000,  6.00000,
   10.500,  7.00000,
   11.500,  3.00000,
   12.500,  1.10000,
   13.500,  3.00000,
   14.500,  5.00000,
   15.500,  3.00000,
   16.500,  8.00000,
   17.500,  7.00000,
   18.500,  6.00000,
   19.500,  4.00000,
   20.500,  4.00000,
   21.500,  12.0000,
   22.500,  4.00000,
   23.500,  6.00000,
   24.500,  8.00000,
   25.500,  10.0000,
   26.500,  12.0000,
   27.500,  7.00000,
   28.500,  11.0000,
   29.500,  12.0000,
   30.500,  15.0000,
   31.500,  13.0000,
   32.500,  17.0000,
   33.500,  20.0000,
   34.500,  23.0000,
   35.500,  20.0000,
   36.500,  29.0000,
   37.500,  24.0000,
   38.500,  26.0000,
   39.500,  31.0000,
   40.500,  33.0000,
   41.500,  30.0000,
   42.500,  47.0000,
   43.500,  55.0000,
   44.500,  76.0000,
   45.500,  65.0000,
   46.500,  75.0000,
   47.500,  94.0000,
   48.500,  112.000,
   49.500,  133.000,
   50.500,  178.000,
   51.500,  182.000,
   52.500,  236.000,
   53.500,  266.000,
   54.500,  274.000,
   55.500,  280.000,
   56.500,  293.000,
   57.500,  290.000,
   58.500,  295.000,
   59.500,  263.000,
   60.500,  291.000,
   61.500,  239.000,
   62.500,  224.000,
   63.500,  170.000,
   64.500,  174.000,
   65.500,  115.000,
   66.500,  95.0000,
   67.500,  95.0000,
   68.500,  80.0000,
   69.500,  41.0000,
   70.500,  33.0000,
   71.500,  28.0000,
   72.500,  16.0000,
   73.500,  12.0000,
   74.500,  8.00000,
   75.500,  3.00000,
   76.500,  2.00000,
   77.500,  0.00000,
   78.500,  0.00000,
   79.500,  0.00000,
   80.500,  0.00000,
   81.500,  0.00000,
   82.500,  0.00000,
   83.500,  0.00000,
   84.500,  0.00000,
   85.500,  0.00000,
   86.500,  0.00000,
   87.500,  0.00000,
   88.500,  0.00000,
   89.500,  0.00000,
   90.500,  0.00000,
   91.500,  0.00000,
   92.500,  0.00000,
   93.500,  0.00000,
   94.500,  0.00000,
   95.500,  0.00000,
   96.500,  0.00000,
   97.500,  0.00000,
   98.500,  0.00000,
   99.500,  0.00000,
  };

  return dd[i][k];
}

float dat8(int i, int k) {

  float dd[100][2] = {   
   1.0000,  0.00000, //# of data pairs 100
   3.0000,  0.00000,
   5.0000,  12.0000,
   7.0000,  8.00000,
   9.0000,  8.00000,
   11.000,  14.0000,
   13.000,  8.00000,
   15.000,  14.0000,
   17.000,  15.0000,
   19.000,  15.0000,
   21.000,  11.0000,
   23.000,  11.0000,
   25.000,  7.00000,
   27.000,  17.0000,
   29.000,  21.0000,
   31.000,  23.0000,
   33.000,  34.0000,
   35.000,  33.0000,
   37.000,  50.0000,
   39.000,  48.0000,
   41.000,  78.0000,
   43.000,  86.0000,
   45.000,  139.000,
   47.000,  217.000,
   49.000,  273.000,
   51.000,  386.000,
   53.000,  463.000,
   55.000,  589.000,
   57.000,  647.000,
   59.000,  614.000,
   61.000,  519.000,
   63.000,  394.000,
   65.000,  274.000,
   67.000,  182.000,
   69.000,  109.000,
   71.000,  57.0000,
   73.000,  21.0000,
   75.000,  5.00000,
   77.000,  2.00000,
   79.000,  0.00000,
   81.000,  0.00000,
   83.000,  0.00000,
   85.000,  0.00000,
   87.000,  0.00000,
   89.000,  0.00000,
   91.000,  0.00000,
   93.000,  0.00000,
   95.000,  0.00000,
   97.000,  0.00000,
   99.000,  0.00000,
   101.00,  0.00000,
   103.00,  0.00000,
   105.00,  0.00000,
   107.00,  0.00000,
   109.00,  0.00000,
   111.00,  0.00000,
   113.00,  0.00000,
   115.00,  0.00000,
   117.00,  0.00000,
   119.00,  0.00000,
   121.00,  0.00000,
   123.00,  0.00000,
   125.00,  0.00000,
   127.00,  0.00000,
   129.00,  0.00000,
   131.00,  0.00000,
   133.00,  0.00000,
   135.00,  0.00000,
   137.00,  0.00000,
   139.00,  0.00000,
   141.00,  0.00000,
   143.00,  0.00000,
   145.00,  0.00000,
   147.00,  0.00000,
   149.00,  0.00000,
   151.00,  0.00000,
   153.00,  0.00000,
   155.00,  0.00000,
   157.00,  0.00000,
   159.00,  0.00000,
   161.00,  0.00000,
   163.00,  0.00000,
   165.00,  0.00000,
   167.00,  0.00000,
   169.00,  0.00000,
   171.00,  0.00000,
   173.00,  0.00000,
   175.00,  0.00000,
   177.00,  0.00000,
   179.00,  0.00000,
   181.00,  0.00000,
   183.00,  0.00000,
   185.00,  0.00000,
   187.00,  0.00000,
   189.00,  0.00000,
   191.00,  0.00000,
   193.00,  0.00000,
   195.00,  0.00000,
   197.00,  0.00000,
   199.00,  0.00000,
  };

  return dd[i][k];
}


int tb_57_pi_2x2() {
   gROOT->Reset();
   float x[1000], y[1000], ex[1000], ey[1000];
   int n=100;

   c1 = new TCanvas(" tb_57_pi_2x2"," tb_57_pi_2x2",0,0,800,600);
   c1->GetFrame()->SetFillColor(21);
   c1->GetFrame()->SetBorderSize(12);


   // Inside the canvas, we create 4 pads

pad1 = new TPad("pad1", "",0.02,.52,0.48,0.98,21);
pad2 = new TPad("pad2", "",0.52,.52,0.98,0.98,21);
pad3 = new TPad("pad3", "",0.02,.02,0.48,0.48,21);
pad4 = new TPad("pad4", "",0.52,.02,0.98,0.48,21);
pad1->Draw();
pad2->Draw();
pad3->Draw();
pad4->Draw();


// Edit pad1

   pad1->cd();

   TH2F *hr1 = new TH2F("hr1","Etot:All Pions",
		       1,0.,100,1,0,2500);

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

// create graph

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


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

   TText* text = new TText(60,2200,"Blue: Test Beam");
   text->Draw();     

   for (int i=0; i<n; i++) {
     x [i] = dat2(i,0);
     y [i] = dat2(i,1); //note the renormalization
     ex[i] = 0;
     ey[i] = 0;
   }

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

   TText* text = new TText(60,2000,"Green: QFL");
   text->Draw();     

c1->Update();

   // Edit pad2

pad2->cd();

 TH2F *hr2 = new TH2F("hr2","EEM: All Pions",
		       1,0.,60,1,0,800);

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

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

   gr1 = new TGraphErrors(n,x,y,ex,ey); 
//Note that I revert to gr1 for the new pad.  Otherwise, there's a bug!
   gr1->SetMarkerColor(kBlue);
   gr1->SetMarkerStyle(21);
   gr1->Draw("LP");

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

     for (int i=0; i<n; i++) {
     x [i] = dat4(i,0);
     y [i] = dat4(i,1);//notice the renormalization
     ex[i] = 0;
     ey[i] = 0;
   }

   printf("%f",n);
   gr2 = new TGraphErrors(n,x,y,ex,ey);
   gr2->SetMarkerColor(kGreen);
   gr2->SetMarkerStyle(20);
   gr2->Draw("LP");

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

c1->Update();

//done with pad 2
   pad3->cd();

 TH2F *hr3 = new TH2F("hr3","Ehad: All Pions",
		       1,0.,100,1,0,800);

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

    for (int i=0; i<n; i++) {
     x [i] = dat5(i,0);
     y [i] = dat5(i,1);
     ex[i] = 0;
     ey[i] = 0;
   }
   printf("%f",n);
   gr1 = new TGraphErrors(n,x,y,ex,ey);
   gr1->SetMarkerColor(kBlue);
   gr1->SetMarkerStyle(21);
   gr1->Draw("LP");

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


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

   printf("%f",n);
   gr2 = new TGraphErrors(n,x,y,ex,ey);
   gr2->SetMarkerColor(kGreen);
   gr2->SetMarkerStyle(20);
   gr2->Draw("LP");

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


c1->Update();

   // pad4 commands

   pad4->cd();

 TH2F *hr4 = new TH2F("hr4","Emip",
		       1,0.,100,1,0,800);

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

   for (int i=0; i<n; i++) {
     x [i] = dat7(i,0);
     y [i] = dat7(i,1)*2;
     ex[i] = 0;
     ey[i] = 0;
   }
   printf("%f",n);
   gr1 = new TGraphErrors(n,x,y,ex,ey);
   gr1->SetMarkerColor(kBlue);
   gr1->SetMarkerStyle(22);
   gr1->Draw("LP");

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

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

   gr2 = new TGraphErrors(n,x,y,ex,ey);
   gr2->SetMarkerColor(kGreen);
   gr2->SetMarkerStyle(20);
   gr2->Draw("PL");

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

   //   label.Draw();
   c1->Update();

}