{
  gROOT->Reset();
  gROOT->SetStyle("Plain");

  gROOT->ForceStyle();

  gStyle->SetOptStat(0);
  gStyle->SetOptFit(0);

  gStyle->SetTitleX(0.1f);
  gStyle->SetTitleW(0.8f);

  gStyle->SetLineWidth(1.5);
  gStyle->SetTextSize(1.5);
  
  gStyle->SetTitleFont(132,"xy");
  gStyle->SetLabelFont(132,"xy");

  //  gStyle->SetLabelFont(62,"xy");

  //  gStyle->SetLabelSize(0.073,"xy");
  gStyle->SetLabelSize(0.05,"xy");
  gStyle->SetTitleSize(0.055,"xy");

  gStyle->SetTitleOffset(1.0,"x");
  gStyle->SetTitleOffset(1.1,"y");

  gStyle->SetTitleBorderSize(0);


  Double_t current[6] = {50, 100, 200, 300, 400, 500};

  Double_t ex[6];


  Double_t rate1[6]      =  {11.4, 23.6, 47.6, 73, 99.2, 127.1};
  Double_t rate1_ps[6]   =  {0.95, 2.72, 3.93, 6.03, 7.7, 9.62};
  
  Double_t rate2[6]     =  {8.5, 19.0, 32.9, 51.1, 67, 86.8};
  Double_t rate2_ps[6]  =  {1.0, 2.1, 3.69, 5.69, 7.23, 9.52};

  Double_t rate3[6]     =  {6.5, 13.5, 26.7, 40.2, 54, 69};
  Double_t rate3_ps[6]  =  {0.97, 1.93, 3.88, 5.84, 7.71, 9.51};

  //  for(int ii = 0; ii < 6; ii++){
  //    rate1[ii]  =  rate1[ii] - rate1_ps[ii];
  //    rate2[ii]  =  rate2[ii] - rate2_ps[ii];
  //    rate3[ii]  =  rate3[ii] - rate3_ps[ii];
  //  }

  TCanvas *c1 = new TCanvas("c1","c1",200,10, 600, 500);

  gPad->SetBottomMargin(0.19);
  gPad->SetTopMargin(0.1);
  gPad->SetLeftMargin(0.15);
  gPad->SetRightMargin(0.015);
  

  gPad->SetGrid(1);

  TGraph *gr1  = new TGraph(6, current, rate1); 
  gr1->SetMarkerStyle(20);
  gr1->SetMarkerSize(1.2);
  gr1->SetMarkerColor(2); 
  gr1->SetFillColor(4);

  gr1->SetTitle("BCAL & FCAL Trigger");
  gr1->GetXaxis()->SetTitle("Current (nA)");
  gr1->GetYaxis()->SetTitle("Rate  (kHz)");
  gr1->Draw("APL");
  

  TGraph *gr2  = new TGraph(6, current, rate2); 
  gr2->SetMarkerStyle(21);
  gr2->SetMarkerSize(1.2);
  gr2->SetMarkerColor(4); 
  gr2->SetFillColor(4);
  gr2->Draw("PLsame");


  TGraph *gr3  = new TGraph(6, current, rate3); 
  gr3->SetMarkerStyle(22);
  gr3->SetMarkerSize(1.2);
  gr3->SetMarkerColor(6); 
  gr3->SetFillColor(5);
  gr3->Draw("PLsame");

  leg = new TLegend(0.2,0.65,0.75,0.85);
  leg->SetFillColor(0);
  leg->SetTextSize(0.025);

  leg->AddEntry(gr1,"E^{FCAL}_{THR} = 0.5 GeV  (E_{FCAL} + 0.5 E_{BCAL} > 0.5 GeV ) ","p");
  leg->AddEntry(gr2,"E^{FCAL}_{THR} = 0.83 GeV ","p");
  leg->AddEntry(gr3,"E^{FCAL}_{THR} = 1.25 GeV ","p");
  leg->Draw();


}