{
  // ========================================================================
  // --- make plot of CDF expected Mtop uncertainties
  // --- created: Mar-2006 (D.Glenzinski)
  // --- based on extrapolations of the CDF Average Mt using 750/pb of data
  // --- cf. http://www-cdf.fnal.gov/internal/physics/top/run2mass/CDFExtrap/
  // --- for more details
  // ========================================================================
  gROOT->Reset();

  Int_t MAX = 10;
  Float_t Lumi[MAX];
  Float_t rootL[MAX];
  Float_t dampd[MAX];

  // --- luminosities
  Lumi[0] = 110;
  Lumi[1] = 320;
  Lumi[2] = 680;
  Lumi[3] = 1000;
  Lumi[4] = 2000;
  Lumi[5] = 4000;
  Lumi[6] = 6000;
  Lumi[7] = 8000;
  
  // --- rootL extrapolation
  rootL[0] = 6.6;
  rootL[1] = 3.8;
  rootL[2] = 2.7;
  rootL[3] = 2.19;
  rootL[4] = 1.55;
  rootL[5] = 1.09;
  rootL[6] = 0.89;
  rootL[7] = 0.77;

  // --- "damped" extrapolation (ie. stat+jes scales as rootL, rest fixed)
  dampd[0] = 6.6;
  dampd[1] = 3.8;
  dampd[2] = 2.7;
  dampd[3] = 2.47;
  dampd[4] = 2.07;
  dampd[5] = 1.69;
  dampd[6] = 1.57;
  dampd[7] = 1.49;


  // --- make TGraphs to plot
  grLum = new TGraph(8, Lumi, rootL);
  grLum->SetLineColor(4);
  grLum->SetLineWidth(3);
  grLum->SetLineStyle(2);

  grDmp = new TGraph(8, Lumi, dampd);
  grDmp->SetLineColor(4);
  grDmp->SetLineWidth(3);
  grLum->SetLineStyle(4);


  Float_t OLum[3], OSta[3], OTot[3];
  OLum[0] = 110;
  OSta[0] = 4.2;
  OTot[0] = 6.6;
  o110 = new TGraph(1, OLum, OTot);
  o110->SetMarkerStyle(23);
  o110->SetMarkerSize(1.5);
  o110->SetMarkerColor(4);
  o110->SetLineWidth(3);

  OLum[0] = 320;
  OSta[0] = 2.3;
  OTot[0] = 3.8;
  o320 = new TGraph(1, OLum, OTot);
  o320->SetMarkerStyle(23);
  o320->SetMarkerSize(1.5);
  o320->SetMarkerColor(4);
  o320->SetLineWidth(3);

  OLum[0] = 680;
  OSta[0] = 1.6;
  OTot[0] = 2.7;
  o680 = new TGraph(1, OLum, OTot);
  o680->SetMarkerStyle(23);
  o680->SetMarkerSize(1.5);
  o680->SetMarkerColor(4);
  o680->SetLineWidth(3);

  OLum[0] = 2000;
  OSta[0] = 3.0;
  OTot[0] = 3.0;
  tdr = new TGraph(1, OLum, OTot);
  tdr->SetMarkerStyle(29);
  tdr->SetMarkerSize(2.0);
  tdr->SetMarkerColor(2);
  tdr->SetLineWidth(3);

  OLum[0] = 115;
  OTot[0] = 0.8;
  star = new TGraph(1, OLum, OTot);
  star->SetMarkerStyle(29);
  star->SetMarkerSize(2.0);
  star->SetMarkerColor(2);
  star->SetLineWidth(3);

  OLum[0] = 115;
  OTot[0] = 1.4;
  triangle = new TGraph(1, OLum, OTot);
  triangle->SetMarkerStyle(23);
  triangle->SetMarkerSize(2.0);
  triangle->SetMarkerColor(4);
  triangle->SetLineWidth(3);

  // --- plot it
  gROOT->SetStyle("Plain");
  gPad->SetLogx(1);
  gPad->SetLogy(1);
  gStyle->SetOptStat(0);

  TH2F *nullTH2F = new TH2F("nullTH2F", "", 200, 70., 20000., 200, 0., 30.);
  nullTH2F->GetXaxis()->SetTitle("Integrated Luminosity (pb^{-1})");
  nullTH2F->GetXaxis()->SetTitleOffset(1.2);
  nullTH2F->GetXaxis()->CenterTitle();
  nullTH2F->GetYaxis()->SetTitle("#Delta M_{t}(total)  GeV/c^{2}");
  nullTH2F->GetYaxis()->SetTitleOffset(1.3);
  nullTH2F->GetYaxis()->CenterTitle();
  nullTH2F->Draw();

  grLum->Draw();
  grDmp->Draw();
  o110->Draw("P");
  o320->Draw("P");
  o680->Draw("P");
  tdr->Draw("P");
  star->Draw("P");
  triangle->Draw("P");

  // --- legend
  l = new TLine(1.90, -0.30, 2.18, -0.30);
  l->SetLineColor(4);
  l->SetLineWidth(3);
  l->Draw();

  l = new TLine(1.90, -0.60, 2.18, -0.60);
  l->SetLineColor(4);
  l->SetLineWidth(3);
  l->SetLineStyle(4);
  l->Draw();

  // --- line & arrow at DM/M <= 1%
  //l = new TLine(3.30, 0.237, 4.02, 0.237);
  //l->SetLineColor(2);
  //l->SetLineWidth(3);
  //l->SetLineStyle(4);
  //l->Draw();

  l = new TLine(3.25, 0.237, 3.54, 0.237);
  l->SetLineColor(2);
  l->SetLineWidth(3);
  l->Draw();

  //a = new TArrow(3.98, 0.237, 3.98, 0.146, 0.015, ">");
  a = new TArrow(3.45, 0.237, 3.45, 0.146, 0.015, ">");
  a->SetLineColor(2);
  a->SetLineWidth(3);
  a->Draw();

  // --- arrows to mark 1, 2, 4, 8 fb-1
  a = new TArrow(3.00, 0.70, 3.00, 0.60, 0.015, ">");
  a->SetLineColor(4);
  a->SetLineWidth(3);
  a->Draw();

  a = new TArrow(3.30, 0.70, 3.30, 0.60, 0.015, ">");
  a->SetLineColor(4);
  a->SetLineWidth(3);
  a->Draw();

  a = new TArrow(3.60, 0.70, 3.60, 0.60, 0.015, ">");
  a->SetLineColor(4);
  a->SetLineWidth(3);
  a->Draw();

  a = new TArrow(3.90, 0.70, 3.90, 0.60, 0.015, ">");
  a->SetLineColor(4);
  a->SetLineWidth(3);
  a->Draw();

  // --- title
  TLatex t;
  t->SetTextSize(0.05);
  t->DrawLatex(2.25, 1.23, "CDF Top Mass Uncertainty");
  t->SetTextSize(0.035);
  t->DrawLatex(2.39, 1.10, "(l+l and l+j channels combined)");

  // --- legend
  t->SetTextSize(0.035);
  t->DrawLatex(2.23,  0.11, "CDF Results");
  t->DrawLatex(2.23, -0.12, "Run IIa goal (TDR 1996)");
  t->DrawLatex(2.23, -0.35, "Scale #Delta(stat) / #sqrt{L},  Fix #Delta(syst)");
  t->DrawLatex(2.23, -0.47, "(assumes no improvements)");
  t->DrawLatex(2.23, -0.65, "Scale #Delta(total) / #sqrt{L}");
  t->DrawLatex(2.23, -0.77, "(improvements required)");

  // --- luminosities
  t->DrawLatex(2.93, 0.75, "1 fb^{-1}");
  t->DrawLatex(3.23, 0.75, "2 fb^{-1}");
  t->DrawLatex(3.53, 0.75, "4 fb^{-1}");
  t->DrawLatex(3.83, 0.75, "8 fb^{-1}");

  // --- DM/M < 1%
  t->SetTextColor(2);
  t->SetTextSize(0.032);
  t->SetTextAngle(-28);
  t->DrawLatex(3.50, 0.15, "#DeltaM/M < 1%");
  //t->DrawLatex(3.50, 0.15, "#DeltaM/M #approx #Gamma/M < 1%");


}