Skip to content

Instantly share code, notes, and snippets.

@restrepo

restrepo/4VL.ipynb

Created October 2, 2015 05:20
Show Gist options
  • Save restrepo/419a25f93a4c103690d9 to your computer and use it in GitHub Desktop.
Save restrepo/419a25f93a4c103690d9 to your computer and use it in GitHub Desktop.
Download ZIP
Display the source blob
Display the rendered blob
Raw
4VL.ipynb
{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"From http://arxiv.org/pdf/1404.4398v2.pdf eq.(65)\n",
"\\begin{align}\n",
"\\sigma(gg\\to h)=\\frac{8\\pi^2}{m_h^3}\\Gamma(h\\to gg)\n",
"\\end{align}\n",
"Details are referenced to [hep-ph/9504378]"
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"import pandas as pd\n",
"import numpy as np\n",
"import os, sys, inspect\n",
"import commands"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"from hep import *\n",
"class fourVL(THDM):\n",
" MESSAGE=THDM.MESSAGE+\" with vectorlike quarks\"\n",
" def __init__(self):\n",
" super(fourVL, self).__init__() #Get THDM.__init__\n",
" self.MassQ=0.\n",
" self.MassU=0.\n",
" self.MassD=0.\n",
" self.fd=0.\n",
" self.fu=0.\n",
" self.fdt=0.\n",
" self.fut=0.\n",
" def init_lha_blocks(self,LHA):\n",
" super(fourVL,self).init_lha_blocks(LHA) #Get THDM init blocks\n",
" if self.LHA.blocks.has_key('LAMQ'):\n",
" for b in ['LAMQ','LAMD','LAMU']:\n",
" for i in range(1,4):\n",
" self.LHA.blocks[b][i]='%.8E # %s' %(0,b)\n",
" else:\n",
" print \"Use hep.buildSLHA to include proper blocks\" \n",
" \n",
" return self.LHA"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"if 1==1:\n",
" xd=pd.DataFrame()\n",
" h=hep()\n",
" h.MODEL='4VL'\n",
" t=fourVL()\n",
" LHA=h.buildSLHA(['modsel','sminputs','minpar','sphenoinput','EPSUIN','EPSDIN','EPSEIN',\\\n",
" 'LAMQ','LAMD','LAMU'])\n",
" LHA=t.init_lha_blocks(LHA)\n",
" \n",
" t.m_h0=125.;t.m_H0=2000;t.m_A0=2100.26;t.m_Hp=2090.99\n",
" \n",
" t.tanb=-34.2325 #-???\n",
" t.sab=-0.98 #np.sin(np.arctan(t.tanb)-np.arccos(cosa))\n",
" t.lambdas[6]=0.;t.lambdas[7]=0.;t.m12_2=0\n",
" \n",
" t.MassQ=800.;t.MassU=1100.;t.MassD=1100.\n",
" t.fd=0.;t.fu=2.3;t.fdt=0.;t.fut=1.6\n",
" \n",
" #Intialize output table\n",
" ins=pd.Series({})\n",
" #for t.m_h0 in np.linspace(100,2000.):\n",
" for t.m_H0 in [2000]: #np.linspace(t.m_h0+10,2000.):\n",
" #convert into general THDM basis\n",
" perturbativity=t.phys_to_gen()\n",
" #generate LHA input file\n",
" for i in range(1,6):\n",
" LHA.blocks['MINPAR'][i]='%.8E # Lambda%dInput' %(t.lambdas[i],i) \n",
" LHA.blocks['MINPAR'][6]='%.8E # TanBeta' %t.tanb\n",
" LHA.blocks['MINPAR'][7]='%.8E # MQInput' %t.MassQ \n",
" LHA.blocks['MINPAR'][8]='%.8E # MUInput' %t.MassU \n",
" LHA.blocks['MINPAR'][9]='%.8E # MDInput' %t.MassD \n",
" LHA.blocks['MINPAR'][10]='%.8E # fdInput' %t.fd \n",
" LHA.blocks['MINPAR'][11]='%.8E # fuInput' %t.fu \n",
" LHA.blocks['MINPAR'][12]='%.8E # fdtInput' %t.fdt \n",
" LHA.blocks['MINPAR'][13]='%.8E # futInput' %t.fut\n",
" \n",
" SPC=h.runSPheno(LHA) #Set paths and commands in h.config\n",
" #reset branchings\n",
" for p in [h.pdg.h0,h.pdg.H0]:\n",
" h.Gamma[p]=0\n",
" h.Br[p]={}\n",
" h.Br[p][h.pdg.g,h.pdg.g]=0\n",
"\n",
" decays=h.branchings(SPC.decays) # -> h.Gamma[pdg], h.Br[pdg][pdg1,pdg2,...] \n",
" \n",
" od=pd.Series({'$m_{h^0}$':SPC.blocks['MASS'][h.pdg.h0],'$m_{H^0}$':SPC.blocks['MASS'][h.pdg.H0],\\\n",
" r'$\\Gamma(h^0\\to gg)$':h.Br[h.pdg.h0][h.pdg.g,h.pdg.g]*h.Gamma[h.pdg.h0],\\\n",
" r'$\\Gamma(H^0\\to gg)$':h.Br[h.pdg.H0][h.pdg.g,h.pdg.g]*h.Gamma[h.pdg.H0],\\\n",
" r'$\\sigma(gg\\to H^0)$':\\\n",
" 8*np.pi**2*h.Br[h.pdg.H0][h.pdg.g,h.pdg.g]*h.Gamma[h.pdg.H0]/SPC.blocks['MASS'][h.pdg.H0]**3,\\\n",
" 'pert.':perturbativity})\n",
" xs=ins.append(od)\n",
" xd=xd.append(xs.to_dict(),ignore_index=True)\n",
"\n",
" xd.to_csv('thdm.csv',index=False)"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/html": [
"<div style=\"max-height:1000px;max-width:1500px;overflow:auto;\">\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <th></th>\n",
" <th>$\\Gamma(H^0\\to gg)$</th>\n",
" <th>$\\Gamma(h^0\\to gg)$</th>\n",
" <th>$\\sigma(gg\\to H^0)$</th>\n",
" <th>$m_{H^0}$</th>\n",
" <th>$m_{h^0}$</th>\n",
" <th>pert.</th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>0</th>\n",
" <td> 15278.746823</td>\n",
" <td> 0.025717</td>\n",
" <td> 0.000151</td>\n",
" <td> 2000</td>\n",
" <td> 125.000002</td>\n",
" <td> False</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" $\\Gamma(H^0\\to gg)$ $\\Gamma(h^0\\to gg)$ $\\sigma(gg\\to H^0)$ $m_{H^0}$ \\\n",
"0 15278.746823 0.025717 0.000151 2000 \n",
"\n",
" $m_{h^0}$ pert. \n",
"0 125.000002 False "
]
},
"execution_count": 4,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"xd"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"BLOCK MODSEL\r\n",
" 1 0 # 1/0: High/low scale input\r\n",
" 2 1 # Boundary Condition\r\n",
" 6 1 # Generation Mixing\r\n",
" 12 1000. # Renormalization scale\r\n",
"\r\n",
"BLOCK SMINPUTS\r\n",
" 2 1.16637000E-05 # G_F,Fermi constant\r\n",
" 3 1.18700000E-01 # alpha_s(MZ) SM MSbar\r\n",
" 4 9.11887000E+01 # Z-boson pole mass\r\n",
" 5 4.18000000E+00 # m_b(mb) SM MSbar\r\n",
" 6 1.73500000E+02 # m_top(pole)\r\n",
" 7 1.77669000E+00 # m_tau(pole)\r\n",
"\r\n",
"BLOCK MINPAR\r\n",
" 1 3.99759756E+04 # Lambda1Input\r\n",
" 2 1.15198706E+00 # Lambda2Input\r\n",
" 3 5.55467608E+02 # Lambda3Input\r\n",
" 4 -7.60465292E+01 # Lambda4Input\r\n",
" 5 -7.74101866E+01 # Lambda5Input\r\n",
" 6 -3.42325000E+01 # TanBeta\r\n",
" 7 8.00000000E+02 # MQInput\r\n",
" 8 1.10000000E+03 # MUInput\r\n",
" 9 1.10000000E+03 # MDInput\r\n",
" 10 0.00000000E+00 # fdInput\r\n",
" 11 2.30000000E+00 # fuInput\r\n",
" 12 0.00000000E+00 # fdtInput\r\n",
" 13 1.60000000E+00 # futInput\r\n",
"\r\n",
"BLOCK SPHENOINPUT\r\n",
" 1 -1 # error level\r\n",
" 2 0 # SPA conventions\r\n",
" 11 1 # calculate branching ratios\r\n",
" 13 1 # 3-Body decays: none (0), fermion (1), scalar (2), both (3)\r\n",
" 14 0 # Run couplings to scale of decaying particle\r\n",
" 12 1.000E-15 # write only branching ratios larger than this value\r\n",
" 15 1.000E-40 # write only decay if width larger than this value\r\n",
" 31 -1 # fixed GUT scale (-1: dynamical GUT scale)\r\n",
" 32 0 # Strict unification\r\n",
" 34 1.000E-04 # Precision of mass calculation\r\n",
" 35 40 # Maximal number of iterations\r\n",
" 37 1 # Set Yukawa scheme\r\n",
" 38 2 # 1- or 2-Loop RGEs\r\n",
" 50 1 # Majorana phases: use only positive masses\r\n",
" 51 0 # Write Output in CKM basis\r\n",
" 52 1 # Write spectrum in case of tachyonic states\r\n",
" 55 0 # Calculate one loop masses\r\n",
" 57 0 # Calculate low energy constraints\r\n",
" 65 0 # Solution tadpole equation\r\n",
" 75 1 # Write WHIZARD files\r\n",
" 76 1 # Write HiggsBounds file\r\n",
" 86 0. # Maximal width to be counted as invisible in Higgs decays; -1: only LSP\r\n",
" 510 0. # Write tree level values for tadpole solutions\r\n",
" 515 0 # Write parameter values at GUT scale\r\n",
" 525 0. # Write loop contributions to diphoton decay of Higgs\r\n",
"\r\n",
"BLOCK EPSUIN\r\n",
" 1 1 0.00000000E+00 #EPSUIN(1,1)\r\n",
" 1 2 0.00000000E+00 #EPSUIN(1,2)\r\n",
" 1 3 0.00000000E+00 #EPSUIN(1,3)\r\n",
" 2 1 0.00000000E+00 #EPSUIN(2,1)\r\n",
" 2 2 0.00000000E+00 #EPSUIN(2,2)\r\n",
" 2 3 0.00000000E+00 #EPSUIN(2,3)\r\n",
" 3 1 0.00000000E+00 #EPSUIN(3,1)\r\n",
" 3 2 0.00000000E+00 #EPSUIN(3,2)\r\n",
" 3 3 0.00000000E+00 #EPSUIN(3,3)\r\n",
"\r\n",
"BLOCK EPSDIN\r\n",
" 1 1 0.00000000E+00 #EPSDIN(1,1)\r\n",
" 1 2 0.00000000E+00 #EPSDIN(1,2)\r\n",
" 1 3 0.00000000E+00 #EPSDIN(1,3)\r\n",
" 2 1 0.00000000E+00 #EPSDIN(2,1)\r\n",
" 2 2 0.00000000E+00 #EPSDIN(2,2)\r\n",
" 2 3 0.00000000E+00 #EPSDIN(2,3)\r\n",
" 3 1 0.00000000E+00 #EPSDIN(3,1)\r\n",
" 3 2 0.00000000E+00 #EPSDIN(3,2)\r\n",
" 3 3 0.00000000E+00 #EPSDIN(3,3)\r\n",
"\r\n",
"BLOCK EPSEIN\r\n",
" 1 1 0.00000000E+00 #EPSEIN(1,1)\r\n",
" 1 2 0.00000000E+00 #EPSEIN(1,2)\r\n",
" 1 3 0.00000000E+00 #EPSEIN(1,3)\r\n",
" 2 1 0.00000000E+00 #EPSEIN(2,1)\r\n",
" 2 2 0.00000000E+00 #EPSEIN(2,2)\r\n",
" 2 3 0.00000000E+00 #EPSEIN(2,3)\r\n",
" 3 1 0.00000000E+00 #EPSEIN(3,1)\r\n",
" 3 2 0.00000000E+00 #EPSEIN(3,2)\r\n",
" 3 3 0.00000000E+00 #EPSEIN(3,3)\r\n",
"\r\n",
"BLOCK LAMQ\r\n",
" 1 0.00000000E+00 # LAMQ\r\n",
" 2 0.00000000E+00 # LAMQ\r\n",
" 3 0.00000000E+00 # LAMQ\r\n",
"\r\n",
"BLOCK LAMD\r\n",
" 1 0.00000000E+00 # LAMD\r\n",
" 2 0.00000000E+00 # LAMD\r\n",
" 3 0.00000000E+00 # LAMD\r\n",
"\r\n",
"BLOCK LAMU\r\n",
" 1 0.00000000E+00 # LAMU\r\n",
" 2 0.00000000E+00 # LAMU\r\n",
" 3 0.00000000E+00 # LAMU"
]
}
],
"source": [
"cat LesHouches.in.4VL_low"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"# SUSY Les Houches Accord 2 - 4VL Spectrum + Decays + Flavor Observables\r\n",
"# SPheno module generated by SARAH\r\n",
"# ----------------------------------------------------------------------\r\n",
"# SPheno v3.3.6 \r\n",
"# W. Porod, Comput. Phys. Commun. 153 (2003) 275-315, hep-ph/0301101\r\n",
"# W. Porod, F.Staub, Comput.Phys.Commun.183 (2012) 2458-2469, arXiv:1104.1573\r\n",
"# SARAH: 4.5.8 \r\n",
"# F. Staub; arXiv:0806.0538 (online manual)\r\n",
"# F. Staub; Comput. Phys. Commun. 181 (2010) 1077-1086; arXiv:0909.2863\r\n",
"# F. Staub; Comput. Phys. Commun. 182 (2011) 808-833; arXiv:1002.0840\r\n",
"# F. Staub; Comput. Phys. Commun. 184 (2013) 1792-1809; arXiv:1207.0906\r\n",
"# F. Staub; Comput. Phys. Commun. 185 (2014) 1773-1790; arXiv:1309.7223 \r\n",
"# Including the calculation of flavor observables based on the FlavorKit \r\n",
"# W. Porod, F. Staub, A. Vicente; Eur.Phys.J. C74 (2014) 8, 2992; arXiv:1405.1434 \r\n",
"# Two-loop masss corrections to Higgs fields based on \r\n",
"# M. D. Goodsell, K. Nickel, F. Staub; arXiv:1411.0675 \r\n",
"# M. D. Goodsell, K. Nickel, F. Staub; arXiv:1503.03098\r\n",
"# \r\n",
"# in case of problems send email to florian.staub@cern.ch and goodsell@lpthe.jussieu.fr \r\n",
"# ----------------------------------------------------------------------\r\n",
"# Created: 01.10.2015, 23:53\r\n",
"Block SPINFO # Program information\r\n",
" 1 SPhenoSARAH # spectrum calculator\r\n",
" 2 v3.3.6 # version number of SPheno\r\n",
" 9 4.5.8 # version number of SARAH\r\n",
"Block MODSEL # Input parameters\r\n",
" 1 1 # GUT scale input\r\n",
" 2 1 # Boundary conditions \r\n",
" 6 1 # switching on flavour violation\r\n",
"Block MINPAR # Input parameters\r\n",
" 1 3.99759756E+04 # Lambda1Input\r\n",
" 2 1.15198706E+00 # Lambda2Input\r\n",
" 3 5.55467608E+02 # Lambda3Input\r\n",
" 4 -7.60465292E+01 # Lambda4Input\r\n",
" 5 -7.74101866E+01 # Lambda5Input\r\n",
" 6 -3.42325000E+01 # TanBeta\r\n",
" 7 8.00000000E+02 # MQInput\r\n",
" 8 1.10000000E+03 # MUInput\r\n",
" 9 1.10000000E+03 # MDInput\r\n",
" 10 0.00000000E+00 # fdInput\r\n",
" 11 2.30000000E+00 # fuInput\r\n",
" 12 0.00000000E+00 # fdtInput\r\n",
" 13 1.60000000E+00 # futInput\r\n",
"Block gaugeGUT Q= -1.00000000E+00 # (GUT scale)\r\n",
" 1 0.00000000E+00 # g1(Q)^DRbar\r\n",
" 2 0.00000000E+00 # g2(Q)^DRbar\r\n",
" 3 0.00000000E+00 # g3(Q)^DRbar\r\n",
"Block SMINPUTS # SM parameters\r\n",
" 2 1.16637000E-05 # G_mu [GeV^-2]\r\n",
" 3 1.18700000E-01 # alpha_s(MZ)^MSbar\r\n",
" 4 9.11887000E+01 # m_Z(pole)\r\n",
" 5 4.18000000E+00 # m_b(m_b), MSbar\r\n",
" 6 1.73500000E+02 # m_t(pole)\r\n",
" 7 1.77669000E+00 # m_tau(pole)\r\n",
"Block GAUGE Q= 1.00000000E+03 # (Renormalization Scale)\r\n",
" 1 3.61543008E-01 # g1\r\n",
" 2 6.71430691E-01 # g2\r\n",
" 3 1.05673600E+00 # g3\r\n",
"Block HMIX Q= 1.00000000E+03 # (Renormalization Scale)\r\n",
" 35 -7.74101866E+01 # Lam5\r\n",
" 31 3.99759756E+04 # Lam1\r\n",
" 34 -7.60465292E+01 # Lam4\r\n",
" 33 5.55467608E+02 # Lam3\r\n",
" 32 1.15198706E+00 # Lam2\r\n",
" 20 -1.33864231E+07 # M112\r\n",
" 21 -7.53539249E+04 # M222\r\n",
" 102 6.97027906E+00 # vd\r\n",
" 103 -2.38610078E+02 # vu\r\n",
" 3 2.38711864E+02 # v\r\n",
" 10 1.54159263E+00 # betaH\r\n",
" 11 2.97046151E+00 # alphaH\r\n",
"Block FU Q= 1.00000000E+03 # (Renormalization Scale)\r\n",
" 1 2.30000000E+00 # fu\r\n",
"Block FDT Q= 1.00000000E+03 # (Renormalization Scale)\r\n",
" 1 0.00000000E+00 # fdt\r\n",
"Block FD Q= 1.00000000E+03 # (Renormalization Scale)\r\n",
" 1 0.00000000E+00 # fd\r\n",
"Block FUT Q= 1.00000000E+03 # (Renormalization Scale)\r\n",
" 1 1.60000000E+00 # fut\r\n",
"Block MD Q= 1.00000000E+03 # (Renormalization Scale)\r\n",
" 1 1.10000000E+03 # MaD\r\n",
"Block MQ Q= 1.00000000E+03 # (Renormalization Scale)\r\n",
" 1 8.00000000E+02 # MaQ\r\n",
"Block MU Q= 1.00000000E+03 # (Renormalization Scale)\r\n",
" 1 1.10000000E+03 # MaU\r\n",
"Block epsU Q= 1.00000000E+03 # (Renormalization Scale)\r\n",
" 1 1 0.00000000E+00 # Real(epsU(1,1),dp)\r\n",
" 1 2 0.00000000E+00 # Real(epsU(1,2),dp)\r\n",
" 1 3 0.00000000E+00 # Real(epsU(1,3),dp)\r\n",
" 2 1 0.00000000E+00 # Real(epsU(2,1),dp)\r\n",
" 2 2 0.00000000E+00 # Real(epsU(2,2),dp)\r\n",
" 2 3 0.00000000E+00 # Real(epsU(2,3),dp)\r\n",
" 3 1 0.00000000E+00 # Real(epsU(3,1),dp)\r\n",
" 3 2 0.00000000E+00 # Real(epsU(3,2),dp)\r\n",
" 3 3 0.00000000E+00 # Real(epsU(3,3),dp)\r\n",
"Block Yu Q= 1.00000000E+03 # (Renormalization Scale)\r\n",
" 1 1 -6.12582336E-06 # Real(Yu(1,1),dp)\r\n",
" 1 2 -1.41689997E-06 # Real(Yu(1,2),dp)\r\n",
" 1 3 -2.15349541E-08 # Real(Yu(1,3),dp)\r\n",
" 2 1 6.90438568E-04 # Real(Yu(2,1),dp)\r\n",
" 2 2 -2.98312287E-03 # Real(Yu(2,2),dp)\r\n",
" 2 3 -1.26213133E-04 # Real(Yu(2,3),dp)\r\n",
" 3 1 -5.14584556E-03 # Real(Yu(3,1),dp)\r\n",
" 3 2 3.53876517E-02 # Real(Yu(3,2),dp)\r\n",
" 3 3 -8.64558232E-01 # Real(Yu(3,3),dp)\r\n",
"Block Yd Q= 1.00000000E+03 # (Renormalization Scale)\r\n",
" 1 1 4.95981278E-04 # Real(Yd(1,1),dp)\r\n",
" 1 2 0.00000000E+00 # Real(Yd(1,2),dp)\r\n",
" 1 3 0.00000000E+00 # Real(Yd(1,3),dp)\r\n",
" 2 1 0.00000000E+00 # Real(Yd(2,1),dp)\r\n",
" 2 2 9.40654143E-03 # Real(Yd(2,2),dp)\r\n",
" 2 3 0.00000000E+00 # Real(Yd(2,3),dp)\r\n",
" 3 1 0.00000000E+00 # Real(Yd(3,1),dp)\r\n",
" 3 2 0.00000000E+00 # Real(Yd(3,2),dp)\r\n",
" 3 3 4.87429860E-01 # Real(Yd(3,3),dp)\r\n",
"Block Ye Q= 1.00000000E+03 # (Renormalization Scale)\r\n",
" 1 1 9.77192633E-05 # Real(Ye(1,1),dp)\r\n",
" 1 2 0.00000000E+00 # Real(Ye(1,2),dp)\r\n",
" 1 3 0.00000000E+00 # Real(Ye(1,3),dp)\r\n",
" 2 1 0.00000000E+00 # Real(Ye(2,1),dp)\r\n",
" 2 2 2.06537504E-02 # Real(Ye(2,2),dp)\r\n",
" 2 3 0.00000000E+00 # Real(Ye(2,3),dp)\r\n",
" 3 1 0.00000000E+00 # Real(Ye(3,1),dp)\r\n",
" 3 2 0.00000000E+00 # Real(Ye(3,2),dp)\r\n",
" 3 3 3.51113762E-01 # Real(Ye(3,3),dp)\r\n",
"Block epsD Q= 1.00000000E+03 # (Renormalization Scale)\r\n",
" 1 1 0.00000000E+00 # Real(epsD(1,1),dp)\r\n",
" 1 2 0.00000000E+00 # Real(epsD(1,2),dp)\r\n",
" 1 3 0.00000000E+00 # Real(epsD(1,3),dp)\r\n",
" 2 1 0.00000000E+00 # Real(epsD(2,1),dp)\r\n",
" 2 2 0.00000000E+00 # Real(epsD(2,2),dp)\r\n",
" 2 3 0.00000000E+00 # Real(epsD(2,3),dp)\r\n",
" 3 1 0.00000000E+00 # Real(epsD(3,1),dp)\r\n",
" 3 2 0.00000000E+00 # Real(epsD(3,2),dp)\r\n",
" 3 3 0.00000000E+00 # Real(epsD(3,3),dp)\r\n",
"Block epsE Q= 1.00000000E+03 # (Renormalization Scale)\r\n",
" 1 1 0.00000000E+00 # Real(epsE(1,1),dp)\r\n",
" 1 2 0.00000000E+00 # Real(epsE(1,2),dp)\r\n",
" 1 3 0.00000000E+00 # Real(epsE(1,3),dp)\r\n",
" 2 1 0.00000000E+00 # Real(epsE(2,1),dp)\r\n",
" 2 2 0.00000000E+00 # Real(epsE(2,2),dp)\r\n",
" 2 3 0.00000000E+00 # Real(epsE(2,3),dp)\r\n",
" 3 1 0.00000000E+00 # Real(epsE(3,1),dp)\r\n",
" 3 2 0.00000000E+00 # Real(epsE(3,2),dp)\r\n",
" 3 3 0.00000000E+00 # Real(epsE(3,3),dp)\r\n",
"Block LAMD Q= 1.00000000E+03 # (Renormalization Scale)\r\n",
" 1 0.00000000E+00 # Real(lamD(1) ,dp)\r\n",
" 2 0.00000000E+00 # Real(lamD(2) ,dp)\r\n",
" 3 0.00000000E+00 # Real(lamD(3) ,dp)\r\n",
"Block LAMQ Q= 1.00000000E+03 # (Renormalization Scale)\r\n",
" 1 0.00000000E+00 # Real(lamQ(1) ,dp)\r\n",
" 2 0.00000000E+00 # Real(lamQ(2) ,dp)\r\n",
" 3 0.00000000E+00 # Real(lamQ(3) ,dp)\r\n",
"Block LAMU Q= 1.00000000E+03 # (Renormalization Scale)\r\n",
" 1 0.00000000E+00 # Real(lamU(1) ,dp)\r\n",
" 2 0.00000000E+00 # Real(lamU(2) ,dp)\r\n",
" 3 0.00000000E+00 # Real(lamU(3) ,dp)\r\n",
"Block MASS # Mass spectrum\r\n",
"# PDG code mass particle\r\n",
" 25 1.25000002E+02 # hh_1\r\n",
" 35 2.00000000E+03 # hh_2\r\n",
" 36 2.10026000E+03 # Ah_2\r\n",
" 37 2.09099000E+03 # Hm_2\r\n",
" 23 9.11887000E+01 # VZ\r\n",
" 24 8.03497269E+01 # VWm\r\n",
" 1 5.00000000E-03 # Fd_1\r\n",
" 3 9.50000000E-02 # Fd_2\r\n",
" 5 4.18000000E+00 # Fd_3\r\n",
" 7 8.00000000E+02 # Fd_4\r\n",
" 9 1.10000000E+03 # Fd_5\r\n",
" 2 2.50000000E-03 # Fu_1\r\n",
" 4 1.27000000E+00 # Fu_2\r\n",
" 6 1.73500000E+02 # Fu_3\r\n",
" 8 8.44153626E+02 # Fu_4\r\n",
" 10 1.16656488E+03 # Fu_5\r\n",
" 11 5.10998930E-04 # Fe_1\r\n",
" 13 1.05658372E-01 # Fe_2\r\n",
" 15 1.77669000E+00 # Fe_3\r\n",
"Block SCALARMIX Q= 1.00000000E+03 # ()\r\n",
" 1 1 -1.70297076E-01 # ZH(1,1)\r\n",
" 1 2 -9.85392767E-01 # ZH(1,2)\r\n",
" 2 1 -9.85392767E-01 # ZH(2,1)\r\n",
" 2 2 1.70297076E-01 # ZH(2,2)\r\n",
"Block PSEUDOSCALARMIX Q= 1.00000000E+03 # ()\r\n",
" 1 1 2.91995502E-02 # ZA(1,1)\r\n",
" 1 2 -9.99573602E-01 # ZA(1,2)\r\n",
" 2 1 -9.99573602E-01 # ZA(2,1)\r\n",
" 2 2 -2.91995502E-02 # ZA(2,2)\r\n",
"Block CHARGEMIX Q= 1.00000000E+03 # ()\r\n",
" 1 1 2.91995502E-02 # ZP(1,1)\r\n",
" 1 2 -9.99573602E-01 # ZP(1,2)\r\n",
" 2 1 -9.99573602E-01 # ZP(2,1)\r\n",
" 2 2 -2.91995502E-02 # ZP(2,2)\r\n",
"Block UDLMIX Q= 1.00000000E+03 # ()\r\n",
" 1 1 1.00000000E+00 # Real(ZDL(1,1),dp)\r\n",
" 1 2 0.00000000E+00 # Real(ZDL(1,2),dp)\r\n",
" 1 3 0.00000000E+00 # Real(ZDL(1,3),dp)\r\n",
" 1 4 0.00000000E+00 # Real(ZDL(1,4),dp)\r\n",
" 1 5 0.00000000E+00 # Real(ZDL(1,5),dp)\r\n",
" 2 1 0.00000000E+00 # Real(ZDL(2,1),dp)\r\n",
" 2 2 1.00000000E+00 # Real(ZDL(2,2),dp)\r\n",
" 2 3 0.00000000E+00 # Real(ZDL(2,3),dp)\r\n",
" 2 4 0.00000000E+00 # Real(ZDL(2,4),dp)\r\n",
" 2 5 0.00000000E+00 # Real(ZDL(2,5),dp)\r\n",
" 3 1 0.00000000E+00 # Real(ZDL(3,1),dp)\r\n",
" 3 2 0.00000000E+00 # Real(ZDL(3,2),dp)\r\n",
" 3 3 1.00000000E+00 # Real(ZDL(3,3),dp)\r\n",
" 3 4 0.00000000E+00 # Real(ZDL(3,4),dp)\r\n",
" 3 5 0.00000000E+00 # Real(ZDL(3,5),dp)\r\n",
" 4 1 0.00000000E+00 # Real(ZDL(4,1),dp)\r\n",
" 4 2 0.00000000E+00 # Real(ZDL(4,2),dp)\r\n",
" 4 3 0.00000000E+00 # Real(ZDL(4,3),dp)\r\n",
" 4 4 1.00000000E+00 # Real(ZDL(4,4),dp)\r\n",
" 4 5 0.00000000E+00 # Real(ZDL(4,5),dp)\r\n",
" 5 1 0.00000000E+00 # Real(ZDL(5,1),dp)\r\n",
" 5 2 0.00000000E+00 # Real(ZDL(5,2),dp)\r\n",
" 5 3 0.00000000E+00 # Real(ZDL(5,3),dp)\r\n",
" 5 4 0.00000000E+00 # Real(ZDL(5,4),dp)\r\n",
" 5 5 1.00000000E+00 # Real(ZDL(5,5),dp)\r\n",
"Block UDRMIX Q= 1.00000000E+03 # ()\r\n",
" 1 1 1.00000000E+00 # Real(ZDR(1,1),dp)\r\n",
" 1 2 0.00000000E+00 # Real(ZDR(1,2),dp)\r\n",
" 1 3 0.00000000E+00 # Real(ZDR(1,3),dp)\r\n",
" 1 4 0.00000000E+00 # Real(ZDR(1,4),dp)\r\n",
" 1 5 0.00000000E+00 # Real(ZDR(1,5),dp)\r\n",
" 2 1 0.00000000E+00 # Real(ZDR(2,1),dp)\r\n",
" 2 2 1.00000000E+00 # Real(ZDR(2,2),dp)\r\n",
" 2 3 0.00000000E+00 # Real(ZDR(2,3),dp)\r\n",
" 2 4 0.00000000E+00 # Real(ZDR(2,4),dp)\r\n",
" 2 5 0.00000000E+00 # Real(ZDR(2,5),dp)\r\n",
" 3 1 0.00000000E+00 # Real(ZDR(3,1),dp)\r\n",
" 3 2 0.00000000E+00 # Real(ZDR(3,2),dp)\r\n",
" 3 3 1.00000000E+00 # Real(ZDR(3,3),dp)\r\n",
" 3 4 0.00000000E+00 # Real(ZDR(3,4),dp)\r\n",
" 3 5 0.00000000E+00 # Real(ZDR(3,5),dp)\r\n",
" 4 1 0.00000000E+00 # Real(ZDR(4,1),dp)\r\n",
" 4 2 0.00000000E+00 # Real(ZDR(4,2),dp)\r\n",
" 4 3 0.00000000E+00 # Real(ZDR(4,3),dp)\r\n",
" 4 4 1.00000000E+00 # Real(ZDR(4,4),dp)\r\n",
" 4 5 0.00000000E+00 # Real(ZDR(4,5),dp)\r\n",
" 5 1 0.00000000E+00 # Real(ZDR(5,1),dp)\r\n",
" 5 2 0.00000000E+00 # Real(ZDR(5,2),dp)\r\n",
" 5 3 0.00000000E+00 # Real(ZDR(5,3),dp)\r\n",
" 5 4 0.00000000E+00 # Real(ZDR(5,4),dp)\r\n",
" 5 5 1.00000000E+00 # Real(ZDR(5,5),dp)\r\n",
"Block UULMIX Q= 1.00000000E+03 # ()\r\n",
" 1 1 -9.74272160E-01 # Real(ZUL(1,1),dp)\r\n",
" 1 2 -2.25348678E-01 # Real(ZUL(1,2),dp)\r\n",
" 1 3 -3.42499367E-03 # Real(ZUL(1,3),dp)\r\n",
" 1 4 0.00000000E+00 # Real(ZUL(1,4),dp)\r\n",
" 1 5 0.00000000E+00 # Real(ZUL(1,5),dp)\r\n",
" 2 1 2.25296231E-01 # Real(ZUL(2,1),dp)\r\n",
" 2 2 -9.73419462E-01 # Real(ZUL(2,2),dp)\r\n",
" 2 3 -4.11844653E-02 # Real(ZUL(2,3),dp)\r\n",
" 2 4 0.00000000E+00 # Real(ZUL(2,4),dp)\r\n",
" 2 5 0.00000000E+00 # Real(ZUL(2,5),dp)\r\n",
" 3 1 -5.94690932E-03 # Real(ZUL(3,1),dp)\r\n",
" 3 2 4.08965161E-02 # Real(ZUL(3,2),dp)\r\n",
" 3 3 -9.99145690E-01 # Real(ZUL(3,3),dp)\r\n",
" 3 4 0.00000000E+00 # Real(ZUL(3,4),dp)\r\n",
" 3 5 0.00000000E+00 # Real(ZUL(3,5),dp)\r\n",
" 4 1 0.00000000E+00 # Real(ZUL(4,1),dp)\r\n",
" 4 2 0.00000000E+00 # Real(ZUL(4,2),dp)\r\n",
" 4 3 0.00000000E+00 # Real(ZUL(4,3),dp)\r\n",
" 4 4 9.37915637E-01 # Real(ZUL(4,4),dp)\r\n",
" 4 5 -3.46863458E-01 # Real(ZUL(4,5),dp)\r\n",
" 5 1 0.00000000E+00 # Real(ZUL(5,1),dp)\r\n",
" 5 2 0.00000000E+00 # Real(ZUL(5,2),dp)\r\n",
" 5 3 0.00000000E+00 # Real(ZUL(5,3),dp)\r\n",
" 5 4 3.46863458E-01 # Real(ZUL(5,4),dp)\r\n",
" 5 5 9.37915637E-01 # Real(ZUL(5,5),dp)\r\n",
"Block UURMIX Q= 1.00000000E+03 # ()\r\n",
" 1 1 1.00000000E+00 # Real(ZUR(1,1),dp)\r\n",
" 1 2 -2.16840434E-19 # Real(ZUR(1,2),dp)\r\n",
" 1 3 0.00000000E+00 # Real(ZUR(1,3),dp)\r\n",
" 1 4 0.00000000E+00 # Real(ZUR(1,4),dp)\r\n",
" 1 5 0.00000000E+00 # Real(ZUR(1,5),dp)\r\n",
" 2 1 2.16840434E-19 # Real(ZUR(2,1),dp)\r\n",
" 2 2 1.00000000E+00 # Real(ZUR(2,2),dp)\r\n",
" 2 3 -0.00000000E+00 # Real(ZUR(2,3),dp)\r\n",
" 2 4 -0.00000000E+00 # Real(ZUR(2,4),dp)\r\n",
" 2 5 -0.00000000E+00 # Real(ZUR(2,5),dp)\r\n",
" 3 1 0.00000000E+00 # Real(ZUR(3,1),dp)\r\n",
" 3 2 0.00000000E+00 # Real(ZUR(3,2),dp)\r\n",
" 3 3 1.00000000E+00 # Real(ZUR(3,3),dp)\r\n",
" 3 4 0.00000000E+00 # Real(ZUR(3,4),dp)\r\n",
" 3 5 0.00000000E+00 # Real(ZUR(3,5),dp)\r\n",
" 4 1 0.00000000E+00 # Real(ZUR(4,1),dp)\r\n",
" 4 2 0.00000000E+00 # Real(ZUR(4,2),dp)\r\n",
" 4 3 0.00000000E+00 # Real(ZUR(4,3),dp)\r\n",
" 4 4 9.99783125E-01 # Real(ZUR(4,4),dp)\r\n",
" 4 5 -2.08255508E-02 # Real(ZUR(4,5),dp)\r\n",
" 5 1 0.00000000E+00 # Real(ZUR(5,1),dp)\r\n",
" 5 2 0.00000000E+00 # Real(ZUR(5,2),dp)\r\n",
" 5 3 0.00000000E+00 # Real(ZUR(5,3),dp)\r\n",
" 5 4 2.08255508E-02 # Real(ZUR(5,4),dp)\r\n",
" 5 5 9.99783125E-01 # Real(ZUR(5,5),dp)\r\n",
"Block UELMIX Q= 1.00000000E+03 # ()\r\n",
" 1 1 1.00000000E+00 # Real(ZEL(1,1),dp)\r\n",
" 1 2 0.00000000E+00 # Real(ZEL(1,2),dp)\r\n",
" 1 3 0.00000000E+00 # Real(ZEL(1,3),dp)\r\n",
" 2 1 0.00000000E+00 # Real(ZEL(2,1),dp)\r\n",
" 2 2 1.00000000E+00 # Real(ZEL(2,2),dp)\r\n",
" 2 3 0.00000000E+00 # Real(ZEL(2,3),dp)\r\n",
" 3 1 0.00000000E+00 # Real(ZEL(3,1),dp)\r\n",
" 3 2 0.00000000E+00 # Real(ZEL(3,2),dp)\r\n",
" 3 3 1.00000000E+00 # Real(ZEL(3,3),dp)\r\n",
"Block UERMIX Q= 1.00000000E+03 # ()\r\n",
" 1 1 1.00000000E+00 # Real(ZER(1,1),dp)\r\n",
" 1 2 0.00000000E+00 # Real(ZER(1,2),dp)\r\n",
" 1 3 0.00000000E+00 # Real(ZER(1,3),dp)\r\n",
" 2 1 0.00000000E+00 # Real(ZER(2,1),dp)\r\n",
" 2 2 1.00000000E+00 # Real(ZER(2,2),dp)\r\n",
" 2 3 0.00000000E+00 # Real(ZER(2,3),dp)\r\n",
" 3 1 0.00000000E+00 # Real(ZER(3,1),dp)\r\n",
" 3 2 0.00000000E+00 # Real(ZER(3,2),dp)\r\n",
" 3 3 1.00000000E+00 # Real(ZER(3,3),dp)\r\n",
"Block SPheno # SPheno internal parameters \r\n",
" 1 -1.00000000E+00 # ErrorLevel\r\n",
" 2 0.00000000E+00 # SPA_conventions\r\n",
" 11 1.00000000E+00 # Branching ratios\r\n",
" 13 1.00000000E+00 # 3 Body decays\r\n",
" 31 -1.00000000E+00 # GUT scale\r\n",
" 33 1.00000000E+03 # SUSY scale\r\n",
" 34 1.00000000E-04 # Precision\r\n",
" 35 4.00000000E+01 # Iterations\r\n",
" 38 2.00000000E+00 # RGE level\r\n",
" 40 7.29735257E-03 # Alpha\r\n",
" 41 2.49520000E+00 # Gamma_Z\r\n",
" 42 2.06000000E+00 # Gamma_W\r\n",
" 50 1.00000000E+00 # Rotate negative fermion masses\r\n",
" 51 0.00000000E+00 # Switch to SCKM matrix\r\n",
" 52 1.00000000E+00 # Ignore negative masses\r\n",
" 53 0.00000000E+00 # Ignore negative masses at MZ\r\n",
" 55 0.00000000E+00 # Calculate one loop masses\r\n",
" 56 1.00000000E+00 # Calculate two-loop Higgs masses\r\n",
" 57 0.00000000E+00 # Calculate low energy\r\n",
" 60 1.00000000E+00 # Include kinetic mixing\r\n",
" 65 1.00000000E+00 # Solution of tadpole equation\r\n",
"Block HiggsLHC7 # Higgs production cross section at LHC7 [pb] \r\n",
" 1 25 4.42505408E+01 # Gluon fusion \r\n",
" 2 25 1.16124208E+00 # Vector boson fusion \r\n",
" 3 25 5.19844836E-01 # W-H production \r\n",
" 4 25 3.03894630E-01 # Z-H production \r\n",
" 5 25 8.36979352E-02 # t-t-H production \r\n",
"Block HiggsLHC8 # Higgs production cross section at LHC8 [pb] \r\n",
" 1 25 5.60503082E+01 # Gluon fusion \r\n",
" 2 25 1.49728527E+00 # Vector boson fusion \r\n",
" 3 25 6.66408136E-01 # W-H production \r\n",
" 4 25 3.79554684E-01 # Z-H production \r\n",
" 5 25 1.26688259E-01 # t-t-H production \r\n",
"Block HiggsBoundsInputHiggsCouplingsFermions # \r\n",
" 3.40143408E+01 0.00000000E+00 3 25 5 5 # h_1 b b coupling \r\n",
" 3.40143408E+01 0.00000000E+00 3 25 3 3 # h_1 s s coupling \r\n",
" 9.71827499E-01 0.00000000E+00 3 25 6 6 # h_1 t t coupling \r\n",
" 9.71827499E-01 0.00000000E+00 3 25 4 4 # h_1 c c coupling \r\n",
" 3.40143408E+01 0.00000000E+00 3 25 15 15 # h_1 tau tau coupling \r\n",
" 3.40143408E+01 0.00000000E+00 3 25 13 13 # h_1 mu mu coupling \r\n",
" 1.13884972E+03 0.00000000E+00 3 35 5 5 # h_2 b b coupling \r\n",
" 1.13884972E+03 0.00000000E+00 3 35 3 3 # h_2 s s coupling \r\n",
" 2.90258419E-02 0.00000000E+00 3 35 6 6 # h_2 t t coupling \r\n",
" 2.90258419E-02 0.00000000E+00 3 35 4 4 # h_2 c c coupling \r\n",
" 1.13884972E+03 0.00000000E+00 3 35 15 15 # h_2 tau tau coupling \r\n",
" 1.13884972E+03 0.00000000E+00 3 35 13 13 # h_2 mu mu coupling \r\n",
" 0.00000000E+00 1.17186406E+03 3 36 5 5 # A_2 b b coupling \r\n",
" 0.00000000E+00 1.17186406E+03 3 36 3 3 # A_2 s s coupling \r\n",
" 0.00000000E+00 8.53341302E-04 3 36 6 6 # A_2 t t coupling \r\n",
" 0.00000000E+00 8.53341302E-04 3 36 4 4 # A_2 c c coupling \r\n",
" 0.00000000E+00 1.17186406E+03 3 36 15 15 # A_2 tau tau coupling \r\n",
" 0.00000000E+00 1.17186406E+03 3 36 13 13 # A_2 mu mu coupling \r\n",
"Block HiggsBoundsInputHiggsCouplingsBosons # \r\n",
" 9.60400000E-01 3 25 24 24 # h_1 W W coupling \r\n",
" 9.60400000E-01 3 25 23 23 # h_1 Z Z coupling \r\n",
" 0.00000000E+00 3 25 23 22 # h_1 Z gamma coupling \r\n",
" 1.11520458E+01 3 25 22 22 # h_1 gamma gamma coupling \r\n",
" 2.87900483E+00 3 25 21 21 # h_1 g g coupling \r\n",
" 0.00000000E+00 4 25 21 21 23 # h_1 g g Z coupling \r\n",
" 3.95999999E-02 3 35 24 24 # h_2 W W coupling \r\n",
" 3.95999999E-02 3 35 23 23 # h_2 Z Z coupling \r\n",
" 0.00000000E+00 3 35 23 22 # h_2 Z gamma coupling \r\n",
" 1.12663703E+02 3 35 22 22 # h_2 gamma gamma coupling \r\n",
" 2.94096247E+02 3 35 21 21 # h_2 g g coupling \r\n",
" 0.00000000E+00 4 35 21 21 23 # h_2 g g Z coupling \r\n",
" 0.00000000E+00 3 36 24 24 # A_2 W W coupling \r\n",
" 0.00000000E+00 3 36 23 23 # A_2 Z Z coupling \r\n",
" 0.00000000E+00 3 36 23 22 # A_2 Z gamma coupling \r\n",
" 1.25992045E-05 3 36 22 22 # A_2 gamma gamma coupling \r\n",
" 6.68876815E-06 3 36 21 21 # A_2 g g coupling \r\n",
" 0.00000000E+00 4 36 21 21 23 # A_2 g g Z coupling \r\n",
" 0.00000000E+00 3 25 25 23 # h_1 h_1 Z coupling \r\n",
" 0.00000000E+00 3 25 35 23 # h_1 h_2 Z coupling \r\n",
" 4.49760176E-02 3 25 36 23 # h_1 A_2 Z coupling \r\n",
" 0.00000000E+00 3 35 25 23 # h_2 h_1 Z coupling \r\n",
" 0.00000000E+00 3 35 35 23 # h_2 h_2 Z coupling \r\n",
" 1.09078201E+00 3 35 36 23 # h_2 A_2 Z coupling \r\n",
" 0.00000000E+00 3 36 36 23 # A_2 A_2 Z coupling \r\n",
"Block EFFHIGGSCOUPLINGS # values of loop-induced couplings \r\n",
" 25 22 22 0.35385989E-04 # H-Photon-Photon \r\n",
" 25 21 21 0.49706469E-03 # H-Gluon-Gluon \r\n",
" 25 22 23 0.00000000E+00 # H-Photon-Z (not yet calculated by SPheno) \r\n",
" 35 22 22 0.24632007E-03 # H-Photon-Photon \r\n",
" 35 21 21 0.60530869E-02 # H-Gluon-Gluon \r\n",
" 35 22 23 0.00000000E+00 # H-Photon-Z (not yet calculated by SPheno) \r\n",
" 36 22 22 0.93168163E-07 # A-Photon-Photon \r\n",
" 36 21 21 0.13311041E-05 # A-Gluon-Gluon \r\n",
" 36 22 23 0.00000000E+00 # A-Photon-Z (not yet calculated by SPheno) \r\n",
"Block SPhenoLowEnergy # low energy observables \r\n",
" 20 0.00000000E+00 # (g-2)_e\r\n",
" 21 0.00000000E+00 # (g-2)_mu\r\n",
" 22 0.00000000E+00 # (g-2)_tau\r\n",
" 23 0.00000000E+00 # EDM(e)\r\n",
" 24 0.00000000E+00 # EDM(mu)\r\n",
" 25 0.00000000E+00 # EDM(tau)\r\n",
" 39 0.00000000E+00 # delta(rho)\r\n",
"Block FlavorKitQFV # quark flavor violating observables \r\n",
" 200 0.00000000E+00 # BR(B->X_s gamma)\r\n",
" 201 0.00000000E+00 # BR(B->X_s gamma)/BR(B->X_s gamma)_SM\r\n",
" 300 0.00000000E+00 # BR(D->mu nu)\r\n",
" 301 0.00000000E+00 # BR(D->mu nu)/BR(D->mu nu)_SM\r\n",
" 400 0.00000000E+00 # BR(Ds->mu nu)\r\n",
" 401 0.00000000E+00 # BR(Ds->mu nu)/BR(Ds->mu nu)_SM\r\n",
" 402 0.00000000E+00 # BR(Ds->tau nu)\r\n",
" 403 0.00000000E+00 # BR(Ds->tau nu)/BR(Ds->tau nu)_SM\r\n",
" 500 0.00000000E+00 # BR(B->mu nu)\r\n",
" 501 0.00000000E+00 # BR(B->mu nu)/BR(B->mu nu)_SM\r\n",
" 502 0.00000000E+00 # BR(B->tau nu)\r\n",
" 503 0.00000000E+00 # BR(B->tau nu)/BR(B->tau nu)_SM\r\n",
" 600 0.00000000E+00 # BR(K->mu nu)\r\n",
" 601 0.00000000E+00 # BR(K->mu nu)/BR(K->mu nu)_SM\r\n",
" 602 0.00000000E+00 # R_K = BR(K->e nu)/(K->mu nu)\r\n",
" 603 0.00000000E+00 # R_K^SM = BR(K->e nu)_SM/(K->mu nu)_SM\r\n",
" 1900 0.00000000E+00 # Delta(M_Bs)\r\n",
" 1901 0.00000000E+00 # Delta(M_Bs)/Delta(M_Bs)_SM\r\n",
" 1902 0.00000000E+00 # Delta(M_Bd)\r\n",
" 1903 0.00000000E+00 # Delta(M_Bd)/Delta(M_Bd)_SM\r\n",
" 4000 0.00000000E+00 # BR(B^0_d->e e)\r\n",
" 4001 0.00000000E+00 # BR(B^0_d->e e)/BR(B^0_d->e e)_SM\r\n",
" 4002 0.00000000E+00 # BR(B^0_s->e e)\r\n",
" 4003 0.00000000E+00 # BR(B^0_s->e e)/BR(B^0_s->e e)_SM\r\n",
" 4004 0.00000000E+00 # BR(B^0_d->mu mu)\r\n",
" 4005 0.00000000E+00 # BR(B^0_d->mu mu)/BR(B^0_d->mu mu)_SM\r\n",
" 4006 0.00000000E+00 # BR(B^0_s->mu mu)\r\n",
" 4007 0.00000000E+00 # BR(B^0_s->mu mu)/BR(B^0_s->mu mu)_SM\r\n",
" 4008 0.00000000E+00 # BR(B^0_d->tau tau)\r\n",
" 4009 0.00000000E+00 # BR(B^0_d->tau tau)/BR(B^0_d->tau tau)_SM\r\n",
" 4010 0.00000000E+00 # BR(B^0_s->tau tau)\r\n",
" 4011 0.00000000E+00 # BR(B^0_s->tau tau)/BR(B^0_s->tau tau)_SM\r\n",
" 5000 0.00000000E+00 # BR(B-> s e e)\r\n",
" 5001 0.00000000E+00 # BR(B-> s e e)/BR(B-> s e e)_SM\r\n",
" 5002 0.00000000E+00 # BR(B-> s mu mu)\r\n",
" 5003 0.00000000E+00 # BR(B-> s mu mu)/BR(B-> s mu mu)_SM\r\n",
" 6000 0.00000000E+00 # BR(B -> K mu mu)\r\n",
" 6001 0.00000000E+00 # BR(B -> K mu mu)/BR(B -> K mu mu)_SM\r\n",
" 7000 0.00000000E+00 # BR(B->s nu nu)\r\n",
" 7001 0.00000000E+00 # BR(B->s nu nu)/BR(B->s nu nu)_SM\r\n",
" 7002 0.00000000E+00 # BR(B->D nu nu)\r\n",
" 7003 0.00000000E+00 # BR(B->D nu nu)/BR(B->D nu nu)_SM\r\n",
" 8000 0.00000000E+00 # BR(K^+ -> pi^+ nu nu)\r\n",
" 8001 0.00000000E+00 # BR(K^+ -> pi^+ nu nu)/BR(K^+ -> pi^+ nu nu)_SM\r\n",
" 8002 0.00000000E+00 # BR(K_L -> pi^0 nu nu)\r\n",
" 8003 0.00000000E+00 # BR(K_L -> pi^0 nu nu)/BR(K_L -> pi^0 nu nu)_SM\r\n",
" 9100 0.00000000E+00 # Delta(M_K)\r\n",
" 9102 0.00000000E+00 # Delta(M_K)/Delta(M_K)_SM\r\n",
" 9103 0.00000000E+00 # epsilon_K\r\n",
" 9104 0.00000000E+00 # epsilon_K/epsilon_K^SM\r\n",
"Block FlavorKitLFV # lepton flavor violating observables \r\n",
" 701 0.00000000E+00 # BR(mu->e gamma)\r\n",
" 702 0.00000000E+00 # BR(tau->e gamma)\r\n",
" 703 0.00000000E+00 # BR(tau->mu gamma)\r\n",
" 800 0.00000000E+00 # CR(mu-e, Al)\r\n",
" 801 0.00000000E+00 # CR(mu-e, Ti)\r\n",
" 802 0.00000000E+00 # CR(mu-e, Sr)\r\n",
" 803 0.00000000E+00 # CR(mu-e, Sb)\r\n",
" 804 0.00000000E+00 # CR(mu-e, Au)\r\n",
" 805 0.00000000E+00 # CR(mu-e, Pb)\r\n",
" 901 0.00000000E+00 # BR(mu->3e)\r\n",
" 902 0.00000000E+00 # BR(tau->3e)\r\n",
" 903 0.00000000E+00 # BR(tau->3mu)\r\n",
" 904 0.00000000E+00 # BR(tau- -> e- mu+ mu-)\r\n",
" 905 0.00000000E+00 # BR(tau- -> mu- e+ e-)\r\n",
" 906 0.00000000E+00 # BR(tau- -> e+ mu- mu-)\r\n",
" 907 0.00000000E+00 # BR(tau- -> mu+ e- e-)\r\n",
" 1001 0.00000000E+00 # BR(Z->e mu)\r\n",
" 1002 0.00000000E+00 # BR(Z->e tau)\r\n",
" 1003 0.00000000E+00 # BR(Z->mu tau)\r\n",
" 1101 0.00000000E+00 # BR(h->e mu)\r\n",
" 1102 0.00000000E+00 # BR(h->e tau)\r\n",
" 1103 0.00000000E+00 # BR(h->mu tau)\r\n",
" 2001 0.00000000E+00 # BR(tau->e pi)\r\n",
" 2002 0.00000000E+00 # BR(tau->e eta)\r\n",
" 2003 0.00000000E+00 # BR(tau->e eta')\r\n",
" 2004 0.00000000E+00 # BR(tau->mu pi)\r\n",
" 2005 0.00000000E+00 # BR(tau->mu eta)\r\n",
" 2006 0.00000000E+00 # BR(tau->mu eta')\r\n",
"Block FWCOEF Q= 1.60000000E+02 # Wilson coefficients at scale Q \r\n",
" 030 4422 00 0 0.00000000E+00 # coeffC7sm\r\n",
" 030 4422 00 2 0.00000000E+00 # coeffC7\r\n",
" 030 4322 00 2 0.00000000E+00 # coeffC7p\r\n",
" 030 6421 00 0 0.00000000E+00 # coeffC8sm\r\n",
" 030 6421 00 2 0.00000000E+00 # coeffC8\r\n",
" 030 6321 00 2 0.00000000E+00 # coeffC8p\r\n",
" 0305111 4133 00 0 0.00000000E+00 # coeffC9eeSM\r\n",
" 0305111 4133 00 2 0.00000000E+00 # coeffC9ee\r\n",
" 0305111 4233 00 2 0.00000000E+00 # coeffC9Pee\r\n",
" 0305111 4137 00 0 0.00000000E+00 # coeffC10eeSM\r\n",
" 0305111 4137 00 2 0.00000000E+00 # coeffC10ee\r\n",
" 0305111 4237 00 2 0.00000000E+00 # coeffC10Pee\r\n",
" 0305131 4133 00 0 0.00000000E+00 # coeffC9mumuSM\r\n",
" 0305131 4133 00 2 0.00000000E+00 # coeffC9mumu\r\n",
" 0305131 4233 00 2 0.00000000E+00 # coeffC9Pmumu\r\n",
" 0305131 4137 00 0 0.00000000E+00 # coeffC10mumuSM\r\n",
" 0305131 4137 00 2 0.00000000E+00 # coeffC10mumu\r\n",
" 0305131 4237 00 2 0.00000000E+00 # coeffC10Pmumu\r\n",
" 0305121 4137 00 0 0.00000000E+00 # coeffC11nu1nu1SM\r\n",
" 0305121 4137 00 2 0.00000000E+00 # coeffC11nu1nu1\r\n",
" 0305121 4237 00 2 0.00000000E+00 # coeffC11Pnu1nu1\r\n",
" 0305141 4137 00 0 0.00000000E+00 # coeffC11nu2nu2SM\r\n",
" 0305141 4137 00 2 0.00000000E+00 # coeffC11nu2nu2\r\n",
" 0305141 4237 00 2 0.00000000E+00 # coeffC11Pnu2nu2\r\n",
" 0305161 4137 00 0 0.00000000E+00 # coeffC11nu3nu3SM\r\n",
" 0305161 4137 00 2 0.00000000E+00 # coeffC11nu3nu3\r\n",
" 0305161 4237 00 2 0.00000000E+00 # coeffC11Pnu3nu3\r\n",
"Block IMFWCOEF Q= 1.60000000E+02 # Im(Wilson coefficients) at scale Q \r\n",
" 030 4422 00 0 0.00000000E+00 # coeffC7sm\r\n",
" 030 4422 00 2 0.00000000E+00 # coeffC7\r\n",
" 030 4322 00 2 0.00000000E+00 # coeffC7p\r\n",
" 030 6421 00 0 0.00000000E+00 # coeffC8sm\r\n",
" 030 6421 00 2 0.00000000E+00 # coeffC8\r\n",
" 030 6321 00 2 0.00000000E+00 # coeffC8p\r\n",
" 0305111 4133 00 0 0.00000000E+00 # coeffC9eeSM\r\n",
" 0305111 4133 00 2 0.00000000E+00 # coeffC9ee\r\n",
" 0305111 4233 00 2 0.00000000E+00 # coeffC9Pee\r\n",
" 0305111 4137 00 0 0.00000000E+00 # coeffC10eeSM\r\n",
" 0305111 4137 00 2 0.00000000E+00 # coeffC10ee\r\n",
" 0305111 4237 00 2 0.00000000E+00 # coeffC10Pee\r\n",
" 0305131 4133 00 0 0.00000000E+00 # coeffC9mumuSM\r\n",
" 0305131 4133 00 2 0.00000000E+00 # coeffC9mumu\r\n",
" 0305131 4233 00 2 0.00000000E+00 # coeffC9Pmumu\r\n",
" 0305131 4137 00 0 0.00000000E+00 # coeffC10mumuSM\r\n",
" 0305131 4137 00 2 0.00000000E+00 # coeffC10mumu\r\n",
" 0305131 4237 00 2 0.00000000E+00 # coeffC10Pmumu\r\n",
" 0305121 4137 00 0 0.00000000E+00 # coeffC11nu1nu1SM\r\n",
" 0305121 4137 00 2 0.00000000E+00 # coeffC11nu1nu1\r\n",
" 0305121 4237 00 2 0.00000000E+00 # coeffC11Pnu1nu1\r\n",
" 0305141 4137 00 0 0.00000000E+00 # coeffC11nu2nu2SM\r\n",
" 0305141 4137 00 2 0.00000000E+00 # coeffC11nu2nu2\r\n",
" 0305141 4237 00 2 0.00000000E+00 # coeffC11Pnu2nu2\r\n",
" 0305161 4137 00 0 0.00000000E+00 # coeffC11nu3nu3SM\r\n",
" 0305161 4137 00 2 0.00000000E+00 # coeffC11nu3nu3\r\n",
" 0305161 4237 00 2 0.00000000E+00 # coeffC11Pnu3nu3\r\n",
"DECAY 4 4.00058449E-13 # Fu_2\r\n",
"# BR NDA ID1 ID2\r\n",
"# BR NDA ID1 ID2 ID3\r\n",
" 3.05502649E-02 3 2 -1 1 # BR(Fu_2 -> Fu_1 Fd_1^* Fd_1 )\r\n",
" 5.45955119E-01 3 2 -1 3 # BR(Fu_2 -> Fu_1 Fd_1^* Fd_2 )\r\n",
" 1.56486271E-03 3 2 -3 1 # BR(Fu_2 -> Fu_1 Fd_2^* Fd_1 )\r\n",
" 2.79270083E-02 3 2 -3 3 # BR(Fu_2 -> Fu_1 Fd_2^* Fd_2 )\r\n",
" 1.07295210E-02 3 1 -11 12 # BR(Fu_2 -> Fd_1 Fe_1^* Fv_1 )\r\n",
" 1.01645136E-02 3 1 -13 14 # BR(Fu_2 -> Fd_1 Fe_2^* Fv_2 )\r\n",
" 1.91744818E-01 3 3 -11 12 # BR(Fu_2 -> Fd_2 Fe_1^* Fv_1 )\r\n",
" 1.81363892E-01 3 3 -13 14 # BR(Fu_2 -> Fd_2 Fe_2^* Fv_2 )\r\n",
"DECAY 6 1.59604378E+00 # Fu_3\r\n",
"# BR NDA ID1 ID2\r\n",
" 3.54387767E-05 2 1 24 # BR(Fu_3 -> Fd_1 VWm^* )\r\n",
" 1.67597777E-03 2 3 24 # BR(Fu_3 -> Fd_2 VWm^* )\r\n",
" 9.98288583E-01 2 5 24 # BR(Fu_3 -> Fd_3 VWm^* )\r\n",
"# BR NDA ID1 ID2 ID3\r\n",
" 4.54660816E-15 3 2 -1 5 # BR(Fu_3 -> Fu_1 Fd_1^* Fd_3 )\r\n",
" 8.74913774E-14 3 2 -3 5 # BR(Fu_3 -> Fu_1 Fd_2^* Fd_3 )\r\n",
" 5.40174553E-14 3 2 -5 5 # BR(Fu_3 -> Fu_1 Fd_3^* Fd_3 )\r\n",
" 1.63195637E-12 3 4 -3 5 # BR(Fu_3 -> Fu_2 Fd_2^* Fd_3 )\r\n",
" 7.80720122E-12 3 4 -5 5 # BR(Fu_3 -> Fu_2 Fd_3^* Fd_3 )\r\n",
" 4.07482129E-15 3 3 -15 16 # BR(Fu_3 -> Fd_2 Fe_3^* Fv_3 )\r\n",
" 2.76867263E-12 3 5 -13 14 # BR(Fu_3 -> Fd_3 Fe_2^* Fv_2 )\r\n",
" 7.99479340E-10 3 5 -15 16 # BR(Fu_3 -> Fd_3 Fe_3^* Fv_3 )\r\n",
"DECAY 8 5.23129292E-04 # Fu_4\r\n",
"# BR NDA ID1 ID2\r\n",
"# BR NDA ID1 ID2 ID3\r\n",
" 3.17047749E-01 3 2 -1 7 # BR(Fu_4 -> Fu_1 Fd_1^* Fd_4 )\r\n",
" 1.69615550E-02 3 2 -3 7 # BR(Fu_4 -> Fu_1 Fd_2^* Fd_4 )\r\n",
" 3.76923823E-06 3 2 -5 7 # BR(Fu_4 -> Fu_1 Fd_3^* Fd_4 )\r\n",
" 1.68967018E-02 3 4 -1 7 # BR(Fu_4 -> Fu_2 Fd_1^* Fd_4 )\r\n",
" 3.15417407E-01 3 4 -3 7 # BR(Fu_4 -> Fu_2 Fd_2^* Fd_4 )\r\n",
" 5.43101905E-04 3 4 -5 7 # BR(Fu_4 -> Fu_2 Fd_3^* Fd_4 )\r\n",
" 1.11337854E-01 3 7 -11 12 # BR(Fu_4 -> Fd_4 Fe_1^* Fv_1 )\r\n",
" 1.11337302E-01 3 7 -13 14 # BR(Fu_4 -> Fd_4 Fe_2^* Fv_2 )\r\n",
" 1.10454560E-01 3 7 -15 16 # BR(Fu_4 -> Fd_4 Fe_3^* Fv_3 )\r\n",
"DECAY 10 1.55835686E+01 # Fu_5\r\n",
"# BR NDA ID1 ID2\r\n",
" 6.29167005E-01 2 7 24 # BR(Fu_5 -> Fd_4 VWm^* )\r\n",
" 1.69289029E-01 2 8 25 # BR(Fu_5 -> Fu_4 hh_1 )\r\n",
" 2.00055677E-01 2 8 23 # BR(Fu_5 -> Fu_4 VZ )\r\n",
"# BR NDA ID1 ID2 ID3\r\n",
" 8.83741211E-15 3 2 -3 7 # BR(Fu_5 -> Fu_1 Fd_2^* Fd_4 )\r\n",
" 5.47758582E-15 3 2 -5 7 # BR(Fu_5 -> Fu_1 Fd_3^* Fd_4 )\r\n",
" 1.64901952E-13 3 4 -3 7 # BR(Fu_5 -> Fu_2 Fd_2^* Fd_4 )\r\n",
" 7.91968427E-13 3 4 -5 7 # BR(Fu_5 -> Fu_2 Fd_3^* Fd_4 )\r\n",
" 5.27459630E-08 3 6 -1 7 # BR(Fu_5 -> Fu_3 Fd_1^* Fd_4 )\r\n",
" 2.49445618E-06 3 6 -3 7 # BR(Fu_5 -> Fu_3 Fd_2^* Fd_4 )\r\n",
" 1.48573788E-03 3 6 -5 7 # BR(Fu_5 -> Fu_3 Fd_3^* Fd_4 )\r\n",
" 3.71171914E-15 3 8 -1 1 # BR(Fu_5 -> Fu_4 Fd_1^* Fd_1 )\r\n",
" 1.33506950E-12 3 8 -3 3 # BR(Fu_5 -> Fu_4 Fd_2^* Fd_2 )\r\n",
" 3.57509046E-09 3 8 -5 5 # BR(Fu_5 -> Fu_4 Fd_3^* Fd_3 )\r\n",
" 2.14545979E-12 3 8 -13 13 # BR(Fu_5 -> Fu_4 Fe_2^* Fe_2 )\r\n",
" 6.19733557E-10 3 8 -15 15 # BR(Fu_5 -> Fu_4 Fe_3^* Fe_3 )\r\n",
" 4.13030054E-15 3 4 -4 8 # BR(Fu_5 -> Fu_2 Fu_2^* Fu_4 )\r\n",
" 2.79660877E-13 3 7 -13 14 # BR(Fu_5 -> Fd_4 Fe_2^* Fv_2 )\r\n",
" 8.08116086E-11 3 7 -15 16 # BR(Fu_5 -> Fd_4 Fe_3^* Fv_3 )\r\n",
"DECAY 13 1.27670536E-18 # Fe_2\r\n",
"# BR NDA ID1 ID2\r\n",
"# BR NDA ID1 ID2 ID3\r\n",
" 2.63969810E-01 3 11 -12 14 # BR(Fe_2 -> Fe_1 Fv_1^* Fv_2 )\r\n",
" 7.36026659E-01 3 14 -2 1 # BR(Fe_2 -> Fv_2 Fu_1^* Fd_1 )\r\n",
" 3.53089158E-06 3 14 -2 3 # BR(Fe_2 -> Fv_2 Fu_1^* Fd_2 )\r\n",
"DECAY 15 2.27012221E-12 # Fe_3\r\n",
"# BR NDA ID1 ID2\r\n",
"# BR NDA ID1 ID2 ID3\r\n",
" 1.99665325E-01 3 11 -12 16 # BR(Fe_3 -> Fe_1 Fv_1^* Fv_3 )\r\n",
" 1.94367445E-01 3 13 -14 16 # BR(Fe_3 -> Fe_2 Fv_2^* Fv_3 )\r\n",
" 5.68539266E-01 3 16 -2 1 # BR(Fe_3 -> Fv_3 Fu_1^* Fd_1 )\r\n",
" 2.97768013E-02 3 16 -2 3 # BR(Fe_3 -> Fv_3 Fu_1^* Fd_2 )\r\n",
" 4.57228360E-04 3 16 -4 1 # BR(Fe_3 -> Fv_3 Fu_2^* Fd_1 )\r\n",
" 7.19393395E-03 3 16 -4 3 # BR(Fe_3 -> Fv_3 Fu_2^* Fd_2 )\r\n",
"DECAY 1 9.89732412E-27 # Fd_1\r\n",
"# BR NDA ID1 ID2\r\n",
"# BR NDA ID1 ID2 ID3\r\n",
" 1.00000000E+00 3 2 -12 11 # BR(Fd_1 -> Fu_1 Fv_1^* Fe_1 )\r\n",
"DECAY 3 3.77354480E-20 # Fd_2\r\n",
"# BR NDA ID1 ID2\r\n",
"# BR NDA ID1 ID2 ID3\r\n",
" 7.34967003E-01 3 1 -2 2 # BR(Fd_2 -> Fd_1 Fu_1^* Fu_1 )\r\n",
" 2.65032997E-01 3 2 -12 11 # BR(Fd_2 -> Fu_1 Fv_1^* Fe_1 )\r\n",
"DECAY 5 1.78283481E-13 # Fd_3\r\n",
"# BR NDA ID1 ID2\r\n",
"# BR NDA ID1 ID2 ID3\r\n",
" 6.12968728E-03 3 1 -2 2 # BR(Fd_3 -> Fd_1 Fu_1^* Fu_1 )\r\n",
" 4.52976684E-01 3 1 -2 4 # BR(Fd_3 -> Fd_1 Fu_1^* Fu_2 )\r\n",
" 1.67710056E-04 3 1 -4 2 # BR(Fd_3 -> Fd_1 Fu_2^* Fu_1 )\r\n",
" 8.80789452E-03 3 1 -4 4 # BR(Fd_3 -> Fd_1 Fu_2^* Fu_2 )\r\n",
" 3.27022113E-04 3 3 -2 2 # BR(Fd_3 -> Fd_2 Fu_1^* Fu_1 )\r\n",
" 2.41192782E-02 3 3 -2 4 # BR(Fd_3 -> Fd_2 Fu_1^* Fu_2 )\r\n",
" 3.11216734E-03 3 3 -4 2 # BR(Fd_3 -> Fd_2 Fu_2^* Fu_1 )\r\n",
" 1.62620924E-01 3 3 -4 4 # BR(Fd_3 -> Fd_2 Fu_2^* Fu_2 )\r\n",
" 2.15258744E-03 3 2 -12 11 # BR(Fd_3 -> Fu_1 Fv_1^* Fe_1 )\r\n",
" 2.14495777E-03 3 2 -14 13 # BR(Fd_3 -> Fu_1 Fv_2^* Fe_2 )\r\n",
" 5.83533347E-04 3 2 -16 15 # BR(Fd_3 -> Fu_1 Fv_3^* Fe_3 )\r\n",
" 1.59074417E-01 3 4 -12 11 # BR(Fd_3 -> Fu_2 Fv_1^* Fe_1 )\r\n",
" 1.58107801E-01 3 4 -14 13 # BR(Fd_3 -> Fu_2 Fv_2^* Fe_2 )\r\n",
" 1.96753358E-02 3 4 -16 15 # BR(Fd_3 -> Fu_2 Fv_3^* Fe_3 )\r\n",
"DECAY 25 9.23369614E-02 # hh_1\r\n",
"# BR NDA ID1 ID2\r\n",
" 1.31730794E-04 2 22 22 # BR(hh_1 -> VP VP )\r\n",
" 2.78510438E-01 2 21 21 # BR(hh_1 -> VG VG )\r\n",
" 9.85507188E-04 2 23 23 # BR(hh_1 -> VZ VZ )\r\n",
" 9.09858108E-03 2 24 -24 # BR(hh_1 -> VWm^* VWm_virt )\r\n",
" 6.37550446E-07 2 -1 1 # BR(hh_1 -> Fd_1^* Fd_1 )\r\n",
" 2.29321049E-04 2 -3 3 # BR(hh_1 -> Fd_2^* Fd_2 )\r\n",
" 6.13723954E-01 2 -5 5 # BR(hh_1 -> Fd_3^* Fd_3 )\r\n",
" 7.45829122E-09 2 -11 11 # BR(hh_1 -> Fe_1^* Fe_1 )\r\n",
" 3.33176957E-04 2 -13 13 # BR(hh_1 -> Fe_2^* Fe_2 )\r\n",
" 9.61718637E-02 2 -15 15 # BR(hh_1 -> Fe_3^* Fe_3 )\r\n",
" 3.43049426E-09 2 -2 2 # BR(hh_1 -> Fu_1^* Fu_1 )\r\n",
" 8.14778313E-04 2 -4 4 # BR(hh_1 -> Fu_2^* Fu_2 )\r\n",
"DECAY 35 1.73241775E+04 # hh_2\r\n",
"# BR NDA ID1 ID2\r\n",
" 1.39350046E-04 2 22 22 # BR(hh_2 -> VP VP )\r\n",
" 8.81932018E-01 2 21 21 # BR(hh_2 -> VG VG )\r\n",
" 1.82037910E-09 2 -1 1 # BR(hh_2 -> Fd_1^* Fd_1 )\r\n",
" 6.54773687E-07 2 -3 3 # BR(hh_2 -> Fd_2^* Fd_2 )\r\n",
" 1.75814213E-03 2 -5 5 # BR(hh_2 -> Fd_3^* Fd_3 )\r\n",
" 2.12954403E-11 2 -11 11 # BR(hh_2 -> Fe_1^* Fe_1 )\r\n",
" 9.51314459E-07 2 -13 13 # BR(hh_2 -> Fe_2^* Fe_2 )\r\n",
" 2.74928588E-04 2 -15 15 # BR(hh_2 -> Fe_3^* Fe_3 )\r\n",
" 8.73767246E-15 2 -2 2 # BR(hh_2 -> Fu_1^* Fu_1 )\r\n",
" 2.07572188E-09 2 -4 4 # BR(hh_2 -> Fu_2^* Fu_2 )\r\n",
" 1.60885997E-04 2 -6 6 # BR(hh_2 -> Fu_3^* Fu_3 )\r\n",
" 1.22390745E-05 2 -8 8 # BR(hh_2 -> Fu_4^* Fu_4 )\r\n",
" 1.06334869E-01 2 25 25 # BR(hh_2 -> hh_1 hh_1 )\r\n",
" 6.25679785E-03 2 -24 24 # BR(hh_2 -> VWm VWm^* )\r\n",
" 3.12915959E-03 2 23 23 # BR(hh_2 -> VZ VZ )\r\n",
"DECAY 36 1.64056318E+02 # Ah_2\r\n",
"# BR NDA ID1 ID2\r\n",
" 2.43798732E-09 2 22 22 # BR(Ah_2 -> VP VP )\r\n",
" 3.53881539E-06 2 21 21 # BR(Ah_2 -> VG VG )\r\n",
" 2.07718578E-07 2 -1 1 # BR(Ah_2 -> Fd_1^* Fd_1 )\r\n",
" 7.47144776E-05 2 -3 3 # BR(Ah_2 -> Fd_2^* Fd_2 )\r\n",
" 2.00617653E-01 2 -5 5 # BR(Ah_2 -> Fd_3^* Fd_3 )\r\n",
" 2.42996585E-09 2 -11 11 # BR(Ah_2 -> Fe_1^* Fe_1 )\r\n",
" 1.08551955E-04 2 -13 13 # BR(Ah_2 -> Fe_2^* Fe_2 )\r\n",
" 3.13714711E-02 2 -15 15 # BR(Ah_2 -> Fe_3^* Fe_3 )\r\n",
" 2.84863437E-14 2 -2 2 # BR(Ah_2 -> Fu_1^* Fu_1 )\r\n",
" 6.76721342E-09 2 -4 4 # BR(Ah_2 -> Fu_2^* Fu_2 )\r\n",
" 5.36133518E-04 2 -6 6 # BR(Ah_2 -> Fu_3^* Fu_3 )\r\n",
" 1.86969950E-04 2 -8 8 # BR(Ah_2 -> Fu_4^* Fu_4 )\r\n",
" 8.04062804E-04 2 -8 10 # BR(Ah_2 -> Fu_4^* Fu_5 )\r\n",
" 7.65200972E-01 2 25 23 # BR(Ah_2 -> hh_1 VZ )\r\n",
" 1.09571318E-03 2 35 23 # BR(Ah_2 -> hh_2 VZ )\r\n",
"DECAY 37 1.58744375E+02 # Hm_2\r\n",
"# BR NDA ID1 ID2\r\n",
" 1.83411169E-07 2 2 -1 # BR(Hm_2^* -> Fu_1 Fd_1^* )\r\n",
" 3.52942975E-06 2 2 -3 # BR(Hm_2^* -> Fu_1 Fd_2^* )\r\n",
" 2.18914690E-06 2 2 -5 # BR(Hm_2^* -> Fu_1 Fd_3^* )\r\n",
" 1.01273507E-08 2 4 -1 # BR(Hm_2^* -> Fu_2 Fd_1^* )\r\n",
" 6.58617967E-05 2 4 -3 # BR(Hm_2^* -> Fu_2 Fd_2^* )\r\n",
" 3.16535130E-04 2 4 -5 # BR(Hm_2^* -> Fu_2 Fd_3^* )\r\n",
" 1.75120588E-08 2 6 -1 # BR(Hm_2^* -> Fu_3 Fd_1^* )\r\n",
" 9.42509250E-07 2 6 -3 # BR(Hm_2^* -> Fu_3 Fd_2^* )\r\n",
" 1.84231336E-01 2 6 -5 # BR(Hm_2^* -> Fu_3 Fd_3^* )\r\n",
" 9.55295309E-05 2 8 -7 # BR(Hm_2^* -> Fu_4 Fd_4^* )\r\n",
" 1.57765671E-03 2 10 -7 # BR(Hm_2^* -> Fu_5 Fd_4^* )\r\n",
" 2.50019384E-09 2 12 -11 # BR(Hm_2^* -> Fv_1 Fe_1^* )\r\n",
" 1.11689195E-04 2 14 -13 # BR(Hm_2^* -> Fv_2 Fe_2^* )\r\n",
" 3.22781317E-02 2 16 -15 # BR(Hm_2^* -> Fv_3 Fe_3^* )\r\n",
" 7.80092186E-01 2 25 24 # BR(Hm_2^* -> hh_1 VWm^* )\r\n",
" 1.22419839E-03 2 35 24 # BR(Hm_2^* -> hh_2 VWm^* )\r\n"
]
}
],
"source": [
"cat SPheno.spc.4VL"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": []
}
],
"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.9"
}
},
"nbformat": 4,
"nbformat_minor": 0
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment