void pl_eff_XY(Int_t iMode=0)
{
  //  TCanvas *can = new TCanvas("can22","can22");
  //  can->Divide(2,2); 
  //  TCanvas *can = new TCanvas("can","can",48,55,700,900);
  TCanvas *can = new TCanvas("can","can",0,0,800,800);
  can->Divide(3,3,0.01,0.01); 
  //  can->Divide(2,2,0,0); 
  Float_t lsize=0.03;

 gROOT->cd();
 gPad->SetFillColor(0);
 gStyle->SetPalette(1);
 gStyle->SetLabelSize(lsize);

 gStyle->SetOptStat(kFALSE);
 //gROOT->cd();
 //gROOT->SetDirLevel(2);

 TH1 *h;
 TH2 *h2;

 TH2 *hs4;
 TH2 *hsS;
 TH2 *hs0;

 TH2 *heff4;
 TH2 *heffS;
 TH2 *heff0;

 Double_t Npredicted=0;
 Double_t Nfound=0.;
 TString hname="";

 can->cd(1);
   gROOT->cd();
   hname=Form("hXY4D4sel");
   hs4=(TH2 *)gROOT->FindObjectAny(hname);
   if (hs4!=NULL) {
     hs4->UseCurrentStyle(); hs4->GetYaxis()->SetLabelSize(lsize);
     hs4->Draw("colz");
   }else{cout<<"Histogram "<<hname<<" not existing. "<<endl;}

 can->cd(2);
   gROOT->cd();
   hname=Form("hXYSel2D4sel");
   hsS=(TH2 *)gROOT->FindObjectAny(hname);
   if (hsS!=NULL) {
     hsS->UseCurrentStyle(); hsS->GetYaxis()->SetLabelSize(lsize);
     hsS->Draw("colz");
   }else{cout<<"Histogram "<<hname<<" not existing. "<<endl;}

 can->cd(3);
   gROOT->cd();
   hname=Form("hXY0D4sel");
   hs0=(TH2 *)gROOT->FindObjectAny(hname);
   if (hs0!=NULL) {
     hs0->UseCurrentStyle(); hs0->GetYaxis()->SetLabelSize(lsize);
     hs0->Draw("colz");
     Npredicted=hs0->GetEntries();
   }else{cout<<"Histogram "<<hname<<" not existing. "<<endl;}

 can->cd(4);
   gROOT->cd();
   hname=Form("hXY4D4best");
   h2=(TH2 *)gROOT->FindObjectAny(hname);
   if (h2!=NULL) {
     h2->UseCurrentStyle(); h2->GetYaxis()->SetLabelSize(lsize);
     h2->Draw("colz");
     heff4=(TH2 *)h2->Clone();
     heff4->SetName("hXEffRef");
     heff4->SetTitle("Selection efficiency on reference counter");
     heff4->Divide(h2,hs4,1.,1.,"B");
   }else{cout<<"Histogram "<<hname<<" not existing. "<<endl;}

 can->cd(5);
   gROOT->cd();
   hname=Form("hXYSel2D4best");
   h2=(TH2 *)gROOT->FindObjectAny(hname);
   if (h2!=NULL) {
     h2->UseCurrentStyle(); h2->GetYaxis()->SetLabelSize(lsize);
     h2->Draw("colz");
     heffS=(TH2 *)h2->Clone();
     heffS->SetName("hXYEffSel2");
     heffS->SetTitle("Selection efficiency on Sel2 counter");
     heffS->Divide(h2,hsS,1.,1.,"B");   
   }else{cout<<"Histogram "<<hname<<" not existing. "<<endl;}

 can->cd(6);
   gROOT->cd();
   hname=Form("hXY0D4best");
   h2=(TH2 *)gROOT->FindObjectAny(hname);
   if (h2!=NULL) {
     h2->UseCurrentStyle(); h2->GetYaxis()->SetLabelSize(lsize);
     h2->Draw("colz");
     Nfound=h2->GetEntries();
     heff0=(TH2 *)h2->Clone();
     heff0->SetName("hXYEffDut");
     heff0->SetTitle("Efficiency of Dut counter");
     heff0->Divide(h2,hs0,1.,1.,"B"); 
     Double_t dEff=Nfound/Npredicted*100.;
     cout<<"Integral Efficiency: "<<Form("%3.1f",dEff)<<" % "<<endl;
 
   }else{cout<<"Histogram "<<hname<<" not existing. "<<endl;}

 can->cd(7);
   if (heff4!=NULL) {
     heff4->UseCurrentStyle(); heff4->GetYaxis()->SetLabelSize(lsize);
     heff4->Draw("colz");
     heff4->SetMaximum(1.5);
     heff4->GetZaxis()->SetRangeUser(0.5, 1.5);
     can->Update();    
     heff4->Draw("colz");
   }else{cout<<"Histogram "<<hname<<" not existing. "<<endl;}

 can->cd(8);
   if (heffS!=NULL) {
     heffS->UseCurrentStyle(); heffS->GetYaxis()->SetLabelSize(lsize);
     heffS->Draw("colz");
     heffS->SetMaximum(1.5);
     heffS->GetZaxis()->SetRangeUser(0.5, 1.5);
     heffS->Draw("colz");
   }else{cout<<"Histogram "<<hname<<" not existing. "<<endl;}

 can->cd(9);
   if (heff0!=NULL) {
     heff0->Draw("colz");
     heff0->SetMaximum(1.5);
     heff0->GetZaxis()->SetRangeUser(0.5, 1.5);
     heff0->UseCurrentStyle(); heff0->GetYaxis()->SetLabelSize(lsize);
     heff0->Draw("colz");
   }else{cout<<"Histogram "<<hname<<" not existing. "<<endl;}

 can->SaveAs("pl_eff_XY.pdf");

}