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

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

(* CreatedBy='Mathematica 8.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[     14026,        382]
NotebookOptionsPosition[     12033,        310]
NotebookOutlinePosition[     12445,        326]
CellTagsIndexPosition[     12402,        323]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[{
 StyleBox[
  RowBox[{
  "Computer", " ", "tools", " ", "in", " ", "particle", " ", "physics"}], 
  "Title"], "\[IndentingNewLine]", 
 StyleBox[
  RowBox[{
   RowBox[{"Course", " ", "by", " ", "Avelino", " ", "Vicente"}], " ", "-", 
   " ", "IFIC", " ", "-", " ", 
   RowBox[{"May", " ", "2019"}]}], "Subtitle"]}], "Input",
 Evaluatable->False,
 CellChangeTimes->CompressedData["
1:eJxTTMoPSmViYGAQAWIQfWm3WF1C+WtH8f8FjSA68WhJF4jWOz1/Gog2e/li
HojO2Gy/AERb873aBBbfsWsziD73O/ssiK7z3gamJRRMpk0G0hcSzMB0esEE
Ff7W146rHraqgugvjn5/HIB02JdWRkcgzXCq0H1n22vHFYxVHiC66PwGta2d
rx33zrHTBNGbdByqDbtfO56qdALTBh8DXuf3A+0NfACmAzfMzyoG0utOyWWD
6KxTGz7OmvDaMeTbLjA9aWPlL7GJrx015kv9A9FxCwz1xYH0+Uh3ML2rxYxT
YtJrx9AlF8D0itIls0M+vHZ88uIImHZza5pTLPTGUe9B8DwQ/WjN+UUgOixk
/mIQDQBAWadN
  "]],

Cell[TextData[{
 StyleBox["With this notebook we will learn the most common commands in ",
  FontSize->18,
  FontColor->GrayLevel[0],
  Background->None],
 StyleBox["SARAH",
  FontSize->18,
  FontWeight->"Bold",
  FontColor->GrayLevel[0],
  Background->None]
}], "Text",
 CellChangeTimes->{{3.766986458621125*^9, 3.766986480811357*^9}, {
  3.766986572646778*^9, 3.766986759367234*^9}, {3.7669871627960043`*^9, 
  3.766987165578465*^9}, {3.766998870952428*^9, 3.766998871156941*^9}, {
  3.767770408031024*^9, 3.767770416528718*^9}},
 Background->RGBColor[1, 1, 0.85]],

Cell[CellGroupData[{

Cell["Load SARAH", "Title",
 Evaluatable->False,
 CellChangeTimes->{{3.686467940463894*^9, 3.686467960453517*^9}, 
   3.686469580779161*^9, {3.6936400185291033`*^9, 3.693640019013015*^9}, {
   3.6990824216357803`*^9, 3.699082443019319*^9}, {3.699085905044036*^9, 
   3.6990859076625643`*^9}, {3.7677704276155*^9, 3.767770430228077*^9}}],

Cell[BoxData[
 RowBox[{"Get", "[", "\"\<codes/SARAH-4.14.1/SARAH.m\>\"", "]"}]], "Input",
 CellChangeTimes->{{3.543066433311716*^9, 3.543066444646021*^9}, 
   3.612673160960286*^9, {3.642139187375333*^9, 3.642139200640821*^9}, {
   3.642159187821062*^9, 3.6421591886540318`*^9}, {3.6422650341114063`*^9, 
   3.642265034239678*^9}, 3.6436190288853693`*^9, {3.657015732584651*^9, 
   3.657015732966357*^9}, {3.6628079285347633`*^9, 3.662807931293388*^9}, {
   3.662809057105275*^9, 3.6628090574333763`*^9}, 3.663315255740694*^9, {
   3.6634129299059753`*^9, 3.6634129305696487`*^9}, {3.663499193352339*^9, 
   3.663499193939356*^9}, {3.68750482518433*^9, 3.687504825664859*^9}, {
   3.687781720909919*^9, 3.68778172276828*^9}, {3.703999956135502*^9, 
   3.7039999615838833`*^9}, {3.767770713967263*^9, 3.7677707155130653`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Load Scotogenic model", "Title",
 Evaluatable->False,
 CellChangeTimes->{{3.686467940463894*^9, 3.686467960453517*^9}, 
   3.686469580779161*^9, {3.6936400185291033`*^9, 3.693640019013015*^9}, {
   3.6990824216357803`*^9, 3.699082443019319*^9}, {3.699085905044036*^9, 
   3.6990859076625643`*^9}, {3.7677704276155*^9, 3.7677704492640553`*^9}}],

Cell[BoxData[
 RowBox[{"Start", "[", "\"\<Scotogenic\>\"", "]"}]], "Input",
 CellChangeTimes->{{3.61769587968014*^9, 3.617695885302827*^9}, {
   3.617700031416925*^9, 3.617700043583078*^9}, {3.617700994605487*^9, 
   3.617701002081102*^9}, {3.617703853456403*^9, 3.6177038578200903`*^9}, {
   3.617704011996723*^9, 3.617704019516409*^9}, {3.6177729141097927`*^9, 
   3.617772916719343*^9}, {3.617772975360723*^9, 3.617772977286265*^9}, 
   3.641892711213243*^9, 3.642139344190316*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Explore the model", "Title",
 Evaluatable->False,
 CellChangeTimes->{{3.686467940463894*^9, 3.686467960453517*^9}, 
   3.686469580779161*^9, {3.6936400185291033`*^9, 3.693640019013015*^9}, {
   3.6990824216357803`*^9, 3.699082443019319*^9}, {3.699085905044036*^9, 
   3.6990859076625643`*^9}, {3.7677704276155*^9, 3.767770463041563*^9}}],

Cell[CellGroupData[{

Cell["Mass matrices", "Subchapter",
 CellChangeTimes->{{3.767770465346221*^9, 3.7677704682569523`*^9}}],

Cell[BoxData[
 RowBox[{"MassMatrix", "[", "Fe", "]"}]], "Input",
 CellChangeTimes->{{3.6418944585737343`*^9, 3.6418944618362913`*^9}, {
   3.70400038925079*^9, 3.704000391798359*^9}, 3.7040007082282763`*^9}],

Cell[BoxData[
 RowBox[{"MassMatrix", "[", "Chi", "]"}]], "Input",
 CellChangeTimes->{{3.641894757416142*^9, 3.641894757714943*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"MassMatrix", "[", "Chi", "]"}], "/.", 
  RowBox[{
   RowBox[{"Mn", "[", 
    RowBox[{"i_", ",", "j_"}], "]"}], "\[RuleDelayed]", 
   RowBox[{"If", "[", 
    RowBox[{
     RowBox[{"i", ">", "j"}], ",", 
     RowBox[{"Mn", "[", 
      RowBox[{"j", ",", "i"}], "]"}], ",", 
     RowBox[{"Mn", "[", 
      RowBox[{"i", ",", "j"}], "]"}]}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.641895406755147*^9, 3.641895434136468*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Mass", "[", "hh", "]"}], "/.", 
  RowBox[{"Masses", "[", "EWSB", "]"}]}]], "Input",
 CellChangeTimes->{{3.641893857192151*^9, 3.64189385976862*^9}, {
  3.641894277612791*^9, 3.6418942855435877`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Tadpole equations", "Subchapter",
 CellChangeTimes->{{3.76777047678607*^9, 3.767770480429604*^9}}],

Cell[BoxData[
 RowBox[{"TadpoleEquations", "[", "EWSB", "]"}]], "Input",
 CellChangeTimes->{{3.641892918846278*^9, 3.641892924317926*^9}}],

Cell[BoxData[
 RowBox[{"TadpoleEquation", "[", "v", "]"}]], "Input",
 CellChangeTimes->{{3.641892939357213*^9, 3.6418929503732023`*^9}}],

Cell[BoxData[
 RowBox[{"solTadpole", "=", 
  RowBox[{"Solve", "[", 
   RowBox[{
    RowBox[{"TadpoleEquation", "[", "v", "]"}], ",", "mH2"}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.6418932816192017`*^9, 3.641893289530778*^9}, {
  3.6418959581884813`*^9, 3.641895960119982*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"Mass", "[", "hh", "]"}], "/.", 
   RowBox[{"Masses", "[", "EWSB", "]"}]}], "/.", "solTadpole"}]], "Input",
 CellChangeTimes->{{3.641895966681736*^9, 3.641895969857542*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Vertices", "Subchapter",
 CellChangeTimes->{{3.7677704882119617`*^9, 3.767770491426187*^9}}],

Cell[BoxData[
 RowBox[{"Vertex", "[", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"bar", "[", "Fe", "]"}], ",", "Fe", ",", "hh"}], "}"}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.6418966032427073`*^9, 3.64189661309589*^9}}],

Cell[BoxData[
 RowBox[{"Vertex", "[", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"bar", "[", "Fe", "]"}], ",", "Fv", ",", 
    RowBox[{"conj", "[", "VWp", "]"}]}], "}"}], "]"}]], "Input",
 CellChangeTimes->{{3.641999275206624*^9, 3.641999283070676*^9}, {
  3.641999427941495*^9, 3.641999439992008*^9}}],

Cell[BoxData[
 RowBox[{"Vertex", "[", 
  RowBox[{"{", 
   RowBox[{"Fv", ",", "Chi", ",", "hh"}], "}"}], "]"}]], "Input",
 CellChangeTimes->{{3.641896631893976*^9, 3.641896632179007*^9}, {
  3.6418968663984213`*^9, 3.6418968675826683`*^9}, {3.6418968985978413`*^9, 
  3.641896898863695*^9}}],

Cell[BoxData[
 RowBox[{"Vertex", "[", 
  RowBox[{"{", 
   RowBox[{"Fv", ",", "Chi", ",", "etR"}], "}"}], "]"}]], "Input",
 CellChangeTimes->{{3.641896874249325*^9, 3.64189688850025*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["LaTeX output", "Subchapter",
 CellChangeTimes->{{3.7677705006286783`*^9, 3.767770505667264*^9}}],

Cell[BoxData[
 RowBox[{"ModelOutput", "[", "EWSB", "]"}]], "Input",
 CellChangeTimes->{{3.6633157695125732`*^9, 3.663315774439954*^9}}],

Cell[BoxData[
 RowBox[{"MakeTeX", "[", "]"}]], "Input",
 CellChangeTimes->CompressedData["
1:eJxTTMoPSmViYGCQBGIQbST8QbLz0yvHx2ofZEF0WKjh/JlAmnl6OZiefvy2
XOvXV463Hv8D00nKjdwa3145pi3ZAKZXafOetvrxyvHQVHswnfWT+wOINq52
BdOBorbb9jG+dnRenQ6mY6fzqB8E0t9fmIHpj9cWTN/F/Nqx78QxMN1y28lj
H5CWuFULphv4l6+6DaT/pV0F0zWCHXtA9Bm3TWB6D1dMtrXsa0ddwyVgusOw
deoH+deOx+58AdN71iVlfwTS6z70gGmbeeWXeRReO36o2Aimdz80mMKk9Nox
uD0UTE9zbmFnAdIRRevA9Maud8YgOm+2tAmIdpNp04sB0gvstoHpTQbCkteA
tJOEPZieJ7bC6BGQDtp6G0xvuxHCJKT82lGArRJMv12lGwCiNSvCwHTmox2L
tqm+dsw//g9MAwCmJ8OU
  "]]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Input for other codes", "Title",
 Evaluatable->False,
 CellChangeTimes->{{3.686467940463894*^9, 3.686467960453517*^9}, 
   3.686469580779161*^9, {3.6936400185291033`*^9, 3.693640019013015*^9}, {
   3.6990824216357803`*^9, 3.699082443019319*^9}, {3.699085905044036*^9, 
   3.6990859076625643`*^9}, {3.7677704276155*^9, 3.767770463041563*^9}, {
   3.767770516406705*^9, 3.767770519100067*^9}}],

Cell[CellGroupData[{

Cell["SPheno module", "Subchapter",
 CellChangeTimes->{{3.767770521418839*^9, 3.767770529341846*^9}}],

Cell[BoxData[
 RowBox[{"MakeSPheno", "[", 
  RowBox[{"IncludeLoopDecays", "\[Rule]", "False"}], "]"}]], "Input",
 CellChangeTimes->{{3.613931911833357*^9, 3.613931920588901*^9}, {
   3.6167554365343027`*^9, 3.616755439897934*^9}, {3.6168231213917027`*^9, 
   3.6168231371752*^9}, {3.61684007729959*^9, 3.616840081183877*^9}, {
   3.616853246013549*^9, 3.6168532469611177`*^9}, {3.617537504856537*^9, 
   3.617537508304668*^9}, 3.6175972471865063`*^9, {3.617697143417962*^9, 
   3.617697146813134*^9}, {3.618052486143476*^9, 3.618052486584504*^9}, {
   3.618381013831895*^9, 3.618381015912339*^9}, {3.618384099383676*^9, 
   3.61838410119975*^9}, {3.618384349960568*^9, 3.618384359990724*^9}, {
   3.6187177917303753`*^9, 3.618717793059134*^9}, {3.6188234229210577`*^9, 
   3.6188234238916063`*^9}, {3.61978398780939*^9, 3.619783989833467*^9}, {
   3.61980506094975*^9, 3.6198050619234877`*^9}, {3.621619305737402*^9, 
   3.621619308978161*^9}, {3.6216730527531023`*^9, 3.6216730538153067`*^9}, {
   3.621674593213233*^9, 3.621674594876628*^9}, {3.621677598277638*^9, 
   3.621677599550735*^9}, {3.622630329106503*^9, 3.622630330943138*^9}, {
   3.623064749056677*^9, 3.623064751218445*^9}, {3.62307087401628*^9, 
   3.623070875501781*^9}, {3.6258336473923492`*^9, 3.625833664998884*^9}, {
   3.6258412251011047`*^9, 3.625841226453388*^9}, {3.626157473988639*^9, 
   3.6261574753543*^9}, {3.627220633033099*^9, 3.627220634267314*^9}, {
   3.627886533873753*^9, 3.627886535778747*^9}, {3.627890646251831*^9, 
   3.627890647430932*^9}, {3.633328961168841*^9, 3.633328963121683*^9}, {
   3.6333297964098063`*^9, 3.633329797406958*^9}, {3.633333083709848*^9, 
   3.633333084937089*^9}, {3.640625974200976*^9, 3.6406259762394533`*^9}, {
   3.640626006524796*^9, 3.640626008041626*^9}, {3.641270431827936*^9, 
   3.641270432996127*^9}, {3.64127432912197*^9, 3.641274355130992*^9}, {
   3.6421676777786427`*^9, 3.642167680099959*^9}, {3.642169615965905*^9, 
   3.642169617676506*^9}, {3.6422662087640123`*^9, 3.64226620981063*^9}, {
   3.642266518638473*^9, 3.642266521006317*^9}, {3.64365135427291*^9, 
   3.643651355232779*^9}, {3.704035642003871*^9, 3.704035648350484*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Input file for micrOmegas", "Subchapter",
 CellChangeTimes->{{3.767770535038866*^9, 3.7677705416248837`*^9}}],

Cell[BoxData[
 RowBox[{"MakeCHep", "[", "]"}]], "Input",
 CellChangeTimes->{{3.642266468948337*^9, 3.642266497973543*^9}, {
  3.642789021753961*^9, 3.64278902288542*^9}, {3.643548904880966*^9, 
  3.6435489063211527`*^9}, {3.6436513564528093`*^9, 3.643651357391623*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Input UFO files for MadGraph", "Subchapter",
 CellChangeTimes->{{3.767770545640856*^9, 3.767770552724971*^9}}],

Cell[BoxData[
 RowBox[{"MakeUFO", "[", "]"}]], "Input",
 CellChangeTimes->{{3.6421679473240767`*^9, 3.642167948801838*^9}, {
  3.64216960516776*^9, 3.642169606837376*^9}, {3.642266444246963*^9, 
  3.642266451160055*^9}, {3.643548908486803*^9, 3.643548909672386*^9}}]
}, Open  ]]
}, Open  ]]
},
WindowSize->{1280, 816},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}},
FrontEndVersion->"10.4 for Linux x86 (64-bit) (April 11, 2016)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[557, 20, 780, 19, 91, "Input",
 Evaluatable->False],
