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")