(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 9.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[      8054,        257]
NotebookOptionsPosition[      6950,        212]
NotebookOutlinePosition[      7306,        228]
CellTagsIndexPosition[      7263,        225]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["Setup", "Section",
 CellChangeTimes->{{3.5817002707890997`*^9, 3.581700271197218*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"$OldDir", "=", 
   RowBox[{"NotebookDirectory", "[", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"$FeynRulesPath", "=", 
   RowBox[{"SetDirectory", "[", "\"\<~/feynrules-current\>\"", "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"<<", "FeynRules`"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"LoadModel", "[", 
   RowBox[{"$FeynRulesPath", "<>", "\"\</Models/HEL.fr\>\""}], "]"}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.581700209879122*^9, 3.581700260936915*^9}, {
   3.581700376890725*^9, 3.58170037821945*^9}, {3.581702836763228*^9, 
   3.581702837338793*^9}, {3.5818481249015427`*^9, 3.581848143347954*^9}, {
   3.582995204405798*^9, 3.58299521528404*^9}, {3.583038702741356*^9, 
   3.583038709578618*^9}, 3.583076573840444*^9, {3.584182000911727*^9, 
   3.58418200121646*^9}, {3.584182161492667*^9, 3.584182163564756*^9}, {
   3.584182282715509*^9, 3.584182284105041*^9}, {3.584182399637722*^9, 
   3.584182404861476*^9}, 3.584182460107353*^9, {3.584182514313562*^9, 
   3.584182515436211*^9}, {3.5841828579996443`*^9, 3.584182858494486*^9}, {
   3.591369267589684*^9, 3.591369308504643*^9}, {3.6037025183316298`*^9, 
   3.603702520287414*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Strongly interacting Higgs", "Section",
 CellChangeTimes->{{3.592050781311898*^9, 3.592050792370556*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"Slagr", "=", 
   RowBox[{"LSILH", "+", "LSM"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Srules", " ", "=", " ", 
   RowBox[{"FeynmanRules", "[", "Slagr", "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"WriteUFO", "[", 
   RowBox[{"Slagr", ",", 
    RowBox[{"Output", "\[Rule]", "\"\<SILH_UFO\>\""}]}], "]"}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.59205079553372*^9, 3.592050838382057*^9}, {
  3.59205088788684*^9, 3.592050916405327*^9}, {3.592917055277775*^9, 
  3.59291706495752*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Pair of Higgs fields with a pair of leptons or quarks", "Section",
 CellChangeTimes->{{3.592050944690618*^9, 3.592050947760697*^9}, {
  3.592050990614394*^9, 3.5920509911924133`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"Plagr", "=", 
   RowBox[{"LF1", "+", "LSM"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Prules", " ", "=", " ", 
   RowBox[{"FeynmanRules", "[", "Plagr", "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"WriteUFO", "[", 
   RowBox[{"Plagr", ",", 
    RowBox[{"Output", "\[Rule]", "\"\<HHFF_UFO\>\""}]}], "]"}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.592050954724972*^9, 3.592050976242144*^9}, {
  3.592917046397978*^9, 3.592917051726028*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
Single Higgs field with a pair of leptons or quarks and a gauge boson \
\>", "Section",
 CellChangeTimes->{{3.592050989074213*^9, 3.5920509937940474`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"SHlagr", "=", 
   RowBox[{"LF2", "+", "LSM"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"SHrules", " ", "=", " ", 
   RowBox[{"FeynmanRules", "[", "SHlagr", "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"WriteUFO", "[", 
   RowBox[{"SHlagr", ",", 
    RowBox[{"Output", "\[Rule]", "\"\<HFFV_UFO\>\""}]}], "]"}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.5920509998861723`*^9, 3.592051016751165*^9}, {
  3.592917036374291*^9, 3.592917042654159*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Gauge bosons self couplings", "Section",
 CellChangeTimes->{{3.592051030955737*^9, 3.592051039361197*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"Glagr", "=", 
   RowBox[{"LBosons", "+", "LSM"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"Grules", " ", "=", " ", 
   RowBox[{"FeynmanRules", "[", "Glagr", "]"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"WriteUFO", "[", 
   RowBox[{"Glagr", ",", 
    RowBox[{"Output", "\[Rule]", "\"\<GaugeBosons_UFO\>\""}]}], "]"}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.592051040002797*^9, 3.592051070805373*^9}, {
  3.592917030382573*^9, 3.5929170336704893`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["CP-violating interactions", "Section",
 CellChangeTimes->{{3.592051082889255*^9, 3.5920510939613533`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"CPlagr", "=", 
   RowBox[{"LCP", "+", "LSM"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"CPrules", " ", "=", " ", 
   RowBox[{"FeynmanRules", "[", "CPlagr", "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"WriteUFO", "[", 
   RowBox[{"CPlagr", ",", 
    RowBox[{"Output", "\[Rule]", "\"\<CPviolating_UFO\>\""}]}], "]"}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.592051098237432*^9, 3.592051151133205*^9}, {
  3.592917021390731*^9, 3.5929170277827*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Gauge - Gauge - Higgs interactions", "Section",
 CellChangeTimes->{{3.592051370099566*^9, 3.592051380461054*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"GGlagr", "=", 
   RowBox[{"LLOPP", "+", "LSM"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"GGrules", " ", "=", " ", 
   RowBox[{"FeynmanRules", "[", "GGlagr", "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"WriteUFO", "[", 
   RowBox[{"GGlagr", ",", 
    RowBox[{"Output", "\[Rule]", "\"\<GGHiggs_UFO\>\""}]}], "]"}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.592051387684909*^9, 3.592051410737917*^9}, {
  3.59291701185524*^9, 3.592917016335073*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Full Lagrangian", "Section",
 CellChangeTimes->{{3.5817009721675787`*^9, 3.5817009738937607`*^9}, {
  3.591369347439767*^9, 3.591369348156604*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"lagr", "=", "LagHEL"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"feynmanrules", "=", 
   RowBox[{"FeynmanRules", "[", 
    RowBox[{"lagr", ",", 
     RowBox[{"ScreenOutput", "\[Rule]", "False"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"WriteUFO", "[", 
   RowBox[{"lagr", ",", 
    RowBox[{"ScreenOutput", "\[Rule]", "False"}]}], "]"}], ";"}]}], "Input",
 CellChangeTimes->{{3.581848751644443*^9, 3.5818487554853983`*^9}, {
  3.591369405252789*^9, 3.591369406748837*^9}, {3.592051417774378*^9, 
  3.5920514331009893`*^9}, {3.592917000529808*^9, 3.592917003555757*^9}, {
  3.603702543963779*^9, 3.6037025584389343`*^9}}]
}, Open  ]]
},
WindowSize->{1280, 976},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
FrontEndVersion->"9.0 for Mac OS X x86 (32-bit, 64-bit Kernel) (November 20, \
2012)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[579, 22, 92, 1, 80, "Section"],
Cell[674, 25, 1232, 24, 80, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[1943, 54, 111, 1, 80, "Section"],
Cell[2057, 57, 556, 15, 63, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[2650, 77, 189, 2, 80, "Section"],
Cell[2842, 81, 508, 14, 63, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[3387, 100, 165, 3, 80, "Section"],
Cell[3555, 105, 514, 14, 63, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[4106, 124, 112, 1, 80, "Section"],
Cell[4221, 127, 484, 13, 63, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[4742, 145, 112, 1, 80, "Section"],
Cell[4857, 148, 517, 14, 63, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[5411, 167, 119, 1, 80, "Section"],
Cell[5533, 170, 516, 14, 63, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[6086, 189, 153, 2, 80, "Section"],
Cell[6242, 193, 692, 16, 63, "Input"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)
