2009年7月17日金曜日

How to use pyROOT 基礎編 その3

ビルトイン関数で fitting

2. ヒストグラムの fitting

h1 = ROOT.TH1S(...)

・Gaussian fitting
h1.Fit("gaus")

・exponential fitting
h1.Fit("expo")

n次のべき関数 fitting
h1.Fit("poln")



fitting parameters を表示するには
ROOT.gStyle.SetOptFit()




#/bin/usr/env python

import sys, math
import ROOT


if __name__=='__main__':

cv = ROOT.TCanvas("cv", "Histogram Example", 200, 10, 700, 500)

# create histogram ( identify, Title, number of class, xmin, xmax )
h1 = ROOT.TH1S('h', 'px', 100, -5., 5.)

for i in xrange(25000):
px = ROOT.gRandom.Gaus()
h1.Fill(px)

# fitting
h1.Fit("gaus")

# output fitting parameters
ROOT.gStyle.SetOptFit()

h1.Draw()
cv.Update()


# export
cv.Print("histogram5.png")



0 件のコメント:

コメントを投稿