From 5231e2b5fcb0240fd2533610880550d51ba740d0 Mon Sep 17 00:00:00 2001 From: Alexander Puck Neuwirth Date: Tue, 30 Jan 2024 12:12:31 +0100 Subject: [PATCH] Add remnant Closes: #201 --- debug/proton_baryon_meson.ipynb | 115 +++++++++++++++++++++++----- pyfeyn2/render/latex/tikzfeynman.py | 1 + 2 files changed, 96 insertions(+), 20 deletions(-) diff --git a/debug/proton_baryon_meson.ipynb b/debug/proton_baryon_meson.ipynb index ff44b999d..f0629e88b 100644 --- a/debug/proton_baryon_meson.ipynb +++ b/debug/proton_baryon_meson.ipynb @@ -47,15 +47,25 @@ }, { "cell_type": "code", - "execution_count": 5, + "execution_count": 42, "id": "f4126da1-c09a-45e9-860f-342c333ec9af", "metadata": {}, "outputs": [ + { + "name": "stderr", + "output_type": "stream", + "text": [ + "/home/apn/.cache/pypoetry/virtualenvs/pyfeyn2-rvdu2REa-py3.11/lib/python3.11/site-packages/feynml/pdgid.py:99: UserWarning: Type remnant is not consistent with pdgid 2, using fermion instead. Using remnant now.\n", + " warnings.warn(\n", + "/home/apn/.cache/pypoetry/virtualenvs/pyfeyn2-rvdu2REa-py3.11/lib/python3.11/site-packages/feynml/pdgid.py:99: UserWarning: Type remnant is not consistent with pdgid -2, using anti fermion instead. Using remnant now.\n", + " warnings.warn(\n" + ] + }, { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAloAAAGxCAMAAACN7ugEAAAAIGNIUk0AAHomAACAhAAA+gAAAIDoAAB1MAAA6mAAADqYAAAXcJy6UTwAAAHsUExURf///wAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA/wAA/wAA/wAA/wAA/5OTk8nJyQAA/wAA/wAA/wAA/wAA/wAA/wAA/2BgYFpaWgAA/wAA/wAA/wAAAAAA/wAA/wAA/wAA/wAA/wAA/wAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAB3AAAAAAAAAAAPAABpAACvAADkAAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAAAAAAAAAAAAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAENDQ1VVVU5OTigoKAkJCVFRUX9/f56env8AAP8AAG5ubnZ2dv8AAP8AAP8AAEFBQR4eHgAAAAAAAAAAAP///3FxcSkpKUlJSTg4ON7e3gUFBQoKComJiaOjo8DAwN3d3aqqqru7uyIiIoiIiO7u7gAA/2ZmZjMzM5mZmURERDIyMhcXFxEREQD/AP8AAMzMzK+vr3BwcCUlJXd3d29vbwgICIhnBEMAAACBdFJOUwAAIbbWdo7HXKP14yKIl/3U9mYRd5kzInbtZkSI3btV7pH39sSqRPeO5MzHwBF3mTPdzLtV7s3fj6iq9+3aQZC0P6Wu6+faoHdEmcyIIrtmVaoz7hHdUM+2ksHXoBFEiKrdu2YzmVXMIu538Zh2Qc3Q5Pigf/G0qIj0wO/H1+s/+p1MWTgAAAABYktHRACIBR1IAAAACXBIWXMAAABkAAAAZAAPlsXdAAAAB3RJTUUH5wsDCR4SitlqfQAAEARJREFUeNrt3Yl/G2l5wPF9k3gTFg2lyibLJt5soKwKbXVYkgs1V1soxGBYWJZeEJskeOOmpQdtuY+dsAcs0IVkuY8W+o8y80h2LNuS5nrmPeb3/Xwc2RvN6B3rtzOvRor00EMAAOAQE4BY2B4FREBpSVTx/jewrYa0Tp0+s/Lw2XPKG3LQlXwRl301pPU688g5s/J61c1IUorN4baIy7pDabUiHW94wx9E0Rv/UGnt7XQrpKMjX563df7RCxfNYxdsD6OEQ2m96a6Sb79w9+4LLyqt/HEzCckE1tYlc+GyubRqexgl1JHWSy/effk7L+umdWJbtn+5ZZw//4RJ6/JXHWl998VXXvie1sofl1m7Ca+tdI/1hO1BlFFDWi9/9+73tcJK04rNvP1WbPu3W8bqqnnsinnS9jCKqyGtb7+kF5akFZs5bXk927p84fKjVx49b3sYxemn9cpLL7yimVY8d67l+RHx/GPJhMv2IErQT+sH/6MY1iSt+cfE2PbvtzCf91cT+mm9Gv/wR/eU01ow3/LV/asXLtoeQzl1pJV47VW9tOJ0O2IT2GPE+4krl2yPooya0orjH//knlJai6Zase3fb2H3xdXVx2wPpLBDab35pyp+Fk/9XGX1b4kXTePjK766P/VHb7U9koIu1bXX+sVrP9BZ/cJpvPd7rac4r7UkLa2j4d2Ap/F+Hw1NPWnpzeElrQWnTE189crqk17eQZ7P4U0dad27pxjW8lOmFy+tXvGxL8/PPJh6np5WT2vp8dDXvrzmf1qLuooPbyp91SuEtGKT/Wke+qpNAGktaGvORtNXHUJIa15b8cItpy9lQaQ1f7+1DH3pqeFf9OiSf9FzUltx5t8BfakI5V9PH20r9yu16KtqoaR1pK242Eroq0LBpDV9LxFT/iXx9FWNgNJKVfZGNfRVWmBpmahT4croqwzSWoa+CiKtTOgrP9LKjr5yIa2c6Csr0iqCvjIgrcLoazHSKoe+5iKtCtDXSUirKvR1BGlVir4eIK3q0ZcgLSX0RVqaGt0Xaakr2Ze3n01AWvUo3pe3n01AWjUq1Je3n01AWnXL3Zevn01AWlbk6cvXzyYgLXsy9uXrZxOQlmXz+rr4oCVPP5uAtFxwQl+rV/07BM4iLWfM9nXl/v2nPNxVHUJabjno62r6Prle77hIy0FJX2+dvrm3xzsu0nLTk9NPJPjjt3n7/CNpuWl1Utbb/+RP/8zX57dJy01X7j/x1OrlyeHQ09dPkJabLh6dZPnXF2n5xKu+SMs7vvRFWn7yoC/S8pjbfZGW75zti7SC4GJfpBUOx/oircC40xdphciJvkgrWLb7Iq2wWeyLtBrATl+k1RS190VajVJnX6TVPDX1RVoNpd8XaTWZal+k1XhafZEWUgp9kRYOVNsXaWFWZX2RFk5QRV+khXlK9kVaWKh4X6SF5Qr1RVrIKG9fpIU8cvRFWsgtW1+khWKW9kVaKGFRX6SFsub0RVqoxPG+SAvVmemLtFCxaV+XSQsaLl56krSghLSghLSghLSghLSghLSghLSghLSghLSghLSghLSghLSghLSghLSghLSghLSghLSghLSghLSghLSghLSghLSgJIC0WtFU+gNpOSOAtLqdqfQH0nJGAGnNIC1nhJBWt9cfmFZfvictZ4SQVtv0O6YdmW480bc9IKRCSKvbXTNpXaY/SYsdlxNCSCvdY5m15LIjZQ1tDwciiLSiyLRGZpz0labVtj0ciCDS6vQ7vVGvm+y+krLWu7aHAxFEWqbbSiZc6eV6HEe2B4OJMNI6kEzkW7bHgInA0mrFf257CJgKLC3zjnfaHgGmQkvrLzip5YrQ0uKJHmeQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFpSQFjQM2mPSQsUG7Wi0thENSAvVmVbVkU87JS1UYqYqQVoo63hVgrRQwpyqBGmhmEVVCdJCbkurEqTVGL0oVfL3k60qQVrIKEdVgrQaIGqP26VWkLcqQVrh20i66hVduFBVgrSCN1jvdHr5yyhTlSCt4HVG+ZcpWZUgreAN0rSy/1qqqEqQVviicac9yHLFyqoSpNUAreWtVFuVIK3GaXdaMz8rVCVIq3H6cbw+isZpX1pVCdJqnF488a53v0OrKkFajdOZpvWe9+reDmk1ihwBJ2VF2rdFWk3xYF41TMIaZjobUQppFdHNdp7IFUdm6/0adlmGtArpRmZQ+Pneep30GLBXwy7LkFY+Y3k1XdRO/qdfs72lS809s6D3oHAGaeUxandNd9iOkrRi21u6iOr5qqxIK4dW+oK6YWRcTsuJqgRp5TBOvvr99LiYFGZ7S49zpypBWvn0N+TPVm9se0tnOFaVIK1c2sPJnTdw5z50sSpBWnlIWe6c0nK2KkFaOXTWk6xaThwK3a5KkFZ2g7V0h9W3vYkeVCVIK7PW+kYURWs20/KlKkFanvCqKkFa7vOvKkFaTvO0KhFeWtGavL9Bd2P/cvI8/2C4f7kh91N7bXo5nLzAJFrbv5xZvDWaXXy8v9ha1F2wWHJZ5lZlEY+rEqGl9b6/7HeH6d2z0esO0z3YsC0P7Lrrg8FaVx7lycmpzrDbS0+sJz/003s36idLmPR5nO4wPbswardksXTxVlJYspgsvp4snvz1OFksnc/3RvuLn3yra8dutT1aeKvddLG/8r0qEVpaf/1+uUNN2k16r6b3YHKHJndacoeOjNzxyR0q92Dyd3LHJ1GkV0n/W3rHp00lVzlYfJy0lN7xySplsShKrtJKM5KrpKtMczvxVntFbrX7gb/54Ieubdr+VZYVWloy1xqsRfKETHu4IacK9i/7G7Jr2b/sDiM5bO1f7i82Xp9dPBqOosOLD/vr49nF1luztzoqd6tJkJvXPvyRj3zY676CTMu0xtOn+qZnzvcvx4OZy+64NXO5v1jryGKdzuzi7daRxbumkludXdz3vsJMKyTe9kVaXtDua0thnaTlD8W+SGu5kNMSOn2R1nLBpyUq7eujW1tbH9vaerrqQZKWt6rsi73Wcg1KS5Tsa+vpj8veirSWa1paonBfn0i6upZ+Q1rLNTItUaCvZz65uXntWa0BkVZQ8vW1+THNsZBWeDL39UyaltppWNIKVKa+tj6++fQzWiMgrZAt7eujajMt0mqAY309u7n1iU/p3y5pNcO0r79Nq/plQnUCP0FaDbJ57e9+OaVxIusI0mqWzU9N0vp7/dd/kVbDPLslaf2D/usLSatxZMcl3+i+vpC0mifZcT2Yxev1RVpNtHnt6M8KfZEWpqrui7RwWIV9kRaOqaYv0sLJSvdFWligTF+khWUK9kVayCR/X6SF7HL1RVrIKWtfpIUiMvRFWihscV+khXLm9kVaqMBJfZEWqnKkL9JCpR70RVqoXtrXFmlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBCWlBSUhpdddj0bc9EKRCSsv0J2m1TLfXH5gWiVkVVFoDKWtkTNv0O6Yd2R5PswWVlhmlabWTQ2N3zaR1waKw0monZa3LN8kea832aBourLRMMpGXw2AUmdbIjG0Pp9ECSytKJ/GJTr/TG/W6tofTaIGl1Yr/cfJNN3mYSFlWBZaW+fSK7RFgKrS0PsPDQleElhZP9DiDtKCEtKCEtKCEtKCEtKCEtKCEtKCEtKCEtKCEtKCEtKCEtKCEtKCEtKCEtKCEtKCEtKCEtKCEtKCEtKCEtFC569s7n72xQ1qoklR1c+fW59hroTKHqhKkhfKOViVIC6WcWJUgLRQ1vypBWihgSVWCtJBPlqoEaSGzzFUJ0kIW+aoSpIUlClQlSAvzFa1KkBZOVKoqQVo4qnxVgrRwSEVVCdLCRJVVCdJC9VUJ0mo2naoEaTWWYlWCtJpIuypBWg1TS1WCtJqjvqoEaTVCzVUJ0nLRzk5167JRlSAtJ92oZC3WqhKk5aLt3bJrsFuVIC0X7dwqsbADVQnScs7Oc8/d3iu2qCtVCdJyTbrHupF/MaeqEqTlmL1/yj3Vcq8qQVqOufPZPFMtR6sSpOWYWzvG3N57bvkVXa5KkJZjru9u39rd3l58HderEqTlnCSY+Q8Q/ahKkJY3PKpKkJabdnZuHToo+laVIC037f4qcXv3zraXVQnSctOdX0398794WJUgLSdd//y0rApfXVM30nLN9AgoYd32dIclSMshh+dVt/3eZRnScsWx2fqu37ssQ1oOOPkx4B2/d1mGtOzy98xCBqRlSdBVCdKqX/hVCdKqVUOqEqRVlyZVJUirBo2rSpCWrmZWJUhLTYOrEqSloelVCdKqGFXtI63qUNUM0qoEVR1HWmVR1RykVQJVLUJaxVDVUqSVG1VlQ1p5UFUOpJURVeVFWstRVSGktRBVFUda81BVSUGlFU+VXQ9VVSGctCZJxQffFUNVlQklrQddyVeRutyq6uaUG6MpIpC04sNdTX6McyzuVlWBCCOtI/usPG25WtX1O7ufM3ulP1DFoiDSio91laktV6sSt8zutrwds7dCSCue9xXPXcTpqiYjvH7DpHX5K4i0YjO/r2Pcr2oi3WPdsD2IMgJIK57bVXzkmr5UJXZ2zN5Nk+GzCVzlf1qL9lnxwbW8qkps727fuXnnuu1hFBdCWrGZN9dK//Swqonre8mEy/YgSqglrfGG4srjBfOsOPazKrFX6gM37TuU1r/+m4p//48vxLHOqlP/GZtFbflZVeLWzYDey/RNdxX8+jV5wlhj1ROPx+lmxCbHY0T37e3c8Pydvc1MWv/1atV+87//N30tQuWrPvDfi+ZaJt720Oe/OH3T+C/ZHkkZe4fS+vJvq/a7l/5/mlblqz7wlXjxXMtDX/3aNK2v2x5JGdvKB8R7P/lxHQfE+W3ZPioUs73r+0ddmDqm8d/4Zqw+jQ9trjWdbQWTlppub01xCxbPtRRvWNutm6RlVxze8TAMNaT1urNiRWv9cYDHwxDUkNa5FfEtrfXH879itY3CUv4fEM2iuRbsqSGtU6fPrDx89pzeNsxrK9a7SSxXx1zLPHLOrLxecSPmtBUr3iSWqmOvder5pK93a27F3GMi7KljrvXwGWPOnFbdjONtMdGyrY60zp42p54/pbsdR9uiLOvqSOuR0ytnFGfxU2lMseGxoTPqmGs9YtTOac2I46reTwQVqCGth8/a3kjYoJ/Wytkzas/xwGEPAQAAQM3vAYHVHziMraDtAAAAJXRFWHRkYXRlOmNyZWF0ZQAyMDIzLTExLTAzVDA5OjMwOjE4KzAwOjAwT5vmTgAAACV0RVh0ZGF0ZTptb2RpZnkAMjAyMy0xMS0wM1QwOTozMDoxOCswMDowMD7GXvIAAAAodEVYdGRhdGU6dGltZXN0YW1wADIwMjMtMTEtMDNUMDk6MzA6MTgrMDA6MDBp038tAAAAE3RFWHRwZGY6VmVyc2lvbgBQREYtMS41UzZawQAAAEp0RVh0c2lnbmF0dXJlADQ5MmFlNGFmMTliMGUzZDQ0M2E4MjFlMzI1ZjhjMDU3NjY3MWRhMDk2YTc4MzUyYWY2ZDk4ODZhMDE2NzVjMmZFi/y2AAAAAElFTkSuQmCC", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAEZCAMAAABcjom2AAAAIGNIUk0AAHomAACAhAAA+gAAAIDoAAB1MAAA6mAAADqYAAAXcJy6UTwAAAIxUExURf///wAAAAAAAAAAAEtLSwAAABsbG4qKiktLS2pqaq+vr46OjgAAAAAAAAAAAAAAAD09PQAAAAAAAAAAAGpqagAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAIyMjAAAAAAAAAAAAAAAAAAAAAAA/wAA/wAA/66urpKSkgAA/wAA/wAA/3d3dwAA/wAA/wAA/wAA/wAA/wAA/2dnZwAAAAAAAAAA/wAA/wAA/wAA/wAA/wAA//8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAAAAAAAAAAAAAAAASgAAXAAAtAAA6AAADAAAGAAAAAAAAAAAAAAAAAAAAAAAAAAAAP8AAP8AAAAAAAAAAP8AAJsAAAQAAP8AAP8AAAAAAAAAAAAAAAAAAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAAAAAAAAAP8AAP8AALkAAA8AAAAAAEcAAAAAAAAAAP8AAP8AAP8AAAAAAAAABgAAXQAAdAAAmAAA/wAA/wAA/wAA/wAA/wAA/wAA/0dHR1VVVRMTEx4eHkFBQTo6OrKysnl5ec/Pz4ODg01NTY2NjWdnZ3d3d93d3f///5mZmaqqqiIiImZmZlVVVYiIiO7u7ru7u8zMzAAAAERERDIyMhkZGXBwcK+vryAgICgoKF9fXxERESQkJBsbG1JSUgAA//8AADMzMwwMDFdXV6KioggICDg4OGdnZwK8xvAAAACZdFJOUwAAEXetRKA/5tr+mSHHdqPW9dbj9TP6tmb+XI4/7feIVd/vmHeIRMfOEWYi1qpVu90z7qiqmZiCkMdg4qpEVZndEYh3Irsi3Vni5sCXnu7MsKiP0mn6ZjHEM+nVp7u4zWzP7vHkzK1OyeIx4eNcz/2FoPb25MD15rrxWsGFkLjn0+nx26iZzO+A86DGes3y69fa9I/+Iexc0iRX2JAAAAABYktHRACIBR1IAAAACXBIWXMAAABkAAAAZAAPlsXdAAAAB3RJTUUH6AEeCwk71/mH3wAAE7ZJREFUeNrtnYl7I8lZh7da9BLAkmNPZGvHG0qTjJTMOJ7ZGSQIYB0tqSUIICWBCASIgIMJLCyEm3CFM+HOwRFmvdlkd9MJSQjHiPuvo6r60NVHlVRXq/U+j3tsq9WS+52vf1XV1a1nnnkGpB4j9w2q3wIfnkm/DfPZJy/ffKPqd8GFtNswcm965ebm5nPfpPqNcCHdNkhZ3Lz85NXPf7Pqt8KFFNvwy+ILr772+s3Nt6h+OzzQzMZBPp8vgEI+f5i05rws3ri5ef2VN5uq3zoPNLMBwNExWtwqxK+0XBZffMuzRdVvmw/a2Tg8AeD0IHaV1bIoGarfNC+0swGeO8zHydjVsiDoZyN/cjv6wd0tC4J+Ng6dKBs7XRYE7Wyc3b51HPrAjpcFQTcbhTNQeG49N3a/LAia2bh9hhYnpyu/zUJZELSyUTh2UEfj7Ogkv/DLjJQFQSsbIWSmLAha28hSWRA0tpGtsiDoaiNzZUHQ00YGy4KgoY1slgVBOxtZLQuCXjYyXBYEnWxkuiwI2tjIelkQNLGxLwuCDjb2ZeEj3objEfX4vizmCLbhaXDcb9cf35fFEvQ2Dk5vFc7ytynXdpmrIF+rOvZlsQK9jUNw9Dwo3GLYNnaxomPBx74s1mGojYMTpOStDNt2Qr589mURBkNunB0DcHxKuzZwwnXgr31ZRMBgI38KDk4OaNeOUuHsyyIaBhtHp4Vj+hB3InV8674somDIjSPwPP12ow5SKMedfVlEQW/jLE+5IsGJ0oEe2ZdFFNQ2CsfHBbo1CZEHKqxjXxYRUNuAZZbNOjFfwNmXRQTIhglpuOO8Lf92qjUxTkxuAGd3y+JuBVPd9OnIxjtuaPgSHnLC4UtF3IFqbYBkjw+y8c43aPiyOxL7T1/5Ks3aTtyBahdt3Kvcr97fchvn58jGu57Q8DVvYPyfv06zdnxuqN51ArhAKu5u/vTzB5WHLzyssh2p/gX3oqmOVLG5oXrX8efR42r17r2NnuqKqDw4B+RI9W2v0fCvjvNv/061JiZruVF9uMmzFkUQqFu4tTpLqzRrufEI22BpSq2JIAg69+c4GcuNyoPq/Ud0q4aLIIg6Exs3TqVibwnnnCY1YkQQhNnI1IGKgiQRBGGzFGJbVTvPxQsXlapfLVQiCOJsRJ9uygCVp5gXLr79O95NKYIgbgaPk93KQM2rpx7f+V2UIggC51M5ma2M8wfvdl08Zhw/FDm7zclgZfgZ8d1YxgVr71zoXENndXbbTstYCusLVBgPmDchfOYnmB+kdtfFequpwl4YQMKsaCfgsiF5H0khovn6yC8M7/RbiWpjEq8YaLZkvZIkqPoRpgvdIJ/M6zfalrSXEg19h44JmTY6XWkvJRBWERDHZU2/2ui1pb2UGDaqiKINQH9At67UK81SnePf873v+b5NDk05iEqD8m69Um2kNccNfJz5/h+A9qAGTdbpR8hECeJv+slK5F6F2bYkvhg3+kN7UC+PsJJiiVnJEORI8xY9KZe0rlwb6cxx3FkoQxj8zKSkOCr1yb82RTNXro105jgMO8SwVomJbCauK/l68VTmeD3ygM+gxIBmOXEqs2QbynO8DyHs4+EKE/1vhUm5WrJrJQOM4nd2pBL26cay76XQtqS+3DomnpNtOPi/aaIMMDSKcDCEFJtdUWKYsDZgHyWVbUN5jhvkMlAHi8BNnRwZ04vSYrNt21Xy3ve9d0CGSUfMb062DfU57qCWTWloujLiMRltEIrvf5s3aF1nfq70u74oz3FUFsWiXU7ojJVyRbcdtAFGzbXxg8xdRek2lOf4APWNgQ3jSyNXy9UGcMR0PdcC5SG28UPMXUX5d0RqW5JfcAW7hHaPDeO7YrhR2zc3/ywoUh7kO6Z+iXwbqnPcxl1rmHAZo93f9mXKw4UQp1Ui34bqHLfxXoLD+B2zSXyvYKweCimUKLh3m+IcJ81ZGDmAZ+MRQnMo6tXjlSiwoTjH+8Ei9FFUOeW6LfbTAyOVqLivYduS/pLUbNio3YAwJSpsqM7xOOTZIKwoUWFDdY6HY9TwCGEu8YwQf+ZKlNyBVXl/PAwIi3BgD7Zu2m4IVgKV2FDeHw/D1uB+G2ruTty2VP/d63DoYmyNGht65XifjJII62IwoMaGVjluDGBtAGFN9fsAyu6jrlOOk/kDpZwGsaHKhk45LrmLEYeqzxhoW6r/8gC4t6FFjpcHw3qZepKsBFTZ0CLH7T4eIWQ/fS0MZZ+GokOOU87jl4cyGzrkuA5djCXUfVJQ21L4Z+fqqEVr6ND9XkKdDZU5XrRNOLBH+jSmPNTZUJnj+Hoj0NfvFskKP9MM53jPUvLSbDe+lodCG81W83ImvWVFRgiT50OrQeGRqvOB2Uy6DVizB9Dkek4J2gQem1Jlw7qcEXqSX3eARwh1PVCpq42xa0Pyq/7wj4jYaqls8jmdru5I1VVgY/KjPzYa2HVYMnkeqvo1UOTTWFaY4i3pNibtCcDXHZUhtDlaMYY2p5kmkm30Gi54rxAdMl/clRHAz4pRHtHd8SgJyTYaHZcx+akl1caKjIBtrRSHAJRTaWOFH5doI0pGwKZWjJzJSYZ0G1a3NQG9qf9jW1p/I1FGALuVIq9z6rJtjEGrAcYd/0f+4+oTkktr3Rh6GQGJVvolyttOUSO9Nqw2wEJ82hbnFxiTXFotuQ1kBIRaMcoQX4acchukLhZGb+WMq28jI2Bu5e0/8UHvImTO71O6jU4Hx0bT/1HAuHp33Ogs/4aLjADD/Mmf8mQcce1FKrDRaDW6064V/Mz9/DiqtuZ06Td8ZRBMm9i4w7UXqaKFa/VQeMx/5J3jvavV3wiQAfDsH2SDzBbl2LdX29/AtC2um+tOV34hRgbA19U6S6NTHKyot8E5x5s4MxaOfsJkoN1fDx2e2sKKehu8c7zVaIznAgTKSGIDK+ptcM/xxXPtCmUEMFjRwIbAeW46yAigsKKBDTHz3KYTzWQExFjRwYaI/vhkdjXRU0ZAiBUdbIiY59aazT7001rLCJhbMXWwIWC+unU1QzrSYcPHMPta2OCf4+4UiKt06dDjSCUgx9vu/KBr1X8YI3rY4J3jDeJChwt22NDDBu8c/xl0lOrInsbIAT1scM7xD89+dqz6L9oITWwk53iz0+xQ7uLJz31Y9d+zIZrYSMzxxhQVUJNmS5r2wKnQxUZSjl+OaVteKZahjY2kHJ/1QI+qvZpmGdrYSMrxK9Sj61BsJtUy9LGRkOPj8fiaIjbSLUMfG8mpkLhC6mVoZCM2xztN0LhM3ETaZWhkIzbHu83mwhysCFIvQyMbQY5bjfmy5y4n8+Wkt7r0Sb8MnWz4OT6d4mNSe4qaUFb7soOr5hodxSZXOMYb19dI0Pj6ConoTNvW/Ok7IEMnG15M96bgGk8QxY3aZsdC33e71hRHB06OVhNPmJr28OTqqbUwVrILMrSy4eZ4o4N2NlmSKdQztGyQ78lyOsFOpvhxpKsZdEF2QoZWNl78ebxEtdEmtYGXXm3gi6E6Y/Q7rzYuJ93uvDY60+n0F15S/e55oJON8S+S4phe4gC5vlzIDTxKhdKD5AYugua03cNTbtsWwGUB9rXBn+sXyS5tzvCyS5admYWWl22yJONU1zjikSWyxAcq65q0u3ZCh0Y2JlO0S1uNVrvbbjWn0267g5YdvGyh5biNlmO0bLXH6KFWu4l+OZ2ipZfku6BDIxstvFvHHXfZwJfPLi17+LY93tLCS9wdbHTmncId0KGPDWvrc+Pp16GPjTHNgHk8qdehj41rDnM80q5DGxuT6fbbSL0ObWy0+My5CdGxfsWytuhiY/sM91jTsX7Fsr7oYoNDhnus6Fi/YlljdLHBI8M9lnV001IXGE1s8Mlwf2OLOlavWNYaTWxwynAPV8fEcje9dMWy3uhhg1uGe2AdkysvihTdHXwT9LDBL8M9Ju2XrmY6fPwNG3rY4JjhHi/90mw2o5tFrRFa2OCa4e4W8VWYszQ1pwha2OCb4cCXIewu7dzu1L2CFtfE8s5wFNwtV4eYERF+d+r2ME0IbXtol3SwwT3DAT49hS+L/ZCQ98vtTt2+hRrMkY8F0eJIxT/DCY3L2ewjQra87Z261yz4aGCDf4b79H75VwT0+za/U3ekBR8NbHDP8AVEnO9gv1N3ogUf9Tb4Z/giQk4/0d6pm9qCj3obIjJ8AWFnA2N2L7MFH/U2BGV4gCgdzrBWWrkD28YWfJTbEJfhwSuI0UHuTzwiH8i1tQUf5TZEZriHGB22d/du51e3t+Cj2obYDPfgrYPUwh3Xxa/x/LQ61TYEZ7gHJx1LRySIXQz5fvCvahuiM9xjOx2huWAiGXVeH0rjodiG+Az3X2gTHbHpbDo29w9hVmxDQoZ7sOigayMJ+HhTtTakZLgHhQ5uLdVNUWtDToZ7ROtQbsFHnQ38eZiSMtxjTYc2FnzU2ejMph/5dbkv6evQzoKPShuz2W+05M47m/zmb92581H9LPiotSHs1HUUv/07v1seKfqDKVBt49qS/sJ7GyF0FMkAtqI/mAK1NlTIAKQLzfljEzmh1IYSGSAHy/VBbWDSP4Pv0GAM6mz8niIZwID4FFE5ahf3IYR9gBYmMPECZMLG73/MUvTKBDNyF5uOgZQ5eBjKlZEFG+0XFb2wSzEyzA2HLLCIklsqNoRbzWajRpkNfyzd6qIeYE/BbPIBMEr2KGwKJ/7U8NLQJDIIu18b/lj6GH/2u9TRQw84HMG+MQh5BJVFsWiXQd8P+p23EYylW/i7lrLLJMPOGKHmVgnYMCiN3bcxrwb8nbprwmBIQ9cumQayoWAoS5WN+Vh6p4NjQ9VFYbixW4TL6WHjzCbNYNkosrFwPrzRanSnyTeCFoQ5goNabvmjqW08kgWHnGcg0KDIxuL5cKuHwkPFmyAU0VEJFGuLvyJHLyjkWrIE1NiQeT6ciqXORz9YyEaNDRUt2liIDfUDiWpsyD0fTsGgXK4P7O3OfEDch6+lb1a0tDlt1BRruAm13a7EYy39wTZbUGND3pw2NiDDKPs6OdRJrG2+BTNXV3KFsnYZ7rOdDWSihHvtfeatIBGjUT1nKqkN7TLcB+3Mfq5W2zDNhyCHx1KgaTA0jgMRBBU2tMtwH2MwGuSKufpGTy6OyIVnRdukHFJZFkFQYEO/DF9hq2kM+CxWYi8+RARBgQ1dMzxgKxsGNMtxI1xRIgjybWib4QFD9FXM8R/BjRVBkG9D2wwPqEE4qMHBNl3zUmlZZrIIgnwb2mZ4gFEqG1s2d/FFgTYkSihFEKTb0D7DfbaamFDyL1/+6B9QiiBIt6F9hvuUtjn/WvRkfJDtcCfbhv4Z7mPWjHJ9MGCvEHJoIi6GrOdIZNvQP8MDoF0vA4OpuTvPCHynBZu5HSDbhv4ZvoJDud5KWEP2wgDSbaQmwwMoaiOs1VQarRVGJXlLkm2kJsMD0OGmn6tF/T9nab5qZyM9GR5QGoxGOTPkvAWLCIJ2NlKU4QF9PAa43NxlFgHOK5XKw0rlfsJqcm2kLsN9gusL2EUE6FYb6ctwH2xjQxH3KverpCZ0s5G+DPcxBrU//KM//pP7j9ifeoFU3MXfaGYjhRkeUC6i/+XVysULDytsSh49rlbv3qNcWaaNNGb4GqxKqg8ZNi7TRmozfA0GJY+wjSrldiXaSG+Gh0OppPKgSl1I4m34I/1OejM8Bgol57SpIdwGkkDG3cgX7QBc6ghVUr37p8wbEmtjQQVw1ewuC0qQiI8/ffr0MfM2RNpwllQAX84ug5V84qkHS2vKRXRtrH/tPn/2566Nv2B+pjgbTpiK3a8OzPlDYuMvmbuKAm2EVsZuZ0fA3cfIxjlzV1GYDSfiKxPFQcrDD3EGJaJsOFnXAe5eLP5Ep0ScjUgdGbERQqISgTYiciPLOjBxSqht1HMs9xaIPVBl2wYhQgmy8Vd/TcPfOM4nP0W1JsaJi/G9DZd1JcjGO25o+BIe+PuPV16nWvkmNjP2NhZYUoJsvPMNGr7sDsT+53/9N83azj43WPCUVJGNdz2h4Wuujf/53/+jWXufG+zcq1bFHan2ucEOU4p/miXF97nBjqgWbmxuvIlpWxlCWO8vLjduXn7y7FY3kdhV6G0cnN4qnOVvU64dF+Nf+BxKllf2BbIOvY1DcPQ8KNyiXDtunOqNz7/65OWbfYGsw1AbBydIyVtpNxylw3E+c3Pz+muv7gskBIbcODsG4PiUdm3HiTpQld6M28ksBXKQz+cLoJDPH6reXYJhsJE/BQcnB9RbjjhQ4X+Kz77li2wFcoT+J4BbBdV7SzQMNo5OC8e0IY6JO7dhMBbIITpKntL/T0grDLlxBJ5n23bIwWqh48dWIM8d5ndfBoONszz71uNnt7EUSP6EpSzTCrWNwvHxBodtZ6EXGDY+RV0gh87eBg+CWdHhD9MVyNntW8eK95QMFH++OCG5QApnoPDcPjckkVAgt8/Q4oS6q5Ne9LAB4gqkcOygxDo7OtmgGZEytLHB3gfZQTSyATbppO8WetnIeoHoZgNkukA0tJHhAtHSBshqgehqI5sFoq8NkMEC0dpG5gpEcxsgWwWiv40sFUgabIDMFEhKbGSkQFJjA2ShQNJkY/cLJF02wI4XSOps7HSBpNAG2N0CSaeNkALh/xFkCkirDbBaIH+r+u3wIMU2lgrk7/5e9ZvhQaptgHmB/IPqN8KFtNvwC+Szqt8GF9JvA+AC+UfVb4EP/w+nuFJ1N428ngAAACV0RVh0ZGF0ZTpjcmVhdGUAMjAyNC0wMS0zMFQxMTowOTo1OSswMDowMLWxdqAAAAAldEVYdGRhdGU6bW9kaWZ5ADIwMjQtMDEtMzBUMTE6MDk6NTkrMDA6MDDE7M4cAAAAKHRFWHRkYXRlOnRpbWVzdGFtcAAyMDI0LTAxLTMwVDExOjA5OjU5KzAwOjAwk/nvwwAAABN0RVh0cGRmOlZlcnNpb24AUERGLTEuNVM2WsEAAABKdEVYdHNpZ25hdHVyZQAzOTBhZDU1YzNjMDRhNjE4ZjkxNjkxZWM1M2UzODU2NDhjNjlhNWMwOWUxODRhYWZjNzNhMGNmZTc2MjIyM2YzEUvjUwAAAABJRU5ErkJggg==", "text/plain": [ - "" + "" ] }, "metadata": {}, @@ -64,36 +74,101 @@ ], "source": [ "fd = FeynmanDiagram().add(\n", - " v1 := Vertex().with_xy(5,0),\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=23).connect(v1,v2),\n", - " Propagator(pdgid=-2).connect(v1,v3),\n", - " Propagator(pdgid=2).connect(v1,v4),\n", - " Leg(pdgid=2212,target=v3,sense='incoming').with_xy(0,5).with_style(\"double-distance: 15\"),\n", - " Leg(pdgid=-2212,target=v4).with_incoming().with_xy(0,-5).with_style(\"label-side: right\"),\n", - " Leg(pdgid=11,sense='outgoing').with_target(v2).with_xy(15,1),\n", - " Leg(pdgid=-11).with_target(v2).with_outgoing().with_xy(15,-1).with_color('green'),\n", - " Leg(pdgid=2,sense='outgoing').with_target(v3).with_xy(15,+5).with_color('blue'),\n", - " Leg(pdgid=2,sense='outgoing').with_target(v3).with_xy(15,+3).with_color('blue'),\n", - " Leg(pdgid=-2).with_target(v4).with_outgoing().with_xy(15,-5).with_color('red'),\n", - " Leg(pdgid=1).with_target(v4).with_outgoing().with_xy(15,-3).with_color('red'),\n", + " Propagator(pdgid=21).connect(v1,v2).with_color('red'),\n", + " Propagator(pdgid=-2).connect(v1,v3).with_color('red'),\n", + " Propagator(pdgid=2).connect(v1,v4).with_color('red'),\n", + " Leg(pdgid=2212,target=v3,sense='incoming').with_xy(0,5),\n", + " Leg(pdgid=2212,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", + "\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", - "fd.render(auto_position=False,render=\"tikz\")" + "scale_positions(fd,0.5)\n", + "fd.render(auto_position=False,render=\"tikz\",file=\"hvq\")" + ] + }, + { + "cell_type": "code", + "execution_count": 40, + "id": "4aa7709f-b9aa-4723-a6fa-faf25de62455", + "metadata": {}, + "outputs": [ + { + "name": "stderr", + "output_type": "stream", + "text": [ + "/home/apn/.cache/pypoetry/virtualenvs/pyfeyn2-rvdu2REa-py3.11/lib/python3.11/site-packages/feynml/pdgid.py:99: UserWarning: Type remnant is not consistent with pdgid 2, using fermion instead. Using remnant now.\n", + " warnings.warn(\n", + "/home/apn/.cache/pypoetry/virtualenvs/pyfeyn2-rvdu2REa-py3.11/lib/python3.11/site-packages/feynml/pdgid.py:99: UserWarning: Type remnant is not consistent with pdgid -2, using anti fermion instead. Using remnant now.\n", + " warnings.warn(\n" + ] + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAASoAAAEZCAMAAAAXCpgeAAAAIGNIUk0AAHomAACAhAAA+gAAAIDoAAB1MAAA6mAAADqYAAAXcJy6UTwAAAG/UExURf///wAAAAAAAAAAAEtLSwAAABsbG4qKiktLS2pqaq+vr46OjgAAAAAAAAAAAAAAAD09PQAAAAAAAAAAAGpqagAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAIyMjAAAAAAAAAAAAAAAAAAAAK6urpKSknd3d2dnZwAAAAAAAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAAAA/wAA/wAA/wAA/wAA/wAA/wAA/wAA/wAA/wAA/wAA/wAA/wAA//8AAAAA/wAA/wAA/wAA/wAA//8AAAAA//8AAP8AAP8AABwAagAA/9IAHgAA/wAA/wAA//8AAP8AAP8AAP8AAP8AAP8AAEdHR1VVVRMTEx4eHkFBQTo6OrKysnl5ec/Pz4ODg01NTY2NjWdnZ3d3d93d3f///5mZmaqqqiIiImZmZlVVVYiIiO7u7ru7u8zMzAAAAERERDIyMhkZGXBwcK+vryAgICgoKF9fXxERESQkJBsbG1JSUv8AAAAA/zMzMwwMDFdXV6KioggICDg4OGdnZ5Xg6xIAAABzdFJOUwAAEXetRKA/5tr+mSHHdqPW9dbj9TP6tmb+XI4/7feIVd/vmMfO1qiqmapEVZndEYh3uzMiZszu7d/TdmCg2Y/45CJVqt2ImUR3EWYz7rvoXZD9kvP05v0x4enPYKCt8YXA9ZC453rN8uvX2vSP/iHsXNIgVEk9AAAAAWJLR0QAiAUdSAAAAAlwSFlzAAAAZAAAAGQAD5bF3QAAAAd0SU1FB+gBHgsGLNOyHtcAABBdSURBVHja7Z2JexvFGcYzoUlpK8lYBptYofQAWva+VdrSIxDHTgylN73v+26iECBhoS09Re8/uPPNSrK01q5mpZ35vjF9n8eKbK12tT+9886xs5Nz584x43X+vndoOMo581FduHjz1uid6o9jOqrz991/ezQavfQu9YcyG5Uw1OjWzTsvv1v9wQxGNTXUK3devTsavUf58Yih6nS73R7rdbtbq7Y8MdS90eju7QcuKP9sxFAxtt3nDzu9+o0WDfXagxcf0vDJyKHa2mVsr1O7SdlQD5/X8snIoWKXtrp1pFAMJUQPVXd3v/pFJEMJ0UO1lVehwjOUEDlUg/2d/tIXMA0lRA1Vb8B6l05nFbKhhIih2h/wh9290l/RDSVEClWvn/MG1WB7tzv3RwqGEiKFaoloGEqINCoyhhIijIqQoYSooqJlKCGaqKgZSoggKoKGEiKHiqShhGihomooIUqo6BpKiAwq0oYSIoKKuKGEKKCibygh9ajyiapeN8FQQopRTRjlxdPTrxtiKCF5VJ29nd6guy+59ZTUlJP4KbMyxlANUW2x7cust9MQVInVHCyTDNUQVaezy3k90gTVkp+pzDJUQ1Rs0Gesvye7NcuXs4If4wzVFFV3j3V2O7JbV3HKTTRUU1Tbe72+fKrnlazea56hGqLqbLPL8vutKns82HPzDNUQ1aAruWGBqooVf8U8QwlJo+r1+z25LQtUrIaVcYYSkkb16Pua7Dav+WG5cYYS4qguPCqj9+cf6H5QaktQXpNVLDfOUFNUj41k9Dp05yCNpVRX/k71bwwRR/X4PRm9UYwP/O73f5DZOq8rf+aieuKmjN6cjKX88U8yW9dnFfZJr42qQQH8M7SzpQpgbVZhn/TaqD70qoz+kud//ZvUlqAzmlVyG374ySYV/BnNKhW7zfMzmVVK9lvXB8Q+5zWlDNWZK3/qLkPU1oFGSh2q6rE9Q6Xu4lZ+tjyl9DpgfqY8pfaSaX6WPKX46nJevmRqMin1F+LZSdkzGpSG6R35TJaNfbKbSeOkIcfFPtnNpHN+ledjn+1G0okqCLHPdh2FUVw80Ykq8bBPu5nSKOWPnj352FonOBoT7HYW8Ecv8RLGIuaKj51qRWVAsCee5cN3ajmMZSIyPBYIVJ7eabPqgt0JnCDdaA9+lnEPBSngifj+OCofHq1UhFXq6kWlLNjtiJ+Ts+67oyH3TRo4Fv+EdsJ35do2PHfdVCCEAwwTvaiUBbuVrmfZYMjfGFsx/1xBwIa8vgu5pZhjpYFAyB/iLMuCIIt1z1tXFezjhCVZo3eEQwhOz+eA7YBFPgRpBF8meIu5fIcsiTIvCDybJanja5/iryrYwQ+B7GfwMn7iQx8iKILM5pAj/hV6Dng+HfNy7FuZ5waew3zHmXlV990QioI9TdNsdVQlHtT+VgIVwJhBGfN8qN6iFCDFw5A/Bp4XBR5/ZjvJ4rt1o1LXYq/8EuCMfQ+SO0hT7urIh8yOEl7ueHKDt/wMCl0IoQQNqthOlu1HNyo1wR44DE6/LNEjyTzuITuIM6jeII6CNOaPaQZVAbNEofMsz414g8q34+qDaL8dSUmwh44TLphKHISfP3/m+h6g4g1KKHIRVGkBxJo75O+wOSQriuBZvOog2lHxYLddhR2cGJAFoi3kMd6otF34driPIiiGIXgvHccAzPUiK4NnKyEhofI/8tR4rARVDAnjiOSOfEjuDAomBwbfDI9q2GQc8gzI3MiLeBuAJZKQUFDZ7hjUMipxynEGyR3ZULFFUOK4pwBY4ELg+1AYfcuyOCRoDTSChIKqINUeKh/25HsZIHEg2a2iY+IAKkf06pjFo8gPoLxBG4CtAQkFFWvNVb4NocTbQJDZkNzcQaIh6kcwiCKqN3cYQ0OJJ5cHQZ5sdkTtWeVujkq0DSM3E73kSccEyjag4gUNqIRwiNALLM8L2mr06r93OdoAVQqlh59/ApntiooNOiYeoGKZBb0me5zy/pzHEyqz2u0ZaEEVRxPBL362DioHui2pqLci6P7zii2cdUwC0ZD0hwFvG0AHN4s2LGt4qBblf7QRKhtSx86gixfYTtExgRE3GxqSaeYVkRWLlmXmRWunNhFUfujGLImmvw0/JvWuGFrgiWdBX81J5jomMDhpD6EYMhe6/5GVWsNM9cC9HlQpBEs6GyX5+NMrtk8CwMCrrmmXhBc5DwZtfTGkBEMmvGZLfd+KHHeYbTZSTAuVD7X5XHemeihGZA+0rX0IJSuGio0nD89vGJHMoG3AOyZhwvsuDm8EhFog6URVjA2d/Lp0KAbCmrkwNMvGbNIxgRE3H0bi3Ayqt3jMC3LoOWGUyQ7kGYeKRzCPqtngW2koRgxghm7kgJMKqqJjEgOq1IOM84eWw5+mqRW5SNfzNaGyXTuM5sZJpkMxDpy3HUGzMbBj6JjEcPHFTQFYPBRtgChLE9tLHTGihChdjQU/4YF18utkjD3xwkhYruiYRGAweKkY2WXRMIztLIiDyFLXBiCHqqwi2KF68+HyJIy4RdCa8DNPVG9DbrLIjUOLzsV7LFRFsPOKjaPwPaggeRMJKv1w7Dq2ZSWpu/4lUDXCQjUJdi+GCi8dFtVbFNpBljiBpaeh1FBoS3194pPwGIvqLeb9fzvwYl7eQmwg1UJDFX4KktrPhoETejYvbyH1OX14C8h9+ineJOCdEydobURJrfBQBTyYAhWDJaqEh8q06Y6YKzjSaTHJCRGVAdMdF4S5Lqhh89gxURk2jx0TlVnBfgV1YV6Dgv2ZZ6+iojIl2K8cXL12iLzcsxHBzg11dJ1hr4xNP9gLQwnhoqIe7FNDCSGvt0452OcMJYSMim6wLxhKCHsVf53BfnQAOpTYsmwoIWxUFIP9tKGEsFHpCfbrBzcOb8htutRQQtio9AT7s5zTkcyGFYYSQkelI9iPnzs8rCRwompDCaGj0hHsh89LbFRnKCF8VDzYk0DtdfZjQFVb860wlBA+quQzlqLbI0508MzhjePql1caSggbVRJ4Cm6POKUr1SRkDCWEjWrY8u0RTSVnKCFsVPFQP6rjw4le+KykoYSwUU1YaUV19NZEn/u8pKGE0FEVrLQe8fqE1BeOG70NB9X09ggx9xVYyb2trcv21wSp55qRIuAqzuqLdaicQAietjOr+MrBl9YhRQIVi79c02KPUp/5Yha/30YfSFR5V9cghYVq8U4S9nRYuWUCkIpZ6unGbfppG+poDVJYqEp3ktQMxUDd6BZ22tRUJ22o62uQQnNV6U6S2qEYt1hCIQn5QzAexsxfY0KtdKO8UlhZVbqTpG4oJs2KJAt5BZjakFjNZ9U2aJRXCm0y9uKdJDVDMYIUpJSYhczgftyGptrcUKioyneSVI6x20NYHonDiWdeajZ3tA1DCaE1Fkp3klQFe+xNDXUCqEHzqiVDCZFoV4GWB3sytHjz0xOoZn+TLn+tGUqIDKpVY+zODFAq179p01BCZFCtGmMPljyrUbuGEqKDqv7iqX8CKGSr1LqhhOigqr94msiPKigwlBAdVO1cPFVjKCFCqFq4eKrKUEKEUEGwp5O7bpJgUuGlk2aqXSwkzHxx37z4Q7GlP1uPV6GhyKEKQtcKxU3LsZcW60y4biB6gM7kD34WuGLoJs1C0TiFNeGEG5UaihyqF78Ca+JyND6n4ItFYDgFWB4v5gDhjwxunYchGZtvFUPzAgBaIUu++rUXVH88SqiK8asgsoWBEs8JxVq7ruuIW79jzykGrqLAEY5KMxvWi+FUw+zrR1effebtgyor0seZ3N8du5NuXzipG2OraFz5gRtPtgyLP8AK1uzw2tWDK28PVHG06R6uK7UWIVRuG/fBr22tg5Vb0EHltzTVcU1rmYQqbecqH2gda5mEKmvzlu9m1rpycHDw/MHBiom1ZFBtHuolNbSWQa5qJdQXJWOt2TRtc1C1FeolrbTWbJq2OahaDPVF1VtLcpq2EBVUrYZ6STXWkpqmPRERVK2H+qIqrbV6mvaJiKBSEOolVVhrxTTtedFApSjUF7XcWlekB7looFIW6iVxa33jhen84qYDgTRQqQz1RV3/5nR68VtGolIc6ou6MSF1rekbSaBSH+rzmrBqfMWCAiotoT4nwepbjQcfKKDSFeozAaujxuNaFFDpC/WpDp8Tod5s8IEAKq2hPtHxpHfcxFoEUOkN9VOSthY+Kt2hflqS1sJHpT3Ul0nGWvio9If6Uq22FjoqjFCv0AproaNCDvVF1VoLGxV+qJdUbS1sVCRCfVFV1sJGRSTUS1pqLWRUhEJ9UUushYyKVKiXVLYWLipyob6oRWvhoiIY6iXNWQsPFfz32TRDfVEzayH+3xDj6NvfweYgp8JamKjG4++6BP4vLRmBtXBRcZFPq6mOsVFlJqyNXQgZlUGkkFGZRAoXlVGkEFF9zzBSiKi+/wOzSCGi8l7EPveGQkM1HX4prWRFWGiopsMvpZWsCAsL1Wz4pbySFV1hrl918oz2sNVEWKhOhl/KK1mRFRKquTH18kpWZIWEan5MvbSSFVnhoCI+pr5cOKiMaByUhYPKhDH1U0JBRfZCaa1QUFG+UFotDFRGhjoOKiNDHQeVkaGOgsrMUEdBZWaoY6AyNNQxUBka6hioDA11BFSmhjoCKlNDXT8qY0NdPypjQ10/KmNDXTsqc0NdOypzQ103KoNDXTcqg0NdNyqDQ10zKpNDXTMqk0NdLyqjQ10vKqNDXS8qo0NdKyqzQ10Hqnwqs0NdOSpOiOXwr3iOfbKbSS2qOU6s4GawVKLKFzixKTlTpdpVp3+MlTpU+TJOJvtKIaqlnjI4r5Shyit+zLWVKlT52WOlDlUlq/+jWlQlp9xYVtKonrzvfIPd1pY/c1H98Ecy+nGe/+SnUluC8rpcNxfVYyMZvQ493r/fviu18ag2p8xF9fg9Gb1RDA/845//ktk6P5tZ9cRNGb1ZoPr3f/4rs/UZzSpFBfAsZpV8rP+sSayfxayS27BhY6E2q+5vtC8qUtYErcuq0a2bFy9gn3ljyaPq7O30Bt19WVQ15e+Vl3ia3TbOWvKottj2ZdbbkUZV3be59/Kdm7dGxlmrgas6u5zXI7I7rmKV5z8fje6+esc8azXIqkGfsf6eNKq8qvw9/AA0OZpYq9Ptdnus1+1uGYKqu8c6ux3pPVeUP/jnoYsPvtbMWtv8a2I7PUxSTVBt7/X6sqlew0rofENrbfHCvyf/NSGj6myzy832vaQMzrU+m1nr0lYXmVQDVINu873XXzJtYq3ubhND46Lq9ftrREU+1xRd1veTttZWbg6qtTWb3rH8ZTlrDfZ3+siksNdbF1ptrd6A9S6Zk1UqtcJa+wP+sCvdpFMkGqhYnbV6/Zyn5GB7d416pU2RQdW8raVbhFCxdZrxGkULFWlrUUPF6FqLICqq1iKJipG0FlVUBK1FFxWjZi3SqGhZizgqRsha9FGRsZYJqBgNaxmCioK1jEHF0K1lEipka5mFimFayzhUeNYyEBVDspaZqJZY6yHlhzQVFStb6xfKj2cwqgVr/fJXyo9mNCp2Yq1fqz+U6aim1vqN+gOZj4qBtX6r4Sj/A1vNkGLJhFaYAAAAJXRFWHRkYXRlOmNyZWF0ZQAyMDI0LTAxLTMwVDExOjA2OjQ0KzAwOjAw6cdMcwAAACV0RVh0ZGF0ZTptb2RpZnkAMjAyNC0wMS0zMFQxMTowNjo0NCswMDowMJia9M8AAAAodEVYdGRhdGU6dGltZXN0YW1wADIwMjQtMDEtMzBUMTE6MDY6NDQrMDA6MDDPj9UQAAAAE3RFWHRwZGY6VmVyc2lvbgBQREYtMS41UzZawQAAAEp0RVh0c2lnbmF0dXJlAGNkYzBhMWNlODJmMWJmNzJjOGQ3NmJhMjA1OTdhNzc1ZTU1ODk1MDliZDIzN2ZlODYxNDY5YzVjZjAyMTBiYWUjKeYMAAAAAElFTkSuQmCC", + "text/plain": [ + "" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "fd = FeynmanDiagram().add(\n", + " v1 := Vertex().with_xy(7,0),\n", + " vh := Vertex().with_xy(6,2.5),\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=22,label=\"$Z/\\\\gamma$\").connect(v1,v2).with_color('red'),\n", + " Propagator(pdgid=-2).connect(v1,vh).with_color('red'),\n", + " Propagator(pdgid=-2).connect(vh,v3).with_color('red'),\n", + " Propagator(pdgid=2).connect(v1,v4).with_color('red'),\n", + " Leg(pdgid=2212,target=v3,sense='incoming').with_xy(0,5),\n", + " Leg(pdgid=2212,target=v4).with_incoming().with_xy(0,-5).with_style(\"label-side: right\"),\n", + "\n", + " \n", + " Leg(pdgid=21,target=vh,sense='outgoing').with_xy(15,4).with_color('red'),\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", + " Leg(pdgid=11,sense='outgoing').with_target(v2).with_xy(15,2).with_color('blue'),\n", + " Leg(pdgid=-11).with_target(v2).with_outgoing().with_xy(15,-2).with_color('blue'),\n", + ")\n", + "scale_positions(fd,0.5)\n", + "fd.render(auto_position=False,render=\"tikz\",file=\"Zj\")" ] }, { "cell_type": "code", - "execution_count": 6, + "execution_count": 3, "id": "27b1d0cd-dc52-4d75-970e-336e95c41e32", "metadata": {}, "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlEAAAGuCAMAAACHjqO9AAAAIGNIUk0AAHomAACAhAAA+gAAAIDoAAB1MAAA6mAAADqYAAAXcJy6UTwAAAEvUExURf///wAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA/wAA/wAA/wAA/wAA/wAAAAAA/wAA/wAA/wAA/wAA/wAA/wAA/66urgAA/wAA/2dnZwAA/wAA/wAA/wAA/wAA/wAA/wAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA/AAAOAAAAAAAdAACDAADKAAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAAAAAAAAAAAAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAEdHR1VVVf8AAP8AAP8AAP8AAAAAAAAA/////4iIiAD/AP8AAOvr428AAABfdFJOUwAAEWYziHfdVcy7mSIRd5kzIkRmRIjdu1Xux83kqKppqOZ2sPbJqu7AoJiCgJDe2rCnd0QRM2aIIrtVqu6ZzFB2kMHkzfHEEUSIqt27ZjOZVcwi7ndQm6h2es2Qx+/mZirjrgAAAAFiS0dEAIgFHUgAAAAJcEhZcwAAAGQAAABkAA+Wxd0AAAAHdElNRQfnCwMJHhRjus9IAAAPVklEQVR42u3diXbbxhWA4QxIcFccw5KztE3SJVL3JYqtbG26S7KtOHaW7qXS5P2foYMLShYlkcRyhzMD/N8RTdtHQwzN/wAgQMEvvQQAAICajFNJr+92AQiL66KS1PQHvp8ktsh1Ub3UmKHrJzEvuF4MynBdVGqLcvtS5ynlS5gbogpB7EUtKppf3GjKN9dFjWxRY3cPf1kRSYXCVVHpeDgZT2xO09lg5Gz2xfbOtCmpnZfvvGLu3vE9jfocFZX1Tc8U7/H6ibPJX1TUpqTumTu75t6e72nU52yr189m7tZNhXkbk9rZuW/yqKLlrKipyTK3U79a0XJSbpfrWL5+uu97Eg24KiqZmCx1O/WVNcW9ktrbM3dfNa/5nkZtrooaDEx/4nTm85YmtXtn9+VXX97xPY3anK2j7G3mdObzdWspp0t2bOeu3ZnyPYn6nJ8pdmbdblS8Rb0S8dqpYIvK0hjNl7dxy7d4i9o9f/3OK74n0Ygt6o2vY7R+oxdtUrvn1qv3fE+jgZiLauOuuRR1fv763l3fM6kr3qLWJ/VqpL5zvvBd3zOpK96i1u1HRb7VO78f8f65Lep7/4tRi/ej3oz38KaJ+OjBfG1Se3f2XtuN8T3TbsR7UIVoi9p0PGpnd3f3nr1FllW8W7sLERdV4nhUpFlFLd6iXpS04SQMWW1VzEVVOQlDVtsScVHVP3FHVlsQc1EXSV1fQ60/dEBWbkVe1G0bvTLHosjKmaiLMvObPwxT4ZQeWbkQd1E3f76q8sFyslIWe1FmccGDi59Tr/cgZKUn+qIW19Ewza+lQVYqWlCU0Do1TFZNtaUo1Z+7IasGKGoVsqqHotYiq8ooajOyqoKiSiKrkiiqCrLajKIqI6u1KKoeslqFohogq1tQVFNktYyiVJDVJYrSQ1Y5ilLW+awoyoUuZ0VRzjTMKtZL5VOUW/WzivVS+RS1BbWyivVS+RS1LZWzivRS+RS1VVWyivRS+RS1fSWzivRS+RTlyaqsXnuRUJyXyqcon27J6v79+HbGl1CUd0tZ7Zyfn+/Ft2K6gqLCcJHVW3Kp4JhXUxQVEJvV28XVzCNeTVFUWN5cXCD/+z+I9ZwgRYXlfnGF/B++/Vasp5opKig75/ff3NtdbPLi/AQDRQXl5vGn6LKiqBjElBVFRSOSrCgqLuFnRVERCjoriopVqFlRVNQCzIqi4hdWVhTVEsFkRVFtEkJWFNU6nrOiqHbylxVFtZiXrCiq7badFUV1whazoqju2E5WFNUxzrOiqC5ymRVFdZajrCiq2/SzoijoZkVRKGhlRVG4QiErisJ1zbKiKNyqdlYUhdXqZEVR2KBiVhSFMspnRVEorVRWFIVqNmVFUahhTVYUhbpuz4qi0MiNrCgKzV3NiqKgxGb18p29PYqCpp1dioIuioIuioIuioIuioIuioIuioIuioIuioIuioIuioIuioIuioIuioIuioIuioIuioIuioIuioIuioIuioIuioKuiIvqTQpp/geKCkXERS2hqFBQFHTFXFQymPbNbCq/p6hQxFxUz0wz00vNaF4Y+Z4QTNxFJcnQ5FGZoQQ19D0f5GIuKl8/mbyjVIoa+J4OclEXlaZmNjEjM5OiEt/TQS7qorJpNpgMbElTG9TU92wgoi7KJDO7M2XvM1tU3/dkIOIu6tJwPvY9BRRaUlRv3vM9BRRaUlTyI/bLA9GSosw7vieAhbYUxVmYUFAUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUdFEUFPWzEUVBh40p3U9TikJz/SzLBvY2y/9AUWhiKSZBUajpZkyColDdipgERaGSdTEJikJZG2MSFNV6gzSXNXqMcjEJisJ6FWISFNViaW/U6L+YrxqToKj22rc5DWqOrRWToKjW6h/YKJIa42rHJCiqtbLq/yYNYxIU1Vr9/N+k9Fs8jZgERbVXOsp6/RLfpxaToKgWm23ci9KNSVBUZ6TLKywHMQmK6oyD+Xw+kaxcxSQoqjMm84Ufu4pJUFRn/GQR1E/dLoaiukA2cz+TnsZl3v01QVEt92KfaZYHlTpfIEVVkZQ6vhOK6zvgB+5XUIaiKklS06976nW7bn03N3G/gjIUVc5IPrSW9uxLMvT9TDdxemigBIoq9ei9xCTjXmqLmvt+pmv4jklQVAmz/HNr49QEXFQQMQmKKmFkb9NpvvGzYfl+pjeEE5OgqHKm+/LrbDDy/UyvCiwmQVGl9MbFafx+jc9EuhFiTIKiypCggjkUFWxMgqJKyA5sTbMQtndhxyQoarP+MF89TT0/wwhiEhS10exgP03ToceiYolJUFTYoopJUFSw4otJUFSIIo1JtKao/nAoP+PfH8t+tOkNi/vBwUSOIaXFfbJ/IGfgk0lxfzlsuDystzzMpPOlYfa+V2H4mqVmV4eZuGMSbSnqnWE/yV8V+76sP7SvRjZOsvwgUm+c9PL1VzqV8yj5L9P8RZ30krF9Ue0vydi+qEk+TIYn9suYkR2WDx/s51/F8HzY/iD/MiYfNpJh15cqw/OljothaYWl/vwX6X4acUyiLUX9ciQvjLxQNgd5nW0L8kLZ11PqmvRMXlf+PfnrnH9P/jrLsJ7kIMPs90ge9nuK4amRPOz35HVdDJvZ+7wuGT7SWGrvV7/+zbuHh+/5/rdspi1FyX5Ub7ExSRf308UWcFzcJ+OJnEzpDy/ux9PFsPTqcDtsdGX4xbAXw4thV4aXWurB/qal2owfHB4ePjyMOatWFWVmiw3Gxf3FabiL+4urACx2WpK+qTQsWx5WdXjJYZFn1a6i2iTWrCgqaK6zOtJ/SIoKn8OsKGqlNhcl3GRFUSu1viihmtX7R0dHHxwdPVSeI0VFRzMr1lErdago0TCro4cfyrqJolbqWlGidlYf2ZzezX9DUSt1sihRI6sHHx8evuvqgARFtUK1rA4/cDgVimqP0lk9yIs6dDQLimqZUlkdfXj48IGjCVBUG23M6n13ZwspqrVuZHX47m9/95HzxVJUuxVZffL7P3z0u29yDg4XXENRHfDgkz9+s+Bqf/wFiuqGDz8uivrE+eetKKoj3vsoD+pP7j/GR1Gdka+mZMfc7cf4KKo77GrqxY65s6woqks+vHZY00VWFNV5yllRFHJ6WVEULqlkRVFY1jQrisItGmRFUVilXlYUhbUqZ0VR2KxKVhSFkkpmRVGoYnNWFIXK1mZFUahnVVYUhQZuyYqi0NRyVhQFFZdZURT0PHh49AFFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQRdFQVcbipplWfZne5v5nghMO4rK5oXM90Rg2lGUGUpQY5MMpn0zm/qeTre1oqiBFNUzPTPNTC/1PZ1ua0VRyYEN6iAxSTI0eVTwqBVFmaktKt/Y5eunoe/JdFw7iurbovJ3emlqZhMz8j2dTmtHUWY8lwNS2TQbTAaJ79l0WkuK6tn98lwysztTvifTbS0pyvzF9wSw0Jai/up7AlhoS1Gc1wsFRUEXRUEXRUEXRUEXRUEXRUEXRUEXRUEXRUEXRUEXRUEXRUEXRUEXRUEXRUEXRUEXRUEXRUEXRUEXRUHVMUVBy/HJ6aPHTygKCo5PTk7O7O1Ttnpo6kpMgqJQ2/WYBEWhjltjEhSFilbHJCgK5W2ISVAUSikTk6AobFI6JkFRWKNaTIKicLsaMQmKwg11YxIUhasaxSQoCgvNYxIUBbWYBEV1nGZMgqK6Sz0mQVGd5CYmQVFd4zAmQVEd4jomQVHdsJWYBEW13vZiEhTVZluOSVBUSE5P1R7KR0yCooLyWONBvMUkKCokJ08bPoDfmARFheT0rP7YAGISFBWM088+e/a81shQYhIUFYp8/VR9NyqomARFBeL551V3o8KLSVBUIL54VGE3KtCYBEUF4uzUmGfPP9v4fSHHJCgqEMdPT86enpys/ZbQYxIUFQwbysq3enHEJCgqdBHFJCgqKMePH52eHF/+KbKYBEWF5fG3ls3qyy8jjElQVFgefbvwVYQxCYoKyfHJ34qenkWak6GoYCz2mf4uQel9TGr7KMq/qzvgx3GvoAxFeXbz3dzjqFdQhqL8WXFo4DTqFZShKC+iPM5UFkVtV6tjEhS1Ne2PSVDUNnQkJkFRjnUpJkFR7nQuJkFRTnQzJkFR2jock6AoRV2PSVCUDmK6QFGNEdMSimqCmG6iqJqIaQWKqo6Y1qGoSohpI4oqi5jKoagSiKkCilqPmKqiqJWIqRbnRSW9/jaeh25RxFSf66KS1PQHW3geakURU0Oui+qlxgwdP4f5hYaPQ0waXBeV2qKavtLr2ZDyBcyL39YUTEz5RITvidQWe1FFRC+Sqr6sYGISz08XfE+kNtdFjWxRY2ePnge0nFSlfMOKqSVcFZWOh5PxxOY0nQ1GzmY/v9zWVU0q2JiOv3j6qXne9H/z8MhRUVnf9EzxHq+fuJr7LWuoMkkFG5M4M09P5MLBsXK21etnM3frpsLlyql0UmHHVEzx+LHJo4qWs6KmJsvcTn3+oqISSUUQUyFfP6n8p2ieuCoqmZgsdTv1pWMG8+t/cUU0MQn7Nu/5E7P5UvmhclXUYGD6bj9hMl+1hrqaVFwxiZOnJ188+eK4+QN54mwdZW8zpzNfvdErbhHGVDh+bnemfE+ivng/e3C1nus5zU8fPX5yehZdTLnjs2e+p9CILeof/4zRfM1+lP2KMibr+dPPv/3W9yQasUW98XWM1m/03J5KdObsiVwq2Pc0Gom5qHW75qfx+de//7O4Or7vmTQSb1HrkzqJz5df/XdRlO+ZNBJvUev3o3yv++s5edqKrV6se+Yt3I/KP83yOP6ifE+hpk3Ho6J19sT3DBqJuqjVx6N8T67DIi5qzX4URfkTb1ErP8zCKsqrmItioxeiiItamRRB+RRzUSuSoiivoi7q+s9WGfbK/Yu7qFuTgleRF6XwE6DQFXtRS4ejCCoA8ReldyUNaPg/wchW3EDy5LgAAAAldEVYdGRhdGU6Y3JlYXRlADIwMjMtMTEtMDNUMDk6MzA6MjArMDA6MDDy+6/KAAAAJXRFWHRkYXRlOm1vZGlmeQAyMDIzLTExLTAzVDA5OjMwOjIwKzAwOjAwg6YXdgAAACh0RVh0ZGF0ZTp0aW1lc3RhbXAAMjAyMy0xMS0wM1QwOTozMDoyMCswMDowMNSzNqkAAAATdEVYdHBkZjpWZXJzaW9uAFBERi0xLjVTNlrBAAAASnRFWHRzaWduYXR1cmUAMWQ4NTA3M2MxNmJjNGIyMzdmYjgyZjA1NzgxYzhhOWM2OWU5NmFmNWZiYTFlMzk2MzUxMDQ5NjJmZWE3NzA4NgD+JBgAAAAASUVORK5CYII=", + "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+Wxd0AAAAHdElNRQfoAR4KLhSnL2LUAAAPQ0lEQVR42u3dCV/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/dotR9aY92QAAACV0RVh0ZGF0ZTpjcmVhdGUAMjAyNC0wMS0zMFQxMDo0NjoyMCswMDowMLaoy+AAAAAldEVYdGRhdGU6bW9kaWZ5ADIwMjQtMDEtMzBUMTA6NDY6MjArMDA6MDDH9XNcAAAAKHRFWHRkYXRlOnRpbWVzdGFtcAAyMDI0LTAxLTMwVDEwOjQ2OjIwKzAwOjAwkOBSgwAAABN0RVh0cGRmOlZlcnNpb24AUERGLTEuNVM2WsEAAABKdEVYdHNpZ25hdHVyZQA3NzBiNjU4MGI0MjY4MmQwMTQ3N2Q3MzcyNGMwNjM5MWRiYTcyYjc5ZWVjNTE0MWUxYWRjM2M5YjQwZjg2Yzk3s3D9igAAAABJRU5ErkJggg==", "text/plain": [ - "" + "" ] }, "metadata": {}, @@ -146,7 +221,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.11.5" + "version": "3.11.7" } }, "nbformat": 4, diff --git a/pyfeyn2/render/latex/tikzfeynman.py b/pyfeyn2/render/latex/tikzfeynman.py index 5a41b63c3..02ca0e534 100644 --- a/pyfeyn2/render/latex/tikzfeynman.py +++ b/pyfeyn2/render/latex/tikzfeynman.py @@ -43,6 +43,7 @@ "phantom": ["draw=none"], "line": ["plain"], "plain": ["plain"], + "remnant": ["plain", "double=none, double distance=%(double_distance)s"], "baryon": ["fermion", "double=none, double distance=%(double_distance)s"], "anti baryon": ["anti fermion", "double=none, double distance=%(double_distance)s"], "meson": ["double=none,double distance =%(double_distance)s"],