void pl_ana_DYXY(Int_t iDut      = 901,
                 Int_t NewNbinsX = 4,
                 Int_t NewNbinsY = 4,
                 Int_t iOpt      = 0) {
  //  TCanvas *can = new TCanvas("can22","can22");
  //  can->Divide(2,2);
  TCanvas* can = new TCanvas("can", "can", 50, 0, 800, 800);
  can->Divide(2, 2);

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

  gROOT->cd();
  gROOT->SetDirLevel(1);

  TH1* h;
  TH1* h1;
  TH2* h2;
  TH3* h3;

  can->cd(1);
  gROOT->cd();
  TString hname = Form("hDutXYDY_%d", iDut);
  h3f           = (TH3*) gROOT->FindObjectAny(hname);
  if (h3f != NULL) {
    h3f->ProjectionZ()->Draw();
    gPad->SetLogy();
  } else {
    cout << hname << " not found" << endl;
  }

  //can->cd(2);Nxbins

  Int_t Nxbins = h3f->GetNbinsX();
  Int_t Nybins = h3f->GetNbinsY();
  Int_t Nzbins = h3f->GetNbinsZ();

  Int_t iDeltaX = Nxbins / NewNbinsX;
  Int_t iDeltaY = Nybins / NewNbinsY;
  cout << "Nbins: " << Nxbins << "," << Nybins << "," << Nzbins << " -> Dbin"
       << iDeltaX << "," << iDeltaY << endl;

  h2n = new TH2D("hTH3FitMean",
                 "#DeltaY Fit Entries; x (cm); y (cm)",
                 NewNbinsX,
                 h3f->GetXaxis()->GetXmin(),
                 h3f->GetXaxis()->GetXmax(),
                 NewNbinsY,
                 h3f->GetYaxis()->GetXmin(),
                 h3f->GetYaxis()->GetXmax());
  h2m = new TH2D("hTH3FitMean",
                 "#DeltaY Fit Mean; x (cm); y (cm)",
                 NewNbinsX,
                 h3f->GetXaxis()->GetXmin(),
                 h3f->GetXaxis()->GetXmax(),
                 NewNbinsY,
                 h3f->GetYaxis()->GetXmin(),
                 h3f->GetYaxis()->GetXmax());
  h2s = new TH2D("hTH3FitSigma",
                 "#DeltaY Fit Sigma; x (cm); y (cm)",
                 NewNbinsX,
                 h3f->GetXaxis()->GetXmin(),
                 h3f->GetXaxis()->GetXmax(),
                 NewNbinsY,
                 h3f->GetYaxis()->GetXmin(),
                 h3f->GetYaxis()->GetXmax());

  Int_t iCol = 0;
  for (Int_t ix = 0; ix < NewNbinsX; ix++) {
    for (Int_t iy = 0; iy < NewNbinsY; iy++) {
      TH1D* hsel      = h3f->ProjectionZ(Form("DeltaY_%d_%d", ix, iy),
                                    ix * iDeltaX + 1,
                                    (ix + 1) * iDeltaX + 1,
                                    iy * iDeltaY + 1,
                                    (iy + 1) * iDeltaY + 1);
      Double_t dFMean = hsel->GetMean();
      Double_t dFLim  = 2.5 * hsel->GetRMS();
      hsel->SetLineColor(iCol++);
      if (hsel->Integral() > 10) {
        TFitResultPtr fRes =
          hsel->Fit("gaus", "S", "HEsame", dFMean - dFLim, dFMean + dFLim);
        h2n->SetBinContent(ix + 1, iy + 1, hsel->Integral());
        h2m->SetBinContent(ix + 1, iy + 1, fRes->Parameter(1));
        h2s->SetBinContent(ix + 1, iy + 1, fRes->Parameter(2));
      }
    }
  }

  can->cd(2);
  {
    h2n->Draw("colz");
    //gPad->SetLinz();
    gPad->SetMargin(0.23, 0.23, 0.23, 0.23);
  }

  can->cd(3);
  {
    h2m->Draw("colz");
    gPad->SetMargin(0.23, 0.23, 0.23, 0.23);
  }

  can->cd(4);
  {
    h2s->Draw("colz");
    gPad->SetMargin(0.23, 0.23, 0.23, 0.23);
  }
  can->SaveAs("pl_ana_DYXY.pdf");
}