Int_t npeaks = 30;
double fpeaks(double *x, double *par) {
   double result = par[0] + par[1]*x[0];
   for (Int_t p=0;p<npeaks;p++) {
      double norm  = par[3*p+2]; // "height" or "area"
      double mean  = par[3*p+3];
      double sigma = par[3*p+4];
      result += norm*TMath::Gaus(x[0],mean,sigma);
   }
   return result;
}

void peakstext(Int_t np=10) {
   npeaks = TMath::Abs(np);
   TH1F *h = new TH1F("h","test",500,0,1000);
   // Generate n peaks at random
   double par[3000];
   par[0] = 0.8;
   par[1] = -0.6/1000;

   for (int p=0;p<npeaks;p++) {
      par[3*p+2] = 1; // "height"
      par[3*p+3] = 10+gRandom->Rndm()*980; // "mean"
      par[3*p+4] = 3+2*gRandom->Rndm(); // "sigma"
   }

   TF1 *f = new TF1("f",fpeaks,0,1000,2+3*npeaks);
   f->SetNpx(1000);
   f->SetParameters(par);
   auto *c1 = new TCanvas();

   h->FillRandom("f",200000);

   TSpectrum *s = new TSpectrum(2*npeaks);
   Int_t nfound = s->Search(h,2,"",0.10);

   TList *functions = h->GetListOfFunctions();
   TPolyMarker *pm = (TPolyMarker*)functions->FindObject("TPolyMarker");
   double *xp = pm->GetX();
   double *yp = pm->GetY();
   int      n = pm->GetN();
   for (int i = 1; i<n; i++) {
      auto t = new TText(xp[i],yp[i],Form("%g", yp[i]));
      t->Draw();
   }
}