1-WH-Morphing-Setup.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# 1. Morphing Setup for $W^{\\pm}(\\ell^{\\pm}\\nu_{\\ell})H(b\\bar{b})$"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Here we set up the morphing basis for $W^{\\pm}(\\ell^{\\pm}\\nu_{\\ell})H(b\\bar{b})$ with all the relevant operators both in the gauge boson couplings and the gauge-fermion couplings:\n",
    "\n",
    "$C_{H\\square}$, $C_{HD}$, $C_{HW}$, and $C^{(3)}_{Hq}$\n",
    "\n",
    "We omit the operator $C_{dH}$ which appears in the bottom-yukawa coupling, modifying the $H \\to b\\bar{b}$ decay, as it is poorly constrained by this process and well measured elsewhere.\n",
    "\n",
    "As can be seen analytically, the operators $C_{HD}$ and $C_{H\\square}$ only change $W^{\\pm}H$ production by an overall rescaling of the couplings, and always enter in the combination \n",
    "\n",
    "$$C_{H\\square} - \\frac{1}{4} C_{HD}.$$\n",
    "\n",
    "As an extra cross check of our procedure, we'll consider all four operators, then rotate the resulting $4 \\times 4$ Fisher Information to explicitly remove the inherent flat direction, and do most of our analysis on the resulting $3 \\times 3$ matrix.\n",
    "\n",
    "This setup is used for both $W^+H$ and $W^-H$ production, though they must be generated and analyzed seprately (the same is true for the $W$ decays)."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Setup"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "MadMiner version: 0.3.0\n"
     ]
    }
   ],
   "source": [
    "madminer_src_path = '/home/shomiller/madminer'\n",
    "import sys\n",
    "import os\n",
    "sys.path.append(madminer_src_path)\n",
    "\n",
    "from __future__ import absolute_import, division, print_function, unicode_literals\n",
    "import logging\n",
    "import numpy as np\n",
    "import matplotlib\n",
    "from matplotlib import pyplot as plt\n",
    "%matplotlib inline\n",
    "\n",
    "from madminer.core import MadMiner\n",
    "\n",
    "import madminer.__version__\n",
    "print( 'MadMiner version: {}'.format(madminer.__version__) )"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# MadMiner output\n",
    "logging.basicConfig(\n",
    "    format='%(asctime)-5.5s %(name)-20.20s %(levelname)-7.7s %(message)s',\n",
    "    datefmt='%H:%M',\n",
    "    level=logging.INFO\n",
    ")\n",
    "\n",
    "# Output of all other modules (e.g. matplotlib)\n",
    "for key in logging.Logger.manager.loggerDict:\n",
    "    if \"madminer\" not in key:\n",
    "        logging.getLogger(key).setLevel(logging.WARNING)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 1. Setup Morphing Basis"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Here we add the parameters to the morphing setup. Note that we rescale the $C_{HW}$ and $C_{Hq}^{(3)}$ coefficients, so that all the Fisher Information entries will be roughly the same size. This helps with numerical aspects later in the analysis."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "14:30 madminer.core        INFO    Added parameter cHbox (LHA: frblock 24, maximal power in squared ME: (2,), range: (-10.0, 10.0))\n",
      "14:30 madminer.core        INFO    Added parameter cHDD (LHA: frblock 25, maximal power in squared ME: (2,), range: (-10.0, 10.0))\n",
      "14:30 madminer.core        INFO    Added parameter cHW (LHA: frblock 28, maximal power in squared ME: (2,), range: (-10.0, 10.0))\n",
      "14:30 madminer.core        INFO    Added parameter cHq3 (LHA: frblock 49, maximal power in squared ME: (2,), range: (-10.0, 10.0))\n"
     ]
    }
   ],
   "source": [
    "miner = MadMiner()\n",
    "\n",
    "miner.add_parameter(\n",
    "    lha_block='frblock',\n",
    "    lha_id=24,\n",
    "    parameter_name='cHbox',\n",
    "    morphing_max_power=2,\n",
    "    parameter_range=(-10.,10.)\n",
    ")\n",
    "miner.add_parameter(\n",
    "    lha_block='frblock',\n",
    "    lha_id=25,\n",
    "    parameter_name='cHDD',\n",
    "    morphing_max_power=2,\n",
    "    parameter_range=(-10.,10.)\n",
    ")\n",
    "miner.add_parameter(\n",
    "    lha_block='frblock',\n",
    "    lha_id=28,\n",
    "    parameter_name='cHW',\n",
    "    morphing_max_power=2,\n",
    "    param_card_transform='0.1*theta',\n",
    "    parameter_range=(-10.,10.)\n",
    ")\n",
    "miner.add_parameter(\n",
    "    lha_block='frblock',\n",
    "    lha_id=49,\n",
    "    parameter_name='cHq3',\n",
    "    morphing_max_power=2,\n",
    "    param_card_transform='0.01*theta',\n",
    "    parameter_range=(-10.,10.)\n",
    ")\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "For now, we'll only add one benchmark manually at the SM, and let the morphing take care of the rest:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "14:30 madminer.core        INFO    Added benchmark sm: cHbox = 0.00e+00, cHq3 = 0.00e+00, cHDD = 0.00e+00, cHW = 0.00e+00)\n"
     ]
    }
   ],
   "source": [
    "miner.add_benchmark(\n",
    "    {'cHbox':0., 'cHDD':0., 'cHW':0., 'cHq3':0.},\n",
    "    'sm'\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now we set the rest of the benchmarks automatically. To try and reduce any possibility of morphing error, we'll run a number of trials and tests to ensure a decent spread of the weights."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "14:30 madminer.core        INFO    Optimizing basis for morphing\n",
      "14:30 madminer.core        INFO    Set up morphing with 4 parameters, 15 morphing components, 1 predefined basis points, and 14 new basis points\n"
     ]
    }
   ],
   "source": [
    "miner.set_morphing(\n",
    "    include_existing_benchmarks=True,\n",
    "    max_overall_power=2,\n",
    "    n_trials=1000,\n",
    "    n_test_thetas=1000\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Setup Systematics"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Here we add scale and PDF systematics.\n",
    "\n",
    "For scale uncertainties, we vary both $\\mu_R$ and $\\mu_F$ independently by a factor of $1/2$ and $2$.\n",
    "\n",
    "For the PDF variations, we run over all the eigenvectors of the `PDF4LHC15_nlo_30` set."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "miner.set_systematics(\n",
    "    scale_variation=(0.5,1,2),\n",
    "    scales='independent',\n",
    "    pdf_variation='90900'\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Save Setup\n",
    "\n",
    "And finally, save the setup to a `.h5` file. Since the basis is well optimized, we should just use the same setup for all our later $WH$ runs."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "14:30 madminer.core        INFO    Saving setup (including morphing) to data/wh_smeft_setup.h5\n"
     ]
    }
   ],
   "source": [
    "miner.save('data/wh_smeft_setup.h5')"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.13"
  },
  "toc": {
   "base_numbering": 1,
   "nav_menu": {},
   "number_sections": true,
   "sideBar": true,
   "skip_h1_title": false,
   "title_cell": "Table of Contents",
   "title_sidebar": "Contents",
   "toc_cell": false,
   "toc_position": {},
   "toc_section_display": true,
   "toc_window_display": false
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}