DavidMStraub/gist:8206886d672f78e39955284f1eda445c

## gistfile1.txt
{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#  Demo of `flavio` plotting functions"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This still needs to be properly documented ..."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Prerequisites"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import flavio\n",
    "import flavio.plots\n",
    "%matplotlib inline\n",
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## $B\\to K^*e^+e^-$ vs. $B\\to K^* \\mu^+\\mu^-$ central values at low $q^2$"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {
    "collapsed": false,
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.text.Text at 0x7f57d3fae6d8>"
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXsAAAERCAYAAACaUQc3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAG8pJREFUeJzt3XmUVPWd9/H3r9beVxaBlp2ARA2rqGBsA8jEJaJEjclJ\nzjFHHJNH53nmRGPGZeycyUQ0xjg5mZwkosnMuGQw+sTouIBga+ZxGRUNoMEFEAERBLrpvdbf80dV\nawlNV3V1V9+qup/XOXWqbnVV15er91O3f9s11lpERKS4eZwuQEREck9hLyLiAgp7EREXyGnYG2Mm\nGWPWGGOuzeXniIhI/0wmHbTGmFXW2h+kbK8AWoHJ1tq7+3nfROCQtbZt8KWKiEi20p7ZG2NWAitS\ntmcD1lq7Prk961jvtda+D9QbY1YaY6oHX66IiGQjbdgnz9y3pzx1KYmzepLPL4HE2b4x5qKU+2pj\nzApr7Q7gVWDlENcuIiIZ8mX4OpPyuAY4lLJdD2CtffioNxnzmjFmMTAJ+EO2RYqIyOBkGvZZSTbj\nvJ/LzxARkfQyHY2T2ovbAtQlH9cAB4e0IhERGXLZNOOsAeYCG4DJwLpsP9wYo7UaRESyYK016V/1\nqUxG46wA5hpjrkh+wOvJ5xcDLdbaN7IptJe1Vrchut1yyy2O11BMN+1P7ct8vWUj7Zm9TXS8PnzE\nc6uz+jQREXGE48slNDU10dzc7HQZIiJ5r7m5maampqzem9EM2lwxxlgnP7/YNDc309jY6HQZRUP7\nc+hoXw4tYwx2gG32CnsRkQKTTdg73owjIiK553jYq81eRCQzarMXEXERNeOIiEifFPYiIi7geNir\nzV5EJDNqsxcRcRG12YuISJ8U9iIiLqCwFxFxAcfDXh20IiKZUQetiIiLqINWRET6pLAXEXEBhb2I\niAso7EVEXMDxsNdoHBGRzGg0joiIi2g0joiI9ElhLyLiAgp7EREXUNiLiLiAwl5ExAUcD3sNvRQR\nyYyGXoqIuIiGXoqISJ8U9iIiLqCwFxFxAYW9iIgLKOxFRFxAYS8i4gIKexERF1DYi4i4gONhrxm0\nIiKZ0QxaEREX0QxaERHpk8JeRMQFFPYiIi6gsBcRcQGFvYiICyjsRURcQGEvIuICCnsRERdQ2IuI\nuIDCXkTEBRT2IiIu4HjYayE0EZHMaCE0EREX0UJoIiLSJ4W9iIgLKOxFRFxAYS8i4gIKexERF8if\nsH/8cWhtdboKEZGilD9hf9tt8MwzTlchIlKU8ifsAbZudboCEZGipLAXEXGB/Ar7t992ugIRkaKU\nX2G/dSto+QQRkSGXX2Hf0QF79zpdhYhI0cmvsPf51G4vIpID+RX2U6cq7EVEcsCX6w8wxqwEtgGH\nrLVv9PviGTPUSSsikgMZndkbY1Ydsb3CGLM4GeT9vW8lsM5auyFt0ANMn64zexGRHEgb9snAXpGy\nPRuw1tr1ye1Z/bx9LjAn+eUwO201OrMXEcmJtGFvrb0b2J7y1KVA7yI224El8MnZ/kUp99VAC/AM\nsBG4JG01EyfCvn3Q1TWgf4SIiPQv0zb71Mtf1QCHUrbrAay1Dx/1JmNuBZYCFvh12k/xemHKFHjn\nHZjV3x8MIiIyEDntoLXWtgFHfQn0q7cpR2EvIjJkMg371GmtLUBd8nENcHAwBXxypfSdO2l8/XUa\n1UkrIvIZzc3NNDc3D+p3GJvB8gTGmLXW2rOTj2cDc621q40x15EYbZN+pE3fv9d+8vlnnAE//jG8\n/z48+SQ88EA2v1JEpOgZY7DWmvSv/FQmo3FWAHONMVcAWGtfTz6/GGjJNuiPSWf2IiJDLm0zTrLj\n9eEjnls9VAU0NTXR2NhIY+8T06cn2uzjcfDk1wRfEREnDaY5J6NmnFzpsxnnjDNg7Fh4+WU4/njH\nahMRyVc5acZxhJpyRESGVH6GvWbSiogMqZwvhJbOUW32kAj7t95yqCIRkfxUfG32L70EV14JmzY5\nVpuISL4qnjb7efNg1y746COnKxERKQr5GfY+H5x1FjzzjNOViIgUhfwMe4ClSxX2IiJDJD87aAGW\nLIEf/QisBTOgpikRkaJUHB20ixbBrbcmOmghEfKTJsETT8DMmY7VKCKSb4qngxYSZ/NqyhERGRL5\nG/aQaMpZt87pKkRECl5+h/3ixfD88xCJOF2JiEhBy98OWoARI2Dq1MSiaIsWDXNlIiL5pTg7aHtd\nfz2UlMAPfzj8BYqI5KHC76Dta4ilOmlFRAYtv8K+L4sWwZtvwv79TlciIlKw8j/sS0rgwgvhvvuc\nrkREpGA5HvZNTU3pOxwuvxzuuScx0UpExKWam5tpamrK6r351UG7alXfo26shc99Du6/H045ZXiL\nFBHJM4XfQXssxiTO7u+91+lKREQKUmGEPcC3vgVr1kBXl9OViIgUnMIJ+4YGOPVUeOQRpysRESk4\nhRP2AN/+tppyRESyUFhhf/75sHkzbN/udCUiIgWlsMI+GEy03f/iF05XIiJSUBwP+4zG2af63vfg\nd7/TxchFxHWKfpz9P/9zYuTl2LHJJ/7u78Dvh5/+dPiKFRHJE0U7zv6mm2DuXFi7NvnE9dfDb38L\n+/Y5WpeISKEoiLCHxOTZyy+Hm2+G6Ohx8PWvwx13OF2WiEhBKJiwP+ss2LgR/vxn+Pu/B37wg8R6\nOVoNU0QkrYIJe4DRo+GPf4THHoM/bWyAr30NfvITp8sSEcl7BRX2ADU1iSadK6+EPZfflGi7f+cd\np8sSEclrBRf2AAsXwne/C9+8fiyxf7gJvvMdLX8sItKPggx7gBtvhFgMbuu+BlpadHETEZF++Jwu\noKmpicbGRhoH+D6vN5Hv8+d7WfjDBzjzukY45xyor89BlSIizmtubh7YJNQU+TOpauFCuO22PidV\nGQPxeN/XI3/66cT6aK9++WbG2A8TI3RERIpY0U6q6s+yZbByJVz29i1E126Ap55yuiQRkbxT8GEP\niYlWwTIfN5/5XGLm1Z49TpckIpJXiiLsvd7EcMwH/jyeNY2/hMsug2jU6bJERPJGUYQ9wIgR8Oij\n8L+eWc5Lodlwyy1OlyQikjfyK+z76oEdgFmz4He/M1y080523LMBnnxyiAoTESls+RX2Q+Dcc+GG\nm7ycW7qe1m9eA3/5i9MliYg4rujCHuDqq+HsC8r4ysgX6DznYti92+mSREQcVZRhD3DnnTDl1FFc\nULaO7mXL4fBhp0sSEXFM0Ya9xwOrV8PIeeP5aue/EV5+CXR3O12WiIgjijbsITEk89//3RCYNZNL\n3/sRofO/Cl1dTpclIjLsijrsIXGp2t//p4G5c/nKW7fScc4lCnwRcR3Hw76pqSnrhX0yFQzCQ3/w\n0PA3J7L0zbs4tOwy6OzM6WeKiAy15uZmmpqasnpvfi2EdvvtifujXnfshdAGwlr4/rVxnrx3L09N\n+Fsanr4ncfkrEZEC4sqF0AbCGPjJTz1cfuNYFmx/gJfnfAfeftvpskREcs5VYd/re9cafv1AFee3\n3899p/w8cRVzEZEi5sqwBzjvPNjwQin/WPYTrvubTUT+5Ze6tKGIFK2CCfvBttf35cQT4X82l7Fl\n3uU03nA6H1z4v9VxKyJFqWDCPldGjID/eraMr/xgJvOf/icem3k9bNnidFkiIkPK9WEPidm2198c\n4JF1VVzdsYqr5r9K+6p/TQwBEhEpAgr7FAsXGTZtryB83kV8oWk5zfOv0yJqIlIUFPZHqK6Gex+q\n4uf/eRzfePcWvjttLa13rNZZvogUNIX9MZx3gZctO6uIn7+cmTcs58EZP8RuedPpskREsqKw70dt\nLfxqTR0PP1vHqq5rWDrnAFu+dbuWSxaRgqOwz8BpCz28umMEFzTN4UsPXcVVYx5l/533QSzmdGki\nIhlR2GfI74drbqhk654qSpcvY+b15/Gjsb+kfc2TmowlInlPYT9AdXXwswdG89Jb1fz1hAuZ9o35\n3DXl5/Ss05ILIpK/
8ifsC+zseOo0w/3NDax7pZZnay5k6pen8rOp/0rnE88V3L9FRIpfTsPeGLPC\nGLPGGPMrY8y1ufwsp5w0y8ujG8fz2AsjeKH2XCZ/5fP8eNLdtDz4lIZrikjeyCjsjTGrjtheYYxZ\nbIxZmeatr1lrLwHWAL/JssaCMPsUPw+9MpFnN9bw9vilTPnmaVwz4kG23boGenqcLk9EXC5t2CcD\nfUXK9mzAWmvXJ7dnHeu91tr3e99mrW0bXKmFYebJPv7t+Uls2VlFxbKFLLhlGRdUN7P2st8S/0Cz\ncUXEGWnD3lp7N7A95alLgdbk4+3AEvjkbP+ilPuq3ueBgxlVk4ulLR0ydpzh1gcnsrOlmnNvnMV1\n687mhEnd3PWF33LwD8+qiUdEhlWmbfapKVwDHErZrgew1j5srX0k5b73TH4Sn/2ycJXycrjyH4/j\njY/Hce9T43i1ZCFTvjaPy6oeZ/3l9xHfucvpEkXEBXI+Gsdae4dbmnD6YwwsXFrGfS9/jh0fV7Lw\nu7O49rEvMnGy4YaJD7D1tke1lr6I5Iwvw9eljiVsAeqSj2vItInmGD65UvquXTRu3Ejj6acP5tcV\nhNpauPr28Vx9O2x6JcR//NMcvnTLSMbe+A5fO+lNLrl6NOO/eSYEAk6XKiJ5oLm5mebm5kH9DmMz\nGBNujFlrrT07+Xg2MNdau9oYcx2wzlr7RlYfboz95PNPPx3uuCNxf9Trin/oejQKzz3ayu9/9iGP\nvDyOGWzlq7O2cdFVo5jwjUVQUuJ0iSKSJ4wxWGsH1MmZyWicFcBcY8wVANba15PPLwZasg16+Syf\nDxavqOHu/57J3s5qbrhnCpu9s5h71TzmV7zFj0/+PVtufwLb0pr+l4mIHCGjM/ucfbjO7NOKROC5\nP7bwp199yKMvjsTX08n5DRs59xz44tUnEzxxmtMlisgwy+bMPtM2+5xpamqisbGRxmP83K0h38vv\nhyUX17Lk4lr+xcKml7p4/BcncMv/DbDlN6NpLF3PsnkHWfb1eqZ8fQGmssLpkkUkRwbTdp/3Z/bW\nJq4R6/bQ78uBjy1rV+9k7R/aWbtlDCWRdpaM3sziM8Kc9a3xjFo2O/FtISJFJZsze4V9kbAWtvxP\nF+vv3cmG9XGef3884+1OzmzYRuOiGF/8xvGMXPIFjfARKQIF2YwjQ8MYOGlBGSctOIH/Q2J0z2sb\njue5++He5z18+/fHM9bu4Izj3mXR/BALl49k8vKTMbU1TpcuIsPA8bBP12Yv2fH5YMHZ1Sw4u5rv\nkwj/zS+M5b8fjPF4c4wbrjqO8LcjnF65ntNmtLLgzBLmXTyJirnTwet1unwR6YPa7CUru7ZHeHHN\nLl58+jAvbynnLwcbmGx2cMqoncw7sYf5X6rkpAsmEzxhclGtWyRS6NRmL4MSDsOm51t55dEPefWF\nMK++V8277ccxw/MOc0btZs7MELMWlnPyeeOpmD0t8eeDiAw7hb0Mua4u2PzcITY+8REbXwrxxrYq\n3mwdSwO7+UL1Tk6a1M7Js32c1FjPpMWT8YwZrb8CRHKsIDto1Waf38rKYMGX61jw5bpPnotGYetr\nY9n8lGXTi53c/YyfzfeP5GCoguneTcys+4gTxndxwoleZpxWy9QzxxGYNkF9ASKDpDZ7yQtthy1b\nX2zhzQ372PpGD1vf87F1Xw07u0bSYPYwvXwP0487zLTJUaadWMK0BXU0nHY83oYx+mtAZADUjCN5\nKRyG7Vu6eOfP+3j7lTbefTvOu7tLePdQPQfDlUwwHzClfB9TRrYxZUKUydP9TDq5kkkLRlE+cwKU\nljr9TxDJKwXZjCPFLxCAGXPKmDFn0lE/6+qCHVsaeO//+dn2RjvvvRdj7ZN+tv1HDTu7R1Jp25jo\n28rEqoNMHNXFhIY4E6b6mXBiJeNn11P9+QaornbgXyVSWBT24qiyMvj8KeV8/pTJR/0sHod9ewPs\neAV2vlHGzq3dbN4e5/H/CvDBfZXs7ByJ10Y43vMWx5ce5PjadhpGhWloMDRMDjBuegXjTqylesYY\nTH2dmorE1RwPe3XQyrF4PDBmnIcx40Zx+vJRR/3cWmg5ZNn1lmH3pjJ2/bWDXdvCPP++YfcrQfa0\nVbCnu45Y3DDOvMfY4EHGVHQwpraHsaOjjBnr4bgJQcZMKWPMzFpqpo3EjByhjmTJW+qgFelHezvs\nebeLD99sYe877Xy4o4e9u2Ps3edlb0uQve0V7OuppjseZDT7GO07xOjSNkZVdDG6NsyokXFGjfYw\nqiHAqAmljJxUwYipNQTGjYTKSv3FIMNObfYifais7O0zKOv3dT09sG/PWD7aWsb+99rYt6OLfbsj\n7Nwf55XNHvY9H+TjzlI+7qnkQKSacjoZyTZG+A8zsrSd+vIeRlSFqa+NU19vqB/lpX5MgPqGUuqO\nL6duYhUlY2qhpkZ/Pciwy5+w16m7OKykBCZM8TFhyghgRL+vjcehtbWGA7tLOLCtko93dHBwTw8H\n9kY4sD/Gu3sMh/7q5WB7kINdpRwKl3MwUoWPKHXspc5zmLpAO7XBbmrLQtRWRKipilNba6mp9VA7\n0kf1CD81x5VQM7aM6rHlVIypxNRUQ3m5/pqQAcufsBcpIB4P1NVBXV0Jnzu5BBid9j3WQkdHgJYD\nQVp2V9HyQTuH9nTT8lGIlv0RWg7G+WsrtH5oaGn30doV5HAoSGuojNZoOT02SBVtVLOTam8H1b4u\nqgPdVJeEqCyNUlUeparCUlVlqawyVNV4qaz1UVnnp3JEMHEbVUrFqDL8tRWJP3lKSvTF4RIKe5Fh\nYkwiXysrvYyfVAVUDej9kQi0tdVy+EAFbXs7Oby3i7b9PRzeH6LtUJS2lhiHW+LsbYf2/Ya2Tg/t\n3T7aegK0h4O0R0roiJbQHi/HT4RK2qngIyo8XVR4e6jw91DuD1MRiFAejFJeEqOiLE5ZqaW8HMor\nTOJW5aWs0ktZlY/y2gBl1X7KaoOU1QYprQniry7DlJcl5kcEAvoyyROOh/1nRuPofwqRY/L7ob4e\n6uv9ML0GyO5aBNZCd7eXzs4SOlprad/XReehEB0Heug4FKajJUzn4Rgdh2N0tsc51GHZ1QmdBwyd\n3R66Qh66Qj46wz66I366on46oyV0xwN0xhMT4MroopTWxL0JUeoJUeoNU+KNUOqLUOqPUOqPURqI\nURqMURKwlASTt5LEHxwlpYaSUkOw1ENJmYeSci/BMi8lFT6CZV6C5T6CFf5Pb5UBAhUBTEkw8SUT\nTN4X0YJ9xTEa57TT4M47E/cpNBpHpLBEIonJct3d0N0epaslRHdriO7DYbrbwnS3RenpiNLdHqW7\nM0ZPZ5yebkt3V5ye7kRHeSgE3SFDKGzoCXvoCXsJRT30RLz0RH2EYr03Pz1xP2HrJxT3EyFAgBAB\nwgQJEex9bMIETJSAJ0LAE03eYgS8yXtfjIA3jt8bJ+CL4/dZ/D5LwJ+49/stfl/iC/eTWwD8foPP\nb/AHErfex76AB3/Q4PN78Ac9+AIpt6D36PugF68/+bjEh9fvwRv0YXzeRGe+z5eYPJi8zKhG44iI\n4/z+RC5VVwPH+UjETPmwfHY8DuFwkHA4SChUSSgE4e4Y4Y4woY4Ike4ooc7ELdITI9wVJRKKE+qK\nEQnFCfckbpFQnEjYEg5ZImFLKAIdEUskDJEQRDogEjVEohCNmuRjQzRuiMYMkZiHaMxDJHmLWUM0\n7iEa9xCJe4lZD9G4l4hNPrbeT24x6yWKlzhePMTwEcVLjD/dtpkl35+T9b5R2ItI0fB4+KQZ6FNe\noDR5KxzWQizmJRbzEo1CMJh90IPCXkQkLxmTaL3x+RLdD4PlGfyvyC211YuIDF7ehz1okI6IyGA5\n3oyjhdBERDJT1EMv4/FEm1U87kCBIiJ5KJuhlwXRjCMiIoOjsBcRcQGFvYiICyjsRURcQGEvIuIC\
nCnsRERdQ2IuIuIAmVYmIFAhNqhIRcRFNqhIRkT4p7EVEXEBhLyLiAgp7EREXUNiLiLiAwl5ExAXy\nJ+x1/UERkZzJn7AXEZGcya+w18VmRURyIr/CXkREciLvw15N+SIig1cQC6GpdUdEpFgWQjv1VLjr\nrsR9ilgMAoHEvYiIaCE0ERE5BoW9iIgLKOxFRFxAYS8i4gIKexERF1DYi4i4gMJeRMQFFPYiIi6g\nsBcRcQGFvYiICyjsRURcQGEvIuICCnsRERdQ2IuIuEBO17M3xlQDk4FaYLu19v1cfp6IiPQtozN7\nY8yqI7ZXGGMWG2NWpnnrPGAJ0ALUZFeiiIgMVtqwTwb6ipTt2YC11q5Pbs/q5+2vAlOA3wDbB1eq\niIhkK23YW2vv5rNBfSnQmny8ncSZe+/Z/kUp99XAJdbaq5Kv+duhLV2OlO3lyqRv2p9DR/vSeZl2\n0KZe/qoGOJSyXQ9grX3YWvtIyv1hYJsx5ksk2u3XDUnFckw6oIaW9ufQ0b50Xk47aK21G3L5+0VE\nJDOZntmnXpW8BahLPq4BDg5pRSIiMuSMtTb9i4xZa609O/l4NjDXWrvaGHMdsM5a+0ZWH25M+g8X\nEZGjWGtN+ld9Km0zjjFmBTDXGHOFtXa1tfZ1Y8xcY8xioCXboM+mWBERyU5GZ/YiIlLYtFyCCEdP\nHDziZ5lOIpSkNPtzVfJe+3MYDVvYpztgdEANTAb7UwdUho6cOHjEzwYyiVDof38mXWmMeRfYNkwl\nFTRjzMrkrc8v0Eyzc1jCPt0BowNqYDLcXzqgMtTHxMFUfU4ilGNLsz8BrrDWTtPQ7PSSfaPrkvt0\ncnLeUurPM87O4TqzT3fA6IAamEz2lw6oodHnJEIZlLrkmeh1ThdSACbz6fG9PbmdKuPsHK6wT3fA\n6IAamEz2lw4oyUvJUX3rgfojz1Tls6y1d1trVyc355BYbyxVxtmpDtoipQNqyGgS4RBKtj1flNw8\nyNFnqtKHZHPNa4MZ6j5cYZ/ugNEBNTD97i8dUFn5zJyP5EJ+AGv4dP9NBp4ZzqIK2LH25zY+3Yf1\nHH2mKn1bbK39hz6ezzg7hyvs+zxgdEBlLd3+1AE1AKkTB1OefgbAWvt68jWDnkToFmn25wZgafI1\nB7Q/0zPGrLTW3pF8vDh5P+DsHLZJVcn/8DuASb1tUMaYV6y184/1czm2DPZn79C3Sb3/o4hIYUmG\n+xoSZ/C1wMXW2g3ZZKdm0IqIuIA6aEVEXEBhLyLiAgp7EREXUNiLiLiAwl5ExAUU9iIiLpDTC46L\nFLrk5JXexaWWWmuvcrIekWzpzF6kf5eQWEL2YaD1iFmhIgVDZ/Yi/UiuI95rMrDWqVpEBkMzaEUy\nkJy2rqU8pGAp7EXSSC4vi7X2dWPM7N7F0UQKiZpxRJKSF3p5jcSCU5OBXwNTgIeAbcaYOuB65yoU\nyZ7O7EUAY8yvgLXW2keMMZOSj6c5XZfIUNFoHHG9ZLgvttY+knxqDrqmghQZhb1IYhx9argvBdY5\nVItITijsRRIXbG5N2b4EndlLkVEHrbietfZhY8yS5IXZpwAHrbVtTtclMpQU9iKAtfY7AMaYWnRW\nL0VIzTginzUftddLEVLYiyQZY1YCK0kEvkhR0Th7EREX0Jm9iIgLKOxFRFxAYS8i4gIKexERF1DY\ni4i4wP8HRmIaCDiolRAAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f57d3f8b940>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "flavio.plots.q2_plot_th_diff('dBR/dq2(B0->K*ee)', 0, 2, color='red')\n",
    "flavio.plots.q2_plot_th_diff('dBR/dq2(B0->K*mumu)', 0, 2, color='blue')\n",
    "plt.yscale('log')\n",
    "plt.xlabel(r'$q^2$')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Experiment vs. theory for binned observables"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "th_bins = [(0.1, 2), (2.0, 4.0), (4.0, 6.0), (15.0, 19.0)]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXIAAAEECAYAAAA1X7/VAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAD5hJREFUeJzt3cFyHMd9x/HfP+VzCANH62CsS3dCK6RAXbAlMC8QWOYL\nEMrRVS7TkU6CTmEqTFWOsYkHcFimHiAkUaOLyDLFIu4RQR1snwIC5Av8c9gGMFotMLvY3pn573w/\nVVua2ZmdaYx3f2529/SYuwsAENffNV0AAMBsCHIACI4gB4DgCHIACO4nTRdgFma2LemWpDeSvnP3\new0XCQBqZ20YtWJmd939s9L6tqQTST13v3/J537u7t+b2ceSvnX3dzUUFwBapfGmFTPbkbRdWl+T\n5O7+JK1fv+iz7v796ccIcQBd1XjTirvfN7Nflt66Jel/0vKhpJuSDlIt3SVZ+u9jd3+X3n9VZ5kB\noE3mHuSphn3i7q/T+pak5yM1aCstL2nY5n1qRZLc/eEFp1iV9ChfiQEglrk3rbj7S0lbZraaQtxz\nNoO4+z2aVQB02URBbmZ3L9m2bWZbqa17LHffk/QvklbdfX/cLqXlY0nLaXlJ0tEkZQSArqoM8tHO\nyJFtE3VMppr4g7S8Om6X0vIDSb203JP0uKqMANBllUGehv8dXrD5lobDBKXzjskfKIX9fqqZb5nZ\n35e2b0vqm9ntdL6X6f0tScfufjDF3wMAnTNrZ+fYjsmy02Aure+NrD+U9PCyfQAAF2t8HDkAYDaz\nBjkdkwDQsEmbVuwHK2bX3P2thh2TfUn7GnZM/mg8t5k1PwcAAATk7la912SjVn7QGZk8TieZqGPS\n3Xllen3xxReNl2FRXlxLrmebX9OorJH7+M7I9dIyHZMA0CA6OwEguMYnzcJ0BoNB00VYGHVey6Io\nVBTF2fLpuQeDwcL8b7oof0dEc5+P3Mx83ucAIjGzqdtA0T3pe5KnsxMA0G4EOQAER5ADQHAEOQAE\nR5ADQHAEOQAER5ADQHAEOQAER5ADQHAEOQAER5ADQHAEOQAER5ADQHAEOQAER5ADQHAEOQAER5AD\nQHC1Pert9bd/retUP7L64c8aOzcAzBs1cgAIjiAHgOAIcgAIjiAHgOAIcgAIjiAHgOAIcgAIjiAH\ngOAIcgAIjiAHgOAIcgAIjiAHZlAUxUKcA7ER5MAMCHK0AUEOAMHVNo0tsIiKopCZTf25aT6zubk5\n9fHRLZU1cjPbNrMtM9up2H47f/GAdhsMBnL3qV6Sptp/MBg0+0ei9S4NcjNbk+Tu/iStXx+z/TBt\nfz26HQAwf1U18luSTtLyoaSbY/b5t/Tfnrsf5CoYEEEdtWVq5KhSFeRLkt6U1lfKG939paRDM3sj\n6Shz2YDWI8jRBjONWjGza5K+k3Rb0n0z+3mGMgEAplA1auVY0nJaXtKPa92fSvq9u78zsxNJn0j6\n99GD7O7u6vhv7yRJG/0b2uh/NFOhAWDRFEVx5XsG7LQXfezGYWdm3933zOyOpEfufmBm19z9rZn9\n1t3vlfa/7e57I8dwd9frb/96pQLmsPrhzxo7NzDKzHTZ7w6Qzr4nE41TvbRG7u4vzaxvZluSjkud\nmY8lrbv7vRTwryQtj4Y4AGD+Lq2RZzkBNXLgB6iRYxLT1Mi5RR8AgiPIASA4ghwAgiPIASA4Zj8E\n0Kjy+OmiKM7uZB0MBtzVOiFGrQA1Y9TKxbg25xi1AgAdQpADQHAEOQAER5ADQHAEOQAER5ADQHCt\nHkf+7MU3evbiaVp+qo3+DUnMaQ4AZWHGkffW39Ph879c6bOMI0ebMFb6Ylybc4wjB4AOIcgBILhW\nt5E3jTkgAERAG/mEaLtDLnyXLsa1OUcbOQB0CEEOAMHV1rSS4TiN/pOr6fNjcfBduhjX5hxNKwDQ\nIYxaAZBFrofHNPkQmlzqvgmRGjkABEeNHKhB+Z6Ezc1N7e7uSuKeBORBkAM1ILAxTzStAEBwBDkA\nBEeQA8ji2Ytvmi6CpPaUo04EOYAsTh8C07S2lKNOBDkABMeoFQBZPHvxVL3192Y+zqzH+IcPNmYu\nQzQEOYAsNvo39Mff/2mmY8wyXfWp//zDf8z0+YhoWgGA4AhyAFls9G80XQRJ7SlHnQhyAFls9D9q\nugiS2lOOOlW2kZvZtqQTST13vz9m+5qknqTlcdsBAPN1aZCnkHZ3f2JmPTO77u4HI7t97u6/MrPf\nXrAdSRum56x7ek0A81fVtHJLw9q4JB1KulnemGrrf5Ykd79HiANA/aqCfEnSm9L6ysj2dUkrZrZm\nZneylgwAMJEcnZ1H7v5SOquhj3U6F3MTmjw3AMxbVWfnsaTltLwk6Whk+5GGTS7SsAnmQ0kPRw+y\nu7uroijO5mSue17m03MDQFuVHz4yraogfyCpL2lfw5EpjyTJzK65+1tJf5J0WgtfkvR83EF2d3fP\nXgCAHxut5H755ZcTf/bSppVSk8mWpONSZ+bjtP21pJPUpLLs7l9ddKyiKGRmV36lclzpRdMKgEVW\n2Ubu7nvu/sTd90rvrY9sf+jun192nMFgIHe/8iud60ovmlUALDLu7ASA4GoL8iZrxdTIASyyTgR5\nW3TxEVQA5q8TTStt6ezs4iOoAMxfJ4IcABZZJ54QdDr0cVazHqOLj6ACMH+dCPLBYDBz84qZnQ2D\nvKpff/qbmT4PtFmumTWZoXN6nQjy47+9yzKFbBumoQWAUZ1oI2/Lo5/aUg4Ai6UjQd6ORz+1pRwA\nFksnghwAFhlBDgDBEeQAEBxBDgDBdWL4IYD2Kj8ZZ3Nz8+wBNE08TSwqghxAowjs2dG0AgDBEeQA\nEBxBDgDBEeQAEBxBDgDBEeQAEBxBDgDBMY68RkyYD2AeqJEDQHDUyAGEV77NvyiKsztFu3LXqM36\nHMrKE5h5jnPM8szMHI9o662/p8Pnf5npGDStAPOX4/m6bZD+jome+E7TCgAER5ADQHCtbiNneksA\nqNbqICewAaBaq4M8gmcvvtGzF0/T8lNt9G9Ikjb6N7TR/6jJogHoiDCjVmZR16iVqn0YtQLMXxdH\nrXQiyHOY5MuxKF+gcbo+ThdxLMrvkCCfg64HeVlX/k7EtCjfz6zjyM1s28y2zGynYr87kxYQAJDP\npUFuZmuS3N2fpPXrF+y3Jelm/uIBAKpU1chvSTpJy4cirAGgdaqCfEnSm9L6yugOZraWauwTteUA\nAPLKMY78pxmO0UrcWQpMrjyaqevqvhZVQX4saTktL0k6Km9MtfH9tBq/m3gEgQ1MjiA/17YgfyCp\nL2lfUk/SI0kys2vu/lZSz8xWNWxyWTGz6+5+MHqQ05qsRDgCwDjlFoBpVY4jN7Pbkl5LWnX3vfTe\nc3dfL+2zI+l3kj4ZDfJFGUc+iUUZv1qlK38npjMYDPT11183XYxW2NzcvHIon8o6jtzd99z9yWmI\np/fWR/a57+7vj6uNA+iGwWAgd2/8JanxMtTd6sB85AAQXJggn/WfKRjiOmJe6Ps6R438AgRQHlxH\nzAtBfo4gBwBMJcyDJYqikFn7bx5texk3NzebLgKAzMLUyNvSI9723vK29aYDmL8wQQ4AGC9MkFOT\nzIPrCCwenhCU0bzueGzbY9a4sxNttijfTx711pA6vkCzniPHZD45/k4mWMK8dDHIwzStII+2jCNv\nSzmARUCQA0BwYcaRI49c4/FnPQbj2YF8qJF3TI7x+NLs4+VpHwfyIcgBIDiaVjqmLTXhtpQDi6Hr\nz9dl+GFGEYYfLkoZgEXH8EMA6BCCHACCI8gBIDiCHACCI8gBIDhGrcyo7pkJ2zBipA1lABYdsx8u\nsDaEaBvKACw6hh8CQIcQ5AAQHEEOAMER5DVqw8MU2lAGAHkR5DVqQ4i2oQwA8iLIASA4prGtURue\nzsOTeYDFQ428Rm14Ok8X5mYGuoYgB4DgCPIataE23IYyAMiLW/SDacPt8W0oA7DouEUfADqkskZu\
nZtuSTiT13P3+mO07afEX7v7ZmO3UyDNqqjZc9yyPQNdlm/3QzNYkrbr7Vymwn7v7QWn7lqRX7v69\nmT2Q9F/uvj9yDII8I5o1gG7I2bRyS8PauCQdSro5sr1Xeu8wrQMAalR1Q9CSpDel9ZXyxpGmlg8k\n/TFTuQAAE8rS2ZmaYF6Um10AAPWoqpEfS1pOy0uSji7Yb8vdP7/oILu7u2fLdI4BaJs2dOaXyzCt\nSTo7++6+Z2Z3JD1y9wMzu+bub9M+O6dNLGa25e5PRo5BZ2dGdHYC89WW31i2zk53f5kOuCXpuNR0\n8rj0/l0z+87MjiQ1/9cDQMdwZ2cwbaktAIuqLb8x7uwEgA4hyAEgOIIcAIIjyAEgOIIcAIIjyAEg\nOIIcAIIjyAEgOG4ICqAN80AAXRHxhiCCHABKIgY5TSsAEBxBDgDBEeQAFspV5/SOXAaCHMBCIcgB\nAOFUPeoNAEIpikJmEw32uNCsn9/c3Jzp89OiRg5goQwGA7n7lV+SZvq8u9d+fwdBDgDBEeQAFkob\n7nauuwzc2QkAJdzZCQCoHUEOAMER5AAQHEEOAMER5AAQHEEOAMER5AAQHEEOAMER5AAQHEEOAMER\n5AAQHEEOAMER5AAQHEEOAMExjS2AziuK4uyByUVRnM0nPhgMGpvffJppbAlyAGihaYK88uHLZrYt\n6URSz93vT7sdADBfl7aRm9maJHf3J2n9+jTbkd/pP/8wO65lXlzP5lR1dt7SsLYtSYeSbk65HZnx\nY8mHa5kX17M5VUG+JOlNaX1lyu0AgDlj+CEABHfpqBUz+1dJj9x9P3Vqrrr7vUm3p30YsgIAV5Br\n1MoDSX1J+5J6kh5Jkpldc/e3F22/SkEAAFdzadOKu7+UJDPbknTs7gdp0+OK7UArmdndkfVtM9sy\ns52myhTZmOt5N/2X61mjud0QxPjyvMzsrrt/ZmY7XM+rSeHyO3d/P62vadgc+FXa9pzKyORGr2d6\n742kI0n/7O77jRUuoNL/+f3C3T9L702Uo3Pp7GR8+Vx8amb/K+lV0wWJKv0QDktvMXx2BmOupyTd\ndvf3CfHppFaNR+ma9szs42lydF6jVviB5McPJI9ynw3DZ/NbTk1Vd5ouSDA9nefkYVqfOEfnFeT8\nQPLjB4LWc/e9VINcMbOPmy5PFO5+39330uoHkr7VFDnKOPIg+IFkU+4UOpa0nJaXNGzbxRWZ2Y6Z\n/VNaPdKwVokppOaUF9P21cwryPmBZMQPJKty08oDnV/LntJoLEylfD1f6fwarmhYq8R0ttz987Q8\ncY7OK8j5geTFDySDNAKgb2a3JYbPzmrM9dyX9I/p/f/jek4njUi7l5a3JP23JszReQ4/vC3ptYbD\nu/aq9sfl0o9DGnP3LIDYUnA/0LAW/lNJn6Q75ifK0bk/WAIAMF90dgJAcAQ5AARHkANAcAQ5AARH\nkANAcAQ5AARHkANAcAQ5AAT3/9Jq0wQuTLOyAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f57d3c200f0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "flavio.plots.q2_plot_exp('<dBR/dq2>(B+->K*mumu)', c='k')\n",
    "flavio.plots.q2_plot_th_bin('<dBR/dq2>(B+->K*mumu)', bin_list=th_bins)"
   ]
  },
  {
   "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.4.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}
	{
	"cells": [
	{
	"cell_type": "markdown",
	"metadata": {},
	"source": [
	"# Demo of `flavio` plotting functions"
	]
	},
	{
	"cell_type": "markdown",
	"metadata": {},
	"source": [
	"This still needs to be properly documented ..."
	]
	},
	{
	"cell_type": "markdown",
	"metadata": {},
	"source": [
	"## Prerequisites"
	]
	},
	{
	"cell_type": "code",
	"execution_count": 14,
	"metadata": {
	"collapsed": true
	},
	"outputs": [],
	"source": [
	"import flavio\n",
	"import flavio.plots\n",
	"%matplotlib inline\n",
	"import matplotlib.pyplot as plt"
	]
	},
	{
	"cell_type": "markdown",
	"metadata": {},
	"source": [
	"## $B\\to K^e^+e^-$ vs. $B\\to K^ \\mu^+\\mu^-$ central values at low $q^2$"
	]
	},
	{
	"cell_type": "code",
	"execution_count": 20,
	"metadata": {
	"collapsed": false,
	"scrolled": true
	},
	"outputs": [
	{
	"data": {
	"text/plain": [
	"<matplotlib.text.Text at 0x7f57d3fae6d8>"
	]
	},
	"execution_count": 20,
	"metadata": {},
	"output_type": "execute_result"
	},
	{
	"data": {
	"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXsAAAERCAYAAACaUQc3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAG8pJREFUeJzt3XmUVPWd9/H3r9beVxaBlp2ARA2rqGBsA8jEJaJEjclJ\nzjFHHJNH53nmRGPGZeycyUQ0xjg5mZwkosnMuGQw+sTouIBga+ZxGRUNoMEFEAERBLrpvdbf80dV\nawlNV3V1V9+qup/XOXWqbnVV15er91O3f9s11lpERKS4eZwuQEREck9hLyLiAgp7EREXyGnYG2Mm\nGWPWGGOuzeXniIhI/0wmHbTGmFXW2h+kbK8AWoHJ1tq7+3nfROCQtbZt8KWKiEi20p7ZG2NWAitS\ntmcD1lq7Prk961jvtda+D9QbY1YaY6oHX66IiGQjbdgnz9y3pzx1KYmzepLPL4HE2b4x5qKU+2pj\nzApr7Q7gVWDlENcuIiIZ8mX4OpPyuAY4lLJdD2CtffioNxnzmjFmMTAJ+EO2RYqIyOBkGvZZSTbj\nvJ/LzxARkfQyHY2T2ovbAtQlH9cAB4e0IhERGXLZNOOsAeYCG4DJwLpsP9wYo7UaRESyYK016V/1\nqUxG46wA5hpjrkh+wOvJ5xcDLdbaN7IptJe1Vrchut1yyy2O11BMN+1P7ct8vWUj7Zm9TXS8PnzE\nc6uz+jQREXGE48slNDU10dzc7HQZIiJ5r7m5maampqzem9EM2lwxxlgnP7/YNDc309jY6HQZRUP7\nc+hoXw4tYwx2gG32CnsRkQKTTdg73owjIiK553jYq81eRCQzarMXEXERNeOIiEifFPYiIi7geNir\nzV5EJDNqsxcRcRG12YuISJ8U9iIiLqCwFxFxAcfDXh20IiKZUQetiIiLqINWRET6pLAXEXEBhb2I\niAso7EVEXMDxsNdoHBGRzGg0joiIi2g0joiI9ElhLyLiAgp7EREXUNiLiLiAwl5ExAUcD3sNvRQR\nyYyGXoqIuIiGXoqISJ8U9iIiLqCwFxFxAYW9iIgLKOxFRFxAYS8i4gIKexERF1DYi4i4gONhrxm0\nIiKZ0QxaEREX0QxaERHpk8JeRMQFFPYiIi6gsBcRcQGFvYiICyjsRURcQGEvIuICCnsRERdQ2IuI\nuIDCXkTEBRT2IiIu4HjYayE0EZHMaCE0EREX0UJoIiLSJ4W9iIgLKOxFRFxAYS8i4gIKexERF8if\nsH/8cWhtdboKEZGilD9hf9tt8MwzTlchIlKU8ifsAbZudboCEZGipLAXEXGB/Ar7t992ugIRkaKU\nX2G/dSto+QQRkSGXX2Hf0QF79zpdhYhI0cmvsPf51G4vIpID+RX2U6cq7EVEcsCX6w8wxqwEtgGH\nrLVv9PviGTPUSSsikgMZndkbY1Ydsb3CGLM4GeT9vW8lsM5auyFt0ANMn64zexGRHEgb9snAXpGy\nPRuw1tr1ye1Z/bx9LjAn+eUwO201OrMXEcmJtGFvrb0b2J7y1KVA7yI224El8MnZ/kUp99VAC/AM\nsBG4JG01EyfCvn3Q1TWgf4SIiPQv0zb71Mtf1QCHUrbrAay1Dx/1JmNuBZYCFvh12k/xemHKFHjn\nHZjV3x8MIiIyEDntoLXWtgFHfQn0q7cpR2EvIjJkMg371GmtLUBd8nENcHAwBXxypfSdO2l8/XUa\n1UkrIvIZzc3NNDc3D+p3GJvB8gTGmLXW2rOTj2cDc621q40x15EYbZN+pE3fv9d+8vlnnAE//jG8\n/z48+SQ88EA2v1JEpOgZY7DWmvSv/FQmo3FWAHONMVcAWGtfTz6/GGjJNuiPSWf2IiJDLm0zTrLj\n9eEjnls9VAU0NTXR2NhIY+8T06cn2uzjcfDk1wRfEREnDaY5J6NmnFzpsxnnjDNg7Fh4+WU4/njH\nahMRyVc5acZxhJpyRESGVH6GvWbSiogMqZwvhJbOUW32kAj7t95yqCIRkfxUfG32L70EV14JmzY5\nVpuISL4qnjb7efNg1y746COnKxERKQr5GfY+H5x1FjzzjNOViIgUhfwMe4ClSxX2IiJDJD87aAGW\nLIEf/QisBTOgpikRkaJUHB20ixbBrbcmOmghEfKTJsETT8DMmY7VKCKSb4qngxYSZ/NqyhERGRL5\nG/aQaMpZt87pKkRECl5+h/3ixfD88xCJOF2JiEhBy98OWoARI2Dq1MSiaIsWDXNlIiL5pTg7aHtd\nfz2UlMAPfzj8BYqI5KHC76Dta4ilOmlFRAYtv8K+L4sWwZtvwv79TlciIlKw8j/sS0rgwgvhvvuc\nrkREpGA5HvZNTU3pOxwuvxzuuScx0UpExKWam5tpamrK6r351UG7alXfo26shc99Du6/H045ZXiL\nFBHJM4XfQXssxiTO7u+91+lKREQKUmGEPcC3vgVr1kBXl9OViIgUnMIJ+4YGOPVUeOQRpysRESk4\nhRP2AN/+tppyRESyUFhhf/75sHkzbN/udCUiIgWlsMI+GEy03f/iF05XIiJSUBwP+4zG2af63vfg\nd7/TxchFxHWKfpz9P/9zYuTl2LHJJ/7u78Dvh5/+dPiKFRHJE0U7zv6mm2DuXFi7NvnE9dfDb38L\n+/Y5WpeISKEoiLCHxOTZyy+Hm2+G6Ohx8PWvwx13OF2WiEhBKJiwP+ss2LgR/vxn+Pu/B37wg8R6\nOVoNU0QkrYIJe4DRo+GPf4THHoM/bWyAr30NfvITp8sSEcl7BRX2ADU1iSadK6+EPZfflGi7f+cd\np8sSEclrBRf2AAsXwne/C9+8fiyxf7gJvvMdLX8sItKPggx7gBtvhFgMbuu+BlpadHETEZF++Jwu\noKmpicbGRhoH+D6vN5Hv8+d7WfjDBzjzukY45xyor89BlSIizmtubh7YJNQU+TOpauFCuO22PidV\nGQPxeN/XI3/66cT6aK9++WbG2A8TI3RERIpY0U6q6s+yZbByJVz29i1E126Ap55yuiQRkbxT8GEP\niYlWwTIfN5/5XGLm1Z49TpckIpJXiiLsvd7EcMwH/jyeNY2/hMsug2jU6bJERPJGUYQ9wIgR8Oij\n8L+eWc5Lodlwyy1OlyQikjfyK+z76oEdgFmz4He/M1y080523LMBnnxyiAoTESls+RX2Q+Dcc+GG\nm7ycW7qe1m9eA3/5i9MliYg4rujCHuDqq+HsC8r4ysgX6DznYti92+mSREQcVZRhD3DnnTDl1FFc\nULaO7mXL4fBhp0sSEXFM0Ya9xwOrV8PIeeP5aue/EV5+CXR3O12WiIgjijbsITEk89//3RCYNZNL\n3/sRofO/Cl1dTpclIjLsijrsIXGp2t//p4G5c/nKW7fScc4lCnwRcR3Hw76pqSnrhX0yFQzCQ3/w\n0PA3J7L0zbs4tOwy6OzM6WeKiAy15uZmmpqasnpvfi2EdvvtifujXnfshdAGwlr4/rVxnrx3L09N\n+Fsanr4ncfkrEZEC4sqF0AbCGPjJTz1cfuNYFmx/gJfnfAfeftvpskREcs5VYd/re9cafv1AFee3\n3899p/w8cRVzEZEi5sqwBzjvPNjwQin/WPYTrvubTUT+5Ze6tKGIFK2CCfvBttf35cQT4X82l7Fl\n3uU03nA6H1z4v9VxKyJFqWDCPldGjID/eraMr/xgJvOf/icem3k9bNnidFkiIkPK9WEPidm2198c\n4JF1VVzdsYqr5r9K+6p/TQwBEhEpAgr7FAsXGTZtryB83kV8oWk5zfOv0yJqIlIUFPZHqK6Gex+q\n4uf/eRzfePcWvjttLa13rNZZvogUNIX9MZx3gZctO6uIn7+cmTcs58EZP8RuedPpskREsqKw70dt\nLfxqTR0PP1vHqq5rWDrnAFu+dbuWSxaRgqOwz8BpCz28umMEFzTN4UsPXcVVYx5l/533QSzmdGki\nIhlR2GfI74drbqhk654qSpcvY+b15/Gjsb+kfc2TmowlInlPYT9AdXXwswdG89Jb1fz1hAuZ9o35\n3DXl5/Ss05ILIpK/
	8ifsC+zseOo0w/3NDax7pZZnay5k6pen8rOp/0rnE88V3L9FRIpfTsPeGLPC\nGLPGGPMrY8y1ufwsp5w0y8ujG8fz2AsjeKH2XCZ/5fP8eNLdtDz4lIZrikjeyCjsjTGrjtheYYxZ\nbIxZmeatr1lrLwHWAL/JssaCMPsUPw+9MpFnN9bw9vilTPnmaVwz4kG23boGenqcLk9EXC5t2CcD\nfUXK9mzAWmvXJ7dnHeu91tr3e99mrW0bXKmFYebJPv7t+Uls2VlFxbKFLLhlGRdUN7P2st8S/0Cz\ncUXEGWnD3lp7N7A95alLgdbk4+3AEvjkbP+ilPuq3ueBgxlVk4ulLR0ydpzh1gcnsrOlmnNvnMV1\n687mhEnd3PWF33LwD8+qiUdEhlWmbfapKVwDHErZrgew1j5srX0k5b73TH4Sn/2ycJXycrjyH4/j\njY/Hce9T43i1ZCFTvjaPy6oeZ/3l9xHfucvpEkXEBXI+Gsdae4dbmnD6YwwsXFrGfS9/jh0fV7Lw\nu7O49rEvMnGy4YaJD7D1tke1lr6I5Iwvw9eljiVsAeqSj2vItInmGD65UvquXTRu3Ejj6acP5tcV\nhNpauPr28Vx9O2x6JcR//NMcvnTLSMbe+A5fO+lNLrl6NOO/eSYEAk6XKiJ5oLm5mebm5kH9DmMz\nGBNujFlrrT07+Xg2MNdau9oYcx2wzlr7RlYfboz95PNPPx3uuCNxf9Trin/oejQKzz3ayu9/9iGP\nvDyOGWzlq7O2cdFVo5jwjUVQUuJ0iSKSJ4wxWGsH1MmZyWicFcBcY8wVANba15PPLwZasg16+Syf\nDxavqOHu/57J3s5qbrhnCpu9s5h71TzmV7zFj0/+PVtufwLb0pr+l4mIHCGjM/ucfbjO7NOKROC5\nP7bwp199yKMvjsTX08n5DRs59xz44tUnEzxxmtMlisgwy+bMPtM2+5xpamqisbGRxmP83K0h38vv\nhyUX17Lk4lr+xcKml7p4/BcncMv/DbDlN6NpLF3PsnkHWfb1eqZ8fQGmssLpkkUkRwbTdp/3Z/bW\nJq4R6/bQ78uBjy1rV+9k7R/aWbtlDCWRdpaM3sziM8Kc9a3xjFo2O/FtISJFJZsze4V9kbAWtvxP\nF+vv3cmG9XGef3884+1OzmzYRuOiGF/8xvGMXPIFjfARKQIF2YwjQ8MYOGlBGSctOIH/Q2J0z2sb\njue5++He5z18+/fHM9bu4Izj3mXR/BALl49k8vKTMbU1TpcuIsPA8bBP12Yv2fH5YMHZ1Sw4u5rv\nkwj/zS+M5b8fjPF4c4wbrjqO8LcjnF65ntNmtLLgzBLmXTyJirnTwet1unwR6YPa7CUru7ZHeHHN\nLl58+jAvbynnLwcbmGx2cMqoncw7sYf5X6rkpAsmEzxhclGtWyRS6NRmL4MSDsOm51t55dEPefWF\nMK++V8277ccxw/MOc0btZs7MELMWlnPyeeOpmD0t8eeDiAw7hb0Mua4u2PzcITY+8REbXwrxxrYq\n3mwdSwO7+UL1Tk6a1M7Js32c1FjPpMWT8YwZrb8CRHKsIDto1Waf38rKYMGX61jw5bpPnotGYetr\nY9n8lGXTi53c/YyfzfeP5GCoguneTcys+4gTxndxwoleZpxWy9QzxxGYNkF9ASKDpDZ7yQtthy1b\nX2zhzQ372PpGD1vf87F1Xw07u0bSYPYwvXwP0487zLTJUaadWMK0BXU0nHY83oYx+mtAZADUjCN5\nKRyG7Vu6eOfP+3j7lTbefTvOu7tLePdQPQfDlUwwHzClfB9TRrYxZUKUydP9TDq5kkkLRlE+cwKU\nljr9TxDJKwXZjCPFLxCAGXPKmDFn0lE/6+qCHVsaeO//+dn2RjvvvRdj7ZN+tv1HDTu7R1Jp25jo\n28rEqoNMHNXFhIY4E6b6mXBiJeNn11P9+QaornbgXyVSWBT24qiyMvj8KeV8/pTJR/0sHod9ewPs\neAV2vlHGzq3dbN4e5/H/CvDBfZXs7ByJ10Y43vMWx5ce5PjadhpGhWloMDRMDjBuegXjTqylesYY\nTH2dmorE1RwPe3XQyrF4PDBmnIcx40Zx+vJRR/3cWmg5ZNn1lmH3pjJ2/bWDXdvCPP++YfcrQfa0\nVbCnu45Y3DDOvMfY4EHGVHQwpraHsaOjjBnr4bgJQcZMKWPMzFpqpo3EjByhjmTJW+qgFelHezvs\nebeLD99sYe877Xy4o4e9u2Ps3edlb0uQve0V7OuppjseZDT7GO07xOjSNkZVdDG6NsyokXFGjfYw\nqiHAqAmljJxUwYipNQTGjYTKSv3FIMNObfYifais7O0zKOv3dT09sG/PWD7aWsb+99rYt6OLfbsj\n7Nwf55XNHvY9H+TjzlI+7qnkQKSacjoZyTZG+A8zsrSd+vIeRlSFqa+NU19vqB/lpX5MgPqGUuqO\nL6duYhUlY2qhpkZ/Pciwy5+w16m7OKykBCZM8TFhyghgRL+vjcehtbWGA7tLOLCtko93dHBwTw8H\n9kY4sD/Gu3sMh/7q5WB7kINdpRwKl3MwUoWPKHXspc5zmLpAO7XBbmrLQtRWRKipilNba6mp9VA7\n0kf1CD81x5VQM7aM6rHlVIypxNRUQ3m5/pqQAcufsBcpIB4P1NVBXV0Jnzu5BBid9j3WQkdHgJYD\nQVp2V9HyQTuH9nTT8lGIlv0RWg7G+WsrtH5oaGn30doV5HAoSGuojNZoOT02SBVtVLOTam8H1b4u\nqgPdVJeEqCyNUlUeparCUlVlqawyVNV4qaz1UVnnp3JEMHEbVUrFqDL8tRWJP3lKSvTF4RIKe5Fh\nYkwiXysrvYyfVAVUDej9kQi0tdVy+EAFbXs7Oby3i7b9PRzeH6LtUJS2lhiHW+LsbYf2/Ya2Tg/t\n3T7aegK0h4O0R0roiJbQHi/HT4RK2qngIyo8XVR4e6jw91DuD1MRiFAejFJeEqOiLE5ZqaW8HMor\nTOJW5aWs0ktZlY/y2gBl1X7KaoOU1QYprQniry7DlJcl5kcEAvoyyROOh/1nRuPofwqRY/L7ob4e\n6uv9ML0GyO5aBNZCd7eXzs4SOlprad/XReehEB0Heug4FKajJUzn4Rgdh2N0tsc51GHZ1QmdBwyd\n3R66Qh66Qj46wz66I366on46oyV0xwN0xhMT4MroopTWxL0JUeoJUeoNU+KNUOqLUOqPUOqPURqI\nURqMURKwlASTt5LEHxwlpYaSUkOw1ENJmYeSci/BMi8lFT6CZV6C5T6CFf5Pb5UBAhUBTEkw8SUT\nTN4X0YJ9xTEa57TT4M47E/cpNBpHpLBEIonJct3d0N0epaslRHdriO7DYbrbwnS3RenpiNLdHqW7\nM0ZPZ5yebkt3V5ye7kRHeSgE3SFDKGzoCXvoCXsJRT30RLz0RH2EYr03Pz1xP2HrJxT3EyFAgBAB\nwgQJEex9bMIETJSAJ0LAE03eYgS8yXtfjIA3jt8bJ+CL4/dZ/D5LwJ+49/stfl/iC/eTWwD8foPP\nb/AHErfex76AB3/Q4PN78Ac9+AIpt6D36PugF68/+bjEh9fvwRv0YXzeRGe+z5eYPJi8zKhG44iI\n4/z+RC5VVwPH+UjETPmwfHY8DuFwkHA4SChUSSgE4e4Y4Y4woY4Ike4ooc7ELdITI9wVJRKKE+qK\nEQnFCfckbpFQnEjYEg5ZImFLKAIdEUskDJEQRDogEjVEohCNmuRjQzRuiMYMkZiHaMxDJHmLWUM0\n7iEa9xCJe4lZD9G4l4hNPrbeT24x6yWKlzhePMTwEcVLjD/dtpkl35+T9b5R2ItI0fB4+KQZ6FNe\noDR5KxzWQizmJRbzEo1CMJh90IPCXkQkLxmTaL3x+RLdD4PlGfyvyC211YuIDF7ehz1okI6IyGA5\n3oyjhdBERDJT1EMv4/FEm1U87kCBIiJ5KJuhlwXRjCMiIoOjsBcRcQGFvYiICyjsRURcQGEvIuIC\
	nCnsRERdQ2IuIuIAmVYmIFAhNqhIRcRFNqhIRkT4p7EVEXEBhLyLiAgp7EREXUNiLiLiAwl5ExAXy\nJ+x1/UERkZzJn7AXEZGcya+w18VmRURyIr/CXkREciLvw15N+SIig1cQC6GpdUdEpFgWQjv1VLjr\nrsR9ilgMAoHEvYiIaCE0ERE5BoW9iIgLKOxFRFxAYS8i4gIKexERF1DYi4i4gMJeRMQFFPYiIi6g\nsBcRcQGFvYiICyjsRURcQGEvIuICCnsRERdQ2IuIuEBO17M3xlQDk4FaYLu19v1cfp6IiPQtozN7\nY8yqI7ZXGGMWG2NWpnnrPGAJ0ALUZFeiiIgMVtqwTwb6ipTt2YC11q5Pbs/q5+2vAlOA3wDbB1eq\niIhkK23YW2vv5rNBfSnQmny8ncSZe+/Z/kUp99XAJdbaq5Kv+duhLV2OlO3lyqRv2p9DR/vSeZl2\n0KZe/qoGOJSyXQ9grX3YWvtIyv1hYJsx5ksk2u3XDUnFckw6oIaW9ufQ0b50Xk47aK21G3L5+0VE\nJDOZntmnXpW8BahLPq4BDg5pRSIiMuSMtTb9i4xZa609O/l4NjDXWrvaGHMdsM5a+0ZWH25M+g8X\nEZGjWGtN+ld9Km0zjjFmBTDXGHOFtXa1tfZ1Y8xcY8xioCXboM+mWBERyU5GZ/YiIlLYtFyCCEdP\nHDziZ5lOIpSkNPtzVfJe+3MYDVvYpztgdEANTAb7UwdUho6cOHjEzwYyiVDof38mXWmMeRfYNkwl\nFTRjzMrkrc8v0Eyzc1jCPt0BowNqYDLcXzqgMtTHxMFUfU4ilGNLsz8BrrDWTtPQ7PSSfaPrkvt0\ncnLeUurPM87O4TqzT3fA6IAamEz2lw6oodHnJEIZlLrkmeh1ThdSACbz6fG9PbmdKuPsHK6wT3fA\n6IAamEz2lw4oyUvJUX3rgfojz1Tls6y1d1trVyc355BYbyxVxtmpDtoipQNqyGgS4RBKtj1flNw8\nyNFnqtKHZHPNa4MZ6j5cYZ/ugNEBNTD97i8dUFn5zJyP5EJ+AGv4dP9NBp4ZzqIK2LH25zY+3Yf1\nHH2mKn1bbK39hz6ezzg7hyvs+zxgdEBlLd3+1AE1AKkTB1OefgbAWvt68jWDnkToFmn25wZgafI1\nB7Q/0zPGrLTW3pF8vDh5P+DsHLZJVcn/8DuASb1tUMaYV6y184/1czm2DPZn79C3Sb3/o4hIYUmG\n+xoSZ/C1wMXW2g3ZZKdm0IqIuIA6aEVEXEBhLyLiAgp7EREXUNiLiLiAwl5ExAUU9iIiLpDTC46L\nFLrk5JXexaWWWmuvcrIekWzpzF6kf5eQWEL2YaD1iFmhIgVDZ/Yi/UiuI95rMrDWqVpEBkMzaEUy\nkJy2rqU8pGAp7EXSSC4vi7X2dWPM7N7F0UQKiZpxRJKSF3p5jcSCU5OBXwNTgIeAbcaYOuB65yoU\nyZ7O7EUAY8yvgLXW2keMMZOSj6c5XZfIUNFoHHG9ZLgvttY+knxqDrqmghQZhb1IYhx9argvBdY5\nVItITijsRRIXbG5N2b4EndlLkVEHrbietfZhY8yS5IXZpwAHrbVtTtclMpQU9iKAtfY7AMaYWnRW\nL0VIzTginzUftddLEVLYiyQZY1YCK0kEvkhR0Th7EREX0Jm9iIgLKOxFRFxAYS8i4gIKexERF1DY\ni4i4wP8HRmIaCDiolRAAAAAASUVORK5CYII=\n",
	"text/plain": [
	"<matplotlib.figure.Figure at 0x7f57d3f8b940>"
	]
	},
	"metadata": {},
	"output_type": "display_data"
	}
	],
	"source": [
	"flavio.plots.q2_plot_th_diff('dBR/dq2(B0->K*ee)', 0, 2, color='red')\n",
	"flavio.plots.q2_plot_th_diff('dBR/dq2(B0->K*mumu)', 0, 2, color='blue')\n",
	"plt.yscale('log')\n",
	"plt.xlabel(r'$q^2$')"
	]
	},
	{
	"cell_type": "markdown",
	"metadata": {},
	"source": [
	"## Experiment vs. theory for binned observables"
	]
	},
	{
	"cell_type": "code",
	"execution_count": 27,
	"metadata": {
	"collapsed": true
	},
	"outputs": [],
	"source": [
	"th_bins = [(0.1, 2), (2.0, 4.0), (4.0, 6.0), (15.0, 19.0)]"
	]
	},
	{
	"cell_type": "code",
	"execution_count": 32,
	"metadata": {
	"collapsed": false
	},
	"outputs": [
	{
	"data": {
	"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXIAAAEECAYAAAA1X7/VAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAD5hJREFUeJzt3cFyHMd9x/HfP+VzCANH62CsS3dCK6RAXbAlMC8QWOYL\nEMrRVS7TkU6CTmEqTFWOsYkHcFimHiAkUaOLyDLFIu4RQR1snwIC5Av8c9gGMFotMLvY3pn573w/\nVVua2ZmdaYx3f2529/SYuwsAENffNV0AAMBsCHIACI4gB4DgCHIACO4nTRdgFma2LemWpDeSvnP3\new0XCQBqZ20YtWJmd939s9L6tqQTST13v3/J537u7t+b2ceSvnX3dzUUFwBapfGmFTPbkbRdWl+T\n5O7+JK1fv+iz7v796ccIcQBd1XjTirvfN7Nflt66Jel/0vKhpJuSDlIt3SVZ+u9jd3+X3n9VZ5kB\noE3mHuSphn3i7q/T+pak5yM1aCstL2nY5n1qRZLc/eEFp1iV9ChfiQEglrk3rbj7S0lbZraaQtxz\nNoO4+z2aVQB02URBbmZ3L9m2bWZbqa17LHffk/QvklbdfX/cLqXlY0nLaXlJ0tEkZQSArqoM8tHO\nyJFtE3VMppr4g7S8Om6X0vIDSb203JP0uKqMANBllUGehv8dXrD5lobDBKXzjskfKIX9fqqZb5nZ\n35e2b0vqm9ntdL6X6f0tScfufjDF3wMAnTNrZ+fYjsmy02Aure+NrD+U9PCyfQAAF2t8HDkAYDaz\nBjkdkwDQsEmbVuwHK2bX3P2thh2TfUn7GnZM/mg8t5k1PwcAAATk7la912SjVn7QGZk8TieZqGPS\n3Xllen3xxReNl2FRXlxLrmebX9OorJH7+M7I9dIyHZMA0CA6OwEguMYnzcJ0BoNB00VYGHVey6Io\nVBTF2fLpuQeDwcL8b7oof0dEc5+P3Mx83ucAIjGzqdtA0T3pe5KnsxMA0G4EOQAER5ADQHAEOQAE\nR5ADQHAEOQAER5ADQHAEOQAER5ADQHAEOQAER5ADQHAEOQAER5ADQHAEOQAER5ADQHAEOQAER5AD\nQHC1Pert9bd/retUP7L64c8aOzcAzBs1cgAIjiAHgOAIcgAIjiAHgOAIcgAIjiAHgOAIcgAIjiAH\ngOAIcgAIjiAHgOAIcgAIjiAHZlAUxUKcA7ER5MAMCHK0AUEOAMHVNo0tsIiKopCZTf25aT6zubk5\n9fHRLZU1cjPbNrMtM9up2H47f/GAdhsMBnL3qV6Sptp/MBg0+0ei9S4NcjNbk+Tu/iStXx+z/TBt\nfz26HQAwf1U18luSTtLyoaSbY/b5t/Tfnrsf5CoYEEEdtWVq5KhSFeRLkt6U1lfKG939paRDM3sj\n6Shz2YDWI8jRBjONWjGza5K+k3Rb0n0z+3mGMgEAplA1auVY0nJaXtKPa92fSvq9u78zsxNJn0j6\n99GD7O7u6vhv7yRJG/0b2uh/NFOhAWDRFEVx5XsG7LQXfezGYWdm3933zOyOpEfufmBm19z9rZn9\n1t3vlfa/7e57I8dwd9frb/96pQLmsPrhzxo7NzDKzHTZ7w6Qzr4nE41TvbRG7u4vzaxvZluSjkud\nmY8lrbv7vRTwryQtj4Y4AGD+Lq2RZzkBNXLgB6iRYxLT1Mi5RR8AgiPIASA4ghwAgiPIASA4Zj8E\n0Kjy+OmiKM7uZB0MBtzVOiFGrQA1Y9TKxbg25xi1AgAdQpADQHAEOQAER5ADQHAEOQAER5ADQHCt\nHkf+7MU3evbiaVp+qo3+DUnMaQ4AZWHGkffW39Ph879c6bOMI0ebMFb6Ylybc4wjB4AOIcgBILhW\nt5E3jTkgAERAG/mEaLtDLnyXLsa1OUcbOQB0CEEOAMHV1rSS4TiN/pOr6fNjcfBduhjX5hxNKwDQ\nIYxaAZBFrofHNPkQmlzqvgmRGjkABEeNHKhB+Z6Ezc1N7e7uSuKeBORBkAM1ILAxTzStAEBwBDkA\nBEeQA8ji2Ytvmi6CpPaUo04EOYAsTh8C07S2lKNOBDkABMeoFQBZPHvxVL3192Y+zqzH+IcPNmYu\nQzQEOYAsNvo39Mff/2mmY8wyXfWp//zDf8z0+YhoWgGA4AhyAFls9G80XQRJ7SlHnQhyAFls9D9q\nugiS2lOOOlW2kZvZtqQTST13vz9m+5qknqTlcdsBAPN1aZCnkHZ3f2JmPTO77u4HI7t97u6/MrPf\nXrAdSRum56x7ek0A81fVtHJLw9q4JB1KulnemGrrf5Ykd79HiANA/aqCfEnSm9L6ysj2dUkrZrZm\nZneylgwAMJEcnZ1H7v5SOquhj3U6F3MTmjw3AMxbVWfnsaTltLwk6Whk+5GGTS7SsAnmQ0kPRw+y\nu7uroijO5mSue17m03MDQFuVHz4yraogfyCpL2lfw5EpjyTJzK65+1tJf5J0WgtfkvR83EF2d3fP\nXgCAHxut5H755ZcTf/bSppVSk8mWpONSZ+bjtP21pJPUpLLs7l9ddKyiKGRmV36lclzpRdMKgEVW\n2Ubu7nvu/sTd90rvrY9sf+jun192nMFgIHe/8iud60ovmlUALDLu7ASA4GoL8iZrxdTIASyyTgR5\nW3TxEVQA5q8TTStt6ezs4iOoAMxfJ4IcABZZJ54QdDr0cVazHqOLj6ACMH+dCPLBYDBz84qZnQ2D\nvKpff/qbmT4PtFmumTWZoXN6nQjy47+9yzKFbBumoQWAUZ1oI2/Lo5/aUg4Ai6UjQd6ORz+1pRwA\nFksnghwAFhlBDgDBEeQAEBxBDgDBdWL4IYD2Kj8ZZ3Nz8+wBNE08TSwqghxAowjs2dG0AgDBEeQA\nEBxBDgDBEeQAEBxBDgDBEeQAEBxBDgDBMY68RkyYD2AeqJEDQHDUyAGEV77NvyiKsztFu3LXqM36\nHMrKE5h5jnPM8szMHI9o662/p8Pnf5npGDStAPOX4/m6bZD+jome+E7TCgAER5ADQHCtbiNneksA\nqNbqICewAaBaq4M8gmcvvtGzF0/T8lNt9G9Ikjb6N7TR/6jJogHoiDCjVmZR16iVqn0YtQLMXxdH\nrXQiyHOY5MuxKF+gcbo+ThdxLMrvkCCfg64HeVlX/k7EtCjfz6zjyM1s28y2zGynYr87kxYQAJDP\npUFuZmuS3N2fpPXrF+y3Jelm/uIBAKpU1chvSTpJy4cirAGgdaqCfEnSm9L6yugOZraWauwTteUA\nAPLKMY78pxmO0UrcWQpMrjyaqevqvhZVQX4saTktL0k6Km9MtfH9tBq/m3gEgQ1MjiA/17YgfyCp\nL2lfUk/SI0kys2vu/lZSz8xWNWxyWTGz6+5+MHqQ05qsRDgCwDjlFoBpVY4jN7Pbkl5LWnX3vfTe\nc3dfL+2zI+l3kj4ZDfJFGUc+iUUZv1qlK38npjMYDPT11183XYxW2NzcvHIon8o6jtzd99z9yWmI\np/fWR/a57+7vj6uNA+iGwWAgd2/8JanxMtTd6sB85AAQXJggn/WfKRjiOmJe6Ps6R438AgRQHlxH\nzAtBfo4gBwBMJcyDJYqikFn7bx5texk3NzebLgKAzMLUyNvSI9723vK29aYDmL8wQQ4AGC9MkFOT\nzIPrCCwenhCU0bzueGzbY9a4sxNttijfTx711pA6vkCzniPHZD45/k4mWMK8dDHIwzStII+2jCNv\nSzmARUCQA0BwYcaRI49c4/FnPQbj2YF8qJF3TI7x+NLs4+VpHwfyIcgBIDiaVjqmLTXhtpQDi6Hr\nz9dl+GFGEYYfLkoZgEXH8EMA6BCCHACCI8gBIDiCHACCI8gBIDhGrcyo7pkJ2zBipA1lABYdsx8u\nsDaEaBvKACw6hh8CQIcQ5AAQHEEOAMER5DVqw8MU2lAGAHkR5DVqQ4i2oQwA8iLIASA4prGtURue\nzsOTeYDFQ428Rm14Ok8X5mYGuoYgB4DgCPIataE23IYyAMiLW/SDacPt8W0oA7DouEUfADqkskZu\
	nZtuSTiT13P3+mO07afEX7v7ZmO3UyDNqqjZc9yyPQNdlm/3QzNYkrbr7Vymwn7v7QWn7lqRX7v69\nmT2Q9F/uvj9yDII8I5o1gG7I2bRyS8PauCQdSro5sr1Xeu8wrQMAalR1Q9CSpDel9ZXyxpGmlg8k\n/TFTuQAAE8rS2ZmaYF6Um10AAPWoqpEfS1pOy0uSji7Yb8vdP7/oILu7u2fLdI4BaJs2dOaXyzCt\nSTo7++6+Z2Z3JD1y9wMzu+bub9M+O6dNLGa25e5PRo5BZ2dGdHYC89WW31i2zk53f5kOuCXpuNR0\n8rj0/l0z+87MjiQ1/9cDQMdwZ2cwbaktAIuqLb8x7uwEgA4hyAEgOIIcAIIjyAEgOIIcAIIjyAEg\nOIIcAIIjyAEgOG4ICqAN80AAXRHxhiCCHABKIgY5TSsAEBxBDgDBEeQAFspV5/SOXAaCHMBCIcgB\nAOFUPeoNAEIpikJmEw32uNCsn9/c3Jzp89OiRg5goQwGA7n7lV+SZvq8u9d+fwdBDgDBEeQAFkob\n7nauuwzc2QkAJdzZCQCoHUEOAMER5AAQHEEOAMER5AAQHEEOAMER5AAQHEEOAMER5AAQHEEOAMER\n5AAQHEEOAMER5AAQHEEOAMExjS2AziuK4uyByUVRnM0nPhgMGpvffJppbAlyAGihaYK88uHLZrYt\n6URSz93vT7sdADBfl7aRm9maJHf3J2n9+jTbkd/pP/8wO65lXlzP5lR1dt7SsLYtSYeSbk65HZnx\nY8mHa5kX17M5VUG+JOlNaX1lyu0AgDlj+CEABHfpqBUz+1dJj9x9P3Vqrrr7vUm3p30YsgIAV5Br\n1MoDSX1J+5J6kh5Jkpldc/e3F22/SkEAAFdzadOKu7+UJDPbknTs7gdp0+OK7UArmdndkfVtM9sy\ns52myhTZmOt5N/2X61mjud0QxPjyvMzsrrt/ZmY7XM+rSeHyO3d/P62vadgc+FXa9pzKyORGr2d6\n742kI0n/7O77jRUuoNL/+f3C3T9L702Uo3Pp7GR8+Vx8amb/K+lV0wWJKv0QDktvMXx2BmOupyTd\ndvf3CfHppFaNR+ma9szs42lydF6jVviB5McPJI9ynw3DZ/NbTk1Vd5ouSDA9nefkYVqfOEfnFeT8\nQPLjB4LWc/e9VINcMbOPmy5PFO5+39330uoHkr7VFDnKOPIg+IFkU+4UOpa0nJaXNGzbxRWZ2Y6Z\n/VNaPdKwVokppOaUF9P21cwryPmBZMQPJKty08oDnV/LntJoLEylfD1f6fwarmhYq8R0ttz987Q8\ncY7OK8j5geTFDySDNAKgb2a3JYbPzmrM9dyX9I/p/f/jek4njUi7l5a3JP23JszReQ4/vC3ptYbD\nu/aq9sfl0o9DGnP3LIDYUnA/0LAW/lNJn6Q75ifK0bk/WAIAMF90dgJAcAQ5AARHkANAcAQ5AARH\nkANAcAQ5AARHkANAcAQ5AAT3/9Jq0wQuTLOyAAAAAElFTkSuQmCC\n",
	"text/plain": [
	"<matplotlib.figure.Figure at 0x7f57d3c200f0>"
	]
	},
	"metadata": {},
	"output_type": "display_data"
	}
	],
	"source": [
	"flavio.plots.q2_plot_exp('<dBR/dq2>(B+->K*mumu)', c='k')\n",
	"flavio.plots.q2_plot_th_bin('<dBR/dq2>(B+->K*mumu)', bin_list=th_bins)"
	]
	},
	{
	"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.4.3"
	}
	},
	"nbformat": 4,
	"nbformat_minor": 0
	}