(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[         0,          0]
NotebookDataLength[    121719,       2992]
NotebookOptionsPosition[    114005,       2750]
NotebookOutlinePosition[    114557,       2770]
CellTagsIndexPosition[    114514,       2767]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[" Tutorial", "Title",
 CellChangeTimes->{{3.539317450114971*^9, 3.5393174577917767`*^9}, 
   3.5559442666324883`*^9}],

Cell[BoxData[
 RowBox[{"Quit", "[", "]"}]], "Input",
 CellChangeTimes->{{3.5393174647226133`*^9, 3.5393174655928297`*^9}}],

Cell[CellGroupData[{

Cell["Loading FeynRules", "Section",
 CellChangeTimes->{{3.539317474608481*^9, 3.5393174784154987`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"$FeynRulesPath", " ", "=", " ", 
  RowBox[{
  "SetDirectory", "[", "\"\<~/feynrules/branch/feynrules-current\>\"", 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.5393174792335243`*^9, 3.539317505486512*^9}, {
  3.5780999932057743`*^9, 3.578099993939001*^9}, {3.578100105316533*^9, 
  3.5781001104184732`*^9}, {3.579349162263041*^9, 3.5793491976568537`*^9}, {
  3.579349302055018*^9, 3.57934930213306*^9}, {3.614226881066931*^9, 
  3.614226898369299*^9}, {3.6178016287606277`*^9, 3.6178016466541862`*^9}, {
  3.677924476563367*^9, 3.6779244836520643`*^9}}],

Cell[BoxData["\<\"/Users/degrande/feynrules/branch/feynrules-current\"\>"], \
"Output",
 CellChangeTimes->{
  3.5393175059082603`*^9, 3.539317637249185*^9, 3.539317710216969*^9, 
   3.53931836320551*^9, 3.539318857008144*^9, {3.555944274657331*^9, 
   3.5559442982248707`*^9}, 3.5559443774661837`*^9, 3.5559446049920673`*^9, 
   3.555944681188205*^9, 3.577373428262697*^9, 3.577374281796172*^9, {
   3.578100112630537*^9, 3.5781001499397497`*^9}, 3.57934927304145*^9, 
   3.579349304490674*^9, 3.579350283582839*^9, 3.579350512924067*^9, 
   3.579350634423061*^9, 3.579358789644754*^9, {3.57935921014121*^9, 
   3.5793592370286303`*^9}, 3.579359528417985*^9, 3.57936106637228*^9, 
   3.579361288219393*^9, 3.579882686241864*^9, 3.579882984714159*^9, 
   3.579955782353389*^9, 3.614226900065674*^9, 3.614227196689986*^9, 
   3.617162109487749*^9, 3.617801329310379*^9, 3.61780165232563*^9, 
   3.617813329932769*^9, 3.677924495713375*^9, 3.677926813709795*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"<<", "FeynRules`"}]], "Input",
 CellChangeTimes->{{3.5393175083150997`*^9, 3.5393175100537786`*^9}}],

Cell[CellGroupData[{

Cell[BoxData["\<\" - FeynRules - \"\>"], "Print",
 CellChangeTimes->{
  3.53931751083799*^9, 3.539317637583722*^9, 3.539317710401309*^9, 
   3.539318367080615*^9, 3.539318857563046*^9, {3.5559442763044662`*^9, 
   3.555944298931198*^9}, 3.555944377919302*^9, 3.5559446052682056`*^9, 
   3.555944681456799*^9, 3.577373391273698*^9, {3.577373428852648*^9, 
   3.5773734513303223`*^9}, 3.5773742883918943`*^9, {3.578100113697337*^9, 
   3.5781001512116632`*^9}, 3.579349307063857*^9, 3.579350284229608*^9, 
   3.579350513044907*^9, 3.579350634908613*^9, 3.579358789778343*^9, {
   3.579359210417015*^9, 3.579359238091683*^9}, 3.5793595291764708`*^9, 
   3.579361067924972*^9, 3.579361288919224*^9, 3.579882686436294*^9, 
   3.579882985016074*^9, 3.5799557837408743`*^9, 3.6142269017242403`*^9, 
   3.614227200332449*^9, 3.617162111297064*^9, 3.617801329548204*^9, 
   3.61780165248982*^9, 3.617813330217751*^9, 3.677924497275835*^9, 
   3.6779268138730087`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Version: \"\>", "\[InvisibleSpace]", "\<\"2.3.22\"\>", 
   "\[InvisibleSpace]", 
   RowBox[{"\<\" (\"\>", " ", "\<\"04 May 2016\"\>"}], 
   "\[InvisibleSpace]", "\<\").\"\>"}],
  SequenceForm["Version: ", "2.3.22", " (" "04 May 2016", ")."],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.53931751083799*^9, 3.539317637583722*^9, 3.539317710401309*^9, 
   3.539318367080615*^9, 3.539318857563046*^9, {3.5559442763044662`*^9, 
   3.555944298931198*^9}, 3.555944377919302*^9, 3.5559446052682056`*^9, 
   3.555944681456799*^9, 3.577373391273698*^9, {3.577373428852648*^9, 
   3.5773734513303223`*^9}, 3.5773742883918943`*^9, {3.578100113697337*^9, 
   3.5781001512116632`*^9}, 3.579349307063857*^9, 3.579350284229608*^9, 
   3.579350513044907*^9, 3.579350634908613*^9, 3.579358789778343*^9, {
   3.579359210417015*^9, 3.579359238091683*^9}, 3.5793595291764708`*^9, 
   3.579361067924972*^9, 3.579361288919224*^9, 3.579882686436294*^9, 
   3.579882985016074*^9, 3.5799557837408743`*^9, 3.6142269017242403`*^9, 
   3.614227200332449*^9, 3.617162111297064*^9, 3.617801329548204*^9, 
   3.61780165248982*^9, 3.617813330217751*^9, 3.677924497275835*^9, 
   3.677926813880662*^9}],

Cell[BoxData["\<\"Authors: A. Alloul, N. Christensen, C. Degrande, C. Duhr, \
B. Fuks\"\>"], "Print",
 CellChangeTimes->{
  3.53931751083799*^9, 3.539317637583722*^9, 3.539317710401309*^9, 
   3.539318367080615*^9, 3.539318857563046*^9, {3.5559442763044662`*^9, 
   3.555944298931198*^9}, 3.555944377919302*^9, 3.5559446052682056`*^9, 
   3.555944681456799*^9, 3.577373391273698*^9, {3.577373428852648*^9, 
   3.5773734513303223`*^9}, 3.5773742883918943`*^9, {3.578100113697337*^9, 
   3.5781001512116632`*^9}, 3.579349307063857*^9, 3.579350284229608*^9, 
   3.579350513044907*^9, 3.579350634908613*^9, 3.579358789778343*^9, {
   3.579359210417015*^9, 3.579359238091683*^9}, 3.5793595291764708`*^9, 
   3.579361067924972*^9, 3.579361288919224*^9, 3.579882686436294*^9, 
   3.579882985016074*^9, 3.5799557837408743`*^9, 3.6142269017242403`*^9, 
   3.614227200332449*^9, 3.617162111297064*^9, 3.617801329548204*^9, 
   3.61780165248982*^9, 3.617813330217751*^9, 3.677924497275835*^9, 
   3.6779268138872623`*^9}],

Cell[BoxData["\<\" \"\>"], "Print",
 CellChangeTimes->{
  3.53931751083799*^9, 3.539317637583722*^9, 3.539317710401309*^9, 
   3.539318367080615*^9, 3.539318857563046*^9, {3.5559442763044662`*^9, 
   3.555944298931198*^9}, 3.555944377919302*^9, 3.5559446052682056`*^9, 
   3.555944681456799*^9, 3.577373391273698*^9, {3.577373428852648*^9, 
   3.5773734513303223`*^9}, 3.5773742883918943`*^9, {3.578100113697337*^9, 
   3.5781001512116632`*^9}, 3.579349307063857*^9, 3.579350284229608*^9, 
   3.579350513044907*^9, 3.579350634908613*^9, 3.579358789778343*^9, {
   3.579359210417015*^9, 3.579359238091683*^9}, 3.5793595291764708`*^9, 
   3.579361067924972*^9, 3.579361288919224*^9, 3.579882686436294*^9, 
   3.579882985016074*^9, 3.5799557837408743`*^9, 3.6142269017242403`*^9, 
   3.614227200332449*^9, 3.617162111297064*^9, 3.617801329548204*^9, 
   3.61780165248982*^9, 3.617813330217751*^9, 3.677924497275835*^9, 
   3.677926813891303*^9}],

Cell[BoxData["\<\"Please cite:\"\>"], "Print",
 CellChangeTimes->{
  3.53931751083799*^9, 3.539317637583722*^9, 3.539317710401309*^9, 
   3.539318367080615*^9, 3.539318857563046*^9, {3.5559442763044662`*^9, 
   3.555944298931198*^9}, 3.555944377919302*^9, 3.5559446052682056`*^9, 
   3.555944681456799*^9, 3.577373391273698*^9, {3.577373428852648*^9, 
   3.5773734513303223`*^9}, 3.5773742883918943`*^9, {3.578100113697337*^9, 
   3.5781001512116632`*^9}, 3.579349307063857*^9, 3.579350284229608*^9, 
   3.579350513044907*^9, 3.579350634908613*^9, 3.579358789778343*^9, {
   3.579359210417015*^9, 3.579359238091683*^9}, 3.5793595291764708`*^9, 
   3.579361067924972*^9, 3.579361288919224*^9, 3.579882686436294*^9, 
   3.579882985016074*^9, 3.5799557837408743`*^9, 3.6142269017242403`*^9, 
   3.614227200332449*^9, 3.617162111297064*^9, 3.617801329548204*^9, 
   3.61780165248982*^9, 3.617813330217751*^9, 3.677924497275835*^9, 
   3.677926813894746*^9}],

Cell[BoxData["\<\"    - Comput.Phys.Commun.185:2250-2300,2014 \
(arXiv:1310.1921);\"\>"], "Print",
 CellChangeTimes->{
  3.53931751083799*^9, 3.539317637583722*^9, 3.539317710401309*^9, 
   3.539318367080615*^9, 3.539318857563046*^9, {3.5559442763044662`*^9, 
   3.555944298931198*^9}, 3.555944377919302*^9, 3.5559446052682056`*^9, 
   3.555944681456799*^9, 3.577373391273698*^9, {3.577373428852648*^9, 
   3.5773734513303223`*^9}, 3.5773742883918943`*^9, {3.578100113697337*^9, 
   3.5781001512116632`*^9}, 3.579349307063857*^9, 3.579350284229608*^9, 
   3.579350513044907*^9, 3.579350634908613*^9, 3.579358789778343*^9, {
   3.579359210417015*^9, 3.579359238091683*^9}, 3.5793595291764708`*^9, 
   3.579361067924972*^9, 3.579361288919224*^9, 3.579882686436294*^9, 
   3.579882985016074*^9, 3.5799557837408743`*^9, 3.6142269017242403`*^9, 
   3.614227200332449*^9, 3.617162111297064*^9, 3.617801329548204*^9, 
   3.61780165248982*^9, 3.617813330217751*^9, 3.677924497275835*^9, 
   3.677926813898588*^9}],

Cell[BoxData["\<\"    - Comput.Phys.Commun.180:1614-1641,2009 \
(arXiv:0806.4194).\"\>"], "Print",
 CellChangeTimes->{
  3.53931751083799*^9, 3.539317637583722*^9, 3.539317710401309*^9, 
   3.539318367080615*^9, 3.539318857563046*^9, {3.5559442763044662`*^9, 
   3.555944298931198*^9}, 3.555944377919302*^9, 3.5559446052682056`*^9, 
   3.555944681456799*^9, 3.577373391273698*^9, {3.577373428852648*^9, 
   3.5773734513303223`*^9}, 3.5773742883918943`*^9, {3.578100113697337*^9, 
   3.5781001512116632`*^9}, 3.579349307063857*^9, 3.579350284229608*^9, 
   3.579350513044907*^9, 3.579350634908613*^9, 3.579358789778343*^9, {
   3.579359210417015*^9, 3.579359238091683*^9}, 3.5793595291764708`*^9, 
   3.579361067924972*^9, 3.579361288919224*^9, 3.579882686436294*^9, 
   3.579882985016074*^9, 3.5799557837408743`*^9, 3.6142269017242403`*^9, 
   3.614227200332449*^9, 3.617162111297064*^9, 3.617801329548204*^9, 
   3.61780165248982*^9, 3.617813330217751*^9, 3.677924497275835*^9, 
   3.67792681390261*^9}],

Cell[BoxData["\<\" \"\>"], "Print",
 CellChangeTimes->{
  3.53931751083799*^9, 3.539317637583722*^9, 3.539317710401309*^9, 
   3.539318367080615*^9, 3.539318857563046*^9, {3.5559442763044662`*^9, 
   3.555944298931198*^9}, 3.555944377919302*^9, 3.5559446052682056`*^9, 
   3.555944681456799*^9, 3.577373391273698*^9, {3.577373428852648*^9, 
   3.5773734513303223`*^9}, 3.5773742883918943`*^9, {3.578100113697337*^9, 
   3.5781001512116632`*^9}, 3.579349307063857*^9, 3.579350284229608*^9, 
   3.579350513044907*^9, 3.579350634908613*^9, 3.579358789778343*^9, {
   3.579359210417015*^9, 3.579359238091683*^9}, 3.5793595291764708`*^9, 
   3.579361067924972*^9, 3.579361288919224*^9, 3.579882686436294*^9, 
   3.579882985016074*^9, 3.5799557837408743`*^9, 3.6142269017242403`*^9, 
   3.614227200332449*^9, 3.617162111297064*^9, 3.617801329548204*^9, 
   3.61780165248982*^9, 3.617813330217751*^9, 3.677924497275835*^9, 
   3.6779268139063377`*^9}],

Cell[BoxData["\<\"http://feynrules.phys.ucl.ac.be\"\>"], "Print",
 CellChangeTimes->{
  3.53931751083799*^9, 3.539317637583722*^9, 3.539317710401309*^9, 
   3.539318367080615*^9, 3.539318857563046*^9, {3.5559442763044662`*^9, 
   3.555944298931198*^9}, 3.555944377919302*^9, 3.5559446052682056`*^9, 
   3.555944681456799*^9, 3.577373391273698*^9, {3.577373428852648*^9, 
   3.5773734513303223`*^9}, 3.5773742883918943`*^9, {3.578100113697337*^9, 
   3.5781001512116632`*^9}, 3.579349307063857*^9, 3.579350284229608*^9, 
   3.579350513044907*^9, 3.579350634908613*^9, 3.579358789778343*^9, {
   3.579359210417015*^9, 3.579359238091683*^9}, 3.5793595291764708`*^9, 
   3.579361067924972*^9, 3.579361288919224*^9, 3.579882686436294*^9, 
   3.579882985016074*^9, 3.5799557837408743`*^9, 3.6142269017242403`*^9, 
   3.614227200332449*^9, 3.617162111297064*^9, 3.617801329548204*^9, 
   3.61780165248982*^9, 3.617813330217751*^9, 3.677924497275835*^9, 
   3.677926813910321*^9}],

Cell[BoxData["\<\" \"\>"], "Print",
 CellChangeTimes->{
  3.53931751083799*^9, 3.539317637583722*^9, 3.539317710401309*^9, 
   3.539318367080615*^9, 3.539318857563046*^9, {3.5559442763044662`*^9, 
   3.555944298931198*^9}, 3.555944377919302*^9, 3.5559446052682056`*^9, 
   3.555944681456799*^9, 3.577373391273698*^9, {3.577373428852648*^9, 
   3.5773734513303223`*^9}, 3.5773742883918943`*^9, {3.578100113697337*^9, 
   3.5781001512116632`*^9}, 3.579349307063857*^9, 3.579350284229608*^9, 
   3.579350513044907*^9, 3.579350634908613*^9, 3.579358789778343*^9, {
   3.579359210417015*^9, 3.579359238091683*^9}, 3.5793595291764708`*^9, 
   3.579361067924972*^9, 3.579361288919224*^9, 3.579882686436294*^9, 
   3.579882985016074*^9, 3.5799557837408743`*^9, 3.6142269017242403`*^9, 
   3.614227200332449*^9, 3.617162111297064*^9, 3.617801329548204*^9, 
   3.61780165248982*^9, 3.617813330217751*^9, 3.677924497275835*^9, 
   3.677926813913814*^9}],

Cell[BoxData["\<\"The FeynRules palette can be opened using the command \
FRPalette[].\"\>"], "Print",
 CellChangeTimes->{
  3.53931751083799*^9, 3.539317637583722*^9, 3.539317710401309*^9, 
   3.539318367080615*^9, 3.539318857563046*^9, {3.5559442763044662`*^9, 
   3.555944298931198*^9}, 3.555944377919302*^9, 3.5559446052682056`*^9, 
   3.555944681456799*^9, 3.577373391273698*^9, {3.577373428852648*^9, 
   3.5773734513303223`*^9}, 3.5773742883918943`*^9, {3.578100113697337*^9, 
   3.5781001512116632`*^9}, 3.579349307063857*^9, 3.579350284229608*^9, 
   3.579350513044907*^9, 3.579350634908613*^9, 3.579358789778343*^9, {
   3.579359210417015*^9, 3.579359238091683*^9}, 3.5793595291764708`*^9, 
   3.579361067924972*^9, 3.579361288919224*^9, 3.579882686436294*^9, 
   3.579882985016074*^9, 3.5799557837408743`*^9, 3.6142269017242403`*^9, 
   3.614227200332449*^9, 3.617162111297064*^9, 3.617801329548204*^9, 
   3.61780165248982*^9, 3.617813330217751*^9, 3.677924497275835*^9, 
   3.677926813918252*^9}]
}, Open  ]]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Loading the model", "Section",
 CellChangeTimes->{{3.539317514565621*^9, 3.539317517827942*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"SetDirectory", "[", 
  "\"\<~/feynrules/branch/feynrules-current/Models/SM\>\"", "]"}]], "Input",
 CellChangeTimes->{{3.5393175793604097`*^9, 3.539317604969696*^9}, {
  3.5559442808384323`*^9, 3.5559442817000427`*^9}, {3.5773734202527857`*^9, 
  3.577373422619042*^9}, {3.579349331911512*^9, 3.579349394474436*^9}, {
  3.579361615574547*^9, 3.5793616313979177`*^9}, {3.579882727446802*^9, 
  3.5798827505288*^9}, {3.579883380265169*^9, 3.579883386157868*^9}, {
  3.579955791202548*^9, 3.579955816296172*^9}, {3.579958189936842*^9, 
  3.579958196374*^9}, {3.614226921305686*^9, 3.614226934448045*^9}, {
  3.614227154675522*^9, 3.614227186038155*^9}, {3.617801666219927*^9, 
  3.617801678049768*^9}, {3.677924508372898*^9, 3.6779245137278633`*^9}}],

Cell[BoxData["\<\"/Users/degrande/feynrules/branch/feynrules-current/Models/\
SM\"\>"], "Output",
 CellChangeTimes->{
  3.539317605382585*^9, 3.5393176390688963`*^9, 3.539317711322501*^9, 
   3.539318369381093*^9, 3.539318860145067*^9, {3.555944282499632*^9, 
   3.5559443001776247`*^9}, 3.555944378775016*^9, 3.5559446066309633`*^9, 
   3.555944683120018*^9, 3.577373454442539*^9, 3.5773742914690447`*^9, 
   3.578100155167124*^9, {3.579349388978545*^9, 3.579349395200293*^9}, 
   3.579350289156604*^9, 3.5793505173865347`*^9, 3.5793506389325314`*^9, 
   3.5793587910217247`*^9, 3.579359242533073*^9, 3.579359532748375*^9, 
   3.579361072123502*^9, 3.579361293397084*^9, 3.579882754989637*^9, 
   3.579882990732761*^9, 3.5799558183623962`*^9, 3.614226935661278*^9, {
   3.614227187078137*^9, 3.614227211205055*^9}, 3.6171621404829483`*^9, 
   3.617801334122128*^9, 3.617801681023378*^9, 3.6178133354359417`*^9, 
   3.6779245423234262`*^9, 3.6779268170201674`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"LoadModel", "[", 
  RowBox[{"\"\<SM.fr\>\"", ",", "\"\<~/Talks/MC4BSM16/Tutorial.fr\>\""}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.539317608005052*^9, 3.539317616212493*^9}, {
  3.555944290574006*^9, 3.555944292699313*^9}, {3.57737340515932*^9, 
  3.577373409019948*^9}, {3.579349400449449*^9, 3.579349402253819*^9}, {
  3.6142271566600933`*^9, 3.6142271794397287`*^9}, {3.6171621252630157`*^9, 
  3.617162134260457*^9}, {3.677924520505478*^9, 3.677924526246369*^9}}],

Cell[CellGroupData[{

Cell[BoxData["\<\"Merging model-files...\"\>"], "Print",
 CellChangeTimes->{{3.539317616899459*^9, 3.53931763972353*^9}, 
   3.539317711767226*^9, 3.539318370032514*^9, 3.53931886085401*^9, {
   3.555944284318446*^9, 3.555944300911276*^9}, 3.555944379117711*^9, 
   3.555944608282445*^9, 3.555944684535471*^9, 3.577373414696068*^9, {
   3.577373444854903*^9, 3.577373455221212*^9}, 3.577374292331169*^9, 
   3.578100156044696*^9, 3.579349403734391*^9, 3.5793502900702133`*^9, 
   3.579350518187063*^9, 3.579350639663134*^9, 3.579358791096305*^9, 
   3.579359243796385*^9, 3.579359533695651*^9, 3.579361073489555*^9, 
   3.579361294460967*^9, 3.579882757761602*^9, 3.579882991463509*^9, 
   3.579955819623869*^9, 3.614226941669352*^9, 3.614227213388626*^9, 
   3.617162142194365*^9, 3.617801335267908*^9, 3.6178016848792973`*^9, 
   3.617813337096767*^9, 3.677924544986567*^9, 3.677926818074012*^9}],

Cell[BoxData["\<\"This model implementation was created by\"\>"], "Print",
 CellChangeTimes->{{3.539317616899459*^9, 3.53931763972353*^9}, 
   3.539317711767226*^9, 3.539318370032514*^9, 3.53931886085401*^9, {
   3.555944284318446*^9, 3.555944300911276*^9}, 3.555944379117711*^9, 
   3.555944608282445*^9, 3.555944684535471*^9, 3.577373414696068*^9, {
   3.577373444854903*^9, 3.577373455221212*^9}, 3.577374292331169*^9, 
   3.578100156044696*^9, 3.579349403734391*^9, 3.5793502900702133`*^9, 
   3.579350518187063*^9, 3.579350639663134*^9, 3.579358791096305*^9, 
   3.579359243796385*^9, 3.579359533695651*^9, 3.579361073489555*^9, 
   3.579361294460967*^9, 3.579882757761602*^9, 3.579882991463509*^9, 
   3.579955819623869*^9, 3.614226941669352*^9, 3.614227213388626*^9, 
   3.617162142194365*^9, 3.617801335267908*^9, 3.6178016848792973`*^9, 
   3.617813337096767*^9, 3.677924544986567*^9, 3.677926818096286*^9}],

Cell[BoxData["\<\"C. Duhr\"\>"], "Print",
 CellChangeTimes->{{3.539317616899459*^9, 3.53931763972353*^9}, 
   3.539317711767226*^9, 3.539318370032514*^9, 3.53931886085401*^9, {
   3.555944284318446*^9, 3.555944300911276*^9}, 3.555944379117711*^9, 
   3.555944608282445*^9, 3.555944684535471*^9, 3.577373414696068*^9, {
   3.577373444854903*^9, 3.577373455221212*^9}, 3.577374292331169*^9, 
   3.578100156044696*^9, 3.579349403734391*^9, 3.5793502900702133`*^9, 
   3.579350518187063*^9, 3.579350639663134*^9, 3.579358791096305*^9, 
   3.579359243796385*^9, 3.579359533695651*^9, 3.579361073489555*^9, 
   3.579361294460967*^9, 3.579882757761602*^9, 3.579882991463509*^9, 
   3.579955819623869*^9, 3.614226941669352*^9, 3.614227213388626*^9, 
   3.617162142194365*^9, 3.617801335267908*^9, 3.6178016848792973`*^9, 
   3.617813337096767*^9, 3.677924544986567*^9, 3.677926818102749*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Model Version: \"\>", "\[InvisibleSpace]", "\<\"1.0\"\>"}],
  SequenceForm["Model Version: ", "1.0"],
  Editable->False]], "Print",
 CellChangeTimes->{{3.539317616899459*^9, 3.53931763972353*^9}, 
   3.539317711767226*^9, 3.539318370032514*^9, 3.53931886085401*^9, {
   3.555944284318446*^9, 3.555944300911276*^9}, 3.555944379117711*^9, 
   3.555944608282445*^9, 3.555944684535471*^9, 3.577373414696068*^9, {
   3.577373444854903*^9, 3.577373455221212*^9}, 3.577374292331169*^9, 
   3.578100156044696*^9, 3.579349403734391*^9, 3.5793502900702133`*^9, 
   3.579350518187063*^9, 3.579350639663134*^9, 3.579358791096305*^9, 
   3.579359243796385*^9, 3.579359533695651*^9, 3.579361073489555*^9, 
   3.579361294460967*^9, 3.579882757761602*^9, 3.579882991463509*^9, 
   3.579955819623869*^9, 3.614226941669352*^9, 3.614227213388626*^9, 
   3.617162142194365*^9, 3.617801335267908*^9, 3.6178016848792973`*^9, 
   3.617813337096767*^9, 3.677924544986567*^9, 3.6779268181097717`*^9}],

Cell[BoxData["\<\"For more information, type ModelInformation[].\"\>"], \
"Print",
 CellChangeTimes->{{3.539317616899459*^9, 3.53931763972353*^9}, 
   3.539317711767226*^9, 3.539318370032514*^9, 3.53931886085401*^9, {
   3.555944284318446*^9, 3.555944300911276*^9}, 3.555944379117711*^9, 
   3.555944608282445*^9, 3.555944684535471*^9, 3.577373414696068*^9, {
   3.577373444854903*^9, 3.577373455221212*^9}, 3.577374292331169*^9, 
   3.578100156044696*^9, 3.579349403734391*^9, 3.5793502900702133`*^9, 
   3.579350518187063*^9, 3.579350639663134*^9, 3.579358791096305*^9, 
   3.579359243796385*^9, 3.579359533695651*^9, 3.579361073489555*^9, 
   3.579361294460967*^9, 3.579882757761602*^9, 3.579882991463509*^9, 
   3.579955819623869*^9, 3.614226941669352*^9, 3.614227213388626*^9, 
   3.617162142194365*^9, 3.617801335267908*^9, 3.6178016848792973`*^9, 
   3.617813337096767*^9, 3.677924544986567*^9, 3.677926818116287*^9}],

Cell[BoxData["\<\"\"\>"], "Print",
 CellChangeTimes->{{3.539317616899459*^9, 3.53931763972353*^9}, 
   3.539317711767226*^9, 3.539318370032514*^9, 3.53931886085401*^9, {
   3.555944284318446*^9, 3.555944300911276*^9}, 3.555944379117711*^9, 
   3.555944608282445*^9, 3.555944684535471*^9, 3.577373414696068*^9, {
   3.577373444854903*^9, 3.577373455221212*^9}, 3.577374292331169*^9, 
   3.578100156044696*^9, 3.579349403734391*^9, 3.5793502900702133`*^9, 
   3.579350518187063*^9, 3.579350639663134*^9, 3.579358791096305*^9, 
   3.579359243796385*^9, 3.579359533695651*^9, 3.579361073489555*^9, 
   3.579361294460967*^9, 3.579882757761602*^9, 3.579882991463509*^9, 
   3.579955819623869*^9, 3.614226941669352*^9, 3.614227213388626*^9, 
   3.617162142194365*^9, 3.617801335267908*^9, 3.6178016848792973`*^9, 
   3.617813337096767*^9, 3.677924544986567*^9, 3.677926818123283*^9}],

Cell[BoxData["\<\"   - Loading particle classes.\"\>"], "Print",
 CellChangeTimes->{{3.539317616899459*^9, 3.53931763972353*^9}, 
   3.539317711767226*^9, 3.539318370032514*^9, 3.53931886085401*^9, {
   3.555944284318446*^9, 3.555944300911276*^9}, 3.555944379117711*^9, 
   3.555944608282445*^9, 3.555944684535471*^9, 3.577373414696068*^9, {
   3.577373444854903*^9, 3.577373455221212*^9}, 3.577374292331169*^9, 
   3.578100156044696*^9, 3.579349403734391*^9, 3.5793502900702133`*^9, 
   3.579350518187063*^9, 3.579350639663134*^9, 3.579358791096305*^9, 
   3.579359243796385*^9, 3.579359533695651*^9, 3.579361073489555*^9, 
   3.579361294460967*^9, 3.579882757761602*^9, 3.579882991463509*^9, 
   3.579955819623869*^9, 3.614226941669352*^9, 3.614227213388626*^9, 
   3.617162142194365*^9, 3.617801335267908*^9, 3.6178016848792973`*^9, 
   3.617813337096767*^9, 3.677924544986567*^9, 3.677926818127863*^9}],

Cell[BoxData["\<\"   - Loading gauge group classes.\"\>"], "Print",
 CellChangeTimes->{{3.539317616899459*^9, 3.53931763972353*^9}, 
   3.539317711767226*^9, 3.539318370032514*^9, 3.53931886085401*^9, {
   3.555944284318446*^9, 3.555944300911276*^9}, 3.555944379117711*^9, 
   3.555944608282445*^9, 3.555944684535471*^9, 3.577373414696068*^9, {
   3.577373444854903*^9, 3.577373455221212*^9}, 3.577374292331169*^9, 
   3.578100156044696*^9, 3.579349403734391*^9, 3.5793502900702133`*^9, 
   3.579350518187063*^9, 3.579350639663134*^9, 3.579358791096305*^9, 
   3.579359243796385*^9, 3.579359533695651*^9, 3.579361073489555*^9, 
   3.579361294460967*^9, 3.579882757761602*^9, 3.579882991463509*^9, 
   3.579955819623869*^9, 3.614226941669352*^9, 3.614227213388626*^9, 
   3.617162142194365*^9, 3.617801335267908*^9, 3.6178016848792973`*^9, 
   3.617813337096767*^9, 3.677924544986567*^9, 3.677926818259626*^9}],

Cell[BoxData["\<\"   - Loading parameter classes.\"\>"], "Print",
 CellChangeTimes->{{3.539317616899459*^9, 3.53931763972353*^9}, 
   3.539317711767226*^9, 3.539318370032514*^9, 3.53931886085401*^9, {
   3.555944284318446*^9, 3.555944300911276*^9}, 3.555944379117711*^9, 
   3.555944608282445*^9, 3.555944684535471*^9, 3.577373414696068*^9, {
   3.577373444854903*^9, 3.577373455221212*^9}, 3.577374292331169*^9, 
   3.578100156044696*^9, 3.579349403734391*^9, 3.5793502900702133`*^9, 
   3.579350518187063*^9, 3.579350639663134*^9, 3.579358791096305*^9, 
   3.579359243796385*^9, 3.579359533695651*^9, 3.579361073489555*^9, 
   3.579361294460967*^9, 3.579882757761602*^9, 3.579882991463509*^9, 
   3.579955819623869*^9, 3.614226941669352*^9, 3.614227213388626*^9, 
   3.617162142194365*^9, 3.617801335267908*^9, 3.6178016848792973`*^9, 
   3.617813337096767*^9, 3.677924544986567*^9, 3.677926818297143*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"\\nModel \"\>", "\[InvisibleSpace]", "\<\"Tutorial\"\>", 
   "\[InvisibleSpace]", "\<\" loaded.\"\>"}],
  SequenceForm["\nModel ", "Tutorial", " loaded."],
  Editable->False]], "Print",
 CellChangeTimes->{{3.539317616899459*^9, 3.53931763972353*^9}, 
   3.539317711767226*^9, 3.539318370032514*^9, 3.53931886085401*^9, {
   3.555944284318446*^9, 3.555944300911276*^9}, 3.555944379117711*^9, 
   3.555944608282445*^9, 3.555944684535471*^9, 3.577373414696068*^9, {
   3.577373444854903*^9, 3.577373455221212*^9}, 3.577374292331169*^9, 
   3.578100156044696*^9, 3.579349403734391*^9, 3.5793502900702133`*^9, 
   3.579350518187063*^9, 3.579350639663134*^9, 3.579358791096305*^9, 
   3.579359243796385*^9, 3.579359533695651*^9, 3.579361073489555*^9, 
   3.579361294460967*^9, 3.579882757761602*^9, 3.579882991463509*^9, 
   3.579955819623869*^9, 3.614226941669352*^9, 3.614227213388626*^9, 
   3.617162142194365*^9, 3.617801335267908*^9, 3.6178016848792973`*^9, 
   3.617813337096767*^9, 3.677924544986567*^9, 3.677926818361171*^9}]
}, Open  ]]
}, Open  ]],

Cell["\<\
\[LineSeparator]Restriction files, to make the quarks massless, and the CKM \
matrix diagonal:\
\>", "Text",
 CellChangeTimes->{{3.539317826891835*^9, 3.539317839321149*^9}, 
   3.539318861815209*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"LoadRestriction", "[", 
  RowBox[{"\"\<DiagonalCKM.rst\>\"", ",", "\"\<Massless.rst\>\""}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.539317840243554*^9, 3.53931785155518*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Loading restrictions from \"\>", 
   "\[InvisibleSpace]", "\<\"DiagonalCKM.rst\"\>", 
   "\[InvisibleSpace]", "\<\" : \"\>", "\[InvisibleSpace]", 
   DynamicBox[ToBoxes[PRIVATE`FR$restrictionCounter, StandardForm],
    ImageSizeCache->{22., {1., 12.}}], "\[InvisibleSpace]", "\<\" / \"\>", 
   "\[InvisibleSpace]", "3"}],
  SequenceForm["Loading restrictions from ", "DiagonalCKM.rst", " : ", 
   Dynamic[PRIVATE`FR$restrictionCounter], " / ", 3],
  Editable->False]], "Print",
 CellChangeTimes->{3.539317852149496*^9, 3.5393183725017023`*^9, 
  3.5393188635676622`*^9, 3.555944382559485*^9, 3.5559446108120213`*^9, 
  3.5559446867107964`*^9, 3.57737345755239*^9, 3.577374294631625*^9, 
  3.578100158857573*^9, 3.579349409804182*^9, 3.579350295931012*^9, 
  3.5793505209064093`*^9, 3.5793506425014772`*^9, 3.579358794923716*^9, 
  3.579359214019527*^9, 3.57935924743237*^9, 3.5793595386756287`*^9, 
  3.579361075960133*^9, 3.579361297098818*^9, 3.579882761262206*^9, 
  3.579882994106968*^9, 3.579955824463958*^9, 3.61422695186008*^9, 
  3.614227219473942*^9, 3.6171621480969677`*^9, 3.617801336926084*^9, 
  3.617801690218315*^9, 3.617813340805356*^9, 3.677924550893873*^9, 
  3.6779268208606853`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Loading restrictions from \"\>", 
   "\[InvisibleSpace]", "\<\"Massless.rst\"\>", 
   "\[InvisibleSpace]", "\<\" : \"\>", "\[InvisibleSpace]", 
   DynamicBox[ToBoxes[PRIVATE`FR$restrictionCounter, StandardForm],
    ImageSizeCache->{22., {1., 12.}}], "\[InvisibleSpace]", "\<\" / \"\>", 
   "\[InvisibleSpace]", "18"}],
  SequenceForm["Loading restrictions from ", "Massless.rst", " : ", 
   Dynamic[PRIVATE`FR$restrictionCounter], " / ", 18],
  Editable->False]], "Print",
 CellChangeTimes->{3.539317852149496*^9, 3.5393183725017023`*^9, 
  3.5393188635676622`*^9, 3.555944382559485*^9, 3.5559446108120213`*^9, 
  3.5559446867107964`*^9, 3.57737345755239*^9, 3.577374294631625*^9, 
  3.578100158857573*^9, 3.579349409804182*^9, 3.579350295931012*^9, 
  3.5793505209064093`*^9, 3.5793506425014772`*^9, 3.579358794923716*^9, 
  3.579359214019527*^9, 3.57935924743237*^9, 3.5793595386756287`*^9, 
  3.579361075960133*^9, 3.579361297098818*^9, 3.579882761262206*^9, 
  3.579882994106968*^9, 3.579955824463958*^9, 3.61422695186008*^9, 
  3.614227219473942*^9, 3.6171621480969677`*^9, 3.617801336926084*^9, 
  3.617801690218315*^9, 3.617813340805356*^9, 3.677924550893873*^9, 
  3.677926821013029*^9}],

Cell[BoxData["\<\"Restrictions loaded.\"\>"], "Print",
 CellChangeTimes->{3.539317852149496*^9, 3.5393183725017023`*^9, 
  3.5393188635676622`*^9, 3.555944382559485*^9, 3.5559446108120213`*^9, 
  3.5559446867107964`*^9, 3.57737345755239*^9, 3.577374294631625*^9, 
  3.578100158857573*^9, 3.579349409804182*^9, 3.579350295931012*^9, 
  3.5793505209064093`*^9, 3.5793506425014772`*^9, 3.579358794923716*^9, 
  3.579359214019527*^9, 3.57935924743237*^9, 3.5793595386756287`*^9, 
  3.579361075960133*^9, 3.579361297098818*^9, 3.579882761262206*^9, 
  3.579882994106968*^9, 3.579955824463958*^9, 3.61422695186008*^9, 
  3.614227219473942*^9, 3.6171621480969677`*^9, 3.617801336926084*^9, 
  3.617801690218315*^9, 3.617813340805356*^9, 3.677924550893873*^9, 
  3.677926821123783*^9}]
}, Open  ]]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["The Lagrangian", "Section",
 CellChangeTimes->{{3.539317863455777*^9, 3.5393178650954237`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"LScalarKin", "=", 
  RowBox[{
   RowBox[{
    RowBox[{"1", "/", "2"}], 
    RowBox[{"del", "[", 
     RowBox[{"pi1", ",", "mu"}], "]"}], 
    RowBox[{"del", "[", 
     RowBox[{"pi1", ",", "mu"}], "]"}]}], "+", 
   RowBox[{
    RowBox[{"1", "/", "2"}], 
    RowBox[{"del", "[", 
     RowBox[{"pi2", ",", "mu"}], "]"}], 
    RowBox[{"del", "[", 
     RowBox[{"pi2", ",", "mu"}], "]"}]}], "-", 
   RowBox[{
    RowBox[{
     RowBox[{"MM1", "^", "2"}], "/", "2"}], 
    RowBox[{"pi1", "^", "2"}]}], "-", 
   RowBox[{
    RowBox[{
     RowBox[{"MM2", "^", "2"}], "/", "2"}], 
    RowBox[{"pi2", "^", "2"}]}], "-", 
   RowBox[{
    RowBox[{"MM12", "^", "2"}], "pi1", " ", "pi2"}]}]}]], "Input",
 CellChangeTimes->{{3.5393178657972727`*^9, 3.539317921175536*^9}, {
  3.539318477867846*^9, 3.539318481606659*^9}, {3.5393189015055037`*^9, 
  3.539318908479196*^9}, {3.5793591059670963`*^9, 3.5793591107344847`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"-", 
    FractionBox["1", "2"]}], " ", 
   SuperscriptBox["MM1", "2"], " ", 
   SuperscriptBox["pi1", "2"]}], "-", 
  RowBox[{
   SuperscriptBox["MM12", "2"], " ", "pi1", " ", "pi2"}], "-", 
  FractionBox[
   RowBox[{
    SuperscriptBox["MM2", "2"], " ", 
    SuperscriptBox["pi2", "2"]}], "2"], "+", 
  RowBox[{
   FractionBox["1", "2"], " ", 
   SuperscriptBox[
    RowBox[{
     SubscriptBox["\<\"\[PartialD]\"\>", "mu"], "[", "pi1", "]"}], "2"]}], 
  "+", 
  RowBox[{
   FractionBox["1", "2"], " ", 
   SuperscriptBox[
    RowBox[{
     SubscriptBox["\<\"\[PartialD]\"\>", "mu"], "[", "pi2", "]"}], 
    "2"]}]}]], "Output",
 CellChangeTimes->{3.539317921603055*^9, 3.539318374039912*^9, 
  3.53931848187812*^9, 3.539318864769656*^9, 3.539318909041028*^9, 
  3.555944385931329*^9, 3.555944612372089*^9, 3.555944688216881*^9, 
  3.577373459175234*^9, 3.5773742965846777`*^9, 3.578100162546414*^9, 
  3.5793494145868607`*^9, 3.57935029915375*^9, 3.579350524477511*^9, 
  3.5793506450603523`*^9, 3.579358798361402*^9, 3.5793591128689413`*^9, 
  3.5793592513840733`*^9, 3.579359541815214*^9, 3.579361078567272*^9, 
  3.579361300753628*^9, 3.579882765103462*^9, 3.579882999682539*^9, 
  3.579955826749527*^9, 3.61422696297958*^9, 3.614227223957199*^9, 
  3.617162156305806*^9, 3.6178013375726833`*^9, 3.617801693858933*^9, 
  3.6779245576132803`*^9, 3.67792683141957*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"LFermionKin", "=", 
  RowBox[{
   RowBox[{"I", " ", 
    RowBox[{"uvbar", ".", 
     RowBox[{"Ga", "[", "mu", "]"}], ".", 
     RowBox[{"DC", "[", 
      RowBox[{"uv", ",", "mu"}], "]"}]}]}], "-", 
   RowBox[{"Muv", " ", 
    RowBox[{"uvbar", ".", "uv"}]}], "+", 
   RowBox[{"I", " ", 
    RowBox[{"evbar", ".", 
     RowBox[{"Ga", "[", "mu", "]"}], ".", 
     RowBox[{"DC", "[", 
      RowBox[{"ev", ",", "mu"}], "]"}]}]}], "-", 
   RowBox[{"Mev", " ", 
    RowBox[{"evbar", ".", "ev"}]}]}]}]], "Input",
 CellChangeTimes->{{3.53931793307228*^9, 3.539317979893841*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"-", "Mev"}], " ", 
   RowBox[{
    OverscriptBox["ev", "\<\"-\"\>"], ".", "ev"}]}], "-", 
  RowBox[{"Muv", " ", 
   RowBox[{
    OverscriptBox["uv", "\<\"-\"\>"], ".", "uv"}]}], "+", 
  RowBox[{"\[ImaginaryI]", " ", 
   RowBox[{
    OverscriptBox["ev", "\<\"-\"\>"], ".", 
    TemplateBox[{"\[Gamma]","mu"},
     "Superscript"], ".", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"\[ImaginaryI]", " ", "ev", " ", 
       SubscriptBox["g", "1"], " ", 
       SubscriptBox["B", "mu"]}], "+", 
      RowBox[{
       SubscriptBox["\<\"\[PartialD]\"\>", "mu"], "[", "ev", "]"}]}], 
     ")"}]}]}], "+", 
  RowBox[{"\[ImaginaryI]", " ", 
   RowBox[{
    OverscriptBox["uv", "\<\"-\"\>"], ".", 
    TemplateBox[{"\[Gamma]","mu"},
     "Superscript"], ".", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{
       RowBox[{"-", 
        FractionBox["2", "3"]}], " ", "\[ImaginaryI]", " ", 
       SubscriptBox["g", "1"], " ", "uv", " ", 
       SubscriptBox["B", "mu"]}], "+", 
      RowBox[{
       SubscriptBox["\<\"\[PartialD]\"\>", "mu"], "[", "uv", "]"}], "-", 
      RowBox[{"\[ImaginaryI]", " ", 
       SubscriptBox["g", "s"], " ", 
       RowBox[{
        SuperscriptBox["T", "a$915"], ".", "uv"}], " ", 
       SubscriptBox["G", 
        RowBox[{"mu", ",", "a$915"}]]}]}], ")"}]}]}]}]], "Output",
 CellChangeTimes->{3.539317980933099*^9, 3.539318374673298*^9, 
  3.5393188655258207`*^9, 3.55594438840161*^9, 3.555944613159431*^9, 
  3.555944688946253*^9, 3.577373459938418*^9, 3.577374297687181*^9, 
  3.578100163601769*^9, 3.579349419437807*^9, 3.579350300083214*^9, 
  3.5793505254085903`*^9, 3.579350646554434*^9, 3.5793588007912607`*^9, 
  3.579359255802463*^9, 3.579359542895461*^9, 3.579361079781403*^9, 
  3.5793613017508173`*^9, 3.579882767926231*^9, 3.579883001062243*^9, 
  3.579955830269313*^9, 3.614226964940504*^9, 3.614227225237988*^9, 
  3.617162158118114*^9, 3.6178013385711737`*^9, 3.617801694456737*^9, 
  3.6779245584340343`*^9, 3.677926832431203*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"LDecays", "=", 
  RowBox[{
   RowBox[{"lam1", " ", "pi1", " ", 
    RowBox[{
     RowBox[{"uvbar", "[", 
      RowBox[{"ss1", ",", "ii"}], "]"}], ".", 
     RowBox[{"uq", "[", 
      RowBox[{"ss2", ",", "1", ",", "ii"}], "]"}]}], 
    RowBox[{"ProjP", "[", 
     RowBox[{"ss1", ",", "ss2"}], "]"}]}], "+", 
   RowBox[{"lam2", " ", "pi2", " ", 
    RowBox[{
     RowBox[{"uvbar", "[", 
      RowBox[{"ss1", ",", "ii"}], "]"}], ".", 
     RowBox[{"uq", "[", 
      RowBox[{"ss2", ",", "1", ",", "ii"}], "]"}]}], 
    RowBox[{"ProjP", "[", 
     RowBox[{"ss1", ",", "ss2"}], "]"}]}], "+", 
   RowBox[{"lam1p", " ", "pi1", " ", 
    RowBox[{
     RowBox[{"evbar", "[", "ss1", "]"}], ".", 
     RowBox[{"l", "[", 
      RowBox[{"ss2", ",", "1"}], "]"}]}], 
    RowBox[{"ProjP", "[", 
     RowBox[{"ss1", ",", "ss2"}], "]"}]}], "+", 
   RowBox[{"lam2p", " ", "pi2", " ", 
    RowBox[{
     RowBox[{"evbar", "[", "ss1", "]"}], ".", 
     RowBox[{"l", "[", 
      RowBox[{"ss2", ",", "1"}], "]"}]}], 
    RowBox[{"ProjP", "[", 
     RowBox[{"ss1", ",", "ss2"}], "]"}]}]}]}], "\[IndentingNewLine]", 
 RowBox[{"LDecays", "=", 
  RowBox[{"LDecays", "+", 
   RowBox[{"HC", "[", "LDecays", "]"}]}]}]}], "Input",
 CellChangeTimes->{{3.5393179924355803`*^9, 3.5393180812160053`*^9}, {
  3.539318112141931*^9, 3.5393181200908337`*^9}, {3.555944392672977*^9, 
  3.555944544909051*^9}, {3.5798827743589373`*^9, 3.5798828027591953`*^9}, {
  3.579955834704596*^9, 3.579955838470811*^9}, {3.6171621642375603`*^9, 
  3.6171621676048594`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"lam1p", " ", "pi1", " ", 
   RowBox[{
    SubscriptBox[
     OverscriptBox["ev", "\<\"-\"\>"], "ss1"], ".", 
    SubscriptBox["l", 
     RowBox[{"ss2", ",", "1"}]]}], " ", 
   SubscriptBox[
    SubscriptBox["P", "\<\"+\"\>"], 
    RowBox[{"ss1", ",", "ss2"}]]}], "+", 
  RowBox[{"lam2p", " ", "pi2", " ", 
   RowBox[{
    SubscriptBox[
     OverscriptBox["ev", "\<\"-\"\>"], "ss1"], ".", 
    SubscriptBox["l", 
     RowBox[{"ss2", ",", "1"}]]}], " ", 
   SubscriptBox[
    SubscriptBox["P", "\<\"+\"\>"], 
    RowBox[{"ss1", ",", "ss2"}]]}], "+", 
  RowBox[{"lam1", " ", "pi1", " ", 
   RowBox[{
    SubscriptBox[
     OverscriptBox["uv", "\<\"-\"\>"], 
     RowBox[{"ss1", ",", "ii"}]], ".", 
    SubscriptBox["uq", 
     RowBox[{"ss2", ",", "1", ",", "ii"}]]}], " ", 
   SubscriptBox[
    SubscriptBox["P", "\<\"+\"\>"], 
    RowBox[{"ss1", ",", "ss2"}]]}], "+", 
  RowBox[{"lam2", " ", "pi2", " ", 
   RowBox[{
    SubscriptBox[
     OverscriptBox["uv", "\<\"-\"\>"], 
     RowBox[{"ss1", ",", "ii"}]], ".", 
    SubscriptBox["uq", 
     RowBox[{"ss2", ",", "1", ",", "ii"}]]}], " ", 
   SubscriptBox[
    SubscriptBox["P", "\<\"+\"\>"], 
    RowBox[{"ss1", ",", "ss2"}]]}]}]], "Output",
 CellChangeTimes->{3.539318081693873*^9, 3.53931812063649*^9, 
  3.53931837527557*^9, 3.539318866186824*^9, 3.555944615341813*^9, 
  3.555944690785212*^9, 3.577373460709169*^9, 3.5773742984859943`*^9, 
  3.578100164435671*^9, 3.579349420510026*^9, 3.579350301178184*^9, 
  3.5793505263540382`*^9, 3.579350647449521*^9, 3.579358801648719*^9, 
  3.579359257144423*^9, 3.579359544057259*^9, 3.579361081206567*^9, 
  3.579361302645125*^9, 3.579882813206955*^9, 3.579883002139476*^9, 
  3.579955869220771*^9, 3.614226966320546*^9, 3.614227226617613*^9, 
  3.6171621723672523`*^9, 3.617801339219819*^9, 3.617801695057027*^9, 
  3.677924559542886*^9, 3.677926833177044*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"lam1p", " ", "pi1", " ", 
   RowBox[{
    SubscriptBox[
     OverscriptBox["l", "\<\"-\"\>"], 
     RowBox[{"r$921", ",", "1"}]], ".", 
    SubscriptBox["ev", "r$922"]}], " ", 
   SubscriptBox[
    SubscriptBox["P", "\<\"-\"\>"], 
    RowBox[{"r$921", ",", "r$922"}]]}], "+", 
  RowBox[{"lam2p", " ", "pi2", " ", 
   RowBox[{
    SubscriptBox[
     OverscriptBox["l", "\<\"-\"\>"], 
     RowBox[{"r$923", ",", "1"}]], ".", 
    SubscriptBox["ev", "r$924"]}], " ", 
   SubscriptBox[
    SubscriptBox["P", "\<\"-\"\>"], 
    RowBox[{"r$923", ",", "r$924"}]]}], "+", 
  RowBox[{"lam1", " ", "pi1", " ", 
   RowBox[{
    SubscriptBox[
     OverscriptBox["uq", "\<\"-\"\>"], 
     RowBox[{"r$925", ",", "1", ",", "ii"}]], ".", 
    SubscriptBox["uv", 
     RowBox[{"r$926", ",", "ii"}]]}], " ", 
   SubscriptBox[
    SubscriptBox["P", "\<\"-\"\>"], 
    RowBox[{"r$925", ",", "r$926"}]]}], "+", 
  RowBox[{"lam2", " ", "pi2", " ", 
   RowBox[{
    SubscriptBox[
     OverscriptBox["uq", "\<\"-\"\>"], 
     RowBox[{"r$927", ",", "1", ",", "ii"}]], ".", 
    SubscriptBox["uv", 
     RowBox[{"r$928", ",", "ii"}]]}], " ", 
   SubscriptBox[
    SubscriptBox["P", "\<\"-\"\>"], 
    RowBox[{"r$927", ",", "r$928"}]]}], "+", 
  RowBox[{"lam1p", " ", "pi1", " ", 
   RowBox[{
    SubscriptBox[
     OverscriptBox["ev", "\<\"-\"\>"], "ss1"], ".", 
    SubscriptBox["l", 
     RowBox[{"ss2", ",", "1"}]]}], " ", 
   SubscriptBox[
    SubscriptBox["P", "\<\"+\"\>"], 
    RowBox[{"ss1", ",", "ss2"}]]}], "+", 
  RowBox[{"lam2p", " ", "pi2", " ", 
   RowBox[{
    SubscriptBox[
     OverscriptBox["ev", "\<\"-\"\>"], "ss1"], ".", 
    SubscriptBox["l", 
     RowBox[{"ss2", ",", "1"}]]}], " ", 
   SubscriptBox[
    SubscriptBox["P", "\<\"+\"\>"], 
    RowBox[{"ss1", ",", "ss2"}]]}], "+", 
  RowBox[{"lam1", " ", "pi1", " ", 
   RowBox[{
    SubscriptBox[
     OverscriptBox["uv", "\<\"-\"\>"], 
     RowBox[{"ss1", ",", "ii"}]], ".", 
    SubscriptBox["uq", 
     RowBox[{"ss2", ",", "1", ",", "ii"}]]}], " ", 
   SubscriptBox[
    SubscriptBox["P", "\<\"+\"\>"], 
    RowBox[{"ss1", ",", "ss2"}]]}], "+", 
  RowBox[{"lam2", " ", "pi2", " ", 
   RowBox[{
    SubscriptBox[
     OverscriptBox["uv", "\<\"-\"\>"], 
     RowBox[{"ss1", ",", "ii"}]], ".", 
    SubscriptBox["uq", 
     RowBox[{"ss2", ",", "1", ",", "ii"}]]}], " ", 
   SubscriptBox[
    SubscriptBox["P", "\<\"+\"\>"], 
    RowBox[{"ss1", ",", "ss2"}]]}]}]], "Output",
 CellChangeTimes->{3.539318081693873*^9, 3.53931812063649*^9, 
  3.53931837527557*^9, 3.539318866186824*^9, 3.555944615341813*^9, 
  3.555944690785212*^9, 3.577373460709169*^9, 3.5773742984859943`*^9, 
  3.578100164435671*^9, 3.579349420510026*^9, 3.579350301178184*^9, 
  3.5793505263540382`*^9, 3.579350647449521*^9, 3.579358801648719*^9, 
  3.579359257144423*^9, 3.579359544057259*^9, 3.579361081206567*^9, 
  3.579361302645125*^9, 3.579882813206955*^9, 3.579883002139476*^9, 
  3.579955869220771*^9, 3.614226966320546*^9, 3.614227226617613*^9, 
  3.6171621723672523`*^9, 3.617801339219819*^9, 3.617801695057027*^9, 
  3.677924559542886*^9, 3.6779268331817007`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"LNew", "=", 
  RowBox[{"LScalarKin", "+", "LFermionKin", "+", "LDecays"}]}]], "Input",
 CellChangeTimes->{{3.5393180928264513`*^9, 3.539318101951029*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"-", 
    FractionBox["1", "2"]}], " ", 
   SuperscriptBox["MM1", "2"], " ", 
   SuperscriptBox["pi1", "2"]}], "-", 
  RowBox[{
   SuperscriptBox["MM12", "2"], " ", "pi1", " ", "pi2"}], "-", 
  FractionBox[
   RowBox[{
    SuperscriptBox["MM2", "2"], " ", 
    SuperscriptBox["pi2", "2"]}], "2"], "+", 
  RowBox[{
   FractionBox["1", "2"], " ", 
   SuperscriptBox[
    RowBox[{
     SubscriptBox["\<\"\[PartialD]\"\>", "mu"], "[", "pi1", "]"}], "2"]}], 
  "+", 
  RowBox[{
   FractionBox["1", "2"], " ", 
   SuperscriptBox[
    RowBox[{
     SubscriptBox["\<\"\[PartialD]\"\>", "mu"], "[", "pi2", "]"}], "2"]}], 
  "-", 
  RowBox[{"Mev", " ", 
   RowBox[{
    OverscriptBox["ev", "\<\"-\"\>"], ".", "ev"}]}], "-", 
  RowBox[{"Muv", " ", 
   RowBox[{
    OverscriptBox["uv", "\<\"-\"\>"], ".", "uv"}]}], "+", 
  RowBox[{"\[ImaginaryI]", " ", 
   RowBox[{
    OverscriptBox["ev", "\<\"-\"\>"], ".", 
    TemplateBox[{"\[Gamma]","mu"},
     "Superscript"], ".", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"\[ImaginaryI]", " ", "ev", " ", 
       SubscriptBox["g", "1"], " ", 
       SubscriptBox["B", "mu"]}], "+", 
      RowBox[{
       SubscriptBox["\<\"\[PartialD]\"\>", "mu"], "[", "ev", "]"}]}], 
     ")"}]}]}], "+", 
  RowBox[{"\[ImaginaryI]", " ", 
   RowBox[{
    OverscriptBox["uv", "\<\"-\"\>"], ".", 
    TemplateBox[{"\[Gamma]","mu"},
     "Superscript"], ".", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{
       RowBox[{"-", 
        FractionBox["2", "3"]}], " ", "\[ImaginaryI]", " ", 
       SubscriptBox["g", "1"], " ", "uv", " ", 
       SubscriptBox["B", "mu"]}], "+", 
      RowBox[{
       SubscriptBox["\<\"\[PartialD]\"\>", "mu"], "[", "uv", "]"}], "-", 
      RowBox[{"\[ImaginaryI]", " ", 
       SubscriptBox["g", "s"], " ", 
       RowBox[{
        SuperscriptBox["T", "a$915"], ".", "uv"}], " ", 
       SubscriptBox["G", 
        RowBox[{"mu", ",", "a$915"}]]}]}], ")"}]}]}], "+", 
  RowBox[{"lam1p", " ", "pi1", " ", 
   RowBox[{
    SubscriptBox[
     OverscriptBox["l", "\<\"-\"\>"], 
     RowBox[{"r$921", ",", "1"}]], ".", 
    SubscriptBox["ev", "r$922"]}], " ", 
   SubscriptBox[
    SubscriptBox["P", "\<\"-\"\>"], 
    RowBox[{"r$921", ",", "r$922"}]]}], "+", 
  RowBox[{"lam2p", " ", "pi2", " ", 
   RowBox[{
    SubscriptBox[
     OverscriptBox["l", "\<\"-\"\>"], 
     RowBox[{"r$923", ",", "1"}]], ".", 
    SubscriptBox["ev", "r$924"]}], " ", 
   SubscriptBox[
    SubscriptBox["P", "\<\"-\"\>"], 
    RowBox[{"r$923", ",", "r$924"}]]}], "+", 
  RowBox[{"lam1", " ", "pi1", " ", 
   RowBox[{
    SubscriptBox[
     OverscriptBox["uq", "\<\"-\"\>"], 
     RowBox[{"r$925", ",", "1", ",", "ii"}]], ".", 
    SubscriptBox["uv", 
     RowBox[{"r$926", ",", "ii"}]]}], " ", 
   SubscriptBox[
    SubscriptBox["P", "\<\"-\"\>"], 
    RowBox[{"r$925", ",", "r$926"}]]}], "+", 
  RowBox[{"lam2", " ", "pi2", " ", 
   RowBox[{
    SubscriptBox[
     OverscriptBox["uq", "\<\"-\"\>"], 
     RowBox[{"r$927", ",", "1", ",", "ii"}]], ".", 
    SubscriptBox["uv", 
     RowBox[{"r$928", ",", "ii"}]]}], " ", 
   SubscriptBox[
    SubscriptBox["P", "\<\"-\"\>"], 
    RowBox[{"r$927", ",", "r$928"}]]}], "+", 
  RowBox[{"lam1p", " ", "pi1", " ", 
   RowBox[{
    SubscriptBox[
     OverscriptBox["ev", "\<\"-\"\>"], "ss1"], ".", 
    SubscriptBox["l", 
     RowBox[{"ss2", ",", "1"}]]}], " ", 
   SubscriptBox[
    SubscriptBox["P", "\<\"+\"\>"], 
    RowBox[{"ss1", ",", "ss2"}]]}], "+", 
  RowBox[{"lam2p", " ", "pi2", " ", 
   RowBox[{
    SubscriptBox[
     OverscriptBox["ev", "\<\"-\"\>"], "ss1"], ".", 
    SubscriptBox["l", 
     RowBox[{"ss2", ",", "1"}]]}], " ", 
   SubscriptBox[
    SubscriptBox["P", "\<\"+\"\>"], 
    RowBox[{"ss1", ",", "ss2"}]]}], "+", 
  RowBox[{"lam1", " ", "pi1", " ", 
   RowBox[{
    SubscriptBox[
     OverscriptBox["uv", "\<\"-\"\>"], 
     RowBox[{"ss1", ",", "ii"}]], ".", 
    SubscriptBox["uq", 
     RowBox[{"ss2", ",", "1", ",", "ii"}]]}], " ", 
   SubscriptBox[
    SubscriptBox["P", "\<\"+\"\>"], 
    RowBox[{"ss1", ",", "ss2"}]]}], "+", 
  RowBox[{"lam2", " ", "pi2", " ", 
   RowBox[{
    SubscriptBox[
     OverscriptBox["uv", "\<\"-\"\>"], 
     RowBox[{"ss1", ",", "ii"}]], ".", 
    SubscriptBox["uq", 
     RowBox[{"ss2", ",", "1", ",", "ii"}]]}], " ", 
   SubscriptBox[
    SubscriptBox["P", "\<\"+\"\>"], 
    RowBox[{"ss1", ",", "ss2"}]]}]}]], "Output",
 CellChangeTimes->{{3.5393181032144747`*^9, 3.539318124801681*^9}, 
   3.539318376704116*^9, 3.539318868615995*^9, 3.5393189114404373`*^9, 
   3.5559447357518663`*^9, 3.577373461914164*^9, 3.5773742995246153`*^9, 
   3.5781001665407763`*^9, 3.579349422580274*^9, 3.5793503056076612`*^9, 
   3.5793505293823853`*^9, 3.579350649975281*^9, 3.5793588066536427`*^9, 
   3.579359141817561*^9, 3.579359259902369*^9, 3.579359547500091*^9, 
   3.579361082767495*^9, 3.57936130745341*^9, 3.5798828207523727`*^9, 
   3.579883006053547*^9, 3.5799558758662443`*^9, 3.614226967884952*^9, 
   3.614227229262*^9, 3.617162174761636*^9, 3.617801340733911*^9, 
   3.6178016966707*^9, 3.677924561378821*^9, 3.677926834199871*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Checking of the model", "Section",
 CellChangeTimes->{{3.5393178553626547`*^9, 3.539317857239614*^9}, {
  3.539318951729939*^9, 3.539318953814794*^9}, {3.677924795889032*^9, 
  3.6779247975870523`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"CheckHermiticity", "[", "LNew", "]"}]], "Input",
 CellChangeTimes->{{3.677924800555821*^9, 3.677924812514064*^9}}],

Cell[CellGroupData[{

Cell[BoxData["\<\"Checking for hermiticity by calculating the Feynman rules \
contained in L-HC[L].\"\>"], "Print",
 CellChangeTimes->{3.677924816254113*^9}],

Cell[BoxData["\<\"If the lagrangian is hermitian, then the number of vertices \
should be zero.\"\>"], "Print",
 CellChangeTimes->{3.677924816256845*^9}],

Cell[BoxData[
 StyleBox["\<\"Starting Feynman rule calculation.\"\>",
  StripOnInput->False,
  LineColor->RGBColor[1, 0.5, 0],
  FrontFaceColor->RGBColor[1, 0.5, 0],
  BackFaceColor->RGBColor[1, 0.5, 0],
  GraphicsColor->RGBColor[1, 0.5, 0],
  FontWeight->Bold,
  FontColor->RGBColor[1, 0.5, 0]]], "Print",
 CellChangeTimes->{3.677924817454109*^9}],

Cell[BoxData["\<\"Expanding the Lagrangian...\"\>"], "Print",
 CellChangeTimes->{3.677924817458523*^9}],

Cell[BoxData["\<\"Collecting the different structures that enter the \
vertex.\"\>"], "Print",
 CellChangeTimes->{3.677924817522003*^9}],

Cell[BoxData["\<\"No vertices found.\"\>"], "Print",
 CellChangeTimes->{3.677924817526813*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"0", "\[InvisibleSpace]", "\<\" vertices obtained.\"\>"}],
  SequenceForm[0, " vertices obtained."],
  Editable->False]], "Print",
 CellChangeTimes->{3.677924817530117*^9}],

Cell[BoxData["\<\"The lagrangian is hermitian.\"\>"], "Print",
 CellChangeTimes->{3.67792481757064*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"{", "}"}]], "Output",
 CellChangeTimes->{3.677924817575622*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"CheckKineticTermNormalisation", "[", "LNew", "]"}]], "Input",
 CellChangeTimes->{{3.6779248203203163`*^9, 3.6779248371505013`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Neglecting all terms with more than \"\>", 
   "\[InvisibleSpace]", "\<\"2\"\>", 
   "\[InvisibleSpace]", "\<\" particles.\"\>"}],
  SequenceForm["Neglecting all terms with more than ", "2", " particles."],
  Editable->False]], "Print",
 CellChangeTimes->{3.677924838066307*^9}],

Cell[BoxData["\<\"All kinetic terms are diagonal.\"\>"], "Print",
 CellChangeTimes->{3.677924838078898*^9}],

Cell[BoxData["\<\"All kinetic terms are correctly normalized.\"\>"], "Print",
 CellChangeTimes->{3.677924838472795*^9}]
}, Open  ]],

Cell[BoxData["True"], "Output",
 CellChangeTimes->{3.677924838476533*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"CheckMassSpectrum", "[", "LNew", "]"}]], "Input",
 CellChangeTimes->{{3.579350411293748*^9, 3.579350415283794*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Neglecting all terms with more than \"\>", 
   "\[InvisibleSpace]", "\<\"2\"\>", 
   "\[InvisibleSpace]", "\<\" particles.\"\>"}],
  SequenceForm["Neglecting all terms with more than ", "2", " particles."],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.5793504171576357`*^9, 3.5793505336303453`*^9, 3.57935065439296*^9, 
   3.579358811124146*^9, {3.579359119161418*^9, 3.579359145870646*^9}, 
   3.5793592657591133`*^9, 3.5793595516192408`*^9, 3.579361091411455*^9, 
   3.579882864653352*^9, 3.579883010247345*^9, 3.579955882985557*^9, 
   3.614226974878377*^9, 3.614227237152254*^9, 3.677924609960373*^9}],

Cell[BoxData["\<\"All mass terms are diagonal.\"\>"], "Print",
 CellChangeTimes->{
  3.5793504171576357`*^9, 3.5793505336303453`*^9, 3.57935065439296*^9, 
   3.579358811124146*^9, {3.579359119161418*^9, 3.579359145870646*^9}, 
   3.5793592657591133`*^9, 3.5793595516192408`*^9, 3.579361091411455*^9, 
   3.579882864653352*^9, 3.579883010247345*^9, 3.579955882985557*^9, 
   3.614226974878377*^9, 3.614227237152254*^9, 3.677924609972211*^9}],

Cell[BoxData["\<\"Getting mass spectrum.\"\>"], "Print",
 CellChangeTimes->{
  3.5793504171576357`*^9, 3.5793505336303453`*^9, 3.57935065439296*^9, 
   3.579358811124146*^9, {3.579359119161418*^9, 3.579359145870646*^9}, 
   3.5793592657591133`*^9, 3.5793595516192408`*^9, 3.579361091411455*^9, 
   3.579882864653352*^9, 3.579883010247345*^9, 3.579955882985557*^9, 
   3.614226974878377*^9, 3.614227237152254*^9, 3.677924610314842*^9}],

Cell[BoxData["\<\"Checking for less then 0.1% agreement with model file \
values.\"\>"], "Print",
 CellChangeTimes->{
  3.5793504171576357`*^9, 3.5793505336303453`*^9, 3.57935065439296*^9, 
   3.579358811124146*^9, {3.579359119161418*^9, 3.579359145870646*^9}, 
   3.5793592657591133`*^9, 3.5793595516192408`*^9, 3.579361091411455*^9, 
   3.579882864653352*^9, 3.579883010247345*^9, 3.579955882985557*^9, 
   3.614226974878377*^9, 3.614227237152254*^9, 3.677924610319518*^9}]
}, Open  ]],

Cell[BoxData[
 TagBox[GridBox[{
    {"\<\"Particle\"\>", "\<\"Analytic value\"\>", "\<\"Numerical value\"\>", \
"\<\"Model-file value\"\>"},
    {"p1", 
     RowBox[{
      SqrtBox["2"], " ", 
      SqrtBox[
       RowBox[{
        RowBox[{
         FractionBox["1", "2"], " ", 
         SuperscriptBox["MM2", "2"], " ", 
         SuperscriptBox[
          RowBox[{"Cos", "[", "th", "]"}], "2"]}], "-", 
        RowBox[{
         SuperscriptBox["MM12", "2"], " ", 
         RowBox[{"Cos", "[", "th", "]"}], " ", 
         RowBox[{"Sin", "[", "th", "]"}]}], "+", 
        RowBox[{
         FractionBox["1", "2"], " ", 
         SuperscriptBox["MM1", "2"], " ", 
         SuperscriptBox[
          RowBox[{"Sin", "[", "th", "]"}], "2"]}]}]]}], "0.9999968746825869`",
      "0.9999968746825563`"},
    {"p2", 
     RowBox[{
      SqrtBox["2"], " ", 
      SqrtBox[
       RowBox[{
        RowBox[{
         FractionBox["1", "2"], " ", 
         SuperscriptBox["MM1", "2"], " ", 
         SuperscriptBox[
          RowBox[{"Cos", "[", "th", "]"}], "2"]}], "+", 
        RowBox[{
         SuperscriptBox["MM12", "2"], " ", 
         RowBox[{"Cos", "[", "th", "]"}], " ", 
         RowBox[{"Sin", "[", "th", "]"}]}], "+", 
        RowBox[{
         FractionBox["1", "2"], " ", 
         SuperscriptBox["MM2", "2"], " ", 
         SuperscriptBox[
          RowBox[{"Sin", "[", "th", "]"}], "2"]}]}]]}], "100.00000003125314`",
      "100.00000003125312`"},
    {"ev", "Mev", "50.`", "50.`"},
    {"uv", "Muv", "400.`", "400.`"}
   },
   GridBoxAlignment->{
    "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
     "RowsIndexed" -> {}},
   GridBoxSpacings->{"Columns" -> {
       Offset[0.27999999999999997`], {
        Offset[2.0999999999999996`]}, 
       Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
       Offset[0.2], {
        Offset[0.4]}, 
       Offset[0.2]}, "RowsIndexed" -> {}}],
  Function[BoxForm`e$, 
   TableForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.579359266696863*^9, 3.579359552626623*^9, 
  3.5793610923537407`*^9, 3.5798828652033567`*^9, 3.579883010753553*^9, 
  3.579955883507883*^9, 3.6142269754416637`*^9, 3.614227237429613*^9, 
  3.6779246103507147`*^9}]
}, Open  ]]
}, Closed]],

Cell[CellGroupData[{

Cell["The Feynman rules", "Section",
 CellChangeTimes->{{3.5393189582943993`*^9, 3.5393189603102827`*^9}, 
   3.577374302527355*^9, 3.614226983421431*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"vertices", "=", 
  RowBox[{"FeynmanRules", "[", "LNew", "]"}]}]], "Input",
 CellChangeTimes->{{3.5393189611068983`*^9, 3.539318966823897*^9}, {
  3.578100172014543*^9, 3.578100173567128*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 StyleBox["\<\"Starting Feynman rule calculation.\"\>",
  StripOnInput->False,
  LineColor->RGBColor[1, 0.5, 0],
  FrontFaceColor->RGBColor[1, 0.5, 0],
  BackFaceColor->RGBColor[1, 0.5, 0],
  GraphicsColor->RGBColor[1, 0.5, 0],
  FontWeight->Bold,
  FontColor->RGBColor[1, 0.5, 0]]], "Print",
 CellChangeTimes->{3.5393189671984*^9, 3.5559447465349827`*^9, 
  3.5773734655358477`*^9, 3.578100176569536*^9, 3.579359675837954*^9, 
  3.579361096122386*^9, 3.579882875210784*^9, 3.5798830142493563`*^9, 
  3.579955924820272*^9, 3.6142269897167387`*^9, 3.6142272437433357`*^9, 
  3.6171621856857367`*^9, 3.6779248843062687`*^9}],

Cell[BoxData["\<\"Expanding the Lagrangian...\"\>"], "Print",
 CellChangeTimes->{3.5393189671984*^9, 3.5559447465349827`*^9, 
  3.5773734655358477`*^9, 3.578100176569536*^9, 3.579359675837954*^9, 
  3.579361096122386*^9, 3.579882875210784*^9, 3.5798830142493563`*^9, 
  3.579955924820272*^9, 3.6142269897167387`*^9, 3.6142272437433357`*^9, 
  3.6171621856857367`*^9, 3.6779248843113537`*^9}],

Cell[BoxData["\<\"Collecting the different structures that enter the \
vertex.\"\>"], "Print",
 CellChangeTimes->{3.5393189671984*^9, 3.5559447465349827`*^9, 
  3.5773734655358477`*^9, 3.578100176569536*^9, 3.579359675837954*^9, 
  3.579361096122386*^9, 3.579882875210784*^9, 3.5798830142493563`*^9, 
  3.579955924820272*^9, 3.6142269897167387`*^9, 3.6142272437433357`*^9, 
  3.6171621856857367`*^9, 3.677924884743679*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
  "13", "\[InvisibleSpace]", "\<\" possible non-zero vertices have been found \
-> starting the computation: \"\>", "\[InvisibleSpace]", 
   DynamicBox[ToBoxes[FeynRules`FR$FeynmanRules, StandardForm],
    ImageSizeCache->{270., {4., 13.}}], "\[InvisibleSpace]", "\<\" / \"\>", 
   "\[InvisibleSpace]", "13", "\[InvisibleSpace]", "\<\".\"\>"}],
  SequenceForm[
  13, " possible non-zero vertices have been found -> starting the \
computation: ", 
   Dynamic[FeynRules`FR$FeynmanRules], " / ", 13, "."],
  Editable->False]], "Print",
 CellChangeTimes->{3.5393189671984*^9, 3.5559447465349827`*^9, 
  3.5773734655358477`*^9, 3.578100176569536*^9, 3.579359675837954*^9, 
  3.579361096122386*^9, 3.579882875210784*^9, 3.5798830142493563`*^9, 
  3.579955924820272*^9, 3.6142269897167387`*^9, 3.6142272437433357`*^9, 
  3.6171621856857367`*^9, 3.677924884766302*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"QN", "::", "NonConserv"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Warning: non quantum number conserving vertex encountered!\
\"\>"}]], "Message", "MSG",
 CellChangeTimes->{3.57935967692867*^9, 3.57936109709974*^9, 
  3.579882876172476*^9, 3.579883015113233*^9, 3.579955925731978*^9, 
  3.614226990346407*^9, 3.614227244355294*^9, 3.617162186403069*^9, 
  3.677924885310598*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Quantum number \"\>", "\[InvisibleSpace]", "LeptonNumber", 
   "\[InvisibleSpace]", "\<\" not conserved in vertex \"\>", 
   "\[InvisibleSpace]", 
   RowBox[{"{", 
    RowBox[{
     OverscriptBox["e", "\<\"-\"\>"], ",", "ev", ",", "p2"}], "}"}], 
   "\[InvisibleSpace]", "\<\".\"\>"}],
  SequenceForm[
  "Quantum number ", $CellContext`LeptonNumber, 
   " not conserved in vertex ", {$CellContext`ebar, $CellContext`ev, \
$CellContext`p2}, "."],
  Editable->False]], "Print",
 CellChangeTimes->{3.5393189671984*^9, 3.5559447465349827`*^9, 
  3.5773734655358477`*^9, 3.578100176569536*^9, 3.579359675837954*^9, 
  3.579361096122386*^9, 3.579882875210784*^9, 3.5798830142493563`*^9, 
  3.579955924820272*^9, 3.6142269897167387`*^9, 3.6142272437433357`*^9, 
  3.6171621856857367`*^9, 3.677924885316235*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"QN", "::", "NonConserv"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Warning: non quantum number conserving vertex encountered!\
\"\>"}]], "Message", "MSG",
 CellChangeTimes->{3.57935967692867*^9, 3.57936109709974*^9, 
  3.579882876172476*^9, 3.579883015113233*^9, 3.579955925731978*^9, 
  3.614226990346407*^9, 3.614227244355294*^9, 3.617162186403069*^9, 
  3.677924885432951*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Quantum number \"\>", "\[InvisibleSpace]", "Y", 
   "\[InvisibleSpace]", "\<\" not conserved in vertex \"\>", 
   "\[InvisibleSpace]", 
   RowBox[{"{", 
    RowBox[{
     OverscriptBox["e", "\<\"-\"\>"], ",", "ev", ",", "p2"}], "}"}], 
   "\[InvisibleSpace]", "\<\".\"\>"}],
  SequenceForm[
  "Quantum number ", FeynRules`Y, 
   " not conserved in vertex ", {$CellContext`ebar, $CellContext`ev, \
$CellContext`p2}, "."],
  Editable->False]], "Print",
 CellChangeTimes->{3.5393189671984*^9, 3.5559447465349827`*^9, 
  3.5773734655358477`*^9, 3.578100176569536*^9, 3.579359675837954*^9, 
  3.579361096122386*^9, 3.579882875210784*^9, 3.5798830142493563`*^9, 
  3.579955924820272*^9, 3.6142269897167387`*^9, 3.6142272437433357`*^9, 
  3.6171621856857367`*^9, 3.6779248854381866`*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"QN", "::", "NonConserv"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Warning: non quantum number conserving vertex encountered!\
\"\>"}]], "Message", "MSG",
 CellChangeTimes->{3.57935967692867*^9, 3.57936109709974*^9, 
  3.579882876172476*^9, 3.579883015113233*^9, 3.579955925731978*^9, 
  3.614226990346407*^9, 3.614227244355294*^9, 3.617162186403069*^9, 
  3.677924885560487*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"General", "::", "stop"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Further output of \[NoBreak]\\!\\(\\*StyleBox[\\(QN :: \
NonConserv\\), \\\"MessageName\\\"]\\)\[NoBreak] will be suppressed during \
this calculation. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \
ButtonStyle->\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/General/stop\\\", ButtonNote -> \
\\\"General::stop\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{3.57935967692867*^9, 3.57936109709974*^9, 
  3.579882876172476*^9, 3.579883015113233*^9, 3.579955925731978*^9, 
  3.614226990346407*^9, 3.614227244355294*^9, 3.617162186403069*^9, 
  3.67792488572552*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Quantum number \"\>", "\[InvisibleSpace]", "LeptonNumber", 
   "\[InvisibleSpace]", "\<\" not conserved in vertex \"\>", 
   "\[InvisibleSpace]", 
   RowBox[{"{", 
    RowBox[{
     OverscriptBox["e", "\<\"-\"\>"], ",", "ev", ",", "p1"}], "}"}], 
   "\[InvisibleSpace]", "\<\".\"\>"}],
  SequenceForm[
  "Quantum number ", $CellContext`LeptonNumber, 
   " not conserved in vertex ", {$CellContext`ebar, $CellContext`ev, \
$CellContext`p1}, "."],
  Editable->False]], "Print",
 CellChangeTimes->{3.5393189671984*^9, 3.5559447465349827`*^9, 
  3.5773734655358477`*^9, 3.578100176569536*^9, 3.579359675837954*^9, 
  3.579361096122386*^9, 3.579882875210784*^9, 3.5798830142493563`*^9, 
  3.579955924820272*^9, 3.6142269897167387`*^9, 3.6142272437433357`*^9, 
  3.6171621856857367`*^9, 3.677924885731688*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Quantum number \"\>", "\[InvisibleSpace]", "Y", 
   "\[InvisibleSpace]", "\<\" not conserved in vertex \"\>", 
   "\[InvisibleSpace]", 
   RowBox[{"{", 
    RowBox[{
     OverscriptBox["e", "\<\"-\"\>"], ",", "ev", ",", "p1"}], "}"}], 
   "\[InvisibleSpace]", "\<\".\"\>"}],
  SequenceForm[
  "Quantum number ", FeynRules`Y, 
   " not conserved in vertex ", {$CellContext`ebar, $CellContext`ev, \
$CellContext`p1}, "."],
  Editable->False]], "Print",
 CellChangeTimes->{3.5393189671984*^9, 3.5559447465349827`*^9, 
  3.5773734655358477`*^9, 3.578100176569536*^9, 3.579359675837954*^9, 
  3.579361096122386*^9, 3.579882875210784*^9, 3.5798830142493563`*^9, 
  3.579955924820272*^9, 3.6142269897167387`*^9, 3.6142272437433357`*^9, 
  3.6171621856857367`*^9, 3.6779248857344627`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Quantum number \"\>", "\[InvisibleSpace]", "Y", 
   "\[InvisibleSpace]", "\<\" not conserved in vertex \"\>", 
   "\[InvisibleSpace]", 
   RowBox[{"{", 
    RowBox[{
     OverscriptBox["u", "\<\"-\"\>"], ",", "uv", ",", "p2"}], "}"}], 
   "\[InvisibleSpace]", "\<\".\"\>"}],
  SequenceForm[
  "Quantum number ", FeynRules`Y, 
   " not conserved in vertex ", {$CellContext`ubar, $CellContext`uv, \
$CellContext`p2}, "."],
  Editable->False]], "Print",
 CellChangeTimes->{3.5393189671984*^9, 3.5559447465349827`*^9, 
  3.5773734655358477`*^9, 3.578100176569536*^9, 3.579359675837954*^9, 
  3.579361096122386*^9, 3.579882875210784*^9, 3.5798830142493563`*^9, 
  3.579955924820272*^9, 3.6142269897167387`*^9, 3.6142272437433357`*^9, 
  3.6171621856857367`*^9, 3.6779248857391644`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Quantum number \"\>", "\[InvisibleSpace]", "Y", 
   "\[InvisibleSpace]", "\<\" not conserved in vertex \"\>", 
   "\[InvisibleSpace]", 
   RowBox[{"{", 
    RowBox[{
     OverscriptBox["u", "\<\"-\"\>"], ",", "uv", ",", "p1"}], "}"}], 
   "\[InvisibleSpace]", "\<\".\"\>"}],
  SequenceForm[
  "Quantum number ", FeynRules`Y, 
   " not conserved in vertex ", {$CellContext`ubar, $CellContext`uv, \
$CellContext`p1}, "."],
  Editable->False]], "Print",
 CellChangeTimes->{3.5393189671984*^9, 3.5559447465349827`*^9, 
  3.5773734655358477`*^9, 3.578100176569536*^9, 3.579359675837954*^9, 
  3.579361096122386*^9, 3.579882875210784*^9, 3.5798830142493563`*^9, 
  3.579955924820272*^9, 3.6142269897167387`*^9, 3.6142272437433357`*^9, 
  3.6171621856857367`*^9, 3.677924885743248*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Quantum number \"\>", "\[InvisibleSpace]", "LeptonNumber", 
   "\[InvisibleSpace]", "\<\" not conserved in vertex \"\>", 
   "\[InvisibleSpace]", 
   RowBox[{"{", 
    RowBox[{
     OverscriptBox["ev", "\<\"-\"\>"], ",", "e", ",", "p1"}], "}"}], 
   "\[InvisibleSpace]", "\<\".\"\>"}],
  SequenceForm[
  "Quantum number ", $CellContext`LeptonNumber, 
   " not conserved in vertex ", {$CellContext`evbar, 
    FeynRules`e, $CellContext`p1}, "."],
  Editable->False]], "Print",
 CellChangeTimes->{3.5393189671984*^9, 3.5559447465349827`*^9, 
  3.5773734655358477`*^9, 3.578100176569536*^9, 3.579359675837954*^9, 
  3.579361096122386*^9, 3.579882875210784*^9, 3.5798830142493563`*^9, 
  3.579955924820272*^9, 3.6142269897167387`*^9, 3.6142272437433357`*^9, 
  3.6171621856857367`*^9, 3.677924885746491*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Quantum number \"\>", "\[InvisibleSpace]", "Y", 
   "\[InvisibleSpace]", "\<\" not conserved in vertex \"\>", 
   "\[InvisibleSpace]", 
   RowBox[{"{", 
    RowBox[{
     OverscriptBox["ev", "\<\"-\"\>"], ",", "e", ",", "p1"}], "}"}], 
   "\[InvisibleSpace]", "\<\".\"\>"}],
  SequenceForm[
  "Quantum number ", FeynRules`Y, 
   " not conserved in vertex ", {$CellContext`evbar, 
    FeynRules`e, $CellContext`p1}, "."],
  Editable->False]], "Print",
 CellChangeTimes->{3.5393189671984*^9, 3.5559447465349827`*^9, 
  3.5773734655358477`*^9, 3.578100176569536*^9, 3.579359675837954*^9, 
  3.579361096122386*^9, 3.579882875210784*^9, 3.5798830142493563`*^9, 
  3.579955924820272*^9, 3.6142269897167387`*^9, 3.6142272437433357`*^9, 
  3.6171621856857367`*^9, 3.6779248857497263`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Quantum number \"\>", "\[InvisibleSpace]", "LeptonNumber", 
   "\[InvisibleSpace]", "\<\" not conserved in vertex \"\>", 
   "\[InvisibleSpace]", 
   RowBox[{"{", 
    RowBox[{
     OverscriptBox["ev", "\<\"-\"\>"], ",", "e", ",", "p2"}], "}"}], 
   "\[InvisibleSpace]", "\<\".\"\>"}],
  SequenceForm[
  "Quantum number ", $CellContext`LeptonNumber, 
   " not conserved in vertex ", {$CellContext`evbar, 
    FeynRules`e, $CellContext`p2}, "."],
  Editable->False]], "Print",
 CellChangeTimes->{3.5393189671984*^9, 3.5559447465349827`*^9, 
  3.5773734655358477`*^9, 3.578100176569536*^9, 3.579359675837954*^9, 
  3.579361096122386*^9, 3.579882875210784*^9, 3.5798830142493563`*^9, 
  3.579955924820272*^9, 3.6142269897167387`*^9, 3.6142272437433357`*^9, 
  3.6171621856857367`*^9, 3.677924885754506*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Quantum number \"\>", "\[InvisibleSpace]", "Y", 
   "\[InvisibleSpace]", "\<\" not conserved in vertex \"\>", 
   "\[InvisibleSpace]", 
   RowBox[{"{", 
    RowBox[{
     OverscriptBox["ev", "\<\"-\"\>"], ",", "e", ",", "p2"}], "}"}], 
   "\[InvisibleSpace]", "\<\".\"\>"}],
  SequenceForm[
  "Quantum number ", FeynRules`Y, 
   " not conserved in vertex ", {$CellContext`evbar, 
    FeynRules`e, $CellContext`p2}, "."],
  Editable->False]], "Print",
 CellChangeTimes->{3.5393189671984*^9, 3.5559447465349827`*^9, 
  3.5773734655358477`*^9, 3.578100176569536*^9, 3.579359675837954*^9, 
  3.579361096122386*^9, 3.579882875210784*^9, 3.5798830142493563`*^9, 
  3.579955924820272*^9, 3.6142269897167387`*^9, 3.6142272437433357`*^9, 
  3.6171621856857367`*^9, 3.6779248857578*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Quantum number \"\>", "\[InvisibleSpace]", "Y", 
   "\[InvisibleSpace]", "\<\" not conserved in vertex \"\>", 
   "\[InvisibleSpace]", 
   RowBox[{"{", 
    RowBox[{
     OverscriptBox["uv", "\<\"-\"\>"], ",", "u", ",", "p1"}], "}"}], 
   "\[InvisibleSpace]", "\<\".\"\>"}],
  SequenceForm[
  "Quantum number ", FeynRules`Y, 
   " not conserved in vertex ", {$CellContext`uvbar, 
    FeynRules`u, $CellContext`p1}, "."],
  Editable->False]], "Print",
 CellChangeTimes->{3.5393189671984*^9, 3.5559447465349827`*^9, 
  3.5773734655358477`*^9, 3.578100176569536*^9, 3.579359675837954*^9, 
  3.579361096122386*^9, 3.579882875210784*^9, 3.5798830142493563`*^9, 
  3.579955924820272*^9, 3.6142269897167387`*^9, 3.6142272437433357`*^9, 
  3.6171621856857367`*^9, 3.677924885761067*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Quantum number \"\>", "\[InvisibleSpace]", "Y", 
   "\[InvisibleSpace]", "\<\" not conserved in vertex \"\>", 
   "\[InvisibleSpace]", 
   RowBox[{"{", 
    RowBox[{
     OverscriptBox["uv", "\<\"-\"\>"], ",", "u", ",", "p2"}], "}"}], 
   "\[InvisibleSpace]", "\<\".\"\>"}],
  SequenceForm[
  "Quantum number ", FeynRules`Y, 
   " not conserved in vertex ", {$CellContext`uvbar, 
    FeynRules`u, $CellContext`p2}, "."],
  Editable->False]], "Print",
 CellChangeTimes->{3.5393189671984*^9, 3.5559447465349827`*^9, 
  3.5773734655358477`*^9, 3.578100176569536*^9, 3.579359675837954*^9, 
  3.579361096122386*^9, 3.579882875210784*^9, 3.5798830142493563`*^9, 
  3.579955924820272*^9, 3.6142269897167387`*^9, 3.6142272437433357`*^9, 
  3.6171621856857367`*^9, 3.6779248857641993`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"13", "\[InvisibleSpace]", "\<\" vertices obtained.\"\>"}],
  SequenceForm[13, " vertices obtained."],
  Editable->False]], "Print",
 CellChangeTimes->{3.5393189671984*^9, 3.5559447465349827`*^9, 
  3.5773734655358477`*^9, 3.578100176569536*^9, 3.579359675837954*^9, 
  3.579361096122386*^9, 3.579882875210784*^9, 3.5798830142493563`*^9, 
  3.579955924820272*^9, 3.6142269897167387`*^9, 3.6142272437433357`*^9, 
  3.6171621856857367`*^9, 3.677924885767263*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         OverscriptBox["ev", "\<\"-\"\>"], ",", "1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"ev", ",", "2"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"A", ",", "3"}], "}"}]}], "}"}], ",", 
     RowBox[{
      RowBox[{"-", "\[ImaginaryI]"}], " ", "e", " ", 
      TemplateBox[{SubscriptBox["\[Gamma]", 
         RowBox[{
           SubscriptBox["\"s\"", "1"], ",", 
           SubscriptBox["\"s\"", "2"]}]],SubscriptBox["\"\[Mu]\"", "3"]},
       "Superscript"]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         OverscriptBox["uv", "\<\"-\"\>"], ",", "1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"uv", ",", "2"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"A", ",", "3"}], "}"}]}], "}"}], ",", 
     RowBox[{
      FractionBox["2", "3"], " ", "\[ImaginaryI]", " ", "e", " ", 
      TemplateBox[{SubscriptBox["\[Gamma]", 
         RowBox[{
           SubscriptBox["\"s\"", "1"], ",", 
           SubscriptBox["\"s\"", "2"]}]],SubscriptBox["\"\[Mu]\"", "3"]},
       "Superscript"], " ", 
      SubscriptBox["\[Delta]", 
       RowBox[{
        SubscriptBox["\<\"m\"\>", "1"], ",", 
        SubscriptBox["\<\"m\"\>", "2"]}]]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         OverscriptBox["e", "\<\"-\"\>"], ",", "1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"ev", ",", "2"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"p2", ",", "3"}], "}"}]}], "}"}], ",", 
     RowBox[{
      RowBox[{"\[ImaginaryI]", " ", "lam1p", " ", 
       RowBox[{"Cos", "[", "th", "]"}], " ", 
       SubscriptBox[
        SubscriptBox["P", "\<\"-\"\>"], 
        RowBox[{
         SubscriptBox["\<\"s\"\>", "1"], ",", 
         SubscriptBox["\<\"s\"\>", "2"]}]]}], "+", 
      RowBox[{"\[ImaginaryI]", " ", "lam2p", " ", 
       SubscriptBox[
        SubscriptBox["P", "\<\"-\"\>"], 
        RowBox[{
         SubscriptBox["\<\"s\"\>", "1"], ",", 
         SubscriptBox["\<\"s\"\>", "2"]}]], " ", 
       RowBox[{"Sin", "[", "th", "]"}]}]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         OverscriptBox["e", "\<\"-\"\>"], ",", "1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"ev", ",", "2"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"p1", ",", "3"}], "}"}]}], "}"}], ",", 
     RowBox[{
      RowBox[{"\[ImaginaryI]", " ", "lam2p", " ", 
       RowBox[{"Cos", "[", "th", "]"}], " ", 
       SubscriptBox[
        SubscriptBox["P", "\<\"-\"\>"], 
        RowBox[{
         SubscriptBox["\<\"s\"\>", "1"], ",", 
         SubscriptBox["\<\"s\"\>", "2"]}]]}], "-", 
      RowBox[{"\[ImaginaryI]", " ", "lam1p", " ", 
       SubscriptBox[
        SubscriptBox["P", "\<\"-\"\>"], 
        RowBox[{
         SubscriptBox["\<\"s\"\>", "1"], ",", 
         SubscriptBox["\<\"s\"\>", "2"]}]], " ", 
       RowBox[{"Sin", "[", "th", "]"}]}]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         OverscriptBox["u", "\<\"-\"\>"], ",", "1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"uv", ",", "2"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"p2", ",", "3"}], "}"}]}], "}"}], ",", 
     RowBox[{
      RowBox[{"\[ImaginaryI]", " ", "lam1", " ", 
       RowBox[{"Cos", "[", "th", "]"}], " ", 
       SubscriptBox["\[Delta]", 
        RowBox[{
         SubscriptBox["\<\"m\"\>", "1"], ",", 
         SubscriptBox["\<\"m\"\>", "2"]}]], " ", 
       SubscriptBox[
        SubscriptBox["P", "\<\"-\"\>"], 
        RowBox[{
         SubscriptBox["\<\"s\"\>", "1"], ",", 
         SubscriptBox["\<\"s\"\>", "2"]}]]}], "+", 
      RowBox[{"\[ImaginaryI]", " ", "lam2", " ", 
       SubscriptBox["\[Delta]", 
        RowBox[{
         SubscriptBox["\<\"m\"\>", "1"], ",", 
         SubscriptBox["\<\"m\"\>", "2"]}]], " ", 
       SubscriptBox[
        SubscriptBox["P", "\<\"-\"\>"], 
        RowBox[{
         SubscriptBox["\<\"s\"\>", "1"], ",", 
         SubscriptBox["\<\"s\"\>", "2"]}]], " ", 
       RowBox[{"Sin", "[", "th", "]"}]}]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         OverscriptBox["u", "\<\"-\"\>"], ",", "1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"uv", ",", "2"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"p1", ",", "3"}], "}"}]}], "}"}], ",", 
     RowBox[{
      RowBox[{"\[ImaginaryI]", " ", "lam2", " ", 
       RowBox[{"Cos", "[", "th", "]"}], " ", 
       SubscriptBox["\[Delta]", 
        RowBox[{
         SubscriptBox["\<\"m\"\>", "1"], ",", 
         SubscriptBox["\<\"m\"\>", "2"]}]], " ", 
       SubscriptBox[
        SubscriptBox["P", "\<\"-\"\>"], 
        RowBox[{
         SubscriptBox["\<\"s\"\>", "1"], ",", 
         SubscriptBox["\<\"s\"\>", "2"]}]]}], "-", 
      RowBox[{"\[ImaginaryI]", " ", "lam1", " ", 
       SubscriptBox["\[Delta]", 
        RowBox[{
         SubscriptBox["\<\"m\"\>", "1"], ",", 
         SubscriptBox["\<\"m\"\>", "2"]}]], " ", 
       SubscriptBox[
        SubscriptBox["P", "\<\"-\"\>"], 
        RowBox[{
         SubscriptBox["\<\"s\"\>", "1"], ",", 
         SubscriptBox["\<\"s\"\>", "2"]}]], " ", 
       RowBox[{"Sin", "[", "th", "]"}]}]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         OverscriptBox["ev", "\<\"-\"\>"], ",", "1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"e", ",", "2"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"p1", ",", "3"}], "}"}]}], "}"}], ",", 
     RowBox[{
      RowBox[{"\[ImaginaryI]", " ", "lam2p", " ", 
       RowBox[{"Cos", "[", "th", "]"}], " ", 
       SubscriptBox[
        SubscriptBox["P", "\<\"+\"\>"], 
        RowBox[{
         SubscriptBox["\<\"s\"\>", "1"], ",", 
         SubscriptBox["\<\"s\"\>", "2"]}]]}], "-", 
      RowBox[{"\[ImaginaryI]", " ", "lam1p", " ", 
       SubscriptBox[
        SubscriptBox["P", "\<\"+\"\>"], 
        RowBox[{
         SubscriptBox["\<\"s\"\>", "1"], ",", 
         SubscriptBox["\<\"s\"\>", "2"]}]], " ", 
       RowBox[{"Sin", "[", "th", "]"}]}]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         OverscriptBox["ev", "\<\"-\"\>"], ",", "1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"e", ",", "2"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"p2", ",", "3"}], "}"}]}], "}"}], ",", 
     RowBox[{
      RowBox[{"\[ImaginaryI]", " ", "lam1p", " ", 
       RowBox[{"Cos", "[", "th", "]"}], " ", 
       SubscriptBox[
        SubscriptBox["P", "\<\"+\"\>"], 
        RowBox[{
         SubscriptBox["\<\"s\"\>", "1"], ",", 
         SubscriptBox["\<\"s\"\>", "2"]}]]}], "+", 
      RowBox[{"\[ImaginaryI]", " ", "lam2p", " ", 
       SubscriptBox[
        SubscriptBox["P", "\<\"+\"\>"], 
        RowBox[{
         SubscriptBox["\<\"s\"\>", "1"], ",", 
         SubscriptBox["\<\"s\"\>", "2"]}]], " ", 
       RowBox[{"Sin", "[", "th", "]"}]}]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         OverscriptBox["uv", "\<\"-\"\>"], ",", "1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"u", ",", "2"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"p1", ",", "3"}], "}"}]}], "}"}], ",", 
     RowBox[{
      RowBox[{"\[ImaginaryI]", " ", "lam2", " ", 
       RowBox[{"Cos", "[", "th", "]"}], " ", 
       SubscriptBox["\[Delta]", 
        RowBox[{
         SubscriptBox["\<\"m\"\>", "1"], ",", 
         SubscriptBox["\<\"m\"\>", "2"]}]], " ", 
       SubscriptBox[
        SubscriptBox["P", "\<\"+\"\>"], 
        RowBox[{
         SubscriptBox["\<\"s\"\>", "1"], ",", 
         SubscriptBox["\<\"s\"\>", "2"]}]]}], "-", 
      RowBox[{"\[ImaginaryI]", " ", "lam1", " ", 
       SubscriptBox["\[Delta]", 
        RowBox[{
         SubscriptBox["\<\"m\"\>", "1"], ",", 
         SubscriptBox["\<\"m\"\>", "2"]}]], " ", 
       SubscriptBox[
        SubscriptBox["P", "\<\"+\"\>"], 
        RowBox[{
         SubscriptBox["\<\"s\"\>", "1"], ",", 
         SubscriptBox["\<\"s\"\>", "2"]}]], " ", 
       RowBox[{"Sin", "[", "th", "]"}]}]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         OverscriptBox["uv", "\<\"-\"\>"], ",", "1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"u", ",", "2"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"p2", ",", "3"}], "}"}]}], "}"}], ",", 
     RowBox[{
      RowBox[{"\[ImaginaryI]", " ", "lam1", " ", 
       RowBox[{"Cos", "[", "th", "]"}], " ", 
       SubscriptBox["\[Delta]", 
        RowBox[{
         SubscriptBox["\<\"m\"\>", "1"], ",", 
         SubscriptBox["\<\"m\"\>", "2"]}]], " ", 
       SubscriptBox[
        SubscriptBox["P", "\<\"+\"\>"], 
        RowBox[{
         SubscriptBox["\<\"s\"\>", "1"], ",", 
         SubscriptBox["\<\"s\"\>", "2"]}]]}], "+", 
      RowBox[{"\[ImaginaryI]", " ", "lam2", " ", 
       SubscriptBox["\[Delta]", 
        RowBox[{
         SubscriptBox["\<\"m\"\>", "1"], ",", 
         SubscriptBox["\<\"m\"\>", "2"]}]], " ", 
       SubscriptBox[
        SubscriptBox["P", "\<\"+\"\>"], 
        RowBox[{
         SubscriptBox["\<\"s\"\>", "1"], ",", 
         SubscriptBox["\<\"s\"\>", "2"]}]], " ", 
       RowBox[{"Sin", "[", "th", "]"}]}]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         OverscriptBox["uv", "\<\"-\"\>"], ",", "1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"uv", ",", "2"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"G", ",", "3"}], "}"}]}], "}"}], ",", 
     RowBox[{"\[ImaginaryI]", " ", 
      SubscriptBox["g", "s"], " ", 
      TemplateBox[{SubscriptBox["\[Gamma]", 
         RowBox[{
           SubscriptBox["\"s\"", "1"], ",", 
           SubscriptBox["\"s\"", "2"]}]],SubscriptBox["\"\[Mu]\"", "3"]},
       "Superscript"], " ", 
      SubsuperscriptBox["T", 
       RowBox[{
        SubscriptBox["\<\"m\"\>", "1"], ",", 
        SubscriptBox["\<\"m\"\>", "2"]}], 
       SubscriptBox["\<\"a\"\>", "3"]]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         OverscriptBox["ev", "\<\"-\"\>"], ",", "1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"ev", ",", "2"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"Z", ",", "3"}], "}"}]}], "}"}], ",", 
     FractionBox[
      RowBox[{"\[ImaginaryI]", " ", "e", " ", 
       SubscriptBox["s", "w"], " ", 
       TemplateBox[{SubscriptBox["\[Gamma]", 
          RowBox[{
            SubscriptBox["\"s\"", "1"], ",", 
            SubscriptBox["\"s\"", "2"]}]],SubscriptBox["\"\[Mu]\"", "3"]},
        "Superscript"]}], 
      SubscriptBox["c", "w"]]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         OverscriptBox["uv", "\<\"-\"\>"], ",", "1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"uv", ",", "2"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"Z", ",", "3"}], "}"}]}], "}"}], ",", 
     RowBox[{"-", 
      FractionBox[
       RowBox[{"2", " ", "\[ImaginaryI]", " ", "e", " ", 
        SubscriptBox["s", "w"], " ", 
        TemplateBox[{SubscriptBox["\[Gamma]", 
           RowBox[{
             SubscriptBox["\"s\"", "1"], ",", 
             SubscriptBox["\"s\"", "2"]}]],SubscriptBox["\"\[Mu]\"", "3"]},
         "Superscript"], " ", 
        SubscriptBox["\[Delta]", 
         RowBox[{
          SubscriptBox["\<\"m\"\>", "1"], ",", 
          SubscriptBox["\<\"m\"\>", "2"]}]]}], 
       RowBox[{"3", " ", 
        SubscriptBox["c", "w"]}]]}]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.539318968485745*^9, 3.555944748124181*^9, 
  3.577373467002137*^9, 3.578100178320201*^9, 3.579359677559607*^9, 
  3.579361097674306*^9, 3.5798828767189703`*^9, 3.579883015711761*^9, 
  3.579955926297202*^9, 3.614226990639884*^9, 3.614227244578875*^9, 
  3.617162186681904*^9, 3.6779248857736053`*^9}]
}, Open  ]]
}, Closed]],

Cell[CellGroupData[{

Cell["Computing decays", "Section",
 CellChangeTimes->{{3.5781006310351267`*^9, 3.578100633256648*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"decays", "=", 
   RowBox[{"ComputeWidths", "[", "vertices", "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.578100182561545*^9, 3.57810020224467*^9}, 
   3.5781019603591213`*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Flavor expansion of the vertices: \"\>", "\[InvisibleSpace]", 
   DynamicBox[ToBoxes[FeynRules`FR$Count1, StandardForm],
    ImageSizeCache->{22., {1., 12.}}], "\[InvisibleSpace]", "\<\" / \"\>", 
   "\[InvisibleSpace]", "13"}],
  SequenceForm["Flavor expansion of the vertices: ", 
   Dynamic[FeynRules`FR$Count1], " / ", 13],
  Editable->False]], "Print",
 CellChangeTimes->{3.5799559328896837`*^9, 3.6142270008489513`*^9, 
  3.6779251765585423`*^9}],

