{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "ca8fc493",
   "metadata": {},
   "source": [
    "# goftest\n",
    "Example showing usage of goodness of fit tests.\n",
    "\n",
    "Using Anderson-Darling and Kolmogorov-Smirnov goodness of fit tests,\n",
    "1st sample test is performed comparing data with a log-normal distribution\n",
    "and a 2nd sample test is done comparing two gaussian data sets.\n",
    "\n",
    "\n",
    "\n",
    "\n",
    "**Author:** Bartolomeu Rabacal  \n",
    "<i><small>This notebook tutorial was automatically generated with <a href= \"https://github.com/root-project/root/blob/master/documentation/doxygen/converttonotebook.py\">ROOTBOOK-izer</a> from the macro found in the ROOT repository  on Tuesday, May 19, 2026 at 08:25 PM.</small></i>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "24379811",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:25:30.854305Z",
     "iopub.status.busy": "2026-05-19T20:25:30.854126Z",
     "iopub.status.idle": "2026-05-19T20:25:30.861547Z",
     "shell.execute_reply": "2026-05-19T20:25:30.860904Z"
    }
   },
   "outputs": [],
   "source": [
    "%%cpp -d\n",
    "#include <cassert>\n",
    "#include \"TCanvas.h\"\n",
    "#include \"TPaveText.h\"\n",
    "#include \"TH1.h\"\n",
    "#include \"TF1.h\"\n",
    "#include \"Math/GoFTest.h\"\n",
    "#include \"Math/Functor.h\"\n",
    "#include \"TRandom3.h\"\n",
    "#include \"Math/DistFunc.h\""
   ]
  },
  {
   "cell_type": "markdown",
   "id": "22ae87b3",
   "metadata": {},
   "source": [
    " need to use Functor1D\n",
    " "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "366d4041",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:25:30.863190Z",
     "iopub.status.busy": "2026-05-19T20:25:30.863009Z",
     "iopub.status.idle": "2026-05-19T20:25:30.866746Z",
     "shell.execute_reply": "2026-05-19T20:25:30.866034Z"
    }
   },
   "outputs": [],
   "source": [
    "%%cpp -d\n",
    "double landau(double x) {\n",
    "   return ROOT::Math::landau_pdf(x);\n",
    "}"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "14f0cddc",
   "metadata": {},
   "source": [
    "------------------------------------------------------------------------\n",
    "Case 1: Create logNormal random sample"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "7279ec3e",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:25:30.868192Z",
     "iopub.status.busy": "2026-05-19T20:25:30.868027Z",
     "iopub.status.idle": "2026-05-19T20:25:31.200127Z",
     "shell.execute_reply": "2026-05-19T20:25:31.199524Z"
    }
   },
   "outputs": [],
   "source": [
    "UInt_t nEvents1 = 1000;"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "04c36fa9",
   "metadata": {},
   "source": [
    "ROOT::Math::Random<ROOT::Math::GSLRngMT> r;"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "056bc544",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:25:31.202280Z",
     "iopub.status.busy": "2026-05-19T20:25:31.202159Z",
     "iopub.status.idle": "2026-05-19T20:25:31.416272Z",
     "shell.execute_reply": "2026-05-19T20:25:31.415555Z"
    }
   },
   "outputs": [],
   "source": [
    "TF1 * f1 = new TF1(\"logNormal\",\"ROOT::Math::lognormal_pdf(x,[0],[1])\",0,500);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "bf111c83",
   "metadata": {},
   "source": [
    "set the lognormal parameters (m and s)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "cce382ed",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:25:31.418162Z",
     "iopub.status.busy": "2026-05-19T20:25:31.418041Z",
     "iopub.status.idle": "2026-05-19T20:25:31.620404Z",
     "shell.execute_reply": "2026-05-19T20:25:31.619755Z"
    }
   },
   "outputs": [],
   "source": [
    "f1->SetParameters(4.0,1.0);\n",
    "f1->SetNpx(1000);\n",
    "\n",
    "\n",
    "Double_t* sample1 = new Double_t[nEvents1];\n",
    "\n",
    "TH1D* h1smp = new TH1D(\"h1smp\", \"LogNormal distribution histogram\", 100, 0, 500);\n",
    "h1smp->SetStats(kFALSE);\n",
    "\n",
    "for (UInt_t i = 0; i < nEvents1; ++i) {\n",
    "   //Double_t data = f1->GetRandom();\n",
    "   Double_t data = gRandom->Gaus(4,1);\n",
    "   data = TMath::Exp(data);\n",
    "   sample1[i] = data;\n",
    "   h1smp->Fill(data);\n",
    "}"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "9a6f435b",
   "metadata": {},
   "source": [
    "normalize correctly the histogram using the entries inside"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "fd198e80",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:25:31.622314Z",
     "iopub.status.busy": "2026-05-19T20:25:31.622194Z",
     "iopub.status.idle": "2026-05-19T20:25:31.826984Z",
     "shell.execute_reply": "2026-05-19T20:25:31.826280Z"
    }
   },
   "outputs": [],
   "source": [
    "h1smp->Scale( ROOT::Math::lognormal_cdf(500.,4.,1) / nEvents1, \"width\");\n",
    "\n",
    "TCanvas* c = new TCanvas(\"c\",\"1-Sample and 2-Samples GoF Tests\");\n",
    "c->Divide(1, 2);\n",
    "TPad * pad = (TPad *)c->cd(1);\n",
    "h1smp->Draw();\n",
    "h1smp->SetLineColor(kBlue);\n",
    "pad->SetLogy();\n",
    "f1->SetNpx(100); // use same points as histo for drawing\n",
    "f1->SetLineColor(kRed);\n",
    "f1->Draw(\"SAME\");"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "57a6ad0e",
   "metadata": {},
   "source": [
    "-----------------------------------------\n",
    "Create GoFTest object"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "33f91010",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:25:31.829079Z",
     "iopub.status.busy": "2026-05-19T20:25:31.828966Z",
     "iopub.status.idle": "2026-05-19T20:25:32.035521Z",
     "shell.execute_reply": "2026-05-19T20:25:32.034999Z"
    }
   },
   "outputs": [],
   "source": [
    "ROOT::Math::GoFTest* goftest_1 = new ROOT::Math::GoFTest(nEvents1, sample1, ROOT::Math::GoFTest::kLogNormal);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "f21146ae",
   "metadata": {},
   "source": [
    "----------------------------------------------------\n",
    "Possible calls for the Anderson - DarlingTest test\n",
    "a) Returning the Anderson-Darling standardized test statistic"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "f19c203c",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:25:32.037687Z",
     "iopub.status.busy": "2026-05-19T20:25:32.037549Z",
     "iopub.status.idle": "2026-05-19T20:25:32.244100Z",
     "shell.execute_reply": "2026-05-19T20:25:32.243536Z"
    }
   },
   "outputs": [],
   "source": [
    "Double_t A2_1 = goftest_1-> AndersonDarlingTest(\"t\");\n",
    "Double_t A2_2 = (*goftest_1)(ROOT::Math::GoFTest::kAD, \"t\");\n",
    "assert(A2_1 == A2_2);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "450f55fa",
   "metadata": {},
   "source": [
    "b) Returning the p-value for the Anderson-Darling test statistic"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "fda85723",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:25:32.246262Z",
     "iopub.status.busy": "2026-05-19T20:25:32.246150Z",
     "iopub.status.idle": "2026-05-19T20:25:32.452644Z",
     "shell.execute_reply": "2026-05-19T20:25:32.452079Z"
    }
   },
   "outputs": [],
   "source": [
    "Double_t pvalueAD_1 = goftest_1-> AndersonDarlingTest(); // p-value is the default choice\n",
    "Double_t pvalueAD_2 = (*goftest_1)(); // p-value and Anderson - Darling Test are the default choices\n",
    "assert(pvalueAD_1 == pvalueAD_2);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "916a8b64",
   "metadata": {},
   "source": [
    "Rebuild the test using the default 1-sample constructor"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "873fe5fc",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:25:32.454857Z",
     "iopub.status.busy": "2026-05-19T20:25:32.454730Z",
     "iopub.status.idle": "2026-05-19T20:25:32.661357Z",
     "shell.execute_reply": "2026-05-19T20:25:32.660715Z"
    }
   },
   "outputs": [],
   "source": [
    "delete goftest_1;\n",
    "goftest_1 = new ROOT::Math::GoFTest(nEvents1, sample1 ); // User must then input a distribution type option\n",
    "goftest_1->SetDistribution(ROOT::Math::GoFTest::kLogNormal);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "c46a73ba",
   "metadata": {},
   "source": [
    "--------------------------------------------------\n",
    "Possible calls for the Kolmogorov - Smirnov test\n",
    "a) Returning the Kolmogorov-Smirnov standardized test statistic"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "2a64d5a2",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:25:32.663081Z",
     "iopub.status.busy": "2026-05-19T20:25:32.662975Z",
     "iopub.status.idle": "2026-05-19T20:25:32.869540Z",
     "shell.execute_reply": "2026-05-19T20:25:32.869036Z"
    }
   },
   "outputs": [],
   "source": [
    "Double_t Dn_1 = goftest_1-> KolmogorovSmirnovTest(\"t\");\n",
    "Double_t Dn_2 = (*goftest_1)(ROOT::Math::GoFTest::kKS, \"t\");\n",
    "assert(Dn_1 == Dn_2);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "3bbd0bdb",
   "metadata": {},
   "source": [
    "b) Returning the p-value for the Kolmogorov-Smirnov test statistic"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "92557273",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:25:32.871567Z",
     "iopub.status.busy": "2026-05-19T20:25:32.871450Z",
     "iopub.status.idle": "2026-05-19T20:25:33.077857Z",
     "shell.execute_reply": "2026-05-19T20:25:33.077377Z"
    }
   },
   "outputs": [],
   "source": [
    "Double_t pvalueKS_1 = goftest_1-> KolmogorovSmirnovTest();\n",
    "Double_t pvalueKS_2 = (*goftest_1)(ROOT::Math::GoFTest::kKS);\n",
    "assert(pvalueKS_1 == pvalueKS_2);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "fc65b972",
   "metadata": {},
   "source": [
    "Valid but incorrect calls for the 2-samples methods of the 1-samples constructed goftest_1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "id": "a20d940f",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:25:33.079920Z",
     "iopub.status.busy": "2026-05-19T20:25:33.079811Z",
     "iopub.status.idle": "2026-05-19T20:25:33.298839Z",
     "shell.execute_reply": "2026-05-19T20:25:33.296648Z"
    }
   },
   "outputs": [],
   "source": [
    "#ifdef TEST_ERROR_MESSAGE\n",
    " Double_t A2 = (*goftest_1)(ROOT::Math::GoFTest::kAD2s, \"t\");     // Issues error message\n",
    " Double_t pvalueKS = (*goftest_1)(ROOT::Math::GoFTest::kKS2s);    // Issues error message\n",
    " assert(A2 == pvalueKS);\n",
    "#endif\n",
    "\n",
    "TPaveText* pt1 = new TPaveText(0.58, 0.6, 0.88, 0.80, \"brNDC\");\n",
    "Char_t str1[50];\n",
    "sprintf(str1, \"p-value for A-D 1-smp test: %f\", pvalueAD_1);\n",
    "pt1->AddText(str1);\n",
    "pt1->SetFillColor(18);\n",
    "pt1->SetTextFont(20);\n",
    "pt1->SetTextColor(4);\n",
    "Char_t str2[50];\n",
    "sprintf(str2, \"p-value for K-S 1-smp test: %f\", pvalueKS_1);\n",
    "pt1->AddText(str2);\n",
    "pt1->Draw();"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "97b61dfb",
   "metadata": {},
   "source": [
    "------------------------------------------------------------------------\n",
    "Case 2: Create Gaussian random samples"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "id": "885a8616",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:25:33.309987Z",
     "iopub.status.busy": "2026-05-19T20:25:33.309852Z",
     "iopub.status.idle": "2026-05-19T20:25:33.516352Z",
     "shell.execute_reply": "2026-05-19T20:25:33.515802Z"
    }
   },
   "outputs": [],
   "source": [
    "  \n",
    "UInt_t nEvents2 = 2000;\n",
    "\n",
    "Double_t* sample2 = new Double_t[nEvents2];\n",
    "\n",
    "TH1D* h2smps_1 = new TH1D(\"h2smps_1\", \"Gaussian distribution histograms\", 100, 0, 500);\n",
    "h2smps_1->SetStats(kFALSE);\n",
    "\n",
    "TH1D* h2smps_2 = new TH1D(\"h2smps_2\", \"Gaussian distribution histograms\", 100, 0, 500);\n",
    "h2smps_2->SetStats(kFALSE);\n",
    "\n",
    "TRandom3 r;\n",
    "for (UInt_t i = 0; i < nEvents1; ++i) {\n",
    "   Double_t data = r.Gaus(300, 50);\n",
    "   sample1[i] = data;\n",
    "   h2smps_1->Fill(data);\n",
    "}\n",
    "h2smps_1->Scale(1. / nEvents1, \"width\");\n",
    "c->cd(2);\n",
    "h2smps_1->Draw();\n",
    "h2smps_1->SetLineColor(kBlue);\n",
    "\n",
    "for (UInt_t i = 0; i < nEvents2; ++i) {\n",
    "   Double_t data = r.Gaus(300, 50);\n",
    "   sample2[i] = data;\n",
    "   h2smps_2->Fill(data);\n",
    "}\n",
    "h2smps_2->Scale(1. / nEvents2, \"width\");\n",
    "h2smps_2->Draw(\"SAME\");\n",
    "h2smps_2->SetLineColor(kRed);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "27a40591",
   "metadata": {},
   "source": [
    "-----------------------------------------\n",
    "Create GoFTest object"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "id": "ad59b688",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:25:33.518156Z",
     "iopub.status.busy": "2026-05-19T20:25:33.518048Z",
     "iopub.status.idle": "2026-05-19T20:25:33.724527Z",
     "shell.execute_reply": "2026-05-19T20:25:33.724042Z"
    }
   },
   "outputs": [],
   "source": [
    "ROOT::Math::GoFTest* goftest_2 = new ROOT::Math::GoFTest(nEvents1, sample1, nEvents2, sample2);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "f981ece1",
   "metadata": {},
   "source": [
    "----------------------------------------------------\n",
    "Possible calls for the Anderson - DarlingTest test\n",
    "a) Returning the Anderson-Darling standardized test statistic"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "id": "0669cdb8",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:25:33.726393Z",
     "iopub.status.busy": "2026-05-19T20:25:33.726287Z",
     "iopub.status.idle": "2026-05-19T20:25:33.932759Z",
     "shell.execute_reply": "2026-05-19T20:25:33.932242Z"
    }
   },
   "outputs": [],
   "source": [
    "A2_1 = goftest_2->AndersonDarling2SamplesTest(\"t\");\n",
    "A2_2 = (*goftest_2)(ROOT::Math::GoFTest::kAD2s, \"t\");\n",
    "assert(A2_1 == A2_2);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "83a942ec",
   "metadata": {},
   "source": [
    "b) Returning the p-value for the Anderson-Darling test statistic"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "id": "3be8f2f0",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:25:33.934722Z",
     "iopub.status.busy": "2026-05-19T20:25:33.934595Z",
     "iopub.status.idle": "2026-05-19T20:25:34.141019Z",
     "shell.execute_reply": "2026-05-19T20:25:34.140521Z"
    }
   },
   "outputs": [],
   "source": [
    "pvalueAD_1 = goftest_2-> AndersonDarling2SamplesTest(); // p-value is the default choice\n",
    "pvalueAD_2 = (*goftest_2)(ROOT::Math::GoFTest::kAD2s);  // p-value is the default choices\n",
    "assert(pvalueAD_1 == pvalueAD_2);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "5e29d386",
   "metadata": {},
   "source": [
    "--------------------------------------------------\n",
    "Possible calls for the Kolmogorov - Smirnov test\n",
    "a) Returning the Kolmogorov-Smirnov standardized test statistic"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "id": "98410a69",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:25:34.142851Z",
     "iopub.status.busy": "2026-05-19T20:25:34.142743Z",
     "iopub.status.idle": "2026-05-19T20:25:34.349239Z",
     "shell.execute_reply": "2026-05-19T20:25:34.348740Z"
    }
   },
   "outputs": [],
   "source": [
    "Dn_1 = goftest_2-> KolmogorovSmirnov2SamplesTest(\"t\");\n",
    "Dn_2 = (*goftest_2)(ROOT::Math::GoFTest::kKS2s, \"t\");\n",
    "assert(Dn_1 == Dn_2);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "655213bb",
   "metadata": {},
   "source": [
    "b) Returning the p-value for the Kolmogorov-Smirnov test statistic"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "id": "839ae956",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:25:34.351082Z",
     "iopub.status.busy": "2026-05-19T20:25:34.350974Z",
     "iopub.status.idle": "2026-05-19T20:25:34.557405Z",
     "shell.execute_reply": "2026-05-19T20:25:34.556895Z"
    }
   },
   "outputs": [],
   "source": [
    "pvalueKS_1 = goftest_2-> KolmogorovSmirnov2SamplesTest();\n",
    "pvalueKS_2 = (*goftest_2)(ROOT::Math::GoFTest::kKS2s);\n",
    "assert(pvalueKS_1 == pvalueKS_2);\n",
    "\n",
    "#ifdef TEST_ERROR_MESSAGE\n",
    "/* Valid but incorrect calls for the 1-sample methods of the 2-samples constructed goftest_2 */\n",
    "A2 = (*goftest_2)(ROOT::Math::GoFTest::kAD, \"t\");     // Issues error message\n",
    "pvalueKS = (*goftest_2)(ROOT::Math::GoFTest::kKS);    // Issues error message\n",
    "assert(A2 == pvalueKS);\n",
    "#endif\n",
    "\n",
    "TPaveText* pt2 = new TPaveText(0.13, 0.6, 0.43, 0.8, \"brNDC\");\n",
    "sprintf(str1, \"p-value for A-D 2-smps test: %f\", pvalueAD_1);\n",
    "pt2->AddText(str1);\n",
    "pt2->SetFillColor(18);\n",
    "pt2->SetTextFont(20);\n",
    "pt2->SetTextColor(4);\n",
    "sprintf(str2, \"p-value for K-S 2-smps test: %f\", pvalueKS_1);\n",
    "pt2-> AddText(str2);\n",
    "pt2-> Draw();"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "fba77720",
   "metadata": {},
   "source": [
    "------------------------------------------------------------------------\n",
    "Case 3: Create Landau random sample"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "id": "331b606d",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:25:34.559198Z",
     "iopub.status.busy": "2026-05-19T20:25:34.559090Z",
     "iopub.status.idle": "2026-05-19T20:25:34.761072Z",
     "shell.execute_reply": "2026-05-19T20:25:34.760574Z"
    }
   },
   "outputs": [],
   "source": [
    "UInt_t nEvents3 = 1000;\n",
    "\n",
    "Double_t* sample3 = new Double_t[nEvents3];\n",
    "for (UInt_t i = 0; i < nEvents3; ++i) {\n",
    "   Double_t data = r.Landau();\n",
    "   sample3[i] = data;\n",
    "}"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "bc1a03cf",
   "metadata": {},
   "source": [
    "------------------------------------------\n",
    "Create GoFTest objects\n",
    "\n",
    "Possible constructors for the user input distribution"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "9f44cf0d",
   "metadata": {},
   "source": [
    "a) User input PDF"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "id": "9d4ea5a2",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:25:34.763021Z",
     "iopub.status.busy": "2026-05-19T20:25:34.762915Z",
     "iopub.status.idle": "2026-05-19T20:25:34.967373Z",
     "shell.execute_reply": "2026-05-19T20:25:34.966895Z"
    }
   },
   "outputs": [],
   "source": [
    "ROOT::Math::Functor1D f(&landau);\n",
    "double minimum = 3*TMath::MinElement(nEvents3, sample3);\n",
    "double maximum = 3*TMath::MaxElement(nEvents3, sample3);\n",
    "ROOT::Math::GoFTest* goftest_3a = new ROOT::Math::GoFTest(nEvents3, sample3, f,  ROOT::Math::GoFTest::kPDF, minimum,maximum);  // need to specify am interval"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "d7cf6e0b",
   "metadata": {},
   "source": [
    "b) User input CDF"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "id": "a30132d3",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:25:34.969260Z",
     "iopub.status.busy": "2026-05-19T20:25:34.969147Z",
     "iopub.status.idle": "2026-05-19T20:25:35.175591Z",
     "shell.execute_reply": "2026-05-19T20:25:35.175124Z"
    }
   },
   "outputs": [],
   "source": [
    "ROOT::Math::Functor1D fI(&TMath::LandauI);\n",
    "ROOT::Math::GoFTest* goftest_3b = new ROOT::Math::GoFTest(nEvents3, sample3, fI, ROOT::Math::GoFTest::kCDF,minimum,maximum);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "e81476b7",
   "metadata": {},
   "source": [
    "Returning the p-value for the Anderson-Darling test statistic"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "id": "f8df2dce",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:25:35.177546Z",
     "iopub.status.busy": "2026-05-19T20:25:35.177433Z",
     "iopub.status.idle": "2026-05-19T20:25:35.381933Z",
     "shell.execute_reply": "2026-05-19T20:25:35.381448Z"
    }
   },
   "outputs": [],
   "source": [
    "pvalueAD_1 = goftest_3a-> AndersonDarlingTest(); // p-value is the default choice\n",
    "\n",
    "pvalueAD_2 = (*goftest_3b)(); // p-value and Anderson - Darling Test are the default choices"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "d0eacaf6",
   "metadata": {},
   "source": [
    "Checking consistency between both tests"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "id": "39a70518",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:25:35.383800Z",
     "iopub.status.busy": "2026-05-19T20:25:35.383691Z",
     "iopub.status.idle": "2026-05-19T20:25:35.590198Z",
     "shell.execute_reply": "2026-05-19T20:25:35.589824Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      " \n",
      "\n",
      "TEST with LANDAU distribution:\tOK ( pvalues = 0.777278  )\n"
     ]
    }
   ],
   "source": [
    "std::cout << \" \\n\\nTEST with LANDAU distribution:\\t\";\n",
    "if (TMath::Abs(pvalueAD_1 - pvalueAD_2) > 1.E-1 * pvalueAD_2) {\n",
    "   std::cout << \"FAILED \" << std::endl;\n",
    "   Error(\"goftest\",\"Error in comparing testing using Landau and Landau CDF\");\n",
    "   std::cerr << \" pvalues are \" << pvalueAD_1 << \"  \" << pvalueAD_2 << std::endl;\n",
    "}\n",
    "else\n",
    "   std::cout << \"OK ( pvalues = \" << pvalueAD_2 << \"  )\" << std::endl;"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "8f26e7dd",
   "metadata": {},
   "source": [
    "Draw all canvases "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "id": "739abcff",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:25:35.591886Z",
     "iopub.status.busy": "2026-05-19T20:25:35.591781Z",
     "iopub.status.idle": "2026-05-19T20:25:35.809037Z",
     "shell.execute_reply": "2026-05-19T20:25:35.808684Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "\n",
       "\n",
       "<div id=\"root_plot_1779222335799\" style=\"width: 700px; height: 500px; position: relative\">\n",
       "</div>\n",
       "\n",
       "</div>\n",
       "<script>\n",
       "   function process_root_plot_1779222335799() {\n",
       "      function execCode(Core) {\n",
       "         Core.settings.HandleKeys = false;\n",
       "         \n",
       "Core.unzipJSON(39898,'WkwIqS0A2psAeAHtnX2PmzeS4L+KIewfdwAtsPjO5zkc4DjxZG7t2Iidib2DwUB2q21d2q1etZw4s8h3P/yKfCR1t+zrzCTZycBB3FLxtVgki1XFKuq/Zn/d/nixPF+8Xc6G2bP7i/PvF5ffLl8+PV9cXL5Zb2dmdvrN+eo/3y3/+PlssGZ2+tlqe9m+PX75f5evtqTPKPb4Yrtan3fg31fnJ7PBm9nprqXhv4719aEOvA/Oh2Bmpw9X58v767P1ZjZIB59ufzxb7sFvVyfbNw18sDo764VBFnAqbK1WX55uHy02r1fns8HOSfl69frNtaTP1tvt+u3VYs/WF1cTnp+uQMKZ2emL/dd77SsNP7/cLrb0UitlrkD3GkTlB5vF2+V1vEm7NvBduasD2hWdkqcmqX5IGZr8bL05WW6erv7WqXeQ+Gh9smzz+lza54v++dw12r7on8+363svL5+s3i/Pvus1tutj4GxINTH0KxVCVpLtq1xNmA13W8I3V2oxm9/s6xyCUzd/ulKhNfKnfZWrCbtupqFs18/7aHQoN0Fm0loJPuUiyVapkeW5r//iu0anXh8QNCdwNtylASdSUgg25VxT0nVytv7hq8/vd6IfAs+/udAMSr04+P7tLvXL3bd7Ly+vNHTv5eWVtu69vNxXu/fycl/zm/dvdTdA4R/3X9+/XbxvA/rmx93XZ2+W28Vs8DqwN6v+7d7lxfLV9uvFdrVuo/jq3duXy037/mz16rv3+68/tq8P16974sP1633a31ruk8XJk8XqnN1jZqf3N+vLyzeLVW9wBz5Zd050uKxZXA3er+hH65PV6Wp5MhtOF2eXSzM7/cNmdfL+KvjjHrz38vL+er05KP/FyWq7eMl+327e0cCD1fvlyZVxT00/2azerrar75eXN1jew9UlHHXith1cbDaz4c9/MbP1xZYvP5nZ6Rfvl68uZ8P5u7MzMzv9qrHnVzDLZ6staMzk7tPF24uz5Z3F+ckd14HLO39YP7jzbHm5Vcb61bu3TxZny+12YqBQ86vl++3N1M//+PTJw3svZsPs36avZnb6+frdy7PlZ+9OT6fZ/Hq5XazOIWUnxPPL1d+W31xO+S+ugpr79XJxNhscnWt2gyXMXYmZtf3t6vxk/cOz9QWb7BB+cQh3ZrYv8OUS1t3XyA8TG7j/Ztb3+v3FdntjDu5tt+2UY3zPP1tuf1guzzsbvwIpoR9s1m+frS9mg8xZWM9PFlu4pAIvJoBj7l4D5Ccz++7R+vvl44vFf77bLZfvvl5Cl6uJp1+uXr95yBD6iaVLd7F99WYi7ndP36x/+OL75fn26XaxfXe5W6Hf3Xu3XbNGdiUfLc/ffbbYNJgVdO8Va3BX4/Tr5eLk8fnZj1ON029X2zfrd9vD5Tot4S8Xl30BTimHpf587SD/xcQF2OkHxYVvly+VJazOX39IZmBl3D9bXF72/UK5JqQcJlywVGd2cDGa/m+UwRpr7Og0lW9+sJprx6B5LsYx7urYMfWSpOdelu9lKNU4EVPCWPlegnGSRqHDYPq/UWSQ6owUZyS7UdyQk9H/R/GDWGf6v1HCIC6b/m+UOEj0pv8bJQ2Si+n/RsmDs8H0f6OUwblq+r9R6uCCM/3f6KwWlkp9Ozq5CrpBqjdSxUgqo/ODlGSkAPrRhdZvikaiHV0cpHhyjMQ4OrDqSAY7ujxIjEaCM+Lt6MogPuoASxpdHSRH42wCk9HbQbwz2nXyo5dBfAYDI8mO3g3irZFYDA173/oppdWFVK71E/3o4yCSjdYPcfRpEBvAQFH2eRBbjVbwbvRlEBeUrs7W0dc91WMcgx0kW6PDrTIGGZhixdnHMbjBOTFKH0BIVRqS3o4hDEpiZs2FMcRBaQoaUseQBpqpwdQ0hty6gVAhj6EM2qMommOog4RsdBDixmgVCV1iY5RBQjK6ZsfolA4d8EMt0/cw5DB9j0OU6XsaXJ6+58FPXwsLuzdTB2sCCziOyQ7WpNK+s3WqtO9uYFC6rcbEDhJfWw6biNWh1SNAyQ1I9CGtk6Q7yfdetHfdpDImulfA5THTvwI2jBkEAKTYMYOBAjGNedrDzG8GA82xZcxgAFDimBWBGE2SMU872ecxT92LH3MdGobGjsUODUG+yxB1GHx3Q1Ey8J0t3LiLHQurUskA0LYv/dmxsE3SBLB1G1Z2LGxcJQMA21apbexYO+PyBUAaR5IA0JiXVNhX9ZojKQGEBgQH0LgYi8qOtbGxqi1nLZWF9KLfQ+Z71e/OGzuKtYOImBCVKYgVXaHBmeBGsa4xhcr0yijWDyUb0VVDdtBxSEzGA8bBuWJ8NjGMYtMg0ZrMBq2jWFhG0rKSaIrt6Yz3xjlASGISVJVRlLdG9pAbRWBjYrwYB+SUqjCjVEeBt/Kdfy6OAnNlNiE87UpsGLf1LHBXK6ZWI9aOInmAwTEeKaPAXVl3qRqBo9dBSqVjGcXZQQI8KxhJ8Hf4WKNbAXSDVrPZlDqK84Pz2o8dBdZaIseCYZsLvFW8UVbmKJx0qzspDWe468ER4KwMVdiSVUbHlMAS27/RWXiTN/3f6JgTYTz6b3Q2Nn4LWx8dk1K8sgA7OpsH9mYWk2m5DM7Xhpbk0dk6UAtmHUfHlLAKqDw6kSGLaZVHJ075cq89OiaF/qg8OgnUbJVHJ5GarfLomA84umRqj05yw3bqtOwGp71WxXfqluMOdtSQHp1OCdi2fpmSNlYpfnQ6JTpY58PoXNiNVntmUvpoQdrBxfdUdXkab6tcdgOmsp/2MLzNhtFP25jFxLnFTu6neIPbZtZ860bPfp7yFT6QCRRuu5ryUUbPxu7FFSw7CUJBeFsTN6KMwe4FIzsG22QCqgO5XU2gvZwCFDhOOFHGYKFG416ML+hmVt6v4wlsZmBO4+jHwNKxbje+ILbB5Fs3BoQlUOQEVdjRVx9eEK+QDYq/KB4UZTSihxyVFZpYvR2DKKtv4xJYPVVI56Chgh2D45xpQ3J9SprMOAbEoT4kjrTApu1D0iGzPvqQWj6cxOmQWz7SRx8SJHBlGlIrXqcxaWlvAelAC3tRsM3aGPx06GlVr6dem+8xeD14dxX15N3Vgxw6aWNkSfY1Qiuxym6qG7yfejCKdT/5LT+08bQxjpEl2fpRlGNNoNxmbYw1T5C2pWNvAxxj1aEDWTcmyxToqTgm26R06iTLmHXSxmQZsU7amGyT1pEYWIVthymKCfbVdliHiw65T9uYOFPIb9M2JmVhugJbl6zCtgJ1SImjRXdYa068NscQQZCTRRdhL9zoMQ1LDtWH1KV2uB412+7og+NUYTU2cjhdCMzamFiS0/ioB8fq42OVJFhW22Ed7gd/n/6kPKuNT0mtUrrusA63460tkOTqNDz6Rkpvo1Oo7c9e0u+VqOQbUZBIEtI5e0y/N3K077o0+JonRtjQHnNnhJ3sY2bVdRSjjJlF10Eqx914NXM//Qwvd/FGURhzF3DoSQu30bXMMvFA5Q4FHtjpGmUstnECnZSxdCaoTKLY3ZyT09FpvKf01djxHcu0GpXTFNZiX3vax47YyJiNI7cu+iKkMAU5SjtboaRywt7knhGS09lRG27hHJ0YsnVj4Rxty14Llx2kfex4EUKuLkC6B+inQ6POWDhCJ75OszBJVltjxgUe2UEqd0m4o9RXY6fCWPpqnHpqexUIlNAZe0cI2513NaRYjH1sbIuCythJxLYovhGprZux7BkkGHs9L8Cple1cux/UxU/6QSu859nakaooU9VwqKOUoBybdrTdoBy7bb6xBNWU+joJE8NmnGGPDtiFPTqtmR06DCzs0GmFJ3RoKE7YKNAp1I7osaAu9nnRgcTdMtJ2Ixh17lya1tgJOpbYMOoEH0sEI83FFBIVoz0IRp0/FtUh+3lQUjvQ2rSNpemRbdLGoopkZ5ylKZIN27Gkdry3pTCWdHDI0r+qk/tWG32mUatGuavaVMopL3dTyIRQ7qffhFJuC3siSm4Lu3HsktvCngaamzC0o0JGs5k07JLbwp4IgaLZMllRpS3sdvyV0qZtKlrarLVDoxSds3ZKlKIyUEcHjVPXVptcdM7OM5VGpW39qY+29Xd96NbvXdR+9mi1zqknCnRG3fvvwoEeU6Uz6Q40pjjNZufRu/nsPHpqZncA6YLuPFobqp1Fd2DPoaFbPWDRrOAKk+5nsXVj7Wy69VI7n+4jqROf7lNdO59uxKxdZpigJjH0ma6dUffMbl3TNQVSqBsQvv0bK8ySddz+jdXpDLYztLowRI9aLZ6quswDqnEeq0tDtKqBB5T+PISKoi6BRssQikHxDGmsrg4hmWiNhDJWb4cA+zQSZaxehhCMIhPG6t0QvFH2lMfq/RCcQf1OWBLCEKzB4pbcWH0cfDWxYgsbq0+DLyY1E2H1efDZoCJneiiDjwYFONNDxWySkpGcxhrs4L1J2UiuYw0yeIc2j25cgxu8qnxS4liDH1w16MIljzWEwRU1ElTMGhGzVY5YC8eKJS2ZnI1UeshYCHI1UukBZdUUVrKMNdTBOVOaXFfVjmZU9ctjhTlWUyI2qbHCGospyTjxY41YD0wpxkkcawyDJFOxXpWxRlUIMYQ5GSt8MZrqOsSWV/OEo03kO0wuDVK1uaYGJeQqU3OHWDKmoudjqsFmgC3CSR1r8kCC4UJBRUZXrYKx5fqemxrYDBc15QYyyDLWVBiWWEYJiEHWYJpQMGP6NJhKneSxZlEQc5iCmHwMJqEGKoVEXAexEhusF5jAK5yxGFFbAyDWhWYvVDA3kH7JRZk3gpUDasMZAUEjjrVYBbEUKCgNBA1ysV2r0aOBvoGgEcaKQY5chg+IQdaIY/iAqYGgAZgbCBp+rDBICoMGYFXQg4Yfa21YYXsSN1aUqIL1sYONVtibRMYKl8Q6DRqAbQY9aAC2GfSggQ1OF5R40ADUFaXWVQV1SUlgfdexotIHI1ixLLY0jHbALPgCjDEKhjHBWNWNBMwRmu8HwQKdjLOZ8lgwsUhPcGxwmWDMMkZCxSJO+axwBB2FOVVgO8bZSL6aUyWCD7DYofajWEFRUK8xyHUDxvRmKxErXkFw0dygIKh4cuOASQdMFEyACUSwOGLAw544gUVB0NDcOmDNgyqYGJ3VumodBBQ6gp1VaOocSHL9USGZirpqbpUKyZB1XYcZE9IuNki4E2irEa/DIKribrvmqOCCyQjM4V+KTFv+yiPpHgWscm1jpNC/R2LpMP2jhHG1pLcSwFgYm02xgI/HFtth8PEoh80uin3Ser0RUDtps+RmGK9gni1QwxfYMhc3cHCxsHf6BwYfGDz9N34vFhZP/8D0D5Onf04A+g+ew6LZMaFHCBwlDQafEDloBH6e6T8kjiEp0Jz+Q+aQarD2XzjC2m0V9AiVA67B0CNajj+9fkvgE4XDUc2iif6j4yBtMP13mVjvuJiPGIZIfQ50+kcuzkYKigP9xzRwRiqs/eUhgQ/52l8ZWDsKM/5YB664sPuqqTxZLIgNpr8kAzZr8gP9Jb0wbDD9JT9AG27oAv2lMHC1qDD9pTjQl8L0l5IauhVmfCkParVHjtD+ysDlH7Ca41MdKrd/wPSXLbdKDYa+HAsWhOF3dMjBgIKjCfSoMjMkUo5IQhg4RloCfXI8sGlK1WtCsRwQWFlJcNorjIVuSdBuuRagWxK0W+7t2KNWL9vEclA4rhZJYGY5Kli6mkC3HBZ640gJ1hbHBYuRm0/s/JYDg9VICaFbjgyW3z4Bmzy3lpSA/hwbLDBNADEOjtDx4LbAcnSwxLQEiHF4IJ1pAgSabHCaAKZqhOv0EDBVe0inh3J6PUSgKTcWoK7HSKep8nIOEl15lABTjhLYkiaAKYdJ4zNcNZKAQNDnlgS9BFJOUIzAZvUaiK2tCdy6TDdBmqA3NNw7TG1wlcKZohcRxcBPuQ3S1QzM3YpePTC0YmCoeiEEIUrmolT0Qkgv0LOBwXIjJJX1lg0MljshqQwrGxgs10LCPQ6w3hJhxgGhCW6WKoVBWK8hGBIiHjCaC/gkvbHhaohbS8rDgLkb4kpYYfDjckgXc+GSUbgd4pKF8cBwuSByuraLKeCHLQWOz3Uu/WPd60sdhqtXRHD4ag0Ml1sivXMGpn+se3D8KiZDL+wpuhHEwHC5JXJweGDtnxtG1rQzMFzBvgeHB4Y+GFV0m3gDwxUOGK62qjftIs0NTndNMDBc0QOGBR8MDFc4YHQTIfQDx8F51j8qAXAauMNUmP45YHRLJQPDFQ6YQH+Ji1kR33wXmFMYrgQ7ON1g2cBwhSt5ODgw/XMnr/utGBiwoEfodisGBixoEkG1AxPpH11Cd181MGBBm1ABRVVMEfQJtTNYE6EHGoXCzsCQBZ1CYW8C/aNVKBwNDFkiOiqiRTIwZIlYSVUAMzBkiVhvkUSqCeATubNFcBS9WpUIfsDewKCFK3uFg4FBSwQ/8pOBQUsEP+DMZbZIBD/gajz4JPBDW7EGfi0J/ICdgV1LAj/gYDz4JPADjgZmLQn8gLPx0CM1+rGGuKSUBH7WOERj8EngBywGzi0J/IC9cYoP+AFzbc5VK/gBJwPbFkwyCqMNAYMf+UjiwMxvNawZmLjkNr/4ugj45Njz2YPAqcNIJMDoi9RH6gJm/QFzuYsj0WfvTvEl+zf11ps9OFsvtt7NzOxMvcdiNLPvZ8OfqwumOlS/ZKqDtxRTXTXVW1O54vXOVO9N9cFUH031yVQPjyqm+mpqsKYGMTU4U4M3NQRTA8ojPiLwtmJqqKZGa2oUU6MzNXqDL2rlDjwmUyM8sZgaq6nJmprE1ORMTd7UhKNJNDUlUxO8tJiaqqnZmprF1OxMzd7UHEzNqKLJ1JxNzfDhamqxphYxtThTize1BFMLSmoyFV6q/LuaWq2p3EFXZ2r1puLgUqOpFRU2m1rh+3qvjh0GTmzVEmK5gVDLqOUuBhZsERsszNfisGNhuxZea2GwFrnawlot/NSqVgwntbBPi5RsYZwWOcHCIi2SsIU5WjiiRTyw8EILA7RwPYswa+F3FiZnkV0t7M2q4xKMzCIRWFiYhW9ZxE8Lx7IsOYu0aVlsFq5kWWa2Oe9Qg1VlMcxYeI6N1FCrHEKl5Yy3mF8sMqXlhLdIk5aj3WJtsZzpFkuLRSOxnOgW44rlKLeIkpYz3KIgWA5viwnFIkNajm6LpG45sy0yuMWBwCJBWkRqy1FtESAtZ7RFdLTInhaB0SIFWj3aVWjgWLYIS5YD2SI3Wo5iW6mB2Gg5iS0Co2XOsVogNPAHJwfmvHk1MecYIjjv+YPGz5yroIiVQVRCbP4QzLnKh8KcN8cA5lz9pTANcKryRx0p6IM5R+/neOQPfTDnKhA2jyzmXMVBYc5R1jm7+EMN5hwlnPMIQwZ9MOcq/KFQc7DwhxrMuTpsqZyHesyxwB9qMOdqDBTmXP3j1ONFmHN1cEONhVHzhxrMuXrYqRbaPMOYcxXf1ENO5TZVClVgay4kzLmKa6qdqZymepfk8peffvrJ/FrenPFj3pwtWOEjwSXdYxMf7s05BSe/55ZyZ9fAtSgJwINIiOtBEPvwECEo42Z8yKPF5rvl5iDepCUcNNkTdjEUz5bvt/fOX+OUjQMqYMu0cwsNNP9s9fp8hlWkwQftk/1gjat7Ug/jxfvVTcfxe9vtPdJx3z5Zfb+6XK3PL2dDFHok56DBh4uXyykGhv4Ubj0EelD48enp5VKDU+CzPXGHtle8V6++e7g8f01ojZ1bnJx1DqaqOhac0a9X255N7ua7IlP3OCa/+N2MUKfz7xjhf/xuRriboJ85h58tNgehRZ8tNtOiUBdsgrXYpmdPnrY98flm8UOLyWjw44vtPv6jAT0EpAE9CuTxxfbz5m/fwstwgmcT6TZ6fLHtLIFBPL7YPtBwrF70warHBdzwoKeAJp6stsSUTfCz9fpMHehJaPEr99fn2/W7zWWPXbi37ehc45j3tls2sTKpj/AC9zOZAXuF4bfAIwYJRETEXKEvzk++2GzWPZqLna2gFqerB+/OX3W2QCbgARcD7FNILvEwvTDjB+yF2feAB/P9cPl6eX5yGHQDdi31gMPS0D5x6nsKw6OJiS3sCkK/vhTN7PRLYiGWl9d4eE99erF4RSyA9r2LjTsYwy4wrqeB467cVWx2Rafkqei1rrXc9XEfJO4jjr5cXbIgD/EhifY6Osky7Knc1HGjzlR0Su0Fr2FDqUer89Xbd2//Y7lZ70M9yLgSrqgsvkW9PNksT5ebPzzcl27pB4RrCYfDBNPD1P04W+rny9MvZ0O0TNou5dvZkK+mPJ9xD3BQ5EVPeLI4XHxPFlfWFp3vkvY9a9LNMM0ni5MrQ4d2TxYnN8M+nyxOjkR+PlmcsNif78nTU15cSeFs7LFKdLh69V2PVHqyuGjxlc8709glvJhxPTk7ffpqs1yeP1i8Uu4DerC1A/IDsg8Oli1Jh/Mx1TrYPxQB3O8eUqYF1JZVK7N5CxedpXnQQBoSlWfU0nFRFoNy2lD7dop3ouSXAKLj2G5WF58vX63eLs4ud0FFypK7tON2IsPB6LTAteFp2uH4WCSaeDDAHbw7qpqQggDSwoF3tXYsUOu04dRdacYztc9wGvNZ//H8fLn5muFRkq2mzV7Ohj8TSHTnDn/E3RGU/jvhTpngdCfckQTkLAVCL7KHNOkw3dmrxXd5xfa8Yu8EeyfcCVpS+/6V/vwF5rJcnCw3nNYaM6Vk20EPVtsH06KJfdFoJBazuMvRtaRUe7U408rM/v9Zr85JnASB+4uLQ/DZ6u1OnMylVClBOcYf3y5eL2lox+DvL85PzpbfvlldfrfcfL04f90DmFv6Z+v3Pa3NXktVTA5iOP+0Wp+tzqfUHr7Yit5fbV6dXef2PYvgU5A+OACfI3J/8f7i+aHYMyW+OEx8cazklHilJAUfLd5/vnqtAfYswseb7Zv1/cXb5WbRuc9Nde3J4uQXi79jxj4Yf/dkcfIhfQ1NCgodcCzAifn03AP6XTsZb2phN9n1EV79KUr/bnLzaK21IcWUMFvMNGr/bpjHbIu3NZdaXWFiieKP18pT3M2GuzJ3wadSQolSctLiV4Lpc5xLKtb6GjBF9QJTzH8qc2eLtcHWlFwmjHN6EYDAWVvFu1xLcD7mm+8B3M11XlPxNqVkBeOlvhgwNX63hHkqZMcYUou5nx4PmA13U5nbGkIqKUj2geavPhyQ5jXWWmsJIXiuXz72jkBxcwZircfGRGNXHxVwXuZOC6iVTrMnRJ1L85iojNF+nzUb7l7P2r8ZwJsDjKHkHGNMLsVGgP5oQMu+OclTNsdVleQlFIlWCqbg608S3JW55OCzt1KkYFLcv0lA70fWSm++P1kgIZUYik0u+egZ2v6hgbm1ELfWmnMKxUXou393YB6vZOu66e8Y8CwGxIJezkfRaf9IXnu8wM6rkp81H3IzVLSXDOw83Mi58ibCUVSvvJJwBNvdqwnHOt69onCs75uvKtwNc+cL5jr2JZbQGaV4ZaGJzv2dhbsyjzkF63IMvmL1n/3cdxeYpNu9u6DceT1p3J/eXWjvLlzT8G//bsP0UsNf5dBm2cHbPsfwy4XvdyHnN4neZyV9UHpQJflD8gM845eVHw6fsDm27fQwVFGzn383NtyhOsLIGtzVz5ti2C8mg32Uil/K5x+iYdTnNfryeyOXby8OF+DD9euvkOHP7pysLreb1ct3vNV0583qcrt+vVm8pezBDCQlzT8ow/18S/pXr5Zn6JLcCs1OP2AJ72bwDz1K1SnwXk2wk5rX9Zrftekcq+TLFVZ/rs9mp893L/U8PzhCnrcivCZz+mC1uexW0IeL6RvPd7lJA367/Hx1eXG2OHj6BpVsp9TpsuA6Yf8UzaP1ycPFyw5/xJx/u2n68Z96mnZmgp3JoV+MfOSG42CajkySpv3mU/SB+4jbTdHf/qmnCIKqueJ3PkXXrlAO7les7vYvzrebFcaoDj+7fPf2B7XFVe9zQGkrvilmmuVmg58XiakWh3JYfFzeVRMYFXnXK825tQ4B7cEVlTE1y2GnLHWekrO2uMrt/Oz00eI9dmZ0xXb90g3P+wTOl8mqycbplyftZavTp+/egtNV/5mEY436zyjH55Wl5kXj56X6mGry1ccY8vJuNjLn0TRnc0oxSXTLu8m4uU0+Vqm2lFhCS5IgIXlvJYUYpaW5WFy2vtrikuRezrrocrLRRamRNJnn6HCr8DnbWlMrF6KroZTscvVV+5B5yrxOWEOozpVeN8YcXC76IbW1V0vMzkZJNovraclHV6sLJbrqKVbmMcUSvbMpFmc9o63z7H11wdbiS0mpUcCK867g/eFdUAqUueQUanLBpZBC+VC5OLdJUrDR40cgFIvzEqyPIcRqvXcfSHLz7GKoJYccQymK3JHGjs3YsTTe/HLVpWQ9S0sRSfMk1abiS6o2QKZswjxl54qPfKSo4z+WJjfIfryun0cbnLiUY5BadbTHhnY8zfskQXLx2Fg6jbGHZIvN0qYUGt1rSLmEIjXlnJVUx3B2R9r7R9LyPJcsvjo+LGuqGM/uZi3yEbzuoGN9yC3HcXMuy9F+j+FyLE2O0OpYH/Zo6vGSx8ZyrOQxfG5b7ljd26Ydw88eTT3W4rHax6h4nGLHU4+NWh9C+1nlj2F7vOXjqfR2rI3jON++5DSSY61nfWXxypO13IijFd30vrm9+t9eStzfwU0vJ6o8/tnqnCu0Lzabxzz0yCEJ/Pj75eb0bP0DbgYY7zYb5PFbnpT6JGnN3om3Lnin70LNrc3e2lKKlIr91+AyYEMs1lXrSsWNT5NKdCmLHnheNElCCIlDFwfCXjHF6KP1XnywmmRLyTXEkmwN0Wat6IIXF5yNHh9fcCgBXpgsHSfnNC2X6kp0tmZX8F9VLDhxSwou2BB6VVdyqCFbJ9HXqlWDB6eSwKXyGtXc2lC9t97FZIuNvGKFO5N466MvyWNta/gGCT6H7CM+j729Y+VcTFFyql44klq/rsZapEafQvFRB3s0TXwRFyr+dy7y7haPvR5p79h0HUtzkm30weuZlRr5POd1chLEh8LzFPTBg7TYjquUVBoNrqfl+VW6L+9+uK7kHEuV7KRU/Grp49jYjqZJ0XNIQqyeiOJ5FaQ4m11NWTz91nlJyaYoTpJIqR/F5Vp7isvfmeY5coJLNvD+bk5TvxJzidkGqSHmNuXXu7jtMGR+dSvSxbFub5t2jFTH+uAUudnzsTR7dCzHSh7D8bbljtW9bdoxWh/H+liLx2ofo+Jxih1PPTbq6Ww5lncsjbPs5vI7XvJ4ajuTb7ZxHOdjvR0vOY3E6ubyVzv4FRyEn9579MV4oV5If71cfL/UW+X2XPxHbZ8P1KJ+zOyH+bib/c4mO+eh6e/rx4+fDcOjxfbNMJytX5+rIfSvFyen/+O9+bP9i/mz/OV/XjeA6lXqhw2gHNQHl9jXfYk52n+++fO4Re+rC171VbvSyeptc3f76gINXg0Dz3686I7HX12c7p7m/+rzB03AuP9mdfmf7xabXuaG1n5Dr3+y2KgbIX4v1ljsiE8Wm0c8EH4Ic2k1wU8XvDP8YYHFmwPVHg4eYhGPH32KBe980rjDdbWEXEtySdkwnvHBl5xyCJ5XUdB/bUgV3THZ0s5EPNql5JijBGI3KEVYQEiSrIsldSHDuxiCSwmJp3cZLReP0XrJxM6ARsno2IIkZYnjIS3XUKzLPliJqR+62UXvra25+kjYCOVSkhrFxZhKKtIElGRDqBFNzJVcW1qMLlaLcJJyjO3EiTbmXIoLNkoIrV/ksBJccBgAJiHIxRhcQXmW6JqQ4auArRTrhPgbcPHJOifVSrQ129av905KsSnkmktpgpa3KdcaneNav/fris+ZK2NfKK7tuYTYFmKIsRClQR8uIF95lGCbay/nUsjOipMiE85qUinZpeSRW7Su1CTJ11yLFe96mlqbqgvRFWJB6UOyzbG4EGPyPraxSUQic4lb1NSTEASji9m7iMFKq/pasHD4UjNhai2NQInkpUZu3Fua8+LEueRd6ZIrkRiFIIycJHURkrCXGmwNGVtNox5RRESQ5JixtZk6T6FgukJmS7lJO7bSTEFg8zGSVuaxZCk2h4BxScuVOZEuwUustkhylMvzlLKEWp2NFb2eJBc9Jp/io8UwFkxCnkoBe1mwKWgPaR5tZF2UjNWil5OUQ0wuVxsLlrGADSgSgxRqCrhutLQYQwzF1YixMbQ0l1NIFrUgR8wn1LUSg7cB9CJWtWDCPCOzJ+e9y96rjBXm0ddULC4KwWKlo5x3LiOzB4uNoqWJlBgRF0tVk4UKTy6JL4Ewr4xVjbTM1T++Bt7Von34ecw+upByCljvWrkgugOcdz608fq5i6BbUFJKL6aWEi8+RtQOqrp5zQXFIsTqEtYT0kqoNqaAXlIa6bEWpRRRYGJUywvlEvFY0TtJthQllTpwQPxIpJpXEri5r6r3iHc2Ys2jrqvZW+RQX5uE6uYO1uWT5GpT0tl1c5gnu06Kj9KbQ7HzEEZqbOUwDxTxLO9Ys9VhyLxkcSG4zAZvFJU5+zAEKTBOUYrKPMNmPG5FNjtFT+YpCZbQmDwiPijLPFa4m9QQnLThyjySn1E4a57KhWKjNNtfxBpN3aC2RfZZStL78OgKvggMO+qMi66WnBK2bov1krouM7jgXO6LSuYOpTJnH2PtW0jmwgq3vrDEvUrkMpfgOX2SpBhq74LAvxxrxJfKdRKwnZx4qTXbiSzW5VQqLdQSu0JTncdCmpKLIXkM7XUeS8zBFk6linmZNJeY/cApEzAlR1PmNcLAA5tVSmlpKaIQineuQhWK+RQj3MiWXJgN0my1ydsQuLyGs0T0PFyIss0lObG1pUUomnypyYlX9PLcu+piklBt0f1HXVujtRVN0aZWF95SM7ovVmPtlvhnq9aTdqO8+tuhjyVpfzzfLl9v+MkH5CAuE9+dLRAQrvzU0JT+gRgwDdD8ZeS6+7iRPlnw+zfbpQo2wQiCzb2zs33yU34mY/+zFpqBEPRfs3+70J8emfHxv57h93z+2qzOt2Yag5Z9zC/r/G+G065dZxe47l4uX635JSRLnN3f2xCS79SQ/KRXuxNRPy7cXiDdXnTx9kB8PPh1FmTcxebJYsOvOqmZ6k9LXNMPppVb3okc/QdJ7q/fXqwvV+0Hn0j7xT0iXm74ZRszu41q8GTxvbqUf2gt/SP+JULlmw6qOJh4i0jBfzlEonS786Wbx4ykFLhzIqi6OV+Gkpuzpo++EnS8c76sthK6WXwNJXET91z0Z33svDfaIao8dx1QT8YXO2jnotKMj/hg/vGcIChVqt4sTtY/HLjn3l9v8HZfnKz4URG26c4UuiN8V6mUulAWd+/jcY5TZgtydOgqBx7+oNJcxQlyVDcV9SnoThkP1+evl3gnqLvl/ie2Wlji6vxn/ozONS7zcLFdvv/QugDxPkodIPf6s2F2cff7xdm75Z3T9ebOvbuf35G7l28v7myXl9vhzjTjP4cgE30OoyB2QZ9k9kiSHkLCmA+m6uNK6EN+a6hdu7aJ57L28Wb1enW+uxEv6iJsbUE+rPlGxO8/SqN/v/v0X45Gu19jUobHn1/BFPLy7Jfjcb+gD74uSZkNEsO8pBxQ/TB5T/xN5iGjMfNegobnN/7mQ5wjtHlVWZTtTc7l3PYiXeaqoUM79uZDzmhVrqDxsXRfTLyueuG+v3DrnhLsbcf5EuqqQ2/ItKl8tDPFigSD0pJT8fo7c4d+fAzs72SK00x1drFVVjHxvaPB31Pmcaaofis7JrCPy77BFGGm/wxM8TYehAzoNsSYaPPfxRBjsC6iOBSuWOrVzc4+/8sNHvkpzObKb2beCLH8l/4xzA+G2ahJKRbno8cchlGt8cIPBNpgmwyYjHKyuTqvst5vFmgTQpgn52st+BZVZZz7n8d01s+rJRc5VINBp1ycvIKXPK9SosuxBFj1rxhmE1L+3YXZ8EQGj/ikHLONMejkToE0xNkcXSxTAaKA5oRIcfMe0B7UFPBzQm2OReLcLtgm3AjE+S2CbY7h+0uG23DKHAuv0WuEFEtx3nK3UBBsfubvmqqK8+l3Ta/9kOk18CO/a3pdS7v176J2+evVX4nCbp6/s2ECP8XXHP+xbqTnW2uVHGEHV5vsouPmD3IOfyJ6H09zdIcdCuL08fsKqHGXby8urwZ1/WHx7vJytTj/QEiNPtZ0QPVPMTXNFqW2n4N3HRTehcz/ts9RHQRrfIqp+W1eDVPdD+15Z576FFPTXt34V4lOQ2i9Fpumab952NPfHVOj0h3BMG426C3drAXVvJ8NyXITvktws0G4zJurEHnDmeaGu80vEiTDuzq8NRtw9D9w3zLW7JIKefpH5j14pMc89PAFhXqem2v8gmt5/SM0SHO0jm8N5elDWyit6i4mooW99CgZfAeS0RAcmU8RO9wiEmpDxAmfxHYk08J5WspUIkzBNjebqw2JMo2GwBcaaz323DwNsjXYC8eW2kfjGxTaMOJEq04y/dhRpxFJo1OgrrSqneiapARvMwL8e/cxnxaaXsxfXWi7pObOwY283q3rxblquqF911RNaE5Xk0Nx9z/uH1pi8gJSgN83Vj+mlnVYorsBNY+e2HKaK0/W0jwgS93ug521peYFpumxZXcvsNIbOGwntM5Dw4ifzEYu6ThcQay31Ut2h6zJY6j10x3cnO7HTqZDgrRmpZG00a99VxoraRvtgX+FWwgcMlGmbuGDeev48yYuX9HRfq64/P/xwLyVmvLzfTA/haDfLrb5k7j8SVyeHoWDeagW95FXAj6FoKvSc6D3To85/P7EZddC0lVCJgT9UFxWF7Gr4rJe+e5iy6Ukm+b4hPxW0nJt0rKozNZEZseH/tEE/aYCXCZj/61qnSb4yVyDlWUekQAPhWnXknz/ILuHSWvR0MTILlum1szuQ0umuQbo5tZTl1s7VJqYWVu90hrrH71Ir94/ytSmSq+xQbH1ECZxV/vrUr8GXHdBuAv/0hJjG78Oo+sKKuIqhZQ4iRIt4phvTVDpVN5T/CDhdy8T75dTYFyBn0flQ/9owl5uU49bleD02y7UbRKCuxd5d6dXqbDHEvTQA9+yuv9+E0NbUmjCZZc/W5DD7qOVTE1C7uJrF3c71KXe2jBo4aDTRy/Sq/eP0tvs3XZcWg+TdN6hNo4W1+CvjEoqS6rL+W0ckpqKoFRSAqkIrI61+me3pKB0c/2liU56Pn4FmXjnBHcLofi/yftQyjyHWivPIupPgXRrvAascMNMIC8vJjavQTzv0txRoxKKoK804p+j8TI2J4JtkysVbr5z0JErPjlX/A/1Dvhfyv+QC/G9m8M/j/+hw//wcueA6LJkx3vDvxt/m9/CAfEakQgbCzzR/Lsl0icPRDwQYVlCjF12cedhDcvyEpN3RDZEfltqckH08xwJcasuB/Wqm1hckETcD08HFKvuhJOnoQ94xPDQTiRs/cNOiOqtM3kaEvUjrvLKUIkBKX7HCz85IR51BLh+5/93eGbfxob0L+SESKTH01eb1QWxGapmf7l6/eZs9frN9v76/Hz5arv/UYIHq/fLk3andro4u1z+9P8ApW6eDA==').then(json => {\n",
       "   const obj = Core.parse(json);\n",
       "   Core.draw('root_plot_1779222335799', obj, '');\n",
       "});\n",
       "\n",
       "      }\n",
       "      const servers = ['/static/', 'https://root.cern/js/7.11.0/', 'https://jsroot.gsi.de/7.11.0/'],\n",
       "            path = 'build/jsroot';\n",
       "      if (typeof JSROOT !== 'undefined')\n",
       "         execCode(JSROOT);\n",
       "      else if (typeof requirejs !== 'undefined') {\n",
       "         servers.forEach((s,i) => { servers[i] = s + path; });\n",
       "         requirejs.config({ paths: { 'jsroot' : servers } })(['jsroot'],  execCode);\n",
       "      } else {\n",
       "         const config = document.getElementById('jupyter-config-data');\n",
       "         if (config)\n",
       "            servers[0] = (JSON.parse(config.innerHTML || '{}')?.baseUrl || '/') + 'static/';\n",
       "         else\n",
       "            servers.shift();\n",
       "         function loadJsroot() {\n",
       "            return !servers.length ? 0 : import(servers.shift() + path + '.js').catch(loadJsroot).then(() => execCode(JSROOT));\n",
       "         }\n",
       "         loadJsroot();\n",
       "      }\n",
       "   }\n",
       "   process_root_plot_1779222335799();\n",
       "</script>\n"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "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
}