{
"cells": [
{
"cell_type": "markdown",
"id": "9fc7417d",
"metadata": {},
"source": [
"# gr011_graph2d_errorsfit\n",
"\n",
"\n",
"\n",
"**Author:** Olivier Couet, Jamie Gooding \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:37 PM."
]
},
{
"cell_type": "code",
"execution_count": 1,
"id": "03ed2e45",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2026-05-19T20:37:57.647247Z",
"iopub.status.busy": "2026-05-19T20:37:57.647111Z",
"iopub.status.idle": "2026-05-19T20:37:58.705760Z",
"shell.execute_reply": "2026-05-19T20:37:58.704956Z"
}
},
"outputs": [],
"source": [
"from ctypes import c_double\n",
"\n",
"import ROOT\n",
"\n",
"c1 = ROOT.TCanvas(\"c1\")\n",
"\n",
"e = 0.3\n",
"nd = 500"
]
},
{
"cell_type": "markdown",
"id": "9c916766",
"metadata": {},
"source": [
"To generate some random data to put into the graph"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "91f962d2",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2026-05-19T20:37:58.711778Z",
"iopub.status.busy": "2026-05-19T20:37:58.711616Z",
"iopub.status.idle": "2026-05-19T20:37:58.863654Z",
"shell.execute_reply": "2026-05-19T20:37:58.863161Z"
}
},
"outputs": [],
"source": [
"r = ROOT.TRandom()\n",
"f2 = ROOT.TF2(\"f2\", \"1000*(([0]*sin(x)/x)*([1]*sin(y)/y))+200\", -6, 6, -6, 6)\n",
"f2.SetParameters(1, 1)\n",
"\n",
"dte = ROOT.TGraph2DErrors(nd)"
]
},
{
"cell_type": "markdown",
"id": "f8a0403f",
"metadata": {},
"source": [
"Fill the 2D graph. It was created only specifying the number of points, so all\n",
"elements are empty. We now \"fill\" the values and errors with SetPoint and SetPointError.\n",
"Note that the first point has index zero"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "af6cfefb",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2026-05-19T20:37:58.866465Z",
"iopub.status.busy": "2026-05-19T20:37:58.866335Z",
"iopub.status.idle": "2026-05-19T20:37:58.986462Z",
"shell.execute_reply": "2026-05-19T20:37:58.985467Z"
}
},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"Warning in : function:f2 has 24 negative values: abs assumed\n"
]
}
],
"source": [
"x = c_double()\n",
"y = c_double()\n",
"zmax = 0\n",
"for i in range(nd):\n",
" f2.GetRandom2(x, y)\n",
" rnd = r.Uniform(-e, e) # Generate a random number in [-e,e]\n",
" z = f2.Eval(x, y) * (1 + rnd)\n",
" if z > zmax:\n",
" zmax = z\n",
" dte.SetPoint(i, x, y, z)\n",
" ex = 0.05 * r.Rndm()\n",
" ey = 0.05 * r.Rndm()\n",
" ez = abs(z * rnd)\n",
" dte.SetPointError(i, ex, ey, ez)"
]
},
{
"cell_type": "markdown",
"id": "143306cc",
"metadata": {},
"source": [
"If the fit is not needed, just draw dte here and skip the lines below\n",
"dte.Draw(\"A p0\")"
]
},
{
"cell_type": "markdown",
"id": "5fa0e206",
"metadata": {},
"source": [
"To do the fit we use a function, in this example the same f2 from above"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "bd196b6c",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2026-05-19T20:37:58.988074Z",
"iopub.status.busy": "2026-05-19T20:37:58.987945Z",
"iopub.status.idle": "2026-05-19T20:37:59.150754Z",
"shell.execute_reply": "2026-05-19T20:37:59.149761Z"
}
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"****************************************\n",
"Minimizer is Minuit2 / Migrad\n",
"Chi2 = 515.752\n",
"NDf = 498\n",
"Edm = 8.12059e-06\n",
"NCalls = 35\n",
"p0 = 0.576226 +/- 0.19029 \n",
"p1 = 1.72868 +/- 0.57087 \n"
]
}
],
"source": [
"f2.SetParameters(0.5, 1.5)\n",
"dte.Fit(f2)\n",
"fit2 = dte.FindObject(\"f2\")\n",
"fit2.SetTitle(\"Minuit fit result on the Graph2DErrors points\")\n",
"fit2.SetMaximum(zmax)\n",
"ROOT.gStyle.SetHistTopMargin(0)\n",
"fit2.SetLineColor(1)\n",
"fit2.SetLineWidth(1)\n",
"fit2.Draw(\"surf1\")\n",
"dte.Draw(\"same p0\")"
]
},
{
"cell_type": "markdown",
"id": "5a168a17",
"metadata": {},
"source": [
"Draw all canvases "
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "af6b67ce",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2026-05-19T20:37:59.152262Z",
"iopub.status.busy": "2026-05-19T20:37:59.152131Z",
"iopub.status.idle": "2026-05-19T20:37:59.400558Z",
"shell.execute_reply": "2026-05-19T20:37:59.399616Z"
}
},
"outputs": [
{
"data": {
"text/html": [
"\n",
"\n",
"\n",
"
\n",
"\n",
"\n",
"\n"
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"from ROOT import gROOT \n",
"gROOT.GetListOfCanvases().Draw()"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.12.12"
}
},
"nbformat": 4,
"nbformat_minor": 5
}