Cell[BoxData[
 StyleBox["\<\"Computing the squared matrix elements relevant for the 1->2 \
decays: \"\>",
  StripOnInput->False,
  LineColor->RGBColor[1, 0.5, 0],
  FrontFaceColor->RGBColor[1, 0.5, 0],
  BackFaceColor->RGBColor[1, 0.5, 0],
  GraphicsColor->RGBColor[1, 0.5, 0],
  FontWeight->Bold,
  FontColor->RGBColor[1, 0.5, 0]]], "Print",
 CellChangeTimes->{3.5799559328896837`*^9, 3.6142270008489513`*^9, 
  3.6779251766884413`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
   DynamicBox[ToBoxes[FeynRules`FR$DecayCounter, StandardForm],
    ImageSizeCache->{22., {1., 12.}}], "\[InvisibleSpace]", "\<\" / \"\>", 
   "\[InvisibleSpace]", "10"}],
  SequenceForm[
   Dynamic[FeynRules`FR$DecayCounter], " / ", 10],
  Editable->False]], "Print",
 CellChangeTimes->{3.5799559328896837`*^9, 3.6142270008489513`*^9, 
  3.677925176692851*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"x1", "=", 
  RowBox[{"PartialWidth", "[", 
   RowBox[{"{", 
    RowBox[{"uv", ",", "u", ",", "p1"}], "}"}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"x2", "=", 
  RowBox[{"PartialWidth", "[", 
   RowBox[{"{", 
    RowBox[{"uv", ",", "u", ",", "p2"}], "}"}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"x3", "=", 
  RowBox[{"PartialWidth", "[", 
   RowBox[{"{", 
    RowBox[{"p2", ",", "e", ",", "evbar"}], "}"}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"x4", "=", 
  RowBox[{"PartialWidth", "[", 
   RowBox[{"{", 
    RowBox[{"ev", ",", "e", ",", "p1"}], "}"}], "]"}]}]}], "Input",
 CellChangeTimes->{{3.578100256478664*^9, 3.578100389458416*^9}, {
  3.578100496453477*^9, 3.578100496643045*^9}, {3.5781006406864853`*^9, 
  3.578100657875001*^9}, {3.578101047163535*^9, 3.578101063202585*^9}, {
  3.578101099217561*^9, 3.578101110441593*^9}, {3.579955894464841*^9, 
  3.579955895661914*^9}, {3.677925189660521*^9, 3.677925191006901*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{"MPe1", "-", "Muv"}], ")"}], "2"], " ", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{"MPe1", "+", "Muv"}], ")"}], "2"], " ", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"lam2", " ", 
       RowBox[{"Cos", "[", "th", "]"}]}], "-", 
      RowBox[{"lam1", " ", 
       RowBox[{"Sin", "[", "th", "]"}]}]}], ")"}], "2"]}], 
  RowBox[{"32", " ", "\[Pi]", " ", 
   SuperscriptBox[
    RowBox[{"Abs", "[", "Muv", "]"}], "3"]}]]], "Output",
 CellChangeTimes->{{3.578100267712432*^9, 3.578100306961565*^9}, {
   3.578100337491872*^9, 3.578100389938218*^9}, 3.5781004969919157`*^9, {
   3.578100649784704*^9, 3.5781006582507563`*^9}, 3.5781010640586987`*^9, 
   3.578101111331161*^9, 3.5793596917372513`*^9, 3.579882886900486*^9, 
   3.579883025786955*^9, 3.5799558989419003`*^9, 3.579955935260817*^9, 
   3.6142270029337683`*^9, {3.677925177929203*^9, 3.677925192406869*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{"MPe2", "-", "Muv"}], ")"}], "2"], " ", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{"MPe2", "+", "Muv"}], ")"}], "2"], " ", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"lam1", " ", 
       RowBox[{"Cos", "[", "th", "]"}]}], "+", 
      RowBox[{"lam2", " ", 
       RowBox[{"Sin", "[", "th", "]"}]}]}], ")"}], "2"]}], 
  RowBox[{"32", " ", "\[Pi]", " ", 
   SuperscriptBox[
    RowBox[{"Abs", "[", "Muv", "]"}], "3"]}]]], "Output",
 CellChangeTimes->{{3.578100267712432*^9, 3.578100306961565*^9}, {
   3.578100337491872*^9, 3.578100389938218*^9}, 3.5781004969919157`*^9, {
   3.578100649784704*^9, 3.5781006582507563`*^9}, 3.5781010640586987`*^9, 
   3.578101111331161*^9, 3.5793596917372513`*^9, 3.579882886900486*^9, 
   3.579883025786955*^9, 3.5799558989419003`*^9, 3.579955935260817*^9, 
   3.6142270029337683`*^9, {3.677925177929203*^9, 3.6779251924122553`*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{"Mev", "-", "MPe2"}], ")"}], "2"], " ", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{"Mev", "+", "MPe2"}], ")"}], "2"], " ", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"lam1p", " ", 
       RowBox[{"Cos", "[", "th", "]"}]}], "+", 
      RowBox[{"lam2p", " ", 
       RowBox[{"Sin", "[", "th", "]"}]}]}], ")"}], "2"]}], 
  RowBox[{"16", " ", "\[Pi]", " ", 
   SuperscriptBox[
    RowBox[{"Abs", "[", "MPe2", "]"}], "3"]}]]], "Output",
 CellChangeTimes->{{3.578100267712432*^9, 3.578100306961565*^9}, {
   3.578100337491872*^9, 3.578100389938218*^9}, 3.5781004969919157`*^9, {
   3.578100649784704*^9, 3.5781006582507563`*^9}, 3.5781010640586987`*^9, 
   3.578101111331161*^9, 3.5793596917372513`*^9, 3.579882886900486*^9, 
   3.579883025786955*^9, 3.5799558989419003`*^9, 3.579955935260817*^9, 
   3.6142270029337683`*^9, {3.677925177929203*^9, 3.677925192414914*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{"Mev", "-", "MPe1"}], ")"}], "2"], " ", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{"Mev", "+", "MPe1"}], ")"}], "2"], " ", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"lam2p", " ", 
       RowBox[{"Cos", "[", "th", "]"}]}], "-", 
      RowBox[{"lam1p", " ", 
       RowBox[{"Sin", "[", "th", "]"}]}]}], ")"}], "2"]}], 
  RowBox[{"32", " ", "\[Pi]", " ", 
   SuperscriptBox[
    RowBox[{"Abs", "[", "Mev", "]"}], "3"]}]]], "Output",
 CellChangeTimes->{{3.578100267712432*^9, 3.578100306961565*^9}, {
   3.578100337491872*^9, 3.578100389938218*^9}, 3.5781004969919157`*^9, {
   3.578100649784704*^9, 3.5781006582507563`*^9}, 3.5781010640586987`*^9, 
   3.578101111331161*^9, 3.5793596917372513`*^9, 3.579882886900486*^9, 
   3.579883025786955*^9, 3.5799558989419003`*^9, 3.579955935260817*^9, 
   3.6142270029337683`*^9, {3.677925177929203*^9, 3.6779251924173937`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"NumericalValue", "[", "x1", "]"}], "\[IndentingNewLine]", 
 RowBox[{"NumericalValue", "[", "x2", "]"}], "\[IndentingNewLine]", 
 RowBox[{"NumericalValue", "[", "x3", "]"}], "\[IndentingNewLine]", 
 RowBox[{"NumericalValue", "[", "x4", "]"}]}], "Input",
 CellChangeTimes->{{3.5781006598477917`*^9, 3.5781006665443563`*^9}, {
  3.578101066979349*^9, 3.5781010707919703`*^9}, {3.578101188434381*^9, 
  3.57810118966578*^9}}],

Cell[BoxData["0.03978624880756019`"], "Output",
 CellChangeTimes->{3.57810066702802*^9, 3.578101071234214*^9, 
  3.578101189891306*^9, 3.579359698611641*^9, 3.579882889407324*^9, 
  3.5798830277150307`*^9, 3.579955939561494*^9, 3.614227006142589*^9, 
  3.677925196435122*^9}],

Cell[BoxData["0.0349723172508463`"], "Output",
 CellChangeTimes->{3.57810066702802*^9, 3.578101071234214*^9, 
  3.578101189891306*^9, 3.579359698611641*^9, 3.579882889407324*^9, 
  3.5798830277150307`*^9, 3.579955939561494*^9, 3.614227006142589*^9, 
  3.677925196439011*^9}],

Cell[BoxData["0.01119114152936453`"], "Output",
 CellChangeTimes->{3.57810066702802*^9, 3.578101071234214*^9, 
  3.578101189891306*^9, 3.579359698611641*^9, 3.579882889407324*^9, 
  3.5798830277150307`*^9, 3.579955939561494*^9, 3.614227006142589*^9, 
  3.677925196441875*^9}],

Cell[BoxData["0.004969365413133469`"], "Output",
 CellChangeTimes->{3.57810066702802*^9, 3.578101071234214*^9, 
  3.578101189891306*^9, 3.579359698611641*^9, 3.579882889407324*^9, 
  3.5798830277150307`*^9, 3.579955939561494*^9, 3.614227006142589*^9, 
  3.6779251964447193`*^9}]
}, Open  ]]
}, Closed]],

