{ "cells": [ { "cell_type": "code", "execution_count": 14, "metadata": {}, "outputs": [], "source": [ "import ROOT\n", "import array\n", "import matplotlib.pyplot as plt\n", "import numpy as np\n", "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "--2019-04-01 07:35:38-- https://www.hepdata.net/record/resource/406719?view=true\n", "Resolving www.hepdata.net (www.hepdata.net)... 188.184.64.140\n", "Connecting to www.hepdata.net (www.hepdata.net)|188.184.64.140|:443... connected.\n", "HTTP request sent, awaiting response... 200 OK\n", "Length: 953915 (932K) [text/xml]\n", "Saving to: ‘bdt.xml’\n", "\n", "bdt.xml 100%[===================>] 931.56K 605KB/s in 1.5s \n", "\n", "2019-04-01 07:35:40 (605 KB/s) - ‘bdt.xml’ saved [953915/953915]\n", "\n" ] } ], "source": [ "!wget -O bdt.xml https://www.hepdata.net/record/resource/406719?view=true" ] }, { "cell_type": "code", "execution_count": 15, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 15, "metadata": {}, "output_type": "execute_result" }, { "name": "stdout", "output_type": "stream", "text": [ " : Booking \"BDT method\" of type \"BDT\" from bdt.xml.\n", " : Reading weight file: bdt.xml\n", "
DataSetInfo : [Default] : Added class \"Signal\"\n", "
DataSetInfo : [Default] : Added class \"Background\"\n", " : Booked classifier \"BDT\" of type: \"BDT\"\n" ] } ], "source": [ "ROOT.TMVA.Tools.Instance()\n", "reader = ROOT.TMVA.Reader()\n", "\n", "var = [array.array('f',[0]) for i in range(7)]\n", "reader.AddVariable(\"MET\",var[0]);\n", "reader.AddVariable(\"MT\",var[1]);\n", "reader.AddVariable(\"dMT200\",var[2]);\n", "reader.AddVariable(\"m_tophad\",var[3]);\n", "reader.AddVariable(\"m_toplep200\",var[4]);\n", "reader.AddVariable(\"dphi_lep_nu200\",var[5]);\n", "reader.AddVariable(\"dphi_rjet_lep\",var[6]);\n", "\n", "reader.BookMVA(\"BDT method\", \"bdt.xml\")" ] }, { "cell_type": "code", "execution_count": 16, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 16, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAD8CAYAAABzTgP2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzt3Xd4nPWV8P3v0aj3akuWZEu2hSX3XjAkFGdjO8WwaRBagISXTdiUfd5N2E1299nkybvZTZY8IRBYAiSQsJBcKcQhpsWGJeAq9yLZyLKaJVu9WF0zv/ePGdmSPLJHmtHcU87nunQhzX3fM2d0mzn6tfMTYwxKKaXUsAirA1BKKRVYNDEopZQaRRODUkqpUTQxKKWUGkUTg1JKqVE0MSillBpFE4NSSqlRNDEopZQaRRODUkqpUSKtDmAyMjMzTUFBgdVhKKVUUNm/f3+zMSbraucFZWIoKCigtLTU6jCUUiqoiEi1J+dpV5JSSqlRNDEopZQaRRODUkqpUTQxKKWUGkUTg1JKqVE0MSillBpFE4NSSqlRgnIdg1KBoqalh98eqMPdFrkJMZHcu76Q6Ej9+0sFF00MSnnhu9tO8Prx84iMfnw4T8RH27hrXYHf41LKG5oYlJqks+29vHniPH9zwxy+sbF41DFjDJ98chdPvH2az6yaqa0GFVT0X6tSk/Tfe5zVBe5YM/OyYyLC3940l/qOPn53oM7foSnlFU0MSk1C36CdF/fWsqFkOnlp8W7P+eA1WSzOS+Enb59myO7wc4RKTZ4mBqUmYdvRBlq7B7jn2oJxz3G2Goqoae3hD4fq/RecUl7SxKDUJDy3q5o5WQlcOyfjiudtKJlGSU4yj79Vgd1x+cwl5R97Kluoa+uxOoyg4ZPEICIbReSkiFSIyMNujouIPOo6fkRElo849jUROS4ix0TkRRGJ9UVMSk2VQ7XtHK5t555rC5Cx05HGGB5rqGzu5pUj2mqwwpDdwX0/38dXXzrkdlqxupzXiUFEbMDjwCZgPnC7iMwfc9omoMj19QDwhOvaXODLwEpjzELABtzmbUxKTaXnd1WREG3j1mW5Hp2/cUE2RdMSefytChzaavC79xsv0D1gp7S6jV2VLVaHExR80WJYDVQYYyqNMQPAS8CWMedsAZ43TruBVBHJcR2LBOJEJBKIB/TPKhWwWi7088rhBj6xIo+k2CiPromIEB66aS6nzl/g9ePnpjhCNdah2nYAEmMi+fH2CoujCQ6+SAy5QO2In+tcj131HGPMWeAHQA3QAHQYY97wQUxKTYmX9tUyYHdw97pZE7ruo4tnUJiZwI93VGh3hp8drm0nJS6Kr9xcxK7KFkqrWq0OKeD5IjG462Qd+y/f7TkikoazNVEIzAASROROty8i8oCIlIpIaVNTk1cBKzUZQ3YHL+yuZv3cDOZOS5rQtbYI4Ys3zOFEQyfbyxqnKELlzqHadpbkp3LH2pmkJ0Tz6A5tNVyNLxJDHZA/4uc8Lu8OGu+cDcAZY0yTMWYQ+B1wrbsXMcY8ZYxZaYxZmZV11b2slfK5P5c1Ut/Rx92TLHFxy7Jc8tLi+PGO97XV4Cc9A0OcOt/F0vxU4qMj+cL1s3nnVNPF7iXlni8Swz6gSEQKRSQa5+Dx1jHnbAXuds1OWouzy6gBZxfSWhGJF+f0jpuBMh/EpNSk7K9uY+fpZvqH7Jcde35XFbmpcdxcPG1Szx1li+CLN8zlcF0Hb5/UVq8/HK3rwGFgaX4KAHetm0VqfBSP7Xjf4sgCm9e1kowxQyLyEPA6zllFzxpjjovIg67jTwLbgM1ABdAD3Os6tkdEfgMcAIaAg8BT3sak1GQM2h3c9cweegbsxEZFsKognfVzM1k/J5OYqAh2nm7h6xvnEWmb/N9Tn1yRx1PvnObbr5xg3ZwMYqNsPnwHaqzDdc6WweK8VMA5AH3f+kIeefMUx+s7WDAjxe11rd0D1Lf3sjDX/fFQ55MiesaYbTg//Ec+9uSI7w3wpXGu/RfgX3wRh1LeKGvopGfAzheuL2TQbth5upnvvVoOOMcIoiMj+MzK/Ks8y5VFR0bwnVsWctcze3nyf07z1Q3X+CJ0NY7DtR3kpcWRmRhz8bF7ri3gp+9U8tiOCp64c8Vl1+w83cxXXjpER+8gR/7lr8IyeWt1VaVcDlS3AXDv+kJmpMYB0NjZx87TLbxX0UxJTjIZIz5gJuv6oiw+vmQGP3nrNFuW5lKYmeD1cyr3DtW2s3Rm6qjHUuKi+Nz6An68o4JT57u4ZrpzIoHdYXhsRwU/2n6KSFsEA0MOWroHyHX9WwgnWhJDKZeDte1MT44hJ+XS4vtpybHcsiyX739qCfddV+iz1/rWR0uIiYrgn14+pgPRU6Spq5+z7b0sy0+97Nh96wtJiLbxmGuGUlNXP3c/u4cf/vkUW5bm8v1PLgaguavfrzEHCk0MSrkcqGlj+cy0q5a58IVpSbF8/cPzeLeima2HdU3nVDjsmnm0xE1iSEuI5s51s3jlSD0v7a1h86N/obSqjf/4xGIe+fQSZqY7K+a2dGtiUCpsNXX1U9vay7KZl3+ITJXPrpnFkrwUvvNKGR29g3573XBxuK4dW4SwYEay2+NfuH420ZERPPy7oyTHRvKHh9bz6VX5iMjFMYnmrgF/hhwwNDEoBRyscY4vLJ+Z5rfXtEUI3711Ea3d/fzg9ZN+e91wcai2nWumJxEf7X4oNTMxhm99ZD6fu7aArQ9dR3F28qhjAE0XwrPFoIPPSgEHatqJsonfpycuzE3hnmsL+PnOKj6xIo+lbro91MQZYzhc285HFs+44nl3rnVf2iQu2kZCtI3mME0M2mJQCmeLYX5OsiVTE//uQ9cwLSmGb/7+qO705iNVLT109g1dXNg2GZlJMbRc0K4kpcLSkN3BkboOlvmxG2mkpNgo/vfHFnC8vpMX9tRYEkOoOVTr7Bp0N/DsqczEGG0xKBWuys910Tto9+vA81gbF2azYlYaL+7VxOALh2s7iI+2UTTBYocjZSZGa2JQKlxZMfA8lojwscU5lJ/r4nTTBcviCBWHattZlJuCLWLyU48zEmNo1q4kpcLTgZp2spJiyEuzdoXrxoXOvateO6ab+XhjYMjBifpOrwfyMxNjaOsZCMtxH00MKuwdrGljWX6qXxa2XUl2SiwrZqWx7WiDpXEEu7KGTgbsDq/GFwCyEqMxBlp7wq/VoIlBhbWWC/1UtfSwfJZ13UgjbVqYzfH6Tqpbuq0OJWgNV1T1NjGE8yI3TQwqrA1v2OKuno4VNi7MBuBV7U6atEO17WQmxjBjRM2ryRgumBiOA9CaGFRYO1DTRmSEXKzXb7W8tHiW5KXwqnYnTdrh2naW+qBrMDMxGtDEoFTYOVDdTklOMnHRgVNzf9OiHA7XdVDX1mN1KEGns2+Q003dXi1sG5aZ5GwxhOMiN00MKmzZHYbDde2Wrl9wZ5OrO0lnJ03ckdoOwPvxBYCkmEiiIyO0xaBUODl5roueAbul6xfcmZWRwIIZyTo7aRLGbuXpDREhMyE6LAvpaWJQYeugq2xCoLUYADYvyuFATTvnOvqsDiWoHKptZ3ZWAilxUT55vsyk8FzkpolBha0D1e1kJERf3JQlkFzqTtJWg6eMMc6tPH04kSAzMYYWbTEoFT4O1rSxzE87tk3U7KxE5k1PYpuOM3isvqOPpq5+n4wvDAvXekmaGFRYauseoLK5OyC7kYZtWpTNvqpWGru0O8kTeypbAFhVkO6z53S2GAZwOMJrX25NDCosDS9sC7SB55E2L8rBGHj9+HmrQwkKuytbSImLojh78hVVx8pIjGHIYcJu61VNDCosHaxpI0JgcZ5/d2ybiKJpiczJStDFbh7ac6aVNYXpRHhRUXWs4UVuLd3h1Z2kiUGFpQM17RRnJ5MQE7i724oImxflsLuyJSwHQCeivr2X6pYe1s7O8OnzZg3v/Rxm9ZI0MaiwY3c4Z68snxW44wvDNi3MwWHgjRPanXQle844xxd8nRiGVz+H2wC0JgYVdsoaOrnQPxTQ4wvDSnKSKMiI18VuV7H7dKvPxxcAMhLCs16SJgYVdnaUNyIC1xdlWR3KVYkI1xdlcbCmHWPCa2bMROw+08JqH48vAKTFR2OLkLCrl+STxCAiG0XkpIhUiMjDbo6LiDzqOn5ERJaPOJYqIr8RkXIRKRORdb6ISanxbC9vZEleKlmuboJAV5KTzIX+Ieraeq0OJSBN1fgCQESEkJ4QfmsZvE4MImIDHgc2AfOB20Vk/pjTNgFFrq8HgCdGHPsR8JoxphhYApR5G5NS42nq6udIXTs3F0+zOhSPleQ4u0fKGjotjiQwXRpf8N36hZEyE2M0MUzCaqDCGFNpjBkAXgK2jDlnC/C8cdoNpIpIjogkAx8AngEwxgwYY9p9EJNSbr19shFj4KaS4EkM87KTEIGyhi6rQwlIw+MLJdnJU/L8mYnRNGlX0oTlArUjfq5zPebJObOBJuBnInJQRJ4WkQQfxKSUWzvKG8lOjmV+ztR8iEyF+OhICjISKD+nLQZ39kzR+MKwzMQYmru0xTBR7u7G2FGy8c6JBJYDTxhjlgHdwGVjFAAi8oCIlIpIaVNTkzfxqjA1MOTgnVNN3Fg8LSDrI11JcXaSdiW50dDRS9UUjS8My0yMpqW7P6wG/32RGOqA/BE/5wH1Hp5TB9QZY/a4Hv8NzkRxGWPMU8aYlcaYlVlZgT+bRAWevWda6R6wB9X4wrCSnGSqW3vo7h+yOpSAsqeyFZi68QVwthj6Bh10D9in7DUCjS8Swz6gSEQKRSQauA3YOuacrcDdrtlJa4EOY0yDMeYcUCsi81zn3Qyc8EFMSl1me/l5YiIjWD830+pQJqw4Owlj4OR5HWcYaXdlC8mxkRRP0fgCOOslAWHVneR1PQBjzJCIPAS8DtiAZ40xx0XkQdfxJ4FtwGagAugB7h3xFH8LvOBKKpVjjinlE8YYtpc1cu2cjIDa39lTJa4xkfKGrqBYmOcvuytbWF2YgW2KxhfgUr2k5gv9FGSGxxCoTwrFGGO24fzwH/nYkyO+N8CXxrn2ELDSF3EoNZ7TTd3UtPbwhQ/MtjqUSclLiyMpJlLHGUYYHl+4c+2sKX2dzOEWwyRnJj3z7hlmpcezYf50X4Y1pXTlswoLO8qdtYZuCsLxBXCugC7OSdKZSSNcGl+YuoFn4OJCyMmsZRi0O/iP18r55Z5qX4c1pTQxqLCwvayR4uwkclPjrA5l0oqzkylv6Aqr2TFXMjy+UDLFU4/TvaiXdKK+k/4hB7WtPb4Oa0ppYlAhr6NnkNLqNm4OokVt7pTkJNOlpTEu2nOmdcrHFwCibBGkxkdNKjHsr24DoK6tN6gSuiYGFfL+5/0m7A7DTcXB08frTrGWxrjoXEcfZ5q7p3Sa6kjDW3xO1P4aZ2LoH3LQFERlNTQxqJC3o+w86QnRLPXhJvFWmDfdWRqj/JxOWZ2q/RfGk5k4uUJ6B6vbSI2PAqC2NXhaepoYVEgbsjt4+1QTN8zLmvIuh6mWEBPJrPR4bTHgHF9I8sP4wjBnIb2JtRjq23up7+hj86IcAOragmecQRODCmkHa9tp7xnk5iDvRhpWkpOsLQZgd6Vzf2d/JfvJ1Es64OpG2rJkBkBQjQ1pYlAhbXtZI5ERwvXXBN9qZ3eKs5OpaummZyB8S2NcGl/wTzcSOLuSuvqH6Bv0vCzG/uo2YqMiWD4rjczE6KCamRS4O6Er5QM7ys+zujCd5Ngoq0PxiZIcZ2mM8nPWr4AeGHIQHem/vy1P1Hfym/11/OHQWQC/ljYZXuTW0j3g8ZTnA9VtLMlLJcoWQV5aPLXalaSU9Wpbezh1/kLQLmpzZ2RpDCv98M1TrPjOm5ya4tpNzRf6eebdM2z60V/Y/Ohf+OXualYXpvOL+1f7bXwBRqx+9rA7qW/QzvH6TpbPcibv/PT4oBp81haDCllbDzuL/G4oCY3xBXCWxki0uDTGi3tr+NH29wH41svH+NUDa6ekjPkz757h37aVMeQwLMlL4dtbFvCxxTNIcy0486eMxIktcjtS18GQw7DC1arLS4vj1aMN2B0mKCZBaGJQIal/yM7Pd1ZxfVFmSBU+ExGKs60rjfHWyUa+9fIxPnhNFhtKpvFPfzjOy4fOcuuyPJ++znsVzXz3Tye4cd40Ht5UTNH0JJ8+/0Rd7ErycGbS8MK2iy2GtHiGHIZznX1Bsfpeu5JUSHrlcANNXf184frgLJp3JSU51pTGOFrXwZdeOEBxdhKP37GcO9bMYml+Kt/9UxkdvYM+e51zHX18+cWDzM5K5NHbl1meFOBSvSRPF6ntr25jdmbCxXIa+enOZBAsA9CaGFTIMcbw9LtnuGZ6ItcXhcZspJGKc5L8XhqjtrWHe3++j7T4aH72uVUkxkQSESH8n1sW0to9wCNvnPTJ6wzaHTz03wfoHbTz5J3LSYgJjE6N2CgbiTGRHnUlGWM4WNPGshGTA/LS4oHgmbKqiUGFnJ2nWyhr6OTz180Oui08PTE86OqvcYb2ngE+97O9DAzZee6+VUxLjr14bGFuCnetncUvdldz7GyH16/176+WU1rdxvc+sZi506xvKYyUkRjt0SK36pYeWroHWDHrUmKYkRqLiLYYlLLM03+pJDMxmo8vnWF1KFPCn6Ux+gbtPPD8fmpbe/np3Svdflj/3V/NIz0hhm++fAyHY/LdW68ebeDpd89wz7pZfHxJ4N07Txe5DY8vjEwMMZE2spNjg2bKqiYGFVIqGrt462QTd60tIDYq+HZq84Q/S2N8//WT7K1q5QefXsKacRaUpcRF8c2PFHO4tp2X9tVO6nUqmy7w9785wtL8VL75kfnehDxlMhOjaen2IDHUtJEUE0nRtMRRj+elxWlXklJWeObdKmIiI7hz7UyrQ5lS/iiN0dk3yIt7a/jr5blX/Qv+lqW5rClM599fK6dlgsXmegfsfPGFA0TZhMfvWO7XRXMT4Wm9pAPVbSydmUrEmGmp+Wnx1GlXklL+1XKhn98dqOOvl+dd3MA9VPmjNMav99XSM2DnvvWFVz1XRPjOLQvp7h/i318r9/g17A7D1397hJPnu/jhZ5YG9FTOjMQY2noGGLI7xj2nq2+Qk+e7RnUjDctLj6ehs4+BofGvDxSaGFTIeGFPDf1DDu6/rsDqUKbccGmMk1PUarA7DM/vqmZVQRoLc1M8uuaa6Uncf10hvy6t409HGq56vsNh+MZvj/DHw/V8/cPF3DAvsFeoZyVGYwy0do/fajhU244xuE8MaXEY49yrOtBpYlAhoW/QzvO7qrhxXlbAzWaZCpdmJk1NYnirvJGa1h4+d+3VWwsjffnmIlbMSuOhFw/wX/9zety1Fg6H4R9+d5Tf7K/jqxuK+Jsb5vgi7Cl1sSzGFbqT9le3IYLbvT/yXVNWg6E0hiYGFRK2Hqqn+cIAnw/BBW3uDJfGmKoV0D/beYaclFg+vGBi5UQSYiJ54fNr2Lwoh397tZxvvnzssq4Xh8PwzZeP8avSWr5801y+uuEaX4Y+ZTKThhPD+GMo+6vbmDc9iSQ3RRsvLnILgplJmhhU0HMuaKukJCeZa+f4rxSzlYZLY0zFzKRT57t4r6KFu9bNItI28Y+I2CgbP75tGV+8YQ7/vaeG+54rpavPuTLaGMM/bz3Gi3tr+OINc/jah4IjKQBkJFy5XpLDYThU036xDMZY2cmx2CIkKDbs0cSggt5f3m/m1PkLfP66wpBc0DaeqSqN8bP3nDO7bl81+ZldERHC1zcW8++fWMTOimY+9eQuzrb38r+3HueXu2v4fz44m7//8Lygul9XazG833iBrv6hi4Xzxoq0RTAjNVa7kpTyhxf2VJOZGMPHAnBR1FRanJdCV/8QJ3zYamjvGeD3B+u4dVmuT6qYfmbVTH5+72rOtvVy83++zXO7qvnC9YU8vLE4qJICQFJMJNGREeMW0nO3sG2s/CDZl0ETgwpq7T0DvFXexJalMwJ2/vtUGZ7Fs6Os0WfP+at9tfQNOvjc+gKfPed1RZn89ovXMis9gQc/OId/3FwSdEkBnN13WYkx4xbS21/dRkZCNLMy4sd9jmBZ5BYYFaqUmqRtR88xYHdw67Jcq0Pxu6ykGJbkp/Ln8kb+9uYir59vyO7g+V3VrJudQXG2bzfBuWZ6Eq9/7QM+fU4rZF6hXtJw4bwrJb38tHiauvrpG7QH9Mr88PoTS4Wclw+dZe60RBbM8N9uXoFkQ/E0Dte20zTBjerd+XPZec629/q0tRBqMtzUS7I7DI+8eYrK5m7WFKZf8fr89OEqq4HdneSTxCAiG0XkpIhUiMjDbo6LiDzqOn5ERJaPOW4TkYMi8oov4lHhoa6th71nWrl1WW5Qdk34wk0lzu6kt8q970569r0q8tLiQmrHO18bWy+psbOPO57ezaPb3+eTK/K4a92sK16flza8L0Ngdyd5nRhExAY8DmwC5gO3i8jYKlibgCLX1wPAE2OOfwUo8zYWFV7+cMi5dWcgVuL0l/k5yeSkxLK9/LxXz3O8voO9Z1q5Z11BUGw9aZXMxBhaLgzgcBjeq2hm86Pvcqi2ne9/cjE/+NSSq3YPhVOLYTVQYYypNMYMAC8BW8acswV43jjtBlJFJAdARPKAjwBP+yAWFSaMMfz+4FlWF6Rf/J8tHIkINxVP4y/vN9M3aJ/08zy3s4q4KBufXpnvw+hCT2ZiDEMOw3e3lXHnM3tIjY9i60PX8SkPf29ZiTFER0ZQG+AD0L5IDLnAyFq7da7HPD3n/wJfBwK/spQKGMfrO6lovMAtYTjoPNaGkun0DNjZXdkyqeubuvp5+WA9n1iRS0r85St21SUZic4pvM+8e4Zbl+Wy9aH1XDOBrUcjIoS81LiA37DHF4nBXbtz7Iobt+eIyEeBRmPM/qu+iMgDIlIqIqVNTU2TiVOFkJcPniXaFsFHFuVYHYrl1s3JIC7Kxo5JjjP8YlcVgw4H918XHuVEvLEsP405WQn8xycX88inlxIfPfGJnXnp8QE/ZdUXiaEOGNmOygPqPTxnPfBxEanC2QV1k4j80t2LGGOeMsasNMaszMrK8kHYKljZHYY/HK7nxuIs/QsXZwmK9XMz2V7WOOFV0L0Ddn6xu5oNJdMpzEyYoghDx8yMeLb/rxu86nLLT4sL+EVuvkgM+4AiESkUkWjgNmDrmHO2Ane7ZietBTqMMQ3GmH8wxuQZYwpc1+0wxtzpg5hUCNt5upmmrv6wXLswng0l0zjb3svJ8xOrtvqbA3W09QzyhTApPhgI8tPjae8ZvFg/KhB5nRiMMUPAQ8DrOGcW/doYc1xEHhSRB12nbQMqgQrgp8AXvX1dFb5+f/AsybGRAV+/359uKnb+LrZPYBW03WF49t0zLMlLYVXB+GUclG8NT1kN5O4kn6x8NsZsw/nhP/KxJ0d8b4AvXeU53gbe9kU8KnT1DAzx+rFzfHzpjIBeOepv05JjWZyXwvay83zpxrkeXfPnsvOcae7msc8uC9t1IFa4tC9Dz8V9NQKNrnxWQeXNE+fpHrCzZal2I411U/E0Dta2X3G/gJGe/kslualxbFyQPcWRqZGGp1cH8pRVTQwqqLx88CwzUmJZXXDl0gPhaEPJdIzxbBX0odp29lW1cd91hZPac0FNXlp8FPHRtoBe5Kb/IlTQaL7QzzvvN7NlWS4Rujr3MgtmJDM9Ocajaas//UslSbGRfGaVLmjzNxFxlt8O4LIYmhhU0HjlcD12h9HZSONwroKezjunmugfGn8VdG1rD68ebeCzq2eSGKMFlq2Qnx6nLQalfOHVY+cozk6a0ErTcLOhZBrdA3b2nmkd95xn3ztDhIhWUbVQXppzkZuvd9/zFU0MKmhUNnezOC/F6jAC2rVzMomJjBh32mpH7yC/3lfLx5bMICclzs/RqWF5aXFc6B+ivScw1zJoYlBBoXfATlNXPzPDuGCeJ+KibVw3N5Pt5efd/jX64t4augfsfP76QguiU8MuzUwKzO4k7WBUQWG4PzacK6l66qaSaWwvb+RDP3yHsWP0dW29XDsngwUztOVlpZGL3BbnpVoczeU0MaigUNOqicFTH108g/3VbfQOXD4AXTQtiQc/OMeCqNRIF1sMAVplVRODCgrDiUG7kq4uJS6KRz691Oow1BUkx0aREhelXUlKeaO2tZe4KBsZCdFWh6KUT+SlxVHV3ENH78QGoOOjbURN8aLEsEoMj+14nz8ebvDJc4nAfdcV6o5XflLT2sPM9Hit6aNCRkFGAn862sCSf31jQtf9/N5VU15AMqwSQ0ZijM9qzte09vCN3x4hJS6KD2utmSlX19aj4wsqpPy/H57H8lkTr2o7JytxCqIZLawSw+2rZ3L76pk+ea7eATu3/XQ3X37xIC8+sJblM7Vs8VQxxlDT2sO6ORlWh6KUzxRmJnD/dYE5bVjXMUxSXLSNZ+5ZyfTkWD7/XCnVLd1WhxSyWroH6Bmw68CzUn6iicELmYkx/PzeVTiM4XM/20dr94DVIYWk4Sl9w3XslVJTSxODl2ZnJfL03Ss5297LF54vpW9w/OJlanIuTlXN0MSglD9oYvCBlQXp/N/PLOVATRtf+9UhHI7ALIwVrLTFoJR/aWLwkc2Lcvjm5hJePXaOH/75lNXhhJTa1l4yE2OIi9atPJXyB00MPnT/dYVsKJnOS/tqA7acbjByrmHQSqBK+YsmBh8SEW4szqKpq5+qlsBc6h6Matt6dEaSUn6kicHH1hQ659rvqWyxOJLQMGh3UN/eq4vblPIjTQw+NicrgczE6CvuoKU8V9/ei8NoVVWl/EkTg4+JCKsL09mjicEnhjdM1xlJSvmPJoYpsLognbPtvQG92Xew0DUMSvmfJoYpsGa2c5xBu5O8V9vWQ5RNyE6OtToUpcKGJoYpMG96Esmxkeyp1MTgrZrWHnJT47CN3aNSKTVlNDFMgYgI5zjD3ipNDN6qbdVy20r5m08Sg4hsFJGTIlIhIg+7OS4i8qjr+BERWe56PF9E3hKRMhE5LiJf8UU8gWBNYQZnmrtp7OyzOpSgpolBKf/zOjGIiA14HNgEzAduF5H5Y07bBBS5vh50E+9rAAARu0lEQVQAnnA9PgT8L2NMCbAW+JKba4PS6sJ0AJ2d5IXOvkHaegZ1cZtSfuaLFsNqoMIYU2mMGQBeAraMOWcL8Lxx2g2kikiOMabBGHMAwBjTBZQBuT6IyXILZiSTEG1jzxld6DZZw8XzNDEo5V++SAy5QO2In+u4/MP9queISAGwDNjj7kVE5AERKRWR0qamJi9DnnqRtghWFKTrzCQv6BoGpazhi8TgbrrI2ApyVzxHRBKB3wJfNcZ0unsRY8xTxpiVxpiVWVlZkw7Wn9YUpnPq/AXdwGeStMWglDV8kRjqgPwRP+cB9Z6eIyJROJPCC8aY3/kgnoCxxjXOoK2Gyalt6yEpNpKU+CirQ1EqrPgiMewDikSkUESigduArWPO2Qrc7ZqdtBboMMY0iIgAzwBlxphHfBBLQFmcl0pMZISOM0ySs9y2thaU8rdIb5/AGDMkIg8BrwM24FljzHERedB1/ElgG7AZqAB6gHtdl68H7gKOisgh12P/aIzZ5m1cgSA6MoLlM9O0xTBJNa09zJueZHUYSoUdrxMDgOuDfNuYx54c8b0BvuTmundxP/4QMlYXpvPojvfp6B0kJU67RDzlcBjq2nrZUDLd6lCUCju68nmKrZmdjjGwv1pbDRPR2NXPwJBDF7cpZQFNDFNsWX4aUTbRhW4TVKMzkpSyjCaGKRYXbWNJXqoW1Jug4amq+Wm617NS/qaJwQ9WF6Zz7GwH3f1DVocSNGpaexCBXE0MSvmdJgY/WDM7gyGH4UBNm9WhBI3ath6yk2OJibRZHYpSYUcTgx+smJWGLUJ02uoEaFVVpayjicEPEmMiWTgjWccZJkAXtyllHU0MfnJdUSZ7q1r59h9P0DdotzqcgNY3aOd8Z78Wz1PKIj5Z4Kau7qEbi+jqG+LZ987w9qlGHvn0Upbmp1odVkCqa3NWVZ2ZoQPPSllBWwx+Ehdt49tbFvLL+9fQN2DnE0/s5D/fOMnAkMPq0ALOpamq2mJQygqaGPzsuqJMXvvaB7h1WS4/3lHBlsffo6zBbaXxsFXbpovblLKSJgYLJMdG8YNPLeGnd6+kqaufLY+9x9bDYyuVh6+alh5iIiPISoqxOhSlwpImBgt9aP503vjaB1g6M5Uvv3iQn713xuqQAkJtm3OqqrMqu1LK3zQxWCw9IZrn71vNhxdM51//eIL/eK0cZzHa8FXT2qvdSEpZSBNDAIiNsvGTO1Zw++qZ/OTt03zjt0cYsofnoLQxxrm4TUthKGUZna4aIGwRwv9360KmJcXwo+3v03JhgMc+u5y46PAqCXGmuZsL/UPMy062OhSlwpa2GAKIiPC1D13Dd25ZyI6Tjdzx9G46+watDsuvSqud9aRWFaRZHIlS4UsTQwC6a+0sHv/scg7WtvPz96qsDsevSqtaSYmLYk5WotWhKBW2NDEEqM2Lclg+M41Xj52zOhS/Kq1uY8WsNCIidEaSUlbRxBDANi3Mpqyhk6rmbqtD8YvW7gEqm7pZqd1ISllKE0MA27gwGyBsWg37XeMLK2elWxyJUuFNE0MAy0uLZ3FeCq8da7A6FL8orW4lyiYszkuxOhSlwpomhgC3cWE2h+s6ONvea3UoU660qo2FuSnERoXXFF2lAo0mhgC3aWEOAK+FeHdS36Cdo3UdrCrQbiSlrKaJIcAVZiZQnJ0U8t1Jx852MGB3sGKWDjwrZTVNDEFg48JsSqvbaOzqszqUKTO8sE0Tg1LW08QQBDYtzMEYeP34eatDmTKlVa3MzkwgM1FLbStlNZ8kBhHZKCInRaRCRB52c1xE5FHX8SMistzTaxVcMz2R2ZkJIdudZIxhv2thm1LKel4nBhGxAY8Dm4D5wO0iMn/MaZuAItfXA8ATE7g27IkIGxdms7uylbbuAavD8bnTTd209QzqwjalAoQvWgyrgQpjTKUxZgB4Cdgy5pwtwPPGaTeQKiI5Hl6rcHYn2R2GN0+EXnfS/upWAFbowjalAoIvEkMuUDvi5zrXY56c48m1CliYm0xeWhyvhmB3UmlVG2nxUczJSrA6FKUUvkkM7qqdjd2CbLxzPLnW+QQiD4hIqYiUNjU1TTDE4CcibFyQzbsVzSFXittZOC9dt/JUKkD4IjHUAfkjfs4Dxu5sP945nlwLgDHmKWPMSmPMyqysLK+DDkabFuUwaDfsKGu0OhSfab7Qz5lmLZynVCDxRWLYBxSJSKGIRAO3AVvHnLMVuNs1O2kt0GGMafDwWuWyLD+V6ckxIdWddKlwniYGpQKF11t7GmOGROQh4HXABjxrjDkuIg+6jj8JbAM2AxVAD3Dvla71NqZQFRHh7E76VWktPQNDxEcH/86s+6vbiI6MYJEWzlMqYPjkk8UYsw3nh//Ix54c8b0BvuTptWp8Gxfm8Nyuat4+2cTmRTlWh+O10qpWFuemEBOphfOUChS68jnIrC5MJzMxhl/ursaZb4NX36Cdo2c7WKHjC0oFFE0MQcYWITx04xx2nm7hz0E+CH2kroNBu2GVrl9QKqBoYghCd6ydxdxpiXz3TyfoH7JbHc6klV5c2KYtBqUCiSaGIBRli+CfPjqfqpYenttZZXU4k7a/qo05WQmkJURbHYpSagRNDEHqg9dkceO8LH68vYLmC/1WhzNhDoehtLpN93dWKgBpYghi3/rofHoH7fznGyetDmXC3jhxno7eQR14VioAaWIIYnOyErl7XQEv7avleH2H1eF4bEf5eb784kEW5abwkRCYcqtUqNHEEOS+cnMRqXFRfPuPJ4Ji+uqO8vM8+IsDzMtO4pf3ryEhJvgX6SkVajQxBLmU+Cj+7q/msedMK68dO+ez5+3uH8Lh8G2iGZsUUuKjfPr8Sinf0MQQAm5flc+86Ul8d1sZfYPeT1/t7Btk3b9t56+f2MmJ+k4fRKhJQalgookhBETaIvjnj82nrq2XZ9494/Xz7TvTSmffECfPdfGxx97l37aV0TMwNO75fYN2th1t4CdvV/DywbPsqWyhtrWHQbsD0KSgVLDRDt4QsX5uJjfMy+L5XVV88YY5Xu1tsOt0C9GREbz99zfwwzdP8V/vVPKnow1855aF3DhvGuCcbrq3qpWXD57lT0cb6Oq7PHGIwPSkWFq6+ynOTtakoFSQ0MQQQj40fzpvn2yiqqWHwszJ74a2q7KFFTPTmJ4cy/c+sZhbl+Xyj78/yr0/28dHFudQkBHPywfrOdveS3y0jY0Ls/nrZXksyU/hfGc/DR291Lf3cra9j4b2XiJtwsMbSzQpKBUkNDGEkLWzMwDnX/yTTQztPQOcaOjkaxuuufjYmtkZbPvK9fzX/1Ty2FsVDNkdXF+Uxdc3zuND86ePKv+dFBvF3GmJ3r0RpZSlNDGEkNmZCUxLimFXZQufXTNzUs+x50wrxlxKMsNiIm18+eaii8+bmRjjdbxKqcCkiSGEiAjr5mSw83QLxphJjTPsOt1CbFQES/Ldb5yjCUGp0KezkkLM2tkZNHX1c7qpe1LX765sYeWsdN04R6kwpokhxKxzdQHtrmyZ8LUtF/opP9fFujkZVz9ZKRWyNDGEmFkZ8WQnx7JrEolhzxnn/ghjxxeUUuFFE0OIGR5n2FPZMuHaSbtOtxAfbWNxnvvxBaVUeNDEEILWzk6n+cIAFY0XJnTdrsoWVhWkE2XTfxZKhTP9BAhB62ZnAkyoO6mxq4+Kxgs6vqCU0sQQivLT48hNjWPXac8Tw+5K5/jCOh1fUCrsaWIIQSLCmtnp7DnT6nHp7F2nW0iKiWTBjOQpjk4pFeg0MYSodbMzaO0e4FRjl0fn765sYXVhOpE6vqBU2NNPgRA1POV0twfdSec6+jjT3K3jC0opQBNDyMpPjycvLc6jAehdlc2Arl9QSjlpYghh62ZneDTOsOt0C8mxkZTk6PiCUsrLxCAi6SLypoi87/pv2jjnbRSRkyJSISIPj3j8+yJSLiJHROT3IpLqTTxqtLWzM2jvGaT83JXHGXZVtrBmdga2iMlv7qOUCh3ethgeBrYbY4qA7a6fRxERG/A4sAmYD9wuIvNdh98EFhpjFgOngH/wMh41wlrXmMGVupPq2nqobe3VaapKqYu8TQxbgOdc3z8H3OLmnNVAhTGm0hgzALzkug5jzBvGmOE9IXcDeV7Go0bITY1jZnr8FdczDB/TgWel1DBvE8N0Y0wDgOu/09yckwvUjvi5zvXYWPcBr3oZjxpj3ewM9p5pwT7OOMOuyhbS4qOYNz3Jz5EppQLVVRODiPxZRI65+dri4Wu467ge9SklIt8EhoAXrhDHAyJSKiKlTU1NHr60Wjcng86+IcoaOi87Zoxh9+kW1s7OIELHF5RSLlfdwc0Ys2G8YyJyXkRyjDENIpIDNLo5rQ7IH/FzHlA/4jnuAT4K3GyuUA7UGPMU8BTAypUrJ1Y2NIytHbE/w8LcS1VTGzp6efXoOeo7+njwBu1GUkpd4u3WnluBe4Dvuf77Bzfn7AOKRKQQOAvcBnwWnLOVgG8AHzTG9HgZi3IjOyWWwswE3q1oZtnMVHaUN7K9rPHiTKU5WQlsXJBtcZRKqUDibWL4HvBrEbkfqAE+BSAiM4CnjTGbjTFDIvIQ8DpgA541xhx3Xf8YEAO86dqfeLcx5kEvY1JjrJ2dzot7a3n7ZBO2CGHlrDT+YVMxNxVPY+60xEntDa2UCl1eJQZjTAtws5vH64HNI37eBmxzc95cb15feeZz1xZiixDWFGbwgaIsUuKjrA5JKRXAvG0xqCAwLzuJ/3PLIqvDUEoFCS2JoZRSahRNDEoppUbRxKCUUmoUTQxKKaVG0cSglFJqFE0MSimlRtHEoJRSahRNDEoppUaRK9StC1gi0gRUT/LyTKDZh+FYIdjfg8ZvvWB/D8EeP1jzHmYZY7KudlJQJgZviEipMWal1XF4I9jfg8ZvvWB/D8EePwT2e9CuJKWUUqNoYlBKKTVKOCaGp6wOwAeC/T1o/NYL9vcQ7PFDAL+HsBtjUEopdWXh2GJQSil1BSGbGERko4icFJEKEXnYzXERkUddx4+IyHIr4hyPB/HfICIdInLI9fXPVsQ5HhF5VkQaReTYOMcD+vcPHr2HQL8H+SLyloiUichxEfmKm3MC9j54GH+g34NYEdkrIodd7+Ff3ZwTePfAGBNyXzi3ED0NzAaigcPA/DHnbAZeBQRYC+yxOu4Jxn8D8IrVsV7hPXwAWA4cG+d4wP7+J/AeAv0e5ADLXd8nAaeC7P8DT+IP9HsgQKLr+yhgD7A20O9BqLYYVgMVxphKY8wA8BKwZcw5W4DnjdNuIFVEcvwd6Dg8iT+gGWPeAVqvcEog//4Bj95DQDPGNBhjDri+7wLKgNwxpwXsffAw/oDm+r1ecP0Y5foaO7AbcPcgVBNDLlA74uc6Lv8H5ck5VvE0tnWuJuqrIrLAP6H5TCD//iciKO6BiBQAy3D+xTpSUNyHK8QPAX4PRMQmIoeARuBNY0zA34NQ3fNZ3Dw2Nkt7co5VPIntAM7l7RdEZDPwMlA05ZH5TiD//j0VFPdARBKB3wJfNcZ0jj3s5pKAug9XiT/g74Exxg4sFZFU4PcistAYM3LcKuDuQai2GOqA/BE/5wH1kzjHKleNzRjTOdxENcZsA6JEJNN/IXotkH//HgmGeyAiUTg/VF8wxvzOzSkBfR+uFn8w3INhxph24G1g45hDAXcPQjUx7AOKRKRQRKKB24CtY87ZCtztmhGwFugwxjT4O9BxXDV+EckWEXF9vxrnvWzxe6STF8i/f48E+j1wxfYMUGaMeWSc0wL2PngSfxDcgyxXSwERiQM2AOVjTgu4exCSXUnGmCEReQh4HecMn2eNMcdF5EHX8SeBbThnA1QAPcC9VsU7lofxfxL4GxEZAnqB24xrikMgEJEXcc4YyRSROuBfcA68Bfzvf5gH7yGg7wGwHrgLOOrq4wb4R2AmBMV98CT+QL8HOcBzImLDmbR+bYx5JdA/i3Tls1JKqVFCtStJKaXUJGliUEopNYomBqWUUqNoYlBKKTWKJgallFKjaGJQSik1iiYGpZRSo2hiUEopNcr/DxM55HCXw16QAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "dphis = np.linspace(0,np.pi)\n", "vals = []\n", "for v in dphis:\n", " var[0][0] = 100000\n", " var[1][0] = 90000\n", " var[2][0] = 30000\n", " var[3][0] = 50000\n", " var[4][0] = 10000\n", " var[5][0] = 0.1\n", " var[6][0] = v\n", " vals.append(reader.EvaluateMVA(\"BDT method\"))\n", "\n", "plt.plot(dphis,vals)" ] } ], "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.7.1" } }, "nbformat": 4, "nbformat_minor": 2 }