{

  TH1D *hped = new TH1D("hped","hped",500,0.,0.);
  TH1D *hped2 = new TH1D("hped2","hped2",500,0.,0.);

  TF1 *fitgaus = new TF1("fitgaus", "gaus",0.,10000.);

    for(int i1=0;i1<5000;i1++)
      hped2->Fill(gRandom->Gaus(8400.,40.));

    hped2->Fit(fitgaus,"","",8000.,8600.);
    double ped2 = fitgaus->GetParameter(1);
    double sigma2 = fitgaus->GetParameter(2);
    cout << Form("pass1 - ped2: %f sigma2: %f\n", ped2, sigma2);


    for(int i1=0;i1<500;i1++)
      hped->Fill(gRandom->Gaus(8400.,40.));

    hped->Fit(fitgaus,"","",8000.,8600.);
    double ped = fitgaus->GetParameter(1);
    double sigma = fitgaus->GetParameter(2);
    cout << Form("pass1 - ped: %f sigma: %f\n", ped, sigma);


}