void pl_XY_3D(Int_t iNSt=4, Int_t iOpt=0, Int_t i2D=0, Int_t i1D=0, Double_t dYmax=0.) {
  //  TCanvas *can = new TCanvas("can22","can22");
  //  can->Divide(2,2); 
  TCanvas *can = new TCanvas("can","can",50,0,800,800);
  can->Divide(1,1); 

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

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

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

 TString hname;
 TString cOpt;
 TString c2D;
 TString c1D;

 const Int_t ColMap[10]={1,2,3,4,6,7,8,9,41,46};

 Int_t iCan=1;
 gROOT->cd();
 switch(iOpt){
 case 0:
   cOpt="DX";
   break;
 case 1:
   cOpt="DY";
   break;
 case 2:
   cOpt="DT";
   break;
 case 3:
   cOpt="TOT";
   break;
 case 4:
   cOpt="CSZ";
   break;
 }

 switch(i2D){
   case 0:
     c2D="yx";
     break;
   case 1:
     c2D="zx";
     break;
   default:
     cout << "i2D mode not implemented " << endl;
     return;
 }

 can->cd(iCan);

 for (Int_t iSt=0; iSt<iNSt; iSt++){
   hname=Form("hXY_%s_%d",cOpt.Data(),iSt);
   h3f=(TH3 *)gROOT->FindObjectAny(hname);
   cout << hname.Data() <<" with pointer  "<<h3f<<" at iCan = "<<iCan<<endl;
   if (h3f!=NULL) {
     switch(i1D){
     case 0:
       h = (TH1 *)h3f->Project3DProfile(c2D.Data())->ProfileX();
       h->SetTitle(Form("%s",h->GetName()));
       //h->GetXaxis()->SetTitle("???");
       break;
     case 1:
       h = (TH1 *)h3f->Project3DProfile(c2D.Data())->ProfileY();
       h->SetTitle(Form("%s",h->GetName()));
       break;
     }
     h->SetLineColor(ColMap[iSt]);
     if(iSt==0) {
       h->Draw();
       if (dYmax>0.) {
	 h->SetMinimum(-dYmax);
	 h->SetMaximum(dYmax);
       }
     } 
     else {
       h->Draw("same");
     }
     gPad->SetGridx();
     gPad->SetGridy();
     //gPad->SetLogy();
   }else cout << hname << " not found" << endl;
 } 
 can->SaveAs("pl_XY_3D.pdf");
}