{

  TChain *a = new TChain("tree1");
  //  a->Add("2115_0.root");
  //  a->Add("2115_1.root");
  //  a->Add("2115_2.root");
  //  a->Add("2115_3.root");

  //  a->Add("2115_4.root");
  //  a->Add("2115_5.root");
  //  a->Add("2115_6.root");
  //  a->Add("2115_7.root");
  //  a->Add("2115_8.root");
  
  //  a->Add("2115_8_new.root");
  
  
  //  a->Add("hd_root.root");  
  //  a->Add("2115_00.root");
  
  a->Add("2115_0_new.root");
  a->Add("2115_1_new.root");
  a->Add("2115_2_new.root");
  a->Add("2115_3_new.root");
  a->Add("2115_4_new.root");
  a->Add("2115_5_new.root");
  a->Add("2115_6_new.root");
  a->Add("2115_7_new.root");
  a->Add("2115_8_new.root");
  a->Add("2115_9_new.root");

  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->SetLabelSize(0.073,"xy");
  gStyle->SetTitleSize(0.086,"xy");

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


  gStyle->SetOptStat(10);

  a->Draw("en>>h11","r_tr>60&&prob_tr>0.2&&mom_tr>1.0&&ndf>12&&en/mom_tr<0.4&&en>-10&&ntr==1","zcol");
  a->Draw("en1>>h22","r_tr>60&&prob_tr>0.2&&mom_tr>1.05&&ndf>12&&en/mom_tr<0.4&&en1>0&&ntr==1","zcol");
  a->Draw("en2>>h33","r_tr>60&&prob_tr>0.2&&mom_tr>1.0&&ndf>12&&en/mom_tr<0.4&&en2>0&&ntr==1","zcol");
  
  a->Draw("r_tr>>h100");
  a->Draw("en/mom_tr:mom_tr>>h101","en/mom_tr<2&&en/mom_tr>0","zcol");

  TCanvas *c1 = new TCanvas("c1","FCAL energy", 200, 10, 800, 300);
  
  c1->Divide(3,1);
  
  c1->cd(1);

  h11->SetMaximum(30);
  h11->SetLineColor(1);
  h11->Draw();
  h11->SetTitle("");
  h11->GetXaxis()->SetTitle("E  (GeV)");
  h11->GetYaxis()->SetTitle("Events");
  
  TLatex *text1 = new TLatex(1.2, 4, "E_{cell} = 10 MeV");
  text1 -> SetTextFont(132);
  text1->SetTextSize(0.06);
  text1->Draw();


  h22->SetMaximum(30);
  c1->cd(2);
  h22->SetLineColor(1);
  h22->Draw();
  h22->SetTitle("");
  h22->GetXaxis()->SetTitle("E  (GeV)");

  TLatex *text2 = new TLatex(1.2, 4, "E_{cell} = 90 MeV");
  text2 -> SetTextFont(132);
  text2->SetTextSize(0.06);
  text2->Draw(); 

  h33->SetMaximum(30);
  c1->cd(3);
  h33->SetLineColor(1);
  h33->Draw();
  h33->SetTitle("");
  h33->GetXaxis()->SetTitle("E  (GeV)");

  TLatex *text3 = new TLatex(1.2, 4.0, "E_{cell} = 90 MeV");
  text3 -> SetTextFont(132);
  text3->SetTextSize(0.06);
  text3->Draw(); 

  TCanvas *c2 = new TCanvas("c2","SiPM shape", 200, 10, 400, 300);
  h100->SetTitle("");
  h100->Draw("zcol");


}