## From \$ROOTSYS/tutorials/fit/fit2.C

```#include "TF2.h"
#include "TH2.h"
#include "TMath.h"

// Fitting a 2-D histogram
// This tutorial illustrates :
//  - how to create a 2-d function
//  - fill a 2-d histogram randomly from this function
//  - fit the histogram
//  - display the fitted function on top of the histogram
//
// This example can be executed via the interpreter or ACLIC
//   root > .x fit2.C
//   root > .x fit2.C++
//Author: Rene Brun

Double_t g2(Double_t *x, Double_t *par) {
Double_t r1 = Double_t((x[0]-par[1])/par[2]);
Double_t r2 = Double_t((x[1]-par[3])/par[4]);
return par[0]*TMath::Exp(-0.5*(r1*r1+r2*r2));
}
Double_t fun2(Double_t *x, Double_t *par) {
Double_t *p1 = &par[0];
Double_t *p2 = &par[5];
Double_t *p3 = &par[10];
Double_t result = g2(x,p1) + g2(x,p2) + g2(x,p3);
return result;
}

void fit2() {
const Int_t npar = 15;
Double_t f2params[npar] =
{100,-3,3,-3,3,160,0,0.8,0,0.9,40,4,0.7,4,0.7};
TF2 *f2 = new TF2("f2",fun2,-10,10,-10,10, npar);
f2->SetParameters(f2params);

//Create an histogram and fill it randomly with f2
TH2F *h2 = new TH2F("h2","from f2",40,-10,10,40,-10,10);
Int_t nentries = 100000;
h2->FillRandom("f2",nentries);
//Fit h2 with original function f2
Float_t ratio = 4*nentries/100000;
f2params[ 0] *= ratio;
f2params[ 5] *= ratio;
f2params[10] *= ratio;
f2->SetParameters(f2params);
h2->Fit("f2");
f2->Draw("cont1 same");
}
```
fit2.C:1
fit2.C:2
fit2.C:3
fit2.C:4
fit2.C:5
fit2.C:6
fit2.C:7
fit2.C:8
fit2.C:9
fit2.C:10
fit2.C:11
fit2.C:12
fit2.C:13
fit2.C:14
fit2.C:15
fit2.C:16
fit2.C:17
fit2.C:18
fit2.C:19
fit2.C:20
fit2.C:21
fit2.C:22
fit2.C:23
fit2.C:24
fit2.C:25
fit2.C:26
fit2.C:27
fit2.C:28
fit2.C:29
fit2.C:30
fit2.C:31
fit2.C:32
fit2.C:33
fit2.C:34
fit2.C:35
fit2.C:36
fit2.C:37
fit2.C:38
fit2.C:39
fit2.C:40
fit2.C:41
fit2.C:42
fit2.C:43
fit2.C:44
fit2.C:45
fit2.C:46
fit2.C:47
fit2.C:48
fit2.C:49
fit2.C:50