Skip to content

Instantly share code, notes, and snippets.

@DavidMStraub

DavidMStraub/2017 LHCb Bqtautau limits.ipynb

Created March 13, 2017 07:51
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save DavidMStraub/63c9b4ecc66225b070bb961c837e9981 to your computer and use it in GitHub Desktop.
Save DavidMStraub/63c9b4ecc66225b070bb961c837e9981 to your computer and use it in GitHub Desktop.
Download ZIP
Display the source blob
Display the rendered blob
Raw
2017 LHCb Bqtautau limits.ipynb
{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Flavio implementation of LHCb 2017 $B_q\\to\\tau\\tau$ limit likelihoods"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"See [arXiv:1703.02508](https://arxiv.org/abs/1703.02508)"
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": false,
"deletable": true,
"editable": true,
"scrolled": false
},
"outputs": [
{
"data": {
"text/plain": [
"'0.20'"
]
},
"execution_count": 1,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"import flavio\n",
"flavio.__version__"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"collapsed": false,
"deletable": true,
"editable": true
},
"outputs": [],
"source": [
"m = flavio.Measurement['LHCb B->tautau 2017']"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"collapsed": false,
"deletable": true,
"editable": true
},
"outputs": [],
"source": [
"_, dist_s = m._parameters['BR(Bs->tautau)']\n",
"_, dist_d = m._parameters['BR(B0->tautau)']"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"collapsed": true,
"deletable": true,
"editable": true
},
"outputs": [],
"source": [
"def pos_ppf(dist, cl):\n",
" \"\"\"Return the PPF (inverse CDF) for the distribution restricted to positive BRs\"\"\"\n",
" a = dist.cdf(0)\n",
" return dist.ppf(a + cl*(1-a)).item()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### 90% and 95% CL limits for $B_s$ and $B_d$"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {
"collapsed": false,
"deletable": true,
"editable": true
},
"outputs": [
{
"data": {
"text/plain": [
"(0.005779399903578894, 0.0070950609668250595)"
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"pos_ppf(dist_s, 0.9), pos_ppf(dist_s, 0.95)"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"collapsed": false,
"deletable": true,
"editable": true
},
"outputs": [
{
"data": {
"text/plain": [
"(0.0016471999725188003, 0.002022179538456282)"
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"pos_ppf(dist_d, 0.9), pos_ppf(dist_d, 0.95)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Plots of the PDFs"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {
"collapsed": true,
"deletable": true,
"editable": true
},
"outputs": [],
"source": [
"import matplotlib.pyplot as plt\n",
"import flavio.plots\n",
"%matplotlib inline\n",
"%config InlineBackend.figure_format = 'svg'"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.legend.Legend at 0x7fc3d5795b38>"
]
},
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/svg+xml": [
"<?xml version=\"1.0\" encoding=\"utf-8\" standalone=\"no\"?>\n",
"<!DOCTYPE svg PUBLIC \"-//W3C//DTD SVG 1.1//EN\"\n",
" \"http://www.w3.org/Graphics/SVG/1.1/DTD/svg11.dtd\">\n",
"<!-- Created with matplotlib (http://matplotlib.org/) -->\n",
"<svg height=\"267pt\" version=\"1.1\" viewBox=\"0 0 391 267\" width=\"391pt\" xmlns=\"http://www.w3.org/2000/svg\" xmlns:xlink=\"http://www.w3.org/1999/xlink\">\n",
" <defs>\n",
" <style type=\"text/css\">\n",
"*{stroke-linecap:butt;stroke-linejoin:round;}\n",
" </style>\n",
" </defs>\n",
" <g id=\"figure_1\">\n",
" <g id=\"patch_1\">\n",
" <path d=\"M 0 267.593582 \n",
"L 391.041911 267.593582 \n",
"L 391.041911 0 \n",
"L 0 0 \n",
"z\n",
"\" style=\"fill:none;\"/>\n",
" </g>\n",
" <g id=\"axes_1\">\n",
" <g id=\"patch_2\">\n",
" <path d=\"M 45.541911 228.14 \n",
"L 380.341911 228.14 \n",
"L 380.341911 10.7 \n",
"L 45.541911 10.7 \n",
"z\n",
"\" style=\"fill:#ffffff;\"/>\n",
" </g>\n",
" <g id=\"PolyCollection_1\">\n",
" <path clip-path=\"url(#pd8980b1118)\" d=\"M 63.809828 186.67406 \n",
"L 63.809828 218.256364 \n",
"L 64.419775 218.256364 \n",
"L 65.029723 218.256364 \n",
"L 65.63967 218.256364 \n",
"L 66.249617 218.256364 \n",
"L 66.859564 218.256364 \n",
"L 67.469511 218.256364 \n",
"L 68.079458 218.256364 \n",
"L 68.689406 218.256364 \n",
"L 69.299353 218.256364 \n",
"L 69.9093 218.256364 \n",
"L 70.519247 218.256364 \n",
"L 71.129194 218.256364 \n",
"L 71.739141 218.256364 \n",
"L 72.349089 218.256364 \n",
"L 72.959036 218.256364 \n",
"L 73.568983 218.256364 \n",
"L 74.17893 218.256364 \n",
"L 74.788877 218.256364 \n",
"L 75.398824 218.256364 \n",
"L 76.008772 218.256364 \n",
"L 76.618719 218.256364 \n",
"L 77.228666 218.256364 \n",
"L 77.838613 218.256364 \n",
"L 78.44856 218.256364 \n",
"L 79.058507 218.256364 \n",
"L 79.668455 218.256364 \n",
"L 80.278402 218.256364 \n",
"L 80.888349 218.256364 \n",
"L 81.498296 218.256364 \n",
"L 82.108243 218.256364 \n",
"L 82.71819 218.256364 \n",
"L 83.328138 218.256364 \n",
"L 83.938085 218.256364 \n",
"L 84.548032 218.256364 \n",
"L 85.157979 218.256364 \n",
"L 85.767926 218.256364 \n",
"L 86.377873 218.256364 \n",
"L 86.987821 218.256364 \n",
"L 87.597768 218.256364 \n",
"L 88.207715 218.256364 \n",
"L 88.817662 218.256364 \n",
"L 89.427609 218.256364 \n",
"L 90.037556 218.256364 \n",
"L 90.647504 218.256364 \n",
"L 91.257451 218.256364 \n",
"L 91.867398 218.256364 \n",
"L 92.477345 218.256364 \n",
"L 93.087292 218.256364 \n",
"L 93.697239 218.256364 \n",
"L 94.307187 218.256364 \n",
"L 94.917134 218.256364 \n",
"L 95.527081 218.256364 \n",
"L 96.137028 218.256364 \n",
"L 96.746975 218.256364 \n",
"L 97.356922 218.256364 \n",
"L 97.96687 218.256364 \n",
"L 98.576817 218.256364 \n",
"L 99.186764 218.256364 \n",
"L 99.796711 218.256364 \n",
"L 100.406658 218.256364 \n",
"L 101.016605 218.256364 \n",
"L 101.626553 218.256364 \n",
"L 102.2365 218.256364 \n",
"L 102.846447 218.256364 \n",
"L 103.456394 218.256364 \n",
"L 104.066341 218.256364 \n",
"L 104.676288 218.256364 \n",
"L 105.286236 218.256364 \n",
"L 105.896183 218.256364 \n",
"L 106.50613 218.256364 \n",
"L 107.116077 218.256364 \n",
"L 107.726024 218.256364 \n",
"L 108.335971 218.256364 \n",
"L 108.945919 218.256364 \n",
"L 109.555866 218.256364 \n",
"L 110.165813 218.256364 \n",
"L 110.77576 218.256364 \n",
"L 111.385707 218.256364 \n",
"L 111.995654 218.256364 \n",
"L 112.605602 218.256364 \n",
"L 113.215549 218.256364 \n",
"L 113.825496 218.256364 \n",
"L 114.435443 218.256364 \n",
"L 115.04539 218.256364 \n",
"L 115.655337 218.256364 \n",
"L 116.265285 218.256364 \n",
"L 116.875232 218.256364 \n",
"L 117.485179 218.256364 \n",
"L 118.095126 218.256364 \n",
"L 118.705073 218.256364 \n",
"L 119.31502 218.256364 \n",
"L 119.924968 218.256364 \n",
"L 120.534915 218.256364 \n",
"L 121.144862 218.256364 \n",
"L 121.754809 218.256364 \n",
"L 122.364756 218.256364 \n",
"L 122.974703 218.256364 \n",
"L 123.584651 218.256364 \n",
"L 124.194598 218.256364 \n",
"L 124.804545 218.256364 \n",
"L 125.414492 218.256364 \n",
"L 126.024439 218.256364 \n",
"L 126.634386 218.256364 \n",
"L 127.244334 218.256364 \n",
"L 127.854281 218.256364 \n",
"L 128.464228 218.256364 \n",
"L 129.074175 218.256364 \n",
"L 129.684122 218.256364 \n",
"L 130.294069 218.256364 \n",
"L 130.904017 218.256364 \n",
"L 131.513964 218.256364 \n",
"L 132.123911 218.256364 \n",
"L 132.733858 218.256364 \n",
"L 133.343805 218.256364 \n",
"L 133.953752 218.256364 \n",
"L 134.5637 218.256364 \n",
"L 135.173647 218.256364 \n",
"L 135.783594 218.256364 \n",
"L 136.393541 218.256364 \n",
"L 137.003488 218.256364 \n",
"L 137.613435 218.256364 \n",
"L 138.223383 218.256364 \n",
"L 138.83333 218.256364 \n",
"L 139.443277 218.256364 \n",
"L 140.053224 218.256364 \n",
"L 140.663171 218.256364 \n",
"L 141.273118 218.256364 \n",
"L 141.883066 218.256364 \n",
"L 142.493013 218.256364 \n",
"L 143.10296 218.256364 \n",
"L 143.712907 218.256364 \n",
"L 144.322854 218.256364 \n",
"L 144.932801 218.256364 \n",
"L 145.542749 218.256364 \n",
"L 146.152696 218.256364 \n",
"L 146.762643 218.256364 \n",
"L 147.37259 218.256364 \n",
"L 147.982537 218.256364 \n",
"L 148.592484 218.256364 \n",
"L 149.202432 218.256364 \n",
"L 149.812379 218.256364 \n",
"L 150.422326 218.256364 \n",
"L 151.032273 218.256364 \n",
"L 151.64222 218.256364 \n",
"L 152.252167 218.256364 \n",
"L 152.862115 218.256364 \n",
"L 153.472062 218.256364 \n",
"L 154.082009 218.256364 \n",
"L 154.691956 218.256364 \n",
"L 155.301903 218.256364 \n",
"L 155.91185 218.256364 \n",
"L 156.521798 218.256364 \n",
"L 157.131745 218.256364 \n",
"L 157.741692 218.256364 \n",
"L 158.351639 218.256364 \n",
"L 158.961586 218.256364 \n",
"L 159.571533 218.256364 \n",
"L 160.181481 218.256364 \n",
"L 160.791428 218.256364 \n",
"L 161.401375 218.256364 \n",
"L 162.011322 218.256364 \n",
"L 162.621269 218.256364 \n",
"L 163.231216 218.256364 \n",
"L 163.841164 218.256364 \n",
"L 164.451111 218.256364 \n",
"L 165.061058 218.256364 \n",
"L 165.671005 218.256364 \n",
"L 166.280952 218.256364 \n",
"L 166.890899 218.256364 \n",
"L 167.500847 218.256364 \n",
"L 168.110794 218.256364 \n",
"L 168.720741 218.256364 \n",
"L 169.330688 218.256364 \n",
"L 169.940635 218.256364 \n",
"L 170.550582 218.256364 \n",
"L 171.16053 218.256364 \n",
"L 171.770477 218.256364 \n",
"L 172.380424 218.256364 \n",
"L 172.990371 218.256364 \n",
"L 173.600318 218.256364 \n",
"L 174.210265 218.256364 \n",
"L 174.820213 218.256364 \n",
"L 175.43016 218.256364 \n",
"L 176.040107 218.256364 \n",
"L 176.650054 218.256364 \n",
"L 177.260001 218.256364 \n",
"L 177.869948 218.256364 \n",
"L 178.479896 218.256364 \n",
"L 179.089843 218.256364 \n",
"L 179.69979 218.256364 \n",
"L 180.309737 218.256364 \n",
"L 180.919684 218.256364 \n",
"L 181.529631 218.256364 \n",
"L 182.139579 218.256364 \n",
"L 182.749526 218.256364 \n",
"L 183.359473 218.256364 \n",
"L 183.96942 218.256364 \n",
"L 184.579367 218.256364 \n",
"L 185.189314 218.256364 \n",
"L 185.799262 218.256364 \n",
"L 186.409209 218.256364 \n",
"L 187.019156 218.256364 \n",
"L 187.629103 218.256364 \n",
"L 188.23905 218.256364 \n",
"L 188.848997 218.256364 \n",
"L 189.458945 218.256364 \n",
"L 190.068892 218.256364 \n",
"L 190.678839 218.256364 \n",
"L 191.288786 218.256364 \n",
"L 191.898733 218.256364 \n",
"L 192.50868 218.256364 \n",
"L 193.118628 218.256364 \n",
"L 193.728575 218.256364 \n",
"L 194.338522 218.256364 \n",
"L 194.948469 218.256364 \n",
"L 195.558416 218.256364 \n",
"L 196.168363 218.256364 \n",
"L 196.778311 218.256364 \n",
"L 197.388258 218.256364 \n",
"L 197.998205 218.256364 \n",
"L 198.608152 218.256364 \n",
"L 199.218099 218.256364 \n",
"L 199.218099 186.722765 \n",
"L 199.218099 186.722765 \n",
"L 198.608152 186.444403 \n",
"L 197.998205 186.16604 \n",
"L 197.388258 185.887678 \n",
"L 196.778311 185.609315 \n",
"L 196.168363 185.330953 \n",
"L 195.558416 185.05259 \n",
"L 194.948469 184.774228 \n",
"L 194.338522 184.495865 \n",
"L 193.728575 184.217503 \n",
"L 193.118628 183.93914 \n",
"L 192.50868 183.660778 \n",
"L 191.898733 183.381621 \n",
"L 191.288786 183.095664 \n",
"L 190.678839 182.809707 \n",
"L 190.068892 182.52375 \n",
"L 189.458945 182.237793 \n",
"L 188.848997 181.951836 \n",
"L 188.23905 181.665879 \n",
"L 187.629103 181.379922 \n",
"L 187.019156 181.093965 \n",
"L 186.409209 180.808008 \n",
"L 185.799262 180.522051 \n",
"L 185.189314 180.236094 \n",
"L 184.579367 179.950137 \n",
"L 183.96942 179.66418 \n",
"L 183.359473 179.378223 \n",
"L 182.749526 179.093571 \n",
"L 182.139579 178.809795 \n",
"L 181.529631 178.526019 \n",
"L 180.919684 178.242244 \n",
"L 180.309737 177.958468 \n",
"L 179.69979 177.674692 \n",
"L 179.089843 177.390917 \n",
"L 178.479896 177.107141 \n",
"L 177.869948 176.823365 \n",
"L 177.260001 176.539589 \n",
"L 176.650054 176.255814 \n",
"L 176.040107 175.972038 \n",
"L 175.43016 175.688262 \n",
"L 174.820213 175.404486 \n",
"L 174.210265 175.122024 \n",
"L 173.600318 174.852619 \n",
"L 172.990371 174.583214 \n",
"L 172.380424 174.313809 \n",
"L 171.770477 174.044404 \n",
"L 171.16053 173.774999 \n",
"L 170.550582 173.505594 \n",
"L 169.940635 173.236188 \n",
"L 169.330688 172.966783 \n",
"L 168.720741 172.697378 \n",
"L 168.110794 172.427973 \n",
"L 167.500847 172.158568 \n",
"L 166.890899 171.889163 \n",
"L 166.280952 171.619758 \n",
"L 165.671005 171.350353 \n",
"L 165.061058 171.097681 \n",
"L 164.451111 170.856891 \n",
"L 163.841164 170.616102 \n",
"L 163.231216 170.375313 \n",
"L 162.621269 170.134524 \n",
"L 162.011322 169.893735 \n",
"L 161.401375 169.652946 \n",
"L 160.791428 169.412157 \n",
"L 160.181481 169.171368 \n",
"L 159.571533 168.930579 \n",
"L 158.961586 168.68979 \n",
"L 158.351639 168.449001 \n",
"L 157.741692 168.208212 \n",
"L 157.131745 167.967423 \n",
"L 156.521798 167.730047 \n",
"L 155.91185 167.532977 \n",
"L 155.301903 167.335908 \n",
"L 154.691956 167.138838 \n",
"L 154.082009 166.941769 \n",
"L 153.472062 166.744699 \n",
"L 152.862115 166.54763 \n",
"L 152.252167 166.35056 \n",
"L 151.64222 166.153491 \n",
"L 151.032273 165.956421 \n",
"L 150.422326 165.759352 \n",
"L 149.812379 165.562282 \n",
"L 149.202432 165.365213 \n",
"L 148.592484 165.168143 \n",
"L 147.982537 164.971074 \n",
"L 147.37259 164.807472 \n",
"L 146.762643 164.668972 \n",
"L 146.152696 164.530471 \n",
"L 145.542749 164.391971 \n",
"L 144.932801 164.25347 \n",
"L 144.322854 164.11497 \n",
"L 143.712907 163.97647 \n",
"L 143.10296 163.837969 \n",
"L 142.493013 163.699469 \n",
"L 141.883066 163.560968 \n",
"L 141.273118 163.422468 \n",
"L 140.663171 163.283968 \n",
"L 140.053224 163.145467 \n",
"L 139.443277 163.006967 \n",
"L 138.83333 162.873089 \n",
"L 138.223383 162.805945 \n",
"L 137.613435 162.738801 \n",
"L 137.003488 162.671657 \n",
"L 136.393541 162.604513 \n",
"L 135.783594 162.537369 \n",
"L 135.173647 162.470225 \n",
"L 134.5637 162.403081 \n",
"L 133.953752 162.335937 \n",
"L 133.343805 162.268793 \n",
"L 132.733858 162.201649 \n",
"L 132.123911 162.134505 \n",
"L 131.513964 162.067361 \n",
"L 130.904017 162.000217 \n",
"L 130.294069 161.933073 \n",
"L 129.684122 161.910669 \n",
"L 129.074175 161.923683 \n",
"L 128.464228 161.936698 \n",
"L 127.854281 161.949712 \n",
"L 127.244334 161.962727 \n",
"L 126.634386 161.975741 \n",
"L 126.024439 161.988756 \n",
"L 125.414492 162.00177 \n",
"L 124.804545 162.014785 \n",
"L 124.194598 162.027799 \n",
"L 123.584651 162.040814 \n",
"L 122.974703 162.053828 \n",
"L 122.364756 162.066843 \n",
"L 121.754809 162.079857 \n",
"L 121.144862 162.097183 \n",
"L 120.534915 162.193913 \n",
"L 119.924968 162.290642 \n",
"L 119.31502 162.387372 \n",
"L 118.705073 162.484102 \n",
"L 118.095126 162.580832 \n",
"L 117.485179 162.677561 \n",
"L 116.875232 162.774291 \n",
"L 116.265285 162.871021 \n",
"L 115.655337 162.967751 \n",
"L 115.04539 163.06448 \n",
"L 114.435443 163.16121 \n",
"L 113.825496 163.25794 \n",
"L 113.215549 163.35467 \n",
"L 112.605602 163.451399 \n",
"L 111.995654 163.591921 \n",
"L 111.385707 163.769024 \n",
"L 110.77576 163.946127 \n",
"L 110.165813 164.12323 \n",
"L 109.555866 164.300333 \n",
"L 108.945919 164.477436 \n",
"L 108.335971 164.654539 \n",
"L 107.726024 164.831642 \n",
"L 107.116077 165.008745 \n",
"L 106.50613 165.185848 \n",
"L 105.896183 165.36295 \n",
"L 105.286236 165.540053 \n",
"L 104.676288 165.717156 \n",
"L 104.066341 165.894259 \n",
"L 103.456394 166.074069 \n",
"L 102.846447 166.322017 \n",
"L 102.2365 166.569964 \n",
"L 101.626553 166.817911 \n",
"L 101.016605 167.065858 \n",
"L 100.406658 167.313805 \n",
"L 99.796711 167.561752 \n",
"L 99.186764 167.809699 \n",
"L 98.576817 168.057646 \n",
"L 97.96687 168.305593 \n",
"L 97.356922 168.55354 \n",
"L 96.746975 168.801487 \n",
"L 96.137028 169.049434 \n",
"L 95.527081 169.297381 \n",
"L 94.917134 169.545328 \n",
"L 94.307187 169.822885 \n",
"L 93.697239 170.126536 \n",
"L 93.087292 170.430186 \n",
"L 92.477345 170.733836 \n",
"L 91.867398 171.037486 \n",
"L 91.257451 171.341136 \n",
"L 90.647504 171.644787 \n",
"L 90.037556 171.948437 \n",
"L 89.427609 172.252087 \n",
"L 88.817662 172.555737 \n",
"L 88.207715 172.859387 \n",
"L 87.597768 173.163037 \n",
"L 86.987821 173.466688 \n",
"L 86.377873 173.770338 \n",
"L 85.767926 174.074907 \n",
"L 85.157979 174.415417 \n",
"L 84.548032 174.755927 \n",
"L 83.938085 175.096438 \n",
"L 83.328138 175.436948 \n",
"L 82.71819 175.777458 \n",
"L 82.108243 176.117969 \n",
"L 81.498296 176.458479 \n",
"L 80.888349 176.798989 \n",
"L 80.278402 177.1395 \n",
"L 79.668455 177.48001 \n",
"L 79.058507 177.82052 \n",
"L 78.44856 178.161031 \n",
"L 77.838613 178.501541 \n",
"L 77.228666 178.842051 \n",
"L 76.618719 179.191175 \n",
"L 76.008772 179.548303 \n",
"L 75.398824 179.905432 \n",
"L 74.788877 180.262561 \n",
"L 74.17893 180.619689 \n",
"L 73.568983 180.976818 \n",
"L 72.959036 181.333947 \n",
"L 72.349089 181.691076 \n",
"L 71.739141 182.048204 \n",
"L 71.129194 182.405333 \n",
"L 70.519247 182.762462 \n",
"L 69.9093 183.11959 \n",
"L 69.299353 183.476719 \n",
"L 68.689406 183.833848 \n",
"L 68.079458 184.190949 \n",
"L 67.469511 184.545679 \n",
"L 66.859564 184.900409 \n",
"L 66.249617 185.255139 \n",
"L 65.63967 185.609869 \n",
"L 65.029723 185.9646 \n",
"L 64.419775 186.31933 \n",
"L 63.809828 186.67406 \n",
"z\n",
"\" style=\"fill:#b3cde3;\"/>\n",
" </g>\n",
" <g id=\"PolyCollection_2\">\n",
" <path clip-path=\"url(#pd8980b1118)\" d=\"M 124.194598 104.659454 \n",
"L 124.194598 218.256364 \n",
"L 124.804545 218.256364 \n",
"L 125.414492 218.256364 \n",
"L 126.024439 218.256364 \n",
"L 126.634386 218.256364 \n",
"L 127.244334 218.256364 \n",
"L 127.854281 218.256364 \n",
"L 128.464228 218.256364 \n",
"L 129.074175 218.256364 \n",
"L 129.684122 218.256364 \n",
"L 130.294069 218.256364 \n",
"L 130.904017 218.256364 \n",
"L 131.513964 218.256364 \n",
"L 132.123911 218.256364 \n",
"L 132.733858 218.256364 \n",
"L 133.343805 218.256364 \n",
"L 133.953752 218.256364 \n",
"L 134.5637 218.256364 \n",
"L 135.173647 218.256364 \n",
"L 135.783594 218.256364 \n",
"L 136.393541 218.256364 \n",
"L 137.003488 218.256364 \n",
"L 137.613435 218.256364 \n",
"L 138.223383 218.256364 \n",
"L 138.83333 218.256364 \n",
"L 139.443277 218.256364 \n",
"L 140.053224 218.256364 \n",
"L 140.663171 218.256364 \n",
"L 141.273118 218.256364 \n",
"L 141.883066 218.256364 \n",
"L 142.493013 218.256364 \n",
"L 143.10296 218.256364 \n",
"L 143.712907 218.256364 \n",
"L 144.322854 218.256364 \n",
"L 144.932801 218.256364 \n",
"L 145.542749 218.256364 \n",
"L 146.152696 218.256364 \n",
"L 146.762643 218.256364 \n",
"L 147.37259 218.256364 \n",
"L 147.982537 218.256364 \n",
"L 148.592484 218.256364 \n",
"L 149.202432 218.256364 \n",
"L 149.812379 218.256364 \n",
"L 150.422326 218.256364 \n",
"L 151.032273 218.256364 \n",
"L 151.64222 218.256364 \n",
"L 152.252167 218.256364 \n",
"L 152.862115 218.256364 \n",
"L 153.472062 218.256364 \n",
"L 154.082009 218.256364 \n",
"L 154.691956 218.256364 \n",
"L 155.301903 218.256364 \n",
"L 155.91185 218.256364 \n",
"L 156.521798 218.256364 \n",
"L 157.131745 218.256364 \n",
"L 157.741692 218.256364 \n",
"L 158.351639 218.256364 \n",
"L 158.961586 218.256364 \n",
"L 159.571533 218.256364 \n",
"L 160.181481 218.256364 \n",
"L 160.791428 218.256364 \n",
"L 161.401375 218.256364 \n",
"L 162.011322 218.256364 \n",
"L 162.011322 105.442276 \n",
"L 162.011322 105.442276 \n",
"L 161.401375 102.015518 \n",
"L 160.791428 98.588759 \n",
"L 160.181481 95.139162 \n",
"L 159.571533 91.618912 \n",
"L 158.961586 88.098662 \n",
"L 158.351639 84.578413 \n",
"L 157.741692 81.061109 \n",
"L 157.131745 77.56771 \n",
"L 156.521798 74.074312 \n",
"L 155.91185 70.580914 \n",
"L 155.301903 67.087516 \n",
"L 154.691956 63.76663 \n",
"L 154.082009 60.450141 \n",
"L 153.472062 57.133653 \n",
"L 152.862115 53.817164 \n",
"L 152.252167 50.796787 \n",
"L 151.64222 47.832573 \n",
"L 151.032273 44.868359 \n",
"L 150.422326 41.904144 \n",
"L 149.812379 39.319903 \n",
"L 149.202432 36.893895 \n",
"L 148.592484 34.467887 \n",
"L 147.982537 32.041878 \n",
"L 147.37259 30.027876 \n",
"L 146.762643 28.322877 \n",
"L 146.152696 26.617879 \n",
"L 145.542749 24.912881 \n",
"L 144.932801 23.59163 \n",
"L 144.322854 22.765059 \n",
"L 143.712907 21.938488 \n",
"L 143.10296 21.111916 \n",
"L 142.493013 20.583636 \n",
"L 141.883066 20.74385 \n",
"L 141.273118 20.904064 \n",
"L 140.663171 21.064278 \n",
"L 140.053224 21.397333 \n",
"L 139.443277 22.588117 \n",
"L 138.83333 23.778901 \n",
"L 138.223383 24.969685 \n",
"L 137.613435 26.193261 \n",
"L 137.003488 28.373472 \n",
"L 136.393541 30.553684 \n",
"L 135.783594 32.733895 \n",
"L 135.173647 34.914106 \n",
"L 134.5637 37.877981 \n",
"L 133.953752 40.930313 \n",
"L 133.343805 43.982645 \n",
"L 132.733858 47.034977 \n",
"L 132.123911 50.611206 \n",
"L 131.513964 54.349267 \n",
"L 130.904017 58.087328 \n",
"L 130.294069 61.825389 \n",
"L 129.684122 65.849062 \n",
"L 129.074175 70.040886 \n",
"L 128.464228 74.232711 \n",
"L 127.854281 78.424536 \n",
"L 127.244334 82.717598 \n",
"L 126.634386 87.114002 \n",
"L 126.024439 91.510406 \n",
"L 125.414492 95.90681 \n",
"L 124.804545 100.292576 \n",
"L 124.194598 104.659454 \n",
"z\n",
"\" style=\"fill:#ccebc5;\"/>\n",
" </g>\n",
" <g id=\"matplotlib.axis_1\">\n",
" <g id=\"xtick_1\">\n",
" <g id=\"line2d_1\">\n",
" <defs>\n",
" <path d=\"M 0 0 \n",
"L 0 3.5 \n",
"\" id=\"me59544ec8f\" style=\"stroke:#000000;stroke-width:0.8;\"/>\n",
" </defs>\n",
" <g>\n",
" <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"60.760092\" xlink:href=\"#me59544ec8f\" y=\"228.14\"/>\n",
" </g>\n",
" </g>\n",
" <g id=\"text_1\">\n",
" <!-- $-0.004$ -->\n",
" <defs>\n",
" <path d=\"M 65.90625 23 \n",
"C 67.59375 23 69.40625 23 69.40625 25 \n",
"C 69.40625 27 67.59375 27 65.90625 27 \n",
"L 11.796875 27 \n",
"C 10.09375 27 8.296875 27 8.296875 25 \n",
"C 8.296875 23 10.09375 23 11.796875 23 \n",
"z\n",
"\" id=\"CMSY10-0\"/>\n",
" <path d=\"M 42 31.640625 \n",
"C 42 37.75 41.90625 48.125 37.703125 56.109375 \n",
"C 34 63.109375 28.09375 65.59375 22.90625 65.59375 \n",
"C 18.09375 65.59375 12 63.40625 8.203125 56.203125 \n",
"C 4.203125 48.71875 3.796875 39.4375 3.796875 31.640625 \n",
"C 3.796875 25.953125 3.90625 17.28125 7 9.671875 \n",
"C 11.296875 -0.609375 19 -2 22.90625 -2 \n",
"C 27.5 -2 34.5 -0.109375 38.59375 9.375 \n",
"C 41.59375 16.28125 42 24.359375 42 31.640625 \n",
"M 22.90625 -0.40625 \n",
"C 16.5 -0.40625 12.703125 5.078125 11.296875 12.6875 \n",
"C 10.203125 18.5625 10.203125 27.15625 10.203125 32.75 \n",
"C 10.203125 40.4375 10.203125 46.828125 11.5 52.921875 \n",
"C 13.40625 61.390625 19 64 22.90625 64 \n",
"C 27 64 32.296875 61.296875 34.203125 53.125 \n",
"C 35.5 47.4375 35.59375 40.734375 35.59375 32.75 \n",
"C 35.59375 26.25 35.59375 18.265625 34.40625 12.375 \n",
"C 32.296875 1.484375 26.40625 -0.40625 22.90625 -0.40625 \n",
"\" id=\"CMR17-48\"/>\n",
" <path d=\"M 18.40625 4.796875 \n",
"C 18.40625 7.6875 16 9.671875 13.59375 9.671875 \n",
"C 10.703125 9.671875 8.703125 7.28125 8.703125 4.890625 \n",
"C 8.703125 2 11.09375 0 13.5 0 \n",
"C 16.40625 0 18.40625 2.390625 18.40625 4.796875 \n",
"\" id=\"CMMI12-58\"/>\n",
" <path d=\"M 33.59375 64.40625 \n",
"C 33.59375 66.5 33.5 66.5 31.703125 66.5 \n",
"L 2 19.59375 \n",
"L 2 17 \n",
"L 27.796875 17 \n",
"L 27.796875 7.140625 \n",
"C 27.796875 3.5 27.59375 2.5 20.59375 2.5 \n",
"L 18.703125 2.5 \n",
"L 18.703125 0 \n",
"C 21.90625 0 27.296875 0 30.703125 0 \n",
"C 34.09375 0 39.5 0 42.703125 0 \n",
"L 42.703125 2.5 \n",
"L 40.796875 2.5 \n",
"C 33.796875 2.5 33.59375 3.5 33.59375 7.140625 \n",
"L 33.59375 17 \n",
"L 43.796875 17 \n",
"L 43.796875 19.59375 \n",
"L 33.59375 19.59375 \n",
"z\n",
"M 28.09375 57.859375 \n",
"L 28.09375 19.59375 \n",
"L 4 19.59375 \n",
"z\n",
"\" id=\"CMR17-52\"/>\n",
" </defs>\n",
" <g transform=\"translate(45.539346 242.058498)scale(0.1 -0.1)\">\n",
" <use transform=\"scale(0.996264)\" xlink:href=\"#CMSY10-0\"/>\n",
" <use transform=\"translate(77.487468 0)scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" <use transform=\"translate(123.177945 0)scale(0.996264)\" xlink:href=\"#CMMI12-58\"/>\n",
" <use transform=\"translate(150.275122 0)scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" <use transform=\"translate(195.965599 0)scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" <use transform=\"translate(241.656076 0)scale(0.996264)\" xlink:href=\"#CMR17-52\"/>\n",
" </g>\n",
" </g>\n",
" </g>\n",
" <g id=\"xtick_2\">\n",
" <g id=\"line2d_2\">\n",
" <g>\n",
" <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"104.240612\" xlink:href=\"#me59544ec8f\" y=\"228.14\"/>\n",
" </g>\n",
" </g>\n",
" <g id=\"text_2\">\n",
" <!-- $-0.002$ -->\n",
" <defs>\n",
" <path d=\"M 41.703125 15.453125 \n",
"L 39.90625 15.453125 \n",
"C 38.90625 8.375 38.09375 7.171875 37.703125 6.5625 \n",
"C 37.203125 5.765625 30 5.765625 28.59375 5.765625 \n",
"L 9.40625 5.765625 \n",
"C 13 9.671875 20 16.75 28.5 24.953125 \n",
"C 34.59375 30.734375 41.703125 37.53125 41.703125 47.421875 \n",
"C 41.703125 59.21875 32.296875 66 21.796875 66 \n",
"C 10.796875 66 4.09375 56.3125 4.09375 47.34375 \n",
"C 4.09375 43.4375 7 42.9375 8.203125 42.9375 \n",
"C 9.203125 42.9375 12.203125 43.546875 12.203125 47.03125 \n",
"C 12.203125 50.109375 9.59375 51 8.203125 51 \n",
"C 7.59375 51 7 50.90625 6.59375 50.703125 \n",
"C 8.5 59.21875 14.296875 63.40625 20.40625 63.40625 \n",
"C 29.09375 63.40625 34.796875 56.515625 34.796875 47.421875 \n",
"C 34.796875 38.734375 29.703125 31.25 24 24.75 \n",
"L 4.09375 2.28125 \n",
"L 4.09375 0 \n",
"L 39.296875 0 \n",
"z\n",
"\" id=\"CMR17-50\"/>\n",
" </defs>\n",
" <g transform=\"translate(89.019865 242.058498)scale(0.1 -0.1)\">\n",
" <use transform=\"scale(0.996264)\" xlink:href=\"#CMSY10-0\"/>\n",
" <use transform=\"translate(77.487468 0)scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" <use transform=\"translate(123.177945 0)scale(0.996264)\" xlink:href=\"#CMMI12-58\"/>\n",
" <use transform=\"translate(150.275122 0)scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" <use transform=\"translate(195.965599 0)scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" <use transform=\"translate(241.656076 0)scale(0.996264)\" xlink:href=\"#CMR17-50\"/>\n",
" </g>\n",
" </g>\n",
" </g>\n",
" <g id=\"xtick_3\">\n",
" <g id=\"line2d_3\">\n",
" <g>\n",
" <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"147.721131\" xlink:href=\"#me59544ec8f\" y=\"228.14\"/>\n",
" </g>\n",
" </g>\n",
" <g id=\"text_3\">\n",
" <!-- $0.000$ -->\n",
" <g transform=\"translate(136.374757 242.058498)scale(0.1 -0.1)\">\n",
" <use transform=\"scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" <use transform=\"translate(45.690477 0)scale(0.996264)\" xlink:href=\"#CMMI12-58\"/>\n",
" <use transform=\"translate(72.787654 0)scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" <use transform=\"translate(118.478131 0)scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" <use transform=\"translate(164.168608 0)scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" </g>\n",
" </g>\n",
" </g>\n",
" <g id=\"xtick_4\">\n",
" <g id=\"line2d_4\">\n",
" <g>\n",
" <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"191.201651\" xlink:href=\"#me59544ec8f\" y=\"228.14\"/>\n",
" </g>\n",
" </g>\n",
" <g id=\"text_4\">\n",
" <!-- $0.002$ -->\n",
" <g transform=\"translate(179.855276 242.058498)scale(0.1 -0.1)\">\n",
" <use transform=\"scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" <use transform=\"translate(45.690477 0)scale(0.996264)\" xlink:href=\"#CMMI12-58\"/>\n",
" <use transform=\"translate(72.787654 0)scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" <use transform=\"translate(118.478131 0)scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" <use transform=\"translate(164.168608 0)scale(0.996264)\" xlink:href=\"#CMR17-50\"/>\n",
" </g>\n",
" </g>\n",
" </g>\n",
" <g id=\"xtick_5\">\n",
" <g id=\"line2d_5\">\n",
" <g>\n",
" <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"234.68217\" xlink:href=\"#me59544ec8f\" y=\"228.14\"/>\n",
" </g>\n",
" </g>\n",
" <g id=\"text_5\">\n",
" <!-- $0.004$ -->\n",
" <g transform=\"translate(223.335796 242.058498)scale(0.1 -0.1)\">\n",
" <use transform=\"scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" <use transform=\"translate(45.690477 0)scale(0.996264)\" xlink:href=\"#CMMI12-58\"/>\n",
" <use transform=\"translate(72.787654 0)scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" <use transform=\"translate(118.478131 0)scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" <use transform=\"translate(164.168608 0)scale(0.996264)\" xlink:href=\"#CMR17-52\"/>\n",
" </g>\n",
" </g>\n",
" </g>\n",
" <g id=\"xtick_6\">\n",
" <g id=\"line2d_6\">\n",
" <g>\n",
" <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"278.16269\" xlink:href=\"#me59544ec8f\" y=\"228.14\"/>\n",
" </g>\n",
" </g>\n",
" <g id=\"text_6\">\n",
" <!-- $0.006$ -->\n",
" <defs>\n",
" <path d=\"M 10.59375 34 \n",
"C 10.59375 57.65625 21.796875 63 28.296875 63 \n",
"C 30.40625 63 35.5 62.640625 37.5 59 \n",
"C 35.90625 59 32.90625 59 32.90625 55.515625 \n",
"C 32.90625 52.828125 35.09375 51.921875 36.5 51.921875 \n",
"C 37.40625 51.921875 40.09375 52.3125 40.09375 55.625 \n",
"C 40.09375 61.78125 35.09375 65.296875 28.203125 65.296875 \n",
"C 16.296875 65.296875 3.796875 52.96875 3.796875 31 \n",
"C 3.796875 3.953125 15.09375 -2 23.09375 -2 \n",
"C 32.796875 -2 42 6.671875 42 20.046875 \n",
"C 42 32.515625 33.90625 41.59375 23.703125 41.59375 \n",
"C 17.59375 41.59375 13.09375 37.609375 10.59375 30.625 \n",
"z\n",
"M 23.09375 0.390625 \n",
"C 10.796875 0.390625 10.796875 18.75 10.796875 22.4375 \n",
"C 10.796875 29.625 14.203125 40 23.5 40 \n",
"C 25.203125 40 30.09375 40 33.40625 33.125 \n",
"C 35.203125 29.21875 35.203125 25.140625 35.203125 20.140625 \n",
"C 35.203125 14.75 35.203125 10.78125 33.09375 6.78125 \n",
"C 30.90625 2.671875 27.703125 0.390625 23.09375 0.390625 \n",
"\" id=\"CMR17-54\"/>\n",
" </defs>\n",
" <g transform=\"translate(266.816315 242.058498)scale(0.1 -0.1)\">\n",
" <use transform=\"scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" <use transform=\"translate(45.690477 0)scale(0.996264)\" xlink:href=\"#CMMI12-58\"/>\n",
" <use transform=\"translate(72.787654 0)scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" <use transform=\"translate(118.478131 0)scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" <use transform=\"translate(164.168608 0)scale(0.996264)\" xlink:href=\"#CMR17-54\"/>\n",
" </g>\n",
" </g>\n",
" </g>\n",
" <g id=\"xtick_7\">\n",
" <g id=\"line2d_7\">\n",
" <g>\n",
" <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"321.643209\" xlink:href=\"#me59544ec8f\" y=\"228.14\"/>\n",
" </g>\n",
" </g>\n",
" <g id=\"text_7\">\n",
" <!-- $0.008$ -->\n",
" <defs>\n",
" <path d=\"M 27.203125 35.375 \n",
"C 33.5 38.546875 39.90625 43.328125 39.90625 50.984375 \n",
"C 39.90625 60.015625 31.09375 65.296875 23 65.296875 \n",
"C 13.90625 65.296875 5.90625 58.734375 5.90625 49.6875 \n",
"C 5.90625 47.203125 6.5 42.921875 10.40625 39.140625 \n",
"C 11.40625 38.15625 15.59375 35.171875 18.296875 33.28125 \n",
"C 13.796875 30.984375 3.296875 25.53125 3.296875 14.59375 \n",
"C 3.296875 4.359375 13.09375 -2 22.796875 -2 \n",
"C 33.5 -2 42.5 5.65625 42.5 15.78125 \n",
"C 42.5 24.84375 36.40625 29.015625 32.40625 31.703125 \n",
"z\n",
"M 14.09375 44.09375 \n",
"C 13.296875 44.59375 9.296875 47.671875 9.296875 52.359375 \n",
"C 9.296875 58.421875 15.59375 63 22.796875 63 \n",
"C 30.703125 63 36.5 57.4375 36.5 50.96875 \n",
"C 36.5 41.703125 26.09375 36.421875 25.59375 36.421875 \n",
"C 25.5 36.421875 25.40625 36.421875 24.59375 37.03125 \n",
"z\n",
"M 32.5 23.734375 \n",
"C 34 22.640625 38.796875 19.375 38.796875 13.296875 \n",
"C 38.796875 5.953125 31.40625 0.390625 23 0.390625 \n",
"C 13.90625 0.390625 7 6.84375 7 14.6875 \n",
"C 7 22.546875 13.09375 29.09375 20 32.1875 \n",
"z\n",
"\" id=\"CMR17-56\"/>\n",
" </defs>\n",
" <g transform=\"translate(310.296835 242.058498)scale(0.1 -0.1)\">\n",
" <use transform=\"scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" <use transform=\"translate(45.690477 0)scale(0.996264)\" xlink:href=\"#CMMI12-58\"/>\n",
" <use transform=\"translate(72.787654 0)scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" <use transform=\"translate(118.478131 0)scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" <use transform=\"translate(164.168608 0)scale(0.996264)\" xlink:href=\"#CMR17-56\"/>\n",
" </g>\n",
" </g>\n",
" </g>\n",
" <g id=\"xtick_8\">\n",
" <g id=\"line2d_8\">\n",
" <g>\n",
" <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"365.123729\" xlink:href=\"#me59544ec8f\" y=\"228.14\"/>\n",
" </g>\n",
" </g>\n",
" <g id=\"text_8\">\n",
" <!-- $0.010$ -->\n",
" <defs>\n",
" <path d=\"M 26.59375 63.40625 \n",
"C 26.59375 65.5 26.5 65.5 25.09375 65.5 \n",
"C 21.203125 61.1875 15.296875 59.796875 9.703125 59.796875 \n",
"C 9.40625 59.796875 8.90625 59.796875 8.796875 59.5 \n",
"C 8.703125 59.296875 8.703125 59.09375 8.703125 57 \n",
"C 11.796875 57 17 57.59375 21 59.984375 \n",
"L 21 7.203125 \n",
"C 21 3.6875 20.796875 2.5 12.203125 2.5 \n",
"L 9.203125 2.5 \n",
"L 9.203125 0 \n",
"C 14 0 19 0 23.796875 0 \n",
"C 28.59375 0 33.59375 0 38.40625 0 \n",
"L 38.40625 2.5 \n",
"L 35.40625 2.5 \n",
"C 26.796875 2.5 26.59375 3.59375 26.59375 7.15625 \n",
"z\n",
"\" id=\"CMR17-49\"/>\n",
" </defs>\n",
" <g transform=\"translate(353.777354 242.058498)scale(0.1 -0.1)\">\n",
" <use transform=\"scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" <use transform=\"translate(45.690477 0)scale(0.996264)\" xlink:href=\"#CMMI12-58\"/>\n",
" <use transform=\"translate(72.787654 0)scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" <use transform=\"translate(118.478131 0)scale(0.996264)\" xlink:href=\"#CMR17-49\"/>\n",
" <use transform=\"translate(164.168608 0)scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" </g>\n",
" </g>\n",
" </g>\n",
" <g id=\"text_9\">\n",
" <!-- branching ratio -->\n",
" <defs>\n",
" <path d=\"M 15 68.6875 \n",
"L 3.09375 67.59375 \n",
"L 3.09375 65 \n",
"C 9 65 9.90625 64.40625 9.90625 59.640625 \n",
"L 9.90625 -0.203125 \n",
"L 11.703125 -0.203125 \n",
"L 14.59375 6.890625 \n",
"C 17.5 2.203125 21.796875 -1 27.59375 -1 \n",
"C 37.59375 -1 47.703125 7.8125 47.703125 21.40625 \n",
"C 47.703125 34.1875 38.703125 43.59375 28.59375 43.59375 \n",
"C 22.40625 43.59375 18 40.296875 15 36.40625 \n",
"z\n",
"M 15.203125 31.3125 \n",
"C 15.203125 33.09375 15.203125 33.40625 16.59375 35.5 \n",
"C 20.5 41.515625 26.09375 42 28 42 \n",
"C 31 42 41.296875 40.390625 41.296875 21.5 \n",
"C 41.296875 1.6875 29.5 0.59375 27.09375 0.59375 \n",
"C 24 0.59375 19.5 1.78125 16.296875 7.59375 \n",
"C 15.203125 9.5 15.203125 9.703125 15.203125 11.484375 \n",
"z\n",
"\" id=\"CMR17-98\"/>\n",
" <path d=\"M 15 23.390625 \n",
"C 15 33 19.09375 42 26.59375 42 \n",
"C 27.296875 42 28 41.90625 28.703125 41.59375 \n",
"C 28.703125 41.59375 26.5 40.890625 26.5 38.3125 \n",
"C 26.5 35.90625 28.40625 34.90625 29.90625 34.90625 \n",
"C 31.09375 34.90625 33.296875 35.609375 33.296875 38.40625 \n",
"C 33.296875 41.609375 30.09375 43.59375 26.703125 43.59375 \n",
"C 19.09375 43.59375 15.796875 36.203125 14.796875 32.703125 \n",
"L 14.703125 32.703125 \n",
"L 14.703125 43.59375 \n",
"L 3.09375 42.5 \n",
"L 3.09375 39.90625 \n",
"C 9 39.90625 9.90625 39.296875 9.90625 34.5 \n",
"L 9.90625 6.6875 \n",
"C 9.90625 2.875 9.5 2.390625 3.09375 2.390625 \n",
"L 3.09375 -0.203125 \n",
"C 5.5 0 10.09375 0 12.703125 0 \n",
"C 15.59375 0 20.703125 0 23.40625 -0.203125 \n",
"L 23.40625 2.390625 \n",
"C 16.203125 2.390625 15 2.390625 15 6.875 \n",
"z\n",
"\" id=\"CMR17-114\"/>\n",
" <path d=\"M 36 25.828125 \n",
"C 36 32.375 36 35.84375 31.796875 39.734375 \n",
"C 28.09375 43 23.796875 44 20.40625 44 \n",
"C 12.5 44 6.796875 37.609375 6.796875 30.796875 \n",
"C 6.796875 27 9.796875 27 10.40625 27 \n",
"C 11.703125 27 14 27.78125 14 30.53125 \n",
"C 14 33 12.09375 34.078125 10.40625 34.078125 \n",
"C 10 34.078125 9.5 33.96875 9.203125 33.875 \n",
"C 11.296875 40.5 16.703125 42.40625 20.203125 42.40625 \n",
"C 25.203125 42.40625 30.703125 37.953125 30.703125 29.453125 \n",
"L 30.703125 25 \n",
"C 24.796875 25 17.703125 24.203125 12.09375 21.203125 \n",
"C 5.796875 17.703125 4 12.703125 4 8.90625 \n",
"C 4 1.203125 13 -1 18.296875 -1 \n",
"C 23.796875 -1 28.90625 2.171875 31.09375 7.984375 \n",
"C 31.296875 3.640625 34.09375 -0.28125 38.5 -0.28125 \n",
"C 40.59375 -0.28125 45.90625 1.09375 45.90625 8.859375 \n",
"L 45.90625 14.375 \n",
"L 44.09375 14.375 \n",
"L 44.09375 8.78125 \n",
"C 44.09375 2.78125 41.40625 2 40.09375 2 \n",
"C 36 2 36 7.15625 36 11.53125 \n",
"z\n",
"M 30.703125 13.59375 \n",
"C 30.703125 4.953125 24.5 0.59375 19 0.59375 \n",
"C 14 0.59375 10.09375 4.296875 10.09375 8.90625 \n",
"C 10.09375 11.890625 11.40625 17.1875 17.203125 20.390625 \n",
"C 22 23.09375 27.5 23.5 30.703125 23.5 \n",
"z\n",
"\" id=\"CMR17-97\"/>\n",
" <path d=\"M 41.5 30.203125 \n",
"C 41.5 35.296875 40.5 43.59375 28.703125 43.59375 \n",
"C 20.796875 43.59375 16.703125 37.5 15.203125 33.5 \n",
"L 15.09375 33.5 \n",
"L 15.09375 43.59375 \n",
"L 3.296875 42.5 \n",
"L 3.296875 39.90625 \n",
"C 9.203125 39.90625 10.09375 39.296875 10.09375 34.5 \n",
"L 10.09375 6.6875 \n",
"C 10.09375 2.875 9.703125 2.390625 3.296875 2.390625 \n",
"L 3.296875 -0.203125 \n",
"C 5.703125 0 10.09375 0 12.703125 0 \n",
"C 15.296875 0 19.796875 0 22.203125 -0.203125 \n",
"L 22.203125 2.390625 \n",
"C 15.796875 2.390625 15.40625 2.78125 15.40625 6.6875 \n",
"L 15.40625 25.90625 \n",
"C 15.40625 35.109375 21 42 28 42 \n",
"C 35.40625 42 36.203125 35.40625 36.203125 30.59375 \n",
"L 36.203125 6.6875 \n",
"C 36.203125 2.875 35.796875 2.390625 29.40625 2.390625 \n",
"L 29.40625 -0.203125 \n",
"C 31.796875 0 36.203125 0 38.796875 0 \n",
"C 41.40625 0 45.90625 0 48.296875 -0.203125 \n",
"L 48.296875 2.390625 \n",
"C 41.90625 2.390625 41.5 2.78125 41.5 6.6875 \n",
"z\n",
"\" id=\"CMR17-110\"/>\n",
" <path d=\"M 34.90625 35 \n",
"C 33 35 30.203125 35 30.203125 31.515625 \n",
"C 30.203125 28.71875 32.5 27.921875 33.796875 27.921875 \n",
"C 34.5 27.921875 37.40625 28.234375 37.40625 31.625 \n",
"C 37.40625 38.546875 30.59375 43.796875 22.90625 43.796875 \n",
"C 12.296875 43.796875 3.296875 34.03125 3.296875 21.296875 \n",
"C 3.296875 8.0625 12.703125 -1 22.90625 -1 \n",
"C 35.296875 -1 38.40625 10.578125 38.40625 11.6875 \n",
"C 38.40625 12.09375 38.296875 12.390625 37.59375 12.390625 \n",
"C 36.90625 12.390625 36.796875 12.296875 36.40625 10.984375 \n",
"C 33.796875 2.8125 28.09375 0.796875 23.796875 0.796875 \n",
"C 17.40625 0.796875 9.703125 6.671875 9.703125 21.390625 \n",
"C 9.703125 36.53125 17.09375 42 23 42 \n",
"C 26.90625 42 32.703125 40.203125 34.90625 35 \n",
"\" id=\"CMR17-99\"/>\n",
" <path d=\"M 41.5 30.203125 \n",
"C 41.5 35.296875 40.5 43.59375 28.703125 43.59375 \n",
"C 20.5 43.59375 16.5 37 15.296875 33.703125 \n",
"L 15.203125 33.703125 \n",
"L 15.203125 68.6875 \n",
"L 3.296875 67.59375 \n",
"L 3.296875 65 \n",
"C 9.203125 65 10.09375 64.40625 10.09375 59.640625 \n",
"L 10.09375 6.65625 \n",
"C 10.09375 2.875 9.703125 2.390625 3.296875 2.390625 \n",
"L 3.296875 -0.203125 \n",
"C 5.703125 0 10.09375 0 12.703125 0 \n",
"C 15.296875 0 19.796875 0 22.203125 -0.203125 \n",
"L 22.203125 2.390625 \n",
"C 15.796875 2.390625 15.40625 2.78125 15.40625 6.6875 \n",
"L 15.40625 25.90625 \n",
"C 15.40625 35.109375 21 42 28 42 \n",
"C 35.40625 42 36.203125 35.40625 36.203125 30.59375 \n",
"L 36.203125 6.6875 \n",
"C 36.203125 2.875 35.796875 2.390625 29.40625 2.390625 \n",
"L 29.40625 -0.203125 \n",
"C 31.796875 0 36.203125 0 38.796875 0 \n",
"C 41.40625 0 45.90625 0 48.296875 -0.203125 \n",
"L 48.296875 2.390625 \n",
"C 41.90625 2.390625 41.5 2.78125 41.5 6.6875 \n",
"z\n",
"\" id=\"CMR17-104\"/>\n",
" <path d=\"M 15.5 61.078125 \n",
"C 15.5 63.265625 13.703125 65.15625 11.40625 65.15625 \n",
"C 9.203125 65.15625 7.296875 63.375 7.296875 61.078125 \n",
"C 7.296875 58.890625 9.09375 57 11.40625 57 \n",
"C 13.59375 57 15.5 58.796875 15.5 61.078125 \n",
"M 3.796875 42.28125 \n",
"L 3.796875 39.703125 \n",
"C 9.40625 39.703125 10.203125 39.09375 10.203125 34.328125 \n",
"L 10.203125 6.671875 \n",
"C 10.203125 2.875 9.796875 2.390625 3.40625 2.390625 \n",
"L 3.40625 -0.203125 \n",
"C 5.796875 0 10.09375 0 12.59375 0 \n",
"C 15 0 19.09375 0 21.40625 -0.203125 \n",
"L 21.40625 2.390625 \n",
"C 15.5 2.390625 15.296875 2.984375 15.296875 6.5625 \n",
"L 15.296875 43.375 \n",
"z\n",
"\" id=\"CMR17-105\"/>\n",
" <path d=\"M 11.09375 18.1875 \n",
"C 13 16.6875 16.296875 15 20.296875 15 \n",
"C 28.203125 15 35 21.203125 35 29.5 \n",
"C 35 32.09375 34.203125 36.09375 31.203125 39.203125 \n",
"C 34 42.0625 38.09375 43 40.40625 43 \n",
"C 40.796875 43 41.40625 43 41.90625 42.71875 \n",
"C 41.5 42.625 40.59375 42.234375 40.59375 40.703125 \n",
"C 40.59375 39.484375 41.5 38.609375 42.796875 38.609375 \n",
"C 44.296875 38.609375 45.09375 39.578125 45.09375 40.84375 \n",
"C 45.09375 42.671875 43.59375 44.59375 40.40625 44.59375 \n",
"C 36.40625 44.59375 32.90625 42.78125 30.09375 40.3125 \n",
"C 27.09375 43.015625 23.5 44 20.296875 44 \n",
"C 12.40625 44 5.59375 37.796875 5.59375 29.5 \n",
"C 5.59375 23.78125 8.90625 20.1875 9.90625 19.1875 \n",
"C 6.90625 15.6875 6.90625 11.609375 6.90625 11.109375 \n",
"C 6.90625 8.515625 7.90625 4.71875 11.296875 2.625 \n",
"C 6.09375 1.328125 2 -2.640625 2 -7.59375 \n",
"C 2 -14.734375 11.5 -20 22.90625 -20 \n",
"C 33.90625 -20 43.796875 -14.9375 43.796875 -7.484375 \n",
"C 43.796875 6 29 6 21.296875 6 \n",
"C 19 6 14.90625 6 14.40625 6.109375 \n",
"C 11.296875 6.609375 9.203125 9.40625 9.203125 12.90625 \n",
"C 9.203125 13.796875 9.203125 16 11.09375 18.1875 \n",
"M 20.296875 16.6875 \n",
"C 11.5 16.6875 11.5 27.5 11.5 29.515625 \n",
"C 11.5 31.5 11.5 42.3125 20.296875 42.3125 \n",
"C 29.09375 42.3125 29.09375 31.5 29.09375 29.515625 \n",
"C 29.09375 27.5 29.09375 16.6875 20.296875 16.6875 \n",
"M 22.90625 -18.296875 \n",
"C 12.90625 -18.296875 6 -13.046875 6 -7.59375 \n",
"C 6 -4.515625 7.703125 -1.640625 9.90625 -0.0625 \n",
"C 12.40625 1.625 13.40625 1.625 20.203125 1.625 \n",
"C 28.40625 1.625 39.796875 1.625 39.796875 -7.59375 \n",
"C 39.796875 -13.046875 32.90625 -18.296875 22.90625 -18.296875 \n",
"\" id=\"CMR17-103\"/>\n",
" <path d=\"M 15.09375 40 \n",
"L 29.09375 40 \n",
"L 29.09375 42.578125 \n",
"L 15.09375 42.578125 \n",
"L 15.09375 61 \n",
"L 13.296875 61 \n",
"C 13.09375 50.703125 9.59375 41.6875 1.09375 41.6875 \n",
"L 1.09375 40 \n",
"L 9.796875 40 \n",
"L 9.796875 12.046875 \n",
"C 9.796875 10.15625 9.796875 -1 21.40625 -1 \n",
"C 27.296875 -1 30.703125 4.78125 30.703125 12.140625 \n",
"L 30.703125 17.8125 \n",
"L 28.90625 17.8125 \n",
"L 28.90625 12.234375 \n",
"C 28.90625 5.375 26.203125 0.796875 22 0.796875 \n",
"C 19.09375 0.796875 15.09375 2.78125 15.09375 11.84375 \n",
"z\n",
"\" id=\"CMR17-116\"/>\n",
" <path d=\"M 43.09375 21.1875 \n",
"C 43.09375 34 33.796875 44 22.90625 44 \n",
"C 12 44 2.703125 34 2.703125 21.1875 \n",
"C 2.703125 8.609375 12 -1 22.90625 -1 \n",
"C 33.796875 -1 43.09375 8.609375 43.09375 21.1875 \n",
"M 22.90625 0.796875 \n",
"C 18.203125 0.796875 14.203125 3.59375 11.90625 7.5 \n",
"C 9.40625 11.984375 9.09375 17.59375 9.09375 22 \n",
"C 9.09375 26.203125 9.296875 31.40625 11.90625 35.90625 \n",
"C 13.90625 39.203125 17.796875 42.40625 22.90625 42.40625 \n",
"C 27.40625 42.40625 31.203125 39.90625 33.59375 36.40625 \n",
"C 36.703125 31.703125 36.703125 25.109375 36.703125 22 \n",
"C 36.703125 18.09375 36.5 12.109375 33.796875 7.296875 \n",
"C 31 2.703125 26.703125 0.796875 22.90625 0.796875 \n",
"\" id=\"CMR17-111\"/>\n",
" </defs>\n",
" <g transform=\"translate(167.228973 257.681553)scale(0.14 -0.14)\">\n",
" <use transform=\"scale(0.996264)\" xlink:href=\"#CMR17-98\"/>\n",
" <use transform=\"translate(50.89537 0)scale(0.996264)\" xlink:href=\"#CMR17-114\"/>\n",
" <use transform=\"translate(86.176138 0)scale(0.996264)\" xlink:href=\"#CMR17-97\"/>\n",
" <use transform=\"translate(131.866615 0)scale(0.996264)\" xlink:href=\"#CMR17-110\"/>\n",
" <use transform=\"translate(182.761984 0)scale(0.996264)\" xlink:href=\"#CMR17-99\"/>\n",
" <use transform=\"translate(220.645191 0)scale(0.996264)\" xlink:href=\"#CMR17-104\"/>\n",
" <use transform=\"translate(271.540561 0)scale(0.996264)\" xlink:href=\"#CMR17-105\"/>\n",
" <use transform=\"translate(296.411619 0)scale(0.996264)\" xlink:href=\"#CMR17-110\"/>\n",
" <use transform=\"translate(347.306989 0)scale(0.996264)\" xlink:href=\"#CMR17-103\"/>\n",
" <use transform=\"translate(423.073433 0)scale(0.996264)\" xlink:href=\"#CMR17-114\"/>\n",
" <use transform=\"translate(458.3542 0)scale(0.996264)\" xlink:href=\"#CMR17-97\"/>\n",
" <use transform=\"translate(504.044677 0)scale(0.996264)\" xlink:href=\"#CMR17-116\"/>\n",
" <use transform=\"translate(539.325445 0)scale(0.996264)\" xlink:href=\"#CMR17-105\"/>\n",
" <use transform=\"translate(564.196504 0)scale(0.996264)\" xlink:href=\"#CMR17-111\"/>\n",
" </g>\n",
" </g>\n",
" </g>\n",
" <g id=\"matplotlib.axis_2\">\n",
" <g id=\"ytick_1\">\n",
" <g id=\"line2d_9\">\n",
" <defs>\n",
" <path d=\"M 0 0 \n",
"L -3.5 0 \n",
"\" id=\"mb2f2f96498\" style=\"stroke:#000000;stroke-width:0.8;\"/>\n",
" </defs>\n",
" <g>\n",
" <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"45.541911\" xlink:href=\"#mb2f2f96498\" y=\"218.256364\"/>\n",
" </g>\n",
" </g>\n",
" <g id=\"text_10\">\n",
" <!-- $0$ -->\n",
" <g transform=\"translate(33.560575 221.715613)scale(0.1 -0.1)\">\n",
" <use transform=\"scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" </g>\n",
" </g>\n",
" </g>\n",
" <g id=\"ytick_2\">\n",
" <g id=\"line2d_10\">\n",
" <g>\n",
" <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"45.541911\" xlink:href=\"#mb2f2f96498\" y=\"175.381851\"/>\n",
" </g>\n",
" </g>\n",
" <g id=\"text_11\">\n",
" <!-- $100$ -->\n",
" <g transform=\"translate(23.597905 178.8411)scale(0.1 -0.1)\">\n",
" <use transform=\"scale(0.996264)\" xlink:href=\"#CMR17-49\"/>\n",
" <use transform=\"translate(45.690477 0)scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" <use transform=\"translate(91.380954 0)scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" </g>\n",
" </g>\n",
" </g>\n",
" <g id=\"ytick_3\">\n",
" <g id=\"line2d_11\">\n",
" <g>\n",
" <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"45.541911\" xlink:href=\"#mb2f2f96498\" y=\"132.507339\"/>\n",
" </g>\n",
" </g>\n",
" <g id=\"text_12\">\n",
" <!-- $200$ -->\n",
" <g transform=\"translate(23.597905 135.966588)scale(0.1 -0.1)\">\n",
" <use transform=\"scale(0.996264)\" xlink:href=\"#CMR17-50\"/>\n",
" <use transform=\"translate(45.690477 0)scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" <use transform=\"translate(91.380954 0)scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" </g>\n",
" </g>\n",
" </g>\n",
" <g id=\"ytick_4\">\n",
" <g id=\"line2d_12\">\n",
" <g>\n",
" <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"45.541911\" xlink:href=\"#mb2f2f96498\" y=\"89.632826\"/>\n",
" </g>\n",
" </g>\n",
" <g id=\"text_13\">\n",
" <!-- $300$ -->\n",
" <defs>\n",
" <path d=\"M 22.09375 33.703125 \n",
"C 31 33.703125 34.90625 25.96875 34.90625 17.03125 \n",
"C 34.90625 5 28.5 0.390625 22.703125 0.390625 \n",
"C 17.40625 0.390625 8.796875 3.03125 6.09375 10.828125 \n",
"C 6.59375 10.625 7.09375 10.625 7.59375 10.625 \n",
"C 10 10.625 11.796875 12.21875 11.796875 14.8125 \n",
"C 11.796875 17.703125 9.59375 19 7.59375 19 \n",
"C 5.90625 19 3.296875 18.203125 3.296875 14.46875 \n",
"C 3.296875 5.234375 12.296875 -2 22.90625 -2 \n",
"C 34 -2 42.5 6.71875 42.5 16.953125 \n",
"C 42.5 26.671875 34.5 33.703125 25 34.796875 \n",
"C 32.59375 36.375 39.90625 43.0625 39.90625 52.015625 \n",
"C 39.90625 59.6875 32 65.296875 23 65.296875 \n",
"C 13.90625 65.296875 5.90625 59.828125 5.90625 52.015625 \n",
"C 5.90625 48.59375 8.5 48 9.796875 48 \n",
"C 11.90625 48 13.703125 49.296875 13.703125 51.890625 \n",
"C 13.703125 54.46875 11.90625 55.765625 9.796875 55.765625 \n",
"C 9.40625 55.765625 8.90625 55.765625 8.5 55.578125 \n",
"C 11.40625 61.859375 19.296875 63 22.796875 63 \n",
"C 26.296875 63 32.90625 61.328125 32.90625 51.875 \n",
"C 32.90625 49.109375 32.5 44.1875 29.09375 39.84375 \n",
"C 26.09375 36 22.703125 36 19.40625 35.6875 \n",
"C 18.90625 35.6875 16.59375 35.453125 16.203125 35.453125 \n",
"C 15.5 35.359375 15.09375 35.265625 15.09375 34.5 \n",
"C 15.09375 33.796875 15.203125 33.703125 17.203125 33.703125 \n",
"z\n",
"\" id=\"CMR17-51\"/>\n",
" </defs>\n",
" <g transform=\"translate(23.597905 93.092075)scale(0.1 -0.1)\">\n",
" <use transform=\"scale(0.996264)\" xlink:href=\"#CMR17-51\"/>\n",
" <use transform=\"translate(45.690477 0)scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" <use transform=\"translate(91.380954 0)scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" </g>\n",
" </g>\n",
" </g>\n",
" <g id=\"ytick_5\">\n",
" <g id=\"line2d_13\">\n",
" <g>\n",
" <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"45.541911\" xlink:href=\"#mb2f2f96498\" y=\"46.758314\"/>\n",
" </g>\n",
" </g>\n",
" <g id=\"text_14\">\n",
" <!-- $400$ -->\n",
" <g transform=\"translate(23.597905 50.217563)scale(0.1 -0.1)\">\n",
" <use transform=\"scale(0.996264)\" xlink:href=\"#CMR17-52\"/>\n",
" <use transform=\"translate(45.690477 0)scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" <use transform=\"translate(91.380954 0)scale(0.996264)\" xlink:href=\"#CMR17-48\"/>\n",
" </g>\n",
" </g>\n",
" </g>\n",
" <g id=\"text_15\">\n",
" <!-- likelihood -->\n",
" <defs>\n",
" <path d=\"M 15.296875 68.6875 \n",
"L 3.40625 67.59375 \n",
"L 3.40625 65 \n",
"C 9.296875 65 10.203125 64.40625 10.203125 59.640625 \n",
"L 10.203125 6.65625 \n",
"C 10.203125 2.875 9.796875 2.390625 3.40625 2.390625 \n",
"L 3.40625 -0.203125 \n",
"C 5.796875 0 10.203125 0 12.703125 0 \n",
"C 15.296875 0 19.703125 0 22.09375 -0.203125 \n",
"L 22.09375 2.390625 \n",
"C 15.703125 2.390625 15.296875 2.78125 15.296875 6.65625 \n",
"z\n",
"\" id=\"CMR17-108\"/>\n",
" <path d=\"M 26.09375 26.140625 \n",
"C 26 26.34375 25.5 26.9375 25.5 27.234375 \n",
"C 25.5 27.4375 32.59375 33.921875 33.5 34.71875 \n",
"C 39.203125 39.5 42.796875 40 45.40625 40 \n",
"L 45.40625 42.703125 \n",
"C 44.796875 42.703125 44.59375 42.703125 44.5 42.5 \n",
"L 39.59375 42.5 \n",
"C 36.796875 42.5 32.5 42.5 29.796875 42.703125 \n",
"L 29.796875 40.09375 \n",
"C 31.203125 39.984375 32.203125 39.1875 32.203125 37.796875 \n",
"C 32.203125 35.796875 30.203125 33.90625 30.09375 33.90625 \n",
"L 15 20.140625 \n",
"L 15 68.6875 \n",
"L 3.09375 67.59375 \n",
"L 3.09375 65 \n",
"C 9 65 9.90625 64.40625 9.90625 59.640625 \n",
"L 9.90625 6.65625 \n",
"C 9.90625 2.875 9.5 2.390625 3.09375 2.390625 \n",
"L 3.09375 -0.203125 \n",
"C 5.5 0 9.90625 0 12.40625 0 \n",
"C 14.90625 0 19.296875 0 21.703125 -0.203125 \n",
"L 21.703125 2.390625 \n",
"C 15.296875 2.390625 14.90625 2.78125 14.90625 6.6875 \n",
"L 14.90625 17.65625 \n",
"L 21.796875 23.9375 \n",
"L 31.90625 8.96875 \n",
"C 33.09375 7.171875 34.09375 5.765625 34.09375 4.46875 \n",
"C 34.09375 2.6875 32.296875 2.390625 31.09375 2.390625 \n",
"L 31.09375 -0.203125 \n",
"C 33.5 0 37.796875 0 40.296875 0 \n",
"C 42.09375 0 45.796875 0 47.5 0 \n",
"L 47.5 2.390625 \n",
"C 42.796875 2.390625 41.703125 3.28125 38.703125 7.671875 \n",
"z\n",
"\" id=\"CMR17-107\"/>\n",
" <path d=\"M 38.09375 23 \n",
"C 38.5 23.40625 38.5 23.609375 38.5 24.625 \n",
"C 38.5 34.96875 33.09375 44 21.703125 44 \n",
"C 11.09375 44 2.703125 33.890625 2.703125 21.609375 \n",
"C 2.703125 8.609375 12.203125 -1 22.796875 -1 \n",
"C 34 -1 38.40625 9.484375 38.40625 11.5625 \n",
"C 38.40625 12.25 37.796875 12.25 37.59375 12.25 \n",
"C 36.90625 12.25 36.796875 12.046875 36.40625 10.875 \n",
"C 34.203125 4.140625 28.703125 0.796875 23.5 0.796875 \n",
"C 19.203125 0.796875 14.90625 3.171875 12.203125 7.5 \n",
"C 9.09375 12.546875 9.09375 18.375 9.09375 23 \n",
"z\n",
"M 9.203125 24.5 \n",
"C 9.90625 39.140625 17.59375 42.40625 21.59375 42.40625 \n",
"C 28.40625 42.40625 33 35.890625 33.09375 24.5 \n",
"z\n",
"\" id=\"CMR17-101\"/>\n",
" <path d=\"M 29.203125 67.59375 \n",
"L 29.203125 65 \n",
"C 35.09375 65 36 64.40625 36 59.6875 \n",
"L 36 36.640625 \n",
"C 35.59375 37.140625 31.5 43.59375 23.40625 43.59375 \n",
"C 13.203125 43.59375 3.296875 34.5 3.296875 21.3125 \n",
"C 3.296875 8.203125 12.59375 -1 22.40625 -1 \n",
"C 30.90625 -1 35.296875 5.265625 35.796875 5.9375 \n",
"L 35.796875 -0.984375 \n",
"L 47.90625 -0.984375 \n",
"L 47.90625 2 \n",
"C 42 2 41.09375 2.59375 41.09375 7.390625 \n",
"L 41.09375 68.6875 \n",
"z\n",
"M 35.796875 11.796875 \n",
"C 35.796875 8.796875 34 6.09375 31.703125 4.09375 \n",
"C 28.296875 1.09375 24.90625 0.59375 23 0.59375 \n",
"C 20.09375 0.59375 9.703125 2.09375 9.703125 21.1875 \n",
"C 9.703125 40.8125 21.296875 42 23.90625 42 \n",
"C 28.5 42 32.203125 39.40625 34.5 35.796875 \n",
"C 35.796875 33.703125 35.796875 33.40625 35.796875 31.609375 \n",
"z\n",
"\" id=\"CMR17-100\"/>\n",
" </defs>\n",
" <g transform=\"translate(16.885876 148.251396)rotate(-90)scale(0.14 -0.14)\">\n",
" <use transform=\"scale(0.996264)\" xlink:href=\"#CMR17-108\"/>\n",
" <use transform=\"translate(24.871058 0)scale(0.996264)\" xlink:href=\"#CMR17-105\"/>\n",
" <use transform=\"translate(49.742117 0)scale(0.996264)\" xlink:href=\"#CMR17-107\"/>\n",
" <use transform=\"translate(95.432579 0)scale(0.996264)\" xlink:href=\"#CMR17-101\"/>\n",
" <use transform=\"translate(135.918239 0)scale(0.996264)\" xlink:href=\"#CMR17-108\"/>\n",
" <use transform=\"translate(160.789298 0)scale(0.996264)\" xlink:href=\"#CMR17-105\"/>\n",
" <use transform=\"translate(185.660356 0)scale(0.996264)\" xlink:href=\"#CMR17-104\"/>\n",
" <use transform=\"translate(236.555726 0)scale(0.996264)\" xlink:href=\"#CMR17-111\"/>\n",
" <use transform=\"translate(284.848642 0)scale(0.996264)\" xlink:href=\"#CMR17-111\"/>\n",
" <use transform=\"translate(333.141558 0)scale(0.996264)\" xlink:href=\"#CMR17-100\"/>\n",
" </g>\n",
" </g>\n",
" </g>\n",
" <g id=\"line2d_14\">\n",
" <path clip-path=\"url(#pd8980b1118)\" d=\"M 60.760092 188.447711 \n",
"L 83.938085 175.096438 \n",
"L 87.597768 173.163037 \n",
"L 95.527081 169.297381 \n",
"L 104.066341 165.894259 \n",
"L 112.605602 163.451399 \n",
"L 121.754809 162.079857 \n",
"L 130.294069 161.933073 \n",
"L 139.443277 163.006967 \n",
"L 148.592484 165.168143 \n",
"L 157.741692 168.208212 \n",
"L 167.500847 172.158568 \n",
"L 179.089843 177.390917 \n",
"L 202.877782 188.34275 \n",
"L 213.246884 192.707251 \n",
"L 221.176197 195.779756 \n",
"L 230.325405 198.988831 \n",
"L 238.864665 201.66061 \n",
"L 248.013873 204.174955 \n",
"L 257.16308 206.353689 \n",
"L 266.312288 208.224493 \n",
"L 275.461495 209.819096 \n",
"L 284.610703 211.169885 \n",
"L 294.369857 212.377971 \n",
"L 304.129012 213.382212 \n",
"L 313.888167 214.21417 \n",
"L 324.257269 214.943197 \n",
"L 335.846265 215.587145 \n",
"L 354.14468 216.353302 \n",
"L 365.123729 216.69646 \n",
"L 365.123729 216.69646 \n",
"\" style=\"fill:none;stroke:#377eb8;stroke-linecap:square;stroke-width:0.6;\"/>\n",
" </g>\n",
" <g id=\"line2d_15\">\n",
" <path clip-path=\"url(#pd8980b1118)\" d=\"M 60.760092 218.256364 \n",
"L 84.548032 218.256364 \n",
"L 85.157979 216.37999 \n",
"L 88.207715 215.511071 \n",
"L 90.647504 214.552261 \n",
"L 92.477345 213.655101 \n",
"L 94.917134 212.078457 \n",
"L 97.356922 209.990024 \n",
"L 99.796711 207.241598 \n",
"L 102.2365 203.62761 \n",
"L 104.676288 198.965301 \n",
"L 106.50613 194.538606 \n",
"L 107.726024 191.297355 \n",
"L 110.165813 183.506683 \n",
"L 112.605602 173.984479 \n",
"L 115.04539 162.595271 \n",
"L 118.095126 145.54522 \n",
"L 120.534915 130.126817 \n",
"L 124.194598 104.659454 \n",
"L 129.684122 65.849062 \n",
"L 131.513964 54.349267 \n",
"L 132.733858 47.034977 \n",
"L 135.173647 34.914106 \n",
"L 137.613435 26.193261 \n",
"L 140.053224 21.397333 \n",
"L 140.663171 21.064278 \n",
"L 142.493013 20.583636 \n",
"L 143.10296 21.111916 \n",
"L 144.932801 23.59163 \n",
"L 145.542749 24.912881 \n",
"L 147.982537 32.041878 \n",
"L 150.422326 41.904144 \n",
"L 153.472062 57.133653 \n",
"L 157.131745 77.56771 \n",
"L 163.841164 115.42697 \n",
"L 166.280952 128.044179 \n",
"L 168.720741 139.68071 \n",
"L 171.16053 150.238759 \n",
"L 173.600318 159.682956 \n",
"L 176.040107 168.031779 \n",
"L 178.479896 175.339456 \n",
"L 180.919684 181.68291 \n",
"L 183.359473 187.151648 \n",
"L 185.799262 191.84025 \n",
"L 188.848997 196.693955 \n",
"L 191.288786 199.960507 \n",
"L 193.728575 202.729049 \n",
"L 196.168363 205.070983 \n",
"L 199.218099 207.477367 \n",
"L 201.657888 209.078786 \n",
"L 204.097677 210.430679 \n",
"L 207.147412 211.825438 \n",
"L 209.587201 212.749627 \n",
"L 210.197148 212.962515 \n",
"L 210.807095 218.256364 \n",
"L 365.123729 218.256364 \n",
"L 365.123729 218.256364 \n",
"\" style=\"fill:none;stroke:#4daf4a;stroke-linecap:square;stroke-width:0.6;\"/>\n",
" </g>\n",
" <g id=\"patch_3\">\n",
" <path d=\"M 45.541911 228.14 \n",
"L 45.541911 10.7 \n",
"\" style=\"fill:none;stroke:#000000;stroke-linecap:square;stroke-linejoin:miter;stroke-width:0.8;\"/>\n",
" </g>\n",
" <g id=\"patch_4\">\n",
" <path d=\"M 380.341911 228.14 \n",
"L 380.341911 10.7 \n",
"\" style=\"fill:none;stroke:#000000;stroke-linecap:square;stroke-linejoin:miter;stroke-width:0.8;\"/>\n",
" </g>\n",
" <g id=\"patch_5\">\n",
" <path d=\"M 45.541911 228.14 \n",
"L 380.341911 228.14 \n",
"\" style=\"fill:none;stroke:#000000;stroke-linecap:square;stroke-linejoin:miter;stroke-width:0.8;\"/>\n",
" </g>\n",
" <g id=\"patch_6\">\n",
" <path d=\"M 45.541911 10.7 \n",
"L 380.341911 10.7 \n",
"\" style=\"fill:none;stroke:#000000;stroke-linecap:square;stroke-linejoin:miter;stroke-width:0.8;\"/>\n",
" </g>\n",
" <g id=\"legend_1\">\n",
" <g id=\"patch_7\">\n",
" <path d=\"M 303.480503 46.574357 \n",
"L 373.341911 46.574357 \n",
"Q 375.341911 46.574357 375.341911 44.574357 \n",
"L 375.341911 17.7 \n",
"Q 375.341911 15.7 373.341911 15.7 \n",
"L 303.480503 15.7 \n",
"Q 301.480503 15.7 301.480503 17.7 \n",
"L 301.480503 44.574357 \n",
"Q 301.480503 46.574357 303.480503 46.574357 \n",
"z\n",
"\" style=\"fill:#ffffff;opacity:0.8;stroke:#cccccc;stroke-linejoin:miter;\"/>\n",
" </g>\n",
" <g id=\"line2d_16\">\n",
" <path d=\"M 305.480503 23.2 \n",
"L 325.480503 23.2 \n",
"\" style=\"fill:none;stroke:#377eb8;stroke-linecap:square;stroke-width:0.6;\"/>\n",
" </g>\n",
" <g id=\"line2d_17\"/>\n",
" <g id=\"text_16\">\n",
" <!-- $B_s\\to\\tau\\tau$ -->\n",
" <defs>\n",
" <path d=\"M 36.59375 61.265625 \n",
"C 37.5 65 37.90625 65 41.796875 65 \n",
"L 54.796875 65 \n",
"C 66.09375 65 66.09375 55.375 66.09375 54.453125 \n",
"C 66.09375 46.34375 58 36 44.796875 36 \n",
"L 30.40625 36 \n",
"z\n",
"M 53.5 35.546875 \n",
"C 64.40625 37.53125 74.296875 45.09375 74.296875 54.265625 \n",
"C 74.296875 62.03125 67.40625 68 56.09375 68 \n",
"L 24 68 \n",
"C 22.09375 68 21.203125 68 21.203125 66.03125 \n",
"C 21.203125 65 22.09375 65 23.59375 65 \n",
"C 29.703125 65 29.703125 64.25 29.703125 63.1875 \n",
"C 29.703125 63 29.703125 62.390625 29.296875 60.890625 \n",
"L 15.796875 7.375 \n",
"C 14.90625 3.890625 14.703125 2.890625 7.703125 2.890625 \n",
"C 5.796875 2.890625 4.796875 2.890625 4.796875 1.09375 \n",
"C 4.796875 0 5.40625 0 7.40625 0 \n",
"L 41.703125 0 \n",
"C 57 0 68.796875 11.546875 68.796875 21.609375 \n",
"C 68.796875 29.78125 61.59375 34.75 53.5 35.546875 \n",
"M 39.296875 2.890625 \n",
"L 25.796875 2.890625 \n",
"C 24.40625 2.890625 24.203125 2.890625 23.59375 3 \n",
"C 22.5 3.09375 22.40625 3.296875 22.40625 4.078125 \n",
"C 22.40625 4.765625 22.59375 5.359375 22.796875 6.234375 \n",
"L 29.796875 34 \n",
"L 48.59375 34 \n",
"C 60.40625 34 60.40625 23.171875 60.40625 22.390625 \n",
"C 60.40625 12.9375 51.703125 2.890625 39.296875 2.890625 \n",
"\" id=\"CMMI12-66\"/>\n",
" <path d=\"M 22.796875 20 \n",
"C 24.5 19.703125 27.203125 19.109375 27.796875 19 \n",
"C 29.09375 18.625 33.59375 17.03125 33.59375 12.234375 \n",
"C 33.59375 9.15625 30.796875 1.09375 19.203125 1.09375 \n",
"C 17.09375 1.09375 9.59375 1.09375 7.59375 6.8125 \n",
"C 11.59375 6.3125 13.59375 9.421875 13.59375 11.640625 \n",
"C 13.59375 13.75 12.203125 14.84375 10.203125 14.84375 \n",
"C 8 14.84375 5.09375 13.140625 5.09375 8.625 \n",
"C 5.09375 2.703125 11.09375 -1 19.09375 -1 \n",
"C 34.296875 -1 38.796875 10.171875 38.796875 15.375 \n",
"C 38.796875 16.859375 38.796875 19.65625 35.59375 22.84375 \n",
"C 33.09375 25.25 30.703125 25.75 25.296875 26.84375 \n",
"C 22.59375 27.4375 18.296875 28.34375 18.296875 32.828125 \n",
"C 18.296875 34.828125 20.09375 42.015625 29.59375 42.015625 \n",
"C 33.796875 42.015625 37.90625 40.453125 38.90625 37 \n",
"C 34.5 37 34.296875 33.203125 34.296875 33.109375 \n",
"C 34.296875 31 36.203125 30.40625 37.09375 30.40625 \n",
"C 38.5 30.40625 41.296875 31.484375 41.296875 35.640625 \n",
"C 41.296875 39.765625 37.5 44 29.703125 44 \n",
"C 16.59375 44 13.09375 33.75 13.09375 29.671875 \n",
"C 13.09375 22.09375 20.5 20.5 22.796875 20 \n",
"\" id=\"CMMI12-115\"/>\n",
" <path d=\"M 83.40625 23 \n",
"C 77.90625 18.796875 75.203125 14.703125 74.40625 13.40625 \n",
"C 69.90625 6.5 69.09375 0.203125 69.09375 0.09375 \n",
"C 69.09375 -1.09375 70.296875 -1.09375 71.09375 -1.09375 \n",
"C 72.796875 -1.09375 72.90625 -0.90625 73.296875 0.90625 \n",
"C 75.59375 10.703125 81.5 19.09375 92.796875 23.703125 \n",
"C 94 24.09375 94.296875 24.296875 94.296875 25 \n",
"C 94.296875 25.703125 93.703125 26 93.5 26.09375 \n",
"C 89.09375 27.796875 77 32.796875 73.203125 49.59375 \n",
"C 72.90625 50.796875 72.796875 51.09375 71.09375 51.09375 \n",
"C 70.296875 51.09375 69.09375 51.09375 69.09375 49.890625 \n",
"C 69.09375 49.6875 70 43.390625 74.203125 36.703125 \n",
"C 76.203125 33.703125 79.09375 30.203125 83.40625 27 \n",
"L 9.09375 27 \n",
"C 7.296875 27 5.5 27 5.5 25 \n",
"C 5.5 23 7.296875 23 9.09375 23 \n",
"z\n",
"\" id=\"CMSY10-33\"/>\n",
" <path d=\"M 28.703125 37.625 \n",
"L 45.296875 37.625 \n",
"C 46.59375 37.625 49.90625 37.625 49.90625 40.8125 \n",
"C 49.90625 43 48 43 46.203125 43 \n",
"L 18.703125 43 \n",
"C 16.40625 43 13 43 8.40625 38.125 \n",
"C 5.796875 35.234375 2.59375 29.953125 2.59375 29.359375 \n",
"C 2.59375 28.765625 3.09375 28.5625 3.703125 28.5625 \n",
"C 4.40625 28.5625 4.5 28.859375 5 29.453125 \n",
"C 10.203125 37.625 15.40625 37.625 17.90625 37.625 \n",
"L 26.203125 37.625 \n",
"L 15.796875 3.484375 \n",
"C 15.296875 1.984375 15.296875 1.78125 15.296875 1.484375 \n",
"C 15.296875 0.390625 16 -1 18 -1 \n",
"C 21.09375 -1 21.59375 1.6875 21.90625 3.171875 \n",
"z\n",
"\" id=\"CMMI12-28\"/>\n",
" </defs>\n",
" <g transform=\"translate(333.480503 26.7)scale(0.1 -0.1)\">\n",
" <use transform=\"translate(0 14.943915)scale(0.996264)\" xlink:href=\"#CMMI12-66\"/>\n",
" <use transform=\"translate(74.10802 0)scale(0.697382)\" xlink:href=\"#CMMI12-115\"/>\n",
" <use transform=\"translate(134.44505 14.943915)scale(0.996264)\" xlink:href=\"#CMSY10-33\"/>\n",
" <use transform=\"translate(261.745733 14.943915)scale(0.996264)\" xlink:href=\"#CMMI12-28\"/>\n",
" <use transform=\"translate(315.236898 14.943915)scale(0.996264)\" xlink:href=\"#CMMI12-28\"/>\n",
" </g>\n",
" </g>\n",
" <g id=\"line2d_18\">\n",
" <path d=\"M 305.480503 37.137178 \n",
"L 325.480503 37.137178 \n",
"\" style=\"fill:none;stroke:#4daf4a;stroke-linecap:square;stroke-width:0.6;\"/>\n",
" </g>\n",
" <g id=\"line2d_19\"/>\n",
" <g id=\"text_17\">\n",
" <!-- $B^0\\to\\tau\\tau$ -->\n",
" <g transform=\"translate(333.480503 40.637178)scale(0.1 -0.1)\">\n",
" <use transform=\"scale(0.996264)\" xlink:href=\"#CMMI12-66\"/>\n",
" <use transform=\"translate(79.135523 36.153639)scale(0.697382)\" xlink:href=\"#CMR17-48\"/>\n",
" <use transform=\"translate(139.290846 0)scale(0.996264)\" xlink:href=\"#CMSY10-33\"/>\n",
" <use transform=\"translate(266.591529 0)scale(0.996264)\" xlink:href=\"#CMMI12-28\"/>\n",
" <use transform=\"translate(320.082695 0)scale(0.996264)\" xlink:href=\"#CMMI12-28\"/>\n",
" </g>\n",
" </g>\n",
" </g>\n",
" </g>\n",
" </g>\n",
" <defs>\n",
" <clipPath id=\"pd8980b1118\">\n",
" <rect height=\"217.44\" width=\"334.8\" x=\"45.541911\" y=\"10.7\"/>\n",
" </clipPath>\n",
" </defs>\n",
"</svg>\n"
],
"text/plain": [
"<matplotlib.figure.Figure at 0x7fc3d7b108d0>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"flavio.plots.pdf_plot(dist_s, -4e-3, 10e-3, col=1, label=r'$B_s\\to\\tau\\tau$')\n",
"flavio.plots.pdf_plot(dist_d, -4e-3, 10e-3, col=2, label=r'$B^0\\to\\tau\\tau$')\n",
"plt.xlabel('branching ratio')\n",
"plt.ylabel('likelihood')\n",
"plt.legend()"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"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.5.2+"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment