void pl_over_Mat04D4best(Int_t iLog=0){
  //  TCanvas *can = new TCanvas("can22","can22");
  //  can->Divide(2,2); 
  TCanvas *can = new TCanvas("can","can",50,0,1000,800);
  can->Divide(5,5); 

gPad->SetFillColor(0);
gStyle->SetPalette(1);
gStyle->SetOptStat(kTRUE);
gStyle->SetOptFit(kTRUE);

 gROOT->cd();
 gROOT->SetDirLevel(1);
 // cout << " DirLevel "<< gROOT->GetDirLevel()<< endl;

 TH1 *h;
 TH2 *h2;
 Double_t NEvents;
 Double_t NFinalHits;
 Double_t dTRes;
 Double_t dTRMS;

 // if (hPla!=NULL) hPla->Delete();

can->cd(1);
 gROOT->cd();
 TString hname="hDXDY04D4best";
 h2=(TH2 *)gROOT->FindObjectAny(hname);
 if (h2!=NULL) {
  h2->Draw("colz");
  if(0==h2->GetEntries()) return;
  gPad->SetLogz();
  can->cd(2);
  h2px=h2->ProjectionX(); 
  TFitResultPtr fResx=h2px->Fit("gaus","S");
  dSigX=fResx->Parameter(2);
  h2py=h2->ProjectionY(); 
  TFitResultPtr fResy=h2py->Fit("gaus","S");
  dSigY=fResy->Parameter(2);
 }else {  cout << hname << " not found" << endl;  }

can->cd(2);
 gROOT->cd();
 TString hname="hDXDT04D4best";
 h2=(TH2 *)gROOT->FindObjectAny(hname);
 if (h2!=NULL) {
  h2->Draw("colz");
  gPad->SetLogz();
 }else  { cout << hname << " not found" << endl; }

can->cd(3);
 gROOT->cd();
 TString hname="hDYDT04D4best";
 h2=(TH2 *)gROOT->FindObjectAny(hname);
 if (h2!=NULL) {
  h2->Draw("colz");
  gPad->SetLogz();
 }else  { cout << hname << " not found" << endl; }

can->cd(4);
 gROOT->cd();
 TString hname="hDTD4DT04D4best";
 h2=(TH2 *)gROOT->FindObjectAny(hname);
 if (h2!=NULL) {
  h2->Draw("colz");
  gPad->SetLogz();
  TProfile *h2pfx=h2->ProfileX();
  h2pfx->Draw("same");

  TString hname2="hDTD4DT04D4best_pfx_px";
  h2pfxpx=(TH2 *)gROOT->FindObjectAny(hname2);
  if (h2pfxpx!=NULL) h2pfxpx->Draw("same");

  Int_t iNbinsX=h2->GetNbinsX();
  cout<< " Histo "<<hname<<" has "<< iNbinsX << " Xbins"<<endl;
  const char *cnamP=Form("%s_py",(char *)hname);
  TH1D *hDT04=h2->ProjectionY(cnamP,0,iNbinsX-1);

  /*
  Double_t BRes=1.E8;
  Int_t BL=0;
  Double_t BEntries=0.;
  Double_t Brms=0.;
  Double_t BNorm=0.;
  for (Int_t iBinL=0; iBinL<iNbinsX/2-1; iBinL++){
      const char *cnam=Form("%s_py%d",(char *)hname,iBinL);
      TH1D *hDT04tmp=h2->ProjectionY(cnam,iBinL,iNbinsX-iBinL,"");
      TFitResultPtr fRes=hDT04tmp->Fit("gaus","S");
      cout <<Form("Fit with iBinL = %d : %6.1f, RMS : %6.1f ",iBinL,fRes->Parameter(2),hDT04tmp->GetRMS());
      cout<<endl;
      if(fRes->Parameter(2)<BRes){
	BRes=fRes->Parameter(2);
	BNorm=fRes->Parameter(0);
	Brms=hDT04tmp->GetRMS();
	BEntries=hDT04tmp->GetSum();
	BL=iBinL;
      }
      if (hDT04tmp->GetEntries()<1000) break;
    }
    //    gPad->SetOptFit(1);
    */
 }else  { cout << hname << " not found" << endl; }


can->cd(5);
 gROOT->cd();
 TString hname="hChi04D4best";
 h1=(TH1 *)gROOT->FindObjectAny(hname);
 if (h1!=NULL) {
  h1->Draw();
  NFinalHits=h1->GetEntries();
 }else  { cout << hname << " not found" << endl; }

can->cd(6);
 gROOT->cd();
 TString hname="hNMatchD4sel";  // normalization
 h1=(TH1 *)gROOT->FindObjectAny(hname);
 if (h1!=NULL) {
  h1->Draw();
  gPad->SetLogy();
  NEvents=h1->GetEntries();
 }else  { cout << hname << " not found" << endl; }

can->cd(7);
 gROOT->cd();
 TFitResultPtr fRes=hDT04->Fit("gaus","S");
 // fRes->Print();
 dTMean=fRes->Parameter(1);
 dTRes=fRes->Parameter(2);
 dTRMS=hDT04->GetRMS();

can->cd(8);
 gROOT->cd();
 TString hname="hDTD4sel";
 h1=(TH1 *)gROOT->FindObjectAny(hname);
 if (h1!=NULL) {
  h1->Draw();
  //  Double_t NSel=h1->Integral(BL,h1->GetNbinsX(),"");
 }else  { cout << hname << " not found" << endl; }

can->cd(9);
 gROOT->cd();
 // TString hname="hDistDT04D4best";
 TString hname="hTexpDT04D4best";
 h2=(TH2 *)gROOT->FindObjectAny(hname);
 if (h2!=NULL) {
  h2->Draw("colz");
  h2->ProfileX()->Draw("same");
  gPad->SetLogz();
 }else  { cout << hname << " not found" << endl; }

can->cd(10);
 gROOT->cd();
 TString hname="hDTMul4D4best";
 h2=(TH2 *)gROOT->FindObjectAny(hname);
 if (h2!=NULL) {
  h2->Draw("colz");
  gPad->SetLogz();
 }else  { cout << hname << " not found" << endl; }

can->cd(11);
 gROOT->cd();
 TString hname="hDTX4D4best";
 h2=(TH2 *)gROOT->FindObjectAny(hname);
 if (h2!=NULL) {
  h2->Draw("colz");
  gPad->SetLogz();
  TProfile *h2pfx=h2->ProfileX();
  h2pfx->Draw("same");
 }else  { cout << hname << " not found" << endl; }

can->cd(12);
 gROOT->cd();
 TString hname="hDTY4D4best";
 h2=(TH2 *)gROOT->FindObjectAny(hname);
 if (h2!=NULL) {
  h2->Draw("colz");
  TProfile *h2pfxx=h2->ProfileX();
  h2pfxx->Draw("same");
  gPad->SetLogz();
 }else  { cout << hname << " not found" << endl; }

can->cd(13);
 gROOT->cd();
 TString hname="hDTMul0D4best";
 h2=(TH2 *)gROOT->FindObjectAny(hname);
 if (h2!=NULL) {
  h2->Draw("colz");
  gPad->SetLogz();
 }else  { cout << hname << " not found" << endl; }

can->cd(14);
 gROOT->cd();
 TString hname="hDTX0D4best";
 h2=(TH2 *)gROOT->FindObjectAny(hname);
 if (h2!=NULL) {
  h2->Draw("colz");
  gPad->SetLogz();
 }else  { cout << hname << " not found" << endl; }

can->cd(15);
 gROOT->cd();
 TString hname="hDTY0D4best";
 h2=(TH2 *)gROOT->FindObjectAny(hname);
 if (h2!=NULL) {
  h2->Draw("colz");
  gPad->SetLogz();
 }else  { cout << hname << " not found" << endl; }

can->cd(16);
 gROOT->cd();
 TString hname="hCluSize0DT04D4best";
 h2=(TH2 *)gROOT->FindObjectAny(hname);
 if (h2!=NULL) {
  h2->Draw("colz");
  gPad->SetLogz();
  TH1D* h2clu=h2->ProjectionX("_px",0,h2->GetNbinsY()-1);
  dMeanCluSize0=h2clu->GetMean();
 }else  { cout << hname << " not found" << endl; }

can->cd(17);
 gROOT->cd();
 TString hname="hTot0DT04D4best";
 h2=(TH2 *)gROOT->FindObjectAny(hname);
 if (h2!=NULL) {
  h2->Draw("colz");
  TProfile *h2pfxt=h2->ProfileX();
  h2pfxt->Draw("same");
  gPad->SetLogz();
 }else  { cout << hname << " not found" << endl; }

can->cd(18);
 gROOT->cd();
 TString hname="hDXX4D4best";
 h2=(TH2 *)gROOT->FindObjectAny(hname);
 if (h2!=NULL) {
  h2->Draw("colz");
  gPad->SetLogz();
 }else  { cout << hname << " not found" << endl; }

can->cd(19);
 gROOT->cd();
 TString hname="hDXY4D4best";
 h2=(TH2 *)gROOT->FindObjectAny(hname);
 if (h2!=NULL) {
  h2->Draw("colz");
  gPad->SetLogz();
 }else  { cout << hname << " not found" << endl; }

can->cd(20);
 gROOT->cd();
 TString hname="hDYX4D4best";
 h2=(TH2 *)gROOT->FindObjectAny(hname);
 if (h2!=NULL) {
  h2->Draw("colz");
  gPad->SetLogz();
 }else  { cout << hname << " not found" << endl; }

can->cd(21);
 gROOT->cd();
 TString hname="hDYY4D4best";
 h2=(TH2 *)gROOT->FindObjectAny(hname);
 if (h2!=NULL) {
  h2->Draw("colz");
  gPad->SetLogz();
 }else  { cout << hname << " not found" << endl; }

can->cd(22);
 gROOT->cd();
 TString hname="hDXX0D4best";
 h2=(TH2 *)gROOT->FindObjectAny(hname);
 if (h2!=NULL) {
  h2->Draw("colz");
  gPad->SetLogz();
 }else  { cout << hname << " not found" << endl; }

can->cd(23);
 gROOT->cd();
 TString hname="hDXY0D4best";
 h2=(TH2 *)gROOT->FindObjectAny(hname);
 if (h2!=NULL) {
  h2->Draw("colz");
  gPad->SetLogz();
 }else  { cout << hname << " not found" << endl; }

can->cd(24);
 gROOT->cd();
 TString hname="hDYX0D4best";
 h2=(TH2 *)gROOT->FindObjectAny(hname);
 if (h2!=NULL) {
  h2->Draw("colz");
  gPad->SetLogz();
 }else  { cout << hname << " not found" << endl; }

can->cd(25);
 gROOT->cd();
 TString hname="hDYY0D4best";
 h2=(TH2 *)gROOT->FindObjectAny(hname);
 if (h2!=NULL) {
  h2->Draw("colz");
  gPad->SetLogz();
 }else  { cout << hname << " not found" << endl; }

 if (iLog != 0) {
   .> "CalLog.out" ;
 }
 //report summary
 Double_t dEff=NFinalHits/NEvents;
 cout << Form("<IRes> Efficiency: %6.3f of %7.0f events, mean Dt %6.1f ps, 1-sig time res: %6.1f ps, RMS: %6.1f ps, SigX %5.2f, SigY %5.2f, CluSize %4.1f",
	      dEff,NEvents, dTMean, dTRes,dTRMS,dSigX,dSigY,dMeanCluSize0) <<endl; 

 // cout << Form("    D4best 1-sigma timing resolution at lower index %d: %6.1f ps, RMS: %6.1f ps from  %6.0f entries of %6.0f (%6.3f)",BL,BRes,Brms,BEntries,NSel,BEntries/NSel) <<endl; 

 can->SaveAs("pl_over_Mat04D4best.pdf");
}