void pl_Eff_DTLH(Int_t iDut=910, Double_t dEffMin=0.5, Double_t dEffMax=1., TString sysinfo=""){
  gROOT->LoadMacro(((TString)gSystem->Getenv("VMCWORKDIR") + "/macro/beamtime/pl_Datime.C").Data());
  //  TCanvas *can = new TCanvas("can22","can22");
  //  can->Divide(2,2); 
  TCanvas *can = new TCanvas("can","can",48,55,600,600);
  can->Divide(2,2); 

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

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

 TH1 *h;
 TH1 *h1;
 TH1 *h1f;
 TH1 *h1m;
 TH1 *h1all;
 TH2 *h2;
 // if (hPla!=NULL) hPla->Delete();
 TString hname="";
 TProfile *h2pfx=NULL;

can->cd(1);
 Double_t Nfound=0.;
 Double_t Nmissed=0.;
 gROOT->cd();
 hname=Form("hDutDTLH_Found_%d",iDut);
 cout << " Look for histo " << hname << endl; 
 h1=(TH1 *)gROOT->FindObjectAny(hname);
 if (h1!=NULL) {
  h1->Draw();
  gPad->SetLogy();
  h1f = (TH1 *)h1->Clone();
  Nfound=h1f->GetEntries();
 }else  { cout << hname << " not found" << endl; return; }

 hname=Form("hDutDTLH_Missed_%d",iDut);
 h1=(TH1 *)gROOT->FindObjectAny(hname);
 if (h1!=NULL) {
  h1m = (TH1 *)h1->Clone();
  Nmissed=h1m->GetEntries();
  h1m->Draw("same");
  h1m->SetLineColor(2);
 }else  { cout << hname << " not found" << endl; return; }

can->cd(2);

 h1all = (TH1 *)h1f->Clone("hDutDTLH_all");
 h1all->Add(h1m,h1f,1.,1.);

 h1all->Draw();
 h1all->SetLineColor(1);
 h1m->Draw("same");
 h1f->Draw("same");

can->cd(3);

 TEfficiency *pEffDut = new TEfficiency(*h1f,*h1all);
 pEffDut->SetTitle("Efficiency of DUT");
 pEffDut->SetName("hDutDTLH_eff");
 pEffDut->Draw("AP");
 gPad->Update(); 
 
 auto graph = pEffDut->GetPaintedGraph(); 
 graph->SetMinimum(dEffMin);
 graph->SetMaximum(dEffMax); 
 graph->GetXaxis()->SetRangeUser(0.,12.);
 graph->Draw();
 /*
 auto heff = pEffDut->GetPaintedHistogram();
 heff->SetMinimum(dEffMin);
 heff->SetMaximum(dEffMax); 
 */
 gPad->Update();

 graph->GetXaxis()->SetRangeUser(0.,12.);
 gPad->SetGridx();
 gPad->SetGridy();

 /*
  TF1 *f1 = new TF1("Restore","1.-0.5*exp(-exp(x)*0.0001)",1.,10.);
  f1->Draw("same");
 */

 Double_t dEff = Nfound/(Nfound+Nmissed);
 cout<< "Average efficiency of Dut: " << dEff << endl;
 Double_t dDt = TMath::Power(10.,h1f->GetMean());

 cout<< "Average time distance of Dut-Hits: " << dDt << " ns => <rate> = "<< 10.E9/dDt <<" Hz / cell " << endl;

can->cd(4);
 TPaveText *pt = new TPaveText(0.1,0.1,0.9,0.9);
pt->AddText(
Form("average efficiency of %d: %5.3f",iDut,dEff));
pt->AddText(Form("average rate: %9.3f kHz/cell ",10.E6/dDt));
pt->SetFillColor(0);
pt->SetTextFont(52);
pt->SetBorderSize(0);
pt->Draw();

 TString FADD=Form("pl_Datime(\"%s\")",sysinfo.Data());
 if(gROOT->IsBatch())
 {
   gInterpreter->ProcessLine(FADD.Data());
 }


 can->SaveAs(Form("pl_Eff_DTLH_%d.pdf",iDut));
}