{ "cells": [ { "cell_type": "markdown", "id": "69f511eb", "metadata": {}, "source": [ "# TestBinomial\n", "Perform a fit to a set of data with binomial errors\n", "like those derived from the division of two histograms.\n", "Three different fits are performed and compared:\n", "\n", " - simple least square fit to the divided histogram obtained\n", " from TH1::Divide with option b\n", " - least square fit to the TGraphAsymmErrors obtained from\n", " TGraphAsymmErrors::BayesDivide\n", " - likelihood fit performed on the dividing histograms using\n", " binomial statistics with the TBinomialEfficiency class\n", "\n", "The first two methods are biased while the last one is statistical correct.\n", "Running the script passing an integer value n larger than 1, n fits are\n", "performed and the bias are also shown.\n", "To run the script :\n", "\n", " to show the bias performing 100 fits for 1000 events per \"experiment\"\n", "\n", "```cpp\n", " root[0]: .x TestBinomial.C+\n", "```\n", "\n", " to show the bias performing 100 fits for 1000 events per \"experiment\"\n", "\n", "```cpp\n", " .x TestBinomial.C+(100, 1000)\n", "```\n", "\n", "\n", "\n", "\n", "**Author:** Rene Brun \n", "This notebook tutorial was automatically generated with ROOTBOOK-izer from the macro found in the ROOT repository on Tuesday, May 19, 2026 at 08:25 PM." ] }, { "cell_type": "markdown", "id": "f9177c2d", "metadata": {}, "source": [ " Arguments are defined. " ] }, { "cell_type": "code", "execution_count": 1, "id": "1d26b8ac", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:25:21.107781Z", "iopub.status.busy": "2026-05-19T20:25:21.107665Z", "iopub.status.idle": "2026-05-19T20:25:21.426455Z", "shell.execute_reply": "2026-05-19T20:25:21.425812Z" } }, "outputs": [], "source": [ "int nloop = 100;\n", "int nevts = 100;\n", "bool plot = false;\n", "bool debug = false;\n", "int seed = 111;" ] }, { "cell_type": "code", "execution_count": 2, "id": "649a4f04", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:25:21.428439Z", "iopub.status.busy": "2026-05-19T20:25:21.428322Z", "iopub.status.idle": "2026-05-19T20:25:21.656032Z", "shell.execute_reply": "2026-05-19T20:25:21.645661Z" } }, "outputs": [], "source": [ "gStyle->SetMarkerStyle(20);\n", "gStyle->SetLineWidth(2.0);\n", "gStyle->SetOptStat(11);\n", "\n", "TObjArray hbiasNorm;\n", "hbiasNorm.Add(new TH1D(\"h0Norm\", \"Bias Histogram fit\",100,-5,5));\n", "hbiasNorm.Add(new TH1D(\"h1Norm\",\"Bias Binomial fit\",100,-5,5));\n", "TObjArray hbiasThreshold;\n", "hbiasThreshold.Add(new TH1D(\"h0Threshold\", \"Bias Histogram fit\",100,-5,5));\n", "hbiasThreshold.Add(new TH1D(\"h1Threshold\",\"Bias Binomial fit\",100,-5,5));\n", "TObjArray hbiasWidth;\n", "hbiasWidth.Add(new TH1D(\"h0Width\", \"Bias Histogram fit\",100,-5,5));\n", "hbiasWidth.Add(new TH1D(\"h1Width\",\"Bias Binomial fit\",100,-5,5));\n", "TH1D* hChisquared = new TH1D(\"hChisquared\",\n", " \"#chi^{2} probability (Baker-Cousins)\", 200, 0.0, 1.0);\n", "\n", "TVirtualFitter::SetDefaultFitter(\"Minuit2\");\n", "ROOT::Math::IntegratorOneDimOptions::SetDefaultIntegrator(\"Gauss\");" ] }, { "cell_type": "markdown", "id": "92496308", "metadata": {}, "source": [ "Note: in order to be able to use TH1::FillRandom() to generate\n", "pseudo-experiments, we use a trick: generate \"selected\"\n", "and \"non-selected\" samples independently. These are\n", "statistically independent and therefore can be safely\n", "added to yield the \"before selection\" sample." ] }, { "cell_type": "markdown", "id": "faf4f564", "metadata": {}, "source": [ "Define (arbitrarily?) a distribution of input events.\n", "Here: assume a x^(-2) distribution. Boundaries: [10, 100]." ] }, { "cell_type": "code", "execution_count": 3, "id": "b5bfaeea", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:25:21.666089Z", "iopub.status.busy": "2026-05-19T20:25:21.665940Z", "iopub.status.idle": "2026-05-19T20:25:21.872358Z", "shell.execute_reply": "2026-05-19T20:25:21.871746Z" } }, "outputs": [], "source": [ "double xmin =10, xmax = 100;\n", "TH1D* hM2D = new TH1D(\"hM2D\", \"x^(-2) denominator distribution\",\n", " 45, xmin, xmax);\n", "TH1D* hM2N = new TH1D(\"hM2N\", \"x^(-2) numerator distribution\",\n", " 45, xmin, xmax);\n", "TH1D* hM2E = new TH1D(\"hM2E\", \"x^(-2) efficiency\",\n", " 45, xmin, xmax);\n", "\n", "TF1* fM2D = new TF1(\"fM2D\", \"(1-[0]/(1+exp(([1]-x)/[2])))/(x*x)\",\n", " xmin, xmax);\n", "TF1* fM2N = new TF1(\"fM2N\", \"[0]/(1+exp(([1]-x)/[2]))/(x*x)\",\n", " xmin, xmax);\n", "TF1* fM2Fit = new TF1(\"fM2Fit\", \"[0]/(1+exp(([1]-x)/[2]))\",\n", " xmin, xmax);\n", "TF1* fM2Fit2 = nullptr;\n", "\n", "TRandom3 rb(seed);" ] }, { "cell_type": "markdown", "id": "bd47a260", "metadata": {}, "source": [ "First try: use a single set of parameters.\n", "For each try, we need to find the overall normalization" ] }, { "cell_type": "code", "execution_count": 4, "id": "6987a914", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:25:21.874152Z", "iopub.status.busy": "2026-05-19T20:25:21.874031Z", "iopub.status.idle": "2026-05-19T20:25:22.081416Z", "shell.execute_reply": "2026-05-19T20:25:22.080852Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "32 68\n" ] } ], "source": [ "double normalization = 0.80;\n", "double threshold = 25.0;\n", "double width = 5.0;\n", "\n", "fM2D->SetParameter(0, normalization);\n", "fM2D->SetParameter(1, threshold);\n", "fM2D->SetParameter(2, width);\n", "fM2N->SetParameter(0, normalization);\n", "fM2N->SetParameter(1, threshold);\n", "fM2N->SetParameter(2, width);\n", "double integralN = fM2N->Integral(xmin, xmax);\n", "double integralD = fM2D->Integral(xmin, xmax);\n", "double fracN = integralN/(integralN+integralD);\n", "int nevtsN = rb.Binomial(nevts, fracN);\n", "int nevtsD = nevts - nevtsN;\n", "\n", "std::cout << nevtsN << \" \" << nevtsD << std::endl;\n", "\n", "gStyle->SetOptFit(1111);" ] }, { "cell_type": "markdown", "id": "863b7c85", "metadata": {}, "source": [ "generate many times to see the bias" ] }, { "cell_type": "code", "execution_count": 5, "id": "9974f8e5", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:25:22.083164Z", "iopub.status.busy": "2026-05-19T20:25:22.083042Z", "iopub.status.idle": "2026-05-19T20:25:22.717376Z", "shell.execute_reply": "2026-05-19T20:25:22.716798Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 0.350714\n", "NDf = 3\n", "Edm = 1.27001e-06\n", "NCalls = 81\n", "p0 = 0.694132 +/- 0.210029 \n", "p1 = 19.3471 +/- 5.85483 \n", "p2 = 5.2245 +/- 5.11013 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 14.3081\n", "Chi2 = 28.6162\n", "NDf = 21\n", "Edm = 4.11934e-08\n", "NCalls = 91\n", "p0 = 0.822966 +/- 0.0807557 \t (limited)\n", "p1 = 21.9819 +/- 2.14579 \n", "p2 = 3.36105 +/- 1.05866 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 1.64738\n", "NDf = 5\n", "Edm = 1.01171e-06\n", "NCalls = 137\n", "p0 = 0.677329 +/- 0.135567 \t (limited)\n", "p1 = 15.5028 +/- 4.11118 \n", "p2 = 4.10441 +/- 2.77044 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 19.363\n", "Chi2 = 38.7261\n", "NDf = 25\n", "Edm = 9.40378e-08\n", "NCalls = 80\n", "p0 = 0.778188 +/- 0.0815568 \t (limited)\n", "p1 = 23.7948 +/- 2.21356 \n", "p2 = 4.58769 +/- 1.24153 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 0.928456\n", "NDf = 6\n", "Edm = 1.46917e-08\n", "NCalls = 75\n", "p0 = 0.578806 +/- 0.128412 \t (limited)\n", "p1 = 19.8914 +/- 2.83678 \n", "p2 = 3.39323 +/- 1.96156 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 11.7144\n", "Chi2 = 23.4287\n", "NDf = 23\n", "Edm = 7.44806e-08\n", "NCalls = 73\n", "p0 = 0.740898 +/- 0.079731 \t (limited)\n", "p1 = 21.5953 +/- 2.12198 \n", "p2 = 2.98575 +/- 1.16957 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 2.56551\n", "NDf = 7\n", "Edm = 1.19556e-08\n", "NCalls = 109\n", "p0 = 0.57261 +/- 0.16552 \t (limited)\n", "p1 = 21.4104 +/- 5.10611 \n", "p2 = 5.58786 +/- 3.30889 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 12.7185\n", "Chi2 = 25.4369\n", "NDf = 23\n", "Edm = 7.33141e-08\n", "NCalls = 74\n", "p0 = 0.861581 +/- 0.082929 \t (limited)\n", "p1 = 28.1952 +/- 3.60519 \n", "p2 = 7.31775 +/- 2.14695 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 7.37095\n", "NDf = 10\n", "Edm = 1.54804e-06\n", "NCalls = 87\n", "p0 = 0.667689 +/- 0.152841 \t (limited)\n", "p1 = 24.0146 +/- 2.40617 \n", "p2 = 3.63777 +/- 1.74672 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 9.05615\n", "Chi2 = 18.1123\n", "NDf = 23\n", "Edm = 2.30749e-07\n", "NCalls = 71\n", "p0 = 0.842186 +/- 0.0921999 \t (limited)\n", "p1 = 26.3581 +/- 2.99562 \n", "p2 = 5.28694 +/- 1.78724 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 0.747776\n", "NDf = 6\n", "Edm = 2.75121e-06\n", "NCalls = 76\n", "p0 = 0.651877 +/- 0.155245 \t (limited)\n", "p1 = 18.9865 +/- 2.68321 \n", "p2 = 2.68989 +/- 1.28756 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 15.9733\n", "Chi2 = 31.9466\n", "NDf = 26\n", "Edm = 5.1941e-07\n", "NCalls = 62\n", "p0 = 0.794215 +/- 0.0873938 \t (limited)\n", "p1 = 25.0012 +/- 3.13976 \n", "p2 = 5.16842 +/- 1.65839 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 0.433039\n", "NDf = 4\n", "Edm = 3.07902e-07\n", "NCalls = 104\n", "p0 = 0.500967 +/- 0.254652 \t (limited)\n", "p1 = 19.4199 +/- 9.30614 \n", "p2 = 6.46947 +/- 8.16624 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 17.807\n", "Chi2 = 35.614\n", "NDf = 24\n", "Edm = 5.43825e-09\n", "NCalls = 98\n", "p0 = 0.748311 +/- 0.0866402 \t (limited)\n", "p1 = 22.4098 +/- 2.06254 \n", "p2 = 3.36062 +/- 0.902421 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 4.49197\n", "NDf = 9\n", "Edm = 4.04441e-06\n", "NCalls = 113\n", "p0 = 0.623236 +/- 0.107197 \t (limited)\n", "p1 = 20.8129 +/- 1.0507 \n", "p2 = 1.22179 +/- 1.09027 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 10.5454\n", "Chi2 = 21.0907\n", "NDf = 22\n", "Edm = 1.84998e-08\n", "NCalls = 101\n", "p0 = 0.7513 +/- 0.0923543 \t (limited)\n", "p1 = 23.3695 +/- 3.23382 \n", "p2 = 3.90837 +/- 1.83804 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 2.54966\n", "NDf = 6\n", "Edm = 1.48922e-07\n", "NCalls = 93\n", "p0 = 0.382687 +/- 0.226798 \t (limited)\n", "p1 = 18.0539 +/- 6.38962 \n", "p2 = 3.0464 +/- 5.58751 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 7.09023\n", "Chi2 = 14.1805\n", "NDf = 22\n", "Edm = 1.88451e-08\n", "NCalls = 149\n", "p0 = 1 +/- 0.132534 \t (limited)\n", "p1 = 29.3106 +/- 2.29361 \n", "p2 = 6.14658 +/- 1.45195 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 1.46531\n", "NDf = 7\n", "Edm = 2.36369e-06\n", "NCalls = 90\n", "p0 = 0.575891 +/- 0.187613 \t (limited)\n", "p1 = 17.6126 +/- 7.68638 \n", "p2 = 2.99185 +/- 6.18408 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 15.9432\n", "Chi2 = 31.8864\n", "NDf = 23\n", "Edm = 9.70221e-07\n", "NCalls = 76\n", "p0 = 0.765903 +/- 0.0805873 \t (limited)\n", "p1 = 23.9216 +/- 2.50269 \n", "p2 = 4.42272 +/- 1.34342 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 3.31728\n", "NDf = 6\n", "Edm = 8.18927e-07\n", "NCalls = 64\n", "p0 = 0.59853 +/- 0.124912 \t (limited)\n", "p1 = 17.9048 +/- 3.49062 \n", "p2 = 3.52241 +/- 2.09942 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 21.6362\n", "Chi2 = 43.2724\n", "NDf = 26\n", "Edm = 2.82218e-06\n", "NCalls = 80\n", "p0 = 0.64315 +/- 0.0777258 \t (limited)\n", "p1 = 19.2211 +/- 1.75093 \n", "p2 = 2.29991 +/- 0.890814 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 0.413866\n", "NDf = 6\n", "Edm = 5.34616e-09\n", "NCalls = 82\n", "p0 = 0.654352 +/- 0.198145 \t (limited)\n", "p1 = 20.6178 +/- 3.66715 \n", "p2 = 4.03019 +/- 2.20627 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 9.09822\n", "Chi2 = 18.1964\n", "NDf = 19\n", "Edm = 2.73834e-07\n", "NCalls = 71\n", "p0 = 0.809718 +/- 0.0969513 \t (limited)\n", "p1 = 22.5716 +/- 2.36934 \n", "p2 = 3.95798 +/- 1.33784 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 1.21002\n", "NDf = 4\n", "Edm = 2.90565e-08\n", "NCalls = 82\n", "p0 = 0.651685 +/- 0.159188 \t (limited)\n", "p1 = 19.3369 +/- 2.66645 \n", "p2 = 2.81842 +/- 2.20413 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 16.871\n", "Chi2 = 33.7421\n", "NDf = 24\n", "Edm = 9.78264e-08\n", "NCalls = 73\n", "p0 = 0.780957 +/- 0.0794324 \t (limited)\n", "p1 = 22.7262 +/- 1.91125 \n", "p2 = 2.78007 +/- 0.877452 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 0.752979\n", "NDf = 7\n", "Edm = 5.34806e-06\n", "NCalls = 101\n", "p0 = 0.643852 +/- 0.170146 \t (limited)\n", "p1 = 21.01 +/- 4.32394 \n", "p2 = 5.36976 +/- 3.12567 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 12.9202\n", "Chi2 = 25.8404\n", "NDf = 23\n", "Edm = 1.28006e-08\n", "NCalls = 100\n", "p0 = 0.794265 +/- 0.0958159 \t (limited)\n", "p1 = 24.1679 +/- 3.59479 \n", "p2 = 4.72785 +/- 2.02176 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 2.5161\n", "NDf = 8\n", "Edm = 6.68524e-10\n", "NCalls = 61\n", "p0 = 0.628793 +/- 0.103749 \t (limited)\n", "p1 = 17.4724 +/- 2.12101 \n", "p2 = 2.59351 +/- 1.42575 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 17.7487\n", "Chi2 = 35.4974\n", "NDf = 24\n", "Edm = 4.85421e-08\n", "NCalls = 84\n", "p0 = 0.739601 +/- 0.111147 \t (limited)\n", "p1 = 23.7592 +/- 3.70376 \n", "p2 = 4.88756 +/- 1.97092 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 1.27043\n", "NDf = 8\n", "Edm = 6.23845e-06\n", "NCalls = 72\n", "p0 = 0.561127 +/- 0.138134 \t (limited)\n", "p1 = 18.1601 +/- 3.22498 \n", "p2 = 3.50338 +/- 1.84312 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 10.3177\n", "Chi2 = 20.6353\n", "NDf = 25\n", "Edm = 4.17867e-06\n", "NCalls = 74\n", "p0 = 0.803115 +/- 0.0973933 \t (limited)\n", "p1 = 24.5223 +/- 3.6204 \n", "p2 = 6.04914 +/- 2.15838 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 2.84972\n", "NDf = 8\n", "Edm = 1.39524e-06\n", "NCalls = 90\n", "p0 = 0.693419 +/- 0.285373 \t (limited)\n", "p1 = 24.7834 +/- 7.96871 \n", "p2 = 5.95505 +/- 3.88805 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 14.2763\n", "Chi2 = 28.5526\n", "NDf = 27\n", "Edm = 1.16259e-07\n", "NCalls = 71\n", "p0 = 0.81816 +/- 0.0960388 \t (limited)\n", "p1 = 27.5551 +/- 3.48223 \n", "p2 = 5.91219 +/- 1.84521 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 3.67016\n", "NDf = 5\n", "Edm = 1.80141e-06\n", "NCalls = 82\n", "p0 = 0.502645 +/- 0.176121 \t (limited)\n", "p1 = 20.1539 +/- 2.79976 \n", "p2 = 2.70329 +/- 1.65981 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 9.70246\n", "Chi2 = 19.4049\n", "NDf = 21\n", "Edm = 4.83549e-06\n", "NCalls = 68\n", "p0 = 0.862268 +/- 0.0732858 \t (limited)\n", "p1 = 22.9201 +/- 1.83831 \n", "p2 = 2.68862 +/- 0.964839 \n", "****************************************\n", " Invalid FitResult (status = 3 )\n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 0.936088\n", "NDf = 1\n", "Edm = 0.106873\n", "NCalls = 83\n", "p0 = 0.41816 +/- 0.669497 \t (limited)\n", "p1 = 18.3049 +/- 3.25491 \n", "p2 = 3.06774 +/- 40.715 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 3.77203\n", "NDf = 8\n", "Edm = 1.56279e-05\n", "NCalls = 84\n", "p0 = 0.642277 +/- 0.195607 \t (limited)\n", "p1 = 22.9222 +/- 4.97062 \n", "p2 = 5.05681 +/- 2.69739 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 16.8304\n", "Chi2 = 33.6608\n", "NDf = 24\n", "Edm = 5.58496e-07\n", "NCalls = 72\n", "p0 = 0.62357 +/- 0.0902387 \t (limited)\n", "p1 = 20.2778 +/- 2.24768 \n", "p2 = 2.93813 +/- 1.24803 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 0.47406\n", "NDf = 6\n", "Edm = 1.7448e-06\n", "NCalls = 88\n", "p0 = 0.635811 +/- 0.173792 \t (limited)\n", "p1 = 21.0478 +/- 3.90109 \n", "p2 = 3.91609 +/- 2.85302 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 10.7469\n", "Chi2 = 21.4938\n", "NDf = 21\n", "Edm = 2.82067e-07\n", "NCalls = 90\n", "p0 = 0.761502 +/- 0.0976432 \t (limited)\n", "p1 = 23.1813 +/- 2.58651 \n", "p2 = 3.04135 +/- 1.22248 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 2.39556\n", "NDf = 8\n", "Edm = 4.51695e-08\n", "NCalls = 81\n", "p0 = 0.564008 +/- 0.188334 \t (limited)\n", "p1 = 19.2305 +/- 5.64705 \n", "p2 = 4.00431 +/- 3.2294 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 16.5802\n", "Chi2 = 33.1604\n", "NDf = 27\n", "Edm = 2.38499e-09\n", "NCalls = 73\n", "p0 = 0.670525 +/- 0.116925 \t (limited)\n", "p1 = 23.7189 +/- 4.27475 \n", "p2 = 5.77511 +/- 2.44309 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 2.71864\n", "NDf = 3\n", "Edm = 2.02651e-07\n", "NCalls = 84\n", "p0 = 0.430955 +/- 0.103386 \t (limited)\n", "p1 = 15.3556 +/- 1.58614 \n", "p2 = 1.18646 +/- 1.05279 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 14.442\n", "Chi2 = 28.884\n", "NDf = 25\n", "Edm = 8.50826e-07\n", "NCalls = 110\n", "p0 = 0.931083 +/- 0.142888 \t (limited)\n", "p1 = 27.3321 +/- 5.16548 \n", "p2 = 5.43478 +/- 2.66681 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 3.91497\n", "NDf = 3\n", "Edm = 9.71484e-06\n", "NCalls = 172\n", "p0 = 0.999991 +/- 0.973989 \t (limited)\n", "p1 = 26.3679 +/- 4.16357 \n", "p2 = 10.2085 +/- 7.97289 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 13.7743\n", "Chi2 = 27.5485\n", "NDf = 26\n", "Edm = 4.38758e-07\n", "NCalls = 97\n", "p0 = 0.94497 +/- 0.0525858 \t (limited)\n", "p1 = 28.5699 +/- 2.31631 \n", "p2 = 4.52203 +/- 1.14841 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 3.69474\n", "NDf = 4\n", "Edm = 1.0286e-07\n", "NCalls = 124\n", "p0 = 0.566565 +/- 0.300627 \t (limited)\n", "p1 = 22.9958 +/- 7.83199 \n", "p2 = 6.59909 +/- 5.25595 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 12.1694\n", "Chi2 = 24.3388\n", "NDf = 24\n", "Edm = 1.18132e-06\n", "NCalls = 93\n", "p0 = 0.902749 +/- 0.0643623 \t (limited)\n", "p1 = 23.1006 +/- 1.54526 \n", "p2 = 3.46387 +/- 0.844882 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 1.01429\n", "NDf = 10\n", "Edm = 2.60252e-06\n", "NCalls = 94\n", "p0 = 0.610818 +/- 0.175926 \t (limited)\n", "p1 = 23.3058 +/- 5.77903 \n", "p2 = 5.09666 +/- 2.59915 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 7.14727\n", "Chi2 = 14.2945\n", "NDf = 19\n", "Edm = 1.36278e-07\n", "NCalls = 95\n", "p0 = 0.900785 +/- 0.229879 \t (limited)\n", "p1 = 28.3188 +/- 10.6921 \n", "p2 = 6.84773 +/- 5.79821 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 3.63097\n", "NDf = 7\n", "Edm = 2.77767e-07\n", "NCalls = 123\n", "p0 = 0.643028 +/- 0.225348 \t (limited)\n", "p1 = 21.4306 +/- 4.46994 \n", "p2 = 3.95327 +/- 2.04398 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 13.9433\n", "Chi2 = 27.8865\n", "NDf = 21\n", "Edm = 3.43188e-07\n", "NCalls = 63\n", "p0 = 0.683396 +/- 0.104286 \t (limited)\n", "p1 = 21.9331 +/- 2.48361 \n", "p2 = 3.36642 +/- 1.25116 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 1.7669\n", "NDf = 7\n", "Edm = 1.75691e-07\n", "NCalls = 98\n", "p0 = 0.73821 +/- 0.220641 \t (limited)\n", "p1 = 23.671 +/- 5.18271 \n", "p2 = 5.53649 +/- 2.72352 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 9.70486\n", "Chi2 = 19.4097\n", "NDf = 21\n", "Edm = 8.73736e-08\n", "NCalls = 108\n", "p0 = 1 +/- 0.129925 \t (limited)\n", "p1 = 30.7473 +/- 2.87702 \n", "p2 = 8.76457 +/- 1.97452 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 2.70826\n", "NDf = 8\n", "Edm = 5.00596e-07\n", "NCalls = 106\n", "p0 = 0.662591 +/- 0.352545 \t (limited)\n", "p1 = 24.2783 +/- 12.3384 \n", "p2 = 7.60846 +/- 6.66009 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 11.7754\n", "Chi2 = 23.5509\n", "NDf = 22\n", "Edm = 8.43281e-07\n", "NCalls = 108\n", "p0 = 0.865161 +/- 0.102972 \t (limited)\n", "p1 = 28.7874 +/- 2.92519 \n", "p2 = 5.66109 +/- 1.45554 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 9.74855\n", "NDf = 4\n", "Edm = 3.46464e-07\n", "NCalls = 124\n", "p0 = 0.606488 +/- 0.231103 \t (limited)\n", "p1 = 22.9147 +/- 5.11615 \n", "p2 = 4.88595 +/- 2.18611 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 14.5599\n", "Chi2 = 29.1197\n", "NDf = 26\n", "Edm = 5.37183e-08\n", "NCalls = 89\n", "p0 = 0.893177 +/- 0.0654554 \t (limited)\n", "p1 = 26.3026 +/- 1.88872 \n", "p2 = 3.66616 +/- 0.981145 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 3.10119\n", "NDf = 5\n", "Edm = 3.71723e-07\n", "NCalls = 234\n", "p0 = 0.999998 +/- 0.501156 \t (limited)\n", "p1 = 28.5973 +/- 4.11085 \n", "p2 = 8.15931 +/- 3.01069 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 17.676\n", "Chi2 = 35.352\n", "NDf = 25\n", "Edm = 4.23826e-07\n", "NCalls = 105\n", "p0 = 0.691753 +/- 0.102495 \t (limited)\n", "p1 = 22.3455 +/- 2.89145 \n", "p2 = 5.80589 +/- 1.86419 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 2.78024\n", "NDf = 7\n", "Edm = 1.5102e-06\n", "NCalls = 82\n", "p0 = 0.65308 +/- 0.206386 \t (limited)\n", "p1 = 23.0571 +/- 4.95138 \n", "p2 = 4.86356 +/- 2.47873 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 9.65754\n", "Chi2 = 19.3151\n", "NDf = 24\n", "Edm = 2.10929e-06\n", "NCalls = 104\n", "p0 = 1 +/- 0.130796 \t (limited)\n", "p1 = 31.0119 +/- 2.70979 \n", "p2 = 7.55857 +/- 1.6782 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 3.05436\n", "NDf = 5\n", "Edm = 1.66357e-05\n", "NCalls = 297\n", "p0 = 0.999547 +/- 0.735932 \t (limited)\n", "p1 = 47.3138 +/- 36.6961 \n", "p2 = 43.5419 +/- 61.6394 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 15.551\n", "Chi2 = 31.1019\n", "NDf = 24\n", "Edm = 3.03617e-07\n", "NCalls = 151\n", "p0 = 1 +/- 0.118805 \t (limited)\n", "p1 = 32.246 +/- 2.88351 \n", "p2 = 8.35109 +/- 1.76088 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 0.139804\n", "NDf = 2\n", "Edm = 1.06491e-06\n", "NCalls = 91\n", "p0 = 0.660157 +/- 0.258707 \t (limited)\n", "p1 = 19.2592 +/- 7.6936 \n", "p2 = 6.44669 +/- 4.78428 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 23.3673\n", "Chi2 = 46.7347\n", "NDf = 30\n", "Edm = 1.84329e-06\n", "NCalls = 94\n", "p0 = 0.74715 +/- 0.110543 \t (limited)\n", "p1 = 28.1656 +/- 7.86858 \n", "p2 = 8.19547 +/- 5.07211 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 2.55321\n", "NDf = 7\n", "Edm = 5.89829e-06\n", "NCalls = 89\n", "p0 = 0.643594 +/- 0.194962 \t (limited)\n", "p1 = 21.9893 +/- 5.24123 \n", "p2 = 5.90243 +/- 2.90213 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 13.6972\n", "Chi2 = 27.3944\n", "NDf = 22\n", "Edm = 5.81988e-08\n", "NCalls = 86\n", "p0 = 0.788092 +/- 0.0892723 \t (limited)\n", "p1 = 25.0025 +/- 2.64366 \n", "p2 = 4.90379 +/- 1.42976 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 2.36128\n", "NDf = 6\n", "Edm = 6.99629e-07\n", "NCalls = 88\n", "p0 = 0.560279 +/- 0.13121 \t (limited)\n", "p1 = 19.8526 +/- 3.03421 \n", "p2 = 2.47661 +/- 2.23047 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 16.199\n", "Chi2 = 32.398\n", "NDf = 22\n", "Edm = 2.10103e-06\n", "NCalls = 62\n", "p0 = 0.693765 +/- 0.0810703 \t (limited)\n", "p1 = 20.5818 +/- 1.69818 \n", "p2 = 3.08945 +/- 0.908341 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 1.33423\n", "NDf = 7\n", "Edm = 2.89896e-06\n", "NCalls = 102\n", "p0 = 0.564562 +/- 0.203287 \t (limited)\n", "p1 = 25.8413 +/- 9.79483 \n", "p2 = 8.35761 +/- 6.07925 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 12.5314\n", "Chi2 = 25.0627\n", "NDf = 26\n", "Edm = 1.61089e-07\n", "NCalls = 97\n", "p0 = 0.835727 +/- 0.083548 \t (limited)\n", "p1 = 30.5594 +/- 4.27567 \n", "p2 = 6.90591 +/- 2.2602 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 1.74047\n", "NDf = 7\n", "Edm = 7.04231e-07\n", "NCalls = 135\n", "p0 = 0.594764 +/- 0.187954 \t (limited)\n", "p1 = 25.3877 +/- 8.57244 \n", "p2 = 6.52667 +/- 4.25157 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 16.4051\n", "Chi2 = 32.8102\n", "NDf = 26\n", "Edm = 3.55957e-06\n", "NCalls = 78\n", "p0 = 0.710631 +/- 0.082398 \t (limited)\n", "p1 = 21.9892 +/- 2.24718 \n", "p2 = 2.94004 +/- 1.13284 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 0.310648\n", "NDf = 7\n", "Edm = 2.45348e-06\n", "NCalls = 169\n", "p0 = 0.539725 +/- 0.760718 \t (limited)\n", "p1 = 12.9065 +/- 72.7676 \n", "p2 = 22.8828 +/- 100.306 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 13.4706\n", "Chi2 = 26.9413\n", "NDf = 25\n", "Edm = 2.38853e-08\n", "NCalls = 153\n", "p0 = 0.812768 +/- 0.131071 \t (limited)\n", "p1 = 24.5484 +/- 4.86163 \n", "p2 = 5.70675 +/- 2.6421 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 2.52738\n", "NDf = 10\n", "Edm = 1.515e-08\n", "NCalls = 119\n", "p0 = 0.631088 +/- 0.21992 \t (limited)\n", "p1 = 26.5773 +/- 9.03434 \n", "p2 = 8.9826 +/- 4.83446 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 11.3338\n", "Chi2 = 22.6677\n", "NDf = 23\n", "Edm = 1.90162e-08\n", "NCalls = 71\n", "p0 = 0.745037 +/- 0.103424 \t (limited)\n", "p1 = 25.6472 +/- 3.41282 \n", "p2 = 7.04987 +/- 2.05878 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 2.51398\n", "NDf = 6\n", "Edm = 6.94471e-07\n", "NCalls = 84\n", "p0 = 0.484246 +/- 0.226988 \t (limited)\n", "p1 = 18.2432 +/- 5.68591 \n", "p2 = 3.93586 +/- 3.23302 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 10.8867\n", "Chi2 = 21.7734\n", "NDf = 23\n", "Edm = 3.78138e-08\n", "NCalls = 84\n", "p0 = 0.836954 +/- 0.108736 \t (limited)\n", "p1 = 24.375 +/- 3.34157 \n", "p2 = 6.16582 +/- 2.0379 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 0.946468\n", "NDf = 5\n", "Edm = 5.45857e-06\n", "NCalls = 74\n", "p0 = 0.584615 +/- 0.133896 \t (limited)\n", "p1 = 19.3728 +/- 3.81107 \n", "p2 = 4.20495 +/- 2.53748 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 12.0156\n", "Chi2 = 24.0312\n", "NDf = 23\n", "Edm = 1.40064e-06\n", "NCalls = 76\n", "p0 = 0.880019 +/- 0.0911312 \t (limited)\n", "p1 = 28.963 +/- 3.37154 \n", "p2 = 6.13976 +/- 1.63654 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 5.46724\n", "NDf = 8\n", "Edm = 1.55331e-06\n", "NCalls = 99\n", "p0 = 0.763745 +/- 0.233374 \t (limited)\n", "p1 = 26.3444 +/- 5.03855 \n", "p2 = 6.07972 +/- 2.87905 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 7.09585\n", "Chi2 = 14.1917\n", "NDf = 22\n", "Edm = 2.55685e-07\n", "NCalls = 100\n", "p0 = 1 +/- 0.137512 \t (limited)\n", "p1 = 31.2245 +/- 2.9769 \n", "p2 = 8.50085 +/- 1.97863 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 0.716283\n", "NDf = 7\n", "Edm = 1.89199e-06\n", "NCalls = 85\n", "p0 = 0.422768 +/- 0.179544 \t (limited)\n", "p1 = 17.3878 +/- 6.43713 \n", "p2 = 3.05887 +/- 4.49587 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 14.7289\n", "Chi2 = 29.4577\n", "NDf = 25\n", "Edm = 2.37579e-07\n", "NCalls = 81\n", "p0 = 0.74917 +/- 0.103618 \t (limited)\n", "p1 = 24.7084 +/- 3.46816 \n", "p2 = 4.65745 +/- 1.83133 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 6.40375\n", "NDf = 8\n", "Edm = 5.62321e-07\n", "NCalls = 104\n", "p0 = 0.507615 +/- 0.198193 \t (limited)\n", "p1 = 22.2924 +/- 5.65071 \n", "p2 = 4.8117 +/- 2.67857 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 14.2056\n", "Chi2 = 28.4112\n", "NDf = 25\n", "Edm = 4.55206e-07\n", "NCalls = 72\n", "p0 = 0.756894 +/- 0.0844728 \t (limited)\n", "p1 = 22.7347 +/- 2.14529 \n", "p2 = 3.73044 +/- 1.05431 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 3.18171\n", "NDf = 9\n", "Edm = 5.5664e-07\n", "NCalls = 83\n", "p0 = 0.694498 +/- 0.134559 \t (limited)\n", "p1 = 21.5399 +/- 2.86559 \n", "p2 = 3.46906 +/- 1.42287 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 14.71\n", "Chi2 = 29.4201\n", "NDf = 29\n", "Edm = 1.37273e-09\n", "NCalls = 63\n", "p0 = 0.708238 +/- 0.0850928 \t (limited)\n", "p1 = 21.4168 +/- 2.75995 \n", "p2 = 3.18025 +/- 1.28681 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 3.13918\n", "NDf = 6\n", "Edm = 4.11187e-06\n", "NCalls = 101\n", "p0 = 0.725096 +/- 0.145408 \t (limited)\n", "p1 = 25.2372 +/- 3.82188 \n", "p2 = 4.88315 +/- 1.80529 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 9.52542\n", "Chi2 = 19.0508\n", "NDf = 25\n", "Edm = 7.91504e-08\n", "NCalls = 72\n", "p0 = 0.879126 +/- 0.0778689 \t (limited)\n", "p1 = 29.5534 +/- 2.95061 \n", "p2 = 5.14138 +/- 1.41346 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 3.15774\n", "NDf = 6\n", "Edm = 2.32136e-06\n", "NCalls = 109\n", "p0 = 0.757722 +/- 0.770168 \t (limited)\n", "p1 = 26.0106 +/- 13.1967 \n", "p2 = 7.01687 +/- 5.25991 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 11.7233\n", "Chi2 = 23.4465\n", "NDf = 21\n", "Edm = 1.90261e-07\n", "NCalls = 89\n", "p0 = 1 +/- 0.073348 \t (limited)\n", "p1 = 30.5684 +/- 2.60029 \n", "p2 = 7.87047 +/- 1.70364 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 1.90008\n", "NDf = 6\n", "Edm = 9.90083e-06\n", "NCalls = 99\n", "p0 = 0.534999 +/- 0.270262 \t (limited)\n", "p1 = 24.0906 +/- 7.8809 \n", "p2 = 5.63935 +/- 3.24361 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 16.1477\n", "Chi2 = 32.2953\n", "NDf = 23\n", "Edm = 2.09463e-06\n", "NCalls = 73\n", "p0 = 0.83612 +/- 0.0936713 \t (limited)\n", "p1 = 25.2247 +/- 3.44688 \n", "p2 = 5.16788 +/- 1.95897 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 2.77855\n", "NDf = 8\n", "Edm = 8.13609e-08\n", "NCalls = 85\n", "p0 = 0.702887 +/- 0.199272 \t (limited)\n", "p1 = 20.9159 +/- 3.68349 \n", "p2 = 4.41504 +/- 2.04434 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 12.5946\n", "Chi2 = 25.1892\n", "NDf = 22\n", "Edm = 4.03445e-07\n", "NCalls = 53\n", "p0 = 0.686397 +/- 0.100168 \t (limited)\n", "p1 = 20.697 +/- 2.5238 \n", "p2 = 4.27225 +/- 1.40442 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 4.37217\n", "NDf = 9\n", "Edm = 2.72083e-07\n", "NCalls = 82\n", "p0 = 0.620959 +/- 0.15843 \t (limited)\n", "p1 = 21.0013 +/- 3.72087 \n", "p2 = 4.56895 +/- 2.37229 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 12.31\n", "Chi2 = 24.62\n", "NDf = 22\n", "Edm = 4.02211e-06\n", "NCalls = 62\n", "p0 = 0.738691 +/- 0.0999045 \t (limited)\n", "p1 = 21.6956 +/- 2.49688 \n", "p2 = 4.68089 +/- 1.47051 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 1.04047\n", "NDf = 4\n", "Edm = 5.64332e-06\n", "NCalls = 102\n", "p0 = 0.593644 +/- 0.174313 \t (limited)\n", "p1 = 21.0676 +/- 4.04434 \n", "p2 = 3.35047 +/- 2.25854 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 8.63265\n", "Chi2 = 17.2653\n", "NDf = 25\n", "Edm = 6.38533e-09\n", "NCalls = 101\n", "p0 = 0.893299 +/- 0.0630243 \t (limited)\n", "p1 = 24.6689 +/- 1.69694 \n", "p2 = 2.34492 +/- 0.886976 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 3.95319\n", "NDf = 10\n", "Edm = 1.32825e-07\n", "NCalls = 96\n", "p0 = 0.574996 +/- 0.152016 \t (limited)\n", "p1 = 20.2114 +/- 5.1893 \n", "p2 = 4.26461 +/- 2.57398 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 11.4512\n", "Chi2 = 22.9024\n", "NDf = 25\n", "Edm = 9.01028e-07\n", "NCalls = 89\n", "p0 = 0.894533 +/- 0.137242 \t (limited)\n", "p1 = 31.8196 +/- 6.42723 \n", "p2 = 9.36301 +/- 3.43315 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 1.2121\n", "NDf = 7\n", "Edm = 3.20149e-06\n", "NCalls = 91\n", "p0 = 0.498388 +/- 0.194573 \t (limited)\n", "p1 = 16.3792 +/- 10.7582 \n", "p2 = 4.04759 +/- 14.7395 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 14.2305\n", "Chi2 = 28.4611\n", "NDf = 27\n", "Edm = 7.65152e-11\n", "NCalls = 112\n", "p0 = 0.568628 +/- 0.0691435 \t (limited)\n", "p1 = 18.8276 +/- 0.571408 \n", "p2 = 0.0507941 +/- 7.32037e-18 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 1.15843\n", "NDf = 4\n", "Edm = 3.66535e-06\n", "NCalls = 105\n", "p0 = 0.570211 +/- 0.236898 \t (limited)\n", "p1 = 17.2578 +/- 15.7293 \n", "p2 = 9.7494 +/- 16.3952 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 20.2423\n", "Chi2 = 40.4846\n", "NDf = 29\n", "Edm = 1.41917e-06\n", "NCalls = 139\n", "p0 = 0.612238 +/- 0.0693685 \t (limited)\n", "p1 = 18.9111 +/- 0.641041 \n", "p2 = 0.0810507 +/- 4.47439 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 4.96865\n", "NDf = 8\n", "Edm = 1.88709e-08\n", "NCalls = 80\n", "p0 = 0.666513 +/- 0.140936 \t (limited)\n", "p1 = 21.2569 +/- 2.26557 \n", "p2 = 3.48413 +/- 2.03653 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 16.7227\n", "Chi2 = 33.4455\n", "NDf = 22\n", "Edm = 1.37912e-07\n", "NCalls = 60\n", "p0 = 0.671333 +/- 0.0830986 \t (limited)\n", "p1 = 20.0899 +/- 1.80751 \n", "p2 = 3.18801 +/- 0.983175 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 2.88079\n", "NDf = 4\n", "Edm = 8.70733e-06\n", "NCalls = 96\n", "p0 = 0.803697 +/- 0.275939 \t (limited)\n", "p1 = 22.1124 +/- 5.26388 \n", "p2 = 5.42701 +/- 3.06391 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 11.2809\n", "Chi2 = 22.5619\n", "NDf = 22\n", "Edm = 1.17703e-06\n", "NCalls = 103\n", "p0 = 1 +/- 0.0550354 \t (limited)\n", "p1 = 28.1892 +/- 1.93212 \n", "p2 = 5.42012 +/- 1.15081 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 3.96392\n", "NDf = 5\n", "Edm = 7.79729e-08\n", "NCalls = 100\n", "p0 = 0.694889 +/- 0.235894 \t (limited)\n", "p1 = 23.9686 +/- 5.40538 \n", "p2 = 5.07476 +/- 3.36335 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 10.6942\n", "Chi2 = 21.3884\n", "NDf = 24\n", "Edm = 4.94209e-06\n", "NCalls = 90\n", "p0 = 0.886412 +/- 0.0741325 \t (limited)\n", "p1 = 26.3919 +/- 2.05984 \n", "p2 = 3.75207 +/- 0.962698 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 0.852791\n", "NDf = 5\n", "Edm = 8.88345e-09\n", "NCalls = 84\n", "p0 = 0.310283 +/- 0.164695 \t (limited)\n", "p1 = 15.0842 +/- 5.86864 \n", "p2 = 2.43104 +/- 3.51103 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 16.7969\n", "Chi2 = 33.5939\n", "NDf = 22\n", "Edm = 1.06214e-06\n", "NCalls = 90\n", "p0 = 0.858009 +/- 0.0896248 \t (limited)\n", "p1 = 26.7658 +/- 2.7434 \n", "p2 = 5.91245 +/- 1.48212 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 1.32394\n", "NDf = 6\n", "Edm = 6.16879e-07\n", "NCalls = 86\n", "p0 = 0.672106 +/- 0.201669 \t (limited)\n", "p1 = 23.2133 +/- 5.10201 \n", "p2 = 5.04859 +/- 2.50327 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 8.61404\n", "Chi2 = 17.2281\n", "NDf = 24\n", "Edm = 6.71438e-08\n", "NCalls = 79\n", "p0 = 0.940369 +/- 0.0595691 \t (limited)\n", "p1 = 28.7077 +/- 2.62934 \n", "p2 = 5.3511 +/- 1.39597 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 2.48447\n", "NDf = 9\n", "Edm = 3.99523e-06\n", "NCalls = 92\n", "p0 = 0.469573 +/- 0.142912 \t (limited)\n", "p1 = 18.7922 +/- 3.85884 \n", "p2 = 2.98044 +/- 2.31127 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 10.8641\n", "Chi2 = 21.7281\n", "NDf = 22\n", "Edm = 2.68775e-06\n", "NCalls = 76\n", "p0 = 0.705796 +/- 0.112524 \t (limited)\n", "p1 = 22.694 +/- 3.20217 \n", "p2 = 4.31371 +/- 1.76006 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 4.05841\n", "NDf = 7\n", "Edm = 3.35768e-06\n", "NCalls = 93\n", "p0 = 0.667164 +/- 0.148637 \t (limited)\n", "p1 = 21.9757 +/- 2.55643 \n", "p2 = 3.34362 +/- 2.015 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 12.3026\n", "Chi2 = 24.6053\n", "NDf = 24\n", "Edm = 6.94447e-09\n", "NCalls = 70\n", "p0 = 0.771226 +/- 0.0842579 \t (limited)\n", "p1 = 22.8924 +/- 2.39483 \n", "p2 = 4.3482 +/- 1.40708 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 3.14114\n", "NDf = 6\n", "Edm = 1.74933e-08\n", "NCalls = 114\n", "p0 = 0.532234 +/- 0.174899 \t (limited)\n", "p1 = 20.0178 +/- 3.71281 \n", "p2 = 3.05935 +/- 2.11129 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 13.6754\n", "Chi2 = 27.3508\n", "NDf = 25\n", "Edm = 3.12686e-06\n", "NCalls = 86\n", "p0 = 0.785653 +/- 0.0837745 \t (limited)\n", "p1 = 23.5391 +/- 1.85462 \n", "p2 = 2.84004 +/- 0.863279 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 2.88405\n", "NDf = 5\n", "Edm = 1.25252e-05\n", "NCalls = 77\n", "p0 = 0.698497 +/- 0.183511 \t (limited)\n", "p1 = 21.6581 +/- 4.18422 \n", "p2 = 4.04207 +/- 1.56968 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 16.2927\n", "Chi2 = 32.5855\n", "NDf = 26\n", "Edm = 9.81736e-09\n", "NCalls = 73\n", "p0 = 0.84928 +/- 0.0880232 \t (limited)\n", "p1 = 26.1083 +/- 2.99956 \n", "p2 = 5.04526 +/- 1.47387 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 0.624725\n", "NDf = 4\n", "Edm = 6.43892e-06\n", "NCalls = 71\n", "p0 = 0.380579 +/- 0.134096 \t (limited)\n", "p1 = 13.9323 +/- 2.56334 \n", "p2 = 1.96554 +/- 1.64705 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 16.2451\n", "Chi2 = 32.4901\n", "NDf = 24\n", "Edm = 3.46148e-07\n", "NCalls = 111\n", "p0 = 1 +/- 0.13976 \t (limited)\n", "p1 = 33.5941 +/- 3.2286 \n", "p2 = 9.3535 +/- 2.03064 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 1.67529\n", "NDf = 7\n", "Edm = 4.64318e-08\n", "NCalls = 80\n", "p0 = 0.616278 +/- 0.1576 \t (limited)\n", "p1 = 18.812 +/- 2.91455 \n", "p2 = 3.21096 +/- 1.42859 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 18.0638\n", "Chi2 = 36.1276\n", "NDf = 26\n", "Edm = 1.22808e-06\n", "NCalls = 66\n", "p0 = 0.675732 +/- 0.0875617 \t (limited)\n", "p1 = 20.9227 +/- 2.1898 \n", "p2 = 3.60151 +/- 1.27833 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 0.374196\n", "NDf = 2\n", "Edm = 8.80185e-06\n", "NCalls = 80\n", "p0 = 0.685904 +/- 0.174442 \t (limited)\n", "p1 = 19.998 +/- 1.62445 \n", "p2 = 1.71179 +/- 1.75684 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 17.4763\n", "Chi2 = 34.9527\n", "NDf = 25\n", "Edm = 1.18689e-07\n", "NCalls = 79\n", "p0 = 0.703138 +/- 0.0910538 \t (limited)\n", "p1 = 20.3182 +/- 1.51335 \n", "p2 = 1.05166 +/- 1.00855 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 2.46931\n", "NDf = 7\n", "Edm = 3.32581e-07\n", "NCalls = 93\n", "p0 = 0.640812 +/- 0.135875 \t (limited)\n", "p1 = 22.3459 +/- 3.96589 \n", "p2 = 4.39764 +/- 2.51185 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 19.573\n", "Chi2 = 39.1461\n", "NDf = 26\n", "Edm = 1.8708e-07\n", "NCalls = 79\n", "p0 = 0.695688 +/- 0.0960438 \t (limited)\n", "p1 = 25.9337 +/- 2.91836 \n", "p2 = 3.8781 +/- 1.39167 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 2.49631\n", "NDf = 3\n", "Edm = 6.6234e-06\n", "NCalls = 102\n", "p0 = 0.611704 +/- 0.35871 \t (limited)\n", "p1 = 20.7057 +/- 8.42105 \n", "p2 = 5.08053 +/- 3.9281 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 19.5603\n", "Chi2 = 39.1206\n", "NDf = 22\n", "Edm = 2.27081e-08\n", "NCalls = 74\n", "p0 = 0.700633 +/- 0.0849755 \t (limited)\n", "p1 = 19.1288 +/- 1.84326 \n", "p2 = 3.20784 +/- 1.0308 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 0.418946\n", "NDf = 2\n", "Edm = 7.44661e-07\n", "NCalls = 53\n", "p0 = 0.60361 +/- 0.106733 \t (limited)\n", "p1 = 15.7806 +/- 2 \n", "p2 = 0.105266 +/- 2 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 19.3993\n", "Chi2 = 38.7986\n", "NDf = 23\n", "Edm = 9.22231e-06\n", "NCalls = 245\n", "p0 = 0.644389 +/- 0.0711313 \t (limited)\n", "p1 = 18.1356 +/- 0.714145 \n", "p2 = 0.0157025 +/- 0.682189 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 3.54456\n", "NDf = 5\n", "Edm = 2.33566e-06\n", "NCalls = 85\n", "p0 = 0.527397 +/- 0.217626 \t (limited)\n", "p1 = 17.9928 +/- 6.38912 \n", "p2 = 5.4498 +/- 3.84602 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 12.2163\n", "Chi2 = 24.4327\n", "NDf = 20\n", "Edm = 6.3676e-06\n", "NCalls = 80\n", "p0 = 0.859519 +/- 0.0906147 \t (limited)\n", "p1 = 25.943 +/- 3.43909 \n", "p2 = 6.71027 +/- 2.02884 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 3.23312\n", "NDf = 8\n", "Edm = 9.46056e-07\n", "NCalls = 70\n", "p0 = 0.515942 +/- 0.0980029 \t (limited)\n", "p1 = 17.2963 +/- 1.44299 \n", "p2 = 1.9009 +/- 1.23849 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 11.3898\n", "Chi2 = 22.7795\n", "NDf = 20\n", "Edm = 6.32872e-09\n", "NCalls = 72\n", "p0 = 0.633638 +/- 0.0838664 \t (limited)\n", "p1 = 18.524 +/- 1.68856 \n", "p2 = 2.16871 +/- 1.02224 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 5.28211\n", "NDf = 11\n", "Edm = 1.87776e-06\n", "NCalls = 71\n", "p0 = 0.612568 +/- 0.109907 \t (limited)\n", "p1 = 19.3198 +/- 2.18618 \n", "p2 = 3.37801 +/- 1.31028 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 9.91351\n", "Chi2 = 19.827\n", "NDf = 22\n", "Edm = 2.00204e-07\n", "NCalls = 62\n", "p0 = 0.636118 +/- 0.097181 \t (limited)\n", "p1 = 19.5971 +/- 2.68815 \n", "p2 = 3.88968 +/- 1.7373 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 2.51155\n", "NDf = 7\n", "Edm = 2.60724e-06\n", "NCalls = 74\n", "p0 = 0.647126 +/- 0.147452 \t (limited)\n", "p1 = 20.8474 +/- 2.85119 \n", "p2 = 3.27715 +/- 1.66359 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 6.94327\n", "Chi2 = 13.8865\n", "NDf = 22\n", "Edm = 4.86511e-06\n", "NCalls = 73\n", "p0 = 0.901211 +/- 0.0749805 \t (limited)\n", "p1 = 24.7879 +/- 2.21116 \n", "p2 = 3.76426 +/- 1.07742 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 2.33728\n", "NDf = 6\n", "Edm = 7.37813e-07\n", "NCalls = 109\n", "p0 = 0.624689 +/- 0.141987 \t (limited)\n", "p1 = 19.3681 +/- 3.57572 \n", "p2 = 3.40927 +/- 2.17685 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 14.2977\n", "Chi2 = 28.5955\n", "NDf = 24\n", "Edm = 3.94253e-10\n", "NCalls = 81\n", "p0 = 0.832764 +/- 0.0867739 \t (limited)\n", "p1 = 24.0485 +/- 2.82721 \n", "p2 = 4.07159 +/- 1.4461 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 3.75642\n", "NDf = 8\n", "Edm = 5.68096e-06\n", "NCalls = 120\n", "p0 = 0.387014 +/- 0.0702296 \t (limited)\n", "p1 = 11.5811 +/- 8.22684 \n", "p2 = 0.313997 +/- 4.44032 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 18.3537\n", "Chi2 = 36.7074\n", "NDf = 24\n", "Edm = 3.14219e-08\n", "NCalls = 122\n", "p0 = 0.766056 +/- 0.130733 \t (limited)\n", "p1 = 26.3265 +/- 5.29894 \n", "p2 = 8.25159 +/- 3.37803 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 1.68054\n", "NDf = 6\n", "Edm = 5.00607e-06\n", "NCalls = 79\n", "p0 = 0.492763 +/- 0.200749 \t (limited)\n", "p1 = 18.8758 +/- 5.209 \n", "p2 = 3.88829 +/- 2.85902 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 12.5237\n", "Chi2 = 25.0474\n", "NDf = 21\n", "Edm = 6.19402e-09\n", "NCalls = 83\n", "p0 = 0.82576 +/- 0.0905004 \t (limited)\n", "p1 = 24.7459 +/- 2.64975 \n", "p2 = 5.53589 +/- 1.50229 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 1.37515\n", "NDf = 4\n", "Edm = 3.40318e-06\n", "NCalls = 106\n", "p0 = 0.759679 +/- 0.234791 \t (limited)\n", "p1 = 24.2145 +/- 5.44887 \n", "p2 = 4.93415 +/- 2.65778 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 9.98776\n", "Chi2 = 19.9755\n", "NDf = 23\n", "Edm = 6.74125e-07\n", "NCalls = 88\n", "p0 = 0.922627 +/- 0.053996 \t (limited)\n", "p1 = 25.2002 +/- 1.76064 \n", "p2 = 3.37452 +/- 0.879146 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 2.35692\n", "NDf = 4\n", "Edm = 1.08773e-07\n", "NCalls = 81\n", "p0 = 0.373315 +/- 0.215845 \t (limited)\n", "p1 = 12.4572 +/- 7.57085 \n", "p2 = 4.00427 +/- 8.97995 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 15.2686\n", "Chi2 = 30.5373\n", "NDf = 24\n", "Edm = 4.23993e-08\n", "NCalls = 154\n", "p0 = 1 +/- 0.165147 \t (limited)\n", "p1 = 30.0737 +/- 2.81991 \n", "p2 = 7.24184 +/- 1.65777 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 1.93137\n", "NDf = 7\n", "Edm = 4.79943e-09\n", "NCalls = 113\n", "p0 = 0.632872 +/- 0.231351 \t (limited)\n", "p1 = 27.1811 +/- 8.5188 \n", "p2 = 6.83092 +/- 3.30519 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 11.4871\n", "Chi2 = 22.9743\n", "NDf = 24\n", "Edm = 3.13018e-07\n", "NCalls = 73\n", "p0 = 0.886287 +/- 0.0822954 \t (limited)\n", "p1 = 26.4994 +/- 2.79725 \n", "p2 = 5.06979 +/- 1.36513 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 3.28016\n", "NDf = 6\n", "Edm = 6.04893e-06\n", "NCalls = 99\n", "p0 = 0.696957 +/- 0.193092 \t (limited)\n", "p1 = 22.9427 +/- 3.44124 \n", "p2 = 3.63739 +/- 2.36972 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 14.002\n", "Chi2 = 28.0041\n", "NDf = 25\n", "Edm = 3.57658e-07\n", "NCalls = 62\n", "p0 = 0.820512 +/- 0.0837157 \t (limited)\n", "p1 = 24.7539 +/- 2.2824 \n", "p2 = 4.0766 +/- 1.13837 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 2.56902\n", "NDf = 7\n", "Edm = 5.20615e-07\n", "NCalls = 94\n", "p0 = 0.68832 +/- 0.19603 \t (limited)\n", "p1 = 24.2988 +/- 4.41739 \n", "p2 = 4.93562 +/- 3.1133 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 10.1565\n", "Chi2 = 20.313\n", "NDf = 23\n", "Edm = 4.51533e-06\n", "NCalls = 73\n", "p0 = 0.83692 +/- 0.100458 \t (limited)\n", "p1 = 26.5351 +/- 2.71955 \n", "p2 = 4.32699 +/- 1.33756 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 3.04459\n", "NDf = 7\n", "Edm = 7.4353e-07\n", "NCalls = 83\n", "p0 = 0.659997 +/- 0.249426 \t (limited)\n", "p1 = 22.9747 +/- 5.95617 \n", "p2 = 5.36022 +/- 2.66473 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 11.1635\n", "Chi2 = 22.327\n", "NDf = 21\n", "Edm = 1.52311e-08\n", "NCalls = 72\n", "p0 = 0.797119 +/- 0.109369 \t (limited)\n", "p1 = 23.7528 +/- 3.02536 \n", "p2 = 5.51075 +/- 1.78264 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 6.27397\n", "NDf = 8\n", "Edm = 5.29961e-07\n", "NCalls = 81\n", "p0 = 0.662433 +/- 0.15714 \t (limited)\n", "p1 = 21.5909 +/- 3.29623 \n", "p2 = 4.09179 +/- 1.74745 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 14.006\n", "Chi2 = 28.012\n", "NDf = 23\n", "Edm = 1.73409e-07\n", "NCalls = 73\n", "p0 = 0.74949 +/- 0.087125 \t (limited)\n", "p1 = 22.4394 +/- 2.25986 \n", "p2 = 3.18063 +/- 1.06939 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 1.94311\n", "NDf = 9\n", "Edm = 1.89133e-07\n", "NCalls = 106\n", "p0 = 0.575157 +/- 0.298322 \t (limited)\n", "p1 = 25.6868 +/- 11.7079 \n", "p2 = 7.54114 +/- 5.74853 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 7.95143\n", "Chi2 = 15.9029\n", "NDf = 22\n", "Edm = 3.34787e-06\n", "NCalls = 100\n", "p0 = 1 +/- 0.0861316 \t (limited)\n", "p1 = 33.4794 +/- 3.24851 \n", "p2 = 9.96621 +/- 2.27222 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 2.75509\n", "NDf = 9\n", "Edm = 9.98946e-08\n", "NCalls = 69\n", "p0 = 0.593401 +/- 0.101469 \t (limited)\n", "p1 = 18.8032 +/- 2.57519 \n", "p2 = 2.53183 +/- 1.7654 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 8.82691\n", "Chi2 = 17.6538\n", "NDf = 19\n", "Edm = 8.94292e-07\n", "NCalls = 66\n", "p0 = 0.707869 +/- 0.0806949 \t (limited)\n", "p1 = 20.75 +/- 1.52055 \n", "p2 = 2.05399 +/- 0.709689 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 2.54952\n", "NDf = 6\n", "Edm = 4.57669e-08\n", "NCalls = 96\n", "p0 = 0.660242 +/- 0.259374 \t (limited)\n", "p1 = 25.0771 +/- 5.98352 \n", "p2 = 5.39543 +/- 2.98249 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 13.093\n", "Chi2 = 26.1861\n", "NDf = 25\n", "Edm = 7.24454e-07\n", "NCalls = 63\n", "p0 = 0.822936 +/- 0.0852676 \t (limited)\n", "p1 = 26.505 +/- 2.8571 \n", "p2 = 4.82335 +/- 1.47173 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 4.62283\n", "NDf = 10\n", "Edm = 1.36401e-08\n", "NCalls = 89\n", "p0 = 0.592713 +/- 0.183874 \t (limited)\n", "p1 = 20.9349 +/- 5.53213 \n", "p2 = 6.39618 +/- 4.21336 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 10.8024\n", "Chi2 = 21.6048\n", "NDf = 22\n", "Edm = 3.61653e-08\n", "NCalls = 73\n", "p0 = 0.718975 +/- 0.133739 \t (limited)\n", "p1 = 24.3254 +/- 4.59911 \n", "p2 = 7.31889 +/- 3.07985 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 5.3127\n", "NDf = 6\n", "Edm = 2.53153e-06\n", "NCalls = 119\n", "p0 = 0.534666 +/- 0.0847463 \t (limited)\n", "p1 = 12.0856 +/- 13.9766 \n", "p2 = 0.685012 +/- 8.82276 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 16.3818\n", "Chi2 = 32.7636\n", "NDf = 22\n", "Edm = 3.8063e-06\n", "NCalls = 100\n", "p0 = 0.789732 +/- 0.103694 \t (limited)\n", "p1 = 24.4115 +/- 2.95524 \n", "p2 = 4.45708 +/- 1.48761 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 0.968364\n", "NDf = 7\n", "Edm = 1.15106e-06\n", "NCalls = 74\n", "p0 = 0.669971 +/- 0.134042 \t (limited)\n", "p1 = 19.9248 +/- 3.66311 \n", "p2 = 4.7186 +/- 3.31403 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 15.3613\n", "Chi2 = 30.7226\n", "NDf = 24\n", "Edm = 5.15702e-06\n", "NCalls = 88\n", "p0 = 0.780567 +/- 0.0894369 \t (limited)\n", "p1 = 25.9084 +/- 2.75967 \n", "p2 = 4.59539 +/- 1.41025 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 2.88345\n", "NDf = 10\n", "Edm = 2.24225e-06\n", "NCalls = 102\n", "p0 = 0.629842 +/- 0.174064 \t (limited)\n", "p1 = 21.9766 +/- 6.0423 \n", "p2 = 6.83408 +/- 3.9968 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 9.27297\n", "Chi2 = 18.5459\n", "NDf = 20\n", "Edm = 1.87394e-06\n", "NCalls = 74\n", "p0 = 0.89538 +/- 0.174255 \t (limited)\n", "p1 = 32.6871 +/- 8.87808 \n", "p2 = 12.1149 +/- 5.18796 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 1.80215\n", "NDf = 6\n", "Edm = 2.68316e-07\n", "NCalls = 87\n", "p0 = 0.495591 +/- 0.211145 \t (limited)\n", "p1 = 19.645 +/- 5.2052 \n", "p2 = 3.59211 +/- 2.51889 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 11.3263\n", "Chi2 = 22.6526\n", "NDf = 23\n", "Edm = 6.21131e-07\n", "NCalls = 85\n", "p0 = 0.921303 +/- 0.0720411 \t (limited)\n", "p1 = 26.5793 +/- 2.38289 \n", "p2 = 5.11338 +/- 1.20119 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 3.79465\n", "NDf = 6\n", "Edm = 4.35871e-08\n", "NCalls = 69\n", "p0 = 0.611171 +/- 0.128277 \t (limited)\n", "p1 = 17.467 +/- 2.93532 \n", "p2 = 3.02848 +/- 1.64432 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 16.2411\n", "Chi2 = 32.4822\n", "NDf = 24\n", "Edm = 5.63919e-07\n", "NCalls = 84\n", "p0 = 0.887021 +/- 0.149113 \t (limited)\n", "p1 = 28.9301 +/- 6.59512 \n", "p2 = 7.99154 +/- 3.69369 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 2.748\n", "NDf = 7\n", "Edm = 5.78153e-09\n", "NCalls = 140\n", "p0 = 1 +/- 0.733966 \t (limited)\n", "p1 = 28.7334 +/- 2.38732 \n", "p2 = 6.33261 +/- 1.85434 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 14.7257\n", "Chi2 = 29.4515\n", "NDf = 25\n", "Edm = 6.17545e-08\n", "NCalls = 98\n", "p0 = 0.793408 +/- 0.149393 \t (limited)\n", "p1 = 25.3941 +/- 5.1633 \n", "p2 = 6.20929 +/- 3.20508 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 4.29305\n", "NDf = 9\n", "Edm = 3.06585e-06\n", "NCalls = 82\n", "p0 = 0.687721 +/- 0.16847 \t (limited)\n", "p1 = 22.2779 +/- 3.60769 \n", "p2 = 4.13806 +/- 1.75728 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 6.51615\n", "Chi2 = 13.0323\n", "NDf = 21\n", "Edm = 2.56493e-06\n", "NCalls = 80\n", "p0 = 0.929011 +/- 0.128044 \t (limited)\n", "p1 = 26.8551 +/- 4.60022 \n", "p2 = 6.16639 +/- 2.6095 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 1.74986\n", "NDf = 2\n", "Edm = 2.69957e-06\n", "NCalls = 256\n", "p0 = 0.999951 +/- 0.988093 \t (limited)\n", "p1 = 33.6961 +/- 9.16268 \n", "p2 = 13.7955 +/- 8.24703 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 21.3855\n", "Chi2 = 42.7711\n", "NDf = 26\n", "Edm = 4.23083e-07\n", "NCalls = 149\n", "p0 = 0.808068 +/- 0.0853178 \t (limited)\n", "p1 = 26.8921 +/- 2.88537 \n", "p2 = 6.0004 +/- 1.56206 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 4.72905\n", "NDf = 5\n", "Edm = 2.89445e-07\n", "NCalls = 67\n", "p0 = 0.417412 +/- 0.107227 \t (limited)\n", "p1 = 14.5796 +/- 2.95636 \n", "p2 = 2.22066 +/- 2.12966 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 20.7746\n", "Chi2 = 41.5492\n", "NDf = 25\n", "Edm = 4.86099e-06\n", "NCalls = 86\n", "p0 = 0.780136 +/- 0.0893488 \t (limited)\n", "p1 = 24.5549 +/- 2.10448 \n", "p2 = 4.36318 +/- 1.12085 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 2.45643\n", "NDf = 6\n", "Edm = 2.63752e-06\n", "NCalls = 99\n", "p0 = 0.554001 +/- 0.226655 \t (limited)\n", "p1 = 23.9121 +/- 6.69531 \n", "p2 = 6.04267 +/- 3.44117 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 13.3134\n", "Chi2 = 26.6268\n", "NDf = 25\n", "Edm = 5.83165e-09\n", "NCalls = 81\n", "p0 = 0.859757 +/- 0.0748063 \t (limited)\n", "p1 = 25.7851 +/- 2.03938 \n", "p2 = 4.58085 +/- 1.11229 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 2.91637\n", "NDf = 4\n", "Edm = 1.377e-05\n", "NCalls = 84\n", "p0 = 0.769946 +/- 0.185015 \t (limited)\n", "p1 = 20.5779 +/- 2.52006 \n", "p2 = 3.50084 +/- 1.35396 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 18.6117\n", "Chi2 = 37.2235\n", "NDf = 24\n", "Edm = 1.80945e-07\n", "NCalls = 72\n", "p0 = 0.788394 +/- 0.104211 \t (limited)\n", "p1 = 24.2646 +/- 3.57577 \n", "p2 = 5.17213 +/- 2.03607 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "Chi2 = 1.21412\n", "NDf = 6\n", "Edm = 1.16053e-08\n", "NCalls = 141\n", "p0 = 0.632519 +/- 0.106789 \t (limited)\n", "p1 = 17.4834 +/- 1.66297 \n", "p2 = 1.44082 +/- 1.10945 \n", "****************************************\n", "Minimizer is Minuit2 / Migrad\n", "MinFCN = 16.0874\n", "Chi2 = 32.1748\n", "NDf = 25\n", "Edm = 1.23857e-08\n", "NCalls = 73\n", "p0 = 0.683032 +/- 0.0852759 \t (limited)\n", "p1 = 20.5028 +/- 2.71991 \n", "p2 = 2.79578 +/- 1.30889 \n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Warning in : Abnormal termination of minimization.\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n", "Warning in : Sum of squares of weights structure already created\n", "Info in : Setting limits for parameter p0 to [0.000000,1.000000]\n", "Info in : Successful Result from Binomial Efficiency fitter of function fM2Fit\n" ] } ], "source": [ "for (int iloop = 0; iloop < nloop; ++iloop) {\n", "\n", " // generate pseudo-experiments\n", " hM2D->Reset();\n", " hM2N->Reset();\n", " hM2D->FillRandom(fM2D->GetName(), nevtsD);\n", " hM2N->FillRandom(fM2N->GetName(), nevtsN);\n", " hM2D->Add(hM2N);\n", "\n", " // construct the \"efficiency\" histogram\n", " hM2N->Sumw2();\n", " hM2E->Divide(hM2N, hM2D, 1, 1, \"b\");\n", "\n", " // Fit twice, using the same fit function.\n", " // In the first (standard) fit, initialize to (arbitrary) values.\n", " // In the second fit, use the results from the first fit (this\n", " // makes it easier for the fit -- but the purpose here is not to\n", " // show how easy or difficult it is to obtain results, but whether\n", " // the CORRECT results are obtained or not!).\n", "\n", " fM2Fit->SetParameter(0, 0.5);\n", " fM2Fit->SetParameter(1, 15.0);\n", " fM2Fit->SetParameter(2, 2.0);\n", " fM2Fit->SetParError(0, 0.1);\n", " fM2Fit->SetParError(1, 1.0);\n", " fM2Fit->SetParError(2, 0.2);\n", " TH1 * hf = fM2Fit->GetHistogram();\n", " // std::cout << \"Function values \" << std::endl;\n", " // for (int i = 1; i <= hf->GetNbinsX(); ++i)\n", " // std::cout << hf->GetBinContent(i) << \" \";\n", " // std::cout << std::endl;\n", "\n", " TCanvas* cEvt;\n", " if (plot) {\n", " cEvt = new TCanvas(Form(\"cEnv%d\",iloop),\n", " Form(\"plots for experiment %d\", iloop),\n", " 1000, 600);\n", " cEvt->Divide(1,2);\n", " cEvt->cd(1);\n", " hM2D->DrawCopy(\"HIST\");\n", " hM2N->SetLineColor(kRed);\n", " hM2N->DrawCopy(\"HIST SAME\");\n", " cEvt->cd(2);\n", " }\n", " for (int fit = 0; fit < 2; ++fit) {\n", " int status = 0;\n", " switch (fit) {\n", " case 0:\n", " {\n", " // TVirtualPad * pad = gPad;\n", " // new TCanvas();\n", " // fM2Fit->Draw();\n", " // gPad = pad;\n", " TString optFit = \"RN\";\n", " if (debug) optFit += TString(\"SV\");\n", " TFitResultPtr res = hM2E->Fit(fM2Fit, optFit);\n", " if (plot) {\n", " hM2E->DrawCopy(\"E\");\n", " fM2Fit->SetLineColor(kBlue);\n", " fM2Fit->DrawCopy(\"SAME\");\n", " }\n", " if (debug) res->Print();\n", " status = res;\n", " break;\n", " }\n", " case 1:\n", " {\n", " // if (fM2Fit2) delete fM2Fit2;\n", " // fM2Fit2 = dynamic_cast(fM2Fit->Clone(\"fM2Fit2\"));\n", " fM2Fit2 = fM2Fit; // do not clone/copy the function\n", " if (fM2Fit2->GetParameter(0) >= 1.0)\n", " fM2Fit2->SetParameter(0, 0.95);\n", " fM2Fit2->SetParLimits(0, 0.0, 1.0);\n", "\n", " // TVirtualPad * pad = gPad;\n", " // new TCanvas();\n", " // fM2Fit2->Draw();\n", " // gPad = pad;\n", "\n", " TBinomialEfficiencyFitter bef(hM2N, hM2D);\n", " TString optFit = \"RI S\";\n", " if (debug) optFit += TString(\"V\");\n", " TFitResultPtr res = bef.Fit(fM2Fit2,optFit);\n", " status = res;\n", " if (status !=0) {\n", " std::cerr << \"Error performing binomial efficiency fit, result = \"\n", " << status << std::endl;\n", " res->Print();\n", " continue;\n", " }\n", " if (plot) {\n", " fM2Fit2->SetLineColor(kRed);\n", " fM2Fit2->DrawCopy(\"SAME\");\n", " \n", " bool confint = (status == 0);\n", " if (confint) {\n", " // compute confidence interval on fitted function\n", " auto htemp = fM2Fit2->GetHistogram();\n", " ROOT::Fit::BinData xdata;\n", " ROOT::Fit::FillData(xdata, fM2Fit2->GetHistogram() );\n", " TGraphErrors gr(fM2Fit2->GetHistogram() );\n", " res->GetConfidenceIntervals(xdata, gr.GetEY(), 0.68, false);\n", " gr.SetFillColor(6);\n", " gr.SetFillStyle(3005);\n", " gr.DrawClone(\"4 same\");\n", " }\n", " }\n", " if (debug) {\n", " res->Print();\n", " }\n", " }\n", " }\n", "\n", " if (status != 0) break;\n", "\n", " double fnorm = fM2Fit->GetParameter(0);\n", " double enorm = fM2Fit->GetParError(0);\n", " double fthreshold = fM2Fit->GetParameter(1);\n", " double ethreshold = fM2Fit->GetParError(1);\n", " double fwidth = fM2Fit->GetParameter(2);\n", " double ewidth = fM2Fit->GetParError(2);\n", " if (fit == 1) {\n", " fnorm = fM2Fit2->GetParameter(0);\n", " enorm = fM2Fit2->GetParError(0);\n", " fthreshold = fM2Fit2->GetParameter(1);\n", " ethreshold = fM2Fit2->GetParError(1);\n", " fwidth = fM2Fit2->GetParameter(2);\n", " ewidth = fM2Fit2->GetParError(2);\n", " hChisquared->Fill(fM2Fit2->GetProb());\n", " }\n", "\n", " TH1D* h = dynamic_cast(hbiasNorm[fit]);\n", " h->Fill((fnorm-normalization)/enorm);\n", " h = dynamic_cast(hbiasThreshold[fit]);\n", " h->Fill((fthreshold-threshold)/ethreshold);\n", " h = dynamic_cast(hbiasWidth[fit]);\n", " h->Fill((fwidth-width)/ewidth);\n", " }\n", "}\n", "\n", "\n", "TCanvas* c1 = new TCanvas(\"c1\",\n", " \"Efficiency fit biases\",10,10,1000,800);\n", "c1->Divide(2,2);\n", "\n", "TH1D *h0, *h1;\n", "c1->cd(1);\n", "h0 = dynamic_cast(hbiasNorm[0]);\n", "h0->Draw(\"HIST\");\n", "h1 = dynamic_cast(hbiasNorm[1]);\n", "h1->SetLineColor(kRed);\n", "h1->Draw(\"HIST SAMES\");\n", "TLegend* l1 = new TLegend(0.1, 0.75, 0.5, 0.9,\n", " \"plateau parameter\", \"ndc\");\n", "l1->AddEntry(h0, Form(\"histogram: mean = %4.2f RMS = \\\n", " %4.2f\", h0->GetMean(), h0->GetRMS()), \"l\");\n", "l1->AddEntry(h1, Form(\"binomial : mean = %4.2f RMS = \\\n", " %4.2f\", h1->GetMean(), h1->GetRMS()), \"l\");\n", "l1->Draw();\n", "\n", "c1->cd(2);\n", "h0 = dynamic_cast(hbiasThreshold[0]);\n", "h0->Draw(\"HIST\");\n", "h1 = dynamic_cast(hbiasThreshold[1]);\n", "h1->SetLineColor(kRed);\n", "h1->Draw(\"HIST SAMES\");\n", "TLegend* l2 = new TLegend(0.1, 0.75, 0.5, 0.9,\n", " \"threshold parameter\", \"ndc\");\n", "l2->AddEntry(h0, Form(\"histogram: mean = %4.2f RMS = \\\n", " %4.2f\", h0->GetMean(), h0->GetRMS()), \"l\");\n", "l2->AddEntry(h1, Form(\"binomial : mean = %4.2f RMS = \\\n", " %4.2f\", h1->GetMean(), h1->GetRMS()), \"l\");\n", "l2->Draw();\n", "\n", "c1->cd(3);\n", "h0 = dynamic_cast(hbiasWidth[0]);\n", "h0->Draw(\"HIST\");\n", "h1 = dynamic_cast(hbiasWidth[1]);\n", "h1->SetLineColor(kRed);\n", "h1->Draw(\"HIST SAMES\");\n", "TLegend* l3 = new TLegend(0.1, 0.75, 0.5, 0.9, \"width parameter\", \"ndc\");\n", "l3->AddEntry(h0, Form(\"histogram: mean = %4.2f RMS = \\\n", " %4.2f\", h0->GetMean(), h0->GetRMS()), \"l\");\n", "l3->AddEntry(h1, Form(\"binomial : mean = %4.2f RMS = \\\n", " %4.2f\", h1->GetMean(), h1->GetRMS()), \"l\");\n", "l3->Draw();\n", "\n", "c1->cd(4);\n", "hChisquared->Draw(\"HIST\");" ] }, { "cell_type": "markdown", "id": "b07b4032", "metadata": {}, "source": [ "Draw all canvases " ] }, { "cell_type": "code", "execution_count": 6, "id": "97d671b5", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:25:22.727380Z", "iopub.status.busy": "2026-05-19T20:25:22.727251Z", "iopub.status.idle": "2026-05-19T20:25:22.952270Z", "shell.execute_reply": "2026-05-19T20:25:22.951571Z" } }, "outputs": [ { "data": { "text/html": [ "\n", "\n", "
\n", "
\n", "\n", "\n", "\n" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "%jsroot on\n", "gROOT->GetListOfCanvases()->Draw()" ] } ], "metadata": { "kernelspec": { "display_name": "ROOT C++", "language": "c++", "name": "root" }, "language_info": { "codemirror_mode": "text/x-c++src", "file_extension": ".C", "mimetype": " text/x-c++src", "name": "c++" } }, "nbformat": 4, "nbformat_minor": 5 }