Cell[CellGroupData[{

Cell["UFO", "Section",
 CellChangeTimes->{{3.539319024076879*^9, 3.53931903003391*^9}, {
  3.578100803596233*^9, 3.578100804390173*^9}, {3.579359706634882*^9, 
  3.579359708801149*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"SetDirectory", "[", "\"\<~/mg5amcnlo/models\>\"", 
  "]"}], "\[IndentingNewLine]", 
 RowBox[{"WriteUFO", "[", 
  RowBox[{
   RowBox[{"LSM", "+", "LNew"}], ",", 
   RowBox[{"Output", "\[Rule]", "\"\<MC4BSM\>\""}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.539319031354458*^9, 3.53931903608384*^9}, {
   3.579359719778757*^9, 3.579359744234468*^9}, {3.579361601368628*^9, 
   3.579361603527499*^9}, {3.579883038689376*^9, 3.579883070017232*^9}, {
   3.5798833617114*^9, 3.5798833680620728`*^9}, {3.579956090474226*^9, 
   3.5799561181857033`*^9}, {3.614227032838636*^9, 3.614227052300777*^9}, 
   3.614227254605687*^9, {3.617162208459701*^9, 3.6171622342337923`*^9}, {
   3.677925868505498*^9, 3.677925871478477*^9}}],

Cell[BoxData["\<\"/Users/degrande/mg5amcnlo/models\"\>"], "Output",
 CellChangeTimes->{3.6142270562894497`*^9, 3.614227255614313*^9, 
  3.617162237645134*^9}],

Cell[CellGroupData[{

Cell[BoxData["\<\" --- Universal FeynRules Output (UFO) v 1.1 ---\"\>"], \
"Print",
 CellChangeTimes->{3.539319036779101*^9, 3.555944767233947*^9, 
  3.5773734727639713`*^9, 3.577374315202577*^9, 3.578100794686623*^9, 
  3.579359751448289*^9, 3.57936111383105*^9, 3.579361316995659*^9, 
  3.579883047778696*^9, 3.579956098171582*^9, 3.614227057826811*^9, 
  3.614227257151779*^9, 3.617162239164522*^9}],

Cell[BoxData[
 StyleBox["\<\"Starting Feynman rule calculation.\"\>",
  StripOnInput->False,
  LineColor->RGBColor[1, 0.5, 0],
  FrontFaceColor->RGBColor[1, 0.5, 0],
  BackFaceColor->RGBColor[1, 0.5, 0],
  GraphicsColor->RGBColor[1, 0.5, 0],
  FontWeight->Bold,
  FontColor->RGBColor[1, 0.5, 0]]], "Print",
 CellChangeTimes->{3.539319036779101*^9, 3.555944767233947*^9, 
  3.5773734727639713`*^9, 3.577374315202577*^9, 3.578100794686623*^9, 
  3.579359751448289*^9, 3.57936111383105*^9, 3.579361316995659*^9, 
  3.579883047778696*^9, 3.579956098171582*^9, 3.614227057826811*^9, 
  3.614227257151779*^9, 3.617162239637867*^9}],

Cell[BoxData["\<\"Expanding the Lagrangian...\"\>"], "Print",
 CellChangeTimes->{3.539319036779101*^9, 3.555944767233947*^9, 
  3.5773734727639713`*^9, 3.577374315202577*^9, 3.578100794686623*^9, 
  3.579359751448289*^9, 3.57936111383105*^9, 3.579361316995659*^9, 
  3.579883047778696*^9, 3.579956098171582*^9, 3.614227057826811*^9, 
  3.614227257151779*^9, 3.617162239638989*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Expanding indices over \"\>", "\[InvisibleSpace]", "4", 
   "\[InvisibleSpace]", "\<\" cores\"\>"}],
  SequenceForm["Expanding indices over ", 4, " cores"],
  Editable->False]], "Print",
 CellChangeTimes->{3.539319036779101*^9, 3.555944767233947*^9, 
  3.5773734727639713`*^9, 3.577374315202577*^9, 3.578100794686623*^9, 
  3.579359751448289*^9, 3.57936111383105*^9, 3.579361316995659*^9, 
  3.579883047778696*^9, 3.579956098171582*^9, 3.614227057826811*^9, 
  3.614227257151779*^9, 3.617162239639653*^9}],

Cell[BoxData["\<\"Collecting the different structures that enter the \
vertex.\"\>"], "Print",
 CellChangeTimes->{3.539319036779101*^9, 3.555944767233947*^9, 
  3.5773734727639713`*^9, 3.577374315202577*^9, 3.578100794686623*^9, 
  3.579359751448289*^9, 3.57936111383105*^9, 3.579361316995659*^9, 
  3.579883047778696*^9, 3.579956098171582*^9, 3.614227057826811*^9, 
  3.614227257151779*^9, 3.6171622421983967`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
  "111", "\[InvisibleSpace]", "\<\" possible non-zero vertices have been \
found -> starting the computation: \"\>", "\[InvisibleSpace]", 
   DynamicBox[ToBoxes[FeynRules`FR$FeynmanRules, StandardForm],
    ImageSizeCache->{162., {4., 12.}}], "\[InvisibleSpace]", "\<\" / \"\>", 
   "\[InvisibleSpace]", "111", "\[InvisibleSpace]", "\<\".\"\>"}],
  SequenceForm[
  111, " possible non-zero vertices have been found -> starting the \
computation: ", 
   Dynamic[FeynRules`FR$FeynmanRules], " / ", 111, "."],
  Editable->False]], "Print",
 CellChangeTimes->{3.539319036779101*^9, 3.555944767233947*^9, 
  3.5773734727639713`*^9, 3.577374315202577*^9, 3.578100794686623*^9, 
  3.579359751448289*^9, 3.57936111383105*^9, 3.579361316995659*^9, 
  3.579883047778696*^9, 3.579956098171582*^9, 3.614227057826811*^9, 
  3.614227257151779*^9, 3.617162242246779*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"106", "\[InvisibleSpace]", "\<\" vertices obtained.\"\>"}],
  SequenceForm[106, " vertices obtained."],
  Editable->False]], "Print",
 CellChangeTimes->{3.539319036779101*^9, 3.555944767233947*^9, 
  3.5773734727639713`*^9, 3.577374315202577*^9, 3.578100794686623*^9, 
  3.579359751448289*^9, 3.57936111383105*^9, 3.579361316995659*^9, 
  3.579883047778696*^9, 3.579956098171582*^9, 3.614227057826811*^9, 
  3.614227257151779*^9, 3.617162244894644*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Flavor expansion of the vertices distributed over \"\>", 
   "\[InvisibleSpace]", "4", "\[InvisibleSpace]", "\<\" cores: \"\>", 
   "\[InvisibleSpace]", 
   DynamicBox[ToBoxes[FeynRules`FR$Count1, StandardForm],
    ImageSizeCache->{97., {3., 12.}}], "\[InvisibleSpace]", "\<\" / \"\>", 
   "\[InvisibleSpace]", "106"}],
  SequenceForm[
  "Flavor expansion of the vertices distributed over ", 4, " cores: ", 
   Dynamic[FeynRules`FR$Count1], " / ", 106],
  Editable->False]], "Print",
 CellChangeTimes->{3.539319036779101*^9, 3.555944767233947*^9, 
  3.5773734727639713`*^9, 3.577374315202577*^9, 3.578100794686623*^9, 
  3.579359751448289*^9, 3.57936111383105*^9, 3.579361316995659*^9, 
  3.579883047778696*^9, 3.579956098171582*^9, 3.614227057826811*^9, 
  3.614227257151779*^9, 3.6171622469681997`*^9}],

Cell[BoxData["\<\"   - Saved vertices in InterfaceRun[ 1 ].\"\>"], "Print",
 CellChangeTimes->{3.539319036779101*^9, 3.555944767233947*^9, 
  3.5773734727639713`*^9, 3.577374315202577*^9, 3.578100794686623*^9, 
  3.579359751448289*^9, 3.57936111383105*^9, 3.579361316995659*^9, 
  3.579883047778696*^9, 3.579956098171582*^9, 3.614227057826811*^9, 
  3.614227257151779*^9, 3.617162249752763*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"QN", "::", "NonConserv"}], "MessageName"], ":", 
  " ", "\<\"Warning: non quantum number conserving vertex \
encountered!\"\>"}]], "Message", "MSG",
 CellChangeTimes->{3.579359768176352*^9, 3.579361130076136*^9, 
  3.5793613334515333`*^9, 3.579883044648332*^9, 3.579956113866609*^9, 
  3.614227056284987*^9, 3.6142272675257683`*^9, 3.617162249862019*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Quantum number \"\>", "\[InvisibleSpace]", "LeptonNumber", 
   "\[InvisibleSpace]", "\<\" not conserved in vertex \"\>", 
   "\[InvisibleSpace]", 
   RowBox[{"{", 
    RowBox[{"p2", ",", 
     OverscriptBox["e", "\<\"-\"\>"], ",", "ev"}], "}"}], 
   "\[InvisibleSpace]", "\<\".\"\>"}],
  SequenceForm[
  "Quantum number ", $CellContext`LeptonNumber, 
   " not conserved in vertex ", {$CellContext`p2, $CellContext`ebar, \
$CellContext`ev}, "."],
  Editable->False]], "Print",
 CellChangeTimes->{3.539319036779101*^9, 3.555944767233947*^9, 
  3.5773734727639713`*^9, 3.577374315202577*^9, 3.578100794686623*^9, 
  3.579359751448289*^9, 3.57936111383105*^9, 3.579361316995659*^9, 
  3.579883047778696*^9, 3.579956098171582*^9, 3.614227057826811*^9, 
  3.614227257151779*^9, 3.617162249863594*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"QN", "::", "NonConserv"}], "MessageName"], ":", 
  " ", "\<\"Warning: non quantum number conserving vertex \
encountered!\"\>"}]], "Message", "MSG",
 CellChangeTimes->{3.579359768176352*^9, 3.579361130076136*^9, 
  3.5793613334515333`*^9, 3.579883044648332*^9, 3.579956113866609*^9, 
  3.614227056284987*^9, 3.6142272675257683`*^9, 3.617162249961705*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Quantum number \"\>", "\[InvisibleSpace]", "Y", 
   "\[InvisibleSpace]", "\<\" not conserved in vertex \"\>", 
   "\[InvisibleSpace]", 
   RowBox[{"{", 
    RowBox[{"p2", ",", 
     OverscriptBox["e", "\<\"-\"\>"], ",", "ev"}], "}"}], 
   "\[InvisibleSpace]", "\<\".\"\>"}],
  SequenceForm[
  "Quantum number ", FeynRules`Y, 
   " not conserved in vertex ", {$CellContext`p2, $CellContext`ebar, \
$CellContext`ev}, "."],
  Editable->False]], "Print",
 CellChangeTimes->{3.539319036779101*^9, 3.555944767233947*^9, 
  3.5773734727639713`*^9, 3.577374315202577*^9, 3.578100794686623*^9, 
  3.579359751448289*^9, 3.57936111383105*^9, 3.579361316995659*^9, 
  3.579883047778696*^9, 3.579956098171582*^9, 3.614227057826811*^9, 
  3.614227257151779*^9, 3.617162249962639*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"QN", "::", "NonConserv"}], "MessageName"], ":", 
  " ", "\<\"Warning: non quantum number conserving vertex \
encountered!\"\>"}]], "Message", "MSG",
 CellChangeTimes->{3.579359768176352*^9, 3.579361130076136*^9, 
  3.5793613334515333`*^9, 3.579883044648332*^9, 3.579956113866609*^9, 
  3.614227056284987*^9, 3.6142272675257683`*^9, 3.6171622500555143`*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"General", "::", "stop"}], "MessageName"], ":", 
  " ", "\<\"Further output of \[NoBreak]\\!\\(\\*StyleBox[\\(QN :: NonConserv\
\\), \\\"MessageName\\\"]\\)\[NoBreak] will be suppressed during this \
calculation. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \
ButtonStyle->\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/General/stop\\\", ButtonNote -> \
\\\"General::stop\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{3.579359768176352*^9, 3.579361130076136*^9, 
  3.5793613334515333`*^9, 3.579883044648332*^9, 3.579956113866609*^9, 
  3.614227056284987*^9, 3.6142272675257683`*^9, 3.617162250080212*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Quantum number \"\>", "\[InvisibleSpace]", "LeptonNumber", 
   "\[InvisibleSpace]", "\<\" not conserved in vertex \"\>", 
   "\[InvisibleSpace]", 
   RowBox[{"{", 
    RowBox[{"p1", ",", 
     OverscriptBox["e", "\<\"-\"\>"], ",", "ev"}], "}"}], 
   "\[InvisibleSpace]", "\<\".\"\>"}],
  SequenceForm[
  "Quantum number ", $CellContext`LeptonNumber, 
   " not conserved in vertex ", {$CellContext`p1, $CellContext`ebar, \
$CellContext`ev}, "."],
  Editable->False]], "Print",
 CellChangeTimes->{3.539319036779101*^9, 3.555944767233947*^9, 
  3.5773734727639713`*^9, 3.577374315202577*^9, 3.578100794686623*^9, 
  3.579359751448289*^9, 3.57936111383105*^9, 3.579361316995659*^9, 
  3.579883047778696*^9, 3.579956098171582*^9, 3.614227057826811*^9, 
  3.614227257151779*^9, 3.617162250081398*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Quantum number \"\>", "\[InvisibleSpace]", "Y", 
   "\[InvisibleSpace]", "\<\" not conserved in vertex \"\>", 
   "\[InvisibleSpace]", 
   RowBox[{"{", 
    RowBox[{"p1", ",", 
     OverscriptBox["e", "\<\"-\"\>"], ",", "ev"}], "}"}], 
   "\[InvisibleSpace]", "\<\".\"\>"}],
  SequenceForm[
  "Quantum number ", FeynRules`Y, 
   " not conserved in vertex ", {$CellContext`p1, $CellContext`ebar, \
$CellContext`ev}, "."],
  Editable->False]], "Print",
 CellChangeTimes->{3.539319036779101*^9, 3.555944767233947*^9, 
  3.5773734727639713`*^9, 3.577374315202577*^9, 3.578100794686623*^9, 
  3.579359751448289*^9, 3.57936111383105*^9, 3.579361316995659*^9, 
  3.579883047778696*^9, 3.579956098171582*^9, 3.614227057826811*^9, 
  3.614227257151779*^9, 3.6171622500821667`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Quantum number \"\>", "\[InvisibleSpace]", "Y", 
   "\[InvisibleSpace]", "\<\" not conserved in vertex \"\>", 
   "\[InvisibleSpace]", 
   RowBox[{"{", 
    RowBox[{"p2", ",", 
     OverscriptBox["u", "\<\"-\"\>"], ",", "uv"}], "}"}], 
   "\[InvisibleSpace]", "\<\".\"\>"}],
  SequenceForm[
  "Quantum number ", FeynRules`Y, 
   " not conserved in vertex ", {$CellContext`p2, $CellContext`ubar, \
$CellContext`uv}, "."],
  Editable->False]], "Print",
 CellChangeTimes->{3.539319036779101*^9, 3.555944767233947*^9, 
  3.5773734727639713`*^9, 3.577374315202577*^9, 3.578100794686623*^9, 
  3.579359751448289*^9, 3.57936111383105*^9, 3.579361316995659*^9, 
  3.579883047778696*^9, 3.579956098171582*^9, 3.614227057826811*^9, 
  3.614227257151779*^9, 3.617162250083528*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Quantum number \"\>", "\[InvisibleSpace]", "Y", 
   "\[InvisibleSpace]", "\<\" not conserved in vertex \"\>", 
   "\[InvisibleSpace]", 
   RowBox[{"{", 
    RowBox[{"p1", ",", 
     OverscriptBox["u", "\<\"-\"\>"], ",", "uv"}], "}"}], 
   "\[InvisibleSpace]", "\<\".\"\>"}],
  SequenceForm[
  "Quantum number ", FeynRules`Y, 
   " not conserved in vertex ", {$CellContext`p1, $CellContext`ubar, \
$CellContext`uv}, "."],
  Editable->False]], "Print",
 CellChangeTimes->{3.539319036779101*^9, 3.555944767233947*^9, 
  3.5773734727639713`*^9, 3.577374315202577*^9, 3.578100794686623*^9, 
  3.579359751448289*^9, 3.57936111383105*^9, 3.579361316995659*^9, 
  3.579883047778696*^9, 3.579956098171582*^9, 3.614227057826811*^9, 
  3.614227257151779*^9, 3.617162250084244*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Quantum number \"\>", "\[InvisibleSpace]", "LeptonNumber", 
   "\[InvisibleSpace]", "\<\" not conserved in vertex \"\>", 
   "\[InvisibleSpace]", 
   RowBox[{"{", 
    RowBox[{"p1", ",", 
     OverscriptBox["ev", "\<\"-\"\>"], ",", "e"}], "}"}], 
   "\[InvisibleSpace]", "\<\".\"\>"}],
  SequenceForm[
  "Quantum number ", $CellContext`LeptonNumber, 
   " not conserved in vertex ", {$CellContext`p1, $CellContext`evbar, 
    FeynRules`e}, "."],
  Editable->False]], "Print",
 CellChangeTimes->{3.539319036779101*^9, 3.555944767233947*^9, 
  3.5773734727639713`*^9, 3.577374315202577*^9, 3.578100794686623*^9, 
  3.579359751448289*^9, 3.57936111383105*^9, 3.579361316995659*^9, 
  3.579883047778696*^9, 3.579956098171582*^9, 3.614227057826811*^9, 
  3.614227257151779*^9, 3.617162250084955*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Quantum number \"\>", "\[InvisibleSpace]", "Y", 
   "\[InvisibleSpace]", "\<\" not conserved in vertex \"\>", 
   "\[InvisibleSpace]", 
   RowBox[{"{", 
    RowBox[{"p1", ",", 
     OverscriptBox["ev", "\<\"-\"\>"], ",", "e"}], "}"}], 
   "\[InvisibleSpace]", "\<\".\"\>"}],
  SequenceForm[
  "Quantum number ", FeynRules`Y, 
   " not conserved in vertex ", {$CellContext`p1, $CellContext`evbar, 
    FeynRules`e}, "."],
  Editable->False]], "Print",
 CellChangeTimes->{3.539319036779101*^9, 3.555944767233947*^9, 
  3.5773734727639713`*^9, 3.577374315202577*^9, 3.578100794686623*^9, 
  3.579359751448289*^9, 3.57936111383105*^9, 3.579361316995659*^9, 
  3.579883047778696*^9, 3.579956098171582*^9, 3.614227057826811*^9, 
  3.614227257151779*^9, 3.61716225008567*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Quantum number \"\>", "\[InvisibleSpace]", "LeptonNumber", 
   "\[InvisibleSpace]", "\<\" not conserved in vertex \"\>", 
   "\[InvisibleSpace]", 
   RowBox[{"{", 
    RowBox[{"p2", ",", 
     OverscriptBox["ev", "\<\"-\"\>"], ",", "e"}], "}"}], 
   "\[InvisibleSpace]", "\<\".\"\>"}],
  SequenceForm[
  "Quantum number ", $CellContext`LeptonNumber, 
   " not conserved in vertex ", {$CellContext`p2, $CellContext`evbar, 
    FeynRules`e}, "."],
  Editable->False]], "Print",
 CellChangeTimes->{3.539319036779101*^9, 3.555944767233947*^9, 
  3.5773734727639713`*^9, 3.577374315202577*^9, 3.578100794686623*^9, 
  3.579359751448289*^9, 3.57936111383105*^9, 3.579361316995659*^9, 
  3.579883047778696*^9, 3.579956098171582*^9, 3.614227057826811*^9, 
  3.614227257151779*^9, 3.617162250086384*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Quantum number \"\>", "\[InvisibleSpace]", "Y", 
   "\[InvisibleSpace]", "\<\" not conserved in vertex \"\>", 
   "\[InvisibleSpace]", 
   RowBox[{"{", 
    RowBox[{"p2", ",", 
     OverscriptBox["ev", "\<\"-\"\>"], ",", "e"}], "}"}], 
   "\[InvisibleSpace]", "\<\".\"\>"}],
  SequenceForm[
  "Quantum number ", FeynRules`Y, 
   " not conserved in vertex ", {$CellContext`p2, $CellContext`evbar, 
    FeynRules`e}, "."],
  Editable->False]], "Print",
 CellChangeTimes->{3.539319036779101*^9, 3.555944767233947*^9, 
  3.5773734727639713`*^9, 3.577374315202577*^9, 3.578100794686623*^9, 
  3.579359751448289*^9, 3.57936111383105*^9, 3.579361316995659*^9, 
  3.579883047778696*^9, 3.579956098171582*^9, 3.614227057826811*^9, 
  3.614227257151779*^9, 3.6171622500871468`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Quantum number \"\>", "\[InvisibleSpace]", "Y", 
   "\[InvisibleSpace]", "\<\" not conserved in vertex \"\>", 
   "\[InvisibleSpace]", 
   RowBox[{"{", 
    RowBox[{"p1", ",", 
     OverscriptBox["uv", "\<\"-\"\>"], ",", "u"}], "}"}], 
   "\[InvisibleSpace]", "\<\".\"\>"}],
  SequenceForm[
  "Quantum number ", FeynRules`Y, 
   " not conserved in vertex ", {$CellContext`p1, $CellContext`uvbar, 
    FeynRules`u}, "."],
  Editable->False]], "Print",
 CellChangeTimes->{3.539319036779101*^9, 3.555944767233947*^9, 
  3.5773734727639713`*^9, 3.577374315202577*^9, 3.578100794686623*^9, 
  3.579359751448289*^9, 3.57936111383105*^9, 3.579361316995659*^9, 
  3.579883047778696*^9, 3.579956098171582*^9, 3.614227057826811*^9, 
  3.614227257151779*^9, 3.61716225008786*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Quantum number \"\>", "\[InvisibleSpace]", "Y", 
   "\[InvisibleSpace]", "\<\" not conserved in vertex \"\>", 
   "\[InvisibleSpace]", 
   RowBox[{"{", 
    RowBox[{"p2", ",", 
     OverscriptBox["uv", "\<\"-\"\>"], ",", "u"}], "}"}], 
   "\[InvisibleSpace]", "\<\".\"\>"}],
  SequenceForm[
  "Quantum number ", FeynRules`Y, 
   " not conserved in vertex ", {$CellContext`p2, $CellContext`uvbar, 
    FeynRules`u}, "."],
  Editable->False]], "Print",
 CellChangeTimes->{3.539319036779101*^9, 3.555944767233947*^9, 
  3.5773734727639713`*^9, 3.577374315202577*^9, 3.578100794686623*^9, 
  3.579359751448289*^9, 3.57936111383105*^9, 3.579361316995659*^9, 
  3.579883047778696*^9, 3.579956098171582*^9, 3.614227057826811*^9, 
  3.614227257151779*^9, 3.6171622500887947`*^9}],

Cell[BoxData[
 StyleBox["\<\"Computing the squared matrix elements relevant for the 1->2 \
decays: \"\>",
  StripOnInput->False,
  LineColor->RGBColor[1, 0.5, 0],
  FrontFaceColor->RGBColor[1, 0.5, 0],
  BackFaceColor->RGBColor[1, 0.5, 0],
  GraphicsColor->RGBColor[1, 0.5, 0],
  FontWeight->Bold,
  FontColor->RGBColor[1, 0.5, 0]]], "Print",
 CellChangeTimes->{3.539319036779101*^9, 3.555944767233947*^9, 
  3.5773734727639713`*^9, 3.577374315202577*^9, 3.578100794686623*^9, 
  3.579359751448289*^9, 3.57936111383105*^9, 3.579361316995659*^9, 
  3.579883047778696*^9, 3.579956098171582*^9, 3.614227057826811*^9, 
  3.614227257151779*^9, 3.6171622500898457`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
   DynamicBox[ToBoxes[PRIVATE`mycounter, StandardForm],
    ImageSizeCache->{184., {4., 13.}}], "\[InvisibleSpace]", "\<\" / \"\>", 
   "\[InvisibleSpace]", "50"}],
  SequenceForm[
   Dynamic[PRIVATE`mycounter], " / ", 50],
  Editable->False]], "Print",
 CellChangeTimes->{3.539319036779101*^9, 3.555944767233947*^9, 
  3.5773734727639713`*^9, 3.577374315202577*^9, 3.578100794686623*^9, 
  3.579359751448289*^9, 3.57936111383105*^9, 3.579361316995659*^9, 
  3.579883047778696*^9, 3.579956098171582*^9, 3.614227057826811*^9, 
  3.614227257151779*^9, 3.617162250125495*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Squared matrix elent compute in \"\>", "\[InvisibleSpace]", 
   "0.65606100000000000527222709933994337916`5.837544134868329", 
   "\[InvisibleSpace]", "\<\" seconds.\"\>"}],
  SequenceForm[
  "Squared matrix elent compute in ", 
   0.65606100000000000527222709933994337916`5.837544134868329, " seconds."],
  Editable->False]], "Print",
 CellChangeTimes->{3.539319036779101*^9, 3.555944767233947*^9, 
  3.5773734727639713`*^9, 3.577374315202577*^9, 3.578100794686623*^9, 
  3.579359751448289*^9, 3.57936111383105*^9, 3.579361316995659*^9, 
  3.579883047778696*^9, 3.579956098171582*^9, 3.614227057826811*^9, 
  3.614227257151779*^9, 3.617162251133788*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
   DynamicBox[ToBoxes[PRIVATE`mycounter, StandardForm],
    ImageSizeCache->{184., {4., 13.}}], "\[InvisibleSpace]", "\<\" / \"\>", 
   "\[InvisibleSpace]", "41"}],
  SequenceForm[
   Dynamic[PRIVATE`mycounter], " / ", 41],
  Editable->False]], "Print",
 CellChangeTimes->{3.539319036779101*^9, 3.555944767233947*^9, 
  3.5773734727639713`*^9, 3.577374315202577*^9, 3.578100794686623*^9, 
  3.579359751448289*^9, 3.57936111383105*^9, 3.579361316995659*^9, 
  3.579883047778696*^9, 3.579956098171582*^9, 3.614227057826811*^9, 
  3.614227257151779*^9, 3.6171622511361303`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Decay widths computed in \"\>", "\[InvisibleSpace]", 
   "0.12118500000000000105249142734464840032`5.104048780466479", 
   "\[InvisibleSpace]", "\<\" seconds.\"\>"}],
  SequenceForm[
  "Decay widths computed in ", 
   0.12118500000000000105249142734464840032`5.104048780466479, " seconds."],
  Editable->False]], "Print",
 CellChangeTimes->{3.539319036779101*^9, 3.555944767233947*^9, 
  3.5773734727639713`*^9, 3.577374315202577*^9, 3.578100794686623*^9, 
  3.579359751448289*^9, 3.57936111383105*^9, 3.579361316995659*^9, 
  3.579883047778696*^9, 3.579956098171582*^9, 3.614227057826811*^9, 
  3.614227257151779*^9, 3.617162251306426*^9}],

