Logo ROOT  
Reference Guide
 
Loading...
Searching...
No Matches
fit2d.C File Reference

Detailed Description

View in nbviewer Open in SWAN
Example illustrating how to fit a 2-d histogram of type y=f(x)

****************************************
Minimizer is Linear / Migrad
Chi2 = 4501.58
NDf = 33
p0 = -0.47055 +/- 0.00713644
p1 = 10.8033 +/- 0.0124078
#include <TCanvas.h>
#include <TCutG.h>
#include <TH2F.h>
#include <TProfile.h>
#include <TRandom.h>
void fit2d()
{
// generate a 2-d histogram using a TCutG
const int n = 6;
float x[n] = {0.092,0.83,0.94,0.81,0.12,0.1};
float y[n] = {0.71,9.4,9,8,0.3,0.71};
TCutG *cut = new TCutG("cut",n,x,y);
TH2F *h2 = new TH2F("h2","h2",40,0,1,40,0,10);
float u,v;
for (int i=0;i<100000;i++) {
u = gRandom->Rndm();
v = 10*gRandom->Rndm();
if (cut->IsInside(u,v)) h2->Fill(u,v);
}
TCanvas *c1 = new TCanvas("c1","show profile",600,900);
c1->Divide(1,2);
c1->cd(1);
h2->Draw();
c1->cd(2);
//use a TProfile to convert the 2-d to 1-d problem
TProfile *prof = h2->ProfileX();
prof->Fit("pol1");
}
R__EXTERN TRandom * gRandom
Definition TRandom.h:62
The Canvas class.
Definition TCanvas.h:23
Graphical cut class.
Definition TCutG.h:20
virtual Int_t IsInside(Double_t x, Double_t y) const
Return 1 if the point (x,y) is inside the polygon defined by the graph vertices 0 otherwise.
Definition TGraph.cxx:1824
virtual TFitResultPtr Fit(const char *formula, Option_t *option="", Option_t *goption="", Double_t xmin=0, Double_t xmax=0)
Fit histogram with function fname.
Definition TH1.cxx:3898
2-D histogram with a float per channel (see TH1 documentation)
Definition TH2.h:307
Profile Histogram.
Definition TProfile.h:32
Double_t Rndm() override
Machine independent random number generator.
Definition TRandom.cxx:559
Double_t y[n]
Definition legend1.C:17
return c1
Definition legend1.C:41
Double_t x[n]
Definition legend1.C:17
const Int_t n
Definition legend1.C:16
Author
Rene Brun

Definition in file fit2d.C.