{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "5ec6101d-be3b-44ba-a904-9e722899bcd4",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "\n",
       "<div id=\"root_plot_1776806457628\" style=\"width: 700px; height: 500px; position: relative\">\n",
       "</div>\n",
       "<script>\n",
       "   function process_root_plot_1776806457628() {\n",
       "      function drawPlot(Core) {\n",
       "         Core.settings.HandleKeys = false;\n",
       "         Core.unzipJSON(21933,'WkwIfB8ArVUAeAHdnGlzG0eSQP+KgrEfOQx0o9HoBmI/6LBG3pUthSWPqXE4NiASkrAmCS0IyfJM+L/ve5nVIEgd4/GMZ+2VhFZnXXlUVlZmHf3Xg//a/vh6ebE4Xx7MDp7eXVy8XVx+s3z+5GLx+vLVentwePDi64vV/7xZfn7vYDYCurPaXubbo+f/vTzZmn5gsUevt6v1RQH+c3VxejAbk7xrafbXD+H6GILxuKnHTUMDD1cXy7vrs/XmYFYV8Mn2xzPoHcBvVqfbVwneX52dlcISKzgUHo2i/PLF9ovF5uUKSkdHpny1evnqRtKd9Xa7Pr9e7On69fWE4xcriahp4tnV6+18teHjy+1iK5a+t8w16HZCVr6/QfgSuk+3aTcY35W7zlAkW3RIHpo0bV8yNnlnvTldbp6s/lKkF5Uz8Yv1KYkK7bg6mP1hfDSdjPwzmXTTftzKQaRPj5pJ11d110z6cWUHHdf09AeKmzyeHlXjUdWP666txiGV7fr288vHq3fLs+8p0HS2sF1/OOFg1tdH054/VJ9W0zqkfa2Fpm6yxLTq2rYVxbOr5j6RC5PVUd2ORs10XFd91QSTX19rXHGQMhC3Dx7M2l6p/OlahWYqhaQNVa4ngDMTBhFs18dI4UPSjgb+drb6NapG0+mka+u67ifjzj65av+ZUn6/G0rz5v6B3Jt9OmSTO0K847Ye9d100k77WmU+Plv/8OW9u6kvz/aB469fR0b0w977N7vUB7s3iLzWEPC1toCvqgFc1fz63bkD+A8amK9/jPfom3fni3ewm6m+1n1/1E7R1lHTjOyap6+W2wVFLP741aq83b58jSn7aoEFS56+fHP+fInFsdjT1cn3NDW8/pivD9cvSyJvV2l/ydzHi9PHi9WFw58W7m7Wl5evFqvS4A58vC6mdH9cSuXNIcngXL1YLTGpLxZnl0tK/HGzOgX/PggRA4iw7q4Z61cpn52utovnGqzt5o0N3GcEnl7je6j8eLM6X21Xb5cQd8NmP1xdOiUM00UBFxsY+/a7w4P1axj+9rufaP6zd8sT6l+8OTsD+jLnl/MfT2KC0eg/XW2l5iBTbm0DtOib88eLs+V2Oxh9Bfjl8t32/dR7nz95/PD2Mxr5t+GVsvfWb+DzzpsXL4YO/IouxxgijcL78SUW8OvLIf/ZdTByv1ouzlAfkUd2wlVzVHeTqcr9DTPc+gdmhePs4x0MPdZKuBjgq4QHS6ebohY/DEbkLtNXGoa7iy25N8R+e7vNmVnBHd9Zbn9YLlH/MIbXoBDq/c36HLqYEI8scXy62GrZA3g2AI6R2wlU9Nf3X6zfLh+9XjDRD1L6/qulcrme+OIB5D+UhTJbhbYuticwkML9/smr9Q+fvV1ebJ8w/b1BB4pefX/7zXatWuxKfrG8eHNnge5YU6W5faLa7Wq8QOinjy7OUOzScd+stq/Wb7b7Gjpo7YPFZdG5IWW/1Lc3JPpPc3E0th91ccASVmB18fJjfo6acfdscXlZhojl0rHaT3itqh6MZvVkclh+8wo14++8jlTfxqSYO5o3kcf7fLKrM5q3paTp01LW927W9Yd1VR12zbz3nSmkrtp5JUJe8zevqlnV14dVx29az6t6Nm0P49+8Gs8qzGv5zatmVtXTw/KbV5NZNRkflt+8amfVtKON+M2rKcMMLPmbV92Maeyw/OZVP6ub+rD85jUkUbDqrQvv1XWwhkIQ9dVh1XbzGqq6FoIFx/MaqsTZTiCFulDVUbiVssm8lqpCZEPuFJopCG4cmHndzaoxIEx27bzuKUtPjGi87+bjEZmUEzV4xshpDPetuOgViGK6qSaANDyGKPF0gNZVVNQVz4S6EFVR1/rNZD6GKNmN+tSFqFGPXCWqno8hqiYX2TI9z8dQNUgdhhqomoJXdvtq3iAqujhoHpOL2tTIRfkIKiookkjYbaBKEdtrdTNvoEqZSkbVzxsUCdb75rBv5w00iUYGmum8gSYxVkHmvIGmBl5koqrnE7pPPVPF5hPk1LTocOgpUoKtAoxnYC/vzWzaDO+T2aQa3iFiOrxPmfNL8U7FLu8978QQtt+OeKdH4t2hAz/xXvNeVTGs5q0jqBqjfdZwEKkdAUwEumkCLUBdJZI2RtK4YAnsMUireSt6KqPL0/lU/AGMmvlUAgSqbjSfSkEAk3Y+Hcaw/TuVgsgZdfOpFAh0DN4ggPe2mk+HkTwGyYC+Gs+nSD4rjObdCEc7KvNezSbBBmaiq2eMDRvl3SGc1gVArQwxmJPDtxRzmLRDHYfuDokDd4fFYRvSpoG+GK4xfTrv0UMQ1oQOAGm8ql7z1Y8jp2rRiXnfJIBaAaQVU6kA0oz1Et1DAI0xIfMOft7RNt5BLxK0AivGQEDpGsSNUahGqB0aSsMNRmyE3mkUoAA7ATye0ctVaI3ZkAEf1aQ9pEMQEDx2h4zOSQOELCajQ9QTnQHUZDCc/NEx1cjhWR+O0aBaUJEctkoVKGzrxDFEs1jWmuiU+KgWQihIVSPR0qy21Xd/9QSYvrE3FbztYl2D4tTnSutKW8ScME7udKaBk58K5rWuCIaGSQNESRiMIKYhbGslry38tNCBcTVkU26dIJKy2mjK+AWkt8gFD1g0rR3liH4c5tCJyEGrKastDFEMdUhImrWue1MAzM8Yj/EPE5+mIKwc8qjpkjDUGusOM26foD3lBwwy7SX4MeN2CiU1AZht+oRlBBVkWgHBPaM7yKqmwLBPrajMXOLcIoqoSZ9QKyuTB00wU2oDQ5P4rAwEReJL+uiQwGdl8qBHMWRtYFgX4YAUJUnGEitKYucGzfPa6U5zlEQDl6ml4LVLklf+I9cuCWYDM50ycMv/ZDuOE7FE0yeF3aw8HfjNytBVGLbyeBjDjios2HgYxiqT85YjOWfwAudgDmUbMVU5nof8gNOoxAwTcI5qyzM8xg7sUjxAqCkeRIDaNlSNH2DDAI/RTl2mLnRpqCqEiEpNIY1cVhRqnE6cUXhXGgi+8NfEYA7bH/w0DuaYZ2O+pXySMPDXoDy2bBv8z2xLT4lJ3AHX4irsNaiPEIWlH/UBsmhAMclZOaDB1EMjqqOWaTgbBjLqRhXfnWiswDsaMyosoR07lvQHdIdkQRYDZtAKyzZd2qgfsrDL15IkS5mP5gJbhx/lu4GlbK4feIrSOEWFqSiMUyRofcWBUwQrAyZ8omQsa8bEO2DBI0rmsp7ioCFEM1Elk/bgZ4JKFtoKfNX1tstyVdKfPALDf9IT/E+YZQZYkidMNIGptDYdoGgreE8GKRmsJ394GnZBcMd7eunib1FF0227RRHjHVm0qKHvsNSqhamBgbTVfAHLljJuNWCwHHA0mfOcZEd+mLDQwGCpVQsVpGyK1qmF9qA0izNIaU7Q1pxZKF24anXbqTiwhSGzbmGseO3RbdTM0VGYc1ahYnTbnJUjBZf8qZIDf2LUYhX+gj5NVuEv4TLxS7+iCpuV/CXslFv4i/yc3sCnNHHTC3tCeunJXUA5PktJ1DF5I0hCGa1Fl/CeEsn3FEe+h2r4Oh0MYZINnKOukAWcoy4lCZgthnyA4K/wC1HT4tpEEuxMi3tjW2JK5pIlwOQuM7vBBko5/l46WckvYKHJTgHKcWCnAOz6XKCQE70OmD0+sNIN2liaSXIKY53uDciz1eLdFCCVUAHDZOdUSgsWBmVawhRHd2UIzSnmKBoFxBohh1AXZNM5j9qMKgmEAhYocOxsEXmpgKVgiRijGUZE5xRqs2ChAeBiJEJTABEP7VrEliCpFBZL0cYdSREzJokWzrFqbhTW10hEZKKMoixEGTGaR1FSyb0SEenAMV1YJGhEK2OkFSjmi2AgyharvWtriA+y6pXNDkQRosi7VYkYY3QnEQSM9F8wF5lhsQcxES/G8A7WiBaTHvkkVrTaruAVOYljR040uiMnqCNaDIMSrRIrJjW2aqxYpAdFwCmhkgacEsquApQiyAiUk92EGsKdJEVRNZqSIjs5Gw6KrkApSn7mXcSQKROAnNCy2wAVT1RTrhFIil+VykByIK1NapIZMot8UhmAg5pdqymfXWZM9EPVDCmHPGLKogPJCGGlLO5IIrL0fRAKsaVg9Nuclf0oPDBKfFk6sLSVZja6g8Kp2IMgDDRBFB037wg1rVqKEmyaNxQl3kyk0RFEnGSE1KgXPlAhx4gzac2CyGhgxR7rcugPOHLo73DE0C/U4CIEiujaYqlJoLV5Vwx1wV+cg5imumKkC5BG0Wr2ZrHR0pfN5BQ9NJNGcVc4jWI01BcTXYArCy2P/Z6JttleI52SE2tfzHRi6YudLpz0g52GpGxrR1K0lN2Vou0102GRsmgx1CXTcEN9yh/ReXafhARRGkv5zh/50YOQwRDtMZWUcxmIsKfHVJJO2GmA1GMpCZSNJFkB67GTTW+gzp4bUDcj6jTwbFqgfkagSGFCUVYERjMqBZ9gwEzSYBADBoxkg2rpHYEBEwnqCL3BgIVskLPBLMLDhSSAmhD4tROgdsZCBBuFrgb0LKkRiRkiEyD2rKixymcAPBVD77IJyxGseLAgwSLK+LCF6mkPVM1Y8jMeJjbuWU0jFibqIwoGwoPoD42FOyjDOrJk4CJBD2WYR5atWDxktRCIgcYaKm2yktazlMYKwZQ2ezHQb+PDTk2GMkwjawksvUYPxzraYYR+YNA4sqLG5FGBQdPYsUoJBAYMI2EcS3s1wXGPXaxYuKNNFgN6DCMBoQthrCP02kUgMATkkI/lCcLofqJ/d8iSCxEpkCOetb+EsIpCBMfmYRaFjPOFoAV8RLAsG/ZYRSCWJQoYxKhMmQs15jIPRmHIEUQdpRW7GKBMCkIQUmMVMkEoEoQIc7GMAYoX+WAaBV0OCxCiBOE7wJAQyxYFhCpzJYM+0TKiRrHeIAhVgpIhCFWC4hWEKsEibS2joGQgfEyjoCsFAUKVoGSYC1X0PJ2cIFQJSgaq6IKcoHgFoUpQ9gWhSlAyBKFKUDJQgQ6qBCVDEKoAWcsNEBMpGagvixwsnCVV6HaCKSuXJCp0RCuJcBgjCWYPsnGeYPYg4ykUEUuJnCoGW4KhUbG6qprizZrLMHXtmh1+VogAIQMQHYEsYehg8RMYuhikjPgCQ5gwlGQ+lDHgXGQaTS0PacLQkjC0CUNMwhDHcMUOsWJkeagDxvAUGPK0itLDotso4n1g6RHGcCJEbFEBWcACDNtpLgurgBLnOiFOrqC0BNgEKCljc1kvIldKAmwFWwlhJW2En8swwFwVsAtQMiKX2Y5cyXCJEQNuXclw4RP7DV7NGQYFsJbI2A9QpOHqAtN8r8j0dRFRwLKot4uItXi9ZGPE7ZKAJTTcXWBQ9NLiHonYwcF+A3Cqf9hI0RuAkedmQSd+zLkLUAGL3yAM/GFF7RIDMe0yaVhVYAwp+AOWHsw64yE3U8SPYXf7IzZXxI9pZwhgyamvNDDujAE3brTgwJh38QtLjwZe/MgDew+MiRe/sPg18uIXh/gx824JxTqm8sDQu6gXsPRg6tF9cPETP8bedTwseqzrcu7ESSrhwM9EKK/OTMoDk4+6Jyw9GH3nU7ffmMeAmQ7B7+ohu07ATIjgC1j8xScOWHlg/FHQyJ/Iv34x9LLZx7aPMJO09AkHPjRPeoQDH6onPmH5Zxqw7133nYiPiSBnQGDxMRVEX5DfiI/JQFlanrkdeDxTNu5AMfMDMyzEJyw+PGVxBSw+XGUXugMOfNNZrNqT1gS+buYOnGViOZ79F+aohMXHrOCeVuQrX6cFVspj34stMBKwKcw3mSDG8JlFSbVY/tdr1ixEgjidHhw0LKK7b0aCdkUpkVAHVg2LaE0ItAyT2HkzIdAyTlRsmWGzjX0DBgqmPRPsWd1oVTdKiNbJInYcrSJad29URiZzsJkApWqjVZjlSHA0lzYyweEM6VFC+TttqGCRIGFOHBjUTAjCNNLwEiUkzMlD7ywS7JNhDS4SpBQvO/ZRZR+zTwKU6piZEJZeZxtrlwmS7jpcKJ4lJN2JJDTPBCllXSTMULQhpU4maWcoIaVOJ6GMYDIBA05C6VvNLD1PguphCUjnQYKEmQDpGAhmxqIfsZXinOJoooT2lBEAbF8DQzi0AWcLGlSIx+DJ2dSNUmCojA10YKiGSmBpApZG5hQsccKS6Oqd+zjmS6Grd6GwA5wrVZEvwbENIUu6eMJGLtLTxo6NW0PMxNGeBhi9iN1Y67OpDYyBDGXGfZQ+9961+PDEdi3oCQ9Ct4Glz7UULT7d4NaZO0S4XwmL39U9LTwy0OC6S4RvmbD4mWA4sgaM66y8XE+JgQAsfiYYNi0yP/BDnxYeGU4Dvz4x+ISVj4sqMUzGhxpcfBomEPHhtovfCSZGDU69+GOCUeGBxe8EE4NIp18Y+vBtKvzeVvxOMFp0YfE7wcSQag/dL8PTYYISH7DyYILx3IJ9qsHF1wEWH7D43ZLXgguL3z35GG8dgYww8ovhBqw8jCRwWNQJAjtgJ0DxAYvfaCIclAgxgaEvPBRg5WFEETBbn0EP9AVM2CZ+o4r0aA41yIzsAreEbsLQF/mEa8qDCSbxEcxJDxNMwG62Sg8TTMK0rzyYYBJuDjXQOEQFZnNVeRBlZP6UHVZh6dMRZYdSephg5Fed0l7jFBWYTVbpYYLJfNqXHiaYhCfscApLn/VpX3kwwSSMJyY9TDC1DqCusfQQcyRMHekh6EgYdyvokT7L4+dJDxNMwjih4ndJJvJpX3kwvySsJy5s//IeTq9w9m/qjDD0Rb5jUBj6AkYnpYfoI2HGhPS4MBP56EjjQSIOwXmW7N/igN7B/bP1YjuuyTiL02O4pwdvOa9H5I59IWhjziE254c9YxwRffNjW5YxQ3zNj3Lw6mlg4mZ+lEP/iYz5UQ5dJ/blRzn0muiWH+XQYeJXbB7lkD8RKj/K4Y0Qg/KjHJ4IUSY/yjEbEEfyoxyubY8D22P/ezxT4kB+lMPaE+bxoxweYY+n0WPfidf4UQ4vjXCMH+XwuHoMOSEWP8rhVRBB8aOctjXsN+WwU8Q//CiHzeg94MJY7xnPPWOUkIWfcxARLA9G3ihWQghZeaOTOC3tg+4hnvCBsSBy4KHZHRlgEBr4sIamFeffhzW0pCPNJ8eleWg4ceB9WENPGA/dhzV0D3DAfVhDq4d/7cMaGjncZx/W0KbhHfuwhiYM55eH7idurg9r6G3ixfqwhqYIJ9WHNehBZjxraHPwNn1YQ6cSZ9KHNfQpR87wOIg8nNpxBX1Yw5UWPD0f1jDYGBlRjHQlR87h+Gk+rKGrjxPmwxp66jhYPqzhZD3Sg8Rh8mENHUicIR/WcHYe6TDi3PiwRjoe1tBZwkHxYQ2nYpwPH9agx3lYwz6nH31QA2n7oAYy8UENKPdBDdr3YQ37nD8+rGGfY4N9WMM+Z8j6sIZ97tIAD2vY5wb/PKxhn4dLGAe8jOwBrWGfG7vzsIZ97uEQHtawz8MLNAjnzRr2uWE2D2vY5+H1GUrzZg37PPw8w2OotYZ9bvTLwxr2efh1cdrF0BZQHPZ5uHOGrLxZwz6PKDS8uAgww32L0DH8tggCw2HLIyT2ebhrEZ2FnxZxFybuu59++unw1zrNiUX8+GnOJ3HrhBIfuxBTTmx6k2Jzof0dzjhnyq1dAzdudty4vXHz4sbV1RCU5EN3WriX8v1ys3dHJhMSXdySKQm7ex9POU99++KlB7A9gCqYl0K4TKAMIv9s9ZJDxihwwnvtm31/7el2+pcTxot3q/fPinNo+bbpFPjydPV2dcntIEpNGBac9CVnr8GHi+fL4f6L+AJODGhpgR+9eHG59Oj0kXa2JO7IRuEV+Mn3D5cXL70ORDEPOUcfDFWDF0+a36xG0h45eYw6GGxqDyY/+91wGN35Czj88++Gw1/ah5w037sOBTQoRRzB9oKZw/Ts8ZMcE/c2ix/yGkbCXG+7uvKRAOfSrZtAufgBcC/P22eOh+AdRDGMyCwmQSaAuPZVhmAAKjcZ752g3yVyicN7cAP8dL0+iwP0JuSVlbto7frNhoHGZMI42xZyblhM0h3EYaQ+YQuYd/4uY+BYkX3vRHAtqUDeiEjos4vTzzabdbmB5sgOMIqL6v6bi5MyDs0U3LNigqULzfUKTCksq4KlsONecK+/Hy5fLi9O9+/ZSF2m3rh8d5U44B6uDtrEvlWKNveMJy0+8C7E8vKGDS+pT14vTrwLEDZ5d59vj4fdZb49Y7Qrd52aXdEhGeS7tD3GI+0m33uJKn1q1AMu99wQhUl7kxVHmGSxlBsQp3SGokNqKXhTENT9YnWxOn9z/uflhstWw2UNq1+7YhkmPm+9PN4sucvzx4dXpTN9T3CZsM+mlO6nXvGZqfeWLx4wG+HWDeVI+eZgNr2egmrGhLWrhDJHAje89vAD3US+S7rCHEnvXy0l+RrraiZpD9+7qkriB26rkqqyQ+kgzJICqXspzo0OtGxcyHEZENdLnBDJDaPxeFESKECggUE62XDh6P7iJKyP5GnW9tgXdBzsqa1J+yIZau2NH4sIlsm4NDsoUKpVltmca0UP2qMmLtKYKCtHuOdJi6wcGZwmSD/mfSdL0s3c9A0+tpvV63vLk9U5t+h2l4rCJN9ObweXu7gMe9xFgRvsRdo+f0oyEvcY3ME7dyOdFB2QvMK8qxXs7HKTHWKR0qb8DO+ywztWav35xcVy85XsWdKhFs3C2rdeyb51y0dV88+X5lY3wC1A1QrVIwsARZErKJL20+vR9eK7vG5EI+bx0vCPh1Dg/pUe38HtA+6GLTfO1nH7PES9g+6vUKyiNJOiNHETy17c5YQuhdROuHFoZXv/P9arCxMHR+Du4vU++HR1vhw8h2nHTWiOTdkfn58vXi6f0NDOwGMyTs+W37xaXeKdf7W4eImddRRk+p31u5KWvZepNrB/zfNPqzUNDqnlFlwWvbvanJzFJfE9a1+yvG8q0Xs5x7rcn717fTwQL9FD4rP9xGcfKjkkXitpwS8W7+6tXsZHAVTCRxtu6N3luueGW7ZhfX7FcC1syvClgpuXJmOOQxM/GK5pNxTQnsESHGxPyd0T342JUXAonIYqb81T0WvyIds69ekZ/793DXnfcshFwmWm+KdL7DPFkJ9k+KTIHlT3PiWwEt8u3y3OX+O30maMOkZgSbqFrm/XL/nSgJl78k3H4R+U8N8f5n55sjzT0Ltk89EwtcSofyOofxfx0R7Lcvi7jmul//nKkNy1LcRzdaX+eHeN/jhLeLmbG+Mbbn/n5LMY3vwcCApuA1rGe6vL12eLH68cDxN3Bje0wlAfpOVmOCpPsF/gT4TaP6+XfvxN91IK6e9bfdjrpV0fZW/5eQNH87++iz6yVvDzuugvv+kuUqBh1HYeW1lX+sQC0W+xi24sb+ytfbh5Qdh9gSvsBx4Y+5Hw9PLNOZ8juAYyrq/BKFw3OmpbTsWPx2M+DMMSNOKy5jvKss8yPmLVm/wx+7KaXJwDQz2+1FFWQErsd5XwJW7aEFhIWFm/yO9JvHhC07Qc5ocVkqvvC7kUIAgHN1Y2PvWJihtFHy/eMuMv+J4RiD9k/2VvbxZTOf7BOUwvoT4KqbL1w86tBwjSZ2Cr86hrpm3Tttyui2XP/LIOF4r4w65yy66ZHqoexbjiwzgdW+tImk2VKcnHVX6ShfSs0bKv4eV92x+yhq/6sMnWsA1Drbpk0X2Jhx2Xhj0T8QxZ4/wWELvfXDwz66b/8vnFymUucp68WvBdjD2/6u56Y5SyOF35MQj7eNeTzzciIKV4FpdDZ+imf3h9msKRWRanVbKbkVm4+C5OXy0dlyDh4Ron3JmLLReVc/gMVC4n4//9Q+r0kIWvdx9TJfne+6LDDcfJWh/nmYVsdseuMW5ze3wPYMaeQsaT/GcgyX8/W4mV50M//5JjMq2gOvdog5N/cWUV29QWruVMunH73o7ML5FFsUm3/v1WjI+IID8hFYgKkaUmDBLYD9t3ivAvEEjzTxHBF8vFxS0E4BYHw3r6ab34/yiBJ9vTW/eWb2+lFNiSr35fMsBRXsdXkMK/vf6IDyPtL/Tz7vyDPcr9gQ8tW1hg50APC2CPFxs+sBNnO1gbpf7kp5+uYQ5U+Q2kMmZRluGbSDGd3mE00zTr71CRVkL40dvl5gUf9MISagtubzZ689fPkLQeLokzJBFY+aWhPEnCuZhDLmgcclaCbXOOgQLxxp4y28W8cNaFf5z24cAQZ3Y4RsTJIO5wcOHDI1tekmDvns19T2W4w8yOsicdKBwb2m5Wu8HsQQBPDvjzPIB7+OzZe1jNz4vUnvpmP7xmtvTAloeuPGhWU9cDXR4qqzlqkOke4DGNOh5icNvf4wVxcE+y4cMNcTex3eF2T92NdI9SxCkMD6pwDsZbFRxGgHJOLbAfzrkH/rFtzl47O/Ls0cs7DYg3ZEJDHDQB2yHn6Q45YxXSa3+FDe7nZ2We/RlrADpF2lWH3a/sE4VZ9qOCfKiq4RRQx8kc5EpP6LGkj8PiGmd/OILRohB6KyjbB0rrEY0nR/VkaksNs5K+yM4jwl3lCwZcdp1wMCj2pXYeESebOBTE+RXuO6dDu/OI+JwgLhTH/fiKACs4UrXziDjM5B/8tQ6Kydr3iGTsF3pEQ08VX2E7fBbt494BuPbmQdT8mmegGMzeW5WOzfX3PCKF/hvyiK4vJf1c7v+vvYAJ3jEfGuReN5ozqm/a5DDS7qasXvsR1Vi49oNquDB+UI3V/BNM8bBnE1/oS+udC6n/C+MWxh0=').then(json => {\n",
       "            const obj = Core.parse(json);\n",
       "            Core.draw('root_plot_1776806457628', obj, '');\n",
       "         });\n",
       "      }\n",
       "      const servers = ['/static/', 'https://root.cern/js/7.10.3/', 'https://jsroot.gsi.de/7.10.3/'],\n",
       "            path = 'build/jsroot';\n",
       "      if (typeof JSROOT !== 'undefined')\n",
       "         drawPlot(JSROOT);\n",
       "      else if (typeof requirejs !== 'undefined') {\n",
       "         servers.forEach((s,i) => { servers[i] = s + path; });\n",
       "         requirejs.config({ paths: { 'jsroot' : servers } })(['jsroot'],  drawPlot);\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(() => drawPlot(JSROOT));\n",
       "         }\n",
       "         loadJsroot();\n",
       "      }\n",
       "   }\n",
       "   process_root_plot_1776806457628();\n",
       "</script>\n"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "from ROOT import TH1D, TCanvas\n",
    "my_canvas = TCanvas(\"mycanvas\",\"canvas title\")\n",
    "example = TH1D(name=\"example\",title=\"example histogram\",nbinsx=100,xlow=-3,xup=3)\n",
    "example.FillRandom(\"gaus\",10000)\n",
    "example.Draw(\"E\")\n",
    "my_canvas.Draw()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "f9be5ebf-c0e3-4f80-8497-8f99b08d0256",
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "UsageError: %%cpp is a cell magic, but the cell body is empty.\n"
     ]
    }
   ],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "01eba49e-0f27-414c-8c92-07f22f221b2c",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.14.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}
