import ROOT ROOT.gROOT.SetStyle("Plain") rndm = ROOT.TRandom3() rndm.SetSeed(123456) f1 = ROOT.TF1("f1", "pol0") chi2 = ROOT.TH1D("chi2", "Chi2 Prob", 100, 0, 1) chi2.SetMinimum(0) for j in range(1, 1000): h1 = ROOT.TH1D("h1", "rnd", 10, 0, 1) for i in range(0, 1000): h1.Fill(rndm.Uniform()) h1.Fit(f1, "Q0") chi2.Fill(ROOT.TMath.Prob(f1.GetChisquare(), f1.GetNDF())) h1 = None c = ROOT.TCanvas() chi2.Draw() c.SaveAs("02-chi2.pdf")