int stmbgr(){

  TFile *f=new TFile("stm_bias_scan0.root");
  TH2F *h=f->Get("h2dadcdac001");
  TF1 *fg=new TF1("fg","gaus",0,100);
  fg->SetNumberFitPoints(100);
  h->FitSlicesX(fg,0,-1,0,"R");
  TCanvas *c1= new TCanvas("c1","c1");
  c1->Divide(1,2);
  c1->cd(1);
  TH1D *h2dadcdac001_1=(TH1D*)gDirectory->Get("h2dadcdac001_1");
  TH1D *h2dadcdac001_2=(TH1D*)gDirectory->Get("h2dadcdac001_2");
  h2dadcdac001_2->Draw();
  c1->cd(2);
  h2dadcdac001_1->Draw();
  for (int i=h2dadcdac001_1->GetNbinsX();i>0;i--){
    printf("%2.1f %2.1f %2.1f\n",
	   h2dadcdac001_1->GetBinCenter(i),
	   h2dadcdac001_1->GetBinContent(i),
	   h2dadcdac001_2->GetBinContent(i)
	   );
  }
  return 0;
}