When calculating the chi2/NDF, empty bin is taken into account in the
max.likelihood method. Least-chi2 option don't have this problem. The
following macro demonstrate this. I tried this in ROOT 5.22. 5.27. They
shows the same results. Is this because I didn't use it properly or this
is a bug? If this is not a bug, how to do the max. likelihood fit
without taking into account the empty bins?
Thanks
--Wei
{
h = new TH1F("h", "", 100, -10, 10); for(int i = 0; i<1000; i++) { float a = gRandom->Gaus(0, 1); h->Fill(a); } h->Sumw2(); h->Scale(0.5); TH1F *h0 = (TH1F*)h->Clone(); h0->SetName("h0"); TH1F *h1 = (TH1F*)h->Clone(); h1->SetName("h1"); TH1F *h2 = (TH1F*)h->Clone(); h2->SetName("h2"); c1 = new TCanvas("c1", "", 900, 900); c1->Divide(2, 2); gStyle->SetOptFit(1); c1->cd(1); h->Fit("gaus", "L", "", -3, 3); c1->cd(2); h0->Fit("gaus", "L", "", -10, 10); c1->cd(3); h1->Fit("gaus", "", "", -3, 3); c1->cd(4); h2->Fit("gaus", "", "", -10, 10);
This archive was generated by hypermail 2.2.0 : Sat Dec 18 2010 - 17:50:01 CET