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

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

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

 TH1 *h;
 TH1 *h1;
 TH1 *hDT04;
 TH1 *h2px;
 TH1 *h2py;
 TH1 *h2pfx;
 TH1 *h2pfxpx;
 TH1 *h2pfxx;
 TH1 *h2pfxt;
 TH1 *h2clu;

 TH2 *h2;
 
 Double_t NEvents;
 Double_t NFinalHits;
 Double_t dTRes;
 Double_t dTRMS;
 Double_t dTMean;
 Double_t dMeanX;
 Double_t dSigX;
 Double_t dMeanY;
 Double_t dSigY;
 Double_t dMeanCluSize0;
 Double_t dMeanCluSize4;

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

can->cd(1);
{
 gROOT->cd();
 TString hname="hDXDY04D4best";
 h2=(TH2 *)gROOT->FindObjectAny(hname);
 if( NULL == h2 ) return;
 if(100 > h2->GetEntries()) {
   cout << "Insufficient number of events " << h2->GetEntries() << endl;
  return;
 }
 if (h2!=NULL) {
  h2->Draw("colz");

  gPad->SetLogz();
  can->cd(2);
  h2px=h2->ProjectionX(); 
  
  TFitResultPtr fResx=h2px->Fit("gaus","S","",-2.*h2px->GetRMS(),2.*h2px->GetRMS());
  if( static_cast<Int_t>(fResx) )
  {
    cout << "res X fit failed" <<endl;
  }
  else
  {
    dMeanX=fResx->Parameter(1);
    dSigX=fResx->Parameter(2);
  }
  h2py=h2->ProjectionY(); 
  TFitResultPtr fResy=h2py->Fit("gaus","S","",-2.*h2py->GetRMS(),2.*h2py->GetRMS());
  if( static_cast<Int_t>(fResy) )
  {
    cout << "res Y fit failed" <<endl;
  }
  else
  {
    dMeanY=fResy->Parameter(1);
    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=(TH1 *)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.Data());
  hDT04=(TH1D *)h2->ProjectionY(cnamP,1,iNbinsX);
 }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();
 hDT04->Draw();
 gPad->SetLogy();
 TFitResultPtr fRes=hDT04->Fit("gaus","S");
 // fRes->Print();
  if( static_cast<Int_t>(fRes) )
  {
    cout << "res T fit failed" <<endl;
  }
  else
  {
    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",1,h2->GetNbinsY());
  dMeanCluSize0=h2clu->GetMean();
  h2->ProfileX()->Draw("same");
 }else  { cout << hname << " not found" << endl; }
}

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

can->cd(18);
{
  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(19);
{
  gROOT->cd();
 TString hname="hTot4DT04D4best";
 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(20);
{
  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(21);
{
  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(22);
{
  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(23);
{
  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(24);
{
  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(25);
{
  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(26);
{
  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(27);
{
  gROOT->cd();
 TString hname="hDYY0D4best";
 h2=(TH2 *)gROOT->FindObjectAny(hname);
 if (h2!=NULL) {
  h2->Draw("colz");
  gPad->SetLogz();
 }else  { cout << hname << " not found" << endl; }
}

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

can->cd(29);
{
  gROOT->cd();
 TString hname="hVelD4best";
 h2=(TH2 *)gROOT->FindObjectAny(hname);
 if (h2!=NULL) {
  h2->Draw("");
  gPad->SetLogy();
 }else  { cout << hname << " not found" << endl; }
}
 //report summary
 Double_t dEff=NFinalHits/NEvents;
 cout << Form("<IRes> Efficiency: %6.3f, %7.0f of %7.0f events, mean Dt %6.4f ns, 1-sig time res: %6.4f ns, RMS: %6.4f ns, DX %5.2f, DY %5.2f, SigX %5.2f, SigY %5.2f, CluSize %4.1f,%4.1f ",
	      dEff,NFinalHits,NEvents,dTMean,dTRes,dTRMS,dMeanX,dMeanY,dSigX,dSigY,dMeanCluSize0,dMeanCluSize4) <<endl; 

 // cout << Form("    D4best 1-sigma timing resolution at lower index %d: %6.1f ns, RMS: %6.1f ns from  %6.0f entries of %6.0f (%6.3f)",BL,BRes,Brms,BEntries,NSel,BEntries/NSel) <<endl; 
 if (iLog != 0) 
 gROOT->ProcessLine(Form(".! echo %f  > Test.res", dTRes)); 

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