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

  gROOT->ForceStyle();

  gStyle->SetOptStat(10);
  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.045,"xy");
  //  gStyle->SetTitleSize(0.05,"xy");
  //  gStyle->SetTitleOffset(1.1,"xy");

  //  gStyle->SetLabelSize(0.059,"xy");
  //  gStyle->SetTitleSize(0.065,"xy");
  //  gStyle->SetTitleOffset(1.,"xy");
  //  gStyle->SetTitleOffset(1,,"y");


  gStyle->SetOptLogz(1);

  //  TFile f("2115_0.root");
  TFile f("ref.root");

  TH2F *e_over_pp1 = (TH2F*)f.Get("e_over_pp1");
  TH2F *e_over_pp = (TH2F*)f.Get("e_over_pp");

  
  TCanvas *c1 = new TCanvas("c1","SiPM shape", 200, 10, 800, 400);


  c1->Divide(2,1);


  c1->cd(1);
  e_over_pp1->Draw("zcol");
  e_over_pp1->GetXaxis()->SetTitle("Track momentum (GeV)");
  e_over_pp1->GetYaxis()->SetTitle("E (cluster) / p");

  c1->cd(2);
  e_over_pp->Draw("zcol");
  
  e_over_pp->GetXaxis()->SetTitle("Track momentum (GeV)");
  e_over_pp->GetYaxis()->SetTitle("E (hit) / p");


  c1->SaveAs("e_over_p.jpg");
  



}