{
  TH1F *h = new TH1F("h", "Beta", 150, 150., 300.);
  h->SetBinContent(52,600);
  h->SetBinContent(53,662);
  h->SetBinContent(54,622);
  h->SetBinContent(55,584);
  h->SetBinContent(56,633);
  h->SetBinContent(57,614);
  h->SetBinContent(58,622);
  h->SetBinContent(59,597);
  h->SetBinContent(60,589);
  h->SetBinContent(61,595);
  h->SetBinContent(62,615);
  h->SetBinContent(63,604);
  h->SetBinContent(64,607);
  h->SetBinContent(65,601);
  h->SetBinContent(66,588);
  h->SetBinContent(67,643);
  h->SetBinContent(68,584);
  h->SetBinContent(69,603);
  h->SetBinContent(70,642);
  h->SetBinContent(71,671);
  h->SetBinContent(72,682);
  h->SetBinContent(73,692);
  h->SetBinContent(74,724);
  h->SetBinContent(75,768);
  h->SetBinContent(76,742);
  h->SetBinContent(77,812);
  h->SetBinContent(78,837);
  h->SetBinContent(79,885);
  h->SetBinContent(80,881);
  h->SetBinContent(81,950);
  h->SetBinContent(82,1007);
  h->SetBinContent(83,1056);
  h->SetBinContent(84,1093);
  h->SetBinContent(85,1161);
  h->SetBinContent(86,1310);
  h->SetBinContent(87,1456);
  h->SetBinContent(88,1546);
  h->SetBinContent(89,1265);
  h->SetBinContent(90,1227);
  h->SetBinContent(91,1213);
  h->SetBinContent(92,1062);
  h->SetBinContent(93,1027);
  h->SetBinContent(94,947);
  h->SetBinContent(95,964);
  h->SetBinContent(96,817);
  h->SetBinContent(97,733);
  h->SetBinContent(98,780);
  h->SetBinContent(99,737);
  h->SetBinContent(100,714);
  h->SetBinContent(101,704);
  h->SetBinContent(102,693);
  h->SetBinContent(103,714);
  h->SetBinContent(104,604);
  h->SetBinContent(105,534);
  h->SetBinContent(106,487);
  h->SetBinContent(107,485);
  h->SetBinContent(108,445);
  h->SetBinContent(109,465);
  h->SetBinContent(110,408);
  h->SetBinContent(111,372);
  h->SetBinContent(112,354);
  h->SetBinContent(113,299);
  h->SetBinContent(114,298);
  h->SetBinContent(115,286);
  h->SetBinContent(116,293);
  h->SetBinContent(117,267);
  h->SetBinContent(118,232);
  h->SetBinContent(119,279);
  h->SetBinContent(120,221);
  h->SetBinContent(121,240);
  h->SetBinContent(122,236);
  h->SetBinContent(123,230);
  h->SetBinContent(124,221);
  h->SetBinContent(125,230);
  h->SetBinContent(126,197);
  h->SetBinContent(127,227);
  h->SetBinContent(128,225);
  h->SetBinContent(129,206);
  h->SetBinContent(130,199);
  h->SetEntries(51015);
  
#if 0 /* 0 or 1 */
  TF1 *f = new TF1("f","gaus(0)+gaus(3)+expo(6)", 201., 280.);
  // note: you MUST set reasonable initial values for all parameters
  // f->SetParameters(900., 235., 10., 200., 250., 10., 10., -0.02);
  f->SetParameters(700., 235., 10., 400., 235., 2., 10., -0.02);
#else /* 0 or 1 */
  TF1 *f = new TF1("f","gausn(0)+gausn(3)+expo(6)", 201., 280.);
  // note: you MUST set reasonable initial values for all parameters
  f->SetParameters(2.e4, 235., 10., 2.5e3, 235., 2., 10., -0.02);
#endif /* 0 or 1 */
  
  // gStyle->SetOptFit();
  
  // f->Print();
  h->Fit(f, "R");
}