//transcribed by Rudolph Magyar from tb_227_pi.top on 7/16/99 // float dat1(int i, int k) { float dd[78][3] = { 5.1075001, 10.7100010, 3.5700002, 8.5125008, 3.5700002, 2.0611405, 11.9174995, 8.3300009, 3.1484442, 15.3225002, 7.1400003, 2.9148932, 18.7275009, 1.1900001, 1.1900001, 25.5374985, 2.3800001, 1.6829141, 28.9424992, 2.3800001, 1.6829141, 32.3474998, 1.1900001, 1.1900001, 35.7524986, 2.3800001, 1.6829141, 39.1575012, 2.3800001, 1.6829141, 42.5625000, 3.5700002, 2.0611405, 45.9674988, 2.3800001, 1.6829141, 52.7775002, 1.1900001, 1.1900001, 56.1824989, 1.1900001, 1.1900001, 59.5874977, 1.1900001, 1.1900001, 62.9925003, 1.1900001, 1.1900001, 69.8024979, 4.7600002, 2.3800001, 73.2074966, 2.3800001, 1.6829141, 80.0175018, 4.7600002, 2.3800001, 83.4225006, 1.1900001, 1.1900001, 86.8274994, 1.1900001, 1.1900001, 90.2324982, 4.7600002, 2.3800001, 93.6374969, 1.1900001, 1.1900001, 97.0425034, 2.3800001, 1.6829141, 100.4475021, 4.7600002, 2.3800001, 103.8525009, 1.1900001, 1.1900001, 107.2574997, 3.5700002, 2.0611405, 110.6625061, 4.7600002, 2.3800001, 114.0674973, 4.7600002, 2.3800001, 117.4724960, 3.5700002, 2.0611405, 120.8775101, 2.3800001, 1.6829141, 124.2825012, 3.5700002, 2.0611405, 127.6875000, 8.3300009, 3.1484442, 131.0924988, 11.9000006, 3.7631106, 134.4974976, 3.5700002, 2.0611405, 137.9024963, 8.3300009, 3.1484442, 141.3074951, 11.9000006, 3.7631106, 144.7124939, 19.0400009, 4.7600002, 148.1174927, 4.7600002, 2.3800001, 151.5225067, 17.8500004, 4.6088505, 154.9275055, 9.5200005, 3.3658283, 158.3325043, 17.8500004, 4.6088505, 161.7374878, 11.9000006, 3.7631106, 165.1425018, 10.7100010, 3.5700002, 168.5475006, 26.1800022, 5.5815949, 171.9524994, 27.3700008, 5.7070398, 175.3574982, 27.3700008, 5.7070398, 178.7625122, 45.2200012, 7.3356528, 182.1674957, 58.3100014, 8.3300009, 185.5724945, 61.8800049, 8.5812120, 188.9774933, 110.6700058, 11.4759455, 192.3825073, 132.0900116, 12.5374289, 195.7875061, 189.2100067, 15.0053291, 199.1925049, 229.6700134, 16.5320091, 202.5975037, 248.7100067, 17.2036304, 206.0025024, 295.1200256, 18.7401390, 209.4075012, 336.7700195, 20.0188980, 212.8125000, 311.7800293, 19.2618351, 216.2174988, 304.6400146, 19.0400009, 219.6224976, 260.6100159, 17.6103935, 223.0274963, 232.0500183, 16.6174469, 226.4324951, 239.1900177, 16.8711624, 229.8374786, 193.9700165, 15.1929035, 233.2424774, 140.4200134, 12.9267092, 236.6474762, 117.8100052, 11.8403502, 240.0524750, 88.0600052, 10.2367678, 243.4574738, 65.4500046, 8.8252764, 246.8624878, 59.5000038, 8.4145708, 250.2674866, 34.5100021, 6.4083462, 253.6724854, 32.1300011, 6.1834216, 257.0774841, 20.2300014, 4.9064960, 260.4824829, 15.4700012, 4.2906060, 263.8874817, 15.4700012, 4.2906060, 267.2924805, 4.7600002, 2.3800001, 270.6974792, 9.5200005, 3.3658283, 274.1024780, 1.1900001, 1.1900001, 277.5075073, 4.7600002, 2.3800001, 284.3174744, 10.7100010, 3.5700002, }; return dd[i][k]; } float dat2(int i, int k) { float dd[300][2] = { 5.00000, 0.00000, 1.50000, 0.00000, 2.50000, 0.00000, 3.50000, 0.00000, 4.50000, 0.00000, 5.50000, 6.00000, 6.50000, 5.00000, 7.50000, 0.00000, 8.50000, 2.00000, 9.50000, 4.00000, 10.5000, 3.00000, 11.5000, 3.00000, 12.5000, 5.00000, 13.5000, 1.00000, 14.5000, 3.00000, 15.5000, 4.00000, 16.5000, 2.00000, 17.5000, 0.00000, 18.5000, 2.00000, 19.5000, 0.00000, 20.5000, 2.00000, 21.5000, 0.00000, 22.5000, 1.00000, 23.5000, 2.00000, 24.5000, 3.00000, 25.5000, 1.00000, 26.5000, 1.00000, 27.5000, 0.00000, 28.5000, 1.00000, 29.5000, 2.00000, 30.5000, 1.00000, 31.5000, 0.00000, 32.5000, 0.00000, 33.5000, 2.00000, 34.5000, 1.00000, 35.5000, 0.00000, 36.5000, 0.00000, 37.5000, 3.00000, 38.5000, 2.00000, 39.5000, 3.00000, 40.5000, 1.00000, 41.5000, 1.00000, 42.5000, 1.00000, 43.5000, 2.00000, 44.5000, 0.00000, 45.5000, 1.00000, 46.5000, 1.00000, 47.5000, 5.00000, 48.5000, 1.00000, 49.5000, 1.00000, 50.5000, 3.00000, 51.5000, 1.00000, 52.5000, 1.00000, 53.5000, 0.00000, 54.5000, 4.00000, 55.5000, 3.00000, 56.5000, 1.00000, 57.5000, 2.00000, 58.5000, 3.00000, 59.5000, 3.00000, 60.5000, 0.00000, 61.5000, 2.00000, 62.5000, 1.00000, 63.5000, 4.00000, 64.5000, 1.00000, 65.5000, 0.00000, 66.5000, 3.00000, 67.5000, 1.00000, 68.5000, 1.00000, 69.5000, 3.00000, 70.5000, 5.00000, 71.5000, 2.00000, 72.5000, 0.00000, 73.5000, 1.00000, 74.5000, 2.00000, 75.5000, 0.00000, 76.5000, 2.00000, 77.5000, 4.00000, 78.5000, 2.00000, 79.5000, 4.00000, 80.5000, 1.00000, 81.5000, 4.00000, 82.5000, 4.00000, 83.5000, 1.00000, 84.5000, 3.00000, 85.5000, 2.00000, 86.5000, 3.00000, 87.5000, 3.00000, 88.5000, 3.00000, 89.5000, 4.00000, 90.5000, 1.00000, 91.5000, 5.00000, 92.5000, 2.00000, 93.5000, 4.00000, 94.5000, 2.00000, 95.5000, 4.00000, 96.5000, 4.00000, 97.5000, 9.00000, 98.5000, 2.00000, 99.5000, 5.00000, 100.500, 2.00000, 101.500, 2.00000, 102.500, 2.00000, 103.500, 5.00000, 104.500, 7.00000, 105.500, 4.00000, 106.500, 2.00000, 107.500, 5.00000, 108.500, 4.00000, 109.500, 8.00000, 110.500, 0.00000, 111.500, 6.00000, 112.500, 7.00000, 113.500, 3.00000, 114.500, 4.00000, 115.500, 2.00000, 116.500, 3.00000, 117.500, 8.00000, 118.500, 3.00000, 119.500, 6.00000, 120.500, 7.00000, 121.500, 7.00000, 122.500, 9.00000, 123.500, 5.00000, 124.500, 4.00000, 125.500, 5.00000, 126.500, 1.00000, 127.500, 5.00000, 128.500, 8.00000, 129.500, 5.00000, 130.500, 4.00000, 131.500, 5.00000, 132.500, 11.0000, 133.500, 8.00000, 134.500, 5.00000, 135.500, 7.00000, 136.500, 8.00000, 137.500, 12.0000, 138.500, 11.0000, 139.500, 8.00000, 140.500, 12.0000, 141.500, 7.00000, 142.500, 9.00000, 143.500, 11.0000, 144.500, 6.00000, 145.500, 13.0000, 146.500, 4.00000, 147.500, 15.0000, 148.500, 7.00000, 149.500, 16.0000, 150.500, 15.0000, 151.500, 5.00000, 152.500, 20.0000, 153.500, 11.0000, 154.500, 9.00000, 155.500, 23.0000, 156.500, 15.0000, 157.500, 11.0000, 158.500, 19.0000, 159.500, 11.0000, 160.500, 15.0000, 161.500, 14.0000, 162.500, 19.0000, 163.500, 17.0000, 164.500, 23.0000, 165.500, 15.0000, 166.500, 24.0000, 167.500, 26.0000, 168.500, 24.0000, 169.500, 19.0000, 170.500, 23.0000, 171.500, 20.0000, 172.500, 28.0000, 173.500, 41.0000, 174.500, 35.0000, 175.500, 34.0000, 176.500, 25.0000, 177.500, 45.0000, 178.500, 45.0000, 179.500, 43.0000, 180.500, 51.0000, 181.500, 48.0000, 182.500, 67.0000, 183.500, 65.0000, 184.500, 78.0000, 185.500, 91.0000, 186.500, 76.0000, 187.500, 104.000, 188.500, 95.0000, 189.500, 89.0000, 190.500, 121.000, 191.500, 115.000, 192.500, 122.000, 193.500, 148.000, 194.500, 152.000, 195.500, 174.000, 196.500, 176.000, 197.500, 201.000, 198.500, 197.000, 199.500, 214.000, 200.500, 208.000, 201.500, 235.000, 202.500, 221.000, 203.500, 243.000, 204.500, 259.000, 205.500, 280.000, 206.500, 265.000, 207.500, 291.000, 208.500, 303.000, 209.500, 279.000, 210.500, 304.000, 211.500, 319.000, 212.500, 276.000, 213.500, 274.000, 214.500, 259.000, 215.500, 291.000, 216.500, 282.000, 217.500, 267.000, 218.500, 266.000, 219.500, 279.000, 220.500, 264.000, 221.500, 257.000, 222.500, 218.000, 223.500, 229.000, 224.500, 222.000, 225.500, 233.000, 226.500, 201.000, 227.500, 218.000, 228.500, 176.000, 229.500, 190.000, 230.500, 180.000, 231.500, 161.000, 232.500, 149.000, 233.500, 157.000, 234.500, 154.000, 235.500, 140.000, 236.500, 158.000, 237.500, 133.000, 238.500, 108.000, 239.500, 130.000, 240.500, 104.000, 241.500, 110.000, 242.500, 102.000, 243.500, 106.000, 244.500, 81.0000, 245.500, 83.0000, 246.500, 82.0000, 247.500, 67.0000, 248.500, 73.0000, 249.500, 52.0000, 250.500, 71.0000, 251.500, 44.0000, 252.500, 49.0000, 253.500, 50.0000, 254.500, 43.0000, 255.500, 42.0000, 256.500, 48.0000, 257.500, 40.0000, 258.500, 49.0000, 259.500, 29.0000, 260.500, 23.0000, 261.500, 19.0000, 262.500, 26.0000, 263.500, 23.0000, 264.500, 24.0000, 265.500, 17.0000, 266.500, 17.0000, 267.500, 11.0000, 268.500, 13.0000, 269.500, 10.0000, 270.500, 9.00000, 271.500, 9.00000, 272.500, 12.0000, 273.500, 8.00000, 274.500, 4.00000, 275.500, 3.00000, 276.500, 6.00000, 277.500, 1.00000, 278.500, 4.00000, 279.500, 1.00000, 280.500, 0.00000, 281.500, 1.00000, 282.500, 1.00000, 283.500, 2.00000, 284.500, 1.00000, 285.500, 1.00000, 286.500, 0.00000, 287.500, 0.00000, 288.500, 0.00000, 289.500, 1.00000, 290.500, 0.00000, 291.500, 0.00000, 292.500, 1.00000, 293.500, 0.00000, 294.500, 0.00000, 295.500, 0.00000, 296.500, 0.00000, 297.500, 0.00000, 298.500, 0.00000, 299.500, 0.00000, }; return dd[i][k]; } int tb_227_pi() { gROOT->Reset(); c1 = new TCanvas("tb_227_pi","tb_227_pi",0,0,600,400); // create a 2D histogram to define the range TH2F *hr = new TH2F("hr","Etot: All Pions", 1,0.,300,1,0,380); 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 = 78; for (int i=0; iSetMarkerColor(kBlue); gr1->SetMarkerStyle(21); gr1->Draw("LP"); TText* text = new TText(40,60,"Blue: 227.0 GeV Testbeam"); text->Draw(); int n2 = 300; for (int i=0; iSetMarkerColor(kRed); gr2->SetMarkerStyle(22); gr2->Draw("LP"); TText* text = new TText(40,100,"Red: QFL"); text->Draw(); }