diff --git a/debug/dis_hvq_Zj.ipynb b/debug/dis_hvq_Zj.ipynb index 18292c0cd..5f2a4e080 100644 --- a/debug/dis_hvq_Zj.ipynb +++ b/debug/dis_hvq_Zj.ipynb @@ -47,7 +47,7 @@ }, { "cell_type": "code", - "execution_count": 2, + "execution_count": 20, "id": "f4126da1-c09a-45e9-860f-342c333ec9af", "metadata": {}, "outputs": [ @@ -62,282 +62,116 @@ ] }, { - "name": "stdout", + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAANwAAAEZCAMAAADL1rsAAAAAIGNIUk0AAHomAACAhAAA+gAAAIDoAAB1MAAA6mAAADqYAAAXcJy6UTwAAAGhUExURf///wAAAAAAAAAAAEtLSwAAABsbG4qKiktLS2pqaq+vr46OjgAAAAAAAAAAAAAAAD09PQAAAAAAAAAAAGpqagAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAIyMjAAAAAAAAAAAAAAAAAAAAK6urpKSknd3dwAA/wAA/wAA/2dnZwAA/wAA/wAAAAAAAAAA/wAA/wAA/wAA//8AAP8AAAAA/wAA/wAA/wAA/wAA//8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAOwAAwMAydkAFo4AKv8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAEdHR1VVVRMTEx4eHkFBQTo6OrKysnl5ec/Pz4ODg01NTY2NjWdnZ3d3d93d3f///5mZmaqqqiIiImZmZlVVVYiIiO7u7ru7u8zMzAAAAERERDIyMhkZGXBwcK+vryAgICgoKF9fXxERESQkJBsbG1JSUgAA//8AADMzMwwMDFdXV6KioggICDg4OGdnZyyq7BkAAABpdFJOUwAAEXetRKA/5tr+mSHHdqPW9dbj9TP6tmb+XI4/7feIVd/vmMfO1kSIu6iZqqqZMxHdzDPdVe5mIneIVZkiZhHud6rMRK271vFchejf8XPAqMfz312wpuCQ/HrPes3y69fa9I/+Iexc0rEUZ+kAAAABYktHRACIBR1IAAAACXBIWXMAAABkAAAAZAAPlsXdAAAAB3RJTUUH6QEPDSUsSm/2WwAACttJREFUeNrt3YmX2sYdB3DPunbTFonAxhAvTtMjiUR3m10WOd4VQlpJlPRu2qT3fXeN48R2lLRNT9P7r+4MApZDEvOTGOkn3vzew2tg4L3PfmdGQmhW165dI6WvvesfiXz8WvlxN25e3h99dBdxe9efeTAajd7+2O7hJqGN7l8+fOfjO4abhfbuw0ePR6NP4MdVFEVRiaoo1U0tr0J7Mho9fvDsDfw4Qmp1+s++mtxoObT3nrt5K7IZOly1QUizkthkNbTn92IaosOR21UlycYZGlKc0jiIf5I7NKS4ahCHg4SGE9c62K9HPgELDSVObRH19vqYA4eGEXfQov80miuPpggNH06tB3QD16o1lIUH04WGDxdRaUPDj8sQGnZcptBQ47KGhhiXPTSsuK2EhhO3pdAQ4rYXGjrcNkPDhdtyaJhwWw8NDU5EaDnhgmnFPS8mtDxwU1UQ/nf9eWGhAXGV5r7aUg44W89sM9nktqoTGBoQVyW1O0TdB9JWdAs8saEBcZVKgwpfgOAibrMSHRoQR1p1QupN3tYkiNaxWw6hQXFKk1QaFd7WcbIgn9CguFpTrfPPJ0Gs7pN5hAbEVWrkDv/7xvVIOqUEeYQGxLUUzoYhLk5Hn8kjNBhOrddVvpYhjiTocggNhnvxU5C3DRJuJMghtBnuxos89engM8pnuVqyChLGHAlyCG2Ge2nEU++z3UQ2D3BVUq9c2wsTiXv5CU99EO7b//4Pf+RpHST1yjxxr1zy1IfTTy5/+jNP6+Qxlw9N02Hd8i9s34KrWyaOuZxwbYr73COe+msQ/O3vXC1ZFTzmtHa7fXh4xL0p+PyrkOm74DHXPiYnmibqk3gQFDrm2uEPUYcZkvYthduE4wrdErQ1oTiSOF+WHhf/aVVkHXeenrZPBXfLWJ3g3I60hTsCj1tG6oKuWFx78Y7Ig7IRuqBriNXlldzsIPPyYUvBOj033PKcOR1uYnWHOeKuvge5mkkE6bROu3N6mtuYiysxOv3u2kOFfD+XTac/fY3uF69RlqaSAnGZdPc0cnJI7q2bseCy6I7p7e7Zeh/Uz7Hgso47/WTp7hnBlFzmcbd4h86T7c7To/VWxX3hH60zez2rv/G1y70yYp4sGhepM23iGJtfeu9sCafFNCvyVI11Xf+CkJ61+ZX68r3zmGaFnoeypnO9SXgba1mj6THNij3JZlVnm4T4m4fcVbFjXEiTW9M51sDkGHLzOuvoR+324UnM00WfHrWaXZdnyM1L1xOfLhq3puMZcvM6R46b6syZzTMBY07DjpvorHEvzUvx46juC+PxEPiiM03TjmPnSTw4Qm3jMWQLQG2v6fph/DyJB2cx29gEvWZDh0SDMye28YWzizjHNSY6yDaAaEccjRDgaA08insd8oqTNkcjHDj6ecD84vhLfE3PdF3T9MMS4Wh9+Stcn81POvq5rp+fcDRFhOM88sA1leDD8em4phKEOD4dz1SCEcelKy0uWafpdKLkmidx4hJ1nbvaOd88iRSXoOPacOPGJeh2ABev6+wALlb3dBdw67pj/YjOlPybb8y4Nd3pXU0716DvghS3qoNOJbhxK7rObuGWdac7hlvS7Rwu1B0ftY/0wx0bczNd55hoGnieLAGO6aAb7vLgSPerX9tdHPn6N3YYN4AdZZc4NCVxEoewJE7iEJbESRzCkjiJQ1gSJ3EIS+IkDmFJnMQhLImTOIQlcRKHsCRO4hCWxJUSZ5vmG980TebjWxxfJpwbLmHyuBfHlwnnhLgB9+L4MuHClXUG/+L4UuG6DOemWBxfBhzxw0WD0MXx5cDZ43E41GCL48uBc741np6XmGLIYceRN96c2kCL40uCe+vbGV6MHbe7u18SJ3E4S+IkDmFJnMQhLImTOIQlcRKHsCRO4hCWxEkcwpI4iUNYEidxCEviJA5hSVxZcW99Z3dx/TdfT/FFf0lwjuH46XVYcc4g/IPcw34GHVbc0DImJmtA0uuw4nzH8Qf0p8FOjkqrw4qzfWfgsb87zu6k1aHEsXOFXMM0etb0j7al1GHE+UN2+mHfNaz5tTXS6TDiDNYpyfJ5v6l0GHGWS0w3PIM0mw4jzvEt1+ubRn/lUbCuYFzP7C2dSOkMDZaXbY2H9uplNeC6YnGDISHe4hV5Bp5jsVNi+35Ea7CuWJxnh5vpq6Jz/5BO/64b1RyqKxY37q9m1DVs23RsI/q0WKCuWNwFzWg+JzqeRzPrW75hxJ7yC9MVi7Nt258POdscGOxO10t4RbQuZqVP4ZuCqyHX9WmndDadhx6li1vpUyjO7E33jKd3h445ID0j+dJK67rYlT6F4txezw0ljsn+4/oXhmFtupDGmi52pU/h3TIs0zTZ7iTdenNcEWtVF7vSBwmOzpk2uwKdx3XdoRVd7EofJDjHd+mHbzqh8DZfzi5mpQ8CnON22UbOv7B8n/dizlPdvA9Hz7AIcJ7pTzbkFmRR6kRnz4KOWemDAEc/b3suga7+ozqbrfhMKgS4rtFlm+CuD3uZ891wrS5enMN+9QPDvLA96MXT7fHGy2AWjBtOTM5guHHTHW1LvgxmwTjD6U4+3figq1mysqaXwUSMcy3S81KtAZ9dBzNpL7voCcUa2n7Xgo63afW85MtgFojzDY92xq5peOm/pOp/7/s/wImjnwIISfM3FxbL+WH8r6ZA3NCOOPaaQhd/5KFAnOObrtV3DfgSaV5dMThnSD8F0PnOG/tuVluCrhjcwHI81hv7wF0uoK6g5Iy+w7ojG3MCdQWNuZ4xcF3HNsA7JiBd/jjHYp+9up6RcOx1S7r8ca5lT3ZI0uxyAXX547rD8FCJ1cv8Vpt0BYw5y3KsAf0Qt+33XdeJxwWzmj9iXhiewXcQL5tOMI6ZmGpyu+JxHXvNrhOLW5BNpWEJiC1KJxIXLMnmVsJ97DWrTnRyqzfbN03+Y68ZdeJwEblNshtsYVeZUycQRyJ1C9OKcJ0wXBBzI6JxCzpRuACDThwuViccN9cJxMWMuRx13LhXr++BbAm9UjxuqqO4H/2Yp34SBD/9GVdLVkHSjJIDLtRR3Esjnnqf7fv+48FjrsajxPGWB26io7iXn/DUB+Gu/T//9W+e1kGxYy7U9SjulUue+jDE/ee//+NpXfSYm+gccd2y2DE3KdCE8nPIhFLwmJvi+BoCNwWJY+4Z0HvlgANWkDTmRvcvb95Ahas099WWcsCLS+iV775NR+WDHOLjx1VJ7Q5R97lx8XtgT955eHl/lEN8gOQqDSp8gfeN43RB8IvR6PGjh3nEBxhzrToh9SY3Lojrlc8/yzYokPgqiqKoRFWUqjCc0iSVRoX7nWN6Jftx6+Zz78Hiq9FfLNlXYTYIrtZU67zzSYJuUnvA+Kp0SDT5f7FgXKVG7sDeO6JnLmy/YfHdripgGwDXUsBvHntQNixIfEoD0mmgOLVeh3b5uWd6UC9it4s7vmogEpe61r8IWSq++FoH+3WMuM21OT61RdTbIsecyNoQ30GL/tPg3sQiw5Gk+NR6QEd7q9YAz2hocPBtX6lwJM2uS4lwW44PG45sMz6EuO3FhxJHthQfVtxW4sOLI9njQ43LGh9yHMkUH35chvjKgCNp4ysJLl18pcGRFPGVCQeOr1w4AouvdDhIfCXEEe74yomLiO/WDuHIany/3C3cUny/+vWu4chVfL+JerLsuFl8v91NHGHx/S7y8f8DsDsJ60ziNo0AAAAldEVYdGRhdGU6Y3JlYXRlADIwMjUtMDEtMTVUMTM6Mzc6NDQrMDA6MDBG9t0WAAAAJXRFWHRkYXRlOm1vZGlmeQAyMDI1LTAxLTE1VDEzOjM3OjQ0KzAwOjAwN6tlqgAAACh0RVh0ZGF0ZTp0aW1lc3RhbXAAMjAyNS0wMS0xNVQxMzozNzo0NCswMDowMGC+RHUAAAATdEVYdHBkZjpWZXJzaW9uAFBERi0xLjVTNlrBAAAASnRFWHRzaWduYXR1cmUAOTVhZTg3NDc2MmZlMzljOGUzZmRlZTc2YjRjM2Q1ZDA4NTE2NTVjMGU0Y2ViNWMwMmIwMmQ2YTQ5ZjNlYTg0ZVWPa50AAAAASUVORK5CYII=", + "text/plain": [ + "" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "fd = FeynmanDiagram().add(\n", + " v3 := Vertex().with_xy(5,+5).with_shape('blob'),\n", + " v4 := Vertex().with_xy(5,-5).with_shape('blob'),\n", + " Leg(pdgid=2212,label=\"$p$\",target=v3,sense='incoming').with_xy(0,5),\n", + " Leg(pdgid=2212,label=\"$p$\",target=v4).with_incoming().with_xy(0,-5).with_style(\"label-side: right\"),\n", + "\n", + " Leg(pdgid=2,type=\"remnant\",label=\"$X$\",sense='outgoing').with_target(v3).with_xy(10,+7),\n", + " Leg(pdgid=-2,type=\"remnant\",label=\"$X$\").with_target(v4).with_outgoing().with_xy(10,-7),\n", + " \n", + " v1 := Vertex().with_xy(8,2),\n", + " v2 := Vertex().with_xy(8,-2),\n", + " Propagator(pdgid=1,label=r\"$q$\").connect(v1,v2).with_color('red'),\n", + " Propagator(pdgid=1,label=r\"$q$\").connect(v3,v1).with_color('red'),\n", + " Propagator(pdgid=21).connect(v4,v2).with_color('red'),\n", + " \n", + " v5 := Vertex().with_xy(11,5),\n", + " Propagator(pdgid=22).connect(v1,v5).with_color('blue'),\n", + " #Propagator(pdgid=-24).connect(v5,v6),\n", + " #Leg(pdgid=11,sense='outgoing').with_target(v6).with_xy(20,5).with_color('blue'),\n", + " #Leg(pdgid=-12,sense='outgoing').with_target(v6).with_xy(20,3),\n", + " #Leg(pdgid=-3,sense='outgoing').with_target(v5).with_xy(20,1),\n", + "\n", + " v7 := Vertex().with_xy(11,-5),\n", + " Propagator(pdgid=1,label=r\"$q$\").connect(v2,v7).with_color('red'),\n", + " #Propagator(pdgid=24).connect(v7,v8),\n", + " #Leg(pdgid=-11,sense='outgoing').with_target(v8).with_xy(20,-5).with_color('blue'),\n", + " #Leg(pdgid=12,sense='outgoing').with_target(v8).with_xy(20,-3),\n", + " #Leg(pdgid=3,sense='outgoing').with_target(v7).with_xy(20,-1),\n", + " \n", + ")\n", + "scale_positions(fd,0.5)\n", + "fd.render(auto_position=False,render=\"tikz\",file=\"dir\")" + ] + }, + { + "cell_type": "code", + "execution_count": 28, + "id": "c9c33e5a-3704-4a44-b474-366661a25c34", + "metadata": {}, + "outputs": [ + { + "name": "stderr", "output_type": "stream", "text": [ - "This is LuaHBTeX, Version 1.16.0 (TeX Live 2023 Gentoo Linux) \n", - " system commands enabled.\n", - "(/tmp/tmpgaexcn32/hvq.tex\n", - "LaTeX2e <2023-11-01> patch level 1\n", - "L3 programming layer <2024-01-22>\n", - " (/usr/share/texmf-dist/tex/latex/standalone/standalone.cls\n", - "Document Class: standalone 2022/10/10 v1.3b Class to compile TeX sub-files stan\n", - "dalone\n", - "(/usr/share/texmf-dist/tex/latex/tools/shellesc.sty)\n", - "(/usr/share/texmf-dist/tex/generic/iftex/ifluatex.sty\n", - "(/usr/share/texmf-dist/tex/generic/iftex/iftex.sty))\n", - "(/usr/share/texmf-dist/tex/latex/xkeyval/xkeyval.sty\n", - "(/usr/share/texmf-dist/tex/generic/xkeyval/xkeyval.tex\n", - "(/usr/share/texmf-dist/tex/generic/xkeyval/xkvutils.tex\n", - "(/usr/share/texmf-dist/tex/generic/xkeyval/keyval.tex))))\n", - "(/usr/share/texmf-dist/tex/latex/standalone/standalone.cfg)\n", - "(/usr/share/texmf-dist/tex/latex/base/article.cls\n", - "Document Class: article 2023/05/17 v1.4n Standard LaTeX document class\n", - "(/usr/share/texmf-dist/tex/latex/base/size10.clo))\n", - "(/usr/share/texmf-dist/tex/generic/luatex85/luatex85.sty)\n", - "(/usr/share/texmf-dist/tex/generic/pdftexcmds/pdftexcmds.sty\n", - "(/usr/share/texmf-dist/tex/generic/infwarerr/infwarerr.sty)\n", - "(/usr/share/texmf-dist/tex/generic/ltxcmds/ltxcmds.sty))\n", - "(/usr/share/texmf-site/tex/latex/pgf/frontendlayer/tikz.sty\n", - "(/usr/share/texmf-site/tex/latex/pgf/basiclayer/pgf.sty\n", - "(/usr/share/texmf-site/tex/latex/pgf/utilities/pgfrcs.sty\n", - "(/usr/share/texmf-site/tex/generic/pgf/utilities/pgfutil-common.tex)\n", - "(/usr/share/texmf-site/tex/generic/pgf/utilities/pgfutil-latex.def)\n", - "(/usr/share/texmf-site/tex/generic/pgf/utilities/pgfrcs.code.tex\n", - "(/usr/share/texmf-site/tex/generic/pgf/pgf.revision.tex)))\n", - "(/usr/share/texmf-site/tex/latex/pgf/basiclayer/pgfcore.sty\n", - "(/usr/share/texmf-dist/tex/latex/graphics/graphicx.sty\n", - "(/usr/share/texmf-dist/tex/latex/graphics/graphics.sty\n", - "(/usr/share/texmf-dist/tex/latex/graphics/trig.sty)\n", - "(/usr/share/texmf-dist/tex/latex/graphics-cfg/graphics.cfg)\n", - "(/usr/share/texmf-dist/tex/latex/graphics-def/luatex.def)))\n", - "(/usr/share/texmf-site/tex/latex/pgf/systemlayer/pgfsys.sty\n", - "(/usr/share/texmf-site/tex/generic/pgf/systemlayer/pgfsys.code.tex\n", - "(/usr/share/texmf-site/tex/generic/pgf/utilities/pgfkeys.code.tex\n", - "(/usr/share/texmf-site/tex/generic/pgf/utilities/pgfkeyslibraryfiltered.code.te\n", - "x)) (/usr/share/texmf-site/tex/generic/pgf/systemlayer/pgf.cfg)\n", - "(/usr/share/texmf-site/tex/generic/pgf/systemlayer/pgfsys-luatex.def\n", - "(/usr/share/texmf-site/tex/generic/pgf/systemlayer/pgfsys-common-pdf.def)))\n", - "(/usr/share/texmf-site/tex/generic/pgf/systemlayer/pgfsyssoftpath.code.tex)\n", - "(/usr/share/texmf-site/tex/generic/pgf/systemlayer/pgfsysprotocol.code.tex))\n", - "(/usr/share/texmf-dist/tex/latex/xcolor/xcolor.sty\n", - "(/usr/share/texmf-dist/tex/latex/graphics-cfg/color.cfg)\n", - "(/usr/share/texmf-dist/tex/latex/graphics/mathcolor.ltx))\n", - "(/usr/share/texmf-site/tex/generic/pgf/basiclayer/pgfcore.code.tex\n", - "(/usr/share/texmf-site/tex/generic/pgf/math/pgfmath.code.tex\n", - "(/usr/share/texmf-site/tex/generic/pgf/math/pgfmathutil.code.tex)\n", - "(/usr/share/texmf-site/tex/generic/pgf/math/pgfmathparser.code.tex)\n", - "(/usr/share/texmf-site/tex/generic/pgf/math/pgfmathfunctions.code.tex)\n", - "(/usr/share/texmf-site/tex/generic/pgf/math/pgfmathfunctions.basic.code.tex)\n", - "(/usr/share/texmf-site/tex/generic/pgf/math/pgfmathfunctions.trigonometric.code\n", - ".tex)\n", - "(/usr/share/texmf-site/tex/generic/pgf/math/pgfmathfunctions.random.code.tex)\n", - "(/usr/share/texmf-site/tex/generic/pgf/math/pgfmathfunctions.comparison.code.te\n", - "x) (/usr/share/texmf-site/tex/generic/pgf/math/pgfmathfunctions.base.code.tex)\n", - "(/usr/share/texmf-site/tex/generic/pgf/math/pgfmathfunctions.round.code.tex)\n", - "(/usr/share/texmf-site/tex/generic/pgf/math/pgfmathfunctions.misc.code.tex)\n", - "(/usr/share/texmf-site/tex/generic/pgf/math/pgfmathfunctions.integerarithmetics\n", - ".code.tex) (/usr/share/texmf-site/tex/generic/pgf/math/pgfmathcalc.code.tex)\n", - "(/usr/share/texmf-site/tex/generic/pgf/math/pgfmathfloat.code.tex))\n", - "(/usr/share/texmf-site/tex/generic/pgf/math/pgfint.code.tex)\n", - "(/usr/share/texmf-site/tex/generic/pgf/basiclayer/pgfcorepoints.code.tex)\n", - "(/usr/share/texmf-site/tex/generic/pgf/basiclayer/pgfcorepathconstruct.code.tex\n", - ") (/usr/share/texmf-site/tex/generic/pgf/basiclayer/pgfcorepathusage.code.tex)\n", - "(/usr/share/texmf-site/tex/generic/pgf/basiclayer/pgfcorescopes.code.tex)\n", - "(/usr/share/texmf-site/tex/generic/pgf/basiclayer/pgfcoregraphicstate.code.tex)\n", - "\n", - "(/usr/share/texmf-site/tex/generic/pgf/basiclayer/pgfcoretransformations.code.t\n", - "ex) (/usr/share/texmf-site/tex/generic/pgf/basiclayer/pgfcorequick.code.tex)\n", - "(/usr/share/texmf-site/tex/generic/pgf/basiclayer/pgfcoreobjects.code.tex)\n", - "(/usr/share/texmf-site/tex/generic/pgf/basiclayer/pgfcorepathprocessing.code.te\n", - "x) (/usr/share/texmf-site/tex/generic/pgf/basiclayer/pgfcorearrows.code.tex)\n", - "(/usr/share/texmf-site/tex/generic/pgf/basiclayer/pgfcoreshade.code.tex)\n", - "(/usr/share/texmf-site/tex/generic/pgf/basiclayer/pgfcoreimage.code.tex)\n", - "(/usr/share/texmf-site/tex/generic/pgf/basiclayer/pgfcoreexternal.code.tex)\n", - "(/usr/share/texmf-site/tex/generic/pgf/basiclayer/pgfcorelayers.code.tex)\n", - "(/usr/share/texmf-site/tex/generic/pgf/basiclayer/pgfcoretransparency.code.tex)\n", - " (/usr/share/texmf-site/tex/generic/pgf/basiclayer/pgfcorepatterns.code.tex)\n", - "(/usr/share/texmf-site/tex/generic/pgf/basiclayer/pgfcorerdf.code.tex)))\n", - "(/usr/share/texmf-site/tex/generic/pgf/modules/pgfmoduleshapes.code.tex)\n", - "(/usr/share/texmf-site/tex/generic/pgf/modules/pgfmoduleplot.code.tex)\n", - "(/usr/share/texmf-site/tex/latex/pgf/compatibility/pgfcomp-version-0-65.sty)\n", - "(/usr/share/texmf-site/tex/latex/pgf/compatibility/pgfcomp-version-1-18.sty))\n", - "(/usr/share/texmf-site/tex/latex/pgf/utilities/pgffor.sty\n", - "(/usr/share/texmf-site/tex/latex/pgf/utilities/pgfkeys.sty\n", - "(/usr/share/texmf-site/tex/generic/pgf/utilities/pgfkeys.code.tex))\n", - "(/usr/share/texmf-site/tex/latex/pgf/math/pgfmath.sty\n", - "(/usr/share/texmf-site/tex/generic/pgf/math/pgfmath.code.tex))\n", - "(/usr/share/texmf-site/tex/generic/pgf/utilities/pgffor.code.tex))\n", - "(/usr/share/texmf-site/tex/generic/pgf/frontendlayer/tikz/tikz.code.tex\n", - "(/usr/share/texmf-site/tex/generic/pgf/libraries/pgflibraryplothandlers.code.te\n", - "x) (/usr/share/texmf-site/tex/generic/pgf/modules/pgfmodulematrix.code.tex)\n", - "(/usr/share/texmf-site/tex/generic/pgf/frontendlayer/tikz/libraries/tikzlibrary\n", - "topaths.code.tex)))) (/usr/share/texmf-dist/tex/latex/base/fontenc.sty\n", - "(/usr/share/texmf-dist/tex/latex/lm/t1lmr.fd))\n", - "(/usr/share/texmf-dist/tex/latex/base/inputenc.sty\n", - "\n", - "Package inputenc Warning: inputenc package ignored with utf8 based engines.\n", - "\n", - ") (/usr/share/texmf-dist/tex/latex/lm/lmodern.sty)\n", - "(/usr/share/texmf-dist/tex/latex/base/textcomp.sty)\n", - "(/usr/share/texmf-dist/tex/latex/lastpage/lastpage.sty\n", - "(/usr/share/texmf-dist/tex/latex/lastpage/lastpage2e.sty\n", - "(/usr/share/texmf-dist/tex/latex/lastpage/lastpagemodern.sty) \n", - ")) (/usr/share/texmf-dist/tex/latex/tikz-feynman/tikz-feynman.sty\n", - "(/usr/share/texmf-dist/tex/latex/pgfopts/pgfopts.sty)\n", - "(/usr/share/texmf-dist/tex/latex/tikz-feynman/tikzlibraryfeynman.code.tex\n", - "(/usr/share/texmf-site/tex/generic/pgf/libraries/pgflibraryarrows.meta.code.tex\n", - ")\n", - "(/usr/share/texmf-site/tex/generic/pgf/frontendlayer/tikz/libraries/tikzlibrary\n", - "calc.code.tex)\n", - "(/usr/share/texmf-site/tex/generic/pgf/frontendlayer/tikz/libraries/tikzlibrary\n", - "decorations.code.tex\n", - "(/usr/share/texmf-site/tex/generic/pgf/modules/pgfmoduledecorations.code.tex))\n", - "(/usr/share/texmf-site/tex/generic/pgf/frontendlayer/tikz/libraries/tikzlibrary\n", - "decorations.markings.code.tex\n", - "(/usr/share/texmf-site/tex/generic/pgf/libraries/decorations/pgflibrarydecorati\n", - "ons.markings.code.tex))\n", - "(/usr/share/texmf-site/tex/generic/pgf/frontendlayer/tikz/libraries/tikzlibrary\n", - "decorations.pathmorphing.code.tex\n", - "(/usr/share/texmf-site/tex/generic/pgf/libraries/decorations/pgflibrarydecorati\n", - "ons.pathmorphing.code.tex))\n", - "(/usr/share/texmf-site/tex/generic/pgf/frontendlayer/tikz/libraries/tikzlibrary\n", - "decorations.pathreplacing.code.tex\n", - "(/usr/share/texmf-site/tex/generic/pgf/libraries/decorations/pgflibrarydecorati\n", - "ons.pathreplacing.code.tex))\n", - "(/usr/share/texmf-site/tex/generic/pgf/frontendlayer/tikz/libraries/graphs/tikz\n", - "librarygraphs.code.tex)\n", - "(/usr/share/texmf-site/tex/generic/pgf/frontendlayer/tikz/libraries/tikzlibrary\n", - "patterns.code.tex\n", - "(/usr/share/texmf-site/tex/generic/pgf/libraries/pgflibrarypatterns.code.tex))\n", - "(/usr/share/texmf-site/tex/generic/pgf/frontendlayer/tikz/libraries/tikzlibrary\n", - "positioning.code.tex)\n", - "(/usr/share/texmf-site/tex/generic/pgf/frontendlayer/tikz/libraries/tikzlibrary\n", - "shapes.geometric.code.tex\n", - "(/usr/share/texmf-site/tex/generic/pgf/libraries/shapes/pgflibraryshapes.geomet\n", - "ric.code.tex))\n", - "(/usr/share/texmf-site/tex/generic/pgf/graphdrawing/tex/tikzlibrarygraphdrawing\n", - ".code.tex\n", - "(/usr/share/texmf-site/tex/generic/pgf/graphdrawing/tex/pgflibrarygraphdrawing.\n", - "code.tex))\n", - "(/usr/share/texmf-dist/tex/latex/tikz-feynman/tikzfeynman.keys.code.tex)))\n", - "(/usr/share/texmf-dist/tex/latex/l3backend/l3backend-luatex.def)\n", - "No file hvq.aux.\n", - "(/usr/share/texmf-dist/tex/latex/base/ts1cmr.fd)\n", - "(/usr/share/texmf-dist/tex/context/base/mkii/supp-pdf.mkii\n", - "[Loading MPS to PDF converter (version 2006.09.02).]\n", - ") (/usr/share/texmf-dist/tex/latex/epstopdf-pkg/epstopdf-base.sty\n", - "(/usr/share/texmf-dist/tex/latex/latexconfig/epstopdf-sys.cfg))\n", - "! Undefined control sequence.\n", - "\\pgfmathsetlength ...\\pgfmath@onquick #2\\pgfmath@ \n", - " {\\begingroup \\pgfmath@sele...\n", - "\n", - "l.20 ...ex (Vertex11) [blob,scale=] at (2.5,2.5) {};\n", - " \n", - "! Argument of \\pgfmath@@onquick has an extra }.\n", - " \n", - "\\par \n", - "l.20 ...ex (Vertex11) [blob,scale=] at (2.5,2.5) {};\n", - " \n", - "Runaway argument?\n", - "{\\begingroup \\pgfmath@selectfont \\pgf@x \\unskip \\expandafter \\endgroup \\ETC.\n", - "! Paragraph ended before \\pgfmath@@onquick was complete.\n", - " \n", - "\\par \n", - "l.20 ...ex (Vertex11) [blob,scale=] at (2.5,2.5) {};\n", - " \n", - "! Undefined control sequence.\n", - "\\pgfmathsetlength ...\\pgfmath@onquick #2\\pgfmath@ \n", - " {\\begingroup \\pgfmath@sele...\n", - "\n", - "l.22 ...x (Vertex12) [blob,scale=] at (2.5,-2.5) {};\n", - " \n", - "! Argument of \\pgfmath@@onquick has an extra }.\n", - " \n", - "\\par \n", - "l.22 ...x (Vertex12) [blob,scale=] at (2.5,-2.5) {};\n", - " \n", - "Runaway argument?\n", - "{\\begingroup \\pgfmath@selectfont \\pgf@x \\unskip \\expandafter \\endgroup \\ETC.\n", - "! Paragraph ended before \\pgfmath@@onquick was complete.\n", - " \n", - "\\par \n", - "l.22 ...x (Vertex12) [blob,scale=] at (2.5,-2.5) {};\n", - " \n", - "(/usr/share/texmf-dist/tex/latex/lm/ot1lmr.fd)\n", - "(/usr/share/texmf-dist/tex/latex/lm/omllmm.fd)\n", - "(/usr/share/texmf-dist/tex/latex/lm/omslmsy.fd)\n", - "(/usr/share/texmf-dist/tex/latex/lm/omxlmex.fd) [1{/home/apn/.texlive/texmf-var\n", - "/fonts/map/pdftex/updmap/pdftex.map}]\n", - "\n", - "enddocument/afterlastpage: lastpage setting LastPage.\n", - "(./hvq.aux)\n", - "\n", - "LaTeX Warning: Label(s) may have changed. Rerun to get cross-references right.\n", - "\n", - "\n", - "Package lastpage Warning: Rerun to get the references right on input line 108.\n", - "\n", - ")\n", - "(see the transcript file for additional information)\n", - " 413 words of node memory still in use:\n", - " 3 hlist, 1 vlist, 1 rule, 2 glue, 3 kern, 1 glyph, 5 attribute, 48 glue_spec\n", - ", 5 attribute_list, 2 write nodes\n", - " avail lists: 1:8,2:4866,3:1320,4:62,5:41,6:2,7:2226,8:3,9:747\n", - "{/usr/share/texmf-dist/fonts/enc/dvips/lm/lm-mathit.enc}{/usr/share/texmf-dist/\n", - "fonts/enc/dvips/lm/lm-rm.enc}{/usr/share/texmf-dist/fonts/enc/dvips/lm/lm-maths\n", - "y.enc}\n", - "Output written on hvq.pdf (1 page, 27694 bytes).\n", - "Transcript written on hvq.log.\n", - "\n" + "/home/apn/.local/lib/python3.13/site-packages/feynml/pdgid.py:111: UserWarning: Type remnant is not consistent with pdgid 2, which is fermion. Using remnant now.\n", + " warnings.warn(\n", + "/home/apn/.local/lib/python3.13/site-packages/feynml/pdgid.py:111: UserWarning: Type remnant is not consistent with pdgid -2, which is anti fermion. Using remnant now.\n", + " warnings.warn(\n" ] }, { - "ename": "CalledProcessError", - "evalue": "Command '['lualatex', '-shell-escape', '--interaction=nonstopmode', '/tmp/tmpgaexcn32/hvq.tex']' returned non-zero exit status 1.", - "output_type": "error", - "traceback": [ - "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m", - "\u001b[0;31mCalledProcessError\u001b[0m Traceback (most recent call last)", - "Cell \u001b[0;32mIn[2], line 31\u001b[0m\n\u001b[1;32m 1\u001b[0m fd \u001b[38;5;241m=\u001b[39m FeynmanDiagram()\u001b[38;5;241m.\u001b[39madd(\n\u001b[1;32m 2\u001b[0m v1 \u001b[38;5;241m:=\u001b[39m Vertex()\u001b[38;5;241m.\u001b[39mwith_xy(\u001b[38;5;241m7\u001b[39m,\u001b[38;5;241m0\u001b[39m),\n\u001b[1;32m 3\u001b[0m v2 \u001b[38;5;241m:=\u001b[39m Vertex()\u001b[38;5;241m.\u001b[39mwith_xy(\u001b[38;5;241m10\u001b[39m,\u001b[38;5;241m0\u001b[39m),\n\u001b[0;32m (...)\u001b[0m\n\u001b[1;32m 28\u001b[0m Leg(pdgid\u001b[38;5;241m=\u001b[39m\u001b[38;5;241m-\u001b[39m\u001b[38;5;241m2\u001b[39m,\u001b[38;5;28mtype\u001b[39m\u001b[38;5;241m=\u001b[39m\u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mremnant\u001b[39m\u001b[38;5;124m\"\u001b[39m,label\u001b[38;5;241m=\u001b[39m\u001b[38;5;124m\"\u001b[39m\u001b[38;5;124m$X$\u001b[39m\u001b[38;5;124m\"\u001b[39m)\u001b[38;5;241m.\u001b[39mwith_target(v4)\u001b[38;5;241m.\u001b[39mwith_outgoing()\u001b[38;5;241m.\u001b[39mwith_xy(\u001b[38;5;241m10\u001b[39m,\u001b[38;5;241m-\u001b[39m\u001b[38;5;241m7\u001b[39m),\n\u001b[1;32m 29\u001b[0m )\n\u001b[1;32m 30\u001b[0m scale_positions(fd,\u001b[38;5;241m0.5\u001b[39m)\n\u001b[0;32m---> 31\u001b[0m \u001b[43mfd\u001b[49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43mrender\u001b[49m\u001b[43m(\u001b[49m\u001b[43mauto_position\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[38;5;28;43;01mFalse\u001b[39;49;00m\u001b[43m,\u001b[49m\u001b[43mrender\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[38;5;124;43m\"\u001b[39;49m\u001b[38;5;124;43mtikz\u001b[39;49m\u001b[38;5;124;43m\"\u001b[39;49m\u001b[43m,\u001b[49m\u001b[43mfile\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[38;5;124;43m\"\u001b[39;49m\u001b[38;5;124;43mhvq\u001b[39;49m\u001b[38;5;124;43m\"\u001b[39;49m\u001b[43m)\u001b[49m\n", - "File \u001b[0;32m~/.local/lib/python3.13/site-packages/feynml/feynmandiagram.py:413\u001b[0m, in \u001b[0;36mFeynmanDiagram.render\u001b[0;34m(self, render, show, file, auto_position, auto_position_legs, deepcopy)\u001b[0m\n\u001b[1;32m 411\u001b[0m fd \u001b[38;5;241m=\u001b[39m auto_bend(fd)\n\u001b[1;32m 412\u001b[0m renderer \u001b[38;5;241m=\u001b[39m renderall\u001b[38;5;241m.\u001b[39mrenderer_from_string(render)\n\u001b[0;32m--> 413\u001b[0m \u001b[43mrenderer\u001b[49m\u001b[43m(\u001b[49m\u001b[43mfd\u001b[49m\u001b[43m)\u001b[49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43mrender\u001b[49m\u001b[43m(\u001b[49m\u001b[43mshow\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mshow\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[43mfile\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mfile\u001b[49m\u001b[43m)\u001b[49m\n", - "File \u001b[0;32m~/.local/lib/python3.13/site-packages/pyfeyn2/render/latex/latex.py:63\u001b[0m, in \u001b[0;36mLatexRender.render\u001b[0;34m(self, file, show, resolution, width, height, clean_up, temp_dir)\u001b[0m\n\u001b[1;32m 61\u001b[0m tfile \u001b[38;5;241m=\u001b[39m re\u001b[38;5;241m.\u001b[39msub(\u001b[38;5;124mr\u001b[39m\u001b[38;5;124m\"\u001b[39m\u001b[38;5;124m\\\u001b[39m\u001b[38;5;124m.pdf$\u001b[39m\u001b[38;5;124m\"\u001b[39m, \u001b[38;5;124m\"\u001b[39m\u001b[38;5;124m\"\u001b[39m, os\u001b[38;5;241m.\u001b[39mpath\u001b[38;5;241m.\u001b[39mbasename(file)\u001b[38;5;241m.\u001b[39mstrip())\n\u001b[1;32m 62\u001b[0m tfile \u001b[38;5;241m=\u001b[39m os\u001b[38;5;241m.\u001b[39mpath\u001b[38;5;241m.\u001b[39mjoin(temp_dir\u001b[38;5;241m.\u001b[39mname, tfile)\n\u001b[0;32m---> 63\u001b[0m \u001b[38;5;28;43mself\u001b[39;49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43mgenerate_pdf\u001b[49m\u001b[43m(\u001b[49m\n\u001b[1;32m 64\u001b[0m \u001b[43m \u001b[49m\u001b[43mtfile\u001b[49m\u001b[43m,\u001b[49m\n\u001b[1;32m 65\u001b[0m \u001b[43m \u001b[49m\u001b[43mclean_tex\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mclean_up\u001b[49m\u001b[43m,\u001b[49m\n\u001b[1;32m 66\u001b[0m \u001b[43m \u001b[49m\u001b[43mcompiler\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[38;5;124;43m\"\u001b[39;49m\u001b[38;5;124;43mlualatex\u001b[39;49m\u001b[38;5;124;43m\"\u001b[39;49m\u001b[43m,\u001b[49m\n\u001b[1;32m 67\u001b[0m \u001b[43m \u001b[49m\u001b[43mcompiler_args\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43m[\u001b[49m\u001b[38;5;124;43m\"\u001b[39;49m\u001b[38;5;124;43m-shell-escape\u001b[39;49m\u001b[38;5;124;43m\"\u001b[39;49m\u001b[43m]\u001b[49m\u001b[43m,\u001b[49m\n\u001b[1;32m 68\u001b[0m \u001b[43m\u001b[49m\u001b[43m)\u001b[49m\n\u001b[1;32m 69\u001b[0m file \u001b[38;5;241m+\u001b[39m\u001b[38;5;241m=\u001b[39m \u001b[38;5;124m\"\u001b[39m\u001b[38;5;124m.pdf\u001b[39m\u001b[38;5;124m\"\u001b[39m\n\u001b[1;32m 70\u001b[0m tfile \u001b[38;5;241m+\u001b[39m\u001b[38;5;241m=\u001b[39m \u001b[38;5;124m\"\u001b[39m\u001b[38;5;124m.pdf\u001b[39m\u001b[38;5;124m\"\u001b[39m\n", - "File \u001b[0;32m~/.local/lib/python3.13/site-packages/pylatex/document.py:276\u001b[0m, in \u001b[0;36mDocument.generate_pdf\u001b[0;34m(self, filepath, clean, clean_tex, compiler, compiler_args, silent)\u001b[0m\n\u001b[1;32m 273\u001b[0m command \u001b[38;5;241m=\u001b[39m [compiler] \u001b[38;5;241m+\u001b[39m arguments \u001b[38;5;241m+\u001b[39m compiler_args \u001b[38;5;241m+\u001b[39m main_arguments\n\u001b[1;32m 275\u001b[0m \u001b[38;5;28;01mtry\u001b[39;00m:\n\u001b[0;32m--> 276\u001b[0m output \u001b[38;5;241m=\u001b[39m \u001b[43msubprocess\u001b[49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43mcheck_output\u001b[49m\u001b[43m(\u001b[49m\n\u001b[1;32m 277\u001b[0m \u001b[43m \u001b[49m\u001b[43mcommand\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[43mstderr\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43msubprocess\u001b[49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43mSTDOUT\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[38;5;241;43m*\u001b[39;49m\u001b[38;5;241;43m*\u001b[39;49m\u001b[43mcheck_output_kwargs\u001b[49m\n\u001b[1;32m 278\u001b[0m \u001b[43m \u001b[49m\u001b[43m)\u001b[49m\n\u001b[1;32m 279\u001b[0m \u001b[38;5;28;01mexcept\u001b[39;00m (\u001b[38;5;167;01mOSError\u001b[39;00m, \u001b[38;5;167;01mIOError\u001b[39;00m) \u001b[38;5;28;01mas\u001b[39;00m e:\n\u001b[1;32m 280\u001b[0m \u001b[38;5;66;03m# Use FileNotFoundError when python 2 is dropped\u001b[39;00m\n\u001b[1;32m 281\u001b[0m os_error \u001b[38;5;241m=\u001b[39m e\n", - "File \u001b[0;32m/usr/lib/python3.13/subprocess.py:472\u001b[0m, in \u001b[0;36mcheck_output\u001b[0;34m(timeout, *popenargs, **kwargs)\u001b[0m\n\u001b[1;32m 469\u001b[0m empty \u001b[38;5;241m=\u001b[39m \u001b[38;5;124mb\u001b[39m\u001b[38;5;124m'\u001b[39m\u001b[38;5;124m'\u001b[39m\n\u001b[1;32m 470\u001b[0m kwargs[\u001b[38;5;124m'\u001b[39m\u001b[38;5;124minput\u001b[39m\u001b[38;5;124m'\u001b[39m] \u001b[38;5;241m=\u001b[39m empty\n\u001b[0;32m--> 472\u001b[0m \u001b[38;5;28;01mreturn\u001b[39;00m \u001b[43mrun\u001b[49m\u001b[43m(\u001b[49m\u001b[38;5;241;43m*\u001b[39;49m\u001b[43mpopenargs\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[43mstdout\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mPIPE\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[43mtimeout\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mtimeout\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[43mcheck\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[38;5;28;43;01mTrue\u001b[39;49;00m\u001b[43m,\u001b[49m\n\u001b[1;32m 473\u001b[0m \u001b[43m \u001b[49m\u001b[38;5;241;43m*\u001b[39;49m\u001b[38;5;241;43m*\u001b[39;49m\u001b[43mkwargs\u001b[49m\u001b[43m)\u001b[49m\u001b[38;5;241m.\u001b[39mstdout\n", - "File \u001b[0;32m/usr/lib/python3.13/subprocess.py:577\u001b[0m, in \u001b[0;36mrun\u001b[0;34m(input, capture_output, timeout, check, *popenargs, **kwargs)\u001b[0m\n\u001b[1;32m 575\u001b[0m retcode \u001b[38;5;241m=\u001b[39m process\u001b[38;5;241m.\u001b[39mpoll()\n\u001b[1;32m 576\u001b[0m \u001b[38;5;28;01mif\u001b[39;00m check \u001b[38;5;129;01mand\u001b[39;00m retcode:\n\u001b[0;32m--> 577\u001b[0m \u001b[38;5;28;01mraise\u001b[39;00m CalledProcessError(retcode, process\u001b[38;5;241m.\u001b[39margs,\n\u001b[1;32m 578\u001b[0m output\u001b[38;5;241m=\u001b[39mstdout, stderr\u001b[38;5;241m=\u001b[39mstderr)\n\u001b[1;32m 579\u001b[0m \u001b[38;5;28;01mreturn\u001b[39;00m CompletedProcess(process\u001b[38;5;241m.\u001b[39margs, retcode, stdout, stderr)\n", - "\u001b[0;31mCalledProcessError\u001b[0m: Command '['lualatex', '-shell-escape', '--interaction=nonstopmode', '/tmp/tmpgaexcn32/hvq.tex']' returned non-zero exit status 1." - ] + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAANwAAAEZCAMAAADL1rsAAAAAIGNIUk0AAHomAACAhAAA+gAAAIDoAAB1MAAA6mAAADqYAAAXcJy6UTwAAAHXUExURf///wAAAAAAAAAAAEtLSwAAABsbG4qKiktLS2pqaq+vr46OjgAAAAAAAAAAAAAAAD09PQAAAAAAAAAAAGpqagAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAIyMjAAAAAAAAAAAAAAAAAAAAK6urpKSkv8AAHd3d/8AAP8AAP8AAP8AAP8AAGdnZ/8AAP8AAP8AAP8AAP8AAP8AAAAAAAAAAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AALEAIuwAAwwAwgAA/wAA/wAA/wAA/wAA/wAA/wAA/wAA/wAA/wAA/wAA/wAA/wAA/wAA/wAA/wAA/wAA/wAA/0dHR1VVVRMTEx4eHkFBQTo6OrKysnl5ec/Pz4ODg01NTY2NjWdnZ3d3d93d3f///5mZmaqqqiIiImZmZlVVVYiIiO7u7ru7u8zMzAAAAERERDIyMhkZGXBwcK+vryAgICgoKF9fXxERESQkJBsbG1JSUv8AAAAA/zMzMwwMDFdXV6KioggICDg4OGdnZwg2JwgAAAB7dFJOUwAAEXetRKA/5tr+mSHHdqPW9dbj9TP6tmb+XI4/7feIVd/vmMfOEdbasMyPu6hEwPHCmaqqmUFwZtYz3fVQ7rhViPjgYHci5+uQrYCFXOjk4f181ajH899dpt/WZhHMM1Uiqoi7md3ud0St1pD0es3y69fa9I/+Iexc0pKxfnkAAAABYktHRACIBR1IAAAACXBIWXMAAABkAAAAZAAPlsXdAAAAB3RJTUUH6QEPDSo5oCoOfwAAC89JREFUeNrt3Yt/m8YdAPCcs2TdJuRabqTGStdua68OZjiFtSVikXVtkOuse3dbt1mJHS997P1+OUrTJinttu4p7/3H7g4kWw9Ad6CDH4zf5yPbkg4++ebH/TgOhE6dOoVyH0unPxT4+qn8486cPbzd/3ARcUunH7rT7/ff+kjxcF7S+rcP77790YLhRkl75+69+/3+x+DjKoqiVFFVUZbntTxJ2oN+//6dh8/AxyG0UqM/VqvRjSaT9u4jZ88FNgOHW64j1KhENplO2qNLIQ3B4dD5ZSXKxpk0oDilvhb+JnfSgOKW3TCcSNJg4pprq7XAN8SSBhJXbaLq+dk+J5w0iLi1Jv1Rb0y9GiNp8HDVmkt3cM2VujL2YrykwcMFRNykwcclSBp0XKKkgcYlTRpgXPKkQcUtJGkwcQtKGkDc4pIGDrfIpMHCLThpkHALTxoYnIykpYRzhxH2vpykpYEbqlz/z9n3pSVNEFdprFabyhpn65FtJPMe0zqJSRPELaOVC6i6Kkib0o3x5CZNEFep1KnwMRFcwGMUspMmiEPNGkK1Bm9r5Abr2COFpInilAaq1Cu8rcNkbjpJE8WtNKo1/nrihuo+nkbSBHGVFXSBf71hWyQtKW4aSRPENRXOhj4uTEffSSNpYrhqrVbla+njUIQuhaSJ4R5/QmS1bsQDuSkkbYQ78zhPfML9pPIprpYs3Ig+h9wUkjbCPdnniffYMJHVAa6I2ipnRmEycU894In3/bH9737/B57WbtRWmRIOY4p7+pAnPhgeufzxTzyto/tcOrb1i2Kb5Z/Z2IJrs4zsc6nY1A2N4j59jyf+4rp//RtXSxYZ9zld27z0jGFy7wo+86xI+c64z2nPqc/ruqwjcdfNtM9tqhb7JWuaIWpsKd2GL3k2ebgMt0qsXtal4lBkvZRsszTZuPCjVZnRsgefvaxqSC4uTCc5b+YVv5ZIxgXr3JZc3OaYTeqkbIDObdlSdbi9MfZM6oyzOz0p67JOIVGH1S09LdxkzRx2t+S6jkk6nWCbZaSIOzkPclJJEuswwgFr0O3NF14c1cmUcAEhZ8skV8drSWa4RDrsmMzQwdNvXFyftmV08jFR7ro2tc0sj7c3ZlpmdGY1kU7ttGZt6pYDBZdI55jTtbJL6yTRZxpmdk48gU63p59vXnppYM42zO6Ef7CO7sGMbvSCuDO1ZECdzBoXqCMWwnb0YthBhEy8srMebMv0Uo1ZXbdNC6ERvVSHLYjG9gR4+1pI00yvQ5nROaaXvKjQGYs4JyUFqxsEIm5GZ9F/ptoVWAHW6fEbcULezfgKoikdNnRiCyzetXfalzUDh7yd9eVR07lrzetyE0F21KiNOGvcjG5Ol5uMz21Hts4cN9SRkc0k/H0Or29Fvp89ztMZg474ghF1EgyO6q4NBprgQl39yvrLoXUSDg5R22Agsgdgu/ud7WfC6yQcnMFsAyK0jB5dJ8HgiGcbtLHIQle2OaoqABx2bE8nsg/A6ypHKwA4GrpJcZcFbOrzPAUIBo4WCHJp8HnepjvrG4RnIAMFR+PlL3Adm2N754tbxOHpoYBwnDMPXHUSHo5Pd6XNXXlA4Xh0WOUvPLBw83XsnHBecdE6nexs89VJmLhInX11/RpfnQSKi9DhTe46CRUXrhtdO5NnXJgOqy8JrggiLlhH6+SgCLhZXYtc3N4ipuBqYOJmdOrV9Q1HF10LUNy0TrROwsZN6vALcWxwceM6rL4YaxVwcSc6Wid5JhVyhfN1LXOzvWWIzmrCx3k6+zl1QxeukznAMd0gVp3MAw61vvTl+DboOPSVryZYGDpOJyWuxIGLbsy9QC5wiaLE5TWg4wrd5wpdLUtcXnGtBONm8LhEUeLyGtBxWOyqqXzhCl0tS1yJK1qUuLxGicswCl1QSlyJAxiFniBKFCUurwEdV+gJolG1DL0HSm5xFiGvfI14nwoPvgdKnnGO/xEm0Yv18oHDPo7t6jhua5AznP/JOvbZao7bGuQO12I4h++2BrnDIdX/0CDHbQ1yiLMGA5Yw4dsa5AKHvz5gewDR2xrkA4de+Yb/W+y2BjnBvXpx+EeMLgceNxp+Cd3WIC+4RFHi8holLsMo9OxXiStxAKPQVzMkihKX14CO+7+YlC1xJQ5QlCcfS1zRAjqu0H2u0NWyxJW4okWJy2tAxxV6DqXQBaXE5RX36jeLi+t+69u7hcVhu3c9vg4qDuv+bSy1bu/6jaLhNMP2ZiwNHfVirwQqTsVYZcepttAd1XOCs1Ssm+y+496z3l5xcGzI5djE7hjDm7bFrCoQcarGLmLrOrZxfFVzPB1EnM02SjR59Wjv+n4xcIaD2Jd/sOsQkwVEHFYNx+wSO8nRDgRch3QmLsfDmk3oL8sYaNbMTuCGaL/LFqdrCE18x6FuYoNdWNkNuqGlcFXJFse+VXRyN01rv0bLvxP4hTuhur2bB3sBu8JscYPudI5atmURbIX0t953gtfTQ72gAWi2uDbNERk9waZJc9Y1VNuOceFoUGSLsyxLPe5yFtFt9qQl/vEWulHe2oeGQ+NdrqXSjRLPvZp5b7bfHeyi3msBTTPFkc5oZOw/1TDRUWfOgcBsVdl/nZJvQcM5nc7wqxWw9x0Ljtq2bWPeR8lmdG/c9JIHDHcShBA2nKR7b47jt95UlnYPEAocegLB0ZrpfbGJGec7gXu33jwI6nJQcFh16ME3+3bYWHEjsMtBwGGnxXZyattQVe7E+X2sN/kUIM4kKmG/DZGPNnpVZXe0MR7cPNiHiaPH2yYbSop9hozqdo+O3oxsAwDXslvsoLsleGf73nePjo5uQsZhNn2n26RtmaJfnk7zRmMfME7zTFjX5u66g21HB4BxNm55Rzeq8NzrLR/3OmCcY6COGeuTxKj3xmtMF3VwnnVBMTRLbRmi/W0YezePjr4HE6faJmZfJ2Sb8S/S3v/+D/Zh4uhRAEJxPrk/EREngTLEadaC5l73AeKwShyj6ySfew2d8csGhzV6FEBHzOZAdZLPBYXO+GWD0w1sEhQy9xpHF3z+LqPM2V3MNsc5X3bOrws+bZ5Rn+vYuuNgK9FJ4fmRPg4bbMq1ZdoLm3v1IyB56eMcw/IGJHGGXBERVFXSx7U0f6rE6CRe1TxdBn3OMLCh04O4Ra93tmbKx7mjOH6FtG3TjjWJN0c3PRKTjGMmpvIeJzyuudfkIRc3JhtK/ZCQNj8mr8aRiXMnZMdWFHvulQM3UVVkZ276YamECMy9JtPJwwXkzcudvoChcqTuZG8uEYcCdWNlRZLu5E9pODfkgWTjxkIWzs1Sd2NPNi5UJx83qioScSF9LkUdN+7Z00tCtoitMoVe51/BSHE//BFP/Nh1f/JTrpYs3KiKklpJobgn+zzxHhv7/v3Ofa7G/cj+libuqQc88b4/tP/HP//F09rNts95sbdHcU8f8sQHPu7f//kvT+us+xyLXk/eZgmjz/EXlJ+JFBQYfY6voeCuILLPPSS0rhRwguFG9bn+7cOzZ0DhKo3ValNZ48VFbJXvvEV75Z0U0sePW0YrF1B1lRsXPgJ78Pbdw9v9FNInkLlKnQof411xmM51f97v3793N430CfS5Zg2hWoMb54ZtlY8+zHYoIumrKIpSRVVFWZaGUxqoUq9wrzlkq2S/zp195F2x9K3Q/1i0WhWzieBWGtUabz2J0HmxJJi+ZdolGvz/scK4ygq6ILbugC1zbP8tlr7zy4qwTQDXVIRXHjop64dI+pS6yEYjiqvWaqKb/LFnOKkXMOziTt+yKxMXO2ZPhEwEX/qaa6s1iLj5MT991SaqnpfZ52TGnPStNemPOvcuFhgORaWvWnNpb2+u1IUrGhic+L4vVzgUZ+iSI9yC0wcNhxaZPoC4xaUPJA4tKH1QcQtJH1wcSp4+0Lik6QOOQ4nSBx+XIH15wKG46csJLl76coNDMdKXJ5xw+vKFQ2Lpyx1OJH05xCHu9OUTF5C+cwXCoen0/aJYuIn0/fJXRcOhk/T9OujNvONG6ftNMXGIpe+3ga//D+2WvY6TsMz+AAAAJXRFWHRkYXRlOmNyZWF0ZQAyMDI1LTAxLTE1VDEzOjQyOjU3KzAwOjAwdKYZ0gAAACV0RVh0ZGF0ZTptb2RpZnkAMjAyNS0wMS0xNVQxMzo0Mjo1NyswMDowMAX7oW4AAAAodEVYdGRhdGU6dGltZXN0YW1wADIwMjUtMDEtMTVUMTM6NDI6NTcrMDA6MDBS7oCxAAAAE3RFWHRwZGY6VmVyc2lvbgBQREYtMS41UzZawQAAAEp0RVh0c2lnbmF0dXJlADM2ZDM3YWUxMTQ2M2IxNTg0MmUyNTE3MzdlZjkzMTU5ZjVjZDY0ZWQ4N2M5YzI2NTIxN2E4N2U3NjNhMmQzOTMjk100AAAAAElFTkSuQmCC", + "text/plain": [ + "" + ] + }, + "metadata": {}, + "output_type": "display_data" } ], "source": [ "fd = FeynmanDiagram().add(\n", - " v1 := Vertex().with_xy(7,0),\n", - " v2 := Vertex().with_xy(10,0),\n", " v3 := Vertex().with_xy(5,+5).with_shape('blob'),\n", " v4 := Vertex().with_xy(5,-5).with_shape('blob'),\n", - " Propagator(pdgid=21).connect(v1,v2).with_color('red'),\n", - " Propagator(pdgid=-2).connect(v3,v1).with_color('red'),\n", - " Propagator(pdgid=2).connect(v4,v1).with_color('red'),\n", - " Leg(pdgid=2212,label=\"$p/A$\",target=v3,sense='incoming').with_xy(0,5),\n", - " Leg(pdgid=2212,label=\"$p/A$\",target=v4).with_incoming().with_xy(0,-5).with_style(\"label-side: right\"),\n", - " v5 := Vertex().with_xy(12,2),\n", - " v6 := Vertex().with_xy(15,4),\n", - " Propagator(pdgid=-4).connect(v2,v5).with_color('red'),\n", - " Propagator(pdgid=-24).connect(v5,v6),\n", - " Leg(pdgid=11,sense='outgoing').with_target(v6).with_xy(20,5).with_color('blue'),\n", - " Leg(pdgid=-12,sense='outgoing').with_target(v6).with_xy(20,3),\n", - " Leg(pdgid=-3,sense='outgoing').with_target(v5).with_xy(20,1),\n", + " Leg(pdgid=2212,label=\"$p$\",target=v3,sense='incoming').with_xy(0,5),\n", + " Leg(pdgid=2212,label=\"$p$\",target=v4).with_incoming().with_xy(0,-5).with_style(\"label-side: right\"),\n", "\n", - " v7 := Vertex().with_xy(12,-2),\n", - " v8 := Vertex().with_xy(15,-4),\n", - " Propagator(pdgid=4).connect(v2,v7).with_color('red'),\n", - " Propagator(pdgid=24).connect(v7,v8),\n", - " Leg(pdgid=-11,sense='outgoing').with_target(v8).with_xy(20,-5).with_color('blue'),\n", - " Leg(pdgid=12,sense='outgoing').with_target(v8).with_xy(20,-3),\n", - " Leg(pdgid=3,sense='outgoing').with_target(v7).with_xy(20,-1),\n", - " \n", " Leg(pdgid=2,type=\"remnant\",label=\"$X$\",sense='outgoing').with_target(v3).with_xy(10,+7),\n", " Leg(pdgid=-2,type=\"remnant\",label=\"$X$\").with_target(v4).with_outgoing().with_xy(10,-7),\n", + " \n", + " v1 := Vertex().with_xy(8,2),\n", + " v2 := Vertex().with_xy(8,-2),\n", + " Propagator(type=\"squark\", label=r\"$\\tilde q$\").connect(v1,v2).with_color('red'),\n", + " Propagator(pdgid=1,label=r\"$q$\").connect(v3,v1).with_color('red'),\n", + " Propagator(pdgid=21).connect(v4,v2).with_color('red'),\n", + " \n", + " v5 := Vertex().with_xy(11,5),\n", + " Propagator(type=\"gaugino\", label=r\"$\\tilde \\chi$\").connect(v1,v5).with_color('red'),\n", + " #Propagator(pdgid=-24).connect(v5,v6),\n", + " #Leg(pdgid=11,sense='outgoing').with_target(v6).with_xy(20,5).with_color('blue'),\n", + " #Leg(pdgid=-12,sense='outgoing').with_target(v6).with_xy(20,3),\n", + " #Leg(pdgid=-3,sense='outgoing').with_target(v5).with_xy(20,1),\n", + "\n", + " v7 := Vertex().with_xy(11,-5),\n", + " Propagator(type=\"squark\", label=r\"$\\tilde q$\").connect(v2,v7).with_color('blue'),\n", + " #Propagator(pdgid=24).connect(v7,v8),\n", + " #Leg(pdgid=-11,sense='outgoing').with_target(v8).with_xy(20,-5).with_color('blue'),\n", + " #Leg(pdgid=12,sense='outgoing').with_target(v8).with_xy(20,-3),\n", + " #Leg(pdgid=3,sense='outgoing').with_target(v7).with_xy(20,-1),\n", + " \n", ")\n", "scale_positions(fd,0.5)\n", - "fd.render(auto_position=False,render=\"tikz\",file=\"hvq\")" + "fd.render(auto_position=False,render=\"tikz\",file=\"sus\")" ] }, { "cell_type": "code", - "execution_count": 24, + "execution_count": 3, "id": "4aa7709f-b9aa-4723-a6fa-faf25de62455", "metadata": {}, "outputs": [ @@ -345,15 +179,15 @@ "name": "stderr", "output_type": "stream", "text": [ - "/home/apn/.cache/pypoetry/virtualenvs/pyfeyn2-rvdu2REa-py3.11/lib/python3.11/site-packages/feynml/pdgid.py:111: UserWarning: Type remnant is not consistent with pdgid 2, which is fermion. Using remnant now.\n", + "/home/apn/.local/lib/python3.13/site-packages/feynml/pdgid.py:111: UserWarning: Type remnant is not consistent with pdgid 2, which is fermion. Using remnant now.\n", " warnings.warn(\n", - "/home/apn/.cache/pypoetry/virtualenvs/pyfeyn2-rvdu2REa-py3.11/lib/python3.11/site-packages/feynml/pdgid.py:111: UserWarning: Type remnant is not consistent with pdgid -2, which is anti fermion. Using remnant now.\n", + "/home/apn/.local/lib/python3.13/site-packages/feynml/pdgid.py:111: UserWarning: Type remnant is not consistent with pdgid -2, which is anti fermion. Using remnant now.\n", " warnings.warn(\n" ] }, { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAASoAAAEZCAMAAAAXCpgeAAAAIGNIUk0AAHomAACAhAAA+gAAAIDoAAB1MAAA6mAAADqYAAAXcJy6UTwAAAG8UExURf///wAAAAAAAAAAAEtLSwAAABsbG4qKiktLS2pqaq+vr46OjgAAAAAAAAAAAAAAAD09PQAAAAAAAAAAAGpqagAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAIyMjAAAAAAAAAAAAAAAAAAAAK6urpKSknd3d2dnZwAAAAAAAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAAAA/wAA/wAA/wAA/wAA/wAA/wAA/wAA/wAA//8AAAAA/wAA/wAA/wAA//8AAAAA/wAA/wAA/wAA/wAA/wAA//8AAP8AABwAagAA/9IAHgAA/wAA/wAA//8AAP8AAP8AAP8AAP8AAEdHR1VVVRMTEx4eHkFBQTo6OrKysnl5ec/Pz4ODg01NTY2NjWdnZ3d3d93d3f///5mZmaqqqiIiImZmZlVVVYiIiO7u7ru7u8zMzAAAAERERDIyMhkZGXBwcK+vryAgICgoKF9fXxERESQkJBsbG1JSUv8AAAAA/zMzMwwMDFdXV6KioggICDg4OGdnZ0Bn/HoAAABydFJOUwAAEXetRKA/5tr+mSHHdqPW9dbj9TP6tmb+XI4/7feIVd/vmMfO1qiqmapEVZndEYh3uzMiZszu6+i00rCY2Y/45N8iVardiJlEdxFwZjPuu8ZdkP2S8+Yx4enPYKCt8VzbuO/3es3y69fa9I/+Iexc0nG6YLsAAAABYktHRACIBR1IAAAACXBIWXMAAABkAAAAZAAPlsXdAAAAB3RJTUUH6AEeDRUYl2TXQgAAERNJREFUeNrtnYl727YZxoN0ybpNomq5lRrLWXe03cr71u52cxvHzuJ02X1v3X1vidK0Tcp2W3equ//h4QMpW5JFEpQI4IOb93ms6JFpkfrpxQvgI8icO3eOaK/zD71Lwl7O6Y/qwsVbtyfvFr8f3VGdf+jhO5PJ5OX3iN+V3qiYoSa3b9195b3id6YxqpmhXr372r3J5H3C94cMVafb7RrE6HZ7dVueGOr+ZHLvziMXhB8bMlSEbPXpw7ZRvdGioV5/9OJjEo4MHaregJBhp3KTZUM9fl7KkaFDRS71ulWklBiKCR+q7mCn/JeKDMWED1UvK0OlzlBM6FCNdrb7K3+h0lBM2FAZI2JcOp1Vig3FhAzVzog+DIZLryo3FBMqVEY/owOq0dagO/ciBkMxoUK1QjgMxYQaFRpDMSFGhchQTFhR4TIUE05U2AzFhBAVQkMxoUOF0lBMuFBhNRRTM1S7Qo8Fr6GYmqGajaIvt38gqA3F1AhVp5jzD7O2DwO5oZgaoerllaTOMOvw/1G98BuKqRGqmalIZvD/UVao7Pc6GIoTVWe4bYy64KfdvDhiGOQSL6qCUZY/Pf17TQzFiapHtnaJsU2fDfMOkAK7XHuabkZqxon9LLPSxlCcqDqdAeUFfV7e/kY9w9jq1r5xAWqJ1RwsnQzFiYqMKKI+tdIOs1IHWmGXD9WKn5n0MhQvqu6QdAa0y8tPz7G86nIMrLLVrOBHO0PxotoaGn1IdXCS0R8Y0BwH9bYq4ZTpaChOVJ2tfD4zGjV646yU1fv1MxQnqlFhIK54miNV8kMR6mcoPlRGv88GUZ2GqMpY0d/oZyg+VDPtNqsqlLY/YKWdoZi4UT3xgSZvm1X8kEw7QzFRVBee4NEHsw91P8y1JSiryCqSaWeoGaonJzx6A6ZzkMZcqmp/p+Y3moiieuo+j97M6wO//8MfebbOqtqfvqievsWjt4payp/+zLN1dVap/tBro2rQAP8C42yuBliZVao/9NqoPvIaj/6aZX/7O9eWoDOaVXwbfvSZJh38Gc0qEW+bZWcyq4S8b9UcUPVnXlPCUJ259ifuRHxlH6ilxKEqr+1pKnHLO7Kz5SmhK2GyM+UpsYuGsrPkKcHrq7LlU6Y6kxK9FG2eld6gJKzay45lWqo/7GaSuMDRdlR/2M0kcy2o66n+tBtJJio/UP1p11EQRvkTmahi9/ip5edCHV9JmNBH1yoOW+pi7JNgj61csWocq2WlPn10Y5ceX0gcdtiJVFQaBHvsmh58p6ZNSMoiwyW59125S/xPgt0LnIjEYWvvbPu2n2z0Dl6aUg/5CeAJ6ftRVB48mgkLq8SRi+ok2BOwdeK39cYWhQ5GWE/hmPom8W2THqEF359jWfDccRKGEHYwjuWiij82e+Z5LilCoA2ZyXpjEX9M/zAyI3owvk/GtL8LqKWIbbJv0aQvkChNU99PI9nX2Hz8E7NncCzuJm+1oGlM4rTRXwRjCE7Xo4Atn4QeBGkIWQXeIg59Q5oPqev7rkXixPakX470yeks2emXSK2+dptZEviBtzXbbko/+NiDCAohsynkkPrbteGrS6b0mDwzdR3ftYln28delYwqnM5YWY4VhEFb4/ckSdJ67LELvb8ZQwcwzX3tQfcWJgApGgf00Xfd0HfpM8teGsnIRzVj5cU0sFp869Kogk/suZDcfpLQnYceZHYY03ZHkxu85aXQ6AIIJRhQRauHewpQTdsfXfk2gY+/LDYjSd0EZgdRCt0bxJGfRPQxSaErICZrdK7pOiEdUHlWVL4TyajcqRBWgW0vtmXWtdLPDw0d+lrqoJR1JSF0aT7EmjOmf2FRSGYYwrOobieSUTFSU7f2sNZXBMh8NhZyCR1UWg7Mp6iPQmiGAXgvmUYAzHFDM4VnnEejANWn2ur2liFBwtgsuUMPkjuFhkmBwfCNRjVsMg1odqVO6IZ0DEDiRl+ZbFSu/+mg9XdlHzlKIblDCzq2EFoc9RQA8x0IfA8ao2eaJoUEo4E1fC0ZlbVQitlcHoSS56aAxIZkN/OJiQ2obDarIyaNIs+H9gZjALJ241dwRXw7NXbPglCiYyDIbEhu6iA2EPVCKKKw7s0ZRzBQosnlQpBvWPFRgGrjUgwbG4ZOymbJxcQEujpARRsaUAmm9PsIXN90Xb+twZuK+yysX2NPoPXQzx9DZjusY4OJiQuoSGrCd2BNEzqfc2lCpWa7tXwVqNapsdvQbSas3wph+k87tuB4YuKzgaQ39unYACa4aSiiuqoCVbNgtyB1rBSmeL5l5xMTqLhZMJBMUjePrIiNLFM3FDdkU3KjE85gj2AEHrsmzNXseG5iAsVJawzNkDgw/Q/NxBynok9pKEFVG+yxDxho1zWbktAm50LR1mMlJSiZ0J4t8TwztJ1xulmlGDWqimBn2QNjaw9CyYygY6PJQ/MbKpIpjA3oxCSI6dzFpoOAQOJBq0G1MtghrIkDpVkyJcXEBCpuHlTinBScGE2diI4B7CBMWyvLI0e1FOysgBk4oQ1OyivJbGISAarEhZGANzZt+jRJzNBRdD5f0f2rZsFuw+e2Qhg2+lYEE5MITr44CQCLxmwMEKZJbLmJzSpKCqUIVRHssRuErMyeT0xCMBj8Kq/sknAcRFbqR35oCizbIEdVBDt0bx6cnoSKWwgTXS91Wfc2piYLnSgw8azKUoUqD3basVEUngvTXTpEgk4/mDq2ZZpx4qx/ClSMVKEqgt2NoMNLxnn3FgaWn8a2b8oZKDWUstsSfuZZeIxY9xbR+b/luxFtb4FqIOVShip4DpLaS8e+HbgWbW+tnRMUJXU3u/zs5+iQgE5ObL+1ipJYqUPl02DykS5FWyl1qFqtscuQwrvN4hkx8UkhKg2WOy5I5T2MNVvHrhKVZuvYVaLSK9j3lN5EXKNgf/6FK0pRzQV7hHnN/97+lasHim9NfxLsCd4rSaihDq8R1Xfxxx/suaGY1KI6Cfa2r45oRzNDMSn+vyGOg73lqyPa0JyhmBSjsmc+avnqiM21YCgm1f/jyPTzxZN2r45YqcN90AHHlsuGYlKOaloUh9u9OmIjnTYUk3pUBat2r45Y0rX96wfX+TZdaSgmBKgKVm1fHTGvFyinQ54NSwzFhAHVVPQJmqMbBwelBE5UbigmFKhEszp4kWOjKkMx4UA1FVu5OgJUlT1fjaGYMKByhZ+N2H/+4PpR+a9rDcWkHpX5BQkDz71yEjyGYlKNKojV1tj5DMWkGhVIdo396KDQzS9yGooJAyrZpZjDtwt96cuchmLCgEp2jf1aQeorR43+DAMq7hp7Wx3lVUbqRjNSOFBVBrudV5LhaTvlrL39r65DCgeqymAPE494bBW/10ZXybq8K2uQQoKqIthjgJSvUk82Xjs7G0MdrkEKCaqKYIcYc3I7bWqqkzHUtTVIIUFVHexOfguFOKAP/nQcEW+NCiD3oLxUSFBVBXuS5klGB/YksSCxmhciGgzKS4UEVUWwM1KQUmwVMoHrcRuaanNDMWFBVRrs1hhuj0ThRMdearZ2tA1DMWFBVRbskTsz1AmgBsOrlgzFhAVVSbDHY5MOP12G6vg17vbXmqGY0KCqK8XYx4ASvvlNm4ZiQoOqrhTjr3hWoXYNxYQHVXUpxjsBFJA6tW4oJjyoqksxMX9VQYChmPCgame5oxhDMSFC1UKNXZShmBChgmBPivNcsV90eEmxjMEqlj567Lp59kK+pXd8P16BhkKHyg8cM2AXLUdukt9nwnF8NgO0ixe81HfYkqwkDdjgFO4Jx9wo1FDoUH3t63BPXIrGoxQ8dhMYSgFujxdRgPAigUvnoSRj0a0iGF4AQDMg8Te+eVP04WFCla9v9EOLGSh27YDda9dxbHbpd+TaeeEq9G3mqCS14H4xlGqQfuvwygvPv3NQpXn62MX13ZFTTPuCom+MzHxw5flOVGwZ5C/AHazJwdUr+3vvDFTRxuuLrwm1FiJUThtLh9a21n7tFnhQeS2dOF3TWjqhanHR+jrW0glV2uYiq2bW2tvf339xf79mYS0aVJuH+pIaWksjV7US6ovisdbxMm19ULUV6kuqtdbxMm19UAm7EqnaWpzLtJmwoGo11JdUYS2uZdqFkKBqPdQXVWqt+mXaJ0KCSkCoL6nEWjXLtOeFA5WgUF/UamvtcRe5cKCSdXkptda3b87WFzctBOJAJTLUF3XtO7PlxW9riUpwqC/qekHqatM/RIFKfKjPq2DV+IwFBlRSQn1OjNV3GxcfMKCSfs8AYHXYuK6FAZW8UJ/p4AYL9WbFBwSopIZ6oaNidtzEWghQyQ31U+K2lnpUskP9tDitpR4VihvB8FhLPSr5ob5S9dZSjkpFqJeoxlrKUSkO9UVVWks1KvWhvqRya6lGhSLUF1VmLdmolm8UhyTUl7TSWrJRLd0oDlGoL2qFtaS7avFGcahCfUnL1pKeVQs3ikMX6otatJZ0VAs3ikMY6kuas5Z0VAs3isMZ6os6tpb8wcLcjeK+933VHPiUW0vluCp67iXVEHgF1lKIKvrBFH1UzelIHapoPNUKlboGCKQeoOIRI/UAFYdyUg9QccgOGSqsM8CVUpZV8Us/fICKU2mchA9Q8YiVX7T675GUocJcflktVaiQl19WSRUq/OWXU1KFSofyy5JUDUG16vtyKUKlX6irQqVhqKtCpWGoq0KlYagrQqVjqCtCpWOoq0GlZairQaVlqKtBpWWoK0GlZ6grQaVnqCs5Ea9nqKtApWmoq0ClaagrQKVrqCtApWuoK1gLqmuoK1hhrGuoy0elbahLR6VvqEtHpW+oy0alcajLRqVxqMtGpXGoS0alc6hLRqVzqMtFpXWoy0WldajLRaV1qEtFpXeoy0CVzaR3qAtHRQmRDP5lz1V/2M0kFtUcJ5Jz01giUWULnMiMnK4S7arTP9pKHKpsFSedfSUQ1UpPaZxXwlBlJT/62koUquzssRKHqpTVA1SLKuWUacuKG9UzD51v8LaV7U9fVD/6MY9+kmU//RnXlqCsKtf1RfXkhEdvwIz3H3fucW08qcwpfVE9dZ9Hb+blgX/+6988W2dnM6uevsWjt3JU//nv/3i2PqNZJagBnsWs4o/1nzeJ9bOYVXwbNhwsVGbVw43eC4uaDUF3ubfMqrJqcvvWxQuqP3ljNUPVLf69XLtlVa6/+jJNszvaWasRqk4//3dYnzZVc8D7r9y9dXuinbUaoert5MSGWWdtVln2i8nk3mt39bNWI1QzU5HMqEeVlbW/xx+BIUcTa3W63a5BjG63hxtVZ7htjLrgp90he8EwyKV6VGXtD/557OKjrzez1hZ8S9sce1WKqke2domxTZ8N8w6QArvM8/VW1arON7RWb0D3W9/qFaPqdOhh9qDPy9vfqGcYW93aNyYr2+Bcf9DMWpd6XcWkeLJqRBH1qZV2mJU60Aq7XKjqTpk2sVZ3sKOYFA+q7pB0BvQrzb9Wllfd+oFVwWpuKLpq7sdtrV6mA6qtodGH4wQnGf2BAc1xwGkrMre8Y/Wv+aw12tnu40fV2crnM6ORqGOot5YxIsYl/Fk1KgzE76PmqrHWDnxLg6FKToQDldHvs+FMRyQqUmUtow8j3tFWg0avBtVMu/xVhTXVdKwlW6r/c5ZFNR/GSxQuVKithQ0VwWsthKiwWgslKoLSWlhRIbQWXlQEm7VQo8JlLeSoCCJr4UeFxlo6oCI4rKUJKgzW0gYVUW4tnVAptpZeqIhKa2mHSp21NERFFFlLT1QrrPWY8F3qioosW+uXwvenMaoFa/3q18L3pjUqcmKt34jfle6oZtb6rfgd6Y+KgLV+J2Ev/wdNm5TrFpZbwAAAACV0RVh0ZGF0ZTpjcmVhdGUAMjAyNC0wMS0zMFQxMzoyMToyNCswMDowMPVEUgcAAAAldEVYdGRhdGU6bW9kaWZ5ADIwMjQtMDEtMzBUMTM6MjE6MjQrMDA6MDCEGeq7AAAAKHRFWHRkYXRlOnRpbWVzdGFtcAAyMDI0LTAxLTMwVDEzOjIxOjI0KzAwOjAw0wzLZAAAABN0RVh0cGRmOlZlcnNpb24AUERGLTEuNVM2WsEAAABKdEVYdHNpZ25hdHVyZQA3Y2FhMjJhNjExMGE2NDBkYzkwYzA5OTkxYTZhZTUyMzFkYWVlZDdjMzYxMDcwYzAyNzc0NjJlYmI1MjQ5NWQ0yGoRtwAAAABJRU5ErkJggg==", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAASoAAAEZCAMAAAAXCpgeAAAAIGNIUk0AAHomAACAhAAA+gAAAIDoAAB1MAAA6mAAADqYAAAXcJy6UTwAAAG8UExURf///wAAAAAAAAAAAEtLSwAAABsbG4qKiktLS2pqaq+vr46OjgAAAAAAAAAAAAAAAD09PQAAAAAAAAAAAGpqagAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAIyMjAAAAAAAAAAAAAAAAAAAAK6urpKSknd3d2dnZwAAAAAAAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAAAA/wAA/wAA/wAA/wAA/wAA/wAA/wAA/wAA//8AAAAA/wAA/wAA/wAA//8AAAAA/wAA/wAA/wAA/wAA/wAA//8AAP8AABwAagAA/9IAHgAA/wAA/wAA//8AAP8AAP8AAP8AAP8AAEdHR1VVVRMTEx4eHkFBQTo6OrKysnl5ec/Pz4ODg01NTY2NjWdnZ3d3d93d3f///5mZmaqqqiIiImZmZlVVVYiIiO7u7ru7u8zMzAAAAERERDIyMhkZGXBwcK+vryAgICgoKF9fXxERESQkJBsbG1JSUv8AAAAA/zMzMwwMDFdXV6KioggICDg4OGdnZ0Bn/HoAAABydFJOUwAAEXetRKA/5tr+mSHHdqPW9dbj9TP6tmb+XI4/7feIVd/vmMfO1qiqmapEVZndEYh3uzMiZszu6+i00rCY2Y/45N8iVardiJlEdxFwZjPuu8ZdkP2S8+Yx4enPYKCt8VzbuO/3es3y69fa9I/+Iexc0nG6YLsAAAABYktHRACIBR1IAAAACXBIWXMAAABkAAAAZAAPlsXdAAAAB3RJTUUH6QEPDR4VKdiRawAAERNJREFUeNrtnYl727YZxoN0ybpNomq5lRrLWXe03cr71u52cxvHzuJ02X1v3X1vidK0Tcp2W3equ//h4QMpW5JFEpQI4IOb93ms6JFpkfrpxQvgI8icO3eOaK/zD71Lwl7O6Y/qwsVbtyfvFr8f3VGdf+jhO5PJ5OX3iN+V3qiYoSa3b9195b3id6YxqpmhXr372r3J5H3C94cMVafb7RrE6HZ7dVueGOr+ZHLvziMXhB8bMlSEbPXpw7ZRvdGioV5/9OJjEo4MHaregJBhp3KTZUM9fl7KkaFDRS71ulWklBiKCR+q7mCn/JeKDMWED1UvK0OlzlBM6FCNdrb7K3+h0lBM2FAZI2JcOp1Vig3FhAzVzog+DIZLryo3FBMqVEY/owOq0dagO/ciBkMxoUK1QjgMxYQaFRpDMSFGhchQTFhR4TIUE05U2AzFhBAVQkMxoUOF0lBMuFBhNRRTM1S7Qo8Fr6GYmqGajaIvt38gqA3F1AhVp5jzD7O2DwO5oZgaoerllaTOMOvw/1G98BuKqRGqmalIZvD/UVao7Pc6GIoTVWe4bYy64KfdvDhiGOQSL6qCUZY/Pf17TQzFiapHtnaJsU2fDfMOkAK7XHuabkZqxon9LLPSxlCcqDqdAeUFfV7e/kY9w9jq1r5xAWqJ1RwsnQzFiYqMKKI+tdIOs1IHWmGXD9WKn5n0MhQvqu6QdAa0y8tPz7G86nIMrLLVrOBHO0PxotoaGn1IdXCS0R8Y0BwH9bYq4ZTpaChOVJ2tfD4zGjV646yU1fv1MxQnqlFhIK54miNV8kMR6mcoPlRGv88GUZ2GqMpY0d/oZyg+VDPtNqsqlLY/YKWdoZi4UT3xgSZvm1X8kEw7QzFRVBee4NEHsw91P8y1JSiryCqSaWeoGaonJzx6A6ZzkMZcqmp/p+Y3moiieuo+j97M6wO//8MfebbOqtqfvqievsWjt4payp/+zLN1dVap/tBro2rQAP8C42yuBliZVao/9NqoPvIaj/6aZX/7O9eWoDOaVXwbfvSZJh38Gc0qEW+bZWcyq4S8b9UcUPVnXlPCUJ259ifuRHxlH6ilxKEqr+1pKnHLO7Kz5SmhK2GyM+UpsYuGsrPkKcHrq7LlU6Y6kxK9FG2eld6gJKzay45lWqo/7GaSuMDRdlR/2M0kcy2o66n+tBtJJio/UP1p11EQRvkTmahi9/ip5edCHV9JmNBH1yoOW+pi7JNgj61csWocq2WlPn10Y5ceX0gcdtiJVFQaBHvsmh58p6ZNSMoiwyW59125S/xPgt0LnIjEYWvvbPu2n2z0Dl6aUg/5CeAJ6ftRVB48mgkLq8SRi+ok2BOwdeK39cYWhQ5GWE/hmPom8W2THqEF359jWfDccRKGEHYwjuWiij82e+Z5LilCoA2ZyXpjEX9M/zAyI3owvk/GtL8LqKWIbbJv0aQvkChNU99PI9nX2Hz8E7NncCzuJm+1oGlM4rTRXwRjCE7Xo4Atn4QeBGkIWQXeIg59Q5oPqev7rkXixPakX470yeks2emXSK2+dptZEviBtzXbbko/+NiDCAohsynkkPrbteGrS6b0mDwzdR3ftYln28delYwqnM5YWY4VhEFb4/ckSdJ67LELvb8ZQwcwzX3tQfcWJgApGgf00Xfd0HfpM8teGsnIRzVj5cU0sFp869Kogk/suZDcfpLQnYceZHYY03ZHkxu85aXQ6AIIJRhQRauHewpQTdsfXfk2gY+/LDYjSd0EZgdRCt0bxJGfRPQxSaErICZrdK7pOiEdUHlWVL4TyajcqRBWgW0vtmXWtdLPDw0d+lrqoJR1JSF0aT7EmjOmf2FRSGYYwrOobieSUTFSU7f2sNZXBMh8NhZyCR1UWg7Mp6iPQmiGAXgvmUYAzHFDM4VnnEejANWn2ur2liFBwtgsuUMPkjuFhkmBwfCNRjVsMg1odqVO6IZ0DEDiRl+ZbFSu/+mg9XdlHzlKIblDCzq2EFoc9RQA8x0IfA8ao2eaJoUEo4E1fC0ZlbVQitlcHoSS56aAxIZkN/OJiQ2obDarIyaNIs+H9gZjALJ241dwRXw7NXbPglCiYyDIbEhu6iA2EPVCKKKw7s0ZRzBQosnlQpBvWPFRgGrjUgwbG4ZOymbJxcQEujpARRsaUAmm9PsIXN90Xb+twZuK+yysX2NPoPXQzx9DZjusY4OJiQuoSGrCd2BNEzqfc2lCpWa7tXwVqNapsdvQbSas3wph+k87tuB4YuKzgaQ39unYACa4aSiiuqoCVbNgtyB1rBSmeL5l5xMTqLhZMJBMUjePrIiNLFM3FDdkU3KjE85gj2AEHrsmzNXseG5iAsVJawzNkDgw/Q/NxBynok9pKEFVG+yxDxho1zWbktAm50LR1mMlJSiZ0J4t8TwztJ1xulmlGDWqimBn2QNjaw9CyYygY6PJQ/MbKpIpjA3oxCSI6dzFpoOAQOJBq0G1MtghrIkDpVkyJcXEBCpuHlTinBScGE2diI4B7CBMWyvLI0e1FOysgBk4oQ1OyivJbGISAarEhZGANzZt+jRJzNBRdD5f0f2rZsFuw+e2Qhg2+lYEE5MITr44CQCLxmwMEKZJbLmJzSpKCqUIVRHssRuErMyeT0xCMBj8Kq/sknAcRFbqR35oCizbIEdVBDt0bx6cnoSKWwgTXS91Wfc2piYLnSgw8azKUoUqD3basVEUngvTXTpEgk4/mDq2ZZpx4qx/ClSMVKEqgt2NoMNLxnn3FgaWn8a2b8oZKDWUstsSfuZZeIxY9xbR+b/luxFtb4FqIOVShip4DpLaS8e+HbgWbW+tnRMUJXU3u/zs5+iQgE5ObL+1ipJYqUPl02DykS5FWyl1qFqtscuQwrvN4hkx8UkhKg2WOy5I5T2MNVvHrhKVZuvYVaLSK9j3lN5EXKNgf/6FK0pRzQV7hHnN/97+lasHim9NfxLsCd4rSaihDq8R1Xfxxx/suaGY1KI6Cfa2r45oRzNDMSn+vyGOg73lqyPa0JyhmBSjsmc+avnqiM21YCgm1f/jyPTzxZN2r45YqcN90AHHlsuGYlKOaloUh9u9OmIjnTYUk3pUBat2r45Y0rX96wfX+TZdaSgmBKgKVm1fHTGvFyinQ54NSwzFhAHVVPQJmqMbBwelBE5UbigmFKhEszp4kWOjKkMx4UA1FVu5OgJUlT1fjaGYMKByhZ+N2H/+4PpR+a9rDcWkHpX5BQkDz71yEjyGYlKNKojV1tj5DMWkGhVIdo396KDQzS9yGooJAyrZpZjDtwt96cuchmLCgEp2jf1aQeorR43+DAMq7hp7Wx3lVUbqRjNSOFBVBrudV5LhaTvlrL39r65DCgeqymAPE494bBW/10ZXybq8K2uQQoKqIthjgJSvUk82Xjs7G0MdrkEKCaqKYIcYc3I7bWqqkzHUtTVIIUFVHexOfguFOKAP/nQcEW+NCiD3oLxUSFBVBXuS5klGB/YksSCxmhciGgzKS4UEVUWwM1KQUmwVMoHrcRuaanNDMWFBVRrs1hhuj0ThRMdearZ2tA1DMWFBVRbskTsz1AmgBsOrlgzFhAVVSbDHY5MOP12G6vg17vbXmqGY0KCqK8XYx4ASvvlNm4ZiQoOqrhTjr3hWoXYNxYQHVXUpxjsBFJA6tW4oJjyoqksxMX9VQYChmPCgame5oxhDMSFC1UKNXZShmBChgmBPivNcsV90eEmxjMEqlj567Lp59kK+pXd8P16BhkKHyg8cM2AXLUdukt9nwnF8NgO0ixe81HfYkqwkDdjgFO4Jx9wo1FDoUH3t63BPXIrGoxQ8dhMYSgFujxdRgPAigUvnoSRj0a0iGF4AQDMg8Te+eVP04WFCla9v9EOLGSh27YDda9dxbHbpd+TaeeEq9G3mqCS14H4xlGqQfuvwygvPv3NQpXn62MX13ZFTTPuCom+MzHxw5flOVGwZ5C/AHazJwdUr+3vvDFTRxuuLrwm1FiJUThtLh9a21n7tFnhQeS2dOF3TWjqhanHR+jrW0glV2uYiq2bW2tvf339xf79mYS0aVJuH+pIaWksjV7US6ovisdbxMm19ULUV6kuqtdbxMm19UAm7EqnaWpzLtJmwoGo11JdUYS2uZdqFkKBqPdQXVWqt+mXaJ0KCSkCoL6nEWjXLtOeFA5WgUF/UamvtcRe5cKCSdXkptda3b87WFzctBOJAJTLUF3XtO7PlxW9riUpwqC/qekHqatM/RIFKfKjPq2DV+IwFBlRSQn1OjNV3GxcfMKCSfs8AYHXYuK6FAZW8UJ/p4AYL9WbFBwSopIZ6oaNidtzEWghQyQ31U+K2lnpUskP9tDitpR4VihvB8FhLPSr5ob5S9dZSjkpFqJeoxlrKUSkO9UVVWks1KvWhvqRya6lGhSLUF1VmLdmolm8UhyTUl7TSWrJRLd0oDlGoL2qFtaS7avFGcahCfUnL1pKeVQs3ikMX6otatJZ0VAs3ikMY6kuas5Z0VAs3isMZ6os6tpb8wcLcjeK+933VHPiUW0vluCp67iXVEHgF1lKIKvrBFH1UzelIHapoPNUKlboGCKQeoOIRI/UAFYdyUg9QccgOGSqsM8CVUpZV8Us/fICKU2mchA9Q8YiVX7T675GUocJcflktVaiQl19WSRUq/OWXU1KFSofyy5JUDUG16vtyKUKlX6irQqVhqKtCpWGoq0KlYagrQqVjqCtCpWOoq0GlZairQaVlqKtBpWWoK0GlZ6grQaVnqCs5Ea9nqKtApWmoq0ClaagrQKVrqCtApWuoK1gLqmuoK1hhrGuoy0elbahLR6VvqEtHpW+oy0alcajLRqVxqMtGpXGoS0alc6hLRqVzqMtFpXWoy0WldajLRaV1qEtFpXeoy0CVzaR3qAtHRQmRDP5lz1V/2M0kFtUcJ5Jz01giUWULnMiMnK4S7arTP9pKHKpsFSedfSUQ1UpPaZxXwlBlJT/62koUquzssRKHqpTVA1SLKuWUacuKG9UzD51v8LaV7U9fVD/6MY9+kmU//RnXlqCsKtf1RfXkhEdvwIz3H3fucW08qcwpfVE9dZ9Hb+blgX/+6988W2dnM6uevsWjt3JU//nv/3i2PqNZJagBnsWs4o/1nzeJ9bOYVXwbNhwsVGbVw43eC4uaDUF3ubfMqrJqcvvWxQuqP3ljNUPVLf69XLtlVa6/+jJNszvaWasRqk4//3dYnzZVc8D7r9y9dXuinbUaoert5MSGWWdtVln2i8nk3mt39bNWI1QzU5HMqEeVlbW/xx+BIUcTa3W63a5BjG63hxtVZ7htjLrgp90he8EwyKV6VGXtD/557OKjrzez1hZ8S9sce1WKqke2domxTZ8N8w6QArvM8/VW1arON7RWb0D3W9/qFaPqdOhh9qDPy9vfqGcYW93aNyYr2+Bcf9DMWpd6XcWkeLJqRBH1qZV2mJU60Aq7XKjqTpk2sVZ3sKOYFA+q7pB0BvQrzb9Wllfd+oFVwWpuKLpq7sdtrV6mA6qtodGH4wQnGf2BAc1xwGkrMre8Y/Wv+aw12tnu40fV2crnM6ORqGOot5YxIsYl/Fk1KgzE76PmqrHWDnxLg6FKToQDldHvs+FMRyQqUmUtow8j3tFWg0avBtVMu/xVhTXVdKwlW6r/c5ZFNR/GSxQuVKithQ0VwWsthKiwWgslKoLSWlhRIbQWXlQEm7VQo8JlLeSoCCJr4UeFxlo6oCI4rKUJKgzW0gYVUW4tnVAptpZeqIhKa2mHSp21NERFFFlLT1QrrPWY8F3qioosW+uXwvenMaoFa/3q18L3pjUqcmKt34jfle6oZtb6rfgd6Y+KgLV+J2Ev/wdNm5TrFpZbwAAAACV0RVh0ZGF0ZTpjcmVhdGUAMjAyNS0wMS0xNVQxMzozMDoyMSswMDowMDB94E8AAAAldEVYdGRhdGU6bW9kaWZ5ADIwMjUtMDEtMTVUMTM6MzA6MjErMDA6MDBBIFjzAAAAKHRFWHRkYXRlOnRpbWVzdGFtcAAyMDI1LTAxLTE1VDEzOjMwOjIxKzAwOjAwFjV5LAAAABN0RVh0cGRmOlZlcnNpb24AUERGLTEuNVM2WsEAAABKdEVYdHNpZ25hdHVyZQA3Y2FhMjJhNjExMGE2NDBkYzkwYzA5OTkxYTZhZTUyMzFkYWVlZDdjMzYxMDcwYzAyNzc0NjJlYmI1MjQ5NWQ0yGoRtwAAAABJRU5ErkJggg==", "text/plain": [ "" ] @@ -389,7 +223,7 @@ }, { "cell_type": "code", - "execution_count": 25, + "execution_count": 4, "id": "1a7fde5f-544c-4a8b-9f90-5d43c781a3c2", "metadata": {}, "outputs": [ @@ -397,13 +231,13 @@ "name": "stderr", "output_type": "stream", "text": [ - "/home/apn/.cache/pypoetry/virtualenvs/pyfeyn2-rvdu2REa-py3.11/lib/python3.11/site-packages/feynml/pdgid.py:111: UserWarning: Type remnant is not consistent with pdgid -2, which is anti fermion. Using remnant now.\n", + "/home/apn/.local/lib/python3.13/site-packages/feynml/pdgid.py:111: UserWarning: Type remnant is not consistent with pdgid -2, which is anti fermion. Using remnant now.\n", " warnings.warn(\n" ] }, { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAASoAAAEXCAQAAACw18HXAAAAIGNIUk0AAHomAACAhAAA+gAAAIDoAAB1MAAA6mAAADqYAAAXcJy6UTwAAAACYktHRAD/h4/MvwAAAAlwSFlzAAAAZAAAAGQAD5bF3QAAAAd0SU1FB+gBHg0VHecOI80AABJYSURBVHja7Z1fiCzZXcc/s2zMEqKm8CVyJexUosaH+LC14B9crtEKxCxBVGofwiCSh74gCEuEaZMHkQSl98EJEiJ0+7JyVwg9aDDkDoFuJBvQQeghIjgg2Z4YbjYq5k4H9cabu7tTPnR1T3V3VU1V9zlVp875fQa2t/9Mz+mZ7/2eX/3O7/wOMTGCIURETQ9BBU80PQABAJ8eU0JmTQ9ECeJUjRMxZEoXr+mBKENE1SBzf+oTNj0QxYioGsI+f7pGRFU7tvrTNSKqWrHZn64RUdWE/f605MmmB2AgaR8ZM1bwjhERAQOetSRlcBPiVFpxyJ+uEaea49Hlggt8Bsre0zV/ukacCoAhHQC6St7NSX9KIaICAi4JCZVclblxfVeMiAoIGSl4F9f96RoRFRAsRbWtIMSfUkigDnDGGREzfM4qf69Ph4gxAyXJBzsQp0rwt3AZ8adM9mJgr+lRtI+FPx1X8CePAAh5N08B3+OTtqYaZPq7mXm8NeMiuf9+fpof2iL/NAP+gqcBuKJnq6SQ6a8EQ+K1r89u+U4Bfb5HzENebPpDaUVEdSPRmqQ6W7yHR4cJIzp8jTf5naY/kmZEVDfiM1sK6nKLpENAnyk9fMDjv/hg0x9IOyKqQubXd19dSiqo9N3X/rTAr/gO7URElUM6P+5vIam0P7mGiCqDzfzTsIKkNv3JNURUK+St34UlJeWyP10jolqyW35c/GlJXclP9SW6Ktl1/S6gQ8gxLywTpG7jvFOJP6lHo6g8enQIjf2F71r/JPFTHhpFpbZEVy3iTzrRJiqVJboqEX/SjzZRqSnRVYv4Uz1odKpdS3RVIv5UI/pSCruU6Kplt/13HhEdZhxzp+HP0R60phS2KdFV+/PFn5pAc54qaOxP0kT8ZOJ1bhNoFlVITMyw1gapdftTSC/5mkgQD9Qkqpi4pl0nTV7f+eJUc+a7aXxtU9QzvLRy/x73+IaWn/ROfpmPVtzfkmaxfjeQ9btdmV/9+dou+59eu/+r/BivavlJP8Uz/BEvb/Gdu13feXSBP9DymdpKbdNfrD1tGNFnUvF7VFzfmZjmbZZaRHVJv6ZrwPJ/XnX58S69Wj5be9Auqimd2rJVPtNSr1ObfxoasWZgElbVUw1L+I769btLIKJTeeq1Fpu2vQd4JZor+orrMwPOiLhgtvRjjw6sXfW6hUVO1WskT9RlshZT9QjtOA1raywSVVRr/LZgSJiqyICACZFhNWR1Y5GoFgs09QbOlwBMIXGngAmhiKrpIRQRbeE9QclrQBUEDAHoEyU7Az3JWpkuKhgR87WK3jOpsTJiIfnFTwwla1Xm6q/TaD3Rv/BzPMc9HvFFPl3qui1KNSjTz6Lwb/ETw8aLEhunjKjOGl1ifSe/DbzF9/nPkpWbPT7U0Fg9PDzjtsrWj/HT35RHfKFCv5UO/cbGGjSU1jANw0UVVM749CSmaRzDRVUdjxETek40FzOV2kRV57Tg02GUXOwLtaO7j3q4TAaEDBQee1aGCXfkSqwRanKq+uu3PaaO57Wbo+UxVX6qUQL25tAsKk/zH9fLWZLxxaea4wnN7z9jrPU6bMY4s+Au4tjiYzoMR7eo5kVsOhlkZrIuUpcIQt1oj6n0F6JkLx+HDJnSl3xVA2gXlf4K7k5B1BYxlfYataO9QceIiGCZUtBRZusV1k+NxKvqRndMFeIRcMZFshFAxxVZXrAOEBjQH8s9NDvVagV3oGkqCnKrLUcSrjeAZlGlK7jDgj/+ruQF61La2wB6p78gKVgbE3HBGb62ArZB5gQoU18zaM+oz5n7yFDbeXfZwbrHhOGW5dByzbg19Z72PmSsLdPd5yyzCiIgIuK4oNlPtCL0xet60h5oa1q+oHxNcbxWlIIdJdvAvFR7Ra/BouT2Y42oijdm5adA/aWQ0tvXO5Ld2gGLRJWfWe8W+M7CwforrxGf2gWLRJWXWS9TrtdfKT6+LinsLc9O9hxvulEFi0QF/czEws17hjtMVmTXTSbLDiGL+EqaWZfHKlFlB+vBDdvgryW1iKP6y++cE+GLT1XAKlHlBesREyY57RhDLpPH/WXflk1X6kkdaQUsE1V+sO7RybgGDJhu+FOWgKTevQL6Kz/r5Th3mpoxYLwmKp8RZ0T06K0E+evpWVnuqYZlTpUXrEPZ3sVRhiyb7XvTPqwTVX5mvV/qCq5X8jEhF9umPzjDy/QVn7DEDunscFz25VTDOqfKC9bLdbjTd/STS1goqrwymCl9yTbVgX3TX17N+oz3MiZkKhGSdix0quIyGE/OktGNjU6VH6wDzPAk8NaLnaLKq1mHeamLpDL1YuX0V7TBdCrld7qx1anyNpiGjMWndGOrqPK6wYAv9QbasXT6g7wymPmRSF2RlkYsFlV+GUxAT7rBaMRiURV3g+lJtkobFouqqC5BehfrxGpR5WfWpXexTiwUVbpje3aw7q3tnhGUYs9p7+mzJWbLyqkBnYyeCAFjWarRiHVOtXq2RHaw7jPN2VsjqMA6Ua2THax7dBkxldP5tGCNqHpERPQzNmEVlcFIm1kdWCKqHh3yJrui/cl9yVZpwApR+Ulv0SizW0t+Zr3cpi2hKlaIqpv0bOnlxE950pHexVqwo0phlpxHH2Y2qs0rgwlATmbXgR2iOsYjpIeXiGud7DIYKSvWhRXTH8wbxub3v8srg5nIti0NWCOqvIhqTn6wHtGXRRvFWCOqHsOCLLmUwdSJNaK6iaL2HBNZslGKM6LKz6x3pBexYpwRVb4fSWGxYuxIKZRjkLtp66LyewlFOORU+WUwHbn6U4lLTpWdWb/gQ/hMGEm+ShkOOVVxGYzPRMpg1OCSUxV3g7lgLJOgIpxyquI+67JpSxWOiSo/s96X0mJlOCaqeWbd2xCQJ+V66rBni1YZfC4Y0ONXNp4JpMGQQpxyKo9zRkx5gy9sPNdjSleu/5TglKgg4BFvEWeIar5jcMRIwvWdcUxUc1nFfDn3+Z4cF7krbuWpAM74NeA9uc8HshK4M845FcCLuWF5WJBzF8ripKjIzaubvmnLo0ufgPlZhYbi3vQ3J3uKM3/TVsRLzM/6Cg1OgjjqVNn4TI13qkWa1uQN+yKqFUKGTI3+gy1WLw1eAXB1+stjzAu8lzE9Y6urFifTX4CpYxRRZXG8cYS3OYzxCImArrHxn0x/mZi6GSKApCbMMzjzL6LKwNzexXErKupFVBmYW64XExNzmdEx0CREVBuEBh/fFqe+RqauUu7FwF7TozCNDh/hGf616WFksJ7qeIOv899ND2oDcapMOoZGVWmnMrdtt4gqE1NbdlxPfYbmqJajFNYIjF2sjbk01p/SoxRRbdA3NQRuQToBRFQZyM6aHZFlmk0ijpseQssRp9rA1CC9NYhTrRMsu7ILWyKiWqcjkx/s6NUy/a0gQfqcgO/scMKYiGoFUzPp9fNB3uT7zHipumvJ2t8qE17YMqIyNV26PQEeEPMf/FLF34k4VQpzM+n18ylifrDlFCiiSmFuJr1uIka8uPV3i6iWSJCuCEkpXCOZdFWIUy3Jy6R7Bu4CNLqVpFud9IrIz6T/LReGbIYKl/IOmTFoeji5iFMl5AXpX85pkdYkvtlOJaKakxek/zWPiCUhWg2Z/uZkB+l/zG/yFhjRX8WjQ8AAj4BjI0aUjzgVkB2kR/wDn+HbxEYE6p3lKEOGyWMeXRMnQkkpQF6Qfswv8jY+Z0iQfoyXnGbvL0uKu5yZWKYjooKicpeIAXeaHh4AM8JE3sHyNsQzRPIrSEw1z0NlCydizMyYor15i1uPkGeXj82aHlQW4lRFmXSzcuwhENLjhURKF8wYGykrCdRzM+m+UWuBHhNWKzJDU1Md4lT5mfSOUTnreUR1sfKIoYkFiamKgvRnK72TTsJk8luE5R6emUE64Pz0l1/uEtFvenAFBPRMzFAlOC6q/Jr0oREpz3ZiuahuOiGhHUF6y7A9UC8+IaEtQXrbsN6pik5IyK9JN7frZxuwXFRFJyS0NUg3Htunv6ITEtqSSW8f1jtVSJ+QLiO6a1OaBOm6sF5UeSckHPCNnNebnAFqBw6IKouAR3wl5zkJ0nfE/pgqi4Axb+efMp+LDF35bxEuiipgxLtyn5UgfXecm/4CLpNe5FkLNBKkK8A1pwoYFUZMkklXgVNOde1SeU4lQboC3HIqf2X9b7PzkgTpanDKqQB8vsL/Jee7rCPlLmpwTlTzTHqH0YaoJEhXhXOium7BuB49SSZdFc6JSspd9OOYqKTcpQbcuvqTcpd6cMyptit3kVirEm7t+6tWk96WZojm4ZRTbRekm94M0TwcEpUE6TXh0vQXFXahWqVNzRDNwyGnqhKkbzZDjBBK4k5KoVqQvtkMURKjpXFn+qvW3WVGtNIMMTC2w4qJODL9VQ/S5yuBHlM8Qjm0rQquTH/VM+npZohn+OJUFXDEqapm0tebIQ7xtz5R2DnccKrq3V2ymiGa0aO4DTjhVFUz6SH9zC4xQjkcEJVk0mvGhelPyl3qxgGnku4uNWO/U0kLxvqx3qmkJr1+LBeVBOkNYPv0J0F6E1juVBKkN4DdTiVBejNY7VQSpDeDxaKSIL0hbJ7+8kPxUIJ0rVjsVHlBuqAZe50qMOZAbeewV1QdmeIaQ+P0d7vBj+VJHqpBNIrqdO3+gxo/VsfUk9CdQJuo9jlfuX9S6xWBBOkNoi+m+jR/krq3n/qvfiRIbxZtTrXuUxBzqO0jrH7lZ9KFGlCxQ3mfz/MBPse7+QU+zysA3ObfUq845KugQb3zd9wjZi+5v7d8XNb2mkOBU90FHnIbOOR+8tjJyrXfCQAPuKt46NdC3by1OKtrOiqc6g/Z5228CnyAdySPPc2ry+ePeMAh8BTvUzbuGNjj2qWybxfOJdSKClF9kyNeAyDgHwE4SHV02uf9fASA31A47j0olBMiqQZRMkuccgLs8zi5vjtNXeedpF6lJlMVF0x6m7cyCdaPElE95IRDzjlI7i/Snoec8zi55rvLAx5vJES3G3I5OS0+nsiqbhT8zvd5SHpR5ogjrQOuditeVT8KRHW05j8q3Ch/uCIr49k9o37Ih/nRVFpTf9Z8L5WfKnMr1I3ykOO21uqEqi61+HgirxpRv/b3MX5e26aCuLJL7Ylf1c9eDHz4Uz/yW6re8O9+8sEPw3u++xOzW/+jerDHz6TyU+Vvr1OlQj3ExPzNaayMcLmsG8T9+FLdG8dxtfxU+lamv1p5EuDbb1fXJfVy+X9n3OH3eZ7n+XF14y2XRU/fikfVzpMArz/1urI3vFy59ybf4u95l7rxVpHT/DY9DQr10MLpb5vJT6a/GmljoL5dqC5uVR/K/xH36WpMKWydp3r+OxrHJazSspmhuqySpZ04nsbd+NbL0si6Blonqq1iqigeJnFZXzxLPy0U1RaLNaO4E/txN57GsXiWflomquqyiom/9L9xHMeXcS/2xbPqoIWiqiYriOk8d75Ib2jzrH1OOeWUQ+CQU04Vb/NoE0aIap/XKuwJrDb5Le4Ft17+3YeTWKtnPVzudryvaY9jOzBCVPBFYq54rWTZTPVQfY6n2bPu8hiAk9r2YhvJ3g2LYzqrOFf52WR71w/4cz5x46sXwy5Oe2anPINbv/fr0cffEQAzBgwIiJJztwd86d/v/RkDZlt+iiv+ivfxMb5Z2+/NQG4SVV0m7nHIE8Ab3OMTJf8k6Z19WXuUiz6YR/Tciwc/E+EBY44ZE9HBBy4Y8Mpfvv7KVmeRnhLw8WSXtrsYMv19q8Lktzr47G1Z5T6U+jjrLvFyT5G7GCGq2zs54maDjiqojLOOOOD+WmsSFzFCVM2jwrMOOQIOuXI7TKcGUZ3s/hY1sZtnHSx3Oz5u0WfWgyZRHSapwFMetiwJuI1nHXLOovvWEQ8V7cRuL5qd6nYr/9XqzmfZjsRUudSRg7cTDaLaX059pxy1PGgVz9oGTU61yNYcLtsLtRnxrGpoEdUBV8v/v2plVLWJeFZ5tIjqZNn787ZlGWbxrDJoEdX9xJ32OW9ZQqEM4lk3oUVUV5xzyulaj2K7EM/KR4OoDpzJUohnZaNBVCe1Hm3UPOJZ62gQ1X0HFynEs9IoFtUBp1xx7miFtnjWHFmmUYx4lohKE257lohKG+56lohKMy56loiqBlzzLBFVbbjjWSKqWrnZs4Kmh6gAEVUDFHjWJz/b9OAUIKJqiEzP+vrVZ9pf0iiiapg1z/rod5sekBJEVI2T8qw/taOgUURlCMGtl+/8c9ODUMP/A1ogL1Nbh0MxAAAAJXRFWHRkYXRlOmNyZWF0ZQAyMDI0LTAxLTMwVDEzOjIxOjI5KzAwOjAwlJMzxwAAACV0RVh0ZGF0ZTptb2RpZnkAMjAyNC0wMS0zMFQxMzoyMToyOSswMDowMOXOi3sAAAAodEVYdGRhdGU6dGltZXN0YW1wADIwMjQtMDEtMzBUMTM6MjE6MjkrMDA6MDCy26qkAAAAE3RFWHRwZGY6VmVyc2lvbgBQREYtMS41UzZawQAAAEp0RVh0c2lnbmF0dXJlADBkMjZlMTkwZGJhYTYxNjU1MmM5ZGI3NWRlMWFlODc0ZmYzM2M0NDgzOGQ0YjI4YjE1Njg2ZTc3NzBiNjE5NTUKh85BAAAAAElFTkSuQmCC", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAASoAAAEXCAQAAACw18HXAAAAIGNIUk0AAHomAACAhAAA+gAAAIDoAAB1MAAA6mAAADqYAAAXcJy6UTwAAAACYktHRAD/h4/MvwAAAAlwSFlzAAAAZAAAAGQAD5bF3QAAAAd0SU1FB+kBDw0eF8fW8EcAABJYSURBVHja7Z1fiCzZXcc/s2zMEqKm8CVyJexUosaH+LC14B9crtEKxCxBVGofwiCSh74gCEuEaZMHkQSl98EJEiJ0+7JyVwg9aDDkDoFuJBvQQeghIjgg2Z4YbjYq5k4H9cabu7tTPnR1T3V3VU1V9zlVp875fQa2t/9Mz+mZ7/2eX/3O7/wOMTGCIURETQ9BBU80PQABAJ8eU0JmTQ9ECeJUjRMxZEoXr+mBKENE1SBzf+oTNj0QxYioGsI+f7pGRFU7tvrTNSKqWrHZn64RUdWE/f605MmmB2AgaR8ZM1bwjhERAQOetSRlcBPiVFpxyJ+uEaea49Hlggt8Bsre0zV/ukacCoAhHQC6St7NSX9KIaICAi4JCZVclblxfVeMiAoIGSl4F9f96RoRFRAsRbWtIMSfUkigDnDGGREzfM4qf69Ph4gxAyXJBzsQp0rwt3AZ8adM9mJgr+lRtI+FPx1X8CePAAh5N08B3+OTtqYaZPq7mXm8NeMiuf9+fpof2iL/NAP+gqcBuKJnq6SQ6a8EQ+K1r89u+U4Bfb5HzENebPpDaUVEdSPRmqQ6W7yHR4cJIzp8jTf5naY/kmZEVDfiM1sK6nKLpENAnyk9fMDjv/hg0x9IOyKqQubXd19dSiqo9N3X/rTAr/gO7URElUM6P+5vIam0P7mGiCqDzfzTsIKkNv3JNURUK+St34UlJeWyP10jolqyW35c/GlJXclP9SW6Ktl1/S6gQ8gxLywTpG7jvFOJP6lHo6g8enQIjf2F71r/JPFTHhpFpbZEVy3iTzrRJiqVJboqEX/SjzZRqSnRVYv4Uz1odKpdS3RVIv5UI/pSCruU6Kplt/13HhEdZhxzp+HP0R60phS2KdFV+/PFn5pAc54qaOxP0kT8ZOJ1bhNoFlVITMyw1gapdftTSC/5mkgQD9Qkqpi4pl0nTV7f+eJUc+a7aXxtU9QzvLRy/x73+IaWn/ROfpmPVtzfkmaxfjeQ9btdmV/9+dou+59eu/+r/BivavlJP8Uz/BEvb/Gdu13feXSBP9DymdpKbdNfrD1tGNFnUvF7VFzfmZjmbZZaRHVJv6ZrwPJ/XnX58S69Wj5be9Auqimd2rJVPtNSr1ObfxoasWZgElbVUw1L+I769btLIKJTeeq1Fpu2vQd4JZor+orrMwPOiLhgtvRjjw6sXfW6hUVO1WskT9RlshZT9QjtOA1raywSVVRr/LZgSJiqyICACZFhNWR1Y5GoFgs09QbOlwBMIXGngAmhiKrpIRQRbeE9QclrQBUEDAHoEyU7Az3JWpkuKhgR87WK3jOpsTJiIfnFTwwla1Xm6q/TaD3Rv/BzPMc9HvFFPl3qui1KNSjTz6Lwb/ETw8aLEhunjKjOGl1ifSe/DbzF9/nPkpWbPT7U0Fg9PDzjtsrWj/HT35RHfKFCv5UO/cbGGjSU1jANw0UVVM749CSmaRzDRVUdjxETek40FzOV2kRV57Tg02GUXOwLtaO7j3q4TAaEDBQee1aGCXfkSqwRanKq+uu3PaaO57Wbo+UxVX6qUQL25tAsKk/zH9fLWZLxxaea4wnN7z9jrPU6bMY4s+Au4tjiYzoMR7eo5kVsOhlkZrIuUpcIQt1oj6n0F6JkLx+HDJnSl3xVA2gXlf4K7k5B1BYxlfYataO9QceIiGCZUtBRZusV1k+NxKvqRndMFeIRcMZFshFAxxVZXrAOEBjQH8s9NDvVagV3oGkqCnKrLUcSrjeAZlGlK7jDgj/+ruQF61La2wB6p78gKVgbE3HBGb62ArZB5gQoU18zaM+oz5n7yFDbeXfZwbrHhOGW5dByzbg19Z72PmSsLdPd5yyzCiIgIuK4oNlPtCL0xet60h5oa1q+oHxNcbxWlIIdJdvAvFR7Ra/BouT2Y42oijdm5adA/aWQ0tvXO5Ld2gGLRJWfWe8W+M7CwforrxGf2gWLRJWXWS9TrtdfKT6+LinsLc9O9hxvulEFi0QF/czEws17hjtMVmTXTSbLDiGL+EqaWZfHKlFlB+vBDdvgryW1iKP6y++cE+GLT1XAKlHlBesREyY57RhDLpPH/WXflk1X6kkdaQUsE1V+sO7RybgGDJhu+FOWgKTevQL6Kz/r5Th3mpoxYLwmKp8RZ0T06K0E+evpWVnuqYZlTpUXrEPZ3sVRhiyb7XvTPqwTVX5mvV/qCq5X8jEhF9umPzjDy/QVn7DEDunscFz25VTDOqfKC9bLdbjTd/STS1goqrwymCl9yTbVgX3TX17N+oz3MiZkKhGSdix0quIyGE/OktGNjU6VH6wDzPAk8NaLnaLKq1mHeamLpDL1YuX0V7TBdCrld7qx1anyNpiGjMWndGOrqPK6wYAv9QbasXT6g7wymPmRSF2RlkYsFlV+GUxAT7rBaMRiURV3g+lJtkobFouqqC5BehfrxGpR5WfWpXexTiwUVbpje3aw7q3tnhGUYs9p7+mzJWbLyqkBnYyeCAFjWarRiHVOtXq2RHaw7jPN2VsjqMA6Ua2THax7dBkxldP5tGCNqHpERPQzNmEVlcFIm1kdWCKqHh3yJrui/cl9yVZpwApR+Ulv0SizW0t+Zr3cpi2hKlaIqpv0bOnlxE950pHexVqwo0phlpxHH2Y2qs0rgwlATmbXgR2iOsYjpIeXiGud7DIYKSvWhRXTH8wbxub3v8srg5nIti0NWCOqvIhqTn6wHtGXRRvFWCOqHsOCLLmUwdSJNaK6iaL2HBNZslGKM6LKz6x3pBexYpwRVb4fSWGxYuxIKZRjkLtp66LyewlFOORU+WUwHbn6U4lLTpWdWb/gQ/hMGEm+ShkOOVVxGYzPRMpg1OCSUxV3g7lgLJOgIpxyquI+67JpSxWOiSo/s96X0mJlOCaqeWbd2xCQJ+V66rBni1YZfC4Y0ONXNp4JpMGQQpxyKo9zRkx5gy9sPNdjSleu/5TglKgg4BFvEWeIar5jcMRIwvWdcUxUc1nFfDn3+Z4cF7krbuWpAM74NeA9uc8HshK4M845FcCLuWF5WJBzF8ripKjIzaubvmnLo0ufgPlZhYbi3vQ3J3uKM3/TVsRLzM/6Cg1OgjjqVNn4TI13qkWa1uQN+yKqFUKGTI3+gy1WLw1eAXB1+stjzAu8lzE9Y6urFifTX4CpYxRRZXG8cYS3OYzxCImArrHxn0x/mZi6GSKApCbMMzjzL6LKwNzexXErKupFVBmYW64XExNzmdEx0CREVBuEBh/fFqe+RqauUu7FwF7TozCNDh/hGf616WFksJ7qeIOv899ND2oDcapMOoZGVWmnMrdtt4gqE1NbdlxPfYbmqJajFNYIjF2sjbk01p/SoxRRbdA3NQRuQToBRFQZyM6aHZFlmk0ijpseQssRp9rA1CC9NYhTrRMsu7ILWyKiWqcjkx/s6NUy/a0gQfqcgO/scMKYiGoFUzPp9fNB3uT7zHipumvJ2t8qE17YMqIyNV26PQEeEPMf/FLF34k4VQpzM+n18ylifrDlFCiiSmFuJr1uIka8uPV3i6iWSJCuCEkpXCOZdFWIUy3Jy6R7Bu4CNLqVpFud9IrIz6T/LReGbIYKl/IOmTFoeji5iFMl5AXpX85pkdYkvtlOJaKakxek/zWPiCUhWg2Z/uZkB+l/zG/yFhjRX8WjQ8AAj4BjI0aUjzgVkB2kR/wDn+HbxEYE6p3lKEOGyWMeXRMnQkkpQF6Qfswv8jY+Z0iQfoyXnGbvL0uKu5yZWKYjooKicpeIAXeaHh4AM8JE3sHyNsQzRPIrSEw1z0NlCydizMyYor15i1uPkGeXj82aHlQW4lRFmXSzcuwhENLjhURKF8wYGykrCdRzM+m+UWuBHhNWKzJDU1Md4lT5mfSOUTnreUR1sfKIoYkFiamKgvRnK72TTsJk8luE5R6emUE64Pz0l1/uEtFvenAFBPRMzFAlOC6q/Jr0oREpz3ZiuahuOiGhHUF6y7A9UC8+IaEtQXrbsN6pik5IyK9JN7frZxuwXFRFJyS0NUg3Htunv6ITEtqSSW8f1jtVSJ+QLiO6a1OaBOm6sF5UeSckHPCNnNebnAFqBw6IKouAR3wl5zkJ0nfE/pgqi4Axb+efMp+LDF35bxEuiipgxLtyn5UgfXecm/4CLpNe5FkLNBKkK8A1pwoYFUZMkklXgVNOde1SeU4lQboC3HIqf2X9b7PzkgTpanDKqQB8vsL/Jee7rCPlLmpwTlTzTHqH0YaoJEhXhXOium7BuB49SSZdFc6JSspd9OOYqKTcpQbcuvqTcpd6cMyptit3kVirEm7t+6tWk96WZojm4ZRTbRekm94M0TwcEpUE6TXh0vQXFXahWqVNzRDNwyGnqhKkbzZDjBBK4k5KoVqQvtkMURKjpXFn+qvW3WVGtNIMMTC2w4qJODL9VQ/S5yuBHlM8Qjm0rQquTH/VM+npZohn+OJUFXDEqapm0tebIQ7xtz5R2DnccKrq3V2ymiGa0aO4DTjhVFUz6SH9zC4xQjkcEJVk0mvGhelPyl3qxgGnku4uNWO/U0kLxvqx3qmkJr1+LBeVBOkNYPv0J0F6E1juVBKkN4DdTiVBejNY7VQSpDeDxaKSIL0hbJ7+8kPxUIJ0rVjsVHlBuqAZe50qMOZAbeewV1QdmeIaQ+P0d7vBj+VJHqpBNIrqdO3+gxo/VsfUk9CdQJuo9jlfuX9S6xWBBOkNoi+m+jR/krq3n/qvfiRIbxZtTrXuUxBzqO0jrH7lZ9KFGlCxQ3mfz/MBPse7+QU+zysA3ObfUq845KugQb3zd9wjZi+5v7d8XNb2mkOBU90FHnIbOOR+8tjJyrXfCQAPuKt46NdC3by1OKtrOiqc6g/Z5228CnyAdySPPc2ry+ePeMAh8BTvUzbuGNjj2qWybxfOJdSKClF9kyNeAyDgHwE4SHV02uf9fASA31A47j0olBMiqQZRMkuccgLs8zi5vjtNXeedpF6lJlMVF0x6m7cyCdaPElE95IRDzjlI7i/Snoec8zi55rvLAx5vJES3G3I5OS0+nsiqbhT8zvd5SHpR5ogjrQOuditeVT8KRHW05j8q3Ch/uCIr49k9o37Ih/nRVFpTf9Z8L5WfKnMr1I3ykOO21uqEqi61+HgirxpRv/b3MX5e26aCuLJL7Ylf1c9eDHz4Uz/yW6re8O9+8sEPw3u++xOzW/+jerDHz6TyU+Vvr1OlQj3ExPzNaayMcLmsG8T9+FLdG8dxtfxU+lamv1p5EuDbb1fXJfVy+X9n3OH3eZ7n+XF14y2XRU/fikfVzpMArz/1urI3vFy59ybf4u95l7rxVpHT/DY9DQr10MLpb5vJT6a/GmljoL5dqC5uVR/K/xH36WpMKWydp3r+OxrHJazSspmhuqySpZ04nsbd+NbL0si6Blonqq1iqigeJnFZXzxLPy0U1RaLNaO4E/txN57GsXiWflomquqyiom/9L9xHMeXcS/2xbPqoIWiqiYriOk8d75Ib2jzrH1OOeWUQ+CQU04Vb/NoE0aIap/XKuwJrDb5Le4Ft17+3YeTWKtnPVzudryvaY9jOzBCVPBFYq54rWTZTPVQfY6n2bPu8hiAk9r2YhvJ3g2LYzqrOFf52WR71w/4cz5x46sXwy5Oe2anPINbv/fr0cffEQAzBgwIiJJztwd86d/v/RkDZlt+iiv+ivfxMb5Z2+/NQG4SVV0m7nHIE8Ab3OMTJf8k6Z19WXuUiz6YR/Tciwc/E+EBY44ZE9HBBy4Y8Mpfvv7KVmeRnhLw8WSXtrsYMv19q8Lktzr47G1Z5T6U+jjrLvFyT5G7GCGq2zs54maDjiqojLOOOOD+WmsSFzFCVM2jwrMOOQIOuXI7TKcGUZ3s/hY1sZtnHSx3Oz5u0WfWgyZRHSapwFMetiwJuI1nHXLOovvWEQ8V7cRuL5qd6nYr/9XqzmfZjsRUudSRg7cTDaLaX059pxy1PGgVz9oGTU61yNYcLtsLtRnxrGpoEdUBV8v/v2plVLWJeFZ5tIjqZNn787ZlGWbxrDJoEdX9xJ32OW9ZQqEM4lk3oUVUV5xzyulaj2K7EM/KR4OoDpzJUohnZaNBVCe1Hm3UPOJZ62gQ1X0HFynEs9IoFtUBp1xx7miFtnjWHFmmUYx4lohKE257lohKG+56lohKMy56loiqBlzzLBFVbbjjWSKqWrnZs4Kmh6gAEVUDFHjWJz/b9OAUIKJqiEzP+vrVZ9pf0iiiapg1z/rod5sekBJEVI2T8qw/taOgUURlCMGtl+/8c9ODUMP/A1ogL1Nbh0MxAAAAJXRFWHRkYXRlOmNyZWF0ZQAyMDI1LTAxLTE1VDEzOjMwOjIzKzAwOjAwp+LxZgAAACV0RVh0ZGF0ZTptb2RpZnkAMjAyNS0wMS0xNVQxMzozMDoyMyswMDowMNa/SdoAAAAodEVYdGRhdGU6dGltZXN0YW1wADIwMjUtMDEtMTVUMTM6MzA6MjMrMDA6MDCBqmgFAAAAE3RFWHRwZGY6VmVyc2lvbgBQREYtMS41UzZawQAAAEp0RVh0c2lnbmF0dXJlADBkMjZlMTkwZGJhYTYxNjU1MmM5ZGI3NWRlMWFlODc0ZmYzM2M0NDgzOGQ0YjI4YjE1Njg2ZTc3NzBiNjE5NTUKh85BAAAAAElFTkSuQmCC", "text/plain": [ "" ] @@ -435,7 +269,7 @@ }, { "cell_type": "code", - "execution_count": 26, + "execution_count": 5, "id": "667e3870-42c0-4a77-955a-4f5c14fe1e8c", "metadata": {}, "outputs": [ @@ -443,13 +277,13 @@ "name": "stderr", "output_type": "stream", "text": [ - "/home/apn/.cache/pypoetry/virtualenvs/pyfeyn2-rvdu2REa-py3.11/lib/python3.11/site-packages/feynml/pdgid.py:111: UserWarning: Type remnant is not consistent with pdgid -2, which is anti fermion. Using remnant now.\n", + "/home/apn/.local/lib/python3.13/site-packages/feynml/pdgid.py:111: UserWarning: Type remnant is not consistent with pdgid -2, which is anti fermion. Using remnant now.\n", " warnings.warn(\n" ] }, { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAASoAAAEWCAQAAAB7ixJyAAAAIGNIUk0AAHomAACAhAAA+gAAAIDoAAB1MAAA6mAAADqYAAAXcJy6UTwAAAACYktHRAD/h4/MvwAAAAlwSFlzAAAAZAAAAGQAD5bF3QAAAAd0SU1FB+gBHg0VJyEC+n8AABJaSURBVHja7Z1fiGTZXcc/HTZmCVFz8SVhJGxfo8aH+LB3wT+4jNG7ELMEUbn7EAaRPNRAQFgjdJs8iCQotQ9OkBChypeVWSFUo8GQaQLVaDagjVBNguCAZKuTMNkkYqYrRCdOZnf7+nBvVd2qure6btU55557zu/TMNXV1V19qus73/O7v/M7v0NKitA4IV3G9IibHogSRFSNkzBgzCFB0wNRhoiqQdzypzkiqoZwz5/miKiM46o/zRFRGcVlf5ojojKE+/40R0RVoDv7YwyIFD6vH/40R0RVIGIIQMhY0TP65E8zHmt6AFYRcwZAl5sKni0hIaLPU0yafmGGEacqMCAGOnQBCLf2Fy/9qYCIqsAFcEgHgC4xAYe1n8O3+KkMEdWMiBEdwvzzITFxrXDdd3+aI6IqJc6nwE39RvypiIiqlIAuMQnJld8p/rSKiGoHxJ9K2UuBvaZH0T5COiSccMTJxj8TEAEx7+Bx4Pt8zNVUg+SpriZLiE44z++/h5/nx7bIP02Av+EJAC7puiopZPrbgAHp0sentnymiB7fJ+UBzzf9orQiorqSZElSnS2eI6DDiCEdvszr/EHTL0kzIqorCZnMBHWxxVVeRI8xXUIg4L95X9MvSDsiqrVk13dfmkmqXu3C3J+mhEqrH2xFRFVBMf8UbiGpoj/5hoiqhNX806CGpFb9yTdEVAtU5cc3XQX02Z/miKhm7JYfF3+aYSr5WXyrTmpkoc0wzY/3txxZRIeYI56bJUj9xnunEn9Sj0ZRBXTpEFv7B9+1vkDipyo0imqQy6l+9aR+xJ90ok1UERfExNYVhog/6UebqOJ8dd8mxJ/MoNGppqKyoSZS/Mkg+lIKZ5yRMCHM99I1x2777wISOkw4UrIX0A+0phTChuMp8adm0Jynihp7S5qIn2y8zm0CzaKKSUkZbLArRR2m/Smmm3+MJIgHDIkqJTW066TJ67tQnCoj200TapuinuSFhft3uMPXtPymt/HrfLDm/pYi0/W7vqzf7Up29bd9K4qreGLp/m/yU7ys5Tf9HE/yZ7y4xU/udn2X9Vv4Ey2vqa0Ym/5S7WnDhB6jmj+j4vrOxjRvsxgR1QU9Q9eAm7+96vLjh3nfBWGKdlGN6RjLVm3aAU9t/mlgxZqBTThVTzXYwHfUr99dAAmd2lOvs7i07T0ioH/ld4WK6zMjzkg4ZzLz44AOLF31+oVDTtVtJE90yGgppsqaEPmMQ6JKjMZvUwbEhYqMrB9fYlkNmWkcEtV0gcZs4HwBwBhyd4oYEYuomh7COpItvCdS1gV9k981AKBHku8MDCRrZbuoYEjKl2t6z8hgZcRU8tPfGEvWapOrv06j9UT/wS/xNHd4yOf4xEbXbUmhQZl+poV/098YN16U2DibiOqs0SXWt/H7wBv8kP/asHKzyzMNjTUgILBuq6x5rJ/+xjzkszX6rXToNTbWqKG0hm1YLqqodsanKzFN41guqvoEDBnR9aK5mK0YE5XJaSGkwzC/2BeMo7uPejxLBsT0N1iZU8mIm3Il1giGnMp8/XbA2PO8dnO0PKaqTjVKwN4cmkUVaH5zg4olmVB8qjnepPn5J5xovQ6bcFJacJdw5PAxHZajW1RZEZtO+qWZrPPCJYJgGu0xlf5ClPLl45gBY3qSr2oA7aLSX8HdWRO1JYylvYZxtDfoGJIQzVIKOspsg7X1U0PxKtPojqliAiLOOM83Aui4IqsK1gEiC/pj+Ydmp1qs4I40TUVRZbXlUML1BtAsqmIFd7zmzd+VqmBdSnsbQO/0F+UFaycknHNGqK2ArV86AcrU1wzaM+oZmY8MtJ13Vx6sB4wYNFwO7SNG1/4GGvfl9SqC9YguY1kHNErLF5TnrI/XpImGSZwR1fqNWZICNYjuIj2TdAgrOtodEpZ2yQvzKbPPeT5BnnCSrxrad4Bce3DIqaoy6+vK9WLS/LGAdJbv78pkuQsutRLKMuurJcvRmjKYMyg8Nt8YmrmU3Udf2otDTlUVrEdro600T3J0uMj9SXqh74j+eiqTnBGUyOeMFxgwqmzHGJBVfZ0RAKG0vN4Vt0RVlVk/4ime4ZxhiawyCWWlhBEQyzS3K66J6qiyuGZCn5NSUYXEHOX3pL2GAlwT1boymLDChUJCJmRBeyCiUoBTgTqsy6yXL+QMZzWpXS5kB44SnBNVVWa9qsf6PCclHVtU4aCoymvWq7adhiWfCbvhoKiqymDG9DxvRW0I1wJ1qArWJ/wMJ8RSBmMAB51qfRlMIGUwunHRqaoy6xkTAtkQrxc3RVWVWQeIQXJRmnFy+lu3wXQsm0t146pTVWXWY07Ep3TjqqiqusFAKFlz7Tg6/UFVZj07EulQpKURh0VV3Q0m27YlGXRdOCyq9d1gpA5dHw6LqnqDqfQu1ovToqrOrEvvYp04KKpiAUt5sB4wEp/ShztbtIpnS0xmG7X6dEo2mEacyFKNRpxzqsWzJcqD9ZBx5d4aYXecE9Uy5cF6wCFDxlLrqQVnRNUlIaG34j/ry2CkzawOHBFVlw5Vk926/ck9yVZpwAlRhXlv0aT0qNvqzHq4Nj0qbIsTojrMD4zsVsRPVdKR3sVacKNKYZJvXi/fLFpVBhOBbHHXgRuiOiIgpktQ0VyjvAxGyop14cT0B9m58L3KR6vKYEaybUsDzoiqKqLKqA7WE3qyaKMYZ0TVZbAmSy5lMCZxRlRX0VvjYyNZslGKN6Kqzqx31sRiwjZ4I6pqP5LCYsW4kVLYjH7lpi3p8qkWj5yqugxG34k5XuKTU5Vn1s95hpARQ8lXKcMjp1pfBhMykjIYNfjkVOu7wZxzIpOgIrxyqnWZddm0pQ7PRFWdWe9JabEyPBNVllkPVgQUSLmeOtzZorUJIef06fIbK49E0mBIIV45VcBdhox5jc+uPNZlzKFc/ynBK1FBxEPeIC0RVbZjcMiwNeG6vVv3PRNVJquUL1Q+3m3ReX+WxoF+5akAzvgt4F2Vj0etWQm09hA5/0QF/8wf8VrFY9mR3O3A3mPkvJv+Mqry6u3YtJXtxra3tNBTUZWzbm3QHqZRn6URFSKqBULG1jvVdDd2bG+9qo8xVTXn3KTD2OoOC0ke89kbUYlTlZIwsra6arokPiK0dYziVGUclR7hbQfZbuwOZ/nJzzYiTlWKrZshosK/do4QRFSl2LwA0oaKehFVCfaW66WkpFyUdAy0CRHVCrHFx7elhY+hrauUeymw1/QobKPDB3iS/2x6GCUspzpe4yv8oOlBrSBOVUrH0qiq6FT2tu0WUZVi67rafOqzNEc1G6WwhL1rgCkX1vpTcZQiqhV6tobALUgngIiqBNlZsyOyTLNKwlHTQ2g54lQr2BqktwZxqmWiWVd2YUtEVMt0ZPJj18Vqmf4WkCA9I+LbO5wwJqJawNZMunnex+v8kAkv1HctWftbZMRzW0ZUtqZLtyciAFK+y6/V/JuIUxWwN5Nuno+T8qMtp0ARVQF7M+mmSRjy/NY/LaKaIUG6IiSlMEcy6aoQp5pRlUkPLNwFaHUrSb866a2jOpP+j5xb0rQjnsk7ZkK/6eFUIk6VUxWkf6GiRVqThHY7lYgqoypI/3sekkpCtB4y/WWUB+l/zu/yBljRtSCgQ0SfgIgjK0ZUjTgVUB6kJ/wrn+RbpFYE6p3ZKGMG+dcCDm2cCCWlAFVB+hG/ypv5tCVB+hFBfpp9OCspPuTMxjIdERWsK3dJ6HOz6eEBMJn1+IxmtzGBJZJfQGKqLA9VLpyEEybWFO1lLW4DYp6afc3Kvi/iVOsy6Xbl2GMgpstzuZTOmXBipawkUK/MpIdWrQUGjFisyIxtTXWIU1Vn0jtW5ayziOp84SuWJhYkploXpD9V65l0EueT3zQsDwjsDNIB76e/6nKXxN7uv0BE18YMVY7noqquSR9YkfJsJ46LKuCQHhEQlpYKtyNIbxmuB+oJLwAB5WFtW4L0tuG8U02jprKG+9U16fZ2/WwDjotqHjWtTmdtDdKtx/Xpb5ouDDmHpe5zbcmktw/nnSqmR8whQw6XpjQJ0nXhvKggICSLrorc4GsV329zBqgdeCCqMiIe8sWKxyRI3xH3Y6oyIk54C18tfSyxdOW/Rfgoqoghb698VIL03fFu+ou4yHuRly3QSJCuAN+cKmK4NmKSTLoKvHKquUtVOZUE6Qrwy6nChfW/1c5LEqSrwSunAgj5Iv+Xn++yjJS7qME7UWWZ9A7DFVFJkK4K70Q1b8G4HD1JJl0V3olKyl3045mopNzFAH5d/Um5ixk8c6rtyl0k1qqFX/v+6tWkt6UZon145VTbBem2N0O0D49EJUG6IXya/pK1XagWaVMzRPvwyKnqBOmrzRAThA3xJ6VQL0hfbYYoidGN8Wf6q9fdZUKy0AwxsrbDio14Mv3VD9KzlcCAMQGxHNpWB1+mv/qZ9GIzxDNCcaoaeOJUdTPpy80QB4RbnyjsHX44Vf3uLmXNEO3oUdwGvHCqupn0mF5plxhhMzwQlWTSDePD9CflLqbxwKmku4th3HcqacFoHuedSmrSzeO4qCRIbwDXpz8J0pvAcaeSIL0B3HYqCdKbwWmnkiC9GRwWlQTpDeHy9FcdiscSpGvFYaeqCtIFzbjrVJE1B2p7h7ui6sgU1xgap7/rDb6sQPJQDaJRVKdL9+8bfFkdW09C9wJtotrn7sL9Y6NXBBKkN4i+mOoT/EXh3n7hX/1IkN4s2pxq2acg5UDbS1j8qM6kCwZQsUN5n8/wXj7NO/gVPsNLAFznG4XvOOBLoEG92TPukbKX39+bfV3W9ppDgVPdBh5wHTjgXv6144Vrv2MA7nNb8dDnQl29dTirazsqnOpP2efNvAy8l7fmX3uCl2eP3+I+B8DjvFvZuFNgj7lLld9OnUswigpRfZ1bvAJAxL8BcKPQ0Wmf9/ABAH5H4bj3YK2cEEk1iJJZ4pRjYJ9H+fXdaeE677jwXWoyVemaSW/1ViZB8ygR1QOOOeAuN/L707TnAXd5lF/z3eY+j1YSotsNeTM5TV+eyMo0Cv7m+zyguChzi1taB1zvVrzKPApEdWvJf1S4UfVwRVbWs3tG/YD385OFtKb+rPleIT+1ya1gGuUhx3Wt1Ql1XWr68kReBlG/9vchflnbpoK0tkvtiV+ZZy8F3v/xn/g9VU/4Tz97/8fhXd/76cm1/1E92KMnC/mpzW/nqVLBDCkp/3CaKiOeLetGaS+9UPfEaVovP1W8lenPKI8BfOst6rqkXsw+O+Mmf8yzPMs71Y13syx68VY8yjiPAbz6+KvKnvBi4d7rfJN/4e3qxltHTtltcRoUzNDC6W+byU+mP4O0MVDfLlQXtzKH8v/EPQ41phS2zlM9+22N4xIWadnMUF9W+dJOmo7Tw/Tai9LI2gCtE9VWMVWSDvK4rCeepZ8WimqLxZph2knD9DAdp6l4ln5aJqr6skpJP/+/aZqmF2k3DcWzTNBCUdWTFaR0nr47TW9o86x9TjnllAPggFNOFW/zaBNWiGqfV2rsCaw3+U3vRdde/MiDUarVsx7Mdjve07THsR1YISr4HCmXvLJh2Uz9UD0j0OxZt3kEwLGxvdhWsnfF4pjOKs5FfjHf3vUj/pqPXvnd02GvT3uWpzyja3/428mH3xoBE/r0iUjyc7f7fP47d/6KPpMtX8Ulf8e7+RBfN/Z3s5CrRGXKxAMOeBPwGnf46IZvSXFnX9ke5XUvLCB5+vkbv5AQACcccUJChxA4p89Lf/vqS1udRXpKxIfzXdr+Ysn0980ak9/i4Mu3ZW32otTHWbdJZ3uK/MUKUV3fyRFXG3TUQWWcdYsb3FtqTeIjVoiqeVR41gG3gAMu/Q7TMSCq492fwhC7edaN2W7HRy16zXrQJKqDPBV4yoOWJQG38awD7jLtvnWLB4p2YrcXzU51vZX/a3Xns1xHYqpKTOTg3USDqPZnU98pt1oetIpnbYMmp5pmaw5m7YXajHhWPbSI6gaXs88vWxlVrSKetTlaRHU86/153bEMs3jWJmgR1b3cnfa527KEwiaIZ12FFlFdcpdTTpd6FLuFeFY1GkR1w5sshXhWORpEdWz0aKPmEc9aRoOo7nm4SCGeVUSxqG5wyiV3Pa3QFs/KkGUaxYhniag04bdniai04a9niag046NniagM4JtniaiM4Y9niaiMcrVnRU0PUQEiqgZY41kf+1TTg1OAiKohSj3rK5efbH9Jo4iqYZY864Pfa3pAShBRNU7Bs/7SjYJGEZUlRNdevPnvTQ9CDf8PnOk1vvoBNh0AAAAldEVYdGRhdGU6Y3JlYXRlADIwMjQtMDEtMzBUMTM6MjE6MzkrMDA6MDBYOTNZAAAAJXRFWHRkYXRlOm1vZGlmeQAyMDI0LTAxLTMwVDEzOjIxOjM5KzAwOjAwKWSL5QAAACh0RVh0ZGF0ZTp0aW1lc3RhbXAAMjAyNC0wMS0zMFQxMzoyMTozOSswMDowMH5xqjoAAAATdEVYdHBkZjpWZXJzaW9uAFBERi0xLjVTNlrBAAAASnRFWHRzaWduYXR1cmUAMDYyZWRiYWUxYWE1YmMxZmZlYjk4Mjg0NTMzMDJkODFlOTBmNjAwNTQxMjA0NTU1NmU5MmJjMmYyNmE5YWVlYs8DDK8AAAAASUVORK5CYII=", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAASoAAAEWCAQAAAB7ixJyAAAAIGNIUk0AAHomAACAhAAA+gAAAIDoAAB1MAAA6mAAADqYAAAXcJy6UTwAAAACYktHRAD/h4/MvwAAAAlwSFlzAAAAZAAAAGQAD5bF3QAAAAd0SU1FB+kBDw0eGrlnjPoAABJaSURBVHja7Z1fiGTZXcc/HTZmCVFz8SVhJGxfo8aH+LB3wT+4jNG7ELMEUbn7EAaRPNRAQFgjdJs8iCQotQ9OkBChypeVWSFUo8GQaQLVaDagjVBNguCAZKuTMNkkYqYrRCdOZnf7+nBvVd2qure6btU55557zu/TMNXV1V19qus73/O7v/M7v0NKitA4IV3G9IibHogSRFSNkzBgzCFB0wNRhoiqQdzypzkiqoZwz5/miKiM46o/zRFRGcVlf5ojojKE+/40R0RVoDv7YwyIFD6vH/40R0RVIGIIQMhY0TP65E8zHmt6AFYRcwZAl5sKni0hIaLPU0yafmGGEacqMCAGOnQBCLf2Fy/9qYCIqsAFcEgHgC4xAYe1n8O3+KkMEdWMiBEdwvzzITFxrXDdd3+aI6IqJc6nwE39RvypiIiqlIAuMQnJld8p/rSKiGoHxJ9K2UuBvaZH0T5COiSccMTJxj8TEAEx7+Bx4Pt8zNVUg+SpriZLiE44z++/h5/nx7bIP02Av+EJAC7puiopZPrbgAHp0sentnymiB7fJ+UBzzf9orQiorqSZElSnS2eI6DDiCEdvszr/EHTL0kzIqorCZnMBHWxxVVeRI8xXUIg4L95X9MvSDsiqrVk13dfmkmqXu3C3J+mhEqrH2xFRFVBMf8UbiGpoj/5hoiqhNX806CGpFb9yTdEVAtU5cc3XQX02Z/miKhm7JYfF3+aYSr5WXyrTmpkoc0wzY/3txxZRIeYI56bJUj9xnunEn9Sj0ZRBXTpEFv7B9+1vkDipyo0imqQy6l+9aR+xJ90ok1UERfExNYVhog/6UebqOJ8dd8mxJ/MoNGppqKyoSZS/Mkg+lIKZ5yRMCHM99I1x2777wISOkw4UrIX0A+0phTChuMp8adm0Jynihp7S5qIn2y8zm0CzaKKSUkZbLArRR2m/Smmm3+MJIgHDIkqJTW066TJ67tQnCoj200TapuinuSFhft3uMPXtPymt/HrfLDm/pYi0/W7vqzf7Up29bd9K4qreGLp/m/yU7ys5Tf9HE/yZ7y4xU/udn2X9Vv4Ey2vqa0Ym/5S7WnDhB6jmj+j4vrOxjRvsxgR1QU9Q9eAm7+96vLjh3nfBWGKdlGN6RjLVm3aAU9t/mlgxZqBTThVTzXYwHfUr99dAAmd2lOvs7i07T0ioH/ld4WK6zMjzkg4ZzLz44AOLF31+oVDTtVtJE90yGgppsqaEPmMQ6JKjMZvUwbEhYqMrB9fYlkNmWkcEtV0gcZs4HwBwBhyd4oYEYuomh7COpItvCdS1gV9k981AKBHku8MDCRrZbuoYEjKl2t6z8hgZcRU8tPfGEvWapOrv06j9UT/wS/xNHd4yOf4xEbXbUmhQZl+poV/098YN16U2DibiOqs0SXWt/H7wBv8kP/asHKzyzMNjTUgILBuq6x5rJ/+xjzkszX6rXToNTbWqKG0hm1YLqqodsanKzFN41guqvoEDBnR9aK5mK0YE5XJaSGkwzC/2BeMo7uPejxLBsT0N1iZU8mIm3Il1giGnMp8/XbA2PO8dnO0PKaqTjVKwN4cmkUVaH5zg4olmVB8qjnepPn5J5xovQ6bcFJacJdw5PAxHZajW1RZEZtO+qWZrPPCJYJgGu0xlf5ClPLl45gBY3qSr2oA7aLSX8HdWRO1JYylvYZxtDfoGJIQzVIKOspsg7X1U0PxKtPojqliAiLOOM83Aui4IqsK1gEiC/pj+Ydmp1qs4I40TUVRZbXlUML1BtAsqmIFd7zmzd+VqmBdSnsbQO/0F+UFaycknHNGqK2ArV86AcrU1wzaM+oZmY8MtJ13Vx6sB4wYNFwO7SNG1/4GGvfl9SqC9YguY1kHNErLF5TnrI/XpImGSZwR1fqNWZICNYjuIj2TdAgrOtodEpZ2yQvzKbPPeT5BnnCSrxrad4Bce3DIqaoy6+vK9WLS/LGAdJbv78pkuQsutRLKMuurJcvRmjKYMyg8Nt8YmrmU3Udf2otDTlUVrEdro600T3J0uMj9SXqh74j+eiqTnBGUyOeMFxgwqmzHGJBVfZ0RAKG0vN4Vt0RVlVk/4ime4ZxhiawyCWWlhBEQyzS3K66J6qiyuGZCn5NSUYXEHOX3pL2GAlwT1boymLDChUJCJmRBeyCiUoBTgTqsy6yXL+QMZzWpXS5kB44SnBNVVWa9qsf6PCclHVtU4aCoymvWq7adhiWfCbvhoKiqymDG9DxvRW0I1wJ1qArWJ/wMJ8RSBmMAB51qfRlMIGUwunHRqaoy6xkTAtkQrxc3RVWVWQeIQXJRmnFy+lu3wXQsm0t146pTVWXWY07Ep3TjqqiqusFAKFlz7Tg6/UFVZj07EulQpKURh0VV3Q0m27YlGXRdOCyq9d1gpA5dHw6LqnqDqfQu1ovToqrOrEvvYp04KKpiAUt5sB4wEp/ShztbtIpnS0xmG7X6dEo2mEacyFKNRpxzqsWzJcqD9ZBx5d4aYXecE9Uy5cF6wCFDxlLrqQVnRNUlIaG34j/ry2CkzawOHBFVlw5Vk926/ck9yVZpwAlRhXlv0aT0qNvqzHq4Nj0qbIsTojrMD4zsVsRPVdKR3sVacKNKYZJvXi/fLFpVBhOBbHHXgRuiOiIgpktQ0VyjvAxGyop14cT0B9m58L3KR6vKYEaybUsDzoiqKqLKqA7WE3qyaKMYZ0TVZbAmSy5lMCZxRlRX0VvjYyNZslGKN6Kqzqx31sRiwjZ4I6pqP5LCYsW4kVLYjH7lpi3p8qkWj5yqugxG34k5XuKTU5Vn1s95hpARQ8lXKcMjp1pfBhMykjIYNfjkVOu7wZxzIpOgIrxyqnWZddm0pQ7PRFWdWe9JabEyPBNVllkPVgQUSLmeOtzZorUJIef06fIbK49E0mBIIV45VcBdhox5jc+uPNZlzKFc/ynBK1FBxEPeIC0RVbZjcMiwNeG6vVv3PRNVJquUL1Q+3m3ReX+WxoF+5akAzvgt4F2Vj0etWQm09hA5/0QF/8wf8VrFY9mR3O3A3mPkvJv+Mqry6u3YtJXtxra3tNBTUZWzbm3QHqZRn6URFSKqBULG1jvVdDd2bG+9qo8xVTXn3KTD2OoOC0ke89kbUYlTlZIwsra6arokPiK0dYziVGUclR7hbQfZbuwOZ/nJzzYiTlWKrZshosK/do4QRFSl2LwA0oaKehFVCfaW66WkpFyUdAy0CRHVCrHFx7elhY+hrauUeymw1/QobKPDB3iS/2x6GCUspzpe4yv8oOlBrSBOVUrH0qiq6FT2tu0WUZVi67rafOqzNEc1G6WwhL1rgCkX1vpTcZQiqhV6tobALUgngIiqBNlZsyOyTLNKwlHTQ2g54lQr2BqktwZxqmWiWVd2YUtEVMt0ZPJj18Vqmf4WkCA9I+LbO5wwJqJawNZMunnex+v8kAkv1HctWftbZMRzW0ZUtqZLtyciAFK+y6/V/JuIUxWwN5Nuno+T8qMtp0ARVQF7M+mmSRjy/NY/LaKaIUG6IiSlMEcy6aoQp5pRlUkPLNwFaHUrSb866a2jOpP+j5xb0rQjnsk7ZkK/6eFUIk6VUxWkf6GiRVqThHY7lYgqoypI/3sekkpCtB4y/WWUB+l/zu/yBljRtSCgQ0SfgIgjK0ZUjTgVUB6kJ/wrn+RbpFYE6p3ZKGMG+dcCDm2cCCWlAFVB+hG/ypv5tCVB+hFBfpp9OCspPuTMxjIdERWsK3dJ6HOz6eEBMJn1+IxmtzGBJZJfQGKqLA9VLpyEEybWFO1lLW4DYp6afc3Kvi/iVOsy6Xbl2GMgpstzuZTOmXBipawkUK/MpIdWrQUGjFisyIxtTXWIU1Vn0jtW5ayziOp84SuWJhYkploXpD9V65l0EueT3zQsDwjsDNIB76e/6nKXxN7uv0BE18YMVY7noqquSR9YkfJsJ46LKuCQHhEQlpYKtyNIbxmuB+oJLwAB5WFtW4L0tuG8U02jprKG+9U16fZ2/WwDjotqHjWtTmdtDdKtx/Xpb5ouDDmHpe5zbcmktw/nnSqmR8whQw6XpjQJ0nXhvKggICSLrorc4GsV329zBqgdeCCqMiIe8sWKxyRI3xH3Y6oyIk54C18tfSyxdOW/Rfgoqoghb698VIL03fFu+ou4yHuRly3QSJCuAN+cKmK4NmKSTLoKvHKquUtVOZUE6Qrwy6nChfW/1c5LEqSrwSunAgj5Iv+Xn++yjJS7qME7UWWZ9A7DFVFJkK4K70Q1b8G4HD1JJl0V3olKyl3045mopNzFAH5d/Um5ixk8c6rtyl0k1qqFX/v+6tWkt6UZon145VTbBem2N0O0D49EJUG6IXya/pK1XagWaVMzRPvwyKnqBOmrzRAThA3xJ6VQL0hfbYYoidGN8Wf6q9fdZUKy0AwxsrbDio14Mv3VD9KzlcCAMQGxHNpWB1+mv/qZ9GIzxDNCcaoaeOJUdTPpy80QB4RbnyjsHX44Vf3uLmXNEO3oUdwGvHCqupn0mF5plxhhMzwQlWTSDePD9CflLqbxwKmku4th3HcqacFoHuedSmrSzeO4qCRIbwDXpz8J0pvAcaeSIL0B3HYqCdKbwWmnkiC9GRwWlQTpDeHy9FcdiscSpGvFYaeqCtIFzbjrVJE1B2p7h7ui6sgU1xgap7/rDb6sQPJQDaJRVKdL9+8bfFkdW09C9wJtotrn7sL9Y6NXBBKkN4i+mOoT/EXh3n7hX/1IkN4s2pxq2acg5UDbS1j8qM6kCwZQsUN5n8/wXj7NO/gVPsNLAFznG4XvOOBLoEG92TPukbKX39+bfV3W9ppDgVPdBh5wHTjgXv6144Vrv2MA7nNb8dDnQl29dTirazsqnOpP2efNvAy8l7fmX3uCl2eP3+I+B8DjvFvZuFNgj7lLld9OnUswigpRfZ1bvAJAxL8BcKPQ0Wmf9/ABAH5H4bj3YK2cEEk1iJJZ4pRjYJ9H+fXdaeE677jwXWoyVemaSW/1ViZB8ygR1QOOOeAuN/L707TnAXd5lF/z3eY+j1YSotsNeTM5TV+eyMo0Cv7m+zyguChzi1taB1zvVrzKPApEdWvJf1S4UfVwRVbWs3tG/YD385OFtKb+rPleIT+1ya1gGuUhx3Wt1Ql1XWr68kReBlG/9vchflnbpoK0tkvtiV+ZZy8F3v/xn/g9VU/4Tz97/8fhXd/76cm1/1E92KMnC/mpzW/nqVLBDCkp/3CaKiOeLetGaS+9UPfEaVovP1W8lenPKI8BfOst6rqkXsw+O+Mmf8yzPMs71Y13syx68VY8yjiPAbz6+KvKnvBi4d7rfJN/4e3qxltHTtltcRoUzNDC6W+byU+mP4O0MVDfLlQXtzKH8v/EPQ41phS2zlM9+22N4xIWadnMUF9W+dJOmo7Tw/Tai9LI2gCtE9VWMVWSDvK4rCeepZ8WimqLxZph2knD9DAdp6l4ln5aJqr6skpJP/+/aZqmF2k3DcWzTNBCUdWTFaR0nr47TW9o86x9TjnllAPggFNOFW/zaBNWiGqfV2rsCaw3+U3vRdde/MiDUarVsx7Mdjve07THsR1YISr4HCmXvLJh2Uz9UD0j0OxZt3kEwLGxvdhWsnfF4pjOKs5FfjHf3vUj/pqPXvnd02GvT3uWpzyja3/428mH3xoBE/r0iUjyc7f7fP47d/6KPpMtX8Ulf8e7+RBfN/Z3s5CrRGXKxAMOeBPwGnf46IZvSXFnX9ke5XUvLCB5+vkbv5AQACcccUJChxA4p89Lf/vqS1udRXpKxIfzXdr+Ysn0980ak9/i4Mu3ZW32otTHWbdJZ3uK/MUKUV3fyRFXG3TUQWWcdYsb3FtqTeIjVoiqeVR41gG3gAMu/Q7TMSCq492fwhC7edaN2W7HRy16zXrQJKqDPBV4yoOWJQG38awD7jLtvnWLB4p2YrcXzU51vZX/a3Xns1xHYqpKTOTg3USDqPZnU98pt1oetIpnbYMmp5pmaw5m7YXajHhWPbSI6gaXs88vWxlVrSKetTlaRHU86/153bEMs3jWJmgR1b3cnfa527KEwiaIZ12FFlFdcpdTTpd6FLuFeFY1GkR1w5sshXhWORpEdWz0aKPmEc9aRoOo7nm4SCGeVUSxqG5wyiV3Pa3QFs/KkGUaxYhniag04bdniai04a9niag046NniagM4JtniaiM4Y9niaiMcrVnRU0PUQEiqgZY41kf+1TTg1OAiKohSj3rK5efbH9Jo4iqYZY864Pfa3pAShBRNU7Bs/7SjYJGEZUlRNdevPnvTQ9CDf8PnOk1vvoBNh0AAAAldEVYdGRhdGU6Y3JlYXRlADIwMjUtMDEtMTVUMTM6MzA6MjYrMDA6MDD12t7BAAAAJXRFWHRkYXRlOm1vZGlmeQAyMDI1LTAxLTE1VDEzOjMwOjI2KzAwOjAwhIdmfQAAACh0RVh0ZGF0ZTp0aW1lc3RhbXAAMjAyNS0wMS0xNVQxMzozMDoyNiswMDowMNOSR6IAAAATdEVYdHBkZjpWZXJzaW9uAFBERi0xLjVTNlrBAAAASnRFWHRzaWduYXR1cmUAMDYyZWRiYWUxYWE1YmMxZmZlYjk4Mjg0NTMzMDJkODFlOTBmNjAwNTQxMjA0NTU1NmU5MmJjMmYyNmE5YWVlYs8DDK8AAAAASUVORK5CYII=", "text/plain": [ "" ] @@ -481,13 +315,13 @@ }, { "cell_type": "code", - "execution_count": 5, + "execution_count": 6, "id": "27b1d0cd-dc52-4d75-970e-336e95c41e32", "metadata": {}, "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlEAAAGuCAMAAACHjqO9AAAAIGNIUk0AAHomAACAhAAA+gAAAIDoAAB1MAAA6mAAADqYAAAXcJy6UTwAAAE4UExURf///wAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA/wAA/wAA/wAA/wAA/wAA/wAA/wAA/wAA/wAAAAAA/wAA/wAA/66urgAA/wAA/2dnZwAA/wAA/wAA/wAA/wAA/wAA/wAA/wAA/wAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA/AAAOAAAAAAAhAACXAADYAAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAAAAAAAAAAAAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAEdHR1VVVf8AAP8AAP8AAP8AAAAAAAAA/////4iIiAD/AP8AACwe9TUAAABidFJOUwAAEVUzmYhm3cy7IncRd5kzIrtmRIhE3VXux83kqKpp1++PsKig9snuqvegmNKAkN7e0uKqiEQRM2a73SJ37lWZzJB2cNfkzfHEEUSIqt27ZjOZVcwi7ndQwMePes2Q75vmN35dUAAAAAFiS0dEAIgFHUgAAAAJcEhZcwAAAGQAAABkAA+Wxd0AAAAHdElNRQfoAR4LIg1uM+8vAAAPQ0lEQVR42u3dCV/bRhqA8Y5syydNm4Bou7vZe0v2PgJNuj223RsIITRput3bdJvv/w1WemUDBmzreMczIz3/X4iTlrHG8RNpkI3y2msAAACoyFgVdbp2NwC/2C4q6plu7PpBYoNsF9VPcxrYfhDTnO3NoAjbRcVpUXaf6iylbAtTQ1Q+CL2oWUXT+QdNuWa7qE5a1NDe3V9URFK+sFVUPByMhqM0p9540rE2+/x4Z5qU1Nbrd94wb951PY3qLBXV6Zq+mcgvu5G1yc8ralJS98ydbXNvx/U0qrN21OsmY3v7pty0iUltbe2aLKpgWSuqZ5LE7tSvVrSYlN3tWpbtn3ZdT6IGW0VFI5NYPrO5tKawd1I7O9ky6i3X06jMVlGTiemOrM582tCktu9sv3739S3X06jM2j4q/Rhbnfl01V7K6pYt23ozXUy5nkR11l8ptmbVMircot4IeO+US4tK4hBNF49xix/hFrV9/vadN1xPopa0qHe+CtHqg16wSW2fp+7ecz2NGkIuqolLcynq/PztnTddz6SqcItandTdQH3jfOabrmdSVbhFrVpHBX7UO98NeH2eFvWt/4Woweuo++Ge3jQBnz2Yrkzq/u7dnXshfs20HfAKKhdsUevOR21t74SYVbhHu7mAiypwPirQrIIWblGXJa15EYasNirkosq8CENWmxJwUeXfcUdWGxByUfOkru+hVp86ICu7Ai/qtoNekXNRZGVN0EWZ6c1vhinxkh5Z2RB2UTe/v6r0yXKyUhZ6UWZ2wYP596lXuxOy0hN8UbPraJj619IgKxUNKEpovTRMVnU1pSjV77shqxooahmyqoaiViKr0ihqPbIqg6IKIquCKKoMslqPokojq5UoqhqyWoaiaiCrW1BUXWS1iKJUkNUFitJDVhmKUtb6rCjKhjZnRVHW1Mwq1EvlU5Rd1bMK9VL5FLUBlbIK9VL5FLUppbMK9FL5FLVRZbIK9FL5FLV5BbMK9FL5FOXIsqzeukwozEvlU5RLt2S1uxveYnwBRTm3kNXW+fn5Tng7pisoyg/zrL4tlwoOeTdFUR5Js/pOfjXzgHdTFOWX+7ML5H/3e6G+JkhRftnNr5D//R/8MNSXminKK1vnu/d3tmeHvDDfwUBRXrl5/im4rCgqBCFlRVHBCCQrigqL/1lRVIC8zoqiQuVrVhQVNA+zoqjw+ZUVRTWEN1lRVJP4kBVFNY7jrCiqmdxlRVEN5iQrimq6TWdFUa2wwawoqj02kxVFtYz1rCiqjWxmRVGtZSkrimo3/awoCrpZURRyWllRFK5QyIqicF29rCgKt6qcFUVhuSpZURTWKJkVRaGI4llRFAorlBVFoZx1WVEUKliRFUWhqtuzoijUciMrikJ9V7OiKCiZZUVR0LS1TVHQRVHQRVHQRVHQRVHQRVHQRVHQRVHQRVHQRVHQRVHQRVHQRVHQRVHQRVHQRVHQRVHQRVHQRVHQRVHQRVHQRVHQFXBR/VEuzn5DUb4IuKgFFOULioKukIuKJr2uGectUZQvQi6qb3qJ6cemM811XE8IJuyiomhgsqjMQIIauJ4PMiEXle2fTNZRLEVNXE8HmaCLiuNsGdUxYykqcj0dZIIuKuklk9EkLamXBtVzPRuIoIsy0ThdTKW3SVpU1/VkIMIu6sJgOnQ9BeQaUlR/2nc9BeQaUlT0I9blnmhIUeZd1xPATFOK4lUYX1AUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUFHX7MUVBRxrTaG8UdygKtc1jGme/oSjUsRCToChUdDMmQVEob0lMgqJQyqqYBEWhqLUxCYpqvEmcSWrdR7GYBEVhtRIxCYpqrijuJ7X+ifmyMQmKaq4HaU6TimMrxSQoqrG6e0kyiSqMqxyToKjGSsr/mdSMSVBUY3WzP5PCX+JpxCQoqrniTtLvFvg8tZgERTXYeO0qSjcmQVGtES/usCzEJCiqNfam0+lIsrIVk6Co1hhNZ35sKyZBUa3xk1lQP7W7GYpqAznM/Ux6Ghb56q8Oimq4yzXTOAsqtr5Biioj6tj+G67p+gJ8z/4OylBUKVHPdO3/Jddw61dzo43MnaKK6Mib1uJ++pQMXD/SdayeGiiAogrdez8y0bAfp0VNXT/SFVzHJCiqgHH2vrVhbDwuyouYBEUVkL2A3+tlB780LNeP9AZ/YhIUVUzvgfw8nnRcP9KrPItJUFQh/WH+Mn63wnsi7fAxJkFRRUhQ3pyK8jYmQVEFJHvpkzf24Xjnd0yCotbrDrLdU8/xIwwgJkFRa433HsRxPHBYVCgxCYryW1AxCYryVngxCYryUaAxicYU1R0M5Hv8u0NZR5v+IL+d7I3kHFKc30YP9uQV+GiU314MGywO6y8OM/F0YVh62y8xfMVWk6vDTNgxiaYU9e6gG2XPSvp1WXeQPhvJMEqyk0j9YdTP9l9xT15HyX7qZU/qqB8N0yc1/Skapk9qlA2T4VH6w5hOOiwbPnmQ/ciHZ8MeTLIfxmTDOjLs+lZleLbVYT4sLrHVn/8i7JhEU4r6ZUeeGHmi0hzkeU5bkCcqfT6lrlHfZHVln5M9z9nnZM+zDOtLDjIs/RzJI/2cfHhsJI/0c7K65sPG6W1WlwzvaGy1/6tf/+bh/sF7rv8s62lKUbKO6s8OJvHstjc7Ag7z22g4khdTuoP57bA3GxZfHZ4O61wZPh92OTwfdmV4oa3uPVi31TTjR4/33/9t0Fk1qigznh0w5rfzl+Hmt/OrAMwWLVHXlBqWLA4rO7zgsMCzalZRTRJqVhTlNdtZ7evfJUX5z2JWFLVUk4sSdrKiqKUaX5RQzeqD/f399/f3HyvPkaKCo5kV+6ilWlSUqJfVe/uPD2TfRFFLta0oUTmrh2lOH2a/oKilWlmUqJDVo48ODj60dUKCohqhXFYH71ucCkU1R+GsHmVFHViaBUU1TKGs9j8+ePzI0gQoqonWZvWBvVcLKaqxbmR18OHvPnlofbMU1Wx5Vp/+/tOHn3ydsXC64BqKaoFHf/jj1zO21uOXKKodPv4oL+pT6++3oqiWeO9hFtSf7L+Nj6JaI9tNycLc7tv4KKo90t3U5cLcWlYU1SYfXzutaSMrimo95awoChm9rCgKF1SyoigsqpsVReEWNbKiKCxTLSuKwkqls6IorFcmK4pCQQWzoiiUsT4rikJpK7OiKFSzLCuKQg23ZEVRqGsxK4qCiousKAp6sqwoCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrqaUNQ4SZI/px9j1xOBaUZRyTSXuJ4ITDOKMgMJamiiSa9rxhwAnWpEURMpqm/6ppeYfux6Ou3WiKKivTSovchE0cBkUcGhRhRlemlRvfQ22z8NXE+m5ZpRVDctKvtKL46zZVTH9XRarRlFmeFU1uNJL5mMJpHr2bRaQ4rqp+vyTDROF1OuJ9NuDSnK/MX1BDDTlKL+6noCmGlKUZzW9AVFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQdUhRUHL4dHxk5OnFAUFeUzHp8846qGuKzEJikJl12MSFIUqbo1JUBRKWh6ToCgUtyYmQVEopEhMgqKwTuGYBEVhhXIxCYrC7SrEJCgKN1SNSVAUrqoVk6AozNSPSVAU1GISFNVymjEJimov9ZgERbWSnZgERbWNxZgERbWI7ZgERbXDRmISFNV4m4tJUFSTbTgmQVE+OT5WuysXMQmK8sqJxp04i0lQlE+OzmregduYBEX55Pi0+lgPYhIU5Y3jzz57/qLSSF9iEhTli2z/VH4Z5VVMgqI88eLzssso/2ISFOWJl09KLKM8jUlQlCdOj415/uKztZ/nc0yCojxxeHZ0enZ0tPJTfI9JUJQ3nqWLqWX/L4yYBEX5LqCYBEV55fDkyfHR4cXvAotJUJRfTl6l0qy++CLAmARF+eXJq5m/BRiToCifHB59mff0/JnrqVRGUZ6YrZn+LkHpvU1q8yjKvasL8MOwd1CGohy7+dXcSdA7KENR7iw5NXD8rMqdeYSiHAjyPFNRFLVZjY5JUNTGND8mQVGb0JKYBEVZ1qaYBEXZ07qYBEVZ0c6YBEVpa3FMgqIUtT0mQVE6iGmOomojpgUUVQcx3URRFRHTEhRVHjGtQlGlENNaFFUUMRVDUQUQUwkUtRoxlUVRSxFTJdaLijrdTTwO3aKIqTrbRUU904038DjUiiKmmmwX1U9zGlh+DNO5mvdDTBpsFxWnRdV9pldLQ8o2MM1/WZE3MR0ezbieSGWhF5VHdJlU+W15E5N4cTzjeiKV2S6qkxY1tHbvWUCLSZXK16+YGsJWUfFwMBqO0px640nH2uynF8e6skl5G9Phy7Nn5sVT19OozlJRna7pm4n8shvZmvste6giSXkbkzg1Z0dy4eBQWTvqdZOxvX1T7mLnVDgpv2PKp3h4YrKogmWtqJ5JErtTn15WVCCpAGLKZfsnlX8UzRFbRUUjk1g+s7lwzmB6/T9cEUxMIv0yL11Grb9Uvq9sFTWZmK7dd5hMl+2hriYVVkzi6Ozo5dOXh/XvyBFr+6j0Y2x15ssPevlHgDHlDl+kiynXk6gu3PceXK3nek7TQGPKHJ4+dz2FWtKi/vHPEE1XrKPSH89c/8FW9OLs81evXE+ilrSod74K0eqDnt2XEq05fSqXCnY9jVpCLmrV0vw4PP/68t+zq+O7nkkt4Ra1Oqmj8Hzxn//OinI9k1rCLWr1Osr1vr+ao7NGHPVCXZk3cB2VvZvlJPyiXE+honXno4J1GvAbD0zgRS0/H+V6ci0WcFEr1lEU5U64RS19Mwu7KKdCLoqDno8CLmppUgTlUshFLUmKopwKuqjr31tlWJW7F3ZRtyYFpwIvSuE7QKEr9KIWTkcRlAfCL0rvShrQ8H8/dotR9aY92QAAACV0RVh0ZGF0ZTpjcmVhdGUAMjAyNC0wMS0zMFQxMTozNDoxMiswMDowMF/gwREAAAAldEVYdGRhdGU6bW9kaWZ5ADIwMjQtMDEtMzBUMTE6MzQ6MTIrMDA6MDAuvXmtAAAAKHRFWHRkYXRlOnRpbWVzdGFtcAAyMDI0LTAxLTMwVDExOjM0OjEyKzAwOjAweahYcgAAABN0RVh0cGRmOlZlcnNpb24AUERGLTEuNVM2WsEAAABKdEVYdHNpZ25hdHVyZQA3NzBiNjU4MGI0MjY4MmQwMTQ3N2Q3MzcyNGMwNjM5MWRiYTcyYjc5ZWVjNTE0MWUxYWRjM2M5YjQwZjg2Yzk3s3D9igAAAABJRU5ErkJggg==", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlEAAAGuCAMAAACHjqO9AAAAIGNIUk0AAHomAACAhAAA+gAAAIDoAAB1MAAA6mAAADqYAAAXcJy6UTwAAAE4UExURf///wAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA/wAA/wAA/wAA/wAA/wAA/wAA/wAA/wAA/wAAAAAA/wAA/wAA/66urgAA/wAA/2dnZwAA/wAA/wAA/wAA/wAA/wAA/wAA/wAA/wAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA/AAAOAAAAAAAhAACXAADYAAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAAAAAAAAAAAAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAEdHR1VVVf8AAP8AAP8AAP8AAAAAAAAA/////4iIiAD/AP8AACwe9TUAAABidFJOUwAAEVUzmYhm3cy7IncRd5kzIrtmRIhE3VXux83kqKpp1++PsKig9snuqvegmNKAkN7e0uKqiEQRM2a73SJ37lWZzJB2cNfkzfHEEUSIqt27ZjOZVcwi7ndQwMePes2Q75vmN35dUAAAAAFiS0dEAIgFHUgAAAAJcEhZcwAAAGQAAABkAA+Wxd0AAAAHdElNRQfpAQ8NHh0nAxlZAAAPQ0lEQVR42u3dCV/bRhqA8Y5syydNm4Bou7vZe0v2PgJNuj223RsIITRput3bdJvv/w1WemUDBmzreMczIz3/X4iTlrHG8RNpkI3y2msAAACoyFgVdbp2NwC/2C4q6plu7PpBYoNsF9VPcxrYfhDTnO3NoAjbRcVpUXaf6iylbAtTQ1Q+CL2oWUXT+QdNuWa7qE5a1NDe3V9URFK+sFVUPByMhqM0p9540rE2+/x4Z5qU1Nbrd94wb951PY3qLBXV6Zq+mcgvu5G1yc8ralJS98ydbXNvx/U0qrN21OsmY3v7pty0iUltbe2aLKpgWSuqZ5LE7tSvVrSYlN3tWpbtn3ZdT6IGW0VFI5NYPrO5tKawd1I7O9ky6i3X06jMVlGTiemOrM582tCktu9sv3739S3X06jM2j4q/Rhbnfl01V7K6pYt23ozXUy5nkR11l8ptmbVMircot4IeO+US4tK4hBNF49xix/hFrV9/vadN1xPopa0qHe+CtHqg16wSW2fp+7ecz2NGkIuqolLcynq/PztnTddz6SqcItandTdQH3jfOabrmdSVbhFrVpHBX7UO98NeH2eFvWt/4Woweuo++Ge3jQBnz2Yrkzq/u7dnXshfs20HfAKKhdsUevOR21t74SYVbhHu7mAiypwPirQrIIWblGXJa15EYasNirkosq8CENWmxJwUeXfcUdWGxByUfOkru+hVp86ICu7Ai/qtoNekXNRZGVN0EWZ6c1vhinxkh5Z2RB2UTe/v6r0yXKyUhZ6UWZ2wYP596lXuxOy0hN8UbPraJj619IgKxUNKEpovTRMVnU1pSjV77shqxooahmyqoaiViKr0ihqPbIqg6IKIquCKKoMslqPokojq5UoqhqyWoaiaiCrW1BUXWS1iKJUkNUFitJDVhmKUtb6rCjKhjZnRVHW1Mwq1EvlU5Rd1bMK9VL5FLUBlbIK9VL5FLUppbMK9FL5FLVRZbIK9FL5FLV5BbMK9FL5FOXIsqzeukwozEvlU5RLt2S1uxveYnwBRTm3kNXW+fn5Tng7pisoyg/zrL4tlwoOeTdFUR5Js/pOfjXzgHdTFOWX+7ML5H/3e6G+JkhRftnNr5D//R/8MNSXminKK1vnu/d3tmeHvDDfwUBRXrl5/im4rCgqBCFlRVHBCCQrigqL/1lRVIC8zoqiQuVrVhQVNA+zoqjw+ZUVRTWEN1lRVJP4kBVFNY7jrCiqmdxlRVEN5iQrimq6TWdFUa2wwawoqj02kxVFtYz1rCiqjWxmRVGtZSkrimo3/awoCrpZURRyWllRFK5QyIqicF29rCgKt6qcFUVhuSpZURTWKJkVRaGI4llRFAorlBVFoZx1WVEUKliRFUWhqtuzoijUciMrikJ9V7OiKCiZZUVR0LS1TVHQRVHQRVHQRVHQRVHQRVHQRVHQRVHQRVHQRVHQRVHQRVHQRVHQRVHQRVHQRVHQRVHQRVHQRVHQRVHQRVHQRVHQFXBR/VEuzn5DUb4IuKgFFOULioKukIuKJr2uGectUZQvQi6qb3qJ6cemM811XE8IJuyiomhgsqjMQIIauJ4PMiEXle2fTNZRLEVNXE8HmaCLiuNsGdUxYykqcj0dZIIuKuklk9EkLamXBtVzPRuIoIsy0ThdTKW3SVpU1/VkIMIu6sJgOnQ9BeQaUlR/2nc9BeQaUlT0I9blnmhIUeZd1xPATFOK4lUYX1AUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUFHX7MUVBRxrTaG8UdygKtc1jGme/oSjUsRCToChUdDMmQVEob0lMgqJQyqqYBEWhqLUxCYpqvEmcSWrdR7GYBEVhtRIxCYpqrijuJ7X+ifmyMQmKaq4HaU6TimMrxSQoqrG6e0kyiSqMqxyToKjGSsr/mdSMSVBUY3WzP5PCX+JpxCQoqrniTtLvFvg8tZgERTXYeO0qSjcmQVGtES/usCzEJCiqNfam0+lIsrIVk6Co1hhNZ35sKyZBUa3xk1lQP7W7GYpqAznM/Ux6Ghb56q8Oimq4yzXTOAsqtr5Biioj6tj+G67p+gJ8z/4OylBUKVHPdO3/Jddw61dzo43MnaKK6Mib1uJ++pQMXD/SdayeGiiAogrdez8y0bAfp0VNXT/SFVzHJCiqgHH2vrVhbDwuyouYBEUVkL2A3+tlB780LNeP9AZ/YhIUVUzvgfw8nnRcP9KrPItJUFQh/WH+Mn63wnsi7fAxJkFRRUhQ3pyK8jYmQVEFJHvpkzf24Xjnd0yCotbrDrLdU8/xIwwgJkFRa433HsRxPHBYVCgxCYryW1AxCYryVngxCYryUaAxicYU1R0M5Hv8u0NZR5v+IL+d7I3kHFKc30YP9uQV+GiU314MGywO6y8OM/F0YVh62y8xfMVWk6vDTNgxiaYU9e6gG2XPSvp1WXeQPhvJMEqyk0j9YdTP9l9xT15HyX7qZU/qqB8N0yc1/Skapk9qlA2T4VH6w5hOOiwbPnmQ/ciHZ8MeTLIfxmTDOjLs+lZleLbVYT4sLrHVn/8i7JhEU4r6ZUeeGHmi0hzkeU5bkCcqfT6lrlHfZHVln5M9z9nnZM+zDOtLDjIs/RzJI/2cfHhsJI/0c7K65sPG6W1WlwzvaGy1/6tf/+bh/sF7rv8s62lKUbKO6s8OJvHstjc7Ag7z22g4khdTuoP57bA3GxZfHZ4O61wZPh92OTwfdmV4oa3uPVi31TTjR4/33/9t0Fk1qigznh0w5rfzl+Hmt/OrAMwWLVHXlBqWLA4rO7zgsMCzalZRTRJqVhTlNdtZ7evfJUX5z2JWFLVUk4sSdrKiqKUaX5RQzeqD/f399/f3HyvPkaKCo5kV+6ilWlSUqJfVe/uPD2TfRFFLta0oUTmrh2lOH2a/oKilWlmUqJDVo48ODj60dUKCohqhXFYH71ucCkU1R+GsHmVFHViaBUU1TKGs9j8+ePzI0gQoqonWZvWBvVcLKaqxbmR18OHvPnlofbMU1Wx5Vp/+/tOHn3ydsXC64BqKaoFHf/jj1zO21uOXKKodPv4oL+pT6++3oqiWeO9hFtSf7L+Nj6JaI9tNycLc7tv4KKo90t3U5cLcWlYU1SYfXzutaSMrimo95awoChm9rCgKF1SyoigsqpsVReEWNbKiKCxTLSuKwkqls6IorFcmK4pCQQWzoiiUsT4rikJpK7OiKFSzLCuKQg23ZEVRqGsxK4qCiousKAp6sqwoCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrooCrqaUNQ4SZI/px9j1xOBaUZRyTSXuJ4ITDOKMgMJamiiSa9rxhwAnWpEURMpqm/6ppeYfux6Ou3WiKKivTSovchE0cBkUcGhRhRlemlRvfQ22z8NXE+m5ZpRVDctKvtKL46zZVTH9XRarRlFmeFU1uNJL5mMJpHr2bRaQ4rqp+vyTDROF1OuJ9NuDSnK/MX1BDDTlKL+6noCmGlKUZzW9AVFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQdUhRUHL4dHxk5OnFAUFeUzHp8846qGuKzEJikJl12MSFIUqbo1JUBRKWh6ToCgUtyYmQVEopEhMgqKwTuGYBEVhhXIxCYrC7SrEJCgKN1SNSVAUrqoVk6AozNSPSVAU1GISFNVymjEJimov9ZgERbWSnZgERbWNxZgERbWI7ZgERbXDRmISFNV4m4tJUFSTbTgmQVE+OT5WuysXMQmK8sqJxp04i0lQlE+OzmregduYBEX55Pi0+lgPYhIU5Y3jzz57/qLSSF9iEhTli2z/VH4Z5VVMgqI88eLzssso/2ISFOWJl09KLKM8jUlQlCdOj415/uKztZ/nc0yCojxxeHZ0enZ0tPJTfI9JUJQ3nqWLqWX/L4yYBEX5LqCYBEV55fDkyfHR4cXvAotJUJRfTl6l0qy++CLAmARF+eXJq5m/BRiToCifHB59mff0/JnrqVRGUZ6YrZn+LkHpvU1q8yjKvasL8MOwd1CGohy7+dXcSdA7KENR7iw5NXD8rMqdeYSiHAjyPFNRFLVZjY5JUNTGND8mQVGb0JKYBEVZ1qaYBEXZ07qYBEVZ0c6YBEVpa3FMgqIUtT0mQVE6iGmOomojpgUUVQcx3URRFRHTEhRVHjGtQlGlENNaFFUUMRVDUQUQUwkUtRoxlUVRSxFTJdaLijrdTTwO3aKIqTrbRUU904038DjUiiKmmmwX1U9zGlh+DNO5mvdDTBpsFxWnRdV9pldLQ8o2MM1/WZE3MR0ezbieSGWhF5VHdJlU+W15E5N4cTzjeiKV2S6qkxY1tHbvWUCLSZXK16+YGsJWUfFwMBqO0px640nH2uynF8e6skl5G9Phy7Nn5sVT19OozlJRna7pm4n8shvZmvste6giSXkbkzg1Z0dy4eBQWTvqdZOxvX1T7mLnVDgpv2PKp3h4YrKogmWtqJ5JErtTn15WVCCpAGLKZfsnlX8UzRFbRUUjk1g+s7lwzmB6/T9cEUxMIv0yL11Grb9Uvq9sFTWZmK7dd5hMl+2hriYVVkzi6Ozo5dOXh/XvyBFr+6j0Y2x15ssPevlHgDHlDl+kiynXk6gu3PceXK3nek7TQGPKHJ4+dz2FWtKi/vHPEE1XrKPSH89c/8FW9OLs81evXE+ilrSod74K0eqDnt2XEq05fSqXCnY9jVpCLmrV0vw4PP/68t+zq+O7nkkt4Ra1Oqmj8Hzxn//OinI9k1rCLWr1Osr1vr+ao7NGHPVCXZk3cB2VvZvlJPyiXE+honXno4J1GvAbD0zgRS0/H+V6ci0WcFEr1lEU5U64RS19Mwu7KKdCLoqDno8CLmppUgTlUshFLUmKopwKuqjr31tlWJW7F3ZRtyYFpwIvSuE7QKEr9KIWTkcRlAfCL0rvShrQ8H8/dotR9aY92QAAACV0RVh0ZGF0ZTpjcmVhdGUAMjAyNS0wMS0xNVQxMzozMDoyOSswMDowMAOSrigAAAAldEVYdGRhdGU6bW9kaWZ5ADIwMjUtMDEtMTVUMTM6MzA6MjkrMDA6MDByzxaUAAAAKHRFWHRkYXRlOnRpbWVzdGFtcAAyMDI1LTAxLTE1VDEzOjMwOjI5KzAwOjAwJdo3SwAAABN0RVh0cGRmOlZlcnNpb24AUERGLTEuNVM2WsEAAABKdEVYdHNpZ25hdHVyZQA3NzBiNjU4MGI0MjY4MmQwMTQ3N2Q3MzcyNGMwNjM5MWRiYTcyYjc5ZWVjNTE0MWUxYWRjM2M5YjQwZjg2Yzk3s3D9igAAAABJRU5ErkJggg==", "text/plain": [ "" ]