validation_demo.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import tensorflow as tf\n",
    "#tf.compat.v1.disable_eager_execution()\n",
    "\n",
    "import numpy as np\n",
    "from matplotlib import pyplot as plt\n",
    "import os\n",
    "\n",
    "from junipr import *"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "dataset = get_dataset('./data/tfrecord/JuniprJets.tfrecord', \n",
    "                      batch_size=500,\n",
    "                      max_padding=True, #Needed for validation code to work\n",
    "                      repeat=False, \n",
    "                      shuffle=0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Using SimpleRNN Network\n"
     ]
    }
   ],
   "source": [
    "junipr = JUNIPR()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "junipr.load_model('./saved_models/Junipr_demo_model.h5')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "PYTHIA_val = junipr.validate(dataset, predict=False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "JUNIPR_val = junipr.validate(dataset, predict=True)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# P_end"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY8AAAERCAYAAACD9ivUAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzt3XmcFdWZ//HPF1BQEYmKS4wGNzbjjkpmEm1QiBqTGRMXXEJwwyjjrpEx6q81k6AhbpO4BE3ECJqooxONuBAQYxxRcUlco0nAxAUExSUqLvD8/qhqvFxu367q7nv73u7v+/XqF7dOnap66pb201WnzjmKCMzMzPLo1tEBmJlZ/XHyMDOz3Jw8zMwsNycPMzPLzcnDzMxyc/IwM7Pcqpo8JO0i6S+SxrZQ7zBJj0maK+kiSapSiGZmlkHVkoek/YFTgLdbqPcF4CLgK8CuwE7A8RUP0MzMMqvmncejEXEo8G4L9Y4CpkfE4ohYDvwC+E7FozMzs8yqljwi4uWMVXcBni9YfhbYRtIa7R+VmZm1Ro+ODqCEDVn50dZbgID1gX8UVpQ0DhgHsNZaa+08aNCgasVoZp3MwneW5qq/YZ9eFYqkuh577LHFEdEv73a1mDwASg24tUqjeURMBiYDDB06NObOnVvpuMysk7pkxgu56p8yckCFIqkuSS+1ZrtafFX3daBvwXJfkmSyqGPCMTOzYrWYPB4FBhYsDwGeiYgPOigeMzMr0uHJQ9L6kh6QtG5adA2wr6T1JHUDxgJXdViAZma2imr289hZ0mxgB2CCpFvTVWsAg4A1ASLiaeB04F7gYeBJ4IpqxWlmZi2rWoN5RDwGNJQo/wfQr6hsGjCtOpGZmVleHf7YyszM6o+Th5mZ5ebkYWZmudVqJ8F2l6UDUKlOPzNmzOCMM87gj3/8I7vvvjvLli1jyZIlHHvssZx44omZj//RRx8xatQo7r//fubNm0f//v1XqTNx4kQ++eQTzjnnnMz7NTPrCF0mebTWyJEjufTSSxk+fDgzZ86kR48ePPPMM+y4444MHjyYkSNHNrutpBWJYvXVV2f27NmUG13+1FNPJaJU53ozs9rix1atsM0227Dtttty9913t+t+e/bsSa9enWO8HDPr3Jw8Wumdd97h9ttvZ4011uDoo48GYNKkSay//vqcdtpp7LPPPgCMHj2ahoYGXnnllRXb3nnnney7775sueWW3HbbbUDyeGzQoEE0NDSsqHf++eczYsQIRowYwX777cerr75avRM0MyvDyaMVZs+ezbx585g6dSrLly+nsbERgJNOOondd9+diy66iLvuuguAX/3qV8yePZtNNtlkxfbvvfce06dPZ+LEiUyYMAFIHo81fW7ymc98hpkzZzJr1iwOOOAAzjzzzOqcoJlZC9zmkcOee+7JsmXL6N69OzfffDO77bYbI0eOZOrUqUyYMIE777xzxR1HOXvvvTcA22+/PfPmzWu23qabbsrw4cNZvnw577zzDh999FG7nYuZWVs4eeTQ1GBeaMyYMZx33nlMmDCBm266iSuvvLLF/fTp0wdI2jg+/vjjknVefPFFDjroIB588EF22WUXZs+ezdixY9t8DmZm7cGPrdro61//Oq+++ir33HMP3bp1o2/fvi1vlMETTzxBnz592GWXXQCaTTJmZh3ByaONevXqxYEHHsgRRxzBwQcfvNK63r178/777zN16lRuueWWXPvdaqutWLJkCS+8kPRPae83u8zM2qLLPLZq7axfTZ0EIWnzOP/889ljjz1WqjNmzBhuvfXWVdo7xo8fz+jRo+nTpw+33HILo0aNApI3sO655x5Gjx4NwKhRozjjjDO44IILWLBgASeccAI/+clPOOussxg1ahTbb789G220EQsWLGDMmDH88pe/bNW5mJm1F3WWTmkdOQ3tc889x09/+lMuv/zyDjm+mbVdF56G9rGIGJp3Oz+2aoMbb7yRTz75hClTpvDtb3+7o8MxM6saJ482ePzxx9lxxx1ZvHgxu+66a0eHY2ZWNV2mzaMSJk2axKRJkzo6DDOzqvOdh5mZ5ebkYWZmuTl5mJlZbk4eZmaWm5OHmZnl1nWSx30TW/4pYdasWTQ0NCCJ7bffngsuuID+/fuzzjrrMGbMmBX1Lr74Yvr3788OO+zAjBkzuPLKKxk0aBB9+/ZdaR6ORx55hGHDhrHRRhtx7LHHAknP9V69ejFw4ECmTJnCU089xbBhw5DEbrvtxpe//GW22WYbvve977F8+fJV1jc0NDBw4EC+9rWvsXDhwsp+j2ZmdKUe5s0kh5UM/89mV0lixowZ7LXXXjQ2NvK73/2OP/zhDyvVaWxspKGhYcWETlOmTOHoo49m33335fbbb19Rb/78+TQ2NjJlypQVZf379+fss89eMbHU/Pnz2XzzzXnxxRfZaqutWLhwIYMHD+bCCy/kmGOOWWX90qVL2XXXXdlpp51W2q+ZZeMe5vl0nTuPDnLiiSdyzz33MG3atDbtZ8MNN2T48OHNDpDYq1cvRo0aRUcN0WJmXYuTR4Vtt912nHvuuZx00km8/vrrbdrXxx9/zGqrrdbs+k8++YTPfe5zbTqGmVkWTh5VcOaZZ7L55pszfvz4Vu/jueeeY+bMmXzzm98suf6NN97g5Zdf5qKLLmr1MczMsvLwJO1M0iplPXr04LrrrmOnnXbi1ltvZaeddsq8v8MOO4yePXuybNkyLrnkEg488MBV1nfv3p3HH3+cU045hW222abN52Bm1hInj1bo2bNnyfnEP/zwQ3r27FlymyFDhnDeeedx/PHHM3369MzHmjZtGltttVWL66+44grGjx/PIYccwnbbbZd5/2ZmreHHVq2w6aabMn/+/FXK582bx6abbtrsdqeffjr9+/fnlFNOafeYjjvuOL7whS9w2WWXtfu+zcyKOXm0wqhRo3jvvfe49tprV5Q988wzLFy4kE022aTZ7bp37851113HI4880u4xSWLChAnccMMNLF68uN33b2ZWqOs8tirThyOrpvaMDTbYgNtuu42zzjqLiy++mHXWWYd1112XyZMnr6h75ZVXrrgL+Oijjxg3bhwAAwcO5Ac/+AF/+tOfVtTdc889WbBgAZMmTaJHjx7svPPOHHPMMUDSpnHqqaeuND/6U089tdL6iRMnMmLECEaPHs25555LQ0MDjY2NHHDAAW0+ZzOzUrpOJ8FWePfdd1l77bV577336N27N4899liuxm4zqx/uJJhP1R5bSeolaYqkOZLmShrVTL21JV0v6VFJj6Tb9K5WnIX22WcfAJ544gl69+7NgAGd4z8WM7O2quZjq0aSO51hkgYAcyQNjojiwZjOAT4PDAMCmAmcDUyoYqwAfPazn2XYsGEsXbqUyZMn07t3h+QwM7OaU5XkIakbcDTwDYCIeEHSE8DhQHGvtm2ARyJiWbrtI8CO1Yiz2E033dQRhzUzq3nVemy1BbAe8HxB2bNAqeds04ERktaStCYwAni48iGamVlW1UoeG6b/vl1Q9hawQXHFiLgcmA38DZgPPA6cX2qnksal7SdzFy1a1J7xmplZGdXu51H8atcqY3lI+h6wPUm7x2Ykdy1jS+4sYnJEDI2Iof369WvnUM3MrDnVSh5Nw8n2LSjrW1Be6ETgZxGxNCKWAlcB369wfGZmlkO1ksdfgTeBgQVlQ4BHS9RdHfi4YPljYO3KhWZmZnlVJXlExHLgauBIAElbAzsA0yQNljRTUve0+u+Ag5UCDgXuq0acZmaWTTXbPBoBSZoD3AgcEhELgHWAQUDTLEfjSdpC5pC8ZdUDGFfFOM3MrAVV6ySYtl+MLVE+B9ikYPl14ODiemZmVjs8qq6ZmeXm5GFmZrk5eZiZWW5OHmZmlpuTh5mZ5dam5CEp9wQiZmZW/5p9VVfS7hm2vxTw1HpmZl1MuX4e9wAL0s8BfA74gGSYkXWB3sBLFY3OzMxqUrnkcUdEHAQg6T+BpyPijqaVkv4N2KXC8ZmZWQ1qts2jKXGkRhYmjnT9b4AvVSowMzOrXVkbzNeTtGthgaQvkjy+MjOzLibr2FYTgJmSFgCLgX4kswAeVHYrMzPrlDIlj4i4S9LmwFeBjYDXgDsj4o1KBmdmZrUp86i6EbEYuK6wTNLJEXFpu0dlZmY1LVPykLQ+cBSwOclMf032JunrYWZmXUjWO4//Bd4G/gQsLSj/sN0jMjOzmpc1eXSPiK8WF0p6vJ3jMTOzOpD1Vd37JW1WotxDk5iZdUFZ7zyGAeMl/Rl4p6B8B+C8do/KzMxqWtbksSHwH0VlAjZu33DMzKweZE0e3y0engRAkgdGNDPrgrJ2ErxD0lrAfsA6JP09BkTEfZUMzszMalPWfh5fBG4HXiHp5zEVuEjStIi4ruzGZmYd5JIZL3R0CJ1W1retLgBGRMQOwMKIeB/YFziyYpGZmVnNypo8IiKeavqcFnwCLK9IVGZmVtOyJo8PJX1bUvemAkn7s3JvczMz6yKyvm11HPBb4GcAkt4B/g58rUJxmZlZDcv6ttXfJG0L7EYyl/k/gIcjwo+tzKxTGPb3yTm3+HFF4qgXeYZkXwb8XwVjMTOzOpG1zcPMzGwFJw8zM8st82MrM7N6k78dw7LKfOchqY+kwyVNSP/tk+dAknpJmiJpjqS5kkaVqfsvkmZK+r2kZyWdkOdYZmZWWVmHJxkK3AUsARYB/YCLJe0dEVknhGoEFBHDJA0A5kgaHBELi461OcnUtvtGxGJJ2+Ce7GZmNSXrY6tJwP4R8YemAkn/ClwEDG9pY0ndgKOBbwBExAuSngAOT/dR6FTg2ohYnNZ9BjgtY5xmVkGVGivqlJEDKrJfq5ysj61UmDgAIuJBkjk9stgCWA94vqDsWWBoibp7AqtLmi7pQUk/ltQr43HMzKwK8rR5fLlo+Uuk41xlsGH679sFZW8BG5So2x/4DjAGaACGkDzGKhXTuLT9ZO6iRYsyhmJmZm2V9bHVGcB0SW/zaZvH2sDeOY9XnGxK3bn0BG5oemwl6VLgdknHF/doj4jJwGSAoUOHZk1kZmbWRlmHJ3lU0pYkk0E1DU/y24h4N+NxXk//7dvM50JLgMJG9JdJEsr6zdQ3M7Mqy/q21RERcS1wQ0HZVEm/jYhfZdjFX4E3gYF8mgCGANNL1H2SlR9n9QM+At7IEquZmVVe1sdW3wKuLSobB8wCWkweEbFc0tUkr9w+IGlrYAfgMEmDgZ8Co9Lxs64GzpV0cTrp1JHA1HSdmXV1903s6AiMFpKHpHPTj/0LPjdZk2Q+86wagaskzUmPe0hELJDUHxgErAYsi4hfS9qCpB/Iu8CfgZNzHMfMzCqspTuPzdN/1yj43ORd4IisB4qIpcDYEuVzgE2KyiYC/vPCzKxGlU0eEXEEgKQTI+K/qxOSmXU1eTofnuIR+WpCpn4eThxmZlbIQ7KbmVluTh5mZpZb1n4e26QDFJqZZZJnLo05m42rYCRWCVnvPG6XtEdFIzEzs7qR9b2Ft4CdJZ0OPAVcHxHPVS4sM+tKcs34t8V6lQvEMsuaPL6SDlR4saTtgCMlfQG4G7gxIjzmlJlZF5I1eQwAFktaPf08iGQSqB7ArpJWAy6MiMcqE6aZmdWSrMnjakkPAAcBrwBTge9ExCsAktYlmaZ2t4pEaWZmNSVr8tgUWAqMbObuYhiwUbtFZWY1KVfbhHVqWd+2uiAiTi5MHJK+Ialp6PSH8F2HmVmXkTV57FWibA3gFoCIWBIRC9otKjMzq2ktDcm+Wfqxl6RNWXna2MfJNyS7mZl1Ei21eczn03nHXypa9y5weXsHZGZmta+lIdm7AUiaHhH7VickM7PmPfQ3z0hdCzK9bdVc4iiY29zM6pWndbVWaDZ5SNo9In6ffi6egrbJWFad29zMzDq5cnceZ0p6PCL+CRxHMhRJsV6VCcvMzGpZs8kjIr5asHhRRPy4uE46UKKZmXUxWaehXSVxpFZrx1jMzKxOlGvz+EWG7fcG3NpmZtbFlGvzGA5MaWH7D9svFDNrF357yqqgXPL4YURcXW5jSa+2czxm1kbuB2HV0GybR0uJI1VqzCszM+vkyrV5TAZOjoj3Jf2tVBVgw4pFZtbJXTLjhcx1Txk5oIKRmOVX7rHVfcAH6ee3gZOL1gu4pBJBmZlZbSvXz+PGgsVvRcTTxXUkfasiUZnZSvLcpQyrYBxmTbKObfW0pL2A0cDGwGvAjRExs5LBmVnCM/hZrcnUSVDSmcAv08WnSR5ZTZX03UoFZmZmtSvrHOaHA9tExJKmAknrAfcDP6pEYGZmVruyTkP7cmHiAIiIN0geX5mZWRdT7lXdzQoWp0v6L+BmYAmwLkn7h9s8zMy6oHKPreaTTEFbOG/5WUV1Arggy4Ek9QKuAgalxz0rIu4tU3814FlgWkQ0ZjmGmZlVR7nkcX9EDC+3saT7chyrEVBEDJM0AJgjaXBELGym/jhggxz7N6srfoPK6lm5No9RGbb/QZaDSOoGHA38HCAiXgCeIGmIL1W/N3AQcHuW/ZuZWXWV6yT4ceGypC2BHYE1C4onAEMyHGcLYD3g+YKyZ4GhzdQ/DfgJsF+GfZuZWZVlelVX0jHAKcBGwB9JJoHaDig15lUpTWNgvV1Q9hYlEo+kfsDuEXGepLLJQ9I4ksdbbLbZZuWqmlWFe4JbV5H1Vd2xJMnijxExPCK+BAwE5uY8XhQtq0Sdc8j4OCwiJkfE0IgY2q9fv5yhmJlZa2VNHv+MiE+A1ZsKIuI1YHDG7V9P/+1bUNa3oBwASVsAm0fErIz7NTOzDpC1h/nHknYBnpZ0K/A7YDeyz2H+V+BNkruVpoQxBJheVG8P4LOSZqfLg4ClkhqAIyJiXsbjmXUIv0FlXUXWO4/vkjximgC8A4wnuXPINKpuRCwHrgaOBJC0NbADME3SYEkzJXWPiGsjYueIaIiIBuBuYEq67MRhZlYjso6q+2zB4thWHqsRuErSnPS4h0TEAkn9Se4wVgOWAUhaHbiXT+88touIb7TyuGZt4kZws1VlfdtKwFEUDckO/CIiihvBS4qIpZRIPBExB9ikqOwjoCHLfs3MrPqytnlcBIwAbgHeANYHTiJptzitMqGZmVmtypo89gR2jYgPmwok/Rh4uCJRmZlZTcuaPF4pTBwAEfGBpAUViMmspvgNKrNVlRuSffeCxd9Lup7ksVXTkOwHAw9VNjyzynAjuFnblLvzuIekYbywF/iXiur8G/D/2jsoMzOrbeWSx/SI+Ga5jSX9TzvHY2ZmdaDZToKlEocSG6Wv7pasY2ZmnV+mHuaS1pJ0DfAB8ArwgaTJktaqaHRmZlaTsg5PchnJiLgjSHp975kuX1ahuMzMrIZlfVV3QEQUvn31IvCgpN9XICYzM6txWZPHapJ6FnUS7Jlje7Oa4r4bZm2T9Zf/rcBcSb8GFgH9SOYYn1KhuMxyydNvA9x3w6ytso6qO0nSa8AY4HPAP4CJEXFjJYMzM7PalHVU3XOBlyNiVIXjMTOzOpD1baujgAcqGYiZmdWPrMljLvCX4kJJJ7dvOGZmVg+yNpg/A9wl6U7g7YLyccCl7R6VmZnVtKzJ4wTgSaB4KtiN2jccs9bxq7dm1ZU1edwQEeOLCyVd2M7xmJlZHcj6qu54WDGX+frA4kicWcngrGvznBtmtSvrwIjrSppGMjDiApKBEadKWrei0ZmZWU3K+rbVdcA/gb2ALwCjSBLJLysUl5mZ1bCsbR4bRsTXCpafI5ma9pEKxGRmZjUu653HK5LWLCxI5/J4qWD5iPYMzMzMalfWO4/FwJNpP48lwLokj65mpkOXAIwFrm33CK3L8uu3ZrUra/LYD7gb6Jv+ADwM9E5/AHq1b2hmZlarsiaPKyLi++UqSDqnHeIxM7M6kKnNo6XEkbWOmZl1Dp4J0KrKHf/MOoesb1uZmZmt4ORhZma5OXmYmVluVWvzkNQLuAoYlB73rIi4t0S9r5EMAd+D5LXgKRHx39WK0yrLfTfMOodqNpg3AoqIYZIGAHMkDY6IhUX1LgEOjIgnJG0EPCPppYj4TRVjNTOzMqry2EpSN+Bo4OcAEfEC8ARweInql0fEE2m9BcB9JL3ZzcysRlSrzWMLYD3g+YKyZ4GhxRUj4pKiol7AosqFZmZmeVUreWyY/ls4//lbwAblNpLUB9iFZsbMkjRO0lxJcxctcn4xM6uWar9tFUXLaqH+hcD5EfFSqZURMTkihkbE0H79+rVLgGZm1rJqNZi/nv7bt5nPq5A0Dvg4Ii6vcGzWRu41btb1VOvO46/Am8DAgrIhwKOlKkvaHxgOnJQub13pAM3MLLuq3HlExHJJVwNHAg+kyWAH4DBJg4GfAqMiYpmk3YETgQOAtSQBnA18uxqxWn7uu2HW9VSzzaMRkKQ5wI3AIemruOuQdBxcLa13I9BAMgHVu+nP56sYp5mZtaBqnQQjYinJbIPF5XOATQqWNymuY2ZmtcVjW5mZWW5OHmZmlpuTh5mZ5eaZBK0k990ws3J852FmZrk5eZiZWW5+bGUlueOfmZXjOw8zM8vNycPMzHJz8jAzs9ycPMzMLDcnDzMzy83Jw8zMcnPyMDOz3NzPowt56Oend3QIZtZJ+M7DzMxyc/IwM7PcnDzMzCw3Jw8zM8vNycPMzHJz8jAzs9ycPMzMLDf386hj7rdhZh3Fdx5mZpabk4eZmeXm5GFmZrk5eZiZWW5OHmZmlpvftqoxl8x4IXPdYRWMw8ysHCePGjPs75M7OgQzsxb5sZWZmeXm5GFmZrk5eZiZWW5Va/OQ1Au4ChiUHvesiLi3mbqnAYemi7+KiEnVibIy3AhuZp1NNRvMGwFFxDBJA4A5kgZHxMLCSpL2Bo4BdkiLnpT0bETcWcVY25Ubwc2ss6nKYytJ3YCjgZ8DRMQLwBPA4SWqHwvcEBFLI2IpMA34TjXiNDOzbKp157EFsB7wfEHZs8DQEnV3AW4sqnd85UJrHY9oa2ZdWbWSx4bpv28XlL0FDGmmbnG9DUrtVNI4YFy6+KGkp9sYZy1bH1jc0UFUUGc+v858btBVz+/oi6ofSWUMbM1G1e4kGEXLylivdKWIycBkAElzI6LUnUyn4POrX5353MDnV+8kzW3NdtV6Vff19N++BWV9C8qL6xbXW1ShuMzMrBWqlTz+CrzJyrdHQ4BHS9R9NGM9MzPrIFVJHhGxHLgaOBJA0tYkr+JOkzRY0kxJ3dPqVwGHSOqV9g05NC1rSWd/H9bnV78687mBz6/eter8FJGpeaHNmuskKGkY8D/AlumruZ2uk6CZWWdTteRhZmadh8e2MjOz3Oo+eaRtI1MkzZE0V9Kojo6pvaTnNbvop3dHx9VWknaR9BdJY4vKd5D0kKQHJd0hab0OCrHVypzbgqLr+F8dFGKrSVpN0slp/Pen12rPgvV1e/0ynFtnuH4npO3LsyQ9JemEgnX5r11E1PUPcAFwXfp5AMlbXRt2dFztdG5TOjqGCpzT/sANwGPA2ILy1YGXgOHp8nnAzR0db3ucW2e5lkB/YB6wTro8EvgnsEm9X79y59aJrt8fgQ3Sz1sDnwC7tvba1fWdR84xs6w2PBoRhwLvFpXvAyyLiPvS5WuAb0jqV9Xo2qa5c+ss3gXOjYi3ASJiBrAU+Bfq//qVO7fO4lsR8TpARLwILCFJmq26dnWdPMg3ZlZdknSVpAck3SVpeEfH01YR8XIzq3ah4DpGxD+A94GdqhFXeyhzbgCD02v4B0mTJa1ftcDaSUS8ERHXNy1LEslfrYuo8+vXwrlB57h+f2r6LOmbJHdW99LKa1fvyaO5MbNKjoVVh54Dfh4RXwbOBe6QtEML29Sr4jHNoHNdy2eAg4Evk5znXekvqHq2B8njjt/T+a5f4blBJ7l+kraV9CxwMTA6It6ildeu3pNHk6xjZtWViLgwIh5NPz8K3EkyZH1nVeq98c5yLY+MiHcieajcSPJX3a4dG1Xrpf22fkjStrM8Le4U16/UuXWW6xcRT0XEEODfgd9IajqH3Neu3pNHnjGzOoO/A5/v6CAqpHhMM+ik1zIi3iN5saMur2X6F/dk4JKIeCwt7hTXr5lzW0m9Xz+AiHiC5I/RE2nltav35JFnzKy6I+m7RUUbAq92RCxVsNKYZpI2BdYkeXOprkkaIWlowfLqwGeo32t5EfBIRNwsqaekzeg812+Vc+sM10/Semk7R6H3gLVo5bWr6+QRZcbM6si42tGpkjYAkLQ58G/A9eU3qVt3AT0k7ZEuHwncGhGdYUTlzYDjCp6Rn0Dyh8/DHRdS60g6k2R4oSlpn6MtSa5V3V+/MufWGa7f2sA5ktYESN+k+ndgJq28dnU/PImaGTOrY6NqH5JOJ+k78AnJXwiXFb4RUo8k7Uzy190OwALg2Yj4RrpuR+AKYBnJa4RjI+KNjoo1r+bOLf3LvJGkH5JIXgs9MX21vG5IGgD8ucSq8yKisZ6vX7lzA35BnV+/9Pfk94C9gI9IHkvdBpwfEctbc+3qPnmYmVn11fVjKzMz6xhOHmZmlpuTh5mZ5ebkYWZmuTl5mJlZbk4eZmaWm5OH5SZpuqSGjo6jrSRtomQSsSgoO1TSNR0ZVzVJulXS0o64npLGSZovaUq1j21t5+RhrTEauD9LRUkNkuZXNpzWiYhXSM6l0K+Bk7NsL6mxEr/4lMwg2Zij/mwVzVyYVdpBc0Frts2r+PuKiMnAlGY3sJrWo6MDsPoTEe90dAyVEhHLSOY5MLMyfOdhuUg6I53PuTFdvkLSW5K+L+lmSS9I+mG6biBwKbBR+tfxzWl5b0m/SCfW+T9J30nLN2t6jCRprKR7JX0o6URJzyuZW/pHkh6TdE9a//p03Q9KxDlHyURaP0kHs2tad5akZyTdDXy9oHxbSU8W3imlj1ZmKZn7eaakIWn5wcBYYO/03L6Xlm+Zxn1/euySM9FJ6ibpyvQ7+L2kayStJekkYG9gbLrfowq+55lp2Y2S+qTlE0mGQ5mQrvtqWv4VJXNS369kTurP5rjGY9Lv7n5JNxQcq9lrna7vLulSJfNj3y3pvPSR2K3NfV+pXpJ+JulxJY9Ee2WN1TpQR8+r65/6+yF51NDXsnEmAAAEwUlEQVRYsDybZHhnARsDHwOfTdc1APOLtr+aT+edXxv4G/CldLk/ydwCY9LlU9N9jiW5I9gsPU7TkNKrk0xa8yGwcbrNYSQTaa2Z1r0JODtdty/wGrBuunxh8r/BithWipdk/pSeBeseKFjXSMHc1iR38s8BR6bL2wGLgbVLfIf7AncVLN8G9C/1/aZlJxYd9/tF3//YguXNScZfGpgujwd+V+Z6zgca0s//msbcL12eBFyT8VofTzJPdq90/c1F3+VK31dB2Ssko9R2A54CDuno/8b90/KP7zysvdwTideAN0iSwCqUzDv/LZLB5oiId4E70rJCv0nXX5zuE+DPEfH3SH7rPAM8HxEfRTIv8yKSX5qQJJpfRcT7ad0bC/Z/IDA9It5Ml3/dwnk9SzKD4wPABcDOZeruRjIS6/Vp7H8i+cW4X4m6S4BtJY1Mv5NDSOZrac7S9E7mfpJ2mnJxHArMjYimgf5uAPaUtHGZbZqMBe6IT0dUvQE4TFpp1rzmrvWBwC0RsbTge8/i4YhYEsko2U/z6XW0GuY2D2svhe0gS0nuCErpB/QEfiTpg7SsL/BkYaWIKJ4WE5K/ppt8UmK56ZifAw7Vp3O+9wKaZrvbmOSv4yZv0gxJ6wC/BY6KiFsk9QfmNVc/PW4AMwp+1/YE1imuGBEPSRoHnEmSSH8GTGwmjgaS0Xq3jYj5aeP42BbiGCJpdkHZSyTzwbxWcovmt+0BLATWI7kjgeav9cYFdaDMd1ukcH8f0vx/O1ZDnDys2haR/IL4j0in2JW0GskjpvbyD2BGRExqKpC0fvrxNZIE1mS9MvsZCPQB7k6XV8tw3I8joqHguGvxaeKioHwdYHZETJe0ZXqMV4BrS+x3V5K7rvk54pgbEV8tON5nWPmXdLlt/xYR4wu2XT8iFpfZpkme79bqnB9bWaW9S5oYJP0E2AT4JSs/pjobGNOOx5wCHNjU8JregfwsXXcTsK+kpl9sh5bZz0skdzS7pct7F61/F1hTidtIJgf6u6Sm+Ul6AP9LMg9Esf2BcQAR8VfgZaB70X7XkjQN+AuwVUHMX2kmjq0lTSJ5XLSbpM+ncWxA0laR5f/3KcBX02TT9NLDHRm2g+S7PUBSr/Qx14HNxNn0fVk96+hGF//U1w9wBkm/gPnAUcCPgLeA54Evkkwos5TkMdQQkrvbWcBDJA2oAnoD1wD/R9Jf5FKSX5zrAnNIHv3MBoakxxyR7v8tksl5xhXEsH/xMdNtTiP5ZT6L5Bf4BgXncBZJW8YM4D8Ljrd9uo+lwM1p3e+kx/ktcEla99503VYkbS8PAWekZU13EfcDD5A2npf4HgeSNDzPAh5Jv4/V03VfTM/3EZLk1i1d/xfgVpLk+xbwo7T+/iQTGT0CDE/LRgEPpuc1CxjWTBy3Fnx3O6dlh6fXYVZ6LgPS8izX+jKSdovpwARgXsGxVvq+0nObn17L44qu66Ed/d+6f8r/eDIoM2sX6Z1et4h4P10+EDg9InYrv6XVIz+2MrP28iXgHFjxVt2hwD0dGpFVjO88zKxdpG+jTSZ5u60nyeOskyPigzKbWZ1y8jAzs9z82MrMzHJz8jAzs9ycPMzMLDcnDzMzy83Jw8zMcvv/5SLiLPS+hpAAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.bar(*p_end(PYTHIA_val['endings'], PYTHIA_val['ending_counts']), label = \"Pythia\", **PLOT_END_SETTINGS)\n",
    "plt.bar(*p_end(JUNIPR_val['endings'], JUNIPR_val['ending_counts']), label = \"JUNIPR\", **PLOT_END_SETTINGS)\n",
    "plt.xticks([0, 5, 10, 15, 20, 25, 30])\n",
    "plt.axis([0, 30, 0, 1])\n",
    "plt.xlabel('intermediate state length')\n",
    "plt.ylabel('probability to end')\n",
    "plt.legend()\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# P_mother"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYsAAAERCAYAAACKHYuuAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzt3Xm8FNWZ//HPF0QQURFFTJSIS1jdAcElBFAYNeIkJmo0LrgMJpq4KyRuaPKLSRyiTuISnAmooEYd84uJqCiLGiMiihrRRBPRuKGiqLigLM/8UXWxafp214XeuPf7fr36dW+dPl3n6b5QT9c5VecoIjAzMyumVa0DMDOz+udkYWZmJTlZmJlZSU4WZmZWkpOFmZmV5GRhZmYlVTVZSOov6R+SRpao9x1Jj0uaI2mcJFUpRDMzK6BqyULSN4AzgPdL1NsRGAf8G7AHsDtwcsUDNDOzRlXzzOKxiDgSWFyi3gnAlIhYGBErgN8C3614dGZm1qiqJYuIeDVj1f7A33K2nwX6SNqg/FGZmVkW69U6gAK6sGpX1XuAgM2BV3IrShoFjALYcMMN+/bs2bNaMZqZNQuPP/74wojoXKpePSYLgEITVq02yB0R44HxAP369Ys5c+ZUOi4zs2ZF0stZ6tXjpbNvAR1ztjuSJI+3axOOmZnVY7J4DOiRs90bmBcRn9QoHjOzFq/myULS5pIektQpLfpv4EBJm0lqBYwErq1ZgGZmVtX7LPpKmgnsCoyRdEf61AZAT6A9QEQ8A5wNTAUeBZ4Erq5WnGZmtrqqDXBHxOPA4ALlrwCd88omA5OrE5mZmZVS824oMzOrf04WZmZWkpOFmZmVVK835ZmZld3l9z1fss4Zw7qvVnbfffdxzjnn8NRTTzFo0CCWL1/OokWLOOmkkzj11FMzt//ZZ58xfPhwHnjgAebPn0+3bt1Wq3PppZeybNkyLrjggsz7rQYnCzOzEoYNG8YVV1zBkCFDmDZtGuuttx7z5s1jt912o1evXgwbNqzR10pamRjWX399Zs6cSbFVF84880wiCk1iUVvuhjIzWwN9+vRhp5124p577inrftu2bUu7du3Kus9ycLIwM1tDH3zwAXfeeScbbLABJ554IgCXXXYZm2++OWeddRYHHHAAAN/+9rcZPHgwr7322srX3nXXXRx44IFsv/32/P73vweS7q6ePXsyePDglfUuueQShg4dytChQznooIN4/fXXq/cGczhZmJmtgZkzZzJ//nwmTZrEihUrGDt2LACnnXYagwYNYty4cdx9990A3HLLLcycOZOtttpq5es/+ugjpkyZwqWXXsqYMWOApLur4fcGm266KdOmTWP69Ol861vfYvTo0dV5g3k8ZmFm1gT77rsvy5cvp3Xr1tx2220MGDCAYcOGMWnSJMaMGcNdd9218oyimP333x+AXXbZhfnz5zdar2vXrgwZMoQVK1bwwQcf8Nlnn5XtvTSFk4WZWRM0DHDnOuaYY7j44osZM2YMt956K9dcc03J/Wy88cZAMkaxdOnSgnVeeOEFDjvsMB5++GH69+/PzJkzGTly5Fq/hzXhbigzs7V08MEH8/rrr3PvvffSqlUrOnbsWPpFGcydO5eNN96Y/v37AzSaVKrBycLMbC21a9eOQw89lOOOO47DDz98lec6dOjAxx9/zKRJk7j99tubtN8ddtiBRYsW8fzzyf0h5b7yqkkiolk8+vbtG2ZmlTB16tTYZZddAohBgwbFzJkzV6vz0EMPxWabbRafffbZKuWjR4+OnXbaKfbee+944403YtiwYQHEgAED4r333osBAwYEEMOGDYupU6dGjx49YpNNNonvf//7ERFx/vnnxzbbbBMHH3xwjBo1Ktq2bRtHH3102d4bMCcyHGMVdXjzx5rwsqpmVkvPPfccv/71r7nqqqtqHUqTSHo8IvqVquduKDOztXDzzTezbNkyJk6cyLHHHlvrcCrGycLMbC088cQT7LbbbixcuJA99tij1uFUjLuhzMxaMHdDmZlZ2ThZmJlZSU4WZmZWkpOFmZmV5GRhZmYlOVmYWcsx49LSjzzTp09n8ODBSGKXXXbhZz/7Gd26dWOTTTbhmGOOWVnvl7/8Jd26dWPXXXflvvvu45prrqFnz5507NhxlTUoZs+ezcCBA9lyyy056aSTgGQm23bt2tGjRw8mTpzIX//6VwYOHIgkBgwYwFe+8hX69OnDeeedx4oVK1Z7fvDgwfTo0YMRI0bw5ptvVuazy3Kb97rw8HQfZlbS9J+WfjQCiPvuuy8iIi666KLYe++9V6tz0UUXxYwZM1ZuT5gwIVq3bh0jRoxYpd78+fPj2GOPXaVsm222ieuuu26VOkC88MILERGxYMGC2HTTTWP8+PEFn//kk09ip512Wm2/pZBxug+fWZiZVdCpp57Kvffey+TJk9dqP126dGHIkCGNTibYrl07hg8fTqXuN3OyMDOroJ133pkLL7yQ0047jbfeemut9rV06VLatGnT6PPLli1j6623Xqs2GuNkYWZWYaNHj2bbbbfllFNOWeN9PPfcc0ybNo1vfvObBZ9/5513ePXVVxk3btwat1GMV8ozMysjSauVrbfeelx//fXsvvvu3HHHHey+++6Z9/ed73yHtm3bsnz5ci6//HIOPfTQ1Z5v3bo1TzzxBGeccQZ9+vRZ6/dQiJOFmVkTtW3btuBa2J9++ilt27Yt+JrevXtz8cUXc/LJJzNlypTMbU2ePJkddtih5PNXX301p5xyCkcccQQ777xz5v1n5W4oM7Mm6tq1Ky+99NJq5fPnz6dr166Nvu7ss8+mW7dunHHGGWWP6Xvf+x477rgjV155Zdn3DU4WZmZNNnz4cD766CMmTJiwsmzevHm8+eabbLXVVo2+rnXr1lx//fXMnj277DFJYsyYMdx0000sXLiw7Puv+f0R5Xr4PgszqyQg7r///pXb9957b/Tt2zd23HHH2HvvvWPEiBHx/PPPr3z+6quvjh49ekSPHj3iN7/5zSr7Gjdu3Cr3QwwdOjTatm0b3bt3jwkTJsTTTz+9crnVPfbYI2655ZZVXp///LRp0yIiYtmyZbHddttFnz594rbbbsv6vrysqpnZ2li8eDEbbbQRH330ER06dODxxx9v0uD0uqDu1rOQ1E7SREmzJM2RNLyRehtJulHSY5Jmp6/pUK04zcwaHHDAAQDMnTuXDh060L179xpHVDvVvBpqLMnKfAMldQdmSeoVEfkTmVwAbAMMBAKYBpwPjKlirGZmfPGLX2TgwIEsWbKE8ePH06FDy/3eWpVkIakVcCJwCEBEPC9pLnAUkH8HSR9gdkQsT187G9itGnGameW69dZbax1C3ahWN9R2wGbA33LKngUK9ZNNAYZK2lBSe2Ao8GjlQzQzs8ZUK1l0SX++n1P2HrBFfsWIuAqYCbwIvAQ8AVxSaKeSRqXjH3PefvvtcsZrZmY5qn2fRf6lV6vdFy/pPGAXknGLL5GclYwsuLOI8RHRLyL6de7cucyhmplZg7VKFpJOz1i1YarFjjllHXPKc50K/CYilkTEEuBa4MdrHqWZma2tTAPckjYHTgC2BdbPeWp/4IoMu/gn8C7Qg88TRG+S8Yl86wNLc7aXAhtlidPMzCoj65nF/wcGkRzwX855fJrlxRGxArgOOB5A0peBXYHJknpJmiapdVr9fuBwpYAjgRkZ4zQzswrIeuls64j4Wn6hpCea0NZY4FpJs9J2j4iIBZK6AT2BNsBy4BTgV8AskjGNfwGjmtCOmZmVWdZk8YCkL0XEv/LKdwf+mGUH6fjDyALls4CtcrbfAg7PGJeZmVVBo8lC0vSczVbAyZKeBz7IKd8VuLhCsZmZWZ0odmbxBeBnOdsT8p5XWsfMzJq5Ysni3Igo2sUk6eUyx2NmZnWo0auhchOFpG6SjpG0Tbq9jaQREeGrlMzMWoCsl87+nGQep4ZLZRcDIyT5ZjkzsxYg69VQXSJi5RVKEfEuMErSg5UJy8zM6knWM4s2+QXptOPVXA/DzMxqJOvB/kFJ9wI3Am8DnYGjSWaHNTOzZi5rsjgfGA1cCGwNvAJMBC6rTFhmZlZPMiWLdNW6n6aPlSR1JFmXwszMmrHMYw7pGMUWJLPCBslNebcAe1UmNDMzqxdZpyj/HvALYMO8p/IXMzIzs2Yo69VQZ5Kshb1BRLRqeJBMXW5mZs1c1m6oZyLisQLlZ5YzGDMzq09Zk8Ujkn4N3MOqs85eQTJNuZmZNWNZk8XFwAIgfwGkLuUNx8zM6lHWZDElIr6ZXyjp+jLHY2ZmdSjTAHehRJHKX+PCzMyaoabcZ7E9sBvQPqd4DNC73EGZmVl9yXqfxX8AZwBbAk+RTCy4M/Bi5UIzM7N6kfU+i5EkyeGpiBgSEfsAPYA5lQrMzMzqR9Zk8WFELCOZ6gOAiHgD6FWRqMzMrK5kHbNYKqk/8IykO4D7gQEUWOfCzMyan6xnFuemP8eQ3JR3CtCRZE0LMzNr5rKeWXwLeDWd8mNk5cIxM7N6lPXM4gTgoUoGYmZm9StrspgD/CO/UNLp5Q3HzMzqUdZuqHnA3ZLuAt7PKR9FMpmgmZk1Y1mTxQ+AJ4FD8sq3LG84ZmZWj7Imi5si4pT8Qkk/L3M8ZmZWh7JOJHgKgBKdJSktH13J4MzMrD5kShaSOkmaDHxCsq7FJ5ImSepU0ejMzKwuZL0a6nrgQ2A/YEdgOEniuCFrQ5LaSZooaZakOZKGF6m7l6Rpkh6U9KykH2Rtx8zMyi/rmMUWETEiZ/s54EFJs5vQ1lhAETFQUndglqReEfFmbiVJ25JcYXVgRCyU1Ac4vgntmJlZmWU9s/iXpNx1LJC0IfBGlhdLagWcCPwPQEQ8D8wFjipQ/UxgQkQsTOvOi4izMsZpZmYV0OiZhaQLczbfAuZIuhdYBHQCDgQeyNjOdsBmwN9yyp4F+hWouy/wD0lTgE2AR4DzI2JJxrbMzKzMip1ZfA/YNn20Bx4lmTxwW5KD+F+Ar2Vsp0v6M/eGvveALQrU7QZ8FzgGGEyyEl/BG/8kjUrHP+a8/fbbGUMxM7OmKjZmcXVE/LjYiyVd0MT2In8XBeq0JbmvY2HaxhXAnZJOjogVq+wsYjwwHqBfv375+zYzszJp9MyiVKLIWif1VvqzY05Zx5zyXIuA3EHvV0kSyOYZ2zIzszLLOsC9tv4JvEuyFGuD3sBjBeo+yardU52Bz4B3KhadmZkVVZVkkXYfXUd6CaykLwO7ApMl9UrvqWidVr8OOCLn6qvjgUkRsbwasZqZ2eqy3mdRDmOBayXNSts9IiIWSOoG9CRZonV5RPxO0nYk92EsBv4OeCp0M7MaUkS2cWFJGwMHA1uTjCPcGREfVDC2JunXr1/MmTOn1mGYma1TJD0eEYVuY1hFpjMLSf2Au0kGn98mGUf4paT9I+KJtYrUzMzqXtZuqMuAb0TEnxsKJO0NjAOGVCIwMzOrH1kHuJWbKAAi4mEK3ydhZmbNTOaroSR9JW97H1a/yc7MzJqhrN1Q5wBTJL3P52MWGwH7VyowMzOrH5mSRUQ8Jml74CCSq6FeAf4UEYsrGZyZmdWHrFdDHRcRE4CbcsomSfpTRNxSsejMzKwuZB2zOLpA2Sh8s5yZWYtQ9MwiZ02LbnnrW0AybfkmFYnKzMzqSqluqG3Tnxvk/N5gMXBc2SMyM7O6UzRZRMRxAJJOjYj/qk5IZmZWbzKNWThRmJm1bNVaz8LMzNZhThZmZlZSpmQhqU+lAzEzs/qV9cziTklfrWgkZmZWt7Imi/eAvpL+KOmnknpVMigzM6svWScS/LeIWEiy4NHOwPGSdgTuAW6OiLcqFqGZmdVc1mTRHVgoaf30954kix6tB+whqQ3w84h4vDJhmplZLWVNFtdJegg4DHgNmAR8NyJeA5DUiWTZ1QEVidLMzGoqa7LoCiwBhjVy9jAQ2LJsUZmZWV3JOsD9s4g4PTdRSDpE0hbp5iP4rMLMrNnKmiz2K1C2AXA7QEQsiogFZYvKzMzqSqkpyr+U/tpOUldAOU8/gacoNzNrEUqNWbwERPr7y3nPLQauKndAZmZWf0pNUd4KQNKUiDiwOiGZmVm9yTpFecFEIcmLH5mZtQCNnllIGhQRD6a/5y+p2mAkMKECcZmZWR0p1g01WtITEfEh8D2SqT3ytatMWGZmVk8aTRYR8bWczXER8Z/5dSSdXZGo1sTiBTDj0tq0PeSHtWnXzKxKso5ZrJYoUm3KGIuZmdWpYmMWv83w+v2BGn2dNzOzaik2ZjEEmFji9Z+WL5S189Gny3jkxXdq0vaeQ2rSrJlZ1RRLFj+NiOuKvVjS61kbktQOuJZkevP1gB9FxNQi9dsAzwKTI2Js1nbMzKz8Gh2zKJUoUoXmjGrMWEARMRA4ErhFUpci9UcBWxR53szMqqTYmMV44PSI+FjSi4WqAMUO9rn7agWcCBwCEBHPS5oLHAWMK1C/A8naGXdm2b+ZmVVWsW6oGcAn6e/vA6fnPS/g8oztbAdsBvwtp+xZoF8j9c8CfgUclHH/ZmZWQcXus7g5Z/PoiHgmv46kozO203AG8n5O2XtA7wL77AwMioiLJRVNFpJGkXRX0aXTRhlDMTOzpsq0Ul5EPCNpP+DbwBeAN4CbI2JaE9uLvG0VqHMB8P8yxjUeGA/Qq9uW+fs2M7MyyXRTnqTRwA3p5jMkB/lJks7N2M5b6c+OOWUdc8ob2tkO2DYipmfcr5mZVUHWNbiPAvpExKKGAkmbAQ8Av8jw+n8C7wI9+DxB9Aam5NX7KvBFSTPT7Z7AEkmDgeMiYn7GeM3MrIyyLqv6am6iAIiId0i6o0qKiBXAdcDxAJK+DOwKTJbUS9I0Sa0jYkJE9I2IwRExmGTywonpthOFmVmNFLt09ks5m1Mk/QS4DVgEdCIZv2jKmMVY4FpJs9J2j4iIBZK6kZxBtAGWp22vD0zl8zOLnSPikCa0ZWZmZVSsG+olkgHp3EHoH+XVCeBnWRqKiCUk61/kl88Ctsor+wwYnGW/ZmZWecW6oR6IiNYR0aqxB/BgtQI1M7PaKZYshmd4faZLXM3MbN1W7Ka8pbnbkrYHdgPa5xSPocCNdWZm1rxkunRW0n8AZwBbAk+RDEbvDBSaM8rMzJqZrJfOjiRJDk9FxJCI2Ifknok5lQrMzMzqR9Zk8WFELAPWbyiIiDeAXhWJyszM6krWO7iXSuoPPCPpDuB+YABeg9vMrEXIemZxLsn9FmOAD4BTSOZ2yjrrrJmZrcOyzjr7bM7myMqEYmZm9SrrrLOSdKKk+yXNS3+eIKnQFONmZtbMZB2zGAcMBW4H3gE2B04jucfirMqEZmZm9SJrstgX2CMiPm0okPSfwKMViWpdM+PS2rQ75Ie1adfMWpysA9yv5SYKgIj4BFhQ/pDMzKzeFJuifFDO5oOSbiTphmqYovxw4JHKhmdmZvWgWDfUvSSLG+UOYu+TV+ffgYvKHZSZmdWXYsliSkR8s9iLJf1vmeMxM7M6VGzW2dUSRXqpbBfgzUgUTSYtxSMvvlOTdvccUpNmzawFynqfxYaS/hv4BHgN+ETSeEkbVjQ6MzOrC1mvhrqSZAnVoSTrYu+bbl9ZobjMzKyOZL3PontE5F4d9QLwsCQvq2pm1gJkPbNoI6ltbkG6nTXZmJnZOizrwf4OYI6k3wFvA52Bw4CJFYrLzMzqSNZZZy+T9AZwDLA18ApwaUTcXMngzMysPmRdg/tC4NWIGF7heMzMrA5lHbM4AXiokoGYmVn9ypos5gD/yC+UdHp5wzEzs3qUdYB7HnC3pLuA93PKRwFXlD0qMzOrK1mTxQ+AJ4FD8sq3LG84ZmZWj7Imi5si4pT8Qkk/L3M8ZmZWhzKNWTQkinQt7s4Na29HxOhKBmdmZvUh60SCnSRNJplIcAHJRIKTJHWqaHRmZlYXsl4NdT3wIbAfsCMwnCRx3FChuMzMrI5kHbPoEhEjcrafI1lqdXbWhiS1A64lmbV2PeBHETG1QL0RJAPq6wEdgYkR8V9Z2zEzs/LLmixek9Q+Ij5uKEjXsng5Z/u4iJhQZB9jAUXEQEndgVmSekXEm3n1LgcOjYi5krYE5kl6OSL+kDFWMzMrs6zJYiHwZHqfxSKgE0lX1LR0KhCAkUDBZCGpFXAi6aW3EfG8pLnAUcC4vOpXRcTctN4CSTPStpwszMxqJGuyOAi4h6RbqGNa9ijQIX0AtCvy+u2AzYC/5ZQ9C/TLrxgRl+cVtSOZ6dbMzGoka7K4OiJ+XKyCpAuKPN0l/Zl79/d7QO8S+9wY6A+sdo9H+vwokrvI6dJpo2K7MjOztZD1PouiiSJrHZKlWHOpRP2fA5dExMuFnoyI8RHRLyL6bbpR+wzNm5nZmsh66ezaeiv92TGnrGNO+WrSs4alEXFVJQMzM7PSqpUs/gm8C/TIKesNPFaosqRvAEOA09LtL1c6QDMza1xVkkVErACuA46HlQf/XYHJknpJmiapdfrcIOBU4PvAhpI6AOdXI04zMyusWmcWkN5nIWkWcDNwREQsADYhuVGvTVrvZmAwyeW6i9PHNlWM08zM8mS9GmqtRcQSknsx8stnAVvlbG+VX8fMzGqrasnCKmDGpbVpd8gPa9OumdVMNbuhzMxsHeVkYWZmJTlZmJlZSU4WZmZWkpOFmZmV5GRhZmYl+dLZddgjL75Tk3b3HFKTZs2shnxmYWZmJTlZmJlZSU4WZmZWkpOFmZmV5GRhZmYlOVmYmVlJThZmZlaSk4WZmZXkZGFmZiX5Dm5rulotugReeMmsRnxmYWZmJTlZmJlZSU4WZmZWkscsrMlqNdsteMZbs1rxmYWZmZXkZGFmZiU5WZiZWUkes7B1S63u8fD9HdbC+czCzMxKcrIwM7OS3A1l65RaXbbrS3atpXOyMMvC82FZC+duKDMzK8lnFmYZ+K51a+mcLMzqnS8XtjpQtWQhqR1wLdAzbfdHETG1kbpnAUemm7dExGXVidKs/nhQ3+pBNc8sxgKKiIGSugOzJPWKiDdzK0naH/gPYNe06ElJz0bEXVWM1cw8qG85qpIsJLUCTgQOAYiI5yXNBY4CxuVVPwm4KSKWpK+dDHwXcLIwq6JajtPw4tk1aXbP7TarSbvrQnKs1pnFdsBmwN9yyp4F+hWo2x+4Oa/eyZULzcwsUbMEWaPk2BTVShZd0p/v55S9B/RupG5+vS0K7VTSKGBUuvnhXieO+/saxrc5sHANX7u2atW233Pzb7eWbfs9rzttb5OlUrWvhoq8bWWsV7hSxHhg/FpFBEiaExGFznIqrlZt+z03/3Zr2bbfc/Nru1o35b2V/uyYU9Yxpzy/bn69tysUl5mZZVCtZPFP4F2gR05Zb+CxAnUfy1jPzMyqpCrJIiJWANcBxwNI+jLJpbGTJfWSNE1S67T6tcARktql92YcmZZV0lp3Za2Dbfs9N/92a9m233Mza1sRmYYH1r6hRm7KkzQQ+F9g+5zLZX1TnplZHalasjAzs3WXZ501M7OSWnSySMdFJkqaJWmOpOFVbLu/pH9IGlnFNttIOl3STEkPSHpE0r5VbP8H6fjUdEl/lfSDarWdtv9lSUslDa5SexPTzzr30aEabaftnyjpYUl/lvS0pK9Woc0Fee93nqQ7K91u2nYvSTPS9ztX0ugqtdtF0u2SHkr/XR9X4fYKHjsk7Zr+n35Y0h8llfV29JY+6+xYMsxXVW6SvgEcyqo3H1bDVsBpwK4R8b6kYcAfJPWIiNeq0P6JwLCIeCu9yOE5SY9GxOwqtA1wCfBZldoCICIGV7O9BpIOBfYFBkXE8vQAtmUVmr4jIlbOuCDpSuCJKrQLcD1wX0Sclx4oX5D0VETcU+F2bwYeiIhvSfoS8Iyk5yJiVrkbauzYIWl94A/AyIiYIelikjHiQ8vVdos9s8iZr+p/IJmvCmiYr6rSHouII4HFVWgr12Lgwoh4HyAi7gOWAHtVqf2jI+KttO0XgEVAt2o0LKk/8CEt556dC4FLImI5QERMiIjfVbrRvETRGhgB3FHpdlN9gEfSON4BXgB2q2SDkrYGhgC/Tdv9F3AfcEKFmmzs2HEAsDwiZqTb/w0cIqlzuRpuscmCps1XVVYR8Wql22ik3Xci4saGbUkC1qdKB9CIeDqn7W+SHLwLTlNfARenj6qSdG3aPXG3pKpM+i1pC6AXsGvaLfOQpJOq0Xae/YBZEVGtL0V3kSQnJG1HkjwerXCbW6U/c6faeJPPZ80uqyLHjv7kHMsi4hXgY2D3crXdkruhmjJfVXP1VeBl4MFqNShpJ+B3wIbAYRHxXhXaPACYFxGvJvmxap4DpkfEY+mZzQxJ+0TEkxVutxvJVDrfIDlgbwHMlvR+RNxS4bZzfQe4sWSt8jkBuFPSP4FOwJkRMb3CbTYcvLcC/pH+viWwaYXbzZc/px4UmVdvTbTkM4sGWeeralbS+15+StLHuaJa7UbEXyOiN/B1kvGSPSrZXtrdOBqo+uIMEfHziHgs/f0xkm++1fiG35bk//avImJ5RLxBctA+vgptAyCpPUn3ZrXOHCHp7nokIrYHdgbOlFS2b9aFpGN99wOnK9GLZKxoSSXbbSycAmVlO5615GTRlPmqmpW0+2k8cHlEPF6LGCJiLsnB89QKN3UkcG9EvFvhdrL4Fxln+FxLi9KfuRdqvApsXYW2G/w78KeGMZNKSw/S+wFXwMpumPuAc6rQ/LeB9sDDwEUk4wUvV6HdXPlz6kGZj2ctOVk0Zb6q5mYcMDsibpPUNr2Co6IkbZaOU+T6iKQ7qpK+AhzUcCknSRfBFZL+UOF2kXRuXlEX4PVKt0sysPsxq3ZBdK5S2w2q3QW1fvpzaU7ZUmDjSjecjgUeHxF7RcS3SQ7Sle7+yrfKnHqSupIksLJ9GWyxyaLYfFW1jKvS0mvP1wMmptf8b091uic2Ai5IuydIr9L4OjCtko1GxEkRsXdEDE4vY10AnB4R/17JdlNnpoPNSNqW5Nt2xQ+gEfEpcAPpFTmSNgQOJ7m0tOIkbQ50rfJZ69+A14Aj0hg2Ag6mCgdtSTem7xlJvYF9SM4uquluYL2ce2mOJ7mMuWwXr7To6T7UyHyoTPcqAAAHL0lEQVRVVWi3L8m3+11JDl7PRsQhVWi3O1BogaiLI2JshdtuB5xH0lXwGcm3r9+TXN5Z8TGTdGzkF8BAkgPLHRFxSYXbPJtkkHkZyRnUlblXo1W47Q1J/m33Ttv/X+CyqMJ/eEknAx0j4qeVbiuv3X7AL0n66TuQjCX8MCKWVbjdS0kS0yKSK/zOjIhnK9RWo8cOSbsBVwPL01hGppcQl6ftlpwszMwsmxbbDWVmZtk5WZiZWUlOFmZmVpKThZmZleRkYWZmJTlZmJlZSU4WtpKkKWuyMJCkvpKelPRSvcRka6+Sf1db9/g+C1tJ0sbA4jW5cSs9oE+MiG71EpOBpG7A/IhYownl1vTvKmks0C0iRmasPxF4qdI3h9qaa8lTlFueiPig1jHkq8eYzFoid0MZAJLOUbJ+8th0+2pJ70n6saTbJD0v6ac59VtLuiJdc3gKyTQaufvrIOm3StZD/ouk76blR0t6Q9KbkvaRdKGkRZLGrW1MBV7fRtJlaft/lnRROo30+unEgiHplLSr659KlqwsFf+XlKzZHpJGSpoq6VNJ3SRtr2ShoUckTZZ0h6SXJJ0k6YX0ff4k3c84Se9IWm1WVEl7pe1OT+M8KOe5Y9L2H0z33yUtvzD9rH4l6SZJf5f0e+CW9PmGNbFbN/a5ZPm7Foj1ovT9zpD0O0lfkHQ4MBLYP23zvJwYp6ePP0n6Ylp+GrA/MDKt3zCnVd/0fT6gZO32nsVisQqLCD/8ICIAJgJjc7ZnkkwjLuALJLN4fjF97mTgKaBd+vwtJN0IDa+9Drg+/X0j4EVgn3R7T5K5a7oChwDfLUdMBV57HslEcq2BNsBfgKNyng/g3PT3w4C/Z4y/W/raY9LtM9NYHiWZiwiSxXDeb4idZMWy94B26XYX4H8aiXs2MCD9fReSbiBIZtB9C+icbp8PTMv7rP5KMtvoJunz3ZL/5tk+l1J/17z99CZZXbKhO/tyYHD6+9iGuHPq/yCn7kjgxiJ/501IVnAcmm5/jWRes1a1/n/SUh8+s7BS7o3EG8A7fL5m9qHA7RGxJJL/zbc2vEDJgkNH8/m6xIuBP6ZlRMQjJAeHG0gO0r8pU0z5RpIc8JdHxFLgtoYYctyT/nwa2DZL/Dn+kD7/S5IpsvcAJqVlrwEPNFSMiCdI1pT4elp0JOm3/gLeBY6W1CUiniI5gAMcQ7JGRMNMohOAoVp1ivn7I+LjiHg/In7SyP5H0vjn0ujftYDFJFO+HyKpDckiU38uUv8VktUCHwROB/oWqXsQ8GGkK91FxF1pWwOKvMYqyMnCSskdM1jC5+sGfIFV1x3OXVyoM8lqbb/Q5+tIfDUta3A+ycHi4fSgVI6Y8m1NMk14QwzfIfk2XWhfS0i+ZWeNn4jIXcbyC+nPxj4TSJLjMenv+9L49OxHkqxH8YSke4DuOe8nd8rpt3PKG+QvrVlIsc+l2N91FZEsMPQ1kkTzCslqhAX/FkqWALgVOCciBpEki/YlYuyU0302k+T9blby3VlFeIDb1tQbJAfVBrn/id8GPgW+H+myouk3z9yDw2EkZxcXSbotIhZUIMZXgJ9ExG1pDK1YfTWxQrLEn++N9GdnkhXxIPlMXsqpMwn4saT9gOei8anZ20bEuZJ+BJxLcgazTfp+cj/zht9fpWmKfS7F/q6rULI2ybMR8XVJW5IsazqaZLW4fLsBHzR8nnyemIvF+Goka5A0tLcxtVmu1PCZha25W4FvSWqXDo4e0fBEehC8gVW7bc4n/VadDmwOjohTgd+RrLtQCROBIyU1fGs+lqS/vqhS8TfympdJxhqOBpC0FbB3Xp3XSbqmbqD4Iki3S2ofyToMD/P5t/6JwIFKF9pJ38/0iPhXgX00WJzG017SGEkDKf65NPp3LWAP4OL0vS0gGVNo2OdioL0Svwf+AWyqZE0VSAa08+NsL2lDSZOBPwGbS+qfxr8hMINkLMNqodaDJn7Ux4NkreIFJN+ETyBZKOg9koWC9iRZVGUJ8CTJwOZ6wJXAM8BU4JL0+TvS/XUgWS3sLyQHyCtIDiTDSA4cT6f1/kgyWPyXtY2pwOvbAD8DZpEcaG4E2qfPTU3bnUVyAJqVbk8tEX+nnLozc9slWXXwz8AjJOMJtwIX5cV0FPBUib/FuWm7M0gGzYfmvX4W8CDJN/kuafmZOZ/VL/L2NxmYQ7Im9QYlPpeif9e8/W6ZvscHSJLaHSSLHgHsAMxLP4tz0rIfp/H9gWScaglwQ/rcnunfdTZwZFrWN/2MH0jf70G1/n/Skh++Kc+sTCR1ioh3c7bvIhmQvian7ABgx4i4rBYxmq0pd0OZlc8VDfcCSOoK7EU6iC2poQvrKJr5Ou/WPHmA26x87gFulPQRSTfWdyPi+fS5gySdAfwxkrELs3WKu6HMzKwkd0OZmVlJThZmZlaSk4WZmZXkZGFmZiU5WZiZWUn/B/Y06wxMUTJ4AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "t = None\n",
    "plt.bar(*p_mother(PYTHIA_val['mothers'], PYTHIA_val['mother_counts'], t=t), label = \"Pythia\", **PLOT_MOTHTER_SETTINGS)\n",
    "plt.bar(*p_mother(JUNIPR_val['mothers'], JUNIPR_val['mother_counts'], t=t), label = \"JUNIPR\", **PLOT_MOTHTER_SETTINGS)\n",
    "plt.xticks(np.arange(0, 11, 1))\n",
    "plt.axis([-0.5, 10.5, 0, 1])\n",
    "plt.xlabel('index in energy sorted state ' + t_label(t))\n",
    "plt.ylabel('probability to branch')\n",
    "plt.legend()\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# P_branch"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## z"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZMAAAEQCAYAAAB1OJkXAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAHfFJREFUeJzt3Xu4VXW97/H3h4uskpTyAkna0pSLeEECpWeTgco66pF82l6OpBEaXhIVPbqTzG1L3W3xkHHqhBJmLnegqR0te1CUQEw9kiK4T6HmJdCjBoIXTAtB+J4/5oAmk7nWmnONOcdksj6v55nPnOM3vuM3vmPMyfoy7ooIzMzM0uhS6wTMzKz+uZiYmVlqLiZmZpaai4mZmaXmYmJmZqm5mJiZWWouJmZmlpqLiZmZpeZiYmZmqXWrdQJZ2X333aOxsbHWaZiZ1ZWnn356TUTs0V5cpykmjY2NLF68uNZpmJnVFUmvlBLn3VxmZpaai4mZmaXmYmJmZqm5mJiZWWouJmZmlpqLiZmZpeZiYmZmqbmYpDRv3jwGDx6MJL70pS8xYsQIBg0axI9+9KOy+lm/fj0jR45EEitWrCgac91113HttddWIGszs8pSZ3kG/NChQ6OtixanzXuh3T4uGd2vaPvChQsZNWoUGzZsoFu3bixbtozDDjuMOXPmMHr06Fb7k8Ty5cvJvzK/WNtmH374IRFBQ0NDu7mamVWCpKcjYmh7cZ3mCvgsDRo0iIMPPpi5c+e2WUzK1aNHj4r1ZWbZKuU/rNXQ2n+CK827uarkvffe47777uNjH/sYEyZMAGDq1KnsvvvuXHrppRx33HEAnHbaaYwcOZLXX399y7Rz5szh+OOP53Of+xz33nsvkNudNmDAAEaOHLkl7pprruGoo47iqKOO4oQTTuCNN97IbgHNzPK4mFTBwoULWb58ObNmzWLTpk00NzcDMGnSJI488khuuOEGHnjgAQB+8YtfsHDhQvr27btl+g8++ID777+f6667jsmTJwMwevToLZ83++QnP8n8+fNZsGABJ598Mpdffnk2C2hmVsC7uSro6KOPZuPGjXTt2pW7776bI444gtGjRzNr1iwmT57MnDlztmyRtOXYY48F4NBDD2X58uWtxu29996MGjWKTZs28d5777F+/fqKLYuZWTlcTCpo/vz5dOu29SodN24cV199NZMnT+auu+7ipptuarefXXbZBcgdI9mwYUPRmBdffJFTTz2Vxx9/nGHDhrFw4ULGjx+fehnMzDrCu7mq7Mtf/jJvvPEGDz74IF26dKFXr14V6Xfp0qXssssuDBs2DKDVomNmlgVvmVRZQ0MDp5xyCmeeeSYzZszYalzPnj3529/+xqxZs2hoaODkk08uud/999+fd955hxdeeIF+/foxd+7cSqduZhU0/NWZNZrz9zOZi7dMUpo3bx4XX3wxkDtm8sgjj2wTM27cONavX7/N8ZKJEydy2mmnMWPGDEaMGEFTUxOQO8Nr7dq1nHbaaQA0NTUxb948pkyZwjPPPMOFF17IkCFDuOKKK2hqauLEE0/k/fffZ+XKlYwbN67KS2xmti1ftJiB5557jh//+MdMnz69JvM3s9p74pbLajLfL3wj3ZZJqRctesukiu644w4++ugjWlpa+PrXv17rdMzMqsbFpIqWLFnCYYcdxpo1azj88MNrnY6ZWdX4AHwVTZ06lalTp9Y6DTOzqvOWiZmZpZbZlomkBmAGMCCZ7xUR8VCRuDHAhUlML6AlIn6UN34ukH/b3Pcj4oRq5m5mZm3LcjdXM7mzx4ZL6gcskjQwIlYVxE0DTomIpZL6AMskvRIRv07Gr4yI8dmlbWZm7clkN5ekLsAE4BaAiHgBWAqcUSR8ekQsTeJWAg8DTVnkaWZmHZPVMZP9gN2A5/PangW2OXc5IqYVNDUAq/OGd5Y0W9Kjkn4l6dCKZ2tmZmXJqpj0Tt7X5rW9C+zZ1kSSdgGGAbfmNb8MfDcivkhuS+cxSZ+uYK4lW7BgwZZH7R566KFMmTKFxsZGdt11162uRP/BD35AY2MjgwcPZt68edx0000MGDCAXr16bfUMkieffJLhw4fTp08fzj33XCB3VX1DQwP9+/enpaWFP/zhDwwfPhxJHHHEEXzxi19k0KBBfOc732HTpk3bjB85ciT9+/dnzJgxrFpVuEfRzKwyMrkCXtI/AY8BPSJifdJ2LfBPEXFUG9PdBPwxIlq9dFzS74F7IuL6IuPOAc4B2GeffT7/yiuvtJ7kw9e1vyCjvt1aDsybN49jjjmG5uZmfvvb3/LYY49tFdPc3MzIkSO3PNyqpaWFCRMmcPzxx3PfffdtiVuxYgXNzc20tLRsaWtsbOTKK6/c8pCtFStWsO+++/Liiy+y//77s2rVKgYOHMj111/P2Wefvc34devWcfjhhzNkyJCt+jWz7PgK+Mp4M3nPv2Vur7z2bSSFYENbhSTxKvDZYiMiYmZEDI2IoXvssUc5+Wbioosu4sEHH2T27Nmp+unduzejRo1q9WaPDQ0NNDU1UavbyZjZji+rYvIy8DbQP6/tQOCpYsGSvgKMAiYlwwck73tKGl8Q3huoy+fVHnLIIVx11VVMmjSJN99sta6WZMOGDXTv3r3V8R999BGf+cxnUs3DzKw1mRSTiNgE3AycBVuKw2BgtqSBkuZL6pqMOxK4CLiA3MH2nsCVSVcfBy6VtHMSO4zcMZU7s1iOarj88svZd999mThxYof7eO6555g/fz4nnXRS0fFvvfUWr732GjfccEOH52Fm1pasrzOZIWlRMt+xEbFSUiO5Cxm7AxuBO4C9gDV5026+r/tK4FfAPEkfATsBX4mIF7NYgLQkbdPWrVs3brvtNoYMGcI999zDkCFDSu7v9NNPp0ePHmzcuJFp06ZxyimnbDO+a9euLFmyhEsuuYRBgwalXgYzs2IyKyYRsQ4YX6R9EdA3b7hvYUxBH/+avLZLPXr0KPos9g8//JAePXoUnebAAw/k6quv5vzzz+f+++8veV6zZ89m//33b3f8jTfeyMSJExk7diyHHHJIyf2bmZXK9+aqsL333psVK1Zs0758+XL23nvvVqe77LLLaGxs5JJLLql4Tt/85jc56KCD+OEPf1jxvs3MwMWk4pqamvjggw+49dZ/XBqzbNkyVq1aRd++rW500bVrV2677TaefPLJiuckicmTJ3P77bezZs2a9icwMyuTi0mFbD4esueee3Lvvfcyffp0Dj74YEaMGMG3v/1tZs78x/Ofb7rpJqZMmcKUKVO2au/fvz/f+973tur36KOPZuXKlUydOnXLRYubH+d7+umnc+edW597UDh+wYIFQO5RwHvttRcjR47kl7/8ZeVXgJl1an5sbwf99a9/5ROf+AQffPABPXv25Omnny7r4LmZdS6+aNGKOu644wBYunQpPXv2pF+/fjXOyMysdvykxQ7aa6+9GD58OOvWrWPmzJn07Nmz1imZmdWMi0kH3XXXXbVOwcxsu+HdXGZmlpqLiZmZpeZiYmZmqbmYmJlZai4mZmaWmouJmZml5mJiZmapuZiYmVlqLiZmZpaai4mZmaXmYmJmZqm5mJiZWWouJmZmlpqLiZmZpeZiYmZmqbmYmJlZai4mZmaWmouJmZml5mJiZmapuZiYmVlqmRUTSQ2SWiQtkrRYUlMrcWMkPSRpgaQlki4qGP8pSfdJeizpa0g2S2BmZq3pluG8mgFFxHBJ/YBFkgZGxKqCuGnAKRGxVFIfYJmkVyLi18n4G4GlEfFdSaOAX0vaPyI+zGxJzMxsK5lsmUjqAkwAbgGIiBeApcAZRcKnR8TSJG4l8DDQlPTzKeCUvH4eBtYDJ1R5EczMrA1Z7ebaD9gNeD6v7VlgaGFgREwraGoAViefhwAfRsSr7fVjZmbZyaqY9E7e1+a1vQvs2dZEknYBhgG35vWztiCs3X7MzKy6sj6bKwqG1U789cA1EfFKG3202o+kc5KD/YtXr15dLMTMzCogq2LyZvLeK6+tV177NiSdA2yIiOkF/exaENpqPxExMyKGRsTQPfbYo/yszcysJFkVk5eBt4H+eW0HAk8VC5b0FWAUMCkZPiAZtQRokLR3Kf2YmVk2MikmEbEJuBk4C7YUh8HAbEkDJc2X1DUZdyRwEXABsLOknsCVST9vAXfn9fMlYCdgThbLYWZmxWV9nckMSYuS+Y6NiJWSGoEBQHdgI3AHsBewJm/aR/I+nw+0SHoM6AqcGBHrqp69mZm1KrNikvzBH1+kfRHQN2+4b2FMQfzbwJcrnZ+ZmXWc781lZmapuZiYmVlqLiZmZpaai4mZmaXmYmJmZqmVXEwkDapmImZmVr/K2TK5L7lI0MzMbCvlFJN3gc9L+o2kf5c0sFpJmZlZfSnnosX/EhFrgB9IOgQ4S9JBwFzgjoho9aaNZma2YyunmPQD1kjaKfk8gNzNGLsBh0vqDlwfEU9XPk0zM9uelVNMbpb0KHAq8DowCzgvIl6HLY/UfQA4ouJZmpnZdq2cYrI3sA4Y3crWx3CgT0WyMjOzulLOAfgpEXFxfiGR9M+SNj8y9wm8VWJm1imVU0yOKdL2MeCXABHxTkSsrEhWZmZWV9rdzSVpn+Tj5icc5j9vfQnbPkbXzMw6mVKOmawAIvn8SsG4vwLTMTOzTq3dYhIRXQAk3R8Rx1c/JTMzqzclHzNprZBIOrNy6ZiZWT1qc8tE0pER8bvk81WthI0Hbq1wXmZmVkfa2811uaQlEfE+8E1yt04p1FD5tMzMrJ60WUwi4r/mDd4QEd8vjJF0WcWzMjOzulLOMZNtCkmie4VyMTOzOtXeMZOfldDHscB1lUnHzMzqUXvHTEYBLe3EfFiZVMzMrF61V0z+PSJubitA0hsVzMfMzOpQm8dM2iskiWL37DIzs06kvWMmM4GLI+Jvkv5cLAToXZXMzMysbrS3m+th4O/J57XAxQXjBUwrZUaSGoAZ5J7Q2A24IiIeaiW2D/AfwBsRMb5g3Fy2vrbl/Yg4oZQczMysOtq7zuSOvMGvRcQfC2Mkfa3EeTUDiojhkvoBiyQNjIhVBf3tB/wEWNNKPysLC4yZmdVWOdeZ/FHSMZJ+KmlO8n50sQJTSFIXYAJwS9LXC8BS4Iwi4e8BY4A/lZqbmZnVVsnFRNLl5HY9AfyR3C6uWZK+VcLk+wG7Ac/ntT0LDC0MjIg1EbGujb52ljRb0qOSfiXp0NKWwMzMqqWcZ8CfAQyKiHc2N0jaDXgE+B/tTLv5IP3avLZ3gQPLmP9mLwM/jYiXJI0BHpPULyL+0oG+zMysAsp5bO9r+YUEICLeAsr5Ix4Fwyoa1VYHEZMj4qXk82/IbeGMKxYr6RxJiyUtXr16dbmzMjOzErV3avA+eYP3S/o34G7gHeBTwGnA/BLm82by3quVz2m8Cny22IiImAnMBBg6dGhhITMzswppbzfXCnJbE/lbEFcUxAQwpZ1+XgbeBvrzjwJyIHB/SVkmJO0JHB8RLXnNvYH/LKcfMzOrrPZ2cz0SEV0joktrL+B37c0kIjYBNwNnAUg6ABgMzJY0UNJ8SV1LyPfjwKWSdk76GQYMA+4sYVozM6uS9rZMmkro43slzqsZmCFpUTLfsRGxUlIjuQsZuwMbk6IyH2gEGiQtBMZFxKvASuBXwDxJHwE7AV+JiBdLzMHMzKqgvYsWN+QPS/occBi5LYTNJlPCWVnJ6b7ji7QvAvrmDW8ERrbRx78mLzMz206UfGqwpLOBS4A+5I5RdAcOAYrds8vMzDqRcq4zGU+ueMyLiFEAkj4NXFuFvMzMKu9hP8evWsopJu9HxEeSdtrcEBF/kTSwCnmZmVXcE39+q9Yp7LDKKSYbkrOn/ijpHuC3wBH4GfBmZp1eOVfAf4vc9SaTyd2McSK5Cw9LvWuwmZntoEreMomIZ/MGx1c+FTMzq1fl3DVYkiZI+q2kZcn7NySVfX8tMzPbsZRzzOQG4Cjgl8BbwO7AJHLXmFxa+dTMzKxelFNMjgYOj4gPNzdI+j7w+4pnZWZmdaWcA/Cv5xcSgIj4O7lbnJiZWSfW3i3oj8wb/J2kn5PbzbX5FvT/DXiieumZmVk9aG8314PkHn6Vf5B9REHMicB3K5mUmZnVl/aKyf0RcVJbAZL+dwXzMTOzOtTmMZNihSQ5RbjP5lOC2ys2Zma24yvnOpOdJf0U+DvwOvB3STM3P6jKzMw6r3LO5vohuUf0HkXuYVZHJ8M/rEJeZmZWR8q5zqRfROSf3fUi8Likdh/ba2ZmO7Zytky6S+qR35AMl1OQzMxsB1ROIbgHWCzpTmA1sAdwKtBShbzMzKyOlHPX4KmS/gKMAz4D/D/guoi4o1rJmZlZfSjnGfBXAa9FRFMV8zEzszpUzjGTbwCPVisRMzOrX+UUk8XAS4WNki6uXDpmZlaPyjkAvwx4QNIcYG1e+znA/6xoVmZmVlfKKSYXAs8A/1zQ3qdy6ZiZWT0qp5jcHhETk3ty7Q6siYiQdH2VcjMzszpRzjGTqyTNIndvrpXk7s31c8DFxMyskyunmLQAHwDHAAcBTcA64D9KmVhSg6QWSYskLZbU6inGyV2JH5LUUmTcpyTdJ+mxpK8hZSyDmZlVQTm7uXpHxJi84efIPX3xyRKnbwYUEcMl9QMWSRoYEavygyTtB/wEWNNKPzcCSyPiu5JGAb+WtH/hI4XNzCw7ZT0DXtLH8xuS28+/kjd8ZrEJJXUBJgC3AETEC8BS4Iwi4e8BY4A/FennU8Apef08DKwHTihjOczMrMLK2TJZAzyTnBq8+RnwTcD85Op4gPHArUWm3Q/YDXg+r+1ZYGhhYESsAUievVVoCPBhRLxapB8/8dHMrEbKKSYnAHOBXskL4PdAz+QF0NDKtL2T9/zrU94FDixj/pv7WVvQ9i6wZ7FgSeeQuw6GffbZp8xZmZlZqcopJjdGxLVtBUj613b6iMJJyph/a3202k9EzARmAgwdOrTYdGZmVgElHzNpr5C0E/Nm8t4rr61XXnup3gR2LWjrSD9mZlZB5RyAT+Nl4G2gf17bgcBTZfazBGiQtHfKfszMrIIyKSYRsQm4GTgLQNIBwGBgtqSBkuZL6lpCP28Bd+f18yVgJ2BOtXI3M7P2ZfnI3WZghqRFyXzHRsRKSY3AAKA7sDEpKvOBRnJbIQuBcXlncJ0PtEh6DOgKnBgR6zJcDjMzK5BZMUn+4I8v0r4I6Js3vBEY2UY/bwNfrnyGZmbWUVkdMzEzsx2Yi4mZmaWW5TETMzMAnrjlslqnYBXmLRMzM0vNxcTMzFJzMTEzs9RcTMzMLDUfgLfiHr6uNvMd9e3azNfMUvGWiZmZpeZiYmZmqXk3lxX1xJ/fqsl8v0CNdq+Bd7GZpeBiYtuVWhUxgC+Mqtmszeqed3OZmVlq3jIpwbR5L9RkvpeM7leT+ZqZlcvFpATDX51Zozl/v0bzNTMrj4uJ2Wa+tsasw1xMtme1+uNmnYd/Y1YhLibbsVqe2dQZ1ex06BqeRebfmFWKz+YyM7PUXEzMzCw17+Yyq7FanXoOMLxmc7YdjbdMzMwsNW+ZmNVY7a5jMqscb5mYmVlqLiZmZpaai4mZmaWW2TETSQ3ADGBAMt8rIuKhVmIvBb6aDP4iIqbmjXsGeDcv/KWImFCdrM3MrBRZHoBvBhQRwyX1AxZJGhgRq/KDJB0LnA0MTpqekfRsRMzZPBwR47NK2szM2pfJbi5JXYAJwC0AEfECsBQ4o0j4ucDtEbEuItYBs4HzssjTzMw6JqtjJvsBuwHP57U9CwwtEjusnbhPS7pX0mOSbpfUWOFczcysTFkVk97J+9q8tneBPVuJbSvuJeC8iBgB/B/gUUkfr2CuZmZWpqzP5oqCYZUY948RERPzjrNMJ3fcZ0yxWEnnSFosafHq1avLTtbMzEqTVTF5M3nvldfWK6+9MLYwrmgliIgAXgM+28r4mRExNCKG7rHHHmUnbWZmpcmqmLwMvA30z2s7EHiqSOxTrcVJOkjS8QXxvYE3KpeqmZmVK5NiEhGbgJuBswAkHUDu1N/ZkgZKmi+paxI+AxgrqSG5NuWrSRvA7sAkSd2Tfk4Cdgbuz2I5zMysuKyvM5khaVEy37ERsTI5G2sA0B3YGBFzJQ0CHk+m+1neNSb/l9yZXo9I2kju2MpxEfF2dothZmaFMismyTUj44u0LwL6FrTdANxQJPZtYFKVUjQzsw7yvbnMzCw1FxMzM0vNxcTMzFJzMTEzs9RcTMzMLDUXEzMzS83FxMzMUnMxMTOz1FxMzMwsNRcTMzNLzcXEzMxSczExM7PUXEzMzCw1FxMzM0vNxcTMzFJzMTEzs9RcTMzMLDUXEzMzS83FxMzMUnMxMTOz1FxMzMwsNRcTMzNLzcXEzMxSczExM7PUXEzMzCw1FxMzM0vNxcTMzFLLrJhIapDUImmRpMWSmtqIvVTS08nrXwrGNUp6WNKjkhZK2rf62ZuZWVu6ZTivZkARMVxSP2CRpIERsSo/SNKxwNnA4KTpGUnPRsScZPgOYGZE3CrpTOBO4PBsFsHMzIrJZMtEUhdgAnALQES8ACwFzigSfi5we0Ssi4h1wGzgvKSfQ8kVmdlJ7GzgYEmfr+4SmJlZW7LazbUfsBvwfF7bs8DQIrHD2ogbBiyPiPUAyftLrfRjZmYZyaqY9E7e1+a1vQvs2Upsa3GF49rqx8zMMpLlMROAKBhWiXHtjSvaj6RzgHOSwfcl/ant9Fq1D/BqB6dty65sWxyrOX2p8aXEtRXTkXG7A2tKyC1rab+javXr7776dozvfsINpf79aq3fz5Y0n4io+gs4gFwR2DOv7X8BvygS+zpwat7wScDK5PPZwLMF8X8Azqty/qur1O/MLKcvNb6UuLZiOjIOWFzN77BW35G/e3/328F3X9Lfr7R5ZbWb62XgbaB/XtuBwFNFYp9qI+4pYF9JOwEk759rpZ9KerdK/f4m4+lLjS8lrq2Yjo7bHlUrX3/3278d5bsv9e9XqryUVKSqkzQF6B0RZ0o6AFgEDAI+CfwYaIqIjcmpwT8AhiSTLgUui+TUYEmLgJsi4jZJXwcuiIhhVc59cUT4IH8VeR13Xv7uqyur9Zv1dSYzkmLQDRgbESslNQIDgO7AxoiYK2kQ8Hgy3c/iH9eYAJwG/EzSBGAjcGoGuc/MYB6dnddx5+XvvroyWb+ZbZmYmdmOy/fmMjOz1DpNMcni3mCSxkh6SNICSUskXVTNZdqeVWp9J+PHSfprskvU6kyZv4U+yb+hlgxTrHulrmNJ4yU9k/zt2vwaUYkcsr7OpJaaqf69waYBp0TEUkl9gGWSXomIX1d52bZHzVRgfScnbvwN6Jld6lZhzZT2W9gP+Anb5zUn27tmSljHiYsjYmGlE+gUWyYZ3htsekQsTeaxEngYaPV/YTuqSq3vxI8j4poqp2xVUuZv4T1gDNDRi4s7pTLXcdV0imJCRvcGi4hpBX01AKvTJl+HKrW+iYjXqpGgZabk30JErEn+Q2HlKeffG8C5kh5JdnGdX6kkOksxyfzeYJJ2IfeH8tYO5FvvKrW+rf6V81uwjilnHa8CHgJGkrus4gJJF1Yiic5STDbL8t5g1wPXRMQrpSS2g6rE+rYdQ6m/Beu4dtdxRDwQEbdGzpvkbmtVka2TzlJM3kzee+W19cprL4wtjFvdyrii/SQ3mNwQEdM7mnCdq9T6tvpXzm/BOibNOn6VUm/k2I7OUkwyuzeYpK8Ao4BJyfABlVmEulKp9W31r5zfgnVMyetY0rcKmnoDb1QiiU5RTCJiE3AzcBZs+QM/GJgtaaCk+ZK6JuEzgLHJedsNwFeTNiLiGeA/gbFJ7FhgWUQ8nfR7JHARcAGws6SewJVZLOP2pFLr2+pfmb8F64Ay1/Fxkr6UxH2c3NmUP69UIp3iRe7MqhZyN5hcTO7GkgDDyd32viEv9lLg6eT1LwX9NAILgEeBhcC+eeNeJ7ffMv+1sNbLXufr+xvJeo6kr6tqvWx+Vee3AHRNvusVwMrk8z61zr8eXmWs4xOA+eQuW3gKmALsVIkcfG8uMzNLrVPs5jIzs+pyMTEzs9RcTMzMLDUXEzMzS83FxMzMUnMxMTOz1FxMzMwsNRcTMzNLzcXEzMxSczExqxFJd0p6P+8VlXq2hFnWfDsVs+2ApP8OfA04OiLernU+ZuXqVusEzDo7SZOAcbiQWB1zMTGrIUkXkLsz8lER8Vat8zHrKBcTsxqRdD5wHjAqItbUOh+zNFxMzGogebTzRHJbJH5MsdU9H4A3qwFJ75J7oNFHec3fjIjKPPXOLGMuJmZmlpqvMzEzs9RcTMzMLDUXEzMzS83FxMzMUnMxMTOz1FxMzMwsNRcTMzNLzcXEzMxSczExM7PU/j92KnA/qnkb+gAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "t = None\n",
    "plt.bar(*p_z(PYTHIA_val['branchings_z'], PYTHIA_val['branchings_counts'], t=t), label = 'Pythia', **PLOT_BRANCH_SETTINGS)\n",
    "plt.bar(*p_z(JUNIPR_val['branchings_z'], JUNIPR_val['branchings_counts'], t=t), label = 'JUNIPR', **PLOT_BRANCH_SETTINGS)\n",
    "plt.xscale('log')\n",
    "plt.xlabel('$z$ ' + t_label(t))\n",
    "plt.ylabel('probability')\n",
    "plt.xticks([0.002, 0.01, 0.1, 0.5], [0.002, 0.01, 0.1, 0.5])\n",
    "plt.legend()\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## theta"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZMAAAERCAYAAAC+ZEqyAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAHeZJREFUeJzt3XmYVNW57/HvTzR0lChxAKNBW6MM4khA8XmI0ip9letwE4erMYegwSEhCfroPaLJMW1yc8RDkJtcB4IxYgJOyY0n5EAwiGLUI0dRzIDm4AB61IDgLIZBeO8ftSFF0UMVe9euLvv3eZ5+qvbaa6/17q7ufntPaykiMDMzS2O7WgdgZmb1z8nEzMxSczIxM7PUnEzMzCw1JxMzM0vNycTMzFJzMjEzs9ScTMzMLDUnEzMzS237WgeQl9133z0aGxtrHYaZWV158sknV0XEHh3V6zLJpLGxkYULF9Y6DDOzuiLppXLq+TSXmZml5mRiZmapOZmYmVlqTiZmZpaak4mZmaXmZGJmZqk5mZiZWWpOJmZmKc2dO5fDDz8cSRx77LEMGzaMgQMH8qMf/aiidtatW8fw4cORxLJly1qtc+211/K9730vg6izpa4yB/zgwYPDDy2aWXsmz13S7vpLR/Rtc938+fNpampi/fr1bL/99ixevJgjjjiCWbNmMWLEiDa3k8TSpUspHqGjtbJN1q5dS0TQ0NDQ4f5kQdKTETG4o3pd5gl4q9CD19am36Yra9OvWcYGDhzIIYccwpw5c9pNJpXq3r17Zm1lyae5zMyq5N1332XmzJl8/OMfZ8yYMQBMnDiR3Xffncsuu4yTTjoJgLPPPpvhw4fz6quvbt521qxZjBw5ks985jPce++9QOF0Wv/+/Rk+fPjmet/97nc57rjjOO644zj55JN57bXX8tvBIk4mZmZVMH/+fJYuXcr06dPZuHEjLS0tAIwbN45jjjmGSZMm8dvf/haAu+66i/nz57P33ntv3n716tXMnj2ba6+9lvHjxwMwYsSIze83+eQnP8m8efN44IEHOOOMM7jiiivy2cESPs1lZpah448/ng0bNtCtWzd+8YtfcNRRRzFixAimT5/O+PHjmTVr1uYjkvaceOKJABx22GEsXbq0zXp9+vShqamJjRs38u6777Ju3brM9qUSTiZmZhmaN28e22+/5Z/WUaNGcc011zB+/Hjuuecebr755g7b2XnnnYHCNZL169e3Wue5557jrLPO4tFHH2XIkCHMnz+f0aNHp96HbeHTXGZmVXbqqafy2muvcd9997HddtvRs2fPTNpdtGgRO++8M0OGDAFoM+nkwUcm1qrHXnyjJv0e3VSTbs2qqqGhgTPPPJPzzjuPKVOmbLGuR48efPDBB0yfPp2GhgbOOOOMsts94IADeOutt1iyZAl9+/Zlzpw5WYdeNh+ZmJmlNHfuXC655BKgcM3koYce2qrOqFGjWLdu3VbXS8aOHcvZZ5/NlClTGDZsGM3NzUDhDq933nmHs88+G4Dm5mbmzp3LhAkTePrpp/nGN77BoEGDuOqqq2hubua0007j/fffZ/ny5YwaNarKe7w1P7RorXrs1str0u/RX/lBTfo1q7Znn32WG264gRtvvLHWoVSk3IcWfWRiZlZFd955Jx9++CHTpk3jy1/+cq3DqZrckomkBknTJC2QtFBScxv1TpH0O0kPSHpK0jdL1u8qaaakR5K2BuWzB2ZmlXvqqac44ogjWLVqFUceeWStw6maPC/At1A4rTZUUl9ggaQBEbGipN5k4MyIWCRpT2CxpJci4tfJ+puARRHxHUlNwK8lHRARa3PbEzOzMk2cOJGJEyfWOoyqy+XIRNJ2wBjgVoCIWAIsAr7USvUbI2JRUm858CDQnLSzK3BmUTsPAuuAk6u8C2Zm1o68TnPtD+wG/KWo7Blgq4s6ETG5pKgBWJm8HwSsjYiXO2rHzMzyk1cy6Z28vlNU9jbQq72NJO0MDAFuK2rnnZJqHbZjZmbVlffdXKX3IauD+tcB342Il9ppo812JF2YXOxfuHLlytaqmJlZBvJKJq8nr8VjCPQsKt+KpAuB9RFRfFP268AuJVXbbCcipkbE4IgYvMcee1QetZmZlSWvZPIC8CbQr6jsIOCJ1ipL+jzQBIxLlg9MVj0FNEjqU047ZmbV9sADD2yeavewww5jwoQJNDY2sssuu2zxJPr1119PY2Mjhx9+OHPnzuXmm2+mf//+9OzZc4s5SB5//HGGDh3KnnvuyUUXXQQUnqpvaGigX79+TJs2jT/96U8MHToUSRx11FF87nOfY+DAgXzrW99i48aNW60fPnw4/fr145RTTmHFitIbaLOR2xPwkiYAvSPivCQ5LAAGAp8EbgCaI2KDpGOAa4AzgE23+94YEV9O2rkLeDYirpF0LDAdODAi1rTXv5+Ar4yfgLcuqaMZRtuZCVQSc+fO5YQTTqClpYX777+fRx55ZIs6LS0tDB8+fPPkVtOmTWPMmDGMHDmSmTNnbq63bNkyWlpamDZt2uayxsZGvv3tb2+eZGvZsmXst99+PPfccxxwwAGsWLGCAQMGcN1113HBBRdstX7NmjUceeSRDBo0aIt2O9IZn4BvASRpAXAncE5y6+8uQH9gh6TencBwYBXwXvK1b1E7XwM+K+kRYAJwWkeJxMyss/rmN7/Jfffdx4wZM1K107t3b5qamtoc7LGhoYHm5maq9U91bg8tJn/wR7dSvgDYu2h579I6JfXfBE7NOj4zs1o49NBDufrqqxk3bhwjRoygV69tvzl1/fr17Ljjjm2u//DDD/n0pz+9ze23x2NzmZnV2BVXXMF+++3H2LFjt7mNZ599lnnz5nH66ae3uv6NN97glVdeYdKkSdvcR3s8n4mZWY6krZ9k2H777bn99tsZNGgQv/rVrxg0qPwhB88991y6d+/Ohg0bmDx5MmeeeeZW67t168ZTTz3FpZdeysCBA1PvQ2ucTMzMMta9e/dW52Jfu3Yt3bt3b3Wbgw46iGuuuYavfe1rzJ49u+y+ZsyYwQEHHNDh+ptuuomxY8dyzjnncOihh5bdfrl8msvMLGN9+vRh2bJlW5UvXbqUPn36bL1B4vLLL6exsZFLL70085i++tWvcvDBB/PDH/4w87bBycTMLHPNzc2sXr2a2267bXPZ4sWLWbFiBXvv3fY9Rt26deP222/n8ccfzzwmSYwfP5477riDVatWZd6+k4mZWUY2XQ/p1asX9957LzfeeCOHHHIIw4YN48orr2Tq1Kmb6958881MmDCBCRMmbFHer18/vv/972/R7vHHH8/y5cuZOHHi5ocWN03ne+6553L33XdvUb90/QMPPAAUpgLea6+9GD58OL/85S+z3XdP22ut8UOLZh177733+MQnPsHq1avp0aMHTz75ZEUXz+tBZ3xo0czsI+Wkk04CYNGiRfTo0YO+ffvWOKLa8d1cZmbbaK+99mLo0KGsWbOGqVOn0qNHj1qHVDNOJmZm2+iee+6pdQidhk9zmZlZak4mZmaWmpOJmZml5mRiZmapOZmYmVlqvpvLOpXJc5fUrO9LR3TdZwTM0vKRiZmZpeZkYmZmqTmZmJlZar5mYp3K0JendlypajzIpNm28pGJmZml5mRiZmapOZmYmVlqTiZmZpaak4mZmaXmZGJmZqk5mZiZWWpOJmZmlpqTiZmZpeZkYmZmqTmZmJlZak4mZmaWmpOJmZml5mRiZmapOZmYmVlqTiZmZpaak4mZmaWW20yLkhqAKUD/pN+rIuJ3bdTdE/gZ8FpEjC5ZNwdoKCp6PyJOrkrQZmZWljyn7W0BFBFDJfUFFkgaEBEriitJ2h/4MbCqjXaWlyYYMzOrrVxOc0naDhgD3AoQEUuARcCXWqn+LnAK8J95xGZmZunldc1kf2A34C9FZc8Ag0srRsSqiFjTTls7SZoh6WFJ/yrpsIxjNTOzCuV1mqt38vpOUdnbwEHb0NYLwE8i4nlJpwCPSOobEX9NG6SZmW2bvO/mipJlVdxAxPiIeD55/xsKRzijWqsr6UJJCyUtXLlyZcXBmplZefI6Mnk9ee3Zxvs0Xgb2bW1FREwFpgIMHjy4NJGZbenBa2vTb9OVtenXLEN5HZm8ALwJ9CsqOwh4opJGJPWSNLqkuDfwWqrozMwslVySSURsBG4BzgeQdCBwODBD0gBJ8yR1K6OpHYHLJO2UtDMEGALcXZ3IzcysHGWf5pI0MCIWp+irBZgiaUHS7zkRsVxSI4UHGXcANiRJZR7QCDRImg+MioiXgeXAvwJzJX0IfAz4fEQ8lyIuMzNLqZJrJjMlnR8RD21LR8ntvqNbKV8A7F20vAEY3k4b/5R8mZlZJ1HJaa63gc9K+o2kf5Y0oFpBmZlZfankyOS/RcQq4HpJhwLnSzoYmAPcGRFZ3JllZmZ1qJJk0hdYJeljyfv+QFPSxpGSdgCui4gnsw/TrPoee/GNmvR7dFNNujXLVCXJ5BZJDwNnAa8C04GLI+JVAEm7Ar8Fjso8SjMz69QqSSZ9gDXAiDaOPoYCe2YSlZmZ1ZVKLsBPiIhLihOJpC9I6pUsPoaPSszMuqRKkskJrZR9HPglQES8FRHLM4nKzMzqSoenuSTtk7xtkNSHLQdnfArYpRqBmZlZ/Sjnmsky/j7a70sl694DbswyIDMzqz8dJpOI2A5A0uyIGFn9kMzMrN6Ufc2krUQi6bzswjEzs3rU7pGJpGMi4vfJ+6vbqDYauC3juMzMrI50dJrrCklPRcT7wFcpDJ1SqiH7sMzMrJ60m0wi4r8XLU6KiB+U1pF0eeZRmZlZXankmslWiSSxQ0axmJlZneromslPy2jjRKBGk2ebmVln0NE1kyZgWgd11mYTipmZ1auOksk/R8Qt7VWQ9FqG8ZiZWR1q95pJR4kk0dqYXWZm1oV0dM1kKnBJRHwg6cXWqgC9qxKZMXnukpr1PbRmPZtZPeroNNeDwN+S9+8Al5SsFzA566DMzKy+dPScyZ1Fi/8QEX8urSPpHzKPyszM6krZMy1GxJ8lnQCcDXwK+CtwZ0TMq1ZwZmZWH8p+aFHSFcDPksU/UzjFNV3SP1YjMDMzqx+VzAH/JWBgRLy1qUDSbsBDwL9kHZiZmdWPSqbtfaU4kQBExBsUTneZmVkX1tGtwfsULc6W9L+BXwBvAbtSuH7iayZmZl1cR6e5llGYsrd43verSuoEMCHDmMzMrM50lEweioim9ipIejDDeKzI0Jen1joEM7OydHTNpLmMNr6fRSBmZla/OnpocX3xsqTPAEcAOxYVjwcOyj40MzOrF2XfGizpAuBSYE/gDxQmxToUaG3MLjMz60IquTV4NIXk8YeIaIqIYUA/YGE1AjMzs/pRSTJ5PyI+BD62qSAi/goMyDwqMzOrK5U8Ab9e0hDgz5J+BdwPHIXngDcz6/IqSSb/CPSgcMF9MjAWeB7wqMFmaTx4be36brqydn3bR0olowY/U7Q4OvtQzMysXlUyarAkjZF0v6TFyetXJKnjrUFSg6RpkhZIWiipzWdYJO0p6XeSprWybldJMyU9krQ1qNx9MDOz6qjkNNck4Djgl8AbwO7AOArPmFxWxvYtgCJiqKS+wAJJAyJiRXElSfsDPwZWtdHOTcCiiPiOpCbg15IOiIi1FeyLWafx2Itv1Kzvo9sd38KsfJUkk+OBI4v/aEv6AfAfHW0oaTtgDPAFgIhYImkRhWHtJ5VUfxc4hcK1mcaSdnYFzgT2S9p5UNI64GTg/1WwL2ZmlqFKbg1+tfS//4j4G7C8jG33B3YD/lJU9gwwuLRiRKyKiDVttDMIWBsRL3fUjpmZ5aejIeiPKVr8vaSfUzjNtWkI+v8JPFZGP72T13eKyt6m8mFYepe0samdXhW2Y2ZmGeroNNd9FCa/Kr7IPqykzmnAd8rsL0qWy7p430EbbbYj6ULgQoB99tmntSpmZpaBjpLJ7Ig4vb0Kksq5VvF68tqzjffleh3YpaSsJ7C4tcoRMRWYCjB48ODWkpCZmWWg3WsmrSWS5BbhPTfdEtxRskm8ALxJYSyvTQ4CnqggVoCngAZJfVK2Y2ZmGarkOZOdJP0E+BvwKvA3SVMl7dTRthGxEbgFOD9p60DgcGCGpAGS5knqVkY7b1CYNnhTO8dSGCtsVrn7YWZm2avkbq4fUrhecRzQn8KtwpGUl6OFwoHNAuBO4JyIWE7htFV/kjG+JHWTNJ/CU/YnSppfMhf914DPSnqEwnTBp7Vz95eZmeWgkudM+kZE8d1dzwGPSvp9ORsnf/BHt1K+ANi7aHkDMLyddt4ETi0vZDMzy0MlRyY7SOpeXJAsV5KQzMzsI6iSRPArYKGku4GVwB7AWcC0KsRlZmZ1pJJRgydK+iswCvg08F/AtRFxZ7WCMzOz+lDJHPBXA69ERJuj/ZqZWddUyTWTrwAPVysQMzOrX5Ukk4UUZlbcgqRLsgvHzMzqUSUX4BcDv5U0iy0HW7wQ+D+ZRmVmZnWlkmTyDeBpkjlJiuyZXThmZlaPKkkmd0TE2GRMrt2BVRERkq6rUmxmZlYnKrlmcrWk6RTG5lpOYWyunwNOJmZmXVwlyWQasBo4ATgYaAbWAD/LPiwzM6snlZzm6h0RpxQtP0th9sXHM47JzMzqTEVzwEvasbggGX7+paLl87IKzMzM6kclRyargKeTW4M3zQHfDMxLno6HwqjAt2UaoZmZdXqVJJOTgTkUpsntmZT9B9Aj+QJoyC40MzOrF5Ukk5si4nvtVZD0TynjMTOzOlT2NZOOEkm5dczM7KPHE1uZdWGT5y6pSb+Xjuhbk36tepxMzLqwoS9PrVHPP6hRv1YtldwabGZm1ionEzMzS83JxMzMUnMyMTOz1JxMzMwsNScTMzNLzcnEzMxSczIxM7PUnEzMzCw1JxMzM0vNycTMzFJzMjEzs9ScTMzMLDUnEzMzS83JxMzMUnMyMTOz1JxMzMwsNScTMzNLLbdpeyU1AFOA/km/V0XE79qoexnwxWTxroiYWLTuaeDtourPR8SY6kRtZmblyHMO+BZAETFUUl9ggaQBEbGiuJKkE4ELgMOToqclPRMRszYtR8TovII2s+w9duvlNen36K947vlqyeU0l6TtgDHArQARsQRYBHypleoXAXdExJqIWAPMAC7OI04zM9s2eV0z2R/YDfhLUdkzwOBW6g7poN6nJN0r6RFJd0hqzDhWMzOrUF7JpHfy+k5R2dtArzbqtlfveeDiiBgG/DvwsKQdM4zVzMwqlPfdXFGyrDLr/X1FxNii6yw3Urjuc0prdSVdKGmhpIUrV66sOFgzMytPXsnk9eS1Z1FZz6Ly0rql9VrNBBERwCvAvm2snxoRgyNi8B577FFx0GZmVp68kskLwJtAv6Kyg4AnWqn7RFv1JB0saWRJ/d7Aa9mFamZmlcolmUTERuAW4HwASQdSuPV3hqQBkuZJ6pZUnwKcI6kheTbli0kZwO7AOEk7JO2cDuwEzM5jP8zMrHV5P2cyRdKCpN9zImJ5cjdWf2AHYENEzJE0EHg02e6nRc+Y/JHCnV4PSdpA4drKSRHxZn67YWZmpXJLJskzI6NbKV8A7F1SNgmY1ErdN4FxVQrRzMy2kcfmMjOz1PI8zWVmVlsPXlu7vpuurF3fOfCRiZmZpeYjkzJMnrukJv0OrUmvZh9dj734Rs36PrqpZl3nwkcmZmaWmo9MyjD05am1DsHMrFPzkYmZmaXmZGJmZqk5mZiZWWpOJmZmlpqTiZmZpea7uczMclCr59UuHdE3l358ZGJmZqn5yMTMLAe1e17tB7n04iMTMzNLzcnEzMxSczIxM7PUnEzMzCw1JxMzM0vNycTMzFJzMjEzs9ScTMzMLDUnEzMzS83JxMzMUnMyMTOz1JxMzMwsNScTMzNLzcnEzMxSczIxM7PUnEzMzCw1JxMzM0vNycTMzFJzMjEzs9ScTMzMLDUnEzMzS83JxMzMUnMyMTOz1HJLJpIaJE2TtEDSQknN7dS9TNKTydf/KlnXKOlBSQ9Lmi9pv+pHb2Zm7dk+x75aAEXEUEl9gQWSBkTEiuJKkk4ELgAOT4qelvRMRMxKlu8EpkbEbZLOA+4GjsxnF8zMrDW5HJlI2g4YA9wKEBFLgEXAl1qpfhFwR0SsiYg1wAzg4qSdwygkmRlJ3RnAIZI+W909MDOz9uR1mmt/YDfgL0VlzwCDW6k7pJ16Q4ClEbEOIHl9vo12zMwsJ3klk97J6ztFZW8Dvdqo21a90nXttWNmZjnJ85oJQJQsq8x6Ha1rtR1JFwIXJovvS/rP9sNr0z7Ay9u4LcAubJ0Es94u67q7A6vKbO+jYls/p2pJ+3NXriz3O21befyuVFK//n9XxkxKG9++ZdWKiKp/AQdSSAK9isr+L3BXK3VfBc4qWj4dWJ68vwB4pqT+n4CLqxz/ypTbT632dlnXBRbm8bPRmb629XOqYjypfu5qsd9p28rjd6WS+h+F35W84svrNNcLwJtAv6Kyg4AnWqn7RDv1ngD2k/QxgOT1M220k6W3U27/mxy2q1bdrqSzfV/S/tyVK8v9TttWHr8rldTvbD8TnZaSzFX9jqQJQO+IOE/SgcACYCDwSeAGoDkiNiS3Bl8PDEo2XQRcHsmtwZIWADdHxO2Svgx8PSKGVDn2hRHRpS7yd8V97mz8GdSHzv455RVf3s+ZTEmSwfbAORGxXFIj0B/YAdgQEXMkDQQeTbb7afz9GROAs4GfShoDbADOyiH2qTn00dl0xX3ubPwZ1IfO/jnlEl9uRyZmZvbR5bG5zMwstS6ZTDIcJ+zpZHywTV8/qX702ajwe7CnpN9JmpZjiB955X4GkpolzU4+gz9K+rGk7nnHa1uTNETS85JG1zqWWsv7OZPOooVsxgl7OiJG5xV0xloo73uwP/BjOvN99PWrhTI+A2AUcH1E3C+pAfgDMA74l3zDtWKSPg+cSed6NqlmutyRSVbjhNWzCr8H7wKnANv6wKe1osLP4LqIuD+ptwb4I9CYT6TWjici4ovAe7UOpD2S9pL0iqTukm6V9JKk95IzKydl1U+XSyZkN04YwKck3SvpEUl3JHem1YOyvwcRsSr5A2bZquQz+NOm95L6AUOBn1c7QGtfRLxS6xjKNBKYQ+FM1H8Bx1J4sv/bwD1Z/d3qiskkq3HCoDDI5MURMQz4d+BhSTtmGGu1VPI9sOqo+DOQ9BjwONASEY9VMTarE8nRxruSNkh6P/naIGmNpBOSaiOB2RGxOiJaImJZRGyMiH8DlgKZjLreFZPJJqnHCYuIsUXnt2+kkPlPySC2vJT7PbDqKfsziIijKQxN9FVJX69qVFYXImIthetn8yKiR0T0oDDiSHNyjW0H4Bhgbum2knoDfYHFWcTSFZPJ68lrz6KynkXlpXVL661srdEoPLDzCuUOilZblXwPrDq26TOIiNeBHwLjqxSX1Z9DKIxRiKSdKAwxtenU6DHAHyJii+s6SZKZAdweEcWnWrdZV0wmmYwTJulgSSNL6vcGXssu1Kqp5Htg1VH2Z5CMfl1sNbBT9UKzOrM5mSTv/xoRbyXLI4HZxZWTmz9+DqwDMjvC7XLJJCI2ArcA5wMk44QdDsyQNEDSPEndkupTgHOS5wEagC8mZVAYdnpckuGRdDqFX/AtPrjOqMLvgVVBhZ/BuGT9pv8oRwPz8o/aOqniZHJY0XsoSSaSROEOwt7A6RGxPqsguvJzJmnHCfsjhTtxHpK0gcK575Mi4s38diOVFsr4HiR/0OZRuBW1QdJ8YFRE5DHPxkddC2V8BhSeJ/mppA8p3IXzZ+CyWgRsf6fCdOGTKPwTMF7SqRHxhZxj2APYg8KdgFD4uVmSrNsP6B4RzxZtcjMwADghIv6WaSwem8vMrD5JOg6YEhF9k+XLKfyj8T8oPNrQPyK+nqzbF1gGrAU+LGrmooiYkToWJxMzs48eSbOBGyIil1PvXe6aiZlZFzEfeDCvznxkYmZmqfnIxMzMUnMyMTOz1JxMzMwsNScTMzNLzcnEzMxSczIxqxFJR0h6VNIHkh6XtE+tYzLbVk4mZjUg6dMUxky6jsIkWS9SmKzIrC45mZjVxiTgloiYmYyRdBeF4S/M6lJXHejRrGYk7QycRmFiok22Azw9stUtJxOz/B1PYUTgPxZGBAegO/DrmkVklpJPc5nlrxGYGRE9N31RGENpTm3DMtt2TiZm+esOfLBpIZl3YjAws2YRmaXkZGKWvyeAYyXtJakPcAfwrTqaWM1sK75mYpa/B4B/ozAj3hvAdRFxS21DMkvHQ9CbmVlqPs1lZmapOZmYmVlqTiZmZpaak4mZmaXmZGJmZqk5mZiZWWpOJmZmlpqTiZmZpeZkYmZmqf1/xQAdpURG9lMAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "t = None\n",
    "plt.bar(*p_theta(PYTHIA_val['branchings_t'], PYTHIA_val['branchings_counts'], t=t), label = 'Pythia', **PLOT_BRANCH_SETTINGS)\n",
    "plt.bar(*p_theta(JUNIPR_val['branchings_t'], JUNIPR_val['branchings_counts'], t=t), label = 'JUNIPR', **PLOT_BRANCH_SETTINGS)\n",
    "plt.xscale('log')\n",
    "plt.xlabel('$\\\\theta$ ' + t_label(t))\n",
    "plt.ylabel('probability')\n",
    "plt.xticks([0.05, 0.1, 0.3, 1, np.pi/2], [0.05, 0.1, 0.3, 1, '$\\pi/2$'])\n",
    "plt.legend()\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## phi"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZMAAAETCAYAAADzrOu5AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAHglJREFUeJzt3Xm4VWX99/H3hwN4nAJFBBMQCRnVUI+oZSAO5JCaZiZihoiomWmaST0/jUYtRbPn+fnjBw6UgZbm1YSShqKmngSESlBxwilQBsUUQYbv88de0Oa49z57s/bggc/rus519rrXd9/r3l4bP2cN91qKCMzMzNJoVesBmJlZy+cwMTOz1BwmZmaWmsPEzMxSc5iYmVlqDhMzM0vNYWJmZqk5TMzMLDWHiZmZpda61gOoll122SW6d+9e62GYmbUYs2fPXhoRHYup3WrCpHv37syaNavWwzAzazEkvVxsrQ9zmZlZag4TMzNLzWFiZmapOUzMzCw1h4mZmaXmMDEzs9QcJmZmlprDxMwshfvvv58BAwYgicGDB3PooYfSv39/fv7zn5fUzwcffMBhhx2GJBYuXJiz5qqrruIHP/hBGUZdftpangHf0NAQnrRoZoVcf/+Cguu/cVSvnO0zZsxgyJAhrFmzhtatWzNv3jz2228/pk6dylFHHZW3P0m89NJLZN+dI1fbBqtXryYiqK+vL+rzpCVpdkQ0FFO71cyAT+XBq2qz3SHfrs12zbZUzf1bXtipmff/Nnf73Bczv2dcDXV19Af26d6RaRO/z1Gtkz9iy/DveZtttkndR6U4TMzMKuCdlav5w2NPc+PvGxl+5ABuAq6542F+csfDfOWz+zP/5TcBOO3YQdS3bc3k75zK7h3bATD1+q8ztfFZnn11KdeedwwnfaY/9896jgv/7x/pvPOOzLj+HAC+/8vpzJj7EgDb1bdhwiUn8fFdPrbpQKr0R2nVzplIqpc0SVKjpFmShhao7SzpPkmTcqzbWdIfJP016Wv/ig7czKxEM+a+yEuL3uJX3zmV9RGM/coRAFz0hU8xaN/ujDv/WO69egQAd1xxGjOuP2djkAC8t2oN91w9gqvOGcqYiX8G4KiGvRgzbPAm29lph22ZPu5sHrhuFKcM2pvLJ0yrzgfMoZp7JmPJnKM5WFIvoFFS34h4I7tIUg/gf4Glefq5EZgTEd+VNAT4vaSeEbG6UgN//MVlleq6oEOo0eE18CE2s81wxKU3s259UNdK3PndYRzUtytHHdCTX90/lzGnD2Zq47McM7B3s/0cfeBeAHyyx268tPitvHVdd23PkEtuYv364J2Vq/lg7dqyfZZSVSVMJLUCRgEnA0TEAklzgDOAcU3K3wGOB8YA3Zv0szPwRWDPpJ8HJX0AfA7IczDTWpRanZ+CmgVocyd9KyXfyWTbfNPHnU3rurpN2s4cuh/f+8UDDD54b26cOovLzjhykz9Qn3zlLRatXL/Je55f9j7vaRmLlv6bNWvXbax/fsm7vPP+Gh5/cRmvvvEWw783hfFjhtFvz848+cyr/PDWaR/64/eQIRX6sE1Ua8+kB9ABeCarbT7woasEImIpZK5oyGF/YHVEvJKjH4fJFqBWe4EAjWtr8z91q57mvl+vvbXjZr1/3r/eSdYvp3XdpmcPdt2tE68uWUHjUwtpJbHjduW5EmvBK2+y/bZt6bdnZwDWrltXln43V7XCZMMlEiuy2t4G+m1GPyuatL0N7JqrWNJoYDRAt27dStxU7dXyf6zV+mvG4OBXJtRoy9fWaLtbl23atGZIQy9+fOs0LvvykZus226bNqz6YC3THp9P2zatObyh+L3FLru259/vreaVxcvp1nlnGp9aWOaRl6baV3M1ndSSc/ejxD7y9hMRE4AJkJlnshnbsq1I7f6nbi3ZE/MW8t93PQzA18fdyTknfor9enfdpOaYQ/rx0JPPccjee27SfvKQAVz5v39i+23b8qPzT+Ci6+4C4MoJU7n+4i9w5YSpAFx03V0M/2wDt937BMtXrOS6KdO55PQjOPO4g7j4+t/Ss2tHOnxse5aveI/v33wvV559TBU++aaqFSZvJr/b53ldSj/tmrS1B+Zt/tDMzDJO6bCw5PcM7N+dgf27F6xpt/22HHFgb1q33vR8yldPGcRXTxm0cfmGS07ZZP3E75y+yfIdPxy5yfLoz3+a0Z//9Mbly8/MP0Gy0qoVJi8Ay4He/CdA+gH3lNjPk0C9pK4R8WpWP5PKMUjLUssT4VYVtTrxD1vPyf/7/vY0hzf0Zupj8zj2kP61Hk5FVWWeSUSsByYCIwEk7QUMACZL6itpuqS6Qn0k/SwD7szqZzDQFphaqbGbmW2uBa+8yYjv38aKd9+nX4/daj2ciqr2PJPxkhqT7Q6LiMWSugN9gDbAuiRUppO5LLhe0gzgzKwruL4KTJL0V6AOODEiVlXxc2wVanny37Z8tdorOrjK2/vaFwfztS8Obr5wC1C1MEn+hz8iR3sjsHvW8jrgsAL9LAdOKP8Izcxsc/kW9GZmlprDxMzMUnOYmJlZar4FvdlWqpaTNBu7ja7Ztq0yHCZmZinMevoVbv3j48xZ8Bo9u3TkyIG9+d1D/+DfK1fxmQE9N85Gv/2+Wdw5fQ47brcNF5wyiNfefJvfTH+S5StWMvkHI+jYfgcA5r+4iOvveJBFS1fwmQE9ufzMo7jw2jv55/Ov07nDx/jysQPps0cnrv7l/cx7cRH99uxMm9Z1vPNeZnujP/9pXvrX0o3rB054hG233ZZFixbRq1cvbrrpJjp1auYhYJvBj+0twuM3f7PMozGzj6IuK2YXXP9auwPyrvvUqHHccMkpHNhvD276/WPMevplxo8ZtknNTb9/jP17d2X/PpnbrUx99Cmu/sV9HLLPnvz0wpM21i1auoKb//A4/zXy6I1tJ18+kRHHHcQJg/bdWPOFMTfxmx+NpEunnVi+4j2GXTGJr57yGU4ctO/G9c899xw9e/Zk1apVDBw4kP33359JkyYV9d+jlMf2+pyJmVkNnXLE/vxt3sv8ufHpVP3s3G57DujTlb/lueFjfX09Q4cOZXP/qG6Ow8TMrIZ6dtmFsz53MD+74wGWv7MyVV9r162jri7//9bXrl1Lly5dUm0jH4eJmVmNnXHMQHbbpR3jJk/f7D4W/msZs55+hSEH7JVz/bJly3jttdcYN67p8wjLwyfgzcyqKNdz/1rXteKKkcdw1vdvY8bs5+i9R85HNOU09qZ7aNO6jvXrg69/aQiHN2z6WODhw4dTV1fHk08+yTe+8Q3696/MDScdJmZmZda2TR1r1q7/UPuateto0zr3PW33/HgHzj7xU1w7+S+Mu+jkorc1dtSxdOm0U971kydPpmfPntx4441ccMEFDBs2jH333bfo/ovlw1xmZmW26047smhZ04fCwr+WrqDTzvkfDXz6ZxvYbZd23PDrGWUf0/nnn8/ee+/NDTfcUPa+wWFiZlZ2A/t3Z9XqNfzpr09tbHvx9aW89c5KOu6UP0zqWrXiv0YezfyXFpd9TJIYM2YMU6ZMYenSpWXv32FiZlYmG86H7Pyx7bjqghO5+8G5nPHdX3De1bcz/u5HNnkS4t0PzuW2e5/gtnuf4HcP/WNj+x6dd+bckz69Sb8XXnsny1e8x5Q/z2Lqo0/xwmtLNj7Sd+xN9/CXJ57ZpD57/fDhw3nggQcAOO200/j4xz/OYYcdxl133VXez+5Ji83zpEUzy+W9VR+wfX1b3l+9hiMu+Dm3XnEGvfco/+zyNA45+9rNfq8nLZqZVcGlP/stAAteeYPttmlD1wInwrd0vprLzGwz7dJ+B8758RQ+WLOWb515FNvVt631kGrGYWJmtpl+eN7xtR7CR4YPc5mZWWoOEzMzS81hYmZmqTlMzMwsNYeJmZml5jAxM7PUHCZmZpaaw8TMzFJzmJiZWWoOEzMzS81hYmZmqVUtTCTVS5okqVHSLElDC9ReKml28nNZk3UnJO0PS5op6djKj97MzAqp5o0ex5J5fsrBknoBjZL6RsQb2UWSjgbOAQYkTXMlzY+IqZI6AL8GToiI+yUdDvxRUteIWF7Fz2JmZlmqsmciqRUwCrgZICIWAHOAM3KUnwtMiYhVEbEKmAycl6zrAdQDjyfLjwPbAb0rN3ozM2tOtQ5z9QA6ANnPlpwP5HqC14EF6uYBLwMnJMvHA+80qTczsyqr1mGuDc+xXJHV9jbQL09t07pdASJipaQTgN9JugLYATg6It7KtVFJo4HRAN26dUv1AczMLL9qX83V9IHzKrIuUyztAkwDzouIvsDZwERJO+bsJGJCRDREREPHjh03d8xmZtaMaoXJm8nv9llt7bPam9Y2rVuSvP4i8HZE3AeQ/A7g9LKO1szMSlKtMHkBWM6mJ8r7ATNz1M4sUNcWWNOkfg3wsfIM08zMNkdVwiQi1gMTgZEAkvYic+nvZEl9JU2XVJeUjweGJfNS6snsdYxP1j0I9JG0T9LPvsDewIxqfA4zM8ut2vNMxktqTLY7LCIWS+oO9AHaAOsiYpqk/sCjyftuiYipABHxD0lnAZMkrSRzAn50ROTawzEzsyqpWpgkc0ZG5GhvBHZv0jYOGJennynAlAoM0czMNpPvzWVmZqk5TMzMLDWHiZmZpeYwMTOz1BwmZmaWmsPEzMxSc5iYmVlqDhMzM0vNYWJmZqk5TMzMLDWHiZmZpeYwMTOz1BwmZmaWmsPEzMxSKzpMkmeMmJmZfUgpeyZ/kDS4YiMxM7MWq5QweRs4QNIfJf1YUt9KDcrMzFqWUp60+NmIWApclzx7faSkvYFpwO0R8WZFRmhmZh95pYRJL2CppLbJ6z7AkKSPgZLaAD+JiNnlH6aZmX2UlRImEyU9ApwKvA78CjgvIl4HkLQzcC9wUNlHaWZmH2mlhElXYBVwVJ69j4OBzmUZlZmZtSilnIC/OiIuzg4SSSdL2jVZfBzvlZiZbZVKCZMjc7RtC9wFEBFvRcTisozKzMxalGYPc0nqlrysl9QVUNbqJ4F2lRiYmZm1HMWcM1kIRPL65Sbr/g38dzkHZGZmLU+zYRIRrQAk3RMRx1Z+SGZm1tIUfc4kX5BIOqt8wzEzs5ao4J6JpEER8XDy+so8ZSOAW8s8LjMza0GaO8x1uaQnI+Jd4Hwyt05pqr78wzIzs5akYJhExHFZi+Mi4tqmNZK+WcyGJNUD48nchqU18J2IuC9P7aXA6cniHRFxTda6tsBVZCZJ1gPvAaf43mBmZrVT9Az4XEGSaFNkF2MBRcTBknoBjZL6RsQb2UWSjgbOAQYkTXMlzY+Iqcny9cD8iLg0qf8fYPtiP4eZmZVfc+dMbimij6PJ7CkU6qcVMAo4GSAiFkiaA5wBjGtSfi4wJSJWJe+dDJwHTJXUCTgeuHBDcUScX8QYzcysgpq7mmsImbklhX5WF7GdHkAH4JmstvlAQ47aAwvUDSYz7+W7kh6VdJ+kzxSxfTMzq6DmDnP9OCImFiqQ9K8ittMp+b0iq+1toF+e2qZ1G+7/1R04BLg7Ij4taRBwn6Q+EdF0QqWZmVVJwT2T5oIkkeueXXm7bLKsnFUfrttgG2Adyaz75LLlOcDwXMWSRkuaJWnWkiVLShimmZmVorlzJhOAiyNipaQXc5Xwn72OQjZcadU+z+umte2zltsDG5LgLWB5RKzJWv8a0CXXRiNiAjABoKGhIV9AmZlZSs0d5noQeD95vQK4uMl6kbm6qjkvAMuB3vwnQPoB9+SonZnUkVU3M3k9F9hZUuuIWJu0dQT+UcQYzMysQpo7zHV7RGz4i/7LEfFQk58ZwJeb20hErAcmAiMBJO1F5tLfyZL6SpouqS4pHw8Mk1SfzE05PWkDeAx4nsxVYEjqTeaE/R3Ff2QzMyu3UuaZPCXpSOA0YDdgEXB7REwvsouxwHhJjcl2h0XEYkndyUxkbAOsi4hpkvoDjybvu2XDHJOIWC/pJGCCpHOT9V+MiOeL/RxmZlZ+RYeJpMuBi8gcmnoK2AX4laTrI+Knzb0/mTcyIkd7I7B7k7ZxfHj+yYZ1z5G5ZNnMzD4iSnkG/BlA/4h4a0ODpA7AQ0CzYWJmZluuUh7b+1p2kABExDIyh7vMzGwr1tylwd2yFu+R9EPgTjKX6O5M5vxJsedMzMxsC9XcYa6FZCYQZk8u/E6TmgCuLuOYzMyshWnuMNdDEVEXEa3y/QAPV2OgZmb20dVcmAwtoo8flWMgZmbWcjX3cKzs25Yg6RPAfsB2Wc1jyH3DRjMz20qUMs/kHOAbQGfg72QmGe4L5Lpnl5mZbUVKuTR4BJnw+HtEDImIQ8ncQ2tWJQZmZmYtRylh8m5yc8W2GxoiYhHQt+yjMjOzFqWUGfBrJB0IPCXpbuAvwEEU/wx4MzPbQpWyZ/ItMvNNxgDvABeQedZIs3cNNjOzLVspdw2en7U4ovxDMTOzlqroPRNljJL0F0nzkt9nS8r36F0zM9tKlHLOZBxwOHAXsIzMLegvIjPH5NLyD83MzFqKUsLkCGBgRKze0CDpWuBvZR+VmZm1KKWcgH89O0gAIuJ9YHF5h2RmZi1Nc7egH5S1+LCk28gc5tpwC/ovAY9XbnhmZtYSNHeY689kHn6VfZL90CY1JwLfLeegzMysZWkuTO6JiC8UKpD02zKOx8zMWqCC50xyBUlyiXDnDZcENxc2Zma25Stlnsn2km4C3gdeB96XNEHS9hUbnZmZtQilXM11A5lH9B4O9CFzqXAk7WZmthUrZZ5Jr4jIvrrrOeBRSX5sr5nZVq6UPZM2krbJbkiWSwkkMzPbApUSBHcDsyT9GlgCdAROBSZVYFxmZtaClHLX4GskLQLOBLoArwJXRcTtlRqcmZm1DKU8A/5K4LWIGFrB8ZiZWQtUyjmTs4FHKjUQMzNruUoJk1nA800bJV1czJsl1UuaJKlR0ixJefdwJF0qaXbyc1memrMkRdGjNzOziinlBPw84F5JU4EVWe2jgZ8V8f6xgCLiYEm9gEZJfSPijewiSUcD5wADkqa5kuZHxNSsmnr8DBUzs4+MUvZMLgS2AU4Gzsr66dzcGyW1AkYBNwNExAJgDnBGjvJzgSkRsSoiVgGTgfNyjOWOEsZuZmYVVEqYTImIIWRmwJ8KHJ4sTyzivT2ADsAzWW3zgYYctQcWqpPUHjgOmFLC2M3MrIJKCZMrJf2KzL25FpO5N9dtwE+KeG+n5Hf24bG3gV3z1Baq+zZwDbC+yHGbmVmFlRImk4D3gCOBvYGhwCrglyX00fSEuXJWfbguUyx1AQZknz8pRNLo5GT/rCVLlpQwTDMzK0UpJ+A7RcTxWctPk3n64hNFvPfN5Hf7PK+b1rbPWm5PZsY9ZE7if6/YAUfEBGACQENDg6/8MjOrkFLC5HVJ20XEyg0Nye3nX85aPisibs3x3heA5UBv/hMg/YB7ctTOTOrIqpuZvN4f6Jk8SqU+2eYM4KGI8NMezcxqpJQwWUrmMt2p/OcZ8EOB6cnseIARwIfCJCLWS5oIjAQekbQXmUt/h0vqC/w/YGhErAPGA9dJuiZ5++nAN5N+9t/Qp6TuwEsRcVgJn8HMzCqglDD5HDCNzGGnDYeh/gbskPxAsreQx1hgvKTGZLvDImJxEgp9gDbAuoiYJqk/8GjyvluaniORNJpMMG3YM7kuIv5QwmcxM7MyKiVMboyIHxQqkHRFvnXJnJEROdobgd2btI0DxhXoa+O5EDMzq72ir+ZqLkiKrTEzsy1PKZcGm5mZ5eQwMTOz1BwmZmaWmsPEzMxSc5iYmVlqDhMzM0vNYWJmZqk5TMzMLDWHiZmZpeYwMTOz1BwmZmaWmsPEzMxSc5iYmVlqDhMzM0vNYWJmZqk5TMzMLDWHiZmZpeYwMTOz1BwmZmaWmsPEzMxSc5iYmVlqDhMzM0vNYWJmZqk5TMzMLDWHiZmZpeYwMTOz1BwmZmaWmsPEzMxSq1qYSKqXNElSo6RZkoYWqL1U0uzk57Ks9o6Srpc0Q9Kjku6R1KM6n8DMzPJpXcVtjQUUEQdL6gU0SuobEW9kF0k6GjgHGJA0zZU0PyKmAscBXYHDI2K9pB8BvwX2q9qnMDOzD6nKnomkVsAo4GaAiFgAzAHOyFF+LjAlIlZFxCpgMnBesm4B8NOIWJ8s3w4MkNSxkuM3M7PCqnWYqwfQAXgmq20+0JCj9sB8dRHxWEQ8kbWuHlgFvFvW0ZqZWUmqFSadkt8rstreBnbNU1tMHcDngFsi4v3UIzQzs81WzXMmANFkWUXWfYikTwAnAEMK1IwGRgN069atyCGamVmpqrVn8mbyu31WW/us9qa1TeuWZBdI2onM+ZdhEZG9F7OJiJgQEQ0R0dCxo0+rmJlVSrXC5AVgOdA7q60fMDNH7cxCdZK2BX4NfCsinpW0q6R25R+ymZkVqyphklx9NREYCSBpLzKX/k6W1FfSdEl1Sfl4YFgyL6UeOD1pI6m5ncxeyXxJOwAn4UuDzcxqqtrzTMZLaky2OywiFkvqDvQB2gDrImKapP7Ao8n7bknmmEAmjE5MfrLlPW9iZmaVV7UwSeaMjMjR3gjs3qRtHDAuR+1EMns4Zmb2EeJ7c5mZWWoOEzMzS81hYmZmqTlMzMwsNYeJmZml5jAxM7PUHCZmZpaaw8TMzFJzmJiZWWoOEzMzS81hYmZmqTlMzMwsNYeJmZml5jAxM7PUHCZmZpaaw8TMzFJzmJiZWWoOEzMzS81hYmZmqTlMzMwsNYeJmZml5jAxM7PUHCZmZpaaw8TMzFJzmJiZWWoOEzMzS81hYmZmqTlMzMwsNYeJmZmlVrUwkVQvaZKkRkmzJA0tUHuppNnJz2VN1nWX9KCkRyTNkLRn5UdvZmaFtK7itsYCioiDJfUCGiX1jYg3soskHQ2cAwxImuZKmh8RU5Pl24EJEXGrpLOAXwMDq/MRzMwsl6rsmUhqBYwCbgaIiAXAHOCMHOXnAlMiYlVErAImA+cl/XySTMhMTmonA/tIOqCyn8DMzAqp1mGuHkAH4JmstvlAQ47aAwvUHQi8FBEfACS/n8/Tj5mZVUm1wqRT8ntFVtvbwK55avPVNV1XqB8zM6uSap4zAYgmyyqyrrl1OfuRNBoYnSy+K+nZwsPLaxdg6Wa+16w5/n5Z5Ywal+b7tUexhdUKkzeT3+3zvG5a2z5ruT2wJM+6Qv0QEROACZsx3k1ImhURPpRmFeHvl1VStb5f1TrM9QKwHOid1dYPmJmjdmaBupnAnpLaAiS/P5GnHzMzq5KqhElErAcmAiMBJO1FclWWpL6SpkuqS8rHA8OSeSn1wOlJGxExF/g7MCypHQbMi4jZ1fgcZmaWW7XnmYyX1Jhsd1hELJbUHegDtAHWRcQ0Sf2BR5P33ZI1xwTgNOAWSaOAdcCpVRh76kNlZgX4+2WVVJXvlyIKnes2MzNrnu/NZWZmqTlMCijlfmJmZlszh0lhY0nuJ0bmQoA7JHUq/BYzs+qStI2kmyW9LOnfkuZKOqaaY3CY5FHi/cTMzGqpNfAqMBhoB/wX8JvkAqeqcJjkV8r9xMyKlvwV+Y6kdZLeTX7WSVol6chaj89anoh4LyLGRsTCiFgfEX8CXgIOqNb3zWGSXyn3EzMrWkSsBi4CpkfEDhGxA5mJvUMj4i+1HZ1tCZLD8b3IzMOryvfNYdK8Yu8nZlaKfYB/AkjansydHP5Z0xHZFkFSGzKP5/hFRGw4slLx75vDJL/se4iR9TrnfcDMSrTxH3fyelFEvFXD8dgWIDnXexvwAfC1rFUV/745TPIr5X5iZqXK/sf9SbxXYilJEpkLhjoBX4iINVmrK/59c5jkUeh+YrUcl7V8kjoCHclc0AGZ2wktqN2IbAvxP0Bf4PiIeH9DY7W+bw6TwsaSCfxGMs+eHxYRi2s7JNsC7AO8kPUP/nXgVEkH1XBM1oJJ2oPMI88HAIuzrtoaTpW+b743l5mZpeY9EzMzS81hYmZmqTlMzMwsNYeJmZml5jAxM7PUHCZmZpaaw8TMzFJzmJiZWWoOE7Mak3SjpCtrPQ6zNBwmZrV3KPDXWg/CLA2HiVkNSGol6VuSFpK5d9Idki6p8bDMNpvDxKw2vg0cB3ydzB1cPwdcI6lrTUdltpkcJmZVJmlH4P8AXwG6AnMi4gngVTK3BzdrcRwmZtV3OPBsRCwk86CiOckT8nbCT/K0FsphYlZ9nck8xRMyz5+YAwwCVgD/qNWgzNJoXesBmG2FngEOkNQD2JtMsEwCvhV+wJC1UH44llmVJc/qvpbMI6HbAc8BP4yI22o6MLMUHCZmNSLpS8BZEXF0rcdilpbPmZjVTm/g2VoPwqwcHCZmtdObzBwTsxbPh7nMzCw175mYmVlqDhMzM0vNYWJmZqk5TMzMLDWHiZmZpeYwMTOz1BwmZmaWmsPEzMxS+/96ArnnJOSzKgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "t = None\n",
    "plt.bar(*p_phi(PYTHIA_val['branchings_p'], PYTHIA_val['branchings_counts'], t=t), label = 'Pythia', **PLOT_BRANCH_SETTINGS)\n",
    "plt.bar(*p_phi(JUNIPR_val['branchings_p'], JUNIPR_val['branchings_counts'], t=t), label = 'JUNIPR', **PLOT_BRANCH_SETTINGS)\n",
    "plt.xlabel('$\\\\phi$ ' + t_label(t))\n",
    "plt.ylabel('probability')\n",
    "plt.xticks([0, np.pi, 2* np.pi], [0,'$\\pi$', '$2\\pi$'])\n",
    "plt.legend()\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## delta"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZMAAAEQCAYAAAB1OJkXAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAHsJJREFUeJzt3Xu0VXW99/H3h4vskpRULmnY9sZFUpFAaQxTQNmPekSH+WiQRqhopeVl6Ek0j23zlHhIfeqIEh5zexLJy8myAQdFFFOPpAj2mKiogT5qIGjipRCE7/PHmuBisdZea+25LuzN5zXGGmut3/zO3/zOORb7y7z85lREYGZmlkaneidgZmbtn4uJmZml5mJiZmapuZiYmVlqLiZmZpaai4mZmaXmYmJmZqm5mJiZWWouJmZmllqXWi1IUgMwDRiQLPeyiHggT9wY4HtJTA+gJSJ+njV9DtCQNcsHEXFcseXvtttu0djYmGodzMy2N08//fTqiOhZLK5mxQRoBhQRwyX1AxZIGhgRK3PirgdOjojFkvoAz0l6NSJ+l0xfERETyl14Y2MjCxcuTJO/mdl2R9KrpcTV5DCXpE7AROAWgIhYCiwGTssTPjUiFidxK4CHgaZa5GlmZm1Tq3MmewO7Ai9ktS0BhuYGRsT1OU0NwKqs7ztKmiHpUUm/lXRQxbM1M7Oy1KqY9E7e12S1vQv0am0mSTsBw4Bbs5pfAX4YEV8hs6fzmKTPVTBXMzMrU62v5sq9372KxF8D/CgiNh+zi4hJEfFy8vn3ZPZwxuebWdLZkhZKWrhq1ap8IWZmVgG1KiZvJe89stp6ZLVvRdLZwPqImFqk79eAL+SbEBHTI2JoRAzt2bPoxQhmZtZGtSomrwDvAP2z2vYHnsoXLOlEYCRwfvJ9v+S9l6QJOeG9gTcrnK+ZmZWhJsUkIjYCNwNnwObiMBiYIWmgpHmSOifTDgfOA75L5mR7d+DypKtPAxdJ2jGJHUbmnMqdtVgPMzPLr5bnTJoBSVoAzATGJZf+7kxmIGPXJG4mMAJYDbyfvDYdxloB/BaYK+kPwL8DJ0bESzVaBzOzrcydO5fBgwcjiSOOOILDDjuMQYMG8fOf/7z4zFnWrVvHiBEjkMTy5cvzxlx99dVcddVVFci6srS9PAN+6NCh4UGLZtaa6+cubXX6haP7FZw2f/58Ro4cyfr16+nSpQvPPfccBx98MLNmzWL06NEF55PEsmXLyL5DR762TT766CMigoaGhq2mVYOkpyNiq2EcuWo5At7MtjHF/nhWS2t/lDuKQYMGccABBzBnzpxWi0m5unXrVrG+Ksk3ejQzq5L33nuP++67j0996lNMnDgRgClTprDbbrtx0UUXccwxxwAwduxYRowYwRtvvLF53lmzZnHssceyzz77cO+99wKZw2kDBgxgxIgRm+N+9KMfMWrUKEaNGsVxxx3Hm2/W53ok75mYbceGvza9Tkv+aZ2WWzvz589n2bJlPP744xx++OE0NzcDcP755/PEE09w7bXXAplDWr/+9a+3OqT14YcfMnv2bO666y4mTZrEiSeeyOjRo5k0aRItLS2b4z772c8yb948JNHS0sIll1zCr371qxqt5SdcTMzMKujII49kw4YNdO7cmbvvvptDDz2U0aNHc/vttzNp0iRmzZq1eY+kNUcffTQABx10EMuWLSsY17dvX0aOHMnGjRt57733WLduXcXWpRwuJmZmFTRv3jy6dNnyT+v48eO58sormTRpEnfddRc33XRT0X522mknIHOOZP369XljXnrpJU455RQef/xxhg0bxvz585kwYULqdWgLnzMxM6uy448/njfffJP777+fTp060aNHj+IzlWDx4sXstNNODBs2DKBg0akFFxMzsypraGjg5JNP5vTTT+drX/vaFtO6d+/O3//+d26//Xbuueeesvrdd999+dvf/sbSpZmr8ubMmVOxnMvlYmJmltLcuXO54IILgMw5k0ceeWSrmPHjx7Nu3bqtzpece+65jB07lmnTpnHYYYfR1JR5fNPYsWNZs2YNY8eOBaCpqYm5c+cyefJknnnmGb73ve8xZMgQLrvsMpqamjjhhBP44IMPWLFiBePH5733bVV50KLZduyJWy6uy3K/fGbHv5or1/PPP88NN9zA1KnF7l27bSl10KL3TMzMqmjmzJl8/PHHtLS08M1vfrPe6VSNi4mZWRUtWrSIgw8+mNWrV3PIIYfUO52q8aXBZmZVNGXKFKZMmVLvNKrOeyZmZpaai4mZmaXmYmJmZqm5mJiZWWouJmZmlpqLiZlZCg899NDmR+0edNBBTJ48mcbGRnbeeectRqJfd911NDY2MnjwYObOnctNN93EgAED6NGjxxbPIHnyyScZPnw4ffr04Vvf+haQGVXf0NBA//79aWlp4dlnn2X48OFI4tBDD+UrX/kKgwYN4gc/+AEbN27cavqIESPo378/Y8aMYeXKlVXZDh4Bb7Yd8wj4HA9f3fr0kZcWnCSJuXPnctRRR9Hc3MyDDz7IY489tkVMc3MzI0aM2Pxwq5aWFiZOnMixxx7Lfffdtzlu+fLlNDc3b/HcksbGRi6//PLND9lavnw5e+21Fy+99BL77rsvK1euZODAgVxzzTWcddZZW01fu3YthxxyCEOGDNmi32I8At7MrB0477zzuP/++5kxY0aqfnr37s3IkSML3uyxoaGBpqYmqvWfahcTM7M6OvDAA7niiis4//zzeeutt1L1tX79erp27Vpw+scff8znP//5VMsoxMXEzKzOLrnkEvbaay/OPffcNvfx/PPPM2/ePE466aS8099++21ef/31zY8LrjTfTsXMrIYkbdXWpUsXbrvtNoYMGcJvfvMbhgwZUnJ/p556Kt26dWPDhg1cf/31nHzyyVtN79y5M4sWLeLCCy9k0KBBqdchHxcTM7MK69atW95nsX/00Ud069Yt7zz7778/V155Jeeccw6zZ88ueVkzZsxg3333LTr9xhtv5Nxzz2XcuHEceOCBJfdfKh/mMjOrsL59+7J8+fKt2pctW0bfvn0LznfxxRfT2NjIhRdeWPGcvvOd7/DFL36Rn/3sZxXvG1xMzMwqrqmpiQ8//JBbb711c9tzzz3HypUr2WOPPQrO17lzZ2677TaefPLJiuckiUmTJnHHHXewevXqivfvYmJmViGbzof06tWLe++9l6lTp3LAAQdw2GGHcemllzJ9+vTNsTfddBOTJ09m8uTJW7T379+fH//4x1v0e+SRR7JixQqmTJmyedDipsf5nnrqqdx5551bxOdOf+ihh4DMo4B33313RowYUfbz5ouuuwctmm2/PGgxnffff5/PfOYzfPjhh3Tv3p2nn366rJPn7YEHLZqZVdkxxxwDwOLFi+nevTv9+vWrc0b146u5zMzaaPfdd2f48OGsXbuW6dOn071793qnVDcuJmZmbXTXXXfVO4Vthg9zmZlZajUrJpIaJLVIWiBpoaSmAnFjJD0g6SFJiySdlzN9F0n3SXos6atjne0yM2uHanmYq5nM1WPDJfUDFkgaGBG5N9e/Hjg5IhZL6gM8J+nViPhdMv1GYHFE/FDSSOB3kvaNiI9qtiZmZraFmuyZSOoETARuAYiIpcBi4LQ84VMjYnEStwJ4GGhK+tkFODmrn4eBdcBxVV4FMzNrRa0Oc+0N7Aq8kNW2BNjq2uWIuD6nqQFYlXweAnwUEa8V68fMzGqnVsWkd/K+JqvtXaBXazNJ2gkYBmy6J0HvnD5a7UfS2cn5mYWrVq3KF2JmZhVQ66u5cofbb30v5i1dA/woIl5tpY+C/UTE9IgYGhFDe/bsWUaaZmZWjloVk02PD+uR1dYjq30rks4G1kfE1Jx+ds4JbbUfMzOrvloVk1eAd4D+WW37A0/lC5Z0IjASOD/5vl8yaRHQICn7Hs4F+zEzs9qoSTGJiI3AzcAZsLk4DAZmSBooaZ6kzsm0w4HzgO8CO0rqDlye9PM2cHdWP0cAOwCzarEeZmaWX63HmUyTtCBZ7riIWCGpERgAdAU2ADOB3YHsG+4/kvX5HKBF0mNAZ+CEiFhb9ezNzKygmhWT5A/+hDztC4A9sr4XfnJMZvo7wPGVzs/MzNrO9+YyM7PUXEzMzCw1FxMzM0vNxcTMzFJzMTEzs9RcTMzMLDUXEzMzS83FxMzMUnMxMTOz1FxMzMwsNRcTMzNLzcXEzMxSczExM7PUXEzMzCw1FxMzM0vNxcTMzFJzMTEzs9RcTMzMLDUXEzMzS83FxMzMUnMxMTOz1FxMzMwsNRcTMzNLzcXEzMxSK7mYSBpUzUTMzKz9KmfP5D5JR1QtEzMza7fKKSbvAl+S9HtJP5E0sFpJmZlZ+9KljNj/FRGrgeskHQicIemLwBxgZkS8VZUMzcxsm1dOMekHrJa0Q/J5ADAy6eMQSV2BayLi6cqnaWZm27JyisnNkh4FTgHeAG4Hvh0RbwBI2gX4b+DQimdpZh3Lw1fXZ7kjL63PcrcD5RSTvsBaYHSBvY/hQJ+KZGVmHdoTf3m7Lsv98si6LHa7UE4xmRwRP8lukPRV4LHkfMkTtLJXIqkBmEbm8FgX4LKIeKBAbB/gP4E3I2JCzrQ5QENW0wcRcVwZ62FmZhVWTjE5CvhJTtungHuAwyPib0XmbwYUEcMl9QMWSBoYESuzgyTtDfwCWF2gnxW5BcbMzOqr6KXBkvaUtCfQIKnvpu9J2yJg5xL66ARMBG4BiIilwGLgtDzh7wFjgBdLXw0zM6unUvZMlgORfH41Z9r7wNQS+tgb2BV4IattCTA0NzC5/BhJhfraUdIMYE/gbeCHEfGnEnIwM7MqKVpMIqITgKTZEXFsG5fTO3lfk9X2LrB/G/p6BfiPiHhZ0hjgMUn9IuKvbczNzMxSKnkEfKFCIun0MpYXOd8L7n60ksekiHg5+fx7Mns44wvkdrakhZIWrlq1qtxFmZlZiVrdM5F0eET8Ifl8RYGwCcCtRZazaXR8jwKf03gN+EK+CRExHZgOMHTo0NxCZmZmFVJsz+QSSd2Tz98B9srzaigwb7ZXgHeA/llt+wNPlZOspF6SJuQ09wbeLKcfMzOrrFb3TCLin7K+XhsRP82NkXRxsYVExEZJNwNnAI9K2g8YDJya3DDyBqApIjYU6erTwEWS7o6IDyUNA4YBZxbLwczMqqfkcSb5Ckmia4ldNAPTJC1IljsuIlZIaiQzkLErsEFSZ2Ae0EjmcuT5wPiIeA1YAfwWmCvpY2AH4MSIeKnU9TAzs8ords7klyX0cTRQ9EY7EbGWzPmV3PYFwB5Z3zcAI1rp41+Sl5mZbSOK7ZmMBFqKxHxUmVTMzKy9KlZMfhIRN7cWIMknv83SqNcddM0qqNgJ+FYLSeIooJQ4M8ujXnfQNaukYudMpgMXRMTfJf0lXwifjG43M7PtVLHDXA8D/0g+rwEuyJku4PpKJ2VmZu1LscNcM7O+fiMi/pwbI+kbFc/KzMzalXLGmfxZ0lHAWOBzwF+BmRExr1rJmZlZ+1DyjR4lXULm6YcAfyZziOt2Sd+vRmJmZtZ+lPOkxdOAQdlPVJS0K/AI8G+VTszMzNqPkvdMgNdzH80bEW+TOdxlZmbbsWKXBu+Z9XW2pH8F7gb+BuxC5vyJz5mYmW3nih3mWk7mgVbZD7G6LCcmgMkVzMnMzNqZYoe5HomIzhHRqdAL+EMtEjUzs21XsWLSVEIfP65EImZm1n4VG7S4Pvu7pH2Ag8k8pGqTSWSemmhmZtupki8NlnQWcCHQB/gTmYdZHQjku2eXtXPXz11al+VeOLpfXZZrZumUM85kApniMTciRgJI+hxwVRXysjob/tr0Oi250AM9zWxbVs44kw8iYtOjcgGIiL8CAyuelZmZtSvl7JmslzQM+LOk3wAPAodS+jPgzcysgypnz+T7ZMabTALeA84FegC+a7CZ2XaunLsGL8n6OqHyqZiZWXtVzl2DJWmipAclPZe8nylJxec2M7OOrJxzJtcCo4B7gLeB3YDzyYwxuajyqZmZWXtRTjE5EjgkIj7a1CDpp8AfK56VmZm1K+UUkzeyCwlARPxD0ooK52Tbs4evrt+yR15av2WbtXPFbkF/eNbXP0j6FZnDXJtuQf814InqpWfbmyf+8nbdlv3lkXVbtFm7V2zP5H4yD7/KPsl+WE7MCcAPK5mUmZm1L8WKyeyIOKm1AEn/VcF8zMysHWr10uB8hSS5RLjPpkuCixUbMzPr+MoZZ7KjpP8A/gG8AfxD0nRJO1YtOzMzaxfKuZ3Kz8g8oncUMIDMpcKRtJuZ2XasnEuD+0VE9tVdLwGPS/Jje83MtnPlFJOukrrlDFrsVmofkhqAaWT2aroAl0XEAwVi+wD/CbwZERNypu0CtJC5NLkLcE5ELCpjPczyq+cYF7N2rpxi8htgoaQ7gVVAT+AUMn/YS9EMKCKGS+oHLJA0MCJWZgdJ2hv4BbC6QD83Aosj4oeSRgK/k7Rv7oBKs3LVc4yLWXtX8jmTiJgCXAMcTuaeXF8Bro6I64rNK6kTMBG4JelrKbAYOC1P+HvAGODFPP3sApyc1c/DwDrguFLXw8zMKq+cZ8BfAbweEU1tWM7ewK7AC1ltS4ChuYERsTpZXr5+hgAfRcRrefrxeBczszop52quM4FH27ic3sn7mqy2d4FebehnTU5bW/oxM7MKKqeYLARezm2UdEEZfUTu7GXMW6iPgv1IOlvSQkkLV61a1YZFmZlZKco5Af8c8N+SZrHl3sHZwP8pMu9byXuPAp9L9Rawc05bjyS3rUTEdGA6wNChQ/MVITMzq4Byisn3gGeAr+a09ylh3leAd4D+fFJA9gdml7F8gEVAg6S+EfH/svppKbMfMzOroHIOc90RESPJjIA/BRiVfL+52IwRsTGJOwNA0n7AYGCGpIGS5knqXEI/bwN3Z/VzBLADMKuM9TAzsworZ8/kCkm3A/8b6Aqsl3Q3mcuES9EMTJO0IFnuuIhYIamRzEDGrsCGpKjMAxrJ7IXMB8ZnXcF1DtAi6TGgM3BCRKwtYz3MzKzCyikmLcCbwFFkngHfE/gGmZHqRcd5JH/wJ+RpXwDskfV9AzCilX7eAY4vI28zM6uycopJ74gYk/X9eTJPX3yywjmZmVk7U845kzckfTq7Ibn9/KtZ30+vVGJmZtZ+lLNnshp4Jrk0eNMz4JuAecnoeMgcxrq1ohmamdk2r5xichwwh8y4jh5J2x+B7skLoKFyqZmZWXtRTjG5MSKuai1A0r+kzMfMzNqhcu4a3GohKTXGzMw6nnJOwJuZmeXlYmJmZqm5mJiZWWouJmZmlpqLiZmZpeZiYmZmqbmYmJlZai4mZmaWmouJmZml5mJiZmapuZiYmVlqLiZmZpaai4mZmaXmYmJmZqm5mJiZWWouJmZmlpqLiZmZpeZiYmZmqbmYmJlZal3qnYAV9sQtF9c7BTOzknjPxMzMUnMxMTOz1FxMzMwsNRcTMzNLzcXEzMxSq9nVXJIagGnAgGS5l0XEAwViLwK+nnz9dURMyZr2DPBuVvjLETGxOlmbmVkpanlpcDOgiBguqR+wQNLAiFiZHSTpaOAsYHDS9IykJRExa9P3iJhQq6TNzKy4mhzmktQJmAjcAhARS4HFwGl5wr8F3BERayNiLTAD+HYt8jQzs7ap1TmTvYFdgRey2pYAQ/PEDisS9zlJ90p6TNIdkhornKuZmZWpVsWkd/K+JqvtXaBXgdjW4l4Gvh0RhwH/Azwq6dMVzNXMzMpU66u5Iue7Soz7ZELEuVnnWaaSOe8zJl+spLMlLZS0cNWqVWUna2ZmpalVMXkree+R1dYjqz03NjcubyWIiABeB75QYPr0iBgaEUN79uxZdtJmZlaaWhWTV4B3gP5ZbfsDT+WJfapQnKQvSjo2J7438GblUjUzs3LVpJhExEbgZuAMAEn7kbn0d4akgZLmSeqchE8DxklqSMamfD1pA9gNOF9S16Sfk4Adgdm1WA8zM8uv1uNMpklakCx3XESsSK7GGgB0BTZExBxJg4DHk/l+mTXG5P+SudLrEUkbyJxbOSYi3qndapiZWa6aFZNkzMiEPO0LgD1y2q4Frs0T+w5wfpVSNDOzNvK9uczMLDUXEzMzS82P7TWz7UY9H4X95TN/Wrdl14L3TMzMLDUXEzMzS83FxMzMUnMxMTOz1FxMzMwsNRcTMzNLzcXEzMxSczExM7PUXEzMzCw1FxMzM0vNxcTMzFJzMTEzs9RcTMzMLDUXEzMzS83FxMzMUnMxMTOz1FxMzMwsNRcTMzNLzcXEzMxSczExM7PUXEzMzCw1FxMzM0vNxcTMzFJzMTEzs9RcTMzMLDUXEzMzS61LvRNoD5645eJ6p2Bmtk3znomZmaXmYmJmZqnVrJhIapDUImmBpIWSmlqJvUjS08nrn3OmNUp6WNKjkuZL2qv62ZuZWWtqec6kGVBEDJfUD1ggaWBErMwOknQ0cBYwOGl6RtKSiJiVfJ8JTI+IWyWdDtwJHFKbVTAzs3xqUkwkdQImAl8FiIilkhYDpwHX5oR/C7gjItYm884Avg3MknQQmSIzI4mdAdwo6UsR8XT118TMrG3qdSHPl8/8aU2WU6vDXHsDuwIvZLUtAYbmiR3WStwwYFlErANI3l8u0I+ZmdVIrYpJ7+R9TVbbu0CvArGF4nKntdaPmZnVSK3HmUTOd5UYV2xa3n4knQ2cnXz9QNKLradX0J7Aa22cd2e2LoBtjSsUU057blup+VVamuVWapu2Nr2S23Q3YHWrmaZXi+1ZSmy527QtbbXYnoXyqPS81f93P/Fa2PpvWDnr9oWSoiKi6i9gPzJFoFdW278Dv84T+wZwStb3k4AVyeezgCU58c8C365y/qtSzDu9UnGFYsppz20rNb8qbNM2L7dS27S16ZXcpsDCjrA9q7FN29JWi+1Zq21aw3/3q8pdbrmvWh3megV4B+if1bY/8FSe2KdaiXsK2EvSDgDJ+z4F+qmkd1PM+/sKxhWKKac9t63U/CotzXIrtU1bm97etmkttmcpseVu0zRt1bYt/EZbiymnPfdvWMW3p5IqVXWSJgO9I+J0SfsBC4BBwGeBG4CmiNiQXBp8HTAkmXUxcHEklwZLWgDcFBG3Sfom8N2IGFbl3BdGhE/yW5v491NZ3p7lq8U2q/U4k2lJMegCjIuIFZIagQFAV2BDRMyRNAh4PJnvl/HJGBOAscAvJU0ENgCn1CD36TVYhnVc/v1Ulrdn+aq+zWq2Z2JmZh2X781lZmapddhiUqt7gRWaV1JPSdcn8zwuabakvSu/plZrlfptJdPHS3o/OdxriTK3cR9JD0hqqWGK25RSt5ekz0j6laSnJD2ZzNO9Ejl05OeZNFPle4EVmfefgL7AqIjYKOnHwH8BB1drha1mmqnAbyu5KOXvQEX+MXcwzZS2jfcGfkFtxp1sy5opYXsB/0Jm3MhwMsM15gGXA5PSJtAh90yy7gV2C2TuBUbmqrDT8oRvvhdYZO4HtuleYBS4F9gBkr5UbF5gKfBvEbEx+T4TGCypZ+XW1GqtUr+txA0R8aMqp9zulLmN3wPGAG0dkNzulbm9BgFPRsSG5G/Tk1ToP7gdsphQu3uBFZw3Iv4nIp7MmtYArAU+aMP62LajUr8tIuL1aiTYAZS8jSNidVKot2fl/CZnA6Mk7Sjp08Ao4I+VSKKjFpNa3QustXlzHUfmMud/tJq5besq9duywsrZxlbG9oqIqcB84C/AcmARUJG9445aTDapxb3Ail5bLWkf4HjgsmKx1m5U4rdlrSt1G1tG0e0l6QfAQWTOm+xJZq9mQiUW3lGLyVvJe4+sth5Z7bmxuXGrCkzL7ae1eQGQ9FkyxzLHRUQ9bqholVWp35YVVs42tvK213nAL7LO400DrqpEEh21mNTqXmCtzYukT5G5+uv7EfGipF6Sdm7rStk2oVK/LSusnG1s5W2vHYD1Wd/XA5+pRBIdspgkVyncDJwBkNwLbDAwQ9JASfMkdU7CpwHjkuu0G4CvJ21ExDPAn4BxSew44Ln45KmOBedN+p9JZq9kSXIt94n40uB2rVK/LSuszG283Stzez0IfE0JMr/JhyuVSId8kbl6qoXMDSUXkrmRJGSur34DaMiKvQh4Onn9c04/jcBDwKNkTlztlTM977xkxhdEnteIem8bv7aZ39aZyW8qkr6uqPe6bSuvUrcx0DnZhsuBFcnnPeud/za8vXqROVryRzKXBd8DfK4SOfjeXGZmllqHPMxlZma15WJiZmapuZiYmVlqLiZmZpaai4mZmaXmYmJmZqm5mJiZWWouJmZ1Iuk0SUskrZH0tqSL652TWVt50KJZHST3aFsFfCkinpW0I7BHZB5sZNbuuJiY1YGkrmRuxPcemVuA3O5CYu2ZD3OZ1ccYMo/0PZzM0+/mSzqgzjmZtZn3TMxqTNIXgFnAQRGxIWmbCbwYEc31zM2srbxnYlZ7pwFPbCokiZ2AD+uUj1lqLiZmtbcnmWd0AyBpF+AIYE7dMjJLycXErPaWAsdJ6impFzAD+G1EPFvnvMzazMXErPamAYuAF8k8oGgRMLGuGZml5BPwZmaWmvdMzMwsNRcTMzNLzcXEzMxSczExM7PUXEzMzCw1FxMzM0vNxcTMzFJzMTEzs9RcTMzMLLX/D4zcmoaEA2WiAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "t = None\n",
    "plt.bar(*p_delta(PYTHIA_val['branchings_d'], PYTHIA_val['branchings_counts'], t=t), label = 'Pythia', **PLOT_BRANCH_SETTINGS)\n",
    "plt.bar(*p_delta(JUNIPR_val['branchings_d'], JUNIPR_val['branchings_counts'], t=t), label = 'JUNIPR', **PLOT_BRANCH_SETTINGS)\n",
    "plt.xscale('log')\n",
    "plt.xlabel('$\\\\delta$ ' + t_label(t))\n",
    "plt.ylabel('probability')\n",
    "plt.xticks([2e-4, 0.01, 0.1, 0.8], [2e-4, 0.01, 0.1, 0.8])\n",
    "plt.legend()\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "tf2",
   "language": "python",
   "name": "tf2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}