ROOT_for_beginners_Day3.pdfVIP

ROOT for beginners Third day Data Fitting 1 Fits • We know how to make a fit by using the graphical interface... • How to fit with the command line? Or within a script ? gStyle-SetOptFit(kTRUE) TH1F *h = new TH1F(hg,Un exemple de fit,100,-2,2) h-FillRandom (gaus,10000) h-Fit(gaus,V,E1,-1,1.5) function name fit options fit limits drawing options 2 Basic fits… 3 Which fitting functions ? • The predefined functions: – gaus = p0*exp(-0.5*pow ((x-p1)/p2),2) – expo = exp(p0+p1*x) – polN = p0 + p1*x + p2*pow (x,2) + p3*... – landau (guess the formula!) • How to obtain the values of the fit parameters ? TF1 *gfit = (TF1 *)h-GetFunction(gaus) gfit-GetParameter(0) gfit-GetParameter(1) ... gfit-GetParError(0) ... double par[3] gfit-GetParameters(par) 4 Creating a user defined function TF1 *fu = new TF1(f1, sin(x)/x, 0, 10) TF1 *fd = new TF1(f2, f1 * 2, 0, 10) fu-Draw () Only thef unction name is known! fd-Draw (same) And many other combinations ! 5 Including parameters TF1 *ft = new TF1(f3, [0]*sqrt([1]*(sin([2]*x)+2)) /([3]+pow (x,2)), 0, 10) ft-SetParameters(1,1,3,5) ft-Draw () index 0 1 2 3 content 1 1 3 5 TH1F *hd = new TH1F(h2,Un exemple, 100, 0, 10) hd-FillRandom(f3,100000) ft-SetParameters (h2-GetMaximum(),1,2.