B[neutrino]:=0;B[lepton]:=0;B[uquark]:=1/3;B[dquark]:=1/3;B[tquark]:=1/3;
B[bquark]:=1/3;
L[neutrino]:=1;L[lepton]:=1;L[uquark]:=0;L[dquark]:=0;L[tquark]:=0;
L[bquark]:=0;
tau3[neutrino]:=1/2;tau3[lepton]:=-1/2;tau3[uquark]:=1/2;tau3[dquark]:=-1/2;
tau3[tquark]:=1/2;tau3[bquark]:=-1/2;
P[tquark]:=1;P[bquark]:=0;

TZ1L[x_]:=-g1 v33 tau3[x] - g2 v31 (B[x]-L[x])/2;
TZ1R[x_]:=-g2 v31 (tau3[x]+(B[x]-L[x])/2);
TZ2L[x_]:=-g1/(1+b1) v33 tau3[x] - g2 v31 (B[x]-L[x])/2+
           gv v32 b1/(1+b1)tau3[x];
TZ2R[x_]:=-g2 (tau3[x]+(B[x]-L[x])/2)(1-b2/(1+b2)P[x]) v31-
          g2/2 b2/(1+b2)(B[x]-L[x])P[x]v31+v32 gv/2 b2/(1+b2)P[x];

TV1L[x_]:=-g1 v23 tau3[x] - g2 v21 (B[x]-L[x])/2;
TV1R[x_]:=-g2 v21 (tau3[x]+(B[x]-L[x])/2);
TV2L[x_]:=-g1/(1+b1) v23 tau3[x] - g2 v21 (B[x]-L[x])/2+
          gv v22 b1/(1+b1)tau3[x];
TV2R[x_]:=-g2 (tau3[x]+(B[x]-L[x])/2)(1-b2/(1+b2)P[x]) v21-
           g2/2 b2/(1+b2)(B[x]-L[x])P[x]v21+v22 gv/2 b2/(1+b2)P[x];

