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



  1. #/bin/usr/env python  
  2.   
  3. import sys, math  
  4. import ROOT  
  5.   
  6.   
  7. if __name__=='__main__':  
  8.   
  9.     cv  = ROOT.TCanvas("cv""Histogram Example"20010700500)  
  10.   
  11.     # create histogram ( identify, Title, number of class, xmin, xmax )  
  12.     h1 = ROOT.TH1S('h''px'100, -5.5.)  
  13.   
  14.     for i in xrange(25000):  
  15.         px = ROOT.gRandom.Gaus()  
  16.         h1.Fill(px)  
  17.   
  18.     # fitting   
  19.     h1.Fit("gaus")  
  20.   
  21.     # output fitting parameters  
  22.     ROOT.gStyle.SetOptFit()  
  23.   
  24.     h1.Draw()  
  25.     cv.Update()  
  26.   
  27.   
  28.     # export  
  29.     cv.Print("histogram5.png")  



0 件のコメント:

コメントを投稿