Cell[BoxData["\<\"Preparing Python output.\"\>"], "Print",
 CellChangeTimes->{3.539319036779101*^9, 3.555944767233947*^9, 
  3.5773734727639713`*^9, 3.577374315202577*^9, 3.578100794686623*^9, 
  3.579359751448289*^9, 3.57936111383105*^9, 3.579361316995659*^9, 
  3.579883047778696*^9, 3.579956098171582*^9, 3.614227057826811*^9, 
  3.614227257151779*^9, 3.617162251308053*^9}],

Cell[BoxData["\<\"    - Splitting vertices into building blocks.\"\>"], \
"Print",
 CellChangeTimes->{3.539319036779101*^9, 3.555944767233947*^9, 
  3.5773734727639713`*^9, 3.577374315202577*^9, 3.578100794686623*^9, 
  3.579359751448289*^9, 3.57936111383105*^9, 3.579361316995659*^9, 
  3.579883047778696*^9, 3.579956098171582*^9, 3.614227057826811*^9, 
  3.614227257151779*^9, 3.6171622514027767`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Splitting of vertices distributed over \"\>", 
   "\[InvisibleSpace]", "4", "\[InvisibleSpace]", "\<\" kernels.\"\>"}],
  SequenceForm["Splitting of vertices distributed over ", 4, " kernels."],
  Editable->False]], "Print",
 CellChangeTimes->{3.539319036779101*^9, 3.555944767233947*^9, 
  3.5773734727639713`*^9, 3.577374315202577*^9, 3.578100794686623*^9, 
  3.579359751448289*^9, 3.57936111383105*^9, 3.579361316995659*^9, 
  3.579883047778696*^9, 3.579956098171582*^9, 3.614227057826811*^9, 
  3.614227257151779*^9, 3.6171622514272833`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"    - Optimizing: \"\>", "\[InvisibleSpace]", 
   DynamicBox[ToBoxes[PRIVATE`PY$SplitVertexCounter, StandardForm],
    ImageSizeCache->{313., {4., 13.}}], "\[InvisibleSpace]", "\<\"/\"\>", 
   "\[InvisibleSpace]", "132", "\[InvisibleSpace]", "\<\" .\"\>"}],
  SequenceForm["    - Optimizing: ", 
   Dynamic[PRIVATE`PY$SplitVertexCounter], "/", 132, " ."],
  Editable->False]], "Print",
 CellChangeTimes->{3.539319036779101*^9, 3.555944767233947*^9, 
  3.5773734727639713`*^9, 3.577374315202577*^9, 3.578100794686623*^9, 
  3.579359751448289*^9, 3.57936111383105*^9, 3.579361316995659*^9, 
  3.579883047778696*^9, 3.579956098171582*^9, 3.614227057826811*^9, 
  3.614227257151779*^9, 3.61716225161374*^9}],

