void plot_lg(){

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

  gROOT->ForceStyle();

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

  gStyle->SetTitleX(0.1f);
  gStyle->SetTitleW(0.8f);
  gStyle->SetTextSize(1.5);
  
  gStyle->SetTitleFont(132,"xy");
  gStyle->SetLabelFont(132,"xy");


  gStyle->SetLabelSize(0.055,"xy");
  gStyle->SetTitleSize(0.065,"xy");

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


  TFile *file = new TFile("test_ref.root","r");

  //  TFile *file = new TFile("test_lg_silicon.root","r");


  // Original
  //  TFile *file = new TFile("test_lg1.root","r");

  TTree *tree1 = (TTree*)file->Get("tree1");


  TFile *file1 = new TFile("test_lg1.root","r");
  //  TFile *file1 = new TFile("test_lg_silicon.root","r");

  TTree *tree2 = (TTree*)file1->Get("tree1");


  const int nbinx = 150;
  const float x_min = 0;
  const float x_max = 150;
  const int nbiny = 250;
  const float y_min = 0;
  const float y_max = 2500;

  TH2F *h1 = new TH2F("h1","h1",nbinx,x_min,x_max,nbiny,y_min,y_max);
  TH2F *h2 = new TH2F("h2","h2",nbinx,x_min,x_max,nbiny,y_min,y_max);
  TH2F *h3 = new TH2F("h3","h3",nbinx,x_min,x_max,nbiny,y_min,y_max);
  TH2F *h4 = new TH2F("h4","h4",nbinx,x_min,x_max,nbiny,y_min,y_max);
  TH2F *h5 = new TH2F("h5","h5",nbinx,x_min,x_max,nbiny,y_min,y_max);
  TH2F *h6 = new TH2F("h6","h6",nbinx,x_min,x_max,nbiny,y_min,y_max);
  TH2F *h7 = new TH2F("h7","h7",nbinx,x_min,x_max,nbiny,y_min,y_max);
  TH2F *h8 = new TH2F("h8","h8",nbinx,x_min,x_max,nbiny,y_min,y_max);
  TH2F *h9 = new TH2F("h9","h9",nbinx,x_min,x_max,nbiny,y_min,y_max);
  TH2F *h10 = new TH2F("h10","h10",nbinx,x_min,x_max,nbiny,y_min,y_max);
  TH2F *h11 = new TH2F("h11","h11",nbinx,x_min,x_max,nbiny,y_min,y_max);
  TH2F *h12 = new TH2F("h12","h12",nbinx,x_min,x_max,nbiny,y_min,y_max);

  tree1->Draw("amp_1:tiler>>h1","amp_1>0","zcol");
  tree1->Draw("amp_2:tiler>>h2","amp_2>0","zcol");
  tree1->Draw("amp_3:tiler>>h3","amp_3>0","zcol");
  tree1->Draw("amp_4:tiler>>h4","amp_4>0","zcol");
  tree1->Draw("amp_5:tiler>>h5","amp_5>0","zcol");
  tree1->Draw("amp_6:tiler>>h6","amp_6>0","zcol");
  tree1->Draw("amp_7:tiler>>h7","amp_7>0","zcol");
  tree1->Draw("amp_8:tiler>>h8","amp_8>0","zcol");
  tree1->Draw("amp_9:tiler>>h9","amp_9>0","zcol");
  tree1->Draw("amp_10:tiler>>h10","amp_10>0","zcol");
  tree1->Draw("amp_11:tiler>>h11","amp_11>0","zcol");
  tree1->Draw("amp_12:tiler>>h12","amp_12>0","zcol");



#if 0

  TCanvas *c1 = new TCanvas("c1","c1", 400, 10, 1200, 800);


  c1->Divide(4,3);

  c1->cd(1);

  gPad->SetLeftMargin(0.19);
  //  gPad->SetBottomMargin(0.19);
  h1->GetYaxis()->SetTitle("FADC peak (counts)");

  h1->SetTitle("Module 1");
  h1->Draw("zcol");

  c1->cd(2);

  h2->SetTitle("Module 2");
  h2->Draw("zcol");


  c1->cd(3);

  h3->SetTitle("Module 3");
  h3->Draw("zcol");

  c1->cd(4);

  h4->SetTitle("Module 4");
  h4->Draw("zcol");


  c1->cd(5);

  h5->GetYaxis()->SetTitle("FADC peak (counts)");
  gPad->SetLeftMargin(0.19);
  //  gPad->SetBottomMargin(0.19);

  h5->SetTitle("Module 5");
  h5->Draw("zcol");

  c1->cd(6);

  h6->SetTitle("Module 6");
  h6->Draw("zcol");


  c1->cd(7);

  h7->SetTitle("Module 7");
  h7->Draw("zcol");

  c1->cd(8);

  h8->SetTitle("Module 8");
  h8->Draw("zcol");



  c1->cd(9);

  gPad->SetLeftMargin(0.19);
  gPad->SetBottomMargin(0.19);

  h9->SetTitle("Module 9");
  h9->Draw("zcol");
  h9->GetXaxis()->SetTitle("PS tile");
  h9->GetYaxis()->SetTitle("FADC peak (counts)");

  c1->cd(10);
  gPad->SetBottomMargin(0.19);
  h10->SetTitle("Module 10");
  h10->Draw("zcol");
  h10->GetXaxis()->SetTitle("PS tile");

  c1->cd(11);
  gPad->SetBottomMargin(0.19);
  h11->SetTitle("Module 11");
  h11->Draw("zcol");
  h11->GetXaxis()->SetTitle("PS tile");

  c1->cd(12);
  gPad->SetBottomMargin(0.19);
  h12->SetTitle("Module 12");
  h12->Draw("zcol");

  h12->GetXaxis()->SetTitle("PS tile");


  //  gPad->SetBottomMargin(0.19);

  //  gr->SetMarkerStyle(20);
  //  gr->SetMarkerSize(0.6);
  //  gr->SetMarkerColor(4);

  //  gr->GetXaxis()->SetTitle("Tagger energy (GeV)");
  //  gr->GetYaxis()->SetTitle("#sigma(E)/E  (%)");

  //  gr->SetTitle("");

  //  gr->SetMinimum(0.2);
  //  gr->SetMaximum(10.3);


  //  gr->Draw("AP");
  
  //  gr->Fit(func1,"R","",3,10.5);

  TCanvas *c2 = new TCanvas("c2","c2", 400, 10, 800, 400);

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

  gStyle->SetLabelSize(0.053,"xy");
  gStyle->SetTitleSize(0.053,"xy");

  gStyle->SetTitleOffset(0.9,"x");
  gStyle->SetTitleOffset(0.9,"y");



  //  TH2F *h55 = new TH2F("h55","h55",150,0,150,125,0,2500);

  TProfile *h5555 = new TProfile("h5555","h5555",150,0,150,0,2500);

  h5555->GetYaxis()->SetTitle("FADC peak (counts)");
  h5555->GetXaxis()->SetTitle("PS tile ");

  tree1->Draw("amp_5:tiler>>h5555","amp_5>0","zcol");

  gPad->SetLeftMargin(0.19);
  gPad->SetBottomMargin(0.19);

  h5555->GetXaxis()->SetRangeUser(60,120);
  
  h5555->SetTitle("PWO Module 5");
  h5555->Draw("zcol");
  
#endif


  TCanvas *c3 = new TCanvas("c3","c3", 400, 10, 1200, 305);


  c3->Divide(3,1);


  c3->cd(1);
  

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

  gStyle->SetLabelSize(0.053,"xy");
  gStyle->SetTitleSize(0.053,"xy");

  gStyle->SetTitleOffset(0.9,"x");
  gStyle->SetTitleOffset(1.6,"y");

  gPad->SetLeftMargin(0.17);
  gPad->SetBottomMargin(0.19);

  TH1F *h9999 = new TH1F("h9999","h9999",150,-0.5,149.5);

  TProfile *h33 = new TProfile("h33","h33",150,-0.5,149.5,0,2500);
  TProfile *h55 = new TProfile("h55","h55",150,-0.5,149.5,0,2500);
  TProfile *h99 = new TProfile("h99","h99",150,-0.5,149.5,0,2500);
  
  TProfile *h333 = new TProfile("h333","h333",150,-0.5,149.5,0,2500);
  TProfile *h555 = new TProfile("h555","h555",150,-0.5,149.5,0,2500);
  TProfile *h999 = new TProfile("h999","h999",150,-0.5,149.5,0,2500);

  tree1->Draw("amp_3:tiler>>h33","amp_3>0","zcol");
  tree1->Draw("amp_5:tiler>>h55","amp_5>0","zcol");
  tree1->Draw("amp_9:tiler>>h99","amp_9>0","zcol");

  tree2->Draw("amp_3:tiler>>h333","amp_3>0","zcol");
  tree2->Draw("amp_5:tiler>>h555","amp_5>0","zcol");
  tree2->Draw("amp_9:tiler>>h999","amp_9>0","zcol");

  h33->GetYaxis()->SetTitle("FADC peak (counts)");
  h33->GetXaxis()->SetTitle("PS tile ");

  h33->GetXaxis()->SetRangeUser(100,120);

  h33->SetTitle("Module 3 (with LG)");

  h33->SetMarkerSize(0.7);
  h33->SetMarkerStyle(20);
  h33->SetMarkerColor(4);
  h33->SetFillColor(4);

  h333->SetMarkerSize(0.7);
  h333->SetMarkerStyle(21);
  h333->SetMarkerColor(2);
  h333->SetFillColor(4);

  h33->Draw("");


  h333->Draw("same");


  c3->cd(2);

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

  gStyle->SetLabelSize(0.053,"xy");
  gStyle->SetTitleSize(0.053,"xy");

  gStyle->SetTitleOffset(0.9,"x");
  gStyle->SetTitleOffset(0.9,"y");


  gPad->SetLeftMargin(0.1);
  gPad->SetBottomMargin(0.19);

  //  h55->GetYaxis()->SetTitle("FADC peak (counts)");
  h55->GetXaxis()->SetTitle("PS tile ");

  h55->GetXaxis()->SetRangeUser(80,100);

  h55->SetTitle("Module 5 (no LG)");


  h55->SetMarkerSize(0.7);
  h55->SetMarkerStyle(20);
  h55->SetMarkerColor(4);
  h55->SetFillColor(4);

  h555->SetMarkerSize(0.7);
  h555->SetMarkerStyle(21);
  h555->SetMarkerColor(2);
  h555->SetFillColor(4);



  h55->Draw("");
  h555->Draw("same");

  c3->cd(3);

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

  gStyle->SetLabelSize(0.053,"xy");
  gStyle->SetTitleSize(0.053,"xy");

  gStyle->SetTitleOffset(0.9,"x");
  gStyle->SetTitleOffset(0.9,"y");


  gPad->SetLeftMargin(0.1);
  gPad->SetBottomMargin(0.19);

  //  h99->GetYaxis()->SetTitle("FADC peak (counts)");
  h99->GetXaxis()->SetTitle("PS tile ");

  h99->GetXaxis()->SetRangeUser(50,70);

  h99->SetTitle("Module 9 (with LG)");

  h99->SetMarkerSize(0.7);
  h99->SetMarkerStyle(20);
  h99->SetMarkerColor(4);
  h99->SetFillColor(4);

  h999->SetMarkerSize(0.7);
  h999->SetMarkerStyle(21);
  h999->SetMarkerColor(2);
  h999->SetFillColor(4);

  h99->Draw("");
  h999->Draw("same");


  TCanvas *c4 = new TCanvas("c4","c4", 400, 10, 600, 400);

  gPad->SetBottomMargin(0.15);
  gPad->SetLeftMargin(0.12);
  gPad->SetRightMargin(0.04);

;



  for(int ii = 1; ii < 150; ii++){
    double tmp    = h99->GetBinContent(ii);
    double tmp_er = h99->GetBinError(ii);
    h9999->SetBinContent(ii,tmp);
    h9999->SetBinError(ii,tmp_er);  
    if(ii == 63) h9999->SetBinContent(ii,1630);
  }

  h9999->GetXaxis()->SetTitle("Pair spectrometer tile ");
  h9999->GetYaxis()->SetTitle("ADC  count");

  h9999->SetMarkerSize(0.9);
  h9999->SetMarkerStyle(20);
  h9999->SetMarkerColor(4);
  h9999->SetFillColor(4);

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

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

  h9999->SetLabelSize(0.055,"xy");
  h9999->SetTitleSize(0.065,"xy");

  h9999->SetTitleOffset(1.0,"y");
  h9999->SetTitleOffset(1.0,"x");

  h9999->SetTitle("");
  h9999->GetXaxis()->SetRangeUser(50,70);

  h9999->SetMaximum(2000);
  h9999->Draw("");

 

  h999->SetMarkerSize(0.9);
  h999->Scale(1.1);
  h999->Draw("same");

  TLegend *leg = new TLegend(0.15,0.65,0.4,0.85);
  leg->SetFillColor(0);
  leg->SetTextFont(132);
  leg->SetTextSize(0.06);

  leg->SetBorderSize(0);

  leg->AddEntry(h9999, " No light guide ","p");
  leg->AddEntry(h999," With light guide ","p");

  leg->Draw();


  

}