void plot(){

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

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

  gStyle->SetTitleOffset(1.2,"x");
  gStyle->SetTitleOffset(1.4,"y");


  gStyle->SetOptStat(0);


  gStyle->SetOptLogy(0);
  gStyle->SetOptLogz(1);

  gStyle->SetPalette(1);


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

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

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

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

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

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

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

  TFile *file = new TFile("9_2_crytur.root","r");
  
  TTree *tree1 = (TTree*)file->Get("tree1");
  
  const int nbinx = 150;
  const float x_min = -0.5;
  const float x_max = 149.5;
  //  const int nbiny = 250;
  const float y_min = 0;
  //  const float y_max = 2500;


  const int nbiny = 320;
  const int y_max = 3200;


  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");



  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);


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

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

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

  h55->GetXaxis()->SetRangeUser(60,120);
  
  h55->SetTitle("PWO Module 5");
  h55->Draw("zcol");
  
}