Cell[1340, 41, 566, 15, 54, "Text"],
Cell[CellGroupData[{
Cell[1931, 60, 336, 5, 93, "Title",
 Evaluatable->False],
Cell[2270, 67, 825, 11, 32, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[3132, 83, 349, 5, 93, "Title",
 Evaluatable->False],
Cell[3484, 90, 484, 7, 32, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[4005, 102, 343, 5, 93, "Title",
 Evaluatable->False],
Cell[CellGroupData[{
Cell[4373, 111, 103, 1, 65, "Subchapter"],
Cell[4479, 114, 207, 3, 32, "Input"],
Cell[4689, 119, 131, 2, 32, "Input"],
Cell[4823, 123, 457, 13, 32, "Input"],
Cell[5283, 138, 234, 5, 32, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[5554, 148, 104, 1, 65, "Subchapter"],
Cell[5661, 151, 138, 2, 32, "Input"],
Cell[5802, 155, 136, 2, 32, "Input"],
Cell[5941, 159, 283, 7, 32, "Input"],
Cell[6227, 168, 219, 5, 32, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[6483, 178, 98, 1, 65, "Subchapter"],
Cell[6584, 181, 221, 6, 32, "Input"],
Cell[6808, 189, 300, 7, 32, "Input"],
Cell[7111, 198, 290, 6, 32, "Input"],
Cell[7404, 206, 186, 4, 32, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[7627, 215, 102, 1, 65, "Subchapter"],
Cell[7732, 218, 135, 2, 32, "Input"],
Cell[7870, 222, 546, 11, 32, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[8465, 239, 397, 6, 93, "Title",
 Evaluatable->False],
Cell[CellGroupData[{
Cell[8887, 249, 101, 1, 65, "Subchapter"],
Cell[8991, 252, 2168, 30, 32, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[11196, 287, 115, 1, 65, "Subchapter"],
Cell[11314, 290, 269, 4, 32, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[11620, 299, 116, 1, 65, "Subchapter"],
Cell[11739, 302, 266, 4, 32, "Input"]
}, Open  ]]
}, Open  ]]
}
]
*)