2009-4-221粒子物理与核物理实验中的数据分析杨振伟清华大学第八讲增补:昀大似然拟合以及Minuit2009-4-222最大似然拟合与Minuit获得拟合参数。求极值,数,利用我们需要自己写似然函时,画出直方图或有的时候,比如不方便法拟合。的参数,选择昀大似然拟合。我们也可以修改函数,默认为的和中直接调用直方图或MinuitTGraph()FitFit()rsTGraphErroTGraphROOT2χhep.tsinghua.edu.cn/~yangzw/CourseDataAna/examples/minuit_fit.C得到拟合参数求出似然函数的极值,用写出似然函数个事例的分布产生按照就是这样的例子:比如Minuit)3)2200011).Cminuit_fit2βxαxf++=2009-4-223minuit_fit.C程序说明()minuit_fitvoid3)iflag)Int_tpar,*Double_tchi2,&Double_tgin,*Double_tnpar,&fcn(Int_tvoid2)data()get_input_void1)数中实现:以上三步分别在如下函类具体参数意义参见。参数个数,拟合次数等的参数,比如分布函数有时可能需要修改数从而获得自己的似然函替换原来的分布函数,用自己期望的分布函数组中即可。实验直接给出,存到数当然,这个分布经常是换成自己需要的分布。少:用户需要修改的地方很TMinuitMinuit()minuit_fitvoid3)iflag)Int_tpar,*Double_tchi2,&Double_tgin,*Double_tnpar,&fcn(Int_tvoid2)data()get_input_void1)