Cell[BoxData["\<\"    - Writing files.\"\>"], "Print",
 CellChangeTimes->{3.539319036779101*^9, 3.555944767233947*^9, 
  3.5773734727639713`*^9, 3.577374315202577*^9, 3.578100794686623*^9, 
  3.579359751448289*^9, 3.57936111383105*^9, 3.579361316995659*^9, 
  3.579883047778696*^9, 3.579956098171582*^9, 3.614227057826811*^9, 
  3.614227257151779*^9, 3.6171622517097607`*^9}],

Cell[BoxData["\<\"Done!\"\>"], "Print",
 CellChangeTimes->{3.539319036779101*^9, 3.555944767233947*^9, 
  3.5773734727639713`*^9, 3.577374315202577*^9, 3.578100794686623*^9, 
  3.579359751448289*^9, 3.57936111383105*^9, 3.579361316995659*^9, 
  3.579883047778696*^9, 3.579956098171582*^9, 3.614227057826811*^9, 
  3.614227257151779*^9, 3.617162251845467*^9}]
}, Closed]]
}, Closed]]
}, Closed]]
}, Open  ]]
},
WindowSize->{1225, 650},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}},
ShowSelection->True,
Magnification:>FEPrivate`If[
  FEPrivate`Equal[FEPrivate`$VersionNumber, 6.], 1.5, 1.5 Inherited],
FrontEndVersion->"10.1 for Mac OS X x86 (32-bit, 64-bit Kernel) (March 23, \
2015)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[422, 15, 122, 2, 137, "Title"],
Cell[547, 19, 122, 2, 43, "Input"],
Cell[CellGroupData[{
Cell[694, 25, 104, 1, 95, "Section"],
Cell[CellGroupData[{
Cell[823, 30, 577, 10, 43, "Input"],
Cell[1403, 42, 959, 14, 43, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[2399, 61, 124, 2, 43, "Input"],
Cell[CellGroupData[{
Cell[2548, 67, 958, 14, 30, "Print"],
Cell[3509, 83, 1227, 21, 30, "Print"],
Cell[4739, 106, 1010, 15, 30, "Print"],
Cell[5752, 123, 942, 14, 30, "Print"],
Cell[6697, 139, 953, 14, 30, "Print"],
Cell[7653, 155, 1005, 15, 30, "Print"],
Cell[8661, 172, 1004, 15, 30, "Print"],
Cell[9668, 189, 944, 14, 30, "Print"],
Cell[10615, 205, 972, 14, 30, "Print"],
Cell[11590, 221, 942, 14, 30, "Print"],
Cell[12535, 237, 1009, 15, 30, "Print"]
}, Open  ]]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[13605, 259, 102, 1, 95, "Section"],
Cell[CellGroupData[{
Cell[13732, 264, 770, 11, 43, "Input"],
Cell[14505, 277, 964, 14, 43, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[15506, 296, 494, 8, 43, "Input"],
Cell[CellGroupData[{
Cell[16025, 308, 898, 12, 30, "Print"],
Cell[16926, 322, 916, 12, 30, "Print"],
Cell[17845, 336, 883, 12, 30, "Print"],
Cell[18731, 350, 1024, 16, 30, "Print"],
Cell[19758, 368, 924, 13, 30, "Print"],
Cell[20685, 383, 876, 12, 30, "Print"],
Cell[21564, 397, 906, 12, 30, "Print"],
Cell[22473, 411, 909, 12, 30, "Print"],
Cell[23385, 425, 907, 12, 30, "Print"],
Cell[24295, 439, 1076, 17, 54, "Print"]
}, Open  ]]
}, Open  ]],
Cell[25398, 460, 210, 5, 76, "Text"],
Cell[CellGroupData[{
Cell[25633, 469, 200, 4, 43, "Input"],
Cell[CellGroupData[{
Cell[25858, 477, 1251, 21, 30, "Print"],
Cell[27112, 500, 1245, 21, 30, "Print"],
Cell[28360, 523, 777, 11, 30, "Print"]
}, Open  ]]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[29198, 541, 101, 1, 95, "Section"],
Cell[CellGroupData[{
Cell[29324, 546, 930, 27, 69, "Input"],
Cell[30257, 575, 1417, 35, 69, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[31711, 615, 593, 17, 43, "Input"],
Cell[32307, 634, 2017, 51, 101, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[34361, 690, 1545, 40, 144, "Input"],
Cell[35909, 732, 1893, 49, 91, "Output"],
Cell[37805, 783, 3123, 87, 169, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[40965, 875, 178, 3, 43, "Input"],
Cell[41146, 880, 5101, 149, 276, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[46296, 1035, 208, 3, 95, "Section"],
Cell[CellGroupData[{
Cell[46529, 1042, 138, 2, 43, "Input"],
Cell[CellGroupData[{
Cell[46692, 1048, 157, 2, 30, "Print"],
Cell[46852, 1052, 153, 2, 30, "Print"],
Cell[47008, 1056, 348, 9, 30, "Print"],
Cell[47359, 1067, 103, 1, 30, "Print"],
Cell[47465, 1070, 136, 2, 30, "Print"],
Cell[47604, 1074, 94, 1, 30, "Print"],
Cell[47701, 1077, 216, 5, 30, "Print"],
Cell[47920, 1084, 103, 1, 30, "Print"]
}, Open  ]],
Cell[48038, 1088, 87, 2, 43, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[48162, 1095, 155, 2, 43, "Input"],
Cell[CellGroupData[{
Cell[48342, 1101, 327, 7, 30, "Print"],
Cell[48672, 1110, 107, 1, 30, "Print"],
Cell[48782, 1113, 119, 1, 30, "Print"]
}, Open  ]],
Cell[48916, 1117, 73, 1, 43, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[49026, 1123, 139, 2, 43, "Input"],
Cell[CellGroupData[{
Cell[49190, 1129, 663, 12, 30, "Print"],
Cell[49856, 1143, 440, 6, 30, "Print"],
Cell[50299, 1151, 434, 6, 30, "Print"],
Cell[50736, 1159, 475, 7, 30, "Print"]
}, Open  ]],
Cell[51226, 1169, 2226, 62, 240, "Output"]
}, Open  ]]
}, Closed]],
Cell[CellGroupData[{
Cell[53501, 1237, 154, 2, 72, "Section"],
Cell[CellGroupData[{
Cell[53680, 1243, 215, 4, 43, "Input"],
Cell[CellGroupData[{
Cell[53920, 1251, 636, 13, 30, "Print"],
Cell[54559, 1266, 391, 5, 30, "Print"],
Cell[54953, 1273, 422, 6, 30, "Print"],
Cell[55378, 1281, 905, 17, 54, "Print"]
}, Open  ]],
Cell[56298, 1301, 434, 10, 34, "Message"],
Cell[56735, 1313, 852, 18, 47, "Print"],
Cell[57590, 1333, 434, 10, 34, "Message"],
Cell[58027, 1345, 829, 18, 47, "Print"],
Cell[58859, 1365, 434, 10, 34, "Message"],
Cell[59296, 1377, 713, 14, 34, "Message"],
Cell[CellGroupData[{
Cell[60034, 1395, 852, 18, 47, "Print"],
Cell[60889, 1415, 829, 18, 47, "Print"],
Cell[61721, 1435, 829, 18, 46, "Print"],
Cell[62553, 1455, 827, 18, 46, "Print"],
Cell[63383, 1475, 852, 18, 47, "Print"],
Cell[64238, 1495, 829, 18, 47, "Print"],
Cell[65070, 1515, 852, 18, 47, "Print"],
Cell[65925, 1535, 825, 18, 47, "Print"],
Cell[66753, 1555, 827, 18, 46, "Print"],
Cell[67583, 1575, 829, 18, 46, "Print"],
Cell[68415, 1595, 504, 9, 30, "Print"]
}, Open  ]],
Cell[68934, 1607, 12378, 356, 538, "Output"]
}, Open  ]]
}, Closed]],
Cell[CellGroupData[{
Cell[81361, 1969, 103, 1, 72, "Section"],
Cell[CellGroupData[{
Cell[81489, 1974, 213, 5, 43, "Input"],
Cell[CellGroupData[{
Cell[81727, 1983, 501, 10, 30, "Print"],
Cell[82231, 1995, 437, 11, 30, "Print"],
Cell[82671, 2008, 406, 10, 30, "Print"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[83126, 2024, 976, 24, 119, "Input"],
Cell[84105, 2050, 982, 24, 74, "Output"],
Cell[85090, 2076, 984, 24, 74, "Output"],
Cell[86077, 2102, 985, 24, 74, "Output"],
Cell[87065, 2128, 986, 24, 74, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[88088, 2157, 446, 7, 119, "Input"],
Cell[88537, 2166, 275, 4, 43, "Output"],
Cell[88815, 2172, 274, 4, 43, "Output"],
Cell[89092, 2178, 275, 4, 43, "Output"],
Cell[89370, 2184, 278, 4, 43, "Output"]
}, Open  ]]
}, Closed]],
Cell[CellGroupData[{
Cell[89697, 2194, 185, 3, 72, "Section"],
Cell[CellGroupData[{
Cell[89907, 2201, 740, 13, 69, "Input"],
Cell[90650, 2216, 158, 2, 64, "Output"],
Cell[CellGroupData[{
Cell[90833, 2222, 402, 6, 45, "Print"],
Cell[91238, 2230, 625, 13, 45, "Print"],
Cell[91866, 2245, 380, 5, 45, "Print"],
Cell[92249, 2252, 554, 10, 45, "Print"],
Cell[92806, 2264, 415, 6, 45, "Print"],
Cell[93224, 2272, 900, 17, 81, "Print"],
Cell[94127, 2291, 497, 9, 45, "Print"],
Cell[94627, 2302, 854, 16, 81, "Print"],
Cell[95484, 2320, 394, 5, 45, "Print"]
}, Open  ]],
Cell[95893, 2328, 402, 8, 88, "Message"],
Cell[96298, 2338, 843, 18, 103, "Print"],
Cell[97144, 2358, 402, 8, 88, "Message"],
Cell[97549, 2368, 818, 18, 69, "Print"],
Cell[98370, 2388, 404, 8, 88, "Message"],
Cell[98777, 2398, 682, 12, 88, "Message"],
Cell[CellGroupData[{
Cell[99484, 2414, 843, 18, 103, "Print"],
Cell[100330, 2434, 820, 18, 69, "Print"],
Cell[101153, 2454, 818, 18, 69, "Print"],
Cell[101974, 2474, 818, 18, 69, "Print"],
Cell[102795, 2494, 843, 18, 103, "Print"],
Cell[103641, 2514, 817, 18, 69, "Print"],
Cell[104461, 2534, 843, 18, 103, "Print"],
Cell[105307, 2554, 820, 18, 69, "Print"],
Cell[106130, 2574, 817, 18, 69, "Print"],
Cell[106950, 2594, 820, 18, 69, "Print"],
Cell[107773, 2614, 663, 14, 81, "Print"],
Cell[108439, 2630, 617, 13, 45, "Print"],
Cell[109059, 2645, 703, 13, 45, "Print"],
Cell[109765, 2660, 619, 13, 45, "Print"],
Cell[110387, 2675, 689, 13, 45, "Print"],
Cell[111079, 2690, 377, 5, 45, "Print"],
Cell[111459, 2697, 403, 6, 45, "Print"],
Cell[111865, 2705, 594, 10, 45, "Print"],
Cell[112462, 2717, 752, 13, 45, "Print"],
Cell[113217, 2732, 375, 5, 45, "Print"],
Cell[113595, 2739, 358, 5, 45, "Print"]
}, Closed]]
}, Closed]]
}, Closed]]
}, Open  ]]
}
]
*)