{ "metadata": { "name": "StandAloneMuonsTrigger" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "code", "collapsed": false, "input": [ "import ROOT\n", "import rootnotes\n", "from DataFormats.FWLite import Events, Handle\n", "import math" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 1 }, { "cell_type": "code", "collapsed": false, "input": [ "handleGen = Handle('std::vector')\n", "labelGen = (\"genParticles\")\n", "\n", "handleL2Muons = Handle('std::vector')\n", "labelL2Muons = (\"hltL2Muons\")\n", "\n", "handleL1Muons = Handle('vector')\n", "labelL1Muons = (\"l1extraParticles\")" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 8 }, { "cell_type": "code", "collapsed": false, "input": [ "histoD0 = ROOT.TH1F(\"histoD0\", \"d_0 distribution\", 100, 0, 100)\n", "histoD0Pass = ROOT.TH1F(\"histoD0Pass\", \"d_0 distribution for L2Muons passing cuts\", 100, 0, 100)" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 9 }, { "cell_type": "code", "collapsed": false, "input": [ "events = Events(\"/afs/cern.ch/user/d/demattia/work/29.0_ZEE+ZEE+DIGI+RECO+ALCAELE/step3.root\")\n", "for event in events:\n", " event.getByLabel(labelGen, handleGen)\n", " genMuons = handleGen.product()\n", " # print \"gen d0 =\", genMuons[0].vertex().x()\n", " # print \"gen pt =\", genMuons[0].pt()\n", "\n", " try:\n", " # L1Muons\n", " event.getByLabel(labelL1Muons, handleL1Muons)\n", " L1Muons = handleL1Muons.product()\n", " # L2Muons\n", " event.getByLabel(labelL2Muons, handleL2Muons)\n", " L2Muons = handleL2Muons.product()\n", "\n", " # print len(L2Muons)\n", " for L2Muon in L2Muons:\n", " # print \"reco pt =\", L2Muon.pt()\n", " histoD0.Fill(L2Muon.d0())\n", " if L2Muon.hitPattern().dtStationsWithValidHits() + L2Muon.hitPattern().cscStationsWithValidHits() >= 2:\n", " histoD0Pass.Fill(L2Muon.d0())\n", " # print L2Muons[0].d0()\n", " except:\n", " # print \"L2Muons not found\"\n", " pass\n", " # print \"\"\n" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "inverted\n", "inverted\n", "inverted\n", "inverted\n", "inverted\n", "not inverted\n", "inverted\n", "inverted\n", "not inverted\n", "inverted\n", "not inverted\n", "not inverted\n", "not inverted\n", "not inverted\n", "not inverted\n", "not inverted\n", "inverted\n", "inverted\n", "not inverted\n", "inverted\n", "not inverted\n", "not inverted\n", "inverted\n", "not inverted\n", "inverted\n", "not inverted\n", "not inverted\n", "not inverted\n", "not inverted\n", "not inverted\n", "inverted\n", "not inverted\n", "not inverted\n", "inverted\n", "inverted\n", "not inverted\n", "inverted\n", "not inverted\n", "inverted\n", "inverted\n", "inverted\n", "inverted\n", "not inverted\n", "inverted\n", "not inverted\n", "not inverted\n", "not inverted\n", "inverted\n", "inverted\n", "not inverted\n", "inverted\n", "inverted\n", "not inverted\n", "inverted\n", "inverted\n", "inverted\n", "not inverted\n", "inverted\n", "inverted\n", "not inverted\n", "not inverted\n", "inverted\n", "inverted\n", "inverted\n", "inverted\n", "not inverted\n", "not inverted\n", "not inverted\n", "not inverted\n", "inverted\n", "inverted\n", "inverted\n", "not inverted\n", "inverted\n", "not inverted\n", "inverted\n", "not inverted\n", "inverted\n", "inverted\n", "inverted\n", "inverted\n", "not inverted\n", "inverted\n", "not inverted\n", "not inverted\n", "not inverted\n", "not inverted\n", "inverted\n", "not inverted\n", "not inverted\n" ] } ], "prompt_number": 17 }, { "cell_type": "code", "collapsed": false, "input": [ "canvasPt = ROOT.TCanvas(\"canvasPt\")\n", "canvasPt = rootnotes.canvas(\"canvasPt\", (400, 400))\n", "histoD0.Draw()\n", "histoD0Pass.Draw(\"same\")\n", "histoD0Pass.SetLineColor(2)\n", "canvasPt" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stderr", "text": [ "TCanvas::Constructor:0: RuntimeWarning: Deleting canvas with same name: canvasPt\n" ] }, { "output_type": "pyout", "png": "iVBORw0KGgoAAAANSUhEUgAAAlQAAAI8CAIAAAC4XaJJAAAABmJLR0QAAAAAAAD5Q7t/AAAgAElE\nQVR4nO3d7dKqOqIuUDy17wvsdV8N7vvqBq7M8yN7UpmAiF8QyBi1atV8BSGg5jEhxMv9fi8AICf/\nb+8CAMDW/mfvAgDk4nK57F2EXDzt1BR+ANtxpWkDa75k6PYEIDvCD4DsCD8AsiP8AMiO8AMgO8IP\ngOwIP4D9VVV1+WP9s4anVFX1s6K9o2may5zqj6Zp1myn67rhzKx/1hoXN50AbONyeVjlxpm3slru\nuu56vb76rAVN09xut7ZtP4/Sqqr6vn+6WlmWXde9upGnJVw4zwM3uQPkbpSjm+n7/lFQLcTn9Xr9\nPOl1ewIcXlmWnzz918lXluVCCac9vV3XxclXlmVd1/EKn7dNhR/AIVVVdf9jofNwd6Fvs+u6obRt\n246ycHQxL/6zruuu65qmiVt7a/pUl7nmB7CRV6/5DZH2dkNnFIqPthPv/ekVtTWlijstH13YG3Vs\nxmfm0RXQ+CkL5Vxzza+4A7CJhSp3VC2PevlChIye0rbtQmX+qJtxtJ2F3sjl3Q3atp0eTrzZacln\nNzhsZ+G44tNS1/XCyXy0aHDOAS9+NwQ4rtmxHgtjQ6YW6sCXtrNcpOB6vS4P2lzYZvxnuKthtM5C\nNn/Y03vaa34Lgb+w9O0nfrL0WJvdZaeOJc2dOpZXl66svuIOw9GiNTX+KELquq7revka29MNLl9j\n6/v+vY7ZuFTDoS3f/BDv9I09Dk4bfgDHFfoAwyCR+PE1oTW6kNY0TdM0XdfF3Ylx0ox2MXQ/DivE\nG4w7G+NOyL7vUx50MyX8ANJS13XcxJle/1svDsuqqoZgWx9U8RbKsoz/bJombrq9Mf3KjhPTCD+A\ntIwi4ZOEuN1uw8RgsxfVnoqbfdPIjB95ox9yx8ai8ANIy0JErQmY6bDMvu9vt9v1en11FtBXw+lA\nPZ/CD+BUqqpa6CkNoz2/OEP0J+Is37gLVPgBnE2YD2U6yHNwu93WbOfVQPokwIbnLmwkbll+OKNb\njuG3MOB4eSzyJ0t/tNm3j+VHRTrWKdqlSD8qz+92mtqrdqZT9GvhOt/9fm/b9sO7HX5hVIDZ8Bt1\n88bh92FLMcfwAzixpmlGv5k3DHh5o7UUP2WaN/EjL22867q49bnQT/toTM2n3aQL92YeV4LHlVqR\nUivPPb0ipVaee3pFSq089+SLtFC8hWr50XRfs48/mjMsiPMpXhQ/JdzJ92iysfg+v1HaxRscLaoj\nsxk5OuRRFi48+OhkLiz9v3WernFE0zO75nz9ukh77XpWauW5p1ek1MpzT69IqZXnnnyRFoq3UE29\nFH73SfCUfyxUhsuBtKZJN5pp86VW4OzsoG88Zc15Hpy223PhpACc2Oh+g/6P+MFRcD79GYflMKvr\n+u0riI9+nOHRPNpFUZRl6ff8ABi7z/0uRBB+Y2EUHqO5WqbC7Giz67Rt+0bylWUZWm+PYizMRzN9\nPPy836u7mzrn7/m9MWf5r6VWpNTKU6RXpNTKU6RXpNTKUyRfpF2K99IIyTW/1ff5rwy+6tU9rjnP\nyb1RviLBDwCAqmkba86zbk8AsnPs8Ou67tWp6gDg2OF3vV73LgIAx3Pg8NPgA+A9/7N3Ad7UNE3f\n92VZfvhL9gBbulwuexeBojhu+N1ut/Cbwno+gQMx2nMDa75hHLLbMxzYgX41EYCkHK/lF6YSWJj5\nJnipbyHp72LTA0m5tABHcLDwC7+CsWZit4U8G93/eIAu+PhY0i8tQPIO1u0ZrvDp8ATgE0dq+YUO\nz3hW09E/3PwAwBpHmmiuqqqFGxvKshyycHlit2m3Z9In4XIZd3umXFrgsdRrm7NYc54P1vKbdnje\nbreiKOq61uwDYKVjfw3puu56vcZtvkDLD0hQ6rXNWfhVBwCYIfwAyM452+C6PYEEzdY2YRz7VFVV\nXxnKELb/aC87mi1Y13XhMta0wAuLRvyS+6qlwg/Yxmxt82iejelohllN09xut+UaryiKtm2TGhUY\nRu+PjnF0KuIyjxYtV9pnG+35kkfvp6RzDsjVdOKqNVkVJr1aXie+NzoRs/etDe25tm3DcV2v11Bj\nh1MRfswg/KRPVVUfHtRpw0/IAQdSVdWj3rxhEo/Q7xd3h8YTfYQHh5Wbpqn+Fm/w867F94Qh+o8W\nFX/uW6uqKg71+NDCj/l84cfs7me0fFyjpamfhFHxEi8t8NhsbROq4rquZ58yTOJfluWo6h7N7x+2\nMK3hwz/atg0bHG1n2O/01wIeFekTYS91Xdd1HQ5qWBQ/MhRmuoVQ/viJU2tq9XPWpMIPSNBC+E2F\nuIrDr23bkBDD0iHJ6roOj8TRFcdhWBo/fdhy/MTwlGG1352KafjdJ8E8BPb0XC1vXPitWir8gG18\nEn6j9dcsnT4SJ9z9T9iMmoxDjv7UbPiNzsBoadu2cdgvbHxNre4+P4D9TWvz5QEvy8M9lp97u90u\nl8vlcglXzsKmhhwNI00ul8vGd0eEMg/noSiKvu/jwwxXPcOi2+1mwAsALxiNLA3/rqrqfr+HGyfC\n47fb7Vv3Gq4RknjYXVmWIfyG0S6jkkwfeYmWH0BewljQUcMujA5tmibugdz+Holhj/F4zsvlcr1e\nR4X5MJWFH8D+hq7I2Pqnr7zvLfRtDk262+0Wt/P6vg+9nUMubnlffLgKGMo2HPvwM65FUVyv13jR\np2V7elXwiJaPqzDgBdjDbG2zXD8/GtIyGrRS/D1uZbr9YQzLMJJz9Ph9xWDLL5od8LKybMv3OdzX\n1eppz+z1LtObAQlKvbZJ3srrfOb2XLU09bej8IOzSL22OQtze87wzgPgtOEn5AB4xGhPALIj/ADI\njvADIDvCD4DsCD8AsiP8AMiO8AMgO8IPgOwIPwCyI/wAyM5ppzcztycAj5w2/IQcAI/o9gQgO8IP\ngOwIPwCyI/wAyI7wAyA7wg+A7Jz2VgeAQ2iaJvyjqqqqqqaPD/84hFDaUZm7ruu6bvr47HNH5+Hp\nxt9zOeX9cJfL0nGNli6vvL/LpYiLN/oTOI7Z2iaekWNY2jTN7XYL/27b9lEYpGYo9qiOjdeZrW+7\nrrter8urDes8PSFranXdngDJGZLvQLqumxZ7SKm2bcuyjB+JhVQry7Ku69ETR+t8y2m7PU1vBhxO\n0zTLfXoL/YfDorjbcHik+Gqf4azZcOr7vvhT8XZd1zTNNNVCIYuotNfrNTxx8P227/2Mlo9rtDT1\nkzAqXuKlBR6brW1CVRxaRWVZ3u/3tm2HR4qiaNs2XnMwPB7Wj4Xt3O/3YbOzT/yisIuh3TYq2P1+\nr+v60X7j1eIjHdYPm52ekEfW1OrnrEmFH5CghfCLY2MUJKGuj6v+UcYMCfdoUdjakDFDNH5LnF7x\n3ofCxB6dmWlkjlL/PgnFR4TfqqXCD9jGQvjFNf5sXT+q94eAHG1t1IqKG5T3KI2+flzDNmfDLwTz\nqDCxOKSHmIxTf/YkLJTnaZkNeAFIQqjlhwtgs67X6+VyuVwuo0tiTdOEx2cvvP10sGjYeFmWw0XH\nYnIUXddVVRUeHJV8WCEc/u12iztpm6YZ1n+08ffsPMq/qqrhwMK5e/qU2ZEsdV3HV3Hd6gAkaOFW\nh9BiG6Ir1IfDoqqqwr9HHYlhbEtYVJZl+DNsJOwo1LFDDTncLfDFSu/R6MK6rkeFGVZ+uvdhtTgj\nYst5saZW33O0Z/yCdV3X9/3lcnnvjpavfBEA2FFc9T2qBoeRnMP94MOivarBuq7jXYesGpI4PBi3\n/AbxSNQQciGx4pbM6DzEG/+03E87Rn9kGL0zPDIMcFp4Vlhn2s09snxchWt+wB5ma5tQFcfXtIq/\nr5+FRUObr67roWMwXi2+Sa74+5rfUGdOx1X+4hjj7Q9Fja/qTdecrjZ7YW9h0Wi15+VcdTQ/EI5w\n9ODTVyW8tB8eufADdvE0/GZTbTTifzD7eFmWceDtHn7xQRV/N29Ga8aH9qiF88XwS+ty19Pu4Lhp\nvLydhXVc8wN2kXpt80uh5/Nbqy072PRm4YBn7wsZhA7fYVzT5XI5ypR3ADlbWVdvVqXv/zUkns90\nNGhzahhWFNrRwyig0VFo+QEJSr22OYvDtPyGHuHb7bYwYGlYdL/fww0lQ2/y9MvC5bHR0u8fDwBp\nS+hryPqfqxiZXinU8gMSlHptcxaHafkFVVWFC36v3q2yZloEABgkFH7Fn95LMQbAT+0WfrPX25Zj\nLwzynI6ICcNeDPsEYKXdwi/0VY6SLPwK8KMBnyHeRr8UHFYe/VoVACzYLfxCaN1ut6Zpws/7xlN9\nDuvEd/JVVRVCLrT/wr2QIQv1lALwgqdzwPzO9NeHR1PaTOf/vE8m+JmdC3T5uArTmwF7SL22OYs1\n5zmJcbfx3N5feYpbHYAEpV7bnMWa83zOV0L4AQlKvbY5i4Pd5wcA29jzx2x/6tG8Zb52AXDa8BNy\nADyi2xOA7Ag/ALIj/ADIjvADIDvCD4DsCD8AsnPaWx0ADmH2d2yqqoqnbxzWiVfuum6Y6HG0cvz4\nZj/3Fso23eOjck4Na47OyfSnC75wUD+ZVXRvy8dVmNga2MNsbfOoch5m7Y9/AyCe/X9ak09/LaCY\n/GDAL4x+byA+zNGihcIsZNOrybWmVj9nTSr8gAQthF9ZlvUfww+Utm17/zvSZhNx2OywTtu2bdsO\n29nguIb9xuUc/gzHFR/UyLB0+Ef8iz3DI4On5Xle5heO7ziEH5CghfAbNYniB2dzbtSiip8Vp8tC\n3nzLULzpnyF9h6wa/RmbPd54g08Db7S1p+ucdsDL5YG9ywXwkXABLPx/aNvFrtfrcM0sVPQ/vexX\nVVUcYLMrxP/o+355zeklw2FpVVVf++ny9Vl6IMvHVWj5AXuYrW2G2njo0xvVz3HnYfGneRQ/Ujxo\nDhabXPCbPZyw39EPki90w44eD3+GBus04J8e1Jpa/Zw1qfADErQcfrG4l28Iv6EDcPjHKPzCytO0\n2ODQRscyemQ21KdPHP0Zh1+cpsWzjtw1h3zabk+AAwmRNiTEo/sfiqLo+364qWB2na7r7vd727ZD\nVGxzt8NwXSlOwZDTfd/3fT/bSRub7dIMhxMOeTgtn3d+Cj+A/Q0XtEJCXK/Xhfo9XDabXhuLRzZU\nVdU0zdO8+ZbZ5CuiK4L3P5ce1xfpp5kt/AASMmTe9XqdLo2v6j3KhrjVuDC65IuGkoySr+u6+Mb2\n2+0Wrzzc0h6vX0xadSHRRw+6yX3e8nEVrvkBe5itbUJVHA/iiEe43P8e+j+6kW50ze9Ro2qbWx1m\nw2Uo7TR04j/j2/tGJ2R4JP7HcpHW1OrnrEmFH5CgleF3j2r8cLv6NE5GA0CGJ04HfP40+e6PEzcs\nHUXj9B7E2ZKP4i3exZob/tbU6pf7g7FGh3a5LB3XaOnyyvu7XIq4eKM/gePYsrYJXY7b7OuLFoq9\n/ojWnOe06/13CT8gQanXNmex5jwb8AJAdk77k0aPZjLztQuA04afkAPgEd2eAGRH+AGQHeEHQHaE\nHwDZEX4AZEf4AZCd097qAJCgR7cgszHhB7AR9x/P2mXWN92eAGRH+AGQndN2e5rbEyApSV3vPG34\nCTmApDyqlncJRd2eAGRH+AGQHeEHQHYOGX5VVV3+qKpq7+IAcDDHG/ASLo2WZVlVVdd1fd9fLpe2\nbaUgACsdLPyapimKoizLruvCI13XXa/XpmmGRwBg2Q6Tynyiqqq+70dlDm3B+MHlyXJGS3eZWecF\nl0sRF2/0J8DB7VIJp13vryP8AI7L3J7vCJf66rreuyAAHMZRw6/rujDas+/7uq7DtUAAWOOo4VcU\nRVmW4R+322062uXy2L0oistl+O9Hje3pfn+zH8hM9OH9v/9Iz0INnEjdmPblrhXCaM+iKOK7HZ70\nIG9yFe1rVxZd84OYT8TpuOb3jqqqwgU/tzoAsNLhw6/4M+ZF+AGw0sHCb7aDWOwB8JKDhV8Y5DIa\n23m73aYPAsAjBxvwMgxvqes6zO0Zki+e8Kww4AVOzCfidMzwssqQf4PpfX7CD07LJ+J0hN9rQlNv\n9scchB+clk/E6Qi/rxF+cFo+EafjPj8A2MLBfs9vvUdT5pyypQuQvqRmeTxt+Ak5gKQ8qpZ3CUXd\nngBkR/gBkB3hB0B2hB8A2RF+AGRH+AGQHeEHQHaEHwDZEX4AZEf4AZCd005vZm5PgKSY23MLQg4g\nKeb2BIA9CT8AsiP8AMiO8AMgO8IPgOwIPwCyI/wAyI7wAyA7wg+A7Jx2hhfTmwEkxfRmWxByAEkx\nvRkA7En4AZAd4QdAdoQfANkRfgBkR/gBkB3hB0B2hB8A2RF+AGRH+AGQndNOb2ZuT4CkmNtzC0IO\nICnm9gSAPQk/ALKzc/hVVXWJNE3z9CmXOWueCADBZcdrY6GftyzLqqqKorjdbuHPruuePmtk9KzL\nZfG4LpciXjr680tGZXhSpMUNbVBaOAyfiNN5v3r8wG4tvxB4dV13Xdc0TdM04eD7vl8Iv7Coruv7\n35bzEgBiu4Vf3/dFUYy6K+u6Lv4k3KywKAQnALxnt1sd2rZ9o7km/AD43J7X/KbC9by2bR/FW1ih\nrutwgbB4cI3QNT84LZ+I09nlml9C4dc0TYi0hSINo13Ksiz+9J1OnyL84LR8Ik5nl/BLZYaXruue\nJt/QwovXqaqq7/uqqkbtv4UpA+5/L/W5AfiupGYym5VEy29o8y10eC4IZ/mFZpaWHxyXT8Tp5HWr\nw6Cqqk+Sr/jTBepuBwBW2rnbM3RaFuahBmBDe7b8QvKVZbky+ZqmmZ3JLMSn+x8AWGm38GuaJiTf\n+u7KeBa0eDvFn55PAFhjt/ALGdb3/cIs1aGpNzTpqqoKIRfW6bpuuF7ogh8A6+0TfstZtTy9WZgC\n7Xa7Xa/Xl3pNASBI4laHNyzPc+ZWBzgtn4jTyX2Gly8SfnBaPhGnk+l9fgCwMeEHQHZSmdvz6x7N\nLHfKbl6A9CU14edpw0/IASTlUbW8Syjq9gQgO8IPgOwIPwCyI/wAyI7wAyA7wg+A7Ag/ALIj/ADI\njvADIDvCD4DsnHZ6M3N7AiTF3J5bEHIASTG3JwDsSfgBkB3hB0B2hB8A2RF+AGRH+AGQHeEHQHaE\nHwDZEX4AZOe0M7yY3gwgKaY324KQA0iK6c0AYE/CD4DsCD8AsiP8AMiO8AMgO8IPgOwIPwCyI/wA\nyI7wAyA7wg+A7Jx2ejNzewIkxdyeWxByAEkxtycA7En4AZAd4QdAdg4ZflVVXSJN0+xdIgCO5HK4\ngSHh0mhZllVVFUVxu93Cn13XxessHdflUsRLR39+r5xxGZ4UaXFDG5QWDsMn4nTerx4/2emxwq+q\nqr7v67qOW3shDtu2DXFYCD84MZ+I09kl/A7W7dn3fVEUo37Ouq6LoohbfgCw4GD3+bVtK+QA+NDB\nuj1n6faEjPhEnM4u3Z4Ha/lNDV2gQ/IFC1MG3P9eei+KIvpz9LSFl2S6i3jl0WYP9umcnr1H5+Hp\n1AyjemrlZpd9azub7WW0qd0r608OLbVjGXn70F56YoLv5MRel6RmMpt17JZf13XX67WYRNRrLb/H\ni5a386Rt963vp7t8z12/0+XyLG8n8XPyxb3Ez02hpfLJoe1+LC+9qdaX8KUnJvhO3v11+YABL69p\nmiYkX9u2e5cFgCM5ardnuOeh+PtSHwCsccjwG5Lv0H22AOzleN2eIfnKspR8ALznYOHXNE1IPnf7\nAfC2g432XBg+G895ZrTnp4z2/O5eUhuJZ7Tnq5v9ZOVv7XT9plJ4j73CaM8nllt72oIArHSwlt9K\nWn6f0vL77l5S+1au5ffqZj9Z+Vs7Xb+pFN5jr9DyA4AtCD8AsnPI+/zWeDQ05pTdvADpS2rCz9OG\nn5ADSMqjanmXUNTtCUB2hB8A2RF+AGRH+AGQHeEHQHaEHwDZEX4AZEf4AZAd4QdAdk47w4vpzQCS\nYnqzLQg5gKSY3gwA9iT8AMiO8AMgO8IPgOwIPwCyI/wAyI7wAyA7wg+A7Ag/ALIj/ADIzmmnNzO3\nJ0BSzO25BSEHkBRzewLAnoQfANkRfgBkR/gBkB3hB0B2hB8A2RF+AGRH+AGQHeEHQHaEHwDZOe30\nZub2BEiKuT23IOQAkmJuTwDYk/ADIDuphF9VVStbvpc5TdP8uIAAnEcS4VdVVd/3n2yh67ovlQWA\n89s//F5KvhBydV3f/yb8AFhvz/BrmuZyufR9X5blyqeEkKuq6nelAuD09gy/kGRt265vtwk/AD63\nc8vvfr+/lGShgzQ0GQNBCMCrLoncDB6Gej4tzDAiNPSUDhcLR0+8XBaP63IpHi39e9HydkZLxyuP\n9rKw02Xf2s6PdrpcnuXtJH5OvriX+LnbvILLPjm03Y/lpTfV+hK+9MQE38m7vy4feFJj/8aRZngZ\nekfj0xTGy1RVNeo7Xbhx4j5Zmsg3AIBzSGoms1kHa/mtfK6W36e0/L67l9S+lWv5vbrZT1b+1k7X\nbyqF99grdmn57X+rw+dCF6i7HQBY6QzhBwAvOVL4hUGe05nMwrAXwz4BWOlI4Rfi7Xa7xQ+GLFx/\nmzwAJB1+oak3NOmqqgohF9p/XddVVRWy0AU/ANZLOvymuq6r67ooitvtdr1ew9RoiQxYBeAoUrnV\n4VXL85y51eFTbnX47l5SG4buVodXN/vJyt/a6fpNpfAee8UutzocNfyWCb9PCb/v7iW1ikn4vbrZ\nT1b+1k7XbyqF99gr3OcHAFsQfgBk50hze77k0cxyp+zmBUhfUhN+njb8hBxAUh5Vy7uEom5PALIj\n/ADIjvADIDvCD4DsCD8AsiP8AMiO8AMgO8IPgOwIPwCyc9oZXkxvBpAU05ttQcgBJMX0ZgCwJ+EH\nQHaEHwDZEX4AZEf4AZAd4QdAdoQfANkRfgBkR/gBkB3hB0B2Tju9mbk9AZJibs8tCDmApJjbEwD2\nJPwAyI7wAyA7wg+A7Ag/ALIj/ADIjvADIDvCD4DsCD8AsiP8AMjOaac3M7cnQFLM7bkFIQeQFHN7\nAsCehB8A2RF+AGTnJOFXVVVSl1IBSNkZwq+qqr7v9y4FAIdx+PCTfAC86sDh1zTN5XLp+74sy73L\nAsCRHDj8uq4riqJt2/APAFjpctybwbuuq6oq/DuMdhmO5XJZPK7LpXi09O9Fy9sZLR2vPNrLwk6X\nfWs7P9rpcnmWt5P4OfniXuLnbvMKLvvk0HY/lpfeVOtL+NITE3wn7/66fOBJjf2jnR43/GLT8FtY\n+V4Uo8X3+E0z2vKjNedW/tqHcHmzy761l+U/V+7xpe1My7C+eAvbedvTvbzk8Xts7ROnPinSJ6/L\n+mN5qfALX7CWy/PS0kcrv30S1jx3TQE+3E7x7rGMvF34Z9uJ/5552uZJlOX0ZgvfMkaPT9p2S+t/\n8V6LhcIvv0U++dS95Cvv1MnZfrjCcvGebucrvvXhfGk73/qu86qV740P343r34FvfyLeLuGHb6pd\n3i0rt/O7L9CLO13oGLtcLts3wg58zQ8A3iP8AMiO8AMgO8IPgOwIPwCyI/wAyI7wAyA7J7nP7xy3\n6gOwDS0/ALIj/ADIzkm6PaceTe+pgxRgF0O1fH82A/MGTht+Qg4gKfFPCDyZNvn3dHsCkB3hB0B2\nhB8A2RF+AGRH+AGQHeEHQHaEHwDZEX4AZEf4AZCd087wYnozgKSY3mwLQg4gKaY3A4A9CT8AsiP8\nAMiO8AMgO8IPgOwIPwCyI/wAyI7wAyA7wg+A7Ag/ALJz2unNzO0JkBRze25ByAEkxdyeALAn4QdA\ndoQfANkRfgBkR/gBkB3hB0B2hB8A2RF+AGRH+AGQndPO8GJ6M4CkmN5sC0IOICmmNwOAPQk/ALIj\n/ADIzs7h1zTN5Y+qqrque/qUy5ymaX5eVgDOYs/wq6rqdrsVRVHXdVmWfd9fr9c1+Tf13rMAyNNu\noz27ruv7vizLIbe6rrter9frdWGgZli5rmtNPQDetlvLL6RXnGFVVZVlufysEH5VVf2sXACc327h\n1/d9MYmx8OdCq074AfC5tEZ7Pk21EJmjYTIbFAyAM0kr/IKno1dut1tZlqGPtO/73afJAeBYjjS9\n2RCK8YiYqqr6vp/eJrGQiNNp5cyFBvBdcTV7X/xzF0cKv6qqpinVdd3lcgndobGlPPt7WjkAvu6v\nanZU65rb8ytC/6db/QBYKa3wM5gTgA3sFn6zzbXl8AuDPKc3QszeNQEAjyR0k3vxLMbC42FGtNF2\nnt4dDwCD3cIvzOcyDNQMrbqiKNq2HdYJDw5ZOEwBE9p/XdcNs4O64AfAente8+u6bpjPOmRY27bL\nvZdd19V1XRTF7Xa7Xq9hdlBDNwF4SRKD/kMb7tWnFI87SC/LNzNcLsXKo/57zfFmR9t56c/VO33N\n23vZ5liWn7iw8ifbedu3tvO7MvzoSLd5J6/f7O9eiIX32HIZvnWKXtrLt7bzrTfV987Ykxr7N5II\nv68Tfi+vKfx+t53flUH4fUj4PX3u+r18UIZdwi+tWx0AYAPCD4DsHGl6s5c8mi/nlN28AOkbqmVz\ne/6QkANIyj26wmpuTwDYmvADIDvCD4DsCD8AsiP8AMiO8AMgO8IPgOwIPwCyI/wAyM5pZ3gxvRlA\nUkxvtgUhB5AU05sBwJ6EHwDZEX4AZEf4AZAd4QdAdoQfANkRfgBkR/gBkB3hB0B2hB8A2Tnt9Gbm\n9gRIirk9tyDkAJJibk8A2JPwAyA7wg+A7Ag/ALIj/ADIjvADIDvCD4DsCD84jN3vCx5JrTyFIrGa\n8AMgO6ed4cX0ZgBJMb3ZFoQcQFJMbwYAexJ+AGRH+AGQnUOGX13Xlz+qqvrvf//7rS3v0vW8vNO3\nlyZ4LLvs9JOlv9jsj8rzyU4T9KNjSe08HO5YzvQeO174VVX1v//7v0VR/Pvf/y7Lsu/7f/7554v5\nB8DpHWy053//+9++78uy7LpueOSff/75559/DO8EYKWDtfyGNt/wyL/+9UGBkf8AAAd6SURBVK+y\nLPcrEQDHc7Dw6/u+KIp//etf8YMh/Oq63qdMABzNwcJvlpYfAC85Q/gFoVEIAE8dbMDLSveiKL4x\nEn20nelmF1Z+unTlE5dLu7CXl07CJ4f29mZHFlZ+6RQ9fdWWd7qytJ8sfe8Fnf453dR7O/3k1X9p\nswtFWj7Stz8RL+30R6doectvf14Wlr5dF7200y9Wj7s4Z/gVz0Z+fmtg6PJ23t7LNuNWPyn8Noe2\nsHIKpygFPyrhj179t7eT4Mud2in64nYSrx6/5TzdngCw0hnCL1ztM+wFgJUOFn4h4UbzuQg/AF5y\nOdbEKGE+l3iGl+LPhdNjHQgAOzpYyy/M59L3fZjPOsxwXRTFf/7zn6IomqaJJ7yOA3JjVVU9GssU\nFg2F3KAYg6Zp9i3Pyt1tXKSg67pHu9uyPKEYU+kUad930ezJmd3vxu+iuPK5XC6zlc9en/19P2if\nVIZr3nvvux/QqIfzP//5T/xgXdfDv9u23bF400Xh8bIs40L+qBjxvobpb8qy3Ks8j3Y3eo02LtJ0\nv2vK/LtiPOq936tIwzsnfhfVdb1XeWZPzvS12+tVS+SNPZyT3cvzSWXYtu3T994nDhl+Qci8IJym\n+N0/nLiNSxVXXqNF4fWLX7xhYrYfFWO2kho+A1uWZ9jd9DWKH9m4SKO9FJPw2748T9+0GxdpWp7R\nI3u9ZOkUaU3ls/1nf/Zr7sbl+bAyfPre+9CBwy8WTtzoq014cLMyxK2r2Rdp/YOfW/OG27I89wcv\nx2h3GxcpGFL5aZWxQXmebnzLIoUzM/pYhe/pw4O7vGTTfcWF3LhI03p8urtEPvubvWqfV4bTrxT3\nB6f6bQe75vdIGPA57fQviuLL3cSPhV7+tm1nu/vDg5sNSW3bdnmm743LE/Y4vL8frVDsMWr3druV\nZTl9n+xVnoU9blyk2+1WTD5WTdN0XRce3OsUxYUJBRgKuXuRplIoUvx6bVCezyvDsM7oUxn+/NZg\njpOE36zNxkoETdPc7/flnU6XhnfA18fmVFU1rc2nddlm5ZkV9j4K6Y2LFC7FL2x8y/KEbVZVFQ8E\nmO5oyyINWw6lejTaZbPyjIS39L6nKJyT2+02bHyI5L2KNBXH3gbl+bwyXCjGt6ZxPnP4BZuN+Vx+\npXccehoM1Vb8HXAXwwiuvu/ruh4Ktn2Rwq6HKzQj25cn7PF2u4VfbA4PXq/XHU9RURRVVV2v177v\n+76/3W5xHu/7rp798rRLkcJb6Hq9hjf27Xar6zqRszQtwwbl+VZl+NMGzPnDL3HbNE+7rgvfkZc7\nHjcrz1Czx9+XH/lRkcI5iXvMVvrdKQqnIlzSCL3EoVZ9epZ+UaSwx5B5o0FS1+t1+bnbvItCC2Dl\ndY2fFmkoQwpv7PAaxXcXNE2zprW0cVfZcgE2+JUe4bezDb6FNU0TaqtHTZyNy1NVVVyzX6/X5Z3+\nqEjhnLyx8d+donBa4tq8qqpQly1X8T991dq2HWqlpmlC/b5jeYaSFK/8hPXvihSiJQzNGL2xdynS\n8NIMPeehJfr0ibu3UOMCbHCJVPhtZK9vVVVVhTZfXIXtWJ7YULMP17o22/VQdXaRsGj4dwqnqPhT\njNkhXRvsd7rTeNDBjqdodjDO7CPblGR0OS3+yrJ9kUbjLVP77K8vwE/z+Mzht/vnc+rRa/mjQlZV\nFSrNRxefNy7Pox2tuSDx3SINV9euf4QqrO/7+BrbZuV5SWpF2rE8j3ax+yna6409CENOQmO0+jPd\nVbyv3U/RcgG2KMZXbpjY3cJ9frtM8jJ7btc/+BWzN67tWJ7ZLY/ufN+sSOFWkJGwo7qu97qJLalT\ntFyeve4WnS3Dmr1v/KqNblPbskjhRszRg6PKcPd30dMHwwkcvcSzN/+9X7CvbGV3sydlx3Sf3fU0\njL97z2ZsOp3K1Jblua+bdGbjIsVm30KZn6LZd1EKL9ny99qNizQ6IXEZhj1uWaTpvqbv7S3L83Zl\nOH3i7Kl+v2Bf2UoKhoZOfH/3Ls2+++PcDY+HtsVQyN8VYNb07bVBee5/z9QX7262bt2mSNPiTb8u\n7HKKRm/j2TjcpkizH6vdX7Knn+4UXrUdz9Lsvmbn9tysLlpZyHiFUWfM066slwv2rQ2lYDRAaK/k\nuz9+vUfjLb/4Qi7sZWR2EsKflmehYNNvmhsXabTf6e6cotHHavfy3Fd06uz+qu37Rprua1oZblae\nRy/WmgKMBql+t2F6qvALdsy89VIr5Mbladv26R6doqROUWrlWcmr9pV1fmrNGftFIQ/2Y7YA8Lkz\n3+oAALOEHwDZEX4AZEf4AZAd4QdAdoQfANkRfgBkR/gBkB3hB0B2hB8A2RF+AGRH+AGQHeEHQHaE\nHwDZEX4AZEf4AZAd4QdAdoQfANkRfgBkR/gBkB3hB0B2hB8A2RF+AGRH+AGQHeEHQHaEHwDZEX4A\nZEf4AZAd4QdAdoQfANkRfgBkR/gBkB3hB0B2hB8A2RF+AGRH+AGQHeEHQHaEHwDZEX4AZEf4AZAd\n4QdAdoQfANkRfgBkR/gBkB3hB0B2hB8A2fn/RAPaSEzlpm8AAAAASUVORK5CYII=\n", "prompt_number": 50, "text": [ "" ] } ], "prompt_number": 50 }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [] } ], "metadata": {} } ] }