{
gROOT->Reset();
gROOT->LoadMacro("~/macros/fitfunction.cc+");
Double_t par[100];
Double_t par_err[100];

par[0] = .6;
par[1] = 0.1;
par[2] = 1;
par[3] = 1.;
par[4] = 150000;
par[5] = 25.1;
par[6] = 30000;
par[7] = 28.3;
par[8] = 15000;
par[9] = 31;
par[10] = 90000;
par[11] = 32.9;
par[12] = 1.05;
par[13] = 0;
par[14] = 0;

TH1F* hDe2 = new TH1F("hDe2","hDe2", 125, 15, 40);
TF1* fun = new TF1("fun", MyQuadroPeaks,15, 40, 15);
fun->SetParameters(par);
for(Int_t i=0;i<100000;i++)
{
  hDe2->Fill(fun->GetRandom(15., 40.));
}

FitMyQuadroPeaks(hDe2, par, par_err, 2., 2., 2., 2., 20., 2., 2., 2., -1, -1, 1, 2, 1, NULL);
}