(* ************************** *)
(* *****  Information   ***** *)
(* ************************** *)
M$ModelName = "SN";

M$Information = {
  Authors      -> {"A. Vicente"}, 
  Version      -> "1.0",
  Date         -> "27. 06. 2022",
  Institutions -> {"IFIC"},
  Emails       -> {"avelino.vicente@ific.uv.es"}
};


(* ************************** *)
(* **** Particle classes **** *)
(* ************************** *)
M$ClassesDescription = {


  F[37] == {
    ClassName        -> Neu,
    SelfConjugate    -> True,
    Mass             -> MN,
    FullName         -> "heavy",
    Width            -> 0
  },

  S[105] == {
    ClassName      -> S,
    Indices        -> {Index[SU2D], Index[Colour]},
    SelfConjugate  -> False,
    Mass           -> MS,
    FullName       -> "heavy",
    QuantumNumbers -> {Y -> 1/6}
  }

};

(* ************************** *)
(* *****   Parameters   ***** *)
(* ************************** *)
M$Parameters = {

  (* External parameters *)

  (* Internal Parameters *)
  YN == { 
    ParameterType     -> Internal, 
    Indices           -> {Index[Generation]},
    ComplexParameter  -> True
  },

  MN == { 
    ParameterType     -> Internal, 
    ComplexParameter  -> False
  },

  YS == { 
    ParameterType     -> Internal, 
    Indices           -> {Index[Generation]},
    ComplexParameter  -> True
  },

  MS == { 
    ParameterType     -> Internal, 
    ComplexParameter  -> False
  },

  lam2 == { 
    ParameterType     -> Internal, 
    ComplexParameter  -> False
  },

  lam3 == { 
    ParameterType     -> Internal, 
    ComplexParameter  -> False
  }

};



(* ************************** *)
(* *****   Lagrangian   ***** *)
(* ************************** *)

gotoBFM={G[a__]->G[a]+GQuantum[a],Wi[a__]->Wi[a]+WiQuantum[a],B[a__]->B[a]+BQuantum[a]};

LN := Block[{mu}, 
  I*(Neubar.Ga[mu].DC[Neu, mu])-MN*Neubar.Neu/2]/.gotoBFM;

LS := Block[{ii,jj,cc,mu},
  DC[Sbar[ii,cc],mu] DC[S[ii,cc],mu] - MS^2 Sbar[ii,cc] S[ii,cc]]/.gotoBFM;

LSPhi := Block[{ii,jj,cc},
  - lam2 Phibar[ii] Phi[ii] Sbar[jj,cc] S[jj,cc] -lam3 Phibar[ii] S[ii,cc] Sbar[jj,cc] Phi[jj]];

LYuk := Block[{sp1,ii,jj,cc,ff1,yuk}, 
  yuk =
   -YN[ff1] Neubar[sp1].LL[sp1,ii,ff1] Phi[jj] Eps[ii,jj]
   -YS[ff1] QLbar[sp1,ii,ff1,cc].Neu[sp1] S[ii,cc];
  yuk+HC[yuk]
 ];

LNP := LN + LS + LSPhi + LYuk;

Ltot := LSM + LNP;