%%%%%%%%             The standard model in unitary gauge     %%%%%%%%%%%

algebraic$
%------------------------- Gauge fields -------------------------------
gaugefields:=
%     No.  Name  SU(n)   notation  couping   gaugefixterm gauge_pa breaking
   '((1    u     1       b         g1        gfix1        infinit   yes    ) 
     (2    su    2       a         g         gfix2        infinit   yes     )
     (3    su    3       gs        g3        gfix3        1         no     ) 
    )$
gfix1:=pd(b(v),v)$ 
gfix2:=pd(a(i1,v),v)$ 
gfix3:=pd(gs(i1,v),v)$ 
%----------------------------------------------------------------------

%------------------------ Matter fields -------------------------------- 
matterinput:={{name,  2*spin,  chiral, "|g",  1,     2,    3}, 
              { hig,   0  ,    rl  ,          1,     2,    0},
%              { hig1,  0  ,    rl  ,         -1,     2,    0},
              { el ,   1,      l  ,          -1,     2,    0}, 
              { er ,   1,      r  ,          -2,     0,    0}, 
              { q1l ,   1,      l  ,          1/3,    2,    3}, 
              { q1dr,   1,      r  ,         -2/3,    0,    3},    
              { q1ur,   1,      r  ,          4/3,    0,    3}
              }$   


gbsymmetrylist:={len,lmuon,ltau};
gbsymmetrydata:={{len,  {el,1},   {er,1}},
                 {lmuon,{mul,1},  {mur,1}},
                 {ltau, {taul,1}, {taur,1}}}$  

pchalist:={{hig,h0(0)}, 
           {el,nue(0),ef(-1)}, {er,ef(-1)}, 
           {q1l,qu(2/3),qd(-1/3)}, {q1ur, qu(2/3)}, {q1dr, qd(-1/3)},
           {a,w(1),w(-1),z(0),p(0)},{b,p(0),z(0)},{gs,gs(0)}};

mixrules:={ a(1,~v) => (w(1,v)+w(-1,v))/2**0.5,
            a(2,~v) => (w(-1,v)-w(1,v))/2**0.5/i,
            a(3,~v) => z(0,v),
            b(~v)   => p(0,v) };

mixfix:={};
% where g1 can be any one of g1,g,sin(theta);
construles:={ g1=>g*sin(theta)/cos(theta) 
            }; 

vancumexpectation:='((hig v0));
bosonmix:='((sxpz p z)); 

matrixlist:={ {sxpz,{cos(theta),sin(theta)},{-sin(theta),cos(theta)}} 
            }; 

%         name  n_group   1_componet 2_componet   ...   
 mdefl:={{hig, {2,        0   ,   (v0+h0)/2**0.5}},  
%         {hig1, {2,        xx2   ,   (vx0+hx0-i*xx3)/2**0.5}},  
         {el , {2,        nue   ,   ef}}, 
         {er , {2,        ef           }}, 
         {q1l , {2,       qu   ,  qd}}, 
         {q1dr, {2,       qd          }},
         {q1ur, {2,       qu            }},  
         {qu,  {3,        qu(1),   qu(2) ,  qu(3)}},
         {qd,  {3,        qd(1)  , qd(2) ,  qd(3)}}
        }$

realfamily:='((q1l q2l q3l) (q1dr q2dr q3dr) (q1ur q2ur q3ur)
              (qu qc qt) (qd qs qb) 
              (el mul taul) (er mur taur)
              (ef mu tau) (nue numu nut)
              );

%CKM Matrix: 

mixlist:='(
%  uncomment the next line and comment out the next next line, quark mixing will be added 
%            (l (s1 (qd qs qb)) (s2 (qu qc qt)) 
            (l (s1 (qd qs qb)) (s1 (qu qc qt))
               (se1 (ef mu tau)) (se1 (nue numu nut)) 
            )
            (r (t1 (qd qs qb)) (t2 (qu qc qt))
               (te1 (ef mu tau))
            ) 
          );  

mixlistcondition:={s2(-1)*s1(1)=>v_ckm(1)};


%%%%%%%%%%% where l=(1-r_5)/2, r=(1+r_5)/2$ %%%%%%%%%%%%%%%%%%%%%%%%%%

phymass:={{w, wm}, {z, zm}, {p, 0}, {h0, hm}, 
          {gs, 0},
          {ef, fme},{nue, 0}, 
          {mu, fmmu},{numu, 0}, 
          {tau, fmtau},{nut, 0}, 
          {qd, fmd}, {qu, fmu},
          {qs, fms}, {qc, fmc},
          {qb, fmb}, {qt, fmt}
          }$

phyinput:={g3, g, sin(theta), wm, hm, 
           fme, 
           fmmu, 
           fmtau, 
           fmd, fmu,
           fms, fmc,
           fmb, fmt
          }$    


charge_def:=tp(2,3)+tp(1,1)/2;


coupling_comment:='( 
                    (g3 " strong interaction ") 
                    (g "electro-weak interaction ")
                   );

proton_like_list:='((proton qu qd qub qdb gs)); 

model_input_comment:='"This is Standard Model in unitary gauge, 
include electro-weak interaction and QCD, 
Quark mixing terms are droped."; 


;end$