Fit example.
TFile** py-fillrandom.root
TFile* py-fillrandom.root
KEY: TFormula form1;1 abs(sin(x)/x)
KEY: TF1 sqroot;1 x*gaus(0) + [3]*form1
KEY: TH1F h1f;1 Test random numbers
Formula based function: sqroot
sqroot : x*gaus(0) + [3]*form1 Ndim= 1, Npar= 4, Number= 0
Formula expression:
x*[p0]*exp(-0.5*((x-[p1])/[p2])*((x-[p1])/[p2]))+[p3]*(abs(sin(x)/x))
FCN=198.935 FROM MIGRAD STATUS=CONVERGED 148 CALLS 149 TOTAL
EDM=2.98567e-07 STRATEGY= 1 ERROR MATRIX ACCURATE
EXT PARAMETER STEP FIRST
NO. NAME VALUE ERROR SIZE DERIVATIVE
1 p0 3.31658e+01 5.45703e-01 3.00376e-03 -1.11540e-03
2 p1 4.00667e+00 1.65304e-02 9.48491e-05 -3.06425e-02
3 p2 9.84663e-01 1.28238e-02 6.05976e-05 -3.04244e-02
4 p3 6.34464e+01 1.33233e+00 8.77483e-03 -3.96109e-04
fit1 : Real Time = 0.27 seconds Cpu Time = 0.27 seconds
3 from ROOT
import TCanvas, TFile, TPaveText
4 from ROOT
import gROOT, gBenchmark
6 c1 =
TCanvas(
'c1',
'The Fit Canvas', 200, 10, 700, 500 )
10 c1.GetFrame().SetBorderMode(-1 )
13 gBenchmark.Start(
'fit1' )
17 fill =
TFile(
'py-fillrandom.root' )
28 sqroot = gROOT.FindObject(
'sqroot' )
34 h1f = gROOT.FindObject(
'h1f' )
35 h1f.SetFillColor( 45 )
41 fitlabel =
TPaveText( 0.6, 0.3, 0.9, 0.80,
'NDC' )
42 fitlabel.SetTextAlign( 12 )
43 fitlabel.SetFillColor( 42 )
44 fitlabel.ReadFile(path.join(gROOT.GetTutorialsDir(),
'pyroot',
'fit1_py.py'))
47 gBenchmark.Show(
'fit1' )
- Author
- Wim Lavrijsen
Definition in file fit1.py.