{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Access TTree in Python using PyROOT and fill a histogram\n", "
" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "First import the ROOT Python module." ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": false }, "outputs": [ { "data": { "application/javascript": [ "\n", "require(['notebook'],\n", " function() {\n", " IPython.CodeCell.config_defaults.highlight_modes['magic_text/x-c++src'] = {'reg':[/^%%cpp/]};\n", " console.log(\"JupyROOT - %%cpp magic configured\");\n", " }\n", ");\n" ] }, "metadata": {}, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "Welcome to JupyROOT 6.07/07\n" ] } ], "source": [ "import ROOT" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Optional: activate the JavaScript visualisation to produce interactive plots." ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": false }, "outputs": [], "source": [ "%jsroot on" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Open a file which is located on the web. No type is to be specified for \"f\"." ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": false }, "outputs": [], "source": [ "f = ROOT.TFile.Open(\"https://root.cern.ch/files/summer_student_tutorial_tracks.root\");" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Loop over the TTree called \"events\" in the file. It is accessed with the dot operator.\n", "Same holds for the access to the branches: no need to set them up - they are just accessed by name, again with the dot operator." ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAArgAAAHYCAIAAAApvgy/AAAABmJLR0QAAAAAAAD5Q7t/AAAgAElE\nQVR4nO3da7KjupKAUanjzgsxmTrDQBrGPZNBjIz+kb2zswTy9gMwiG9FRYWNMQbZ20qnXn6eZwcA\nALDmf759AgAA4LwIFAAAQBWBAgAAqCJQAAAAVQQKAACgikABAABU/efbJwDgNTnnX/cJIex9Aru+\nBIDz8MyjAFxLCGGapsf7DMMQY9zpBLz3zjm+OoCbIKMAXEyM0SYVcs7TNHVdZ3/i83MfwFYIFICL\nCSEsY4IQwn4pBAB3RqAA3ELOuUgzaFriQfrhk+4IzxwfwPkx6gFolvc+hJBz9t73fa8phxCCbBGy\nW/Fcu8/qDsULee81MtBXfPLpAM6MQAFo2TRNfd8757QTQ4xR+kIOwzCO4zAMspttuZD+kl3X2R1q\nlb30bRzHUXeQV3xwfABXMgO4sq7rnHPDMCwfkr/xruuWG+2WcRztbsXd5bOWt8dx1D0lMrDns3pA\nAFdBHwWgccW8C1JtW0WqQH76FwmAYRhyzkVHh2UuYfVFQwgzYymBy6LpAbgXHTQhFX+MUep7JQ0T\ny+hhNUpwi0BEIoxpmrz3xUhOAFdEoADcju2omFKSxos3SHIipbTcLsdMKdGZEbg6AgXgXrz32lFR\nGiDf+9EvLQ7SI6GIA2SoxTzPGjFIgmGDswdwOAIF4EYkJui6bjmtglpNMMiIx6LngTMNDXZspJ1B\nQSKGbc4ewDcQKAA3sjqBUtFvUad6fLCPJQ0QMiRS9uz7vshSvN26AeDrCBSAG5HqP6UkFbl0Ziw6\nGWiSQLoiSu5huZyEPabEAXa4hMQK2l/y11WsAJwWq0cC1ya1+OpykavLPC4jg3EcNR+gOy+7FBQP\nFYe1G1fXt+SrBrgoAgXgdrQbgV1fajl9wupuz7MzRjPqAbguAgUAAFBFHwUAAFBFoAAAAKoIFAAA\nQBWBAgAAqCJQAAAAVSwzDQDYGEt7XMKTwx4JFAAA22Ps/ck9H8zR9AAAAKoIFAAAQBVNDwCAy5M1\nzGoPbTKJuCyb/skphRAeLMR6WmQUAABNmaZpjwVLp2l6KVDIORenkVKyPQNejTy+hYwCAODy7C91\n733XdSepg+1p5Jz7vtf4QCKP8y+ZRkYBANA+XR9VQ4qcs7RKLNsI7EOrR5Md7F1ZJfVxy4LEBJJm\nkD11jdYzI1AAADQu55xSCiGklKRijjH2fS+3U0p6e/nQchihJAZ0/xBC3/dyW17l8cl0Xed+Mg3T\nNF2g18IMAMCmvlu5OOe6rrNbxnEsNi7vDsOwvC1PHMdRto/juHoo2WGe52EY9NolIBgN3bJ8oeM9\n/x7RRwEAcAv2t/u8mA9KMw12zxCC1P26T0pp2QFCuxrYpg2hyQYxDMP5OyUUaHoAANxCUUNLLwTv\n/eM5Cu3oSokSiuEPwzBIC8VqH4Xi1/kFGhoWCBQAALfjvZf+BNIQII0CvxqGQUIEmyeIMc4/jRSr\nfRqujkABAHAvUtmP46jZAp3wYNku4L23LRHyRGeGLcgNGfQoLRpXTBs8QKAAALg1W69LKKDhggYB\ndv8QQtd1KSW5m1KyQy6X+1/ehl0oAQCYzzrqwW6xbQ1d18ldGYMgOys7GkIHLNhXKZot9KVl++Pz\nXJ7qYZ5/j/zMSqAAgE15f4HKRSY7komS3CIZ8FJuwM6p8Oo5vPGsTTz/Hl3gvXxJe71IAOCKGqtc\n2vN8oLDLPAp2HOrxi17w6QSA7+I3W0u2zCjEGLVzR9d1ks/p+/7ImvsS+S4AaBuBwiV8relBxooU\na2tK3HDAiBECBQD4Or6Kz+/592j74ZGSSCgmzZYGiPOvkQUAAKwt+ygUTQ92Ja4NXwUAAByGpgcA\nwMb4Kj4/mh4AAKcTY/RrPjngmadBjDEWFZ+c8PJns87osFyXsvYU3b66lOWWP84/ndvpNwfPOXXA\nFQEAHqt9FcuSzePC46Mt51W0D+nv0rORi7WnJ9Wurlut22UOx2KCyPnnwnW7fYoeSh7SMiwO9aBO\nfL66bK1aJVAAgK97HCi8erQHgcJp6TzQWusXMzpLXa57ak1vi0j3EXq0ohjtbm4xz3QtkHq+SBtc\nFGo1r8WgXgA4Oe99zlm/tHV5RlnTWR6NMUrKXe/apgfZLjQhb495WF3Q971mDsQ0TXaLduYrZnHW\nq5anrDas5JyXB18eaisNBgq1mOjb5wUAcO5nkQVLH+r7Xn4QD8OQUpLmdv0BLU34MrxuHMdihWgJ\nAuTpXddJeCHHtL/dD+jTIMtLPu4i8LjTnjyqheD+XsdS17aW2xqChBB2qeyezDxcRXtXBACXU/sq\nLn4HFzWRW7ToS61vmx6KBL7GEMvmCX26M9n4Z7pEfMieias3PeiFy/7FbnZAgFo2IjyoyrXto3ae\nz1eXu6z1AABAzVz/1fvMz/1iWWdhf3bbjfLjXpIKx4zSl6TI6hlKy4jc7bpOciGSMkkpSabEXp2U\nhhSXNsHYS5jnWbYXwyVCCNM02QmNPtFg0wMAADIy3zmXc5aaO6Xkvd+16UHXRNTVq1NKdulqTWnY\n+j7GOP9kO7SfgbQpaMyh8UTxiiGEcRy1/UU6ZLifGGKTiyKjAABoRFH7aqCgP7jtDMJ7KKIQ+Vmv\np+EWnRb1tj1b99vCy977YRhWsyPSiXLjxMmTTRRX0d4VAcDl1L6KH8+j4BZD+1b7KNgRg3agoDOt\n+M9sP0DtpYshke7vTgl2rGMxPFLuLkdayt3V4n1wbk9eBRmFK/H+3wePzvOfw84EAN6m4xGUDmF4\nwHu/2vZvD9L3vSYMtONk13XaA8BuP5hMoWj7KOgl2z4KzgyIkCuyT5GHdMCnHtyWTFG8nycYWpuO\nu+EJxh9HCYVa0FAchNgCwB42/yp+foaA2p47zTHwqgen8fjMlw/Vtj/p+feotWqVQEEtg4DVIxAr\nANhcw1/FzXj+PaLp4XpWq/ZXwwgAAJ7RYKBQm6Gz7fBWo4dnIoZ5/kNgAQB4RoOBQtsBwfM0FKBx\nAQDwtgYDBRRIHgAA3kagcHZv9GEkMgAAbIUpnE/tvSp/nv/ov9VHPz4vAHhHsQC0kkUQDliI4WyK\nxTOFzgD9YHLG8LflnnaFyQ8RKFzMq9W8DRpqoQMAHGlZq+lSBXfT931RGiEEnXlJVntaPkuWgfj1\nyFtNVk3TwzXsV8HXkhaEFAB2klKyP3a3WrvoWmSBx2KjRADaJV9WhpS1Kord3MNy23ZeKTIKWEdH\nBwB7kBmUbSUXYyymVdYpiovFHu1224Qh60dL+8Xe60NuRVaDLDYuswurz805r661LWKMshjVp6f4\ng0DhpkgYAPgWXbNATNNURAOyBKIsnmQftdu7rtO0xDRNspqzPuX83R1kxcjlRjvC/0HOYJomDZiK\n9ExKaZ7nLaOlJxePuorGrsi5/8q/5l8UQEtqX8XOuWEY7GqQetv9rKBYrA9pd9YlFufF+pDFKot2\nzzN7cKpygfa67LPczyqRumJkccBf18l8vrqkjwIA4FDyY1ea3pftDtJyr7+Si0aK/KNo4Lc/oDfM\nun+F5FRcfeHH2WQdZMhD3/fzPIcQbKJlKzQ9AACOpvVZ0e6gj+pt25bvvdfOfd9aLXpvMca+7yWR\n8GSVrwU4TZOUpw6deDzG8kkNZhTuudYDAFyIVIdSERaBgkQJto6U2/K/fpOfvxfCG6SHwTiOj3sY\nFNW/3l7mZiRo+PCsGgwUCAgA4OSk9kopLRMD8mtYxwTKMMJlWLDVJAGnstrgoi010vLifiIAHSQ5\nTZMUYxFdbdWps8FAAQBwfl3XrbY7SHUojfRCfv7FGFNKmjPWp7c3DYO9duecJBik4tctfd/botg1\nv+Ib+/3tfVNX9JUVIFl2EsCHNvkqlghgda4h/ZG93OE+Prz859+jpqpVR6Bw2RcF0JLGvoqb9Px7\nRNPD6TAlIgDgPBgeeS5ECQCAU9krUHjQsWI5b6VOu1H0SZHtq31VZFrvJofHLH2rCcD7f+XfV14d\nAHAGewUKMrhluV0GiRYDP6SHpy6TpXtqz08dbiu89zIwxvaAbcypFoYmVgCA29o4UJAcwIPKuxj1\n4ZybpmkcR5mS047xkKmpZPs4jjpkVufcyDnrmJltr+LmzhCaAABOYvuMwoOZNWUaajs3ZzEtlx0n\n6kwEIDvI3ZSSPUKxCtlFnS3Jf55kBgDguzYOFKTfwOpPfAkClr0QiqfLxsfDQ+12mbTrg1P+vvPE\nBwCwtxijN55fzmC1RvALF/rpqJ3tauesHfWUPlRsXx5BZnLc5DyPG/Ug81cvtz+5zNfzmYPl5+ZX\nL13IAfgpD6BJMruiLqw8DENKySaP35g+qOu68UfXdUWfttPSdZuK/nlWkWVXy8Uzl/q+32qW64Pm\nUfDeD8NwzPxZ153lg/gAQNskStCKXBuUP6nabXghP6M/POABJALQ2qqWV5BFHJbXoqs81I6/bW17\nREZBL0neTl0HUx61YdGDy16dEhwAcC3F93yMUZLNy9pBusZ771/68rdLTcrLafLYHtnWvrLPOxfz\nLum8b+/Wqr/Vay+eXpAo5Mls/TOOCBSke2PRxKLLgtk9tWuCncfbHsct2iAelxcA4DyGYZBQwNbT\n+oUvvd3lG957L2Pi5Gf3S33R9CCS1ZeWjnEcNQrRtL+IMR5cj+hyVr9OCKRdOoqOCNM0rfbzkDkI\n5nne8qf1vA/nnLzBS/JRsHtqe5VzTh+St01uyzAKuS2xpxzc3tYjbHcRB3Huv/Lv2yey4sznBuC0\nHnwVD8Nga2X9/p9N7bD63W4rDrvdHqE4TlHdyGH1iXr8BxXWTvTytSiWJyBnqz0wbJWtTxnHUepH\ney1SILberJ3Dk2f7/bUeZLlMDe7shEu2p6EWk4yx1PkYuq6jSQIArsI2Ckh/AvkFbPdZDnx77xe/\n5CHsK7qfdnC9UYzSP9Js+ij0fV8UQgjBbpnnWZIHMUa7XTIN8nSpHzfvn7FX08Ncz3sUjTFSFhIZ\nFcVkt9ujyRRM4kIjYS5NZ3pgMCeAtxUtDvJl7h52UHuD7ehgIwxtB3fOyYALeena3D+7si/6YdVu\n+/zJtWvbyurIyVd9P6MgalHFq9txAO//ZYAGgFdJ8/lyPoDlnkUfAufc873zbIZAnmKrYTvgQtYT\nWE7wc4DiWlZPQH5Ur8YQRfWvt4uIR4OGz072gi36j13xik7eD0BPz/779kkBOLXaV7HUO7Y93nZH\ns10KnOmUIPVfrY+CnUdB9tReC0VfB/taenf1sHtbdiywF6u9L+0+ts+B3V92W3bU2LCPwvWq1ccI\nFHZ1oVMF8EUPvoqXv1eLutD2Z1RFt8QHRyuqzOK3u32oiCoOVqxFoNuXEdKyoJaFszz+hoGCny87\nPdEq7693Rdrqf/58/oVOFcAX/fpVrNny2iD51btv2+o426oVwupuy30+vKjnq8vrVauPESjs6kKn\nCuCLrvhVfDfPv0dn6cy4odoEW3xqAQB4VYOBwlUCAsYZAgDOr8FA4RKIEgC07YQL8+I9BApnQas/\ngGZcJbO7h/b6ZxAofBnxAQDgzI5YPRIAAFwUgQIAAKii6QHvYEIFALgJMgr4CMM3AKBtBAp4AfkD\nALgbmh7wGo0VNJdAMwQANIyMAjZDMwQAtKfBQMFXfPu8WkP+AADuoMGmh8amxDqzZTMEAKAxDWYU\nAADAVhrMKJwWP7sBAJdDRuEgRAkAgCsiUPgm+gMCAE6OpoejERwAAC6EjAIAAKgiUAAAAFUECgAA\noIpAAQAAVBEoAACAqgZHPdSWdWBqZwAAXtVgoEBAAADAVmh6AAAAVQQKAACgikABAABUNdhHAV9k\n175ismoAaAAZBeyFBTMBoAF7BQoxxtWNIYQYY87Zbs85P9geQii220NtedI78P5f+fftE9kXyQMA\naNVegUJKydbuOWfvfUpJav2+7/XREELf97KP3S535Xbf9zYmkEPJq9RmTTiD5uMDa57/6L9vnwsA\nYDN+21kHcs4xxmmanHPjOIYQZHsIYZomfS1713uve8r/Eit477uuk9sSNMj+McaUkh7Kez8Mg4YR\n3m98RZ9YBgo3qUT1wm9yvQCgTlUNbWL7jEIIYRiGYuM0TXajBhBSwdu7EmTYR3UHuZtS6rpO99Fg\n4sz4qQ0AuKiNRz1IlwLnnDQNqCK8su0LxdPtRg0glq9ibxevhfMgtQAAV/eFUQ/S7jCOo9y16YEH\nns8c+Ne9fS140q26awBASw4NFGKM3nuJEmrZgs/Nr9vpTG6L/AEANOO4QEHaCIZhmOfZRgm2U8KD\nnME0TfvFFtgc3TIAoA0HzcwovRSXv92lGULvSqCgAYHMo2B3dos2iJzzk40XAADgVQdlFGSoQv6b\n+xnIoKMb7IiGrut0HgU7OEJiDh02OU3T+addAgDgovYa7mlnR5C7y33kpe3ESu7v8RH2WctZGeR2\nkWA41QBW+vw7CgHAzZyqGtrEXk0PRTE9KLUQwjzPRaODPmt1uyYklg8BAIANtRb4nCqU48e0oxAA\n3MypqqFNsHokAACoIlAAAABVBAoAAKDqoHkUjlSbkrmxRiMAAA7QYKBAQAAAwFZoegAAAFUECgAA\noIpAAQAAVDXYRwHnxMxLAHBFZBRwNI0YAADnR0Zhe1SE1jz/oUAA4LoIFDZGpbikbQ0UDgBcDk0P\nu6NJHgBwXWQU9kJ8AABoABkFAABQ1WBGgbUeAADYSoOBAgEBAABboekBAABUNZhRwPnZcZL0+gSA\nMyOjgC9jcgUAODMCBRyH5AEAXA5NDziUjRXIJQDA+ZFRAAAAVQQKAACgiqaHDZBCBwC0iozCp4gS\nAAANI1DYBd37AQBtaLDp4VtrPRAcAADa02CgwFoPAABshaYHAABQRaAAAACqGmx6wOXoyBH6eQDA\n2ZBRwIkw1hQAzoZAAV9D/gAAzm+vQCHGuLoxhLB8KOcs23POy+0hhGL7g0PhWub5j/z79okAANbt\nFSiklIra3XufUpKH7FQHIYS+751zOee+7/VZcldu931vY4LaoQAAwLY2DhQkB7CsvKWan+c55yzz\nHGjFP03TOI4555xz13W6ve/7rutk+ziOEhk8PhQAANjW9hmFEMIwDMXGlFLXdXpXIgD3U8eHEGR7\njHGaJt1NIwDZQe7WDgUAADa38fBI6VLgnNMEgH3I3pYdijpe9tGN9im/HgoAAGzuuHkUarW+TQ88\n8Hzm4I1eC8z6DADAqgYnXKLWBwBgK8fNo1DLB9hOCQ9yBtM01XISAABgJwcFCkXDgQxwcIv2CNlH\nOzqs9mCoHQoAAGzuoEDBDmfIOU/TJEMY7P/u7xENXdfpPAp2cIQcSmIFeygAALC5g/oohBC6rtNu\nhl3XaS5hHMe+73Xkgp1wyXuvTxnH0R5KYwh7KAAAsC1/ZNe/B+MetdHhpe3Lh7w/9IocKx9uZLkc\nFOUJ4IqOr4b21tz1HPUOUbFta3XdSIoUwOW0FyiweuQ7WA15c8QEAHBODc6j8BXUc5/TMiQOA4Dz\nIFD4CPEBAKBtND0AAICqBjMKtbUeGutdAgDAARoMFAgIAADYCk0PAACgikABAABUESgAAIAqAgUA\nAFBFoAAAAKoIFAAAQBWBAgAAqCJQAAAAVQ1OuIRm2NWhWFYDAL6CjAKugSUlAeArGgwUfMW3zwvP\nInkAAOfRYNMDaz00wMYK5BIA4IsazCgAAICtNJhR2A8/bQEAd0NG4VlECQCAGyJQeB997gAAzaPp\n4WXEBwCA+yCjAAAAqggUAABAFYECAACooo8CLkMHntBNBAAOQ0YB18NQVQA4TIOBAms9NIb8AQB8\nUYNND6z10B6NFcglAMDBGswoAACArRAoAACAKgIFAABQRaAAAACqCBQAAEDVcYFCzjmEEEKIMRYP\nxRhXt8tTYow559VDFdsBAMC2DgoUYox938vtlFIIQR/y3qeUZLud7SCEIE/JOfd9rzGB3JXbfd8v\nwwsAALAVf8ysA977YRikUpeafhxHyRaklPQc7G7ee9nHOSf/S6zgve+6Tm7LoewleL/XFTF/8Bnw\nLgA4uf2qoW85rulBswg2nZBS6rpO72oEILGC7hljnKZJd9MsguxAUgEAgJ0cNDNj13WSRXCLIMDG\nDSEEaYYoOh/YjELxFNwTqQUAOMZBgULO2XuvfQskYhC1Wt9mGh7QJIR6Y1mHxtJEt+L9v8QKALCf\ng5oepGPBPM/zPA/DYDsnbm5+3U5ngs0REwDAwY4IFLQPo72rHQtqEYPtlPAgqpimiZaIW5nnP/Lv\n2ycCALfwnQmXtFmhaDjIOctDRd0v+8jcCa7SgwEAAGzuy8Mj7fhGu90thkpqSBFCmKZJnlKMrnRb\nj0tZXdSY37InQX9GACfU3vDIg67HzpLknNMIwP1U/HLbJhiKpxTRgN7WwEIf2uqKVqMER7V0GgQK\nAE6IQOEjmhJY3f7goee37x0oUCedB4ECgBMiUDi7PQIF6qFz4g0CcELtBQoHzaMA7MfmfggaAGBb\nLDONptR6lgAA3kOggKsieQAAB2iw6aE2hXNjjUZwf8cK5BIAYA8NBgoEBAAAbIWmBwAAUEWgAAAA\nqggUAABAFYECAACoIlAAAABVBAoAAKCKQAEAAFQRKAAAgCoCBQAAUEWgAAAAqhqcwpm1Hm5OF31g\n1SgA+FyDGYW54tvnhaOxTBQAfK7BQAH3RP4AAPbQYNMDbktjBXIJALAVMgoAAKCKQAEAAFQRKAAA\ngCoCBQAAUEWgAAAAqggUAABAFYECAACoYh6FEkPwAQBQDQYKn6z1QJTQGNZ9AIAPNdj0sO1aD1Qw\nbSAEBID3NJhR2ATxwaXN8x8iAwDYBIEC2sS6DwCwiQabHgAAwFYIFAAAQBWBAgAAqDo0UIgxhhBi\njE9uzznL9pzzcnsIodgOAAC2dVyg4L1PKTnnUkp2qoPa9hBC3/fOuZxz3/caE8hdud33/TK8AAAA\nW/FvTzDwkhCCc04re+/9OI6SLUgp6Tl474dhkLpf9yme7r3vuk5uS9BgL8H7j66I+Xnaw3sK4Egf\nVkMndND12Fq/2K61vjMBQRFA2ICgOJSNLRyBAhZWh0fy/gLYSXuBwhFNDxIHSP5g2RfBRg8hhGma\nnMk92H1yznqoPc8X7WNyBQB40qF9FKSaX/ZFWN2/67pnDmsTEvpCr3r3mnABJA8A4BPHBQrDMEhK\nQHIy+3VCrK318MBOZ4KTmOc/+u/b5wIAF3NcoGAjA5sGqA1xlDaIx/vIbrREAACwkyMChWLIg1U0\nHOScpcWhqPu1a8LqoQgUAADYyUEZha7rNKOQc56mSe7GGDVzUGx3JgmRUtIuC13X6TwKsgOBAgAA\nOzluFIftM2gHNOpIB7doktCAwDlXTJagt4tRlwyPxGO8xQB21d7wyEOvpza4sZhH4fmnrE7MQKCA\nB5YDI3mvAWyIQOHsCBTwGPMvAdhVe4ECq0fiXogJAOAl//n2CQBH01iB+RkB4FcECs5RYQAAUEHT\nA1ECAABVDWYUams3PN+7hGZsAABEg4HC291NiQ8AACjQ9AAAAKoIFAAAQBWBAgAAqCJQAAAAVQQK\nAACgikABAABUNTg8EniVnXSLUbIAYJFRAP7CTJ0AYBEo4L5IHgDAr2h6wK3ZWIFcAgAsNRgofL7W\nAwAAEA0GCgQEAABshT4KAACgikABAABUESgAAICqBvsoAB/S4Q+MnwQAMgpAFQMmAYBAAfg/5A8A\nYImmB+D/aaxALgEABBkFAABQRaAAAACqCBQAAEBVg30UWOsBAICtNBgoEBAAALAVmh4AAEAVgQIA\nAKhqsOnhGYySBwDgGXcMFIgS8DzWfQBwc19oegghFFtijCGEGGOxPecs23POy+0hhGL7J6gG8Bjx\nJYB7OjpQCCFM02QreO99Ssk5l1KyIxtDCH3fO+dyzn3f61Pkrtzu+34ZXjxvnv/ov7cPgibxkQAA\n4Y8cTKh1/DiOkleIMaaU9By898MwSN3vvdfd5H+JFbz3XdfJbTmgvQTvf78iksl4CR8YAM97phq6\nlkMzCn3fD8Ngt6SUuq7TuxoBSKygjRQxxmmadDfNImi0sdspAwBwa8cFCiEEzRYU2+1tCQiKzgea\nUZDty14OAABgDweNepCUwGrfw1qtbzMND2gSQtWmcF7y/h+50ViaCACArRwRKOScU0rjOB7wWu6J\nWp8mZwAAnnRE04M0N8gYSPm53/e9JhJqQxxtp4QHwyCnaaIlAsfw/l/99+1zAYCDHJFRsBMhhBCk\nA6PU7kXDQc5ZWhy0s4Jud6aRQuZR0EcJFHA87/8lIwXgDr4wisOOe7TjG+W2PlQMldSQQmIIeUox\nutIxPBL7WE0h8PkBsNTe8MgvT+EcQui6TrsfaqbBOTeOY9/3MheTM60POWfvvT7lsK4PuDMbE9Du\nAOBWvhAoFKGWDnp0i6GS8zyvjoesbQcAANtqLUNC0wP2tswo8EECoNprevjColBAY2iMANAwAgXg\nNeQPANzKlzszAleksQK5BADNI6MAAACqGswo1NZ6aKx3CQAAB2gwUCAgAABgKzQ9AACAKgIFAABQ\nRaAAAACqCBQAAEAVgQIAAKhqcNQDcDwWEAHQKjIKwJaYqxFAY26UUeAbHJub5z98rgC07S6BAt/m\n2AnrPgBo232bHmhLBgDgVw1mFGprPTj3X0d8AADAKxoMFFbXeiAtDADAGxoMFIDvslEpGSwAV3ff\nPgrAAUhlAbg6AgVgGyQPADSJpgdgMzZWIJcAoA1kFAAAQBUZBWBfLAMB4NLIKAAHoTECwBURKAC7\nIH8AoA00PQB7YRkIAA0gowAAAKoazCg8XusBAAA8r8FAgbUeAADYCk0PAACgqsGMAnBazKkA4HLI\nKABfQFsYgKsgUAB2R/4AwHXR9AAcgTkVAFzUcRmFnHMIIYQQYyweijGubpenxBhzzquHKrYDF+L9\nv/rv2+cCAFUHBQoxxr7v5XZKyU514L1PKS23hxDkKTnnvu81JpC7crvv+36lClsAAAoeSURBVGV4\nAVwOsQKA0/Krsw5s/zLeD8OglbrejTGmlPQc7G7e+3EcQwjOOflfYgXvfdd1cluCBnsJ3q9fEb3N\ncR6rYQGfTKANtWrouo64nmV1rhW/rfXt9iKAsEewAYRbC0EIFHAhfDKBxrQXKBzR9BBCKEptmiat\n6fWG3J6myf0kD+x22Sjb7VMAAMB+jh71oD0MNAdQq/W7rnvmgDYhIeprPcij/yw3Nhb9AQCwlUPn\nUZD+iV3X7Voxz2uc+6/8qzwKAABWHBQoSHcE59w4jsuxjqtPkTaIx/u4v1sxgOtiqCSAczqo6UES\nCcv6vtiYc5YWB+2soNudaaSQeRT00VqgwHcursj7f+nYCOA8jhv1MI5jsV1mTNLhDLqbVPzFUEkN\nKSSGkKcUgyOc6W5aixL4CsbZLD+rfEqB62pv1MMR1yPVebGxqPiLje7viZXc3/0NbXdFO1TS/RYo\n8P2LM2OoJNAAAoVdFPMoLB96fvsyUOA7F1dBagFoAIHC2REo4LpIgwENaC9QYJlp4CyICQCcEMtM\nAyfCatQAzoaMAgAAqCKjAJyaTS3QNgHgeK31uXiw0ENjV4q2MREIcFHtdWZsMKPAqAc0YJ7/0E0B\nwBk0GCgAbbAxLkEDgG8hUACuhFQZgIMx6gG4JHIMAI5BoABcAPkDAN9C0wNwDczFBOArCBSAq6K/\nAoAD0PQAXB45BgD7IVAALob8AYAj0fQAXM+yvwIzPQPYCRkFoDW0RADYUGtTUrPWA+7mcVhAdgE4\nWHtrPTR3PT/vEB3CcUOrQQN/AsCR2gsUaHoA2kFMAGBzdGYEmsK8TAC2RaAANI5mOACfoOkBuAty\nDADeQEYBaNM8/1lGBky3AOBVrXXO9N4791+7hW9DoJZL4K8D2ByjHgBcTy0g8P5f+Xfw+QC4kNYC\nnyKjwA8mYInpFoD9tJdRaO56mnuHgD0wnyOwk/aqoeauhymcgVfQfQHYFoHC2bX3Dr2BQnAUgnPu\n6UJ4vo/CFaMHPgmOQnDOHVgI7ZV2c9fT3Dv0BgrBUQjOubcK4deg4XKxAp8ERyE45wgUPsCoBwD/\n73JxAIC9tRb4tBfKvYFCcBSCc27TQrjuQAk+CY5CcM6RUfgAMzMCeBPzPAJ30Frg014o9wYKwVEI\nzrmtC+GifR75JDgKwTlHRuEDl7yenHOM0TkXYwwh2Ifee4fe6vN1xFPae6Ezn9thL3Tmc3v+WX/H\nDf8UU6fXmFWwWyiE45/S3gud+dyOfKEzu17TQ8657/uu65xzfd8PwyBBA4Aj2bSB9/88+SwbXhQp\nilPlIQCo6wU+3vuu63LOzrkYY0rJXsKZY8Yzn9thL3Tmczvshc58bp+/0NMtFC8nIT4/t12f1di5\nHfZCZz63I1/ozK53Pd77cRy1xcF7b5MKZ/4onPncDnuhM5/bYS905nPb9YXea6342xvP+udUhfDh\nU9p7oTOf25EvdGYXa3qQRELRLwHAJRStFX/f3XEFy/cO/sazjnlKey905nN7+1ktuVjgIx0UirYG\nbYmQu985MwAAnHPNLS10sYzCrxp7ewB8iJ+DwIdaCBRoiQBQw2AK4EMXW+tBYgJtaLAbAQDA5i4W\nKDjnuq7TMQ5yg0ABAICdXC9QyDlP0+S9996nlMZx1O0hhBBCkW+4idVZp2TmyjtMSKXv/vJiKQR3\np0JQy98PNykE/Rgo++hNCkHULvZWhbCN+ZrGcRzH0d51znVdJzM2DsPwtTP7EuecLZD5p1OnFMh1\n3+hnDMNg3317sfcpBLlGLYSu6/Sh+xSCkou1fxH3KQT75yD0ofsUwly/2FsVwlYaKSb7zSh/J189\nneOM46ifePu1WBRC28FTcXV6926FoB8AiZvl7q0KQWiWUQvkVoVQBAfqzoVw5z+HTTRSoS5/Pdzk\n7R/HcRgG+fQXJWD/TmrfHQ2QWsFu0Yu9cyHo5+E+haDkz9/+RdyqEGrffncrhCLDqtvvUwgbul4f\nhaU7T9coLW2rjW22QEII0zQddlZHCiHMf0+eMU2TXvs9C0Gu+m6FIEIIqwvF3aoQcs7+h+2zdZNC\n0BphtS/CTQphWy0ECqvu2aXRumHkJN+PznTtvFshhBC899M0afrd3akQYozTNP0aN9+B9OLquq7v\ne914q0LQICmlZGfsvVUhbKXZQAF3E0KQ9ceLBMOt5JzneR6Goe/7u8XKOWc7DOq25nkuhoDds2oc\nhiHnLH8RrjIuDE9qNlC459+GdZ96QhMJ4zgWV32fQrDkO1Gv/SaFIFct2Wb5PPR9r98DNymEJe3s\n7G5WCDYysOsB3aoQttJCoMB0jUv2D8M5l3O23xeNkUSC/Iqy2+9TCDHG2nJotyqEYRjkl7QOEZSP\nxH0KYflXoG3w9ymE1RpB3KcQNvbFjpQbsp1XbzU8Urm/e/nabvB2sFx79OoK8/0KQfu624u9TyEU\n3GK8qL3dcCHwSZj/rhFuWwgbaqdCtdHPDd/75VXbSLnhIUASFxb0em9SCPOiHOwAufsUglX8Rdyn\nEPgkCAphQ35uqOfXnXvurNIk253L5FaFUPsTuFUh1NyqEPgkOAphO00FCgAAYFstdGYEAAA7IVAA\nAABVBAoAAKCKQAEA8H0yU5b1TGdDfZZusUtdLBe8qPHeL2dvLGYoWU5TUez55GtdDoECAOAsdCqU\nYRiKBd5qlbSduD3GKDOw6cQJsuDFr/V313UppWJjSqmY2rI2FXSMseGRAYx6AAB833JGxRhjSkkq\nqdp8ixJAaP3tvS+mX3Q/q0Q+ruxyzn3fj+NowxHvvd3i/S81ZrF/M8goAADOyKYTpmmyCYZVxSon\ndrudh0qbCWxzgxzZJgzktn1FzS7Y1o07LDdFoAAAOCOtg2VRhmWqoFBbu8GmHCRLIQ0TwzCklPSY\n0tihzyraHWS9MfeTe5DZHsdxTCk1HysQKAAATmGaJtsJcZqml9YNt9W8W/RqlI0pJVk8zDkXY+y6\nzsYQbtH2YY+mK5Tqs3QFsrb959snAADA/7GRwauN/ct0gtbikjnQlICNAIrwIsaoO6yewDRNRUPG\nSyd5RQQKAIBT0JXB3xNCsCMX7EAJ2V7rC6l3pTHC1dsd7olAAQDQAgkUYozFr/xlj0Wt9YvQQXow\naBODbs85655FV4nlzu2hjwIAoAUhBJkOoRi8YNMMMq2C3l1OsaATKtgUgm2ekCEY+qzl7AvtIVAA\nAJydVM+/5v9zztJ8oH0YdSYG3cE5p4/azoxC7tpeCHJMu4NEG9pHsu10gmPCJQDAJUgdX8QKxYRL\nxc7L/Z95dLnzcrfVI7Q64RJ9FAAAF/BSBfx4588P1V408ABNDwCAC5OGhu+eQ7F8VGNoegAAXNhq\nk8RtT2MPBAoAAKCKpgcAAFBFoAAAAKoIFAAAQBWBAgAAqCJQAAAAVQQKAACgikABAABUESgAAIAq\nAgUAAFBFoAAAAKoIFAAAQNX/Ao465SAjKC7mAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "h = ROOT.TH1F(\"TracksPt\",\"Tracks;Pt [GeV/c];#\",128,0,64)\n", "for event in f.events:\n", " for track in event.tracks:\n", " h.Fill(track.Pt())\n", "c = ROOT.TCanvas()\n", "h.Draw()\n", "c.Draw()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.10" } }, "nbformat": 4, "nbformat_minor": 0 }