Zmst1 = (1 + (as*(-(CF*Mfg^2) - CF*Mt^2 + 2*CF*cthetat*Mfg*Mt*sthetat - CF*Mst1^2*sthetat^2 + CF*Mst2^2*sthetat^2 + CF*Mst1^2*sthetat^4 - 
    CF*Mst2^2*sthetat^4))/(ep*Mst1^2) + as^2*((-3*CA*CF*Me^2)/(8*ep*Mst1^2) + (9*CA*CF*Mfg^2)/(8*ep^2*Mst1^2) - 
   (11*CA*CF*Mfg^2)/(8*ep*Mst1^2) + (3*CF^2*Mfg^2)/(4*ep*Mst1^2) + (3*CA*CF*Mt^2)/(8*ep^2*Mst1^2) + (CF^2*Mt^2)/(2*ep^2*Mst1^2) - 
   (3*CA*CF*Mt^2)/(8*ep*Mst1^2) + (CF^2*Mt^2)/(4*ep*Mst1^2) + (CF^2*Mfg^2*Mt^2)/(ep^2*Mst1^2*(Mst1^2 - Mst2^2)) - 
   (3*CA*CF*cthetat*Mfg*Mt*sthetat)/(2*ep^2*Mst1^2) - (2*CF^2*cthetat*Mfg*Mt*sthetat)/(ep^2*Mst1^2) + 
   (3*CA*CF*cthetat*Mfg*Mt*sthetat)/(2*ep*Mst1^2) - (CF^2*cthetat*Mfg*Mt*sthetat)/(ep*Mst1^2) + (3*CA*CF*sthetat^2)/(8*ep^2) + 
   (CF^2*sthetat^2)/(2*ep^2) - (3*CA*CF*sthetat^2)/(8*ep) + (CF^2*sthetat^2)/(4*ep) - (3*CA*CF*Mst2^2*sthetat^2)/(8*ep^2*Mst1^2) - 
   (CF^2*Mst2^2*sthetat^2)/(2*ep^2*Mst1^2) + (3*CA*CF*Mst2^2*sthetat^2)/(8*ep*Mst1^2) - (CF^2*Mst2^2*sthetat^2)/(4*ep*Mst1^2) - 
   (4*CF^2*Mfg^2*Mt^2*sthetat^2)/(ep^2*Mst1^2*(Mst1^2 - Mst2^2)) + (4*CF^2*cthetat*Mfg*Mt*sthetat^3)/(ep^2*Mst1^2) - 
   (3*CA*CF*sthetat^4)/(8*ep^2) - (3*CF^2*sthetat^4)/(2*ep^2) + (3*CA*CF*sthetat^4)/(8*ep) - (CF^2*sthetat^4)/(4*ep) + 
   (3*CA*CF*Mst2^2*sthetat^4)/(8*ep^2*Mst1^2) + (3*CF^2*Mst2^2*sthetat^4)/(2*ep^2*Mst1^2) - (3*CA*CF*Mst2^2*sthetat^4)/(8*ep*Mst1^2) + 
   (CF^2*Mst2^2*sthetat^4)/(4*ep*Mst1^2) + (4*CF^2*Mfg^2*Mt^2*sthetat^4)/(ep^2*Mst1^2*(Mst1^2 - Mst2^2)) - 
   (4*CF^2*cthetat*Mfg*Mt*sthetat^5)/(ep^2*Mst1^2) + (2*CF^2*sthetat^6)/ep^2 - (2*CF^2*Mst2^2*sthetat^6)/(ep^2*Mst1^2) - 
   (CF^2*sthetat^8)/ep^2 + (CF^2*Mst2^2*sthetat^8)/(ep^2*Mst1^2) + (CF*Me^2*nq*Tf)/(4*ep*Mst1^2) - (3*CF*Mfg^2*nq*Tf)/(4*ep^2*Mst1^2) + 
   (3*CF*Mfg^2*nq*Tf)/(4*ep*Mst1^2) + (CF*Msq^2*nq*Tf)/(2*ep*Mst1^2) - (CF*Mt^2*nq*Tf)/(4*ep^2*Mst1^2) + (CF*Mt^2*nq*Tf)/(4*ep*Mst1^2) + 
   (CF*nt*Tf)/(4*ep) + (CF*Me^2*nt*Tf)/(4*ep*Mst1^2) - (3*CF*Mfg^2*nt*Tf)/(4*ep^2*Mst1^2) + (3*CF*Mfg^2*nt*Tf)/(4*ep*Mst1^2) + 
   (CF*Mst2^2*nt*Tf)/(4*ep*Mst1^2) - (CF*Mt^2*nt*Tf)/(4*ep^2*Mst1^2) - (CF*Mt^2*nt*Tf)/(4*ep*Mst1^2) + 
   (CF*cthetat*Mfg*Mt*nq*sthetat*Tf)/(ep^2*Mst1^2) - (CF*cthetat*Mfg*Mt*nq*sthetat*Tf)/(ep*Mst1^2) + 
   (CF*cthetat*Mfg*Mt*nt*sthetat*Tf)/(ep^2*Mst1^2) - (CF*cthetat*Mfg*Mt*nt*sthetat*Tf)/(ep*Mst1^2) - (CF*nq*sthetat^2*Tf)/(4*ep^2) + 
   (CF*nq*sthetat^2*Tf)/(4*ep) + (CF*Mst2^2*nq*sthetat^2*Tf)/(4*ep^2*Mst1^2) - (CF*Mst2^2*nq*sthetat^2*Tf)/(4*ep*Mst1^2) - 
   (CF*nt*sthetat^2*Tf)/(4*ep^2) + (CF*nt*sthetat^2*Tf)/(4*ep) + (CF*Mst2^2*nt*sthetat^2*Tf)/(4*ep^2*Mst1^2) - 
   (CF*Mst2^2*nt*sthetat^2*Tf)/(4*ep*Mst1^2) + (CF*nq*sthetat^4*Tf)/(4*ep^2) - (CF*nq*sthetat^4*Tf)/(4*ep) - 
   (CF*Mst2^2*nq*sthetat^4*Tf)/(4*ep^2*Mst1^2) + (CF*Mst2^2*nq*sthetat^4*Tf)/(4*ep*Mst1^2) + (CF*nt*sthetat^4*Tf)/(4*ep^2) - 
   (CF*nt*sthetat^4*Tf)/(4*ep) - (CF*Mst2^2*nt*sthetat^4*Tf)/(4*ep^2*Mst1^2) + (CF*Mst2^2*nt*sthetat^4*Tf)/(4*ep*Mst1^2)) + 
 as^3*((9*CA^2*CF*Me^2)/(32*ep^2*Mst1^2) - (15*CA^2*CF*Me^2)/(32*ep*Mst1^2) + (3*CA*CF^2*Me^2)/(16*ep*Mst1^2) - 
   (9*CA^2*CF*Mfg^2)/(8*ep^3*Mst1^2) + (33*CA^2*CF*Mfg^2)/(16*ep^2*Mst1^2) - (3*CA*CF^2*Mfg^2)/(4*ep^2*Mst1^2) - 
   (5*CA^2*CF*Mfg^2)/(4*ep*Mst1^2) + (7*CA*CF^2*Mfg^2)/(8*ep*Mst1^2) - (CF^3*Mfg^2)/(ep*Mst1^2) - (3*CA^2*CF*Mt^2)/(16*ep^3*Mst1^2) - 
   (3*CA*CF^2*Mt^2)/(8*ep^3*Mst1^2) - (CF^3*Mt^2)/(6*ep^3*Mst1^2) + (5*CA^2*CF*Mt^2)/(16*ep^2*Mst1^2) + (CA*CF^2*Mt^2)/(4*ep^2*Mst1^2) - 
   (CF^3*Mt^2)/(4*ep^2*Mst1^2) - (CA^2*CF*Mt^2)/(8*ep*Mst1^2) + (CA*CF^2*Mt^2)/(8*ep*Mst1^2) - (CF^3*Mt^2)/(6*ep*Mst1^2) - 
   (3*CA*CF^2*Mfg^2*Mt^2)/(2*ep^3*Mst1^2*(Mst1^2 - Mst2^2)) - (CF^3*Mfg^2*Mt^2)/(ep^3*Mst1^2*(Mst1^2 - Mst2^2)) + 
   (3*CA*CF^2*Mfg^2*Mt^2)/(2*ep^2*Mst1^2*(Mst1^2 - Mst2^2)) - (CF^3*Mfg^2*Mt^2)/(ep^2*Mst1^2*(Mst1^2 - Mst2^2)) + 
   (9*CA^2*CF*cthetat*Mfg*Mt*sthetat)/(8*ep^3*Mst1^2) + (9*CA*CF^2*cthetat*Mfg*Mt*sthetat)/(4*ep^3*Mst1^2) + 
   (CF^3*cthetat*Mfg*Mt*sthetat)/(ep^3*Mst1^2) - (15*CA^2*CF*cthetat*Mfg*Mt*sthetat)/(8*ep^2*Mst1^2) - 
   (3*CA*CF^2*cthetat*Mfg*Mt*sthetat)/(2*ep^2*Mst1^2) + (3*CF^3*cthetat*Mfg*Mt*sthetat)/(2*ep^2*Mst1^2) + 
   (3*CA^2*CF*cthetat*Mfg*Mt*sthetat)/(4*ep*Mst1^2) - (3*CA*CF^2*cthetat*Mfg*Mt*sthetat)/(4*ep*Mst1^2) + 
   (CF^3*cthetat*Mfg*Mt*sthetat)/(ep*Mst1^2) - (4*CF^3*cthetat*Mfg^3*Mt^3*sthetat)/(ep^3*Mst1^2*(Mst1^2 - Mst2^2)^2) - 
   (3*CA^2*CF*sthetat^2)/(16*ep^3) - (3*CA*CF^2*sthetat^2)/(8*ep^3) - (CF^3*sthetat^2)/(6*ep^3) + (5*CA^2*CF*sthetat^2)/(16*ep^2) + 
   (CA*CF^2*sthetat^2)/(4*ep^2) - (CF^3*sthetat^2)/(4*ep^2) - (CA^2*CF*sthetat^2)/(8*ep) + (CA*CF^2*sthetat^2)/(8*ep) - 
   (CF^3*sthetat^2)/(6*ep) + (3*CA^2*CF*Mst2^2*sthetat^2)/(16*ep^3*Mst1^2) + (3*CA*CF^2*Mst2^2*sthetat^2)/(8*ep^3*Mst1^2) + 
   (CF^3*Mst2^2*sthetat^2)/(6*ep^3*Mst1^2) - (5*CA^2*CF*Mst2^2*sthetat^2)/(16*ep^2*Mst1^2) - (CA*CF^2*Mst2^2*sthetat^2)/(4*ep^2*Mst1^2) + 
   (CF^3*Mst2^2*sthetat^2)/(4*ep^2*Mst1^2) + (CA^2*CF*Mst2^2*sthetat^2)/(8*ep*Mst1^2) - (CA*CF^2*Mst2^2*sthetat^2)/(8*ep*Mst1^2) + 
   (CF^3*Mst2^2*sthetat^2)/(6*ep*Mst1^2) + (6*CA*CF^2*Mfg^2*Mt^2*sthetat^2)/(ep^3*Mst1^2*(Mst1^2 - Mst2^2)) + 
   (10*CF^3*Mfg^2*Mt^2*sthetat^2)/(ep^3*Mst1^2*(Mst1^2 - Mst2^2)) - (6*CA*CF^2*Mfg^2*Mt^2*sthetat^2)/(ep^2*Mst1^2*(Mst1^2 - Mst2^2)) + 
   (4*CF^3*Mfg^2*Mt^2*sthetat^2)/(ep^2*Mst1^2*(Mst1^2 - Mst2^2)) - (9*CA*CF^2*cthetat*Mfg*Mt*sthetat^3)/(2*ep^3*Mst1^2) - 
   (6*CF^3*cthetat*Mfg*Mt*sthetat^3)/(ep^3*Mst1^2) + (9*CA*CF^2*cthetat*Mfg*Mt*sthetat^3)/(2*ep^2*Mst1^2) - 
   (3*CF^3*cthetat*Mfg*Mt*sthetat^3)/(ep^2*Mst1^2) + (16*CF^3*cthetat*Mfg^3*Mt^3*sthetat^3)/(ep^3*Mst1^2*(Mst1^2 - Mst2^2)^2) + 
   (3*CA^2*CF*sthetat^4)/(16*ep^3) + (9*CA*CF^2*sthetat^4)/(8*ep^3) + (7*CF^3*sthetat^4)/(6*ep^3) - (5*CA^2*CF*sthetat^4)/(16*ep^2) - 
   (CA*CF^2*sthetat^4)/ep^2 + (3*CF^3*sthetat^4)/(4*ep^2) + (CA^2*CF*sthetat^4)/(8*ep) - (CA*CF^2*sthetat^4)/(8*ep) + 
   (CF^3*sthetat^4)/(6*ep) - (3*CA^2*CF*Mst2^2*sthetat^4)/(16*ep^3*Mst1^2) - (9*CA*CF^2*Mst2^2*sthetat^4)/(8*ep^3*Mst1^2) - 
   (7*CF^3*Mst2^2*sthetat^4)/(6*ep^3*Mst1^2) + (5*CA^2*CF*Mst2^2*sthetat^4)/(16*ep^2*Mst1^2) + (CA*CF^2*Mst2^2*sthetat^4)/(ep^2*Mst1^2) - 
   (3*CF^3*Mst2^2*sthetat^4)/(4*ep^2*Mst1^2) - (CA^2*CF*Mst2^2*sthetat^4)/(8*ep*Mst1^2) + (CA*CF^2*Mst2^2*sthetat^4)/(8*ep*Mst1^2) - 
   (CF^3*Mst2^2*sthetat^4)/(6*ep*Mst1^2) - (6*CA*CF^2*Mfg^2*Mt^2*sthetat^4)/(ep^3*Mst1^2*(Mst1^2 - Mst2^2)) - 
   (34*CF^3*Mfg^2*Mt^2*sthetat^4)/(ep^3*Mst1^2*(Mst1^2 - Mst2^2)) + (6*CA*CF^2*Mfg^2*Mt^2*sthetat^4)/(ep^2*Mst1^2*(Mst1^2 - Mst2^2)) - 
   (4*CF^3*Mfg^2*Mt^2*sthetat^4)/(ep^2*Mst1^2*(Mst1^2 - Mst2^2)) + (9*CA*CF^2*cthetat*Mfg*Mt*sthetat^5)/(2*ep^3*Mst1^2) + 
   (18*CF^3*cthetat*Mfg*Mt*sthetat^5)/(ep^3*Mst1^2) - (9*CA*CF^2*cthetat*Mfg*Mt*sthetat^5)/(2*ep^2*Mst1^2) + 
   (3*CF^3*cthetat*Mfg*Mt*sthetat^5)/(ep^2*Mst1^2) - (16*CF^3*cthetat*Mfg^3*Mt^3*sthetat^5)/(ep^3*Mst1^2*(Mst1^2 - Mst2^2)^2) - 
   (3*CA*CF^2*sthetat^6)/(2*ep^3) - (4*CF^3*sthetat^6)/ep^3 + (3*CA*CF^2*sthetat^6)/(2*ep^2) - (CF^3*sthetat^6)/ep^2 + 
   (3*CA*CF^2*Mst2^2*sthetat^6)/(2*ep^3*Mst1^2) + (4*CF^3*Mst2^2*sthetat^6)/(ep^3*Mst1^2) - (3*CA*CF^2*Mst2^2*sthetat^6)/(2*ep^2*Mst1^2) + 
   (CF^3*Mst2^2*sthetat^6)/(ep^2*Mst1^2) + (48*CF^3*Mfg^2*Mt^2*sthetat^6)/(ep^3*Mst1^2*(Mst1^2 - Mst2^2)) - 
   (24*CF^3*cthetat*Mfg*Mt*sthetat^7)/(ep^3*Mst1^2) + (3*CA*CF^2*sthetat^8)/(4*ep^3) + (7*CF^3*sthetat^8)/ep^3 - 
   (3*CA*CF^2*sthetat^8)/(4*ep^2) + (CF^3*sthetat^8)/(2*ep^2) - (3*CA*CF^2*Mst2^2*sthetat^8)/(4*ep^3*Mst1^2) - 
   (7*CF^3*Mst2^2*sthetat^8)/(ep^3*Mst1^2) + (3*CA*CF^2*Mst2^2*sthetat^8)/(4*ep^2*Mst1^2) - (CF^3*Mst2^2*sthetat^8)/(2*ep^2*Mst1^2) - 
   (24*CF^3*Mfg^2*Mt^2*sthetat^8)/(ep^3*Mst1^2*(Mst1^2 - Mst2^2)) + (12*CF^3*cthetat*Mfg*Mt*sthetat^9)/(ep^3*Mst1^2) - 
   (6*CF^3*sthetat^10)/ep^3 + (6*CF^3*Mst2^2*sthetat^10)/(ep^3*Mst1^2) + (2*CF^3*sthetat^12)/ep^3 - 
   (2*CF^3*Mst2^2*sthetat^12)/(ep^3*Mst1^2) - (3*CA*CF*Me^2*nq*Tf)/(8*ep^2*Mst1^2) + (CA*CF*Me^2*nq*Tf)/(2*ep*Mst1^2) - 
   (CF^2*Me^2*nq*Tf)/(8*ep*Mst1^2) + (3*CA*CF*Mfg^2*nq*Tf)/(2*ep^3*Mst1^2) - (17*CA*CF*Mfg^2*nq*Tf)/(8*ep^2*Mst1^2) - 
   (CF^2*Mfg^2*nq*Tf)/(2*ep^2*Mst1^2) + (9*CF^2*Mfg^2*nq*Tf)/(4*ep*Mst1^2) - (3*CA*CF*Msq^2*nq*Tf)/(8*ep^2*Mst1^2) + 
   (5*CA*CF*Msq^2*nq*Tf)/(8*ep*Mst1^2) - (CF^2*Msq^2*nq*Tf)/(4*ep*Mst1^2) + (CA*CF*Mt^2*nq*Tf)/(4*ep^3*Mst1^2) + 
   (CF^2*Mt^2*nq*Tf)/(4*ep^3*Mst1^2) - (CA*CF*Mt^2*nq*Tf)/(3*ep^2*Mst1^2) - (CF^2*Mt^2*nq*Tf)/(3*ep^2*Mst1^2) - 
   (CA*CF*Mt^2*nq*Tf)/(24*ep*Mst1^2) + (CF^2*Mt^2*nq*Tf)/(3*ep*Mst1^2) + (CF^2*Mfg^2*Mt^2*nq*Tf)/(ep^3*Mst1^2*(Mst1^2 - Mst2^2)) - 
   (CF^2*Mfg^2*Mt^2*nq*Tf)/(ep^2*Mst1^2*(Mst1^2 - Mst2^2)) - (3*CA*CF*nt*Tf)/(16*ep^2) + (5*CA*CF*nt*Tf)/(16*ep) - (CF^2*nt*Tf)/(8*ep) - 
   (3*CA*CF*Me^2*nt*Tf)/(8*ep^2*Mst1^2) + (CA*CF*Me^2*nt*Tf)/(2*ep*Mst1^2) - (CF^2*Me^2*nt*Tf)/(8*ep*Mst1^2) + 
   (3*CA*CF*Mfg^2*nt*Tf)/(2*ep^3*Mst1^2) - (17*CA*CF*Mfg^2*nt*Tf)/(8*ep^2*Mst1^2) - (CF^2*Mfg^2*nt*Tf)/(2*ep^2*Mst1^2) + 
   (9*CF^2*Mfg^2*nt*Tf)/(4*ep*Mst1^2) - (3*CA*CF*Mst2^2*nt*Tf)/(16*ep^2*Mst1^2) + (5*CA*CF*Mst2^2*nt*Tf)/(16*ep*Mst1^2) - 
   (CF^2*Mst2^2*nt*Tf)/(8*ep*Mst1^2) + (CA*CF*Mt^2*nt*Tf)/(4*ep^3*Mst1^2) + (CF^2*Mt^2*nt*Tf)/(4*ep^3*Mst1^2) + 
   (CA*CF*Mt^2*nt*Tf)/(24*ep^2*Mst1^2) - (CF^2*Mt^2*nt*Tf)/(3*ep^2*Mst1^2) - (2*CA*CF*Mt^2*nt*Tf)/(3*ep*Mst1^2) + 
   (7*CF^2*Mt^2*nt*Tf)/(12*ep*Mst1^2) + (CF^2*Mfg^2*Mt^2*nt*Tf)/(ep^3*Mst1^2*(Mst1^2 - Mst2^2)) - 
   (CF^2*Mfg^2*Mt^2*nt*Tf)/(ep^2*Mst1^2*(Mst1^2 - Mst2^2)) - (3*CA*CF*cthetat*Mfg*Mt*nq*sthetat*Tf)/(2*ep^3*Mst1^2) - 
   (3*CF^2*cthetat*Mfg*Mt*nq*sthetat*Tf)/(2*ep^3*Mst1^2) + (2*CA*CF*cthetat*Mfg*Mt*nq*sthetat*Tf)/(ep^2*Mst1^2) + 
   (2*CF^2*cthetat*Mfg*Mt*nq*sthetat*Tf)/(ep^2*Mst1^2) + (CA*CF*cthetat*Mfg*Mt*nq*sthetat*Tf)/(4*ep*Mst1^2) - 
   (2*CF^2*cthetat*Mfg*Mt*nq*sthetat*Tf)/(ep*Mst1^2) - (3*CA*CF*cthetat*Mfg*Mt*nt*sthetat*Tf)/(2*ep^3*Mst1^2) - 
   (3*CF^2*cthetat*Mfg*Mt*nt*sthetat*Tf)/(2*ep^3*Mst1^2) + (2*CA*CF*cthetat*Mfg*Mt*nt*sthetat*Tf)/(ep^2*Mst1^2) + 
   (2*CF^2*cthetat*Mfg*Mt*nt*sthetat*Tf)/(ep^2*Mst1^2) + (CA*CF*cthetat*Mfg*Mt*nt*sthetat*Tf)/(4*ep*Mst1^2) - 
   (2*CF^2*cthetat*Mfg*Mt*nt*sthetat*Tf)/(ep*Mst1^2) + (CA*CF*nq*sthetat^2*Tf)/(4*ep^3) + (CF^2*nq*sthetat^2*Tf)/(4*ep^3) - 
   (CA*CF*nq*sthetat^2*Tf)/(3*ep^2) - (CF^2*nq*sthetat^2*Tf)/(3*ep^2) - (CA*CF*nq*sthetat^2*Tf)/(24*ep) + (CF^2*nq*sthetat^2*Tf)/(3*ep) - 
   (CA*CF*Mst2^2*nq*sthetat^2*Tf)/(4*ep^3*Mst1^2) - (CF^2*Mst2^2*nq*sthetat^2*Tf)/(4*ep^3*Mst1^2) + 
   (CA*CF*Mst2^2*nq*sthetat^2*Tf)/(3*ep^2*Mst1^2) + (CF^2*Mst2^2*nq*sthetat^2*Tf)/(3*ep^2*Mst1^2) + 
   (CA*CF*Mst2^2*nq*sthetat^2*Tf)/(24*ep*Mst1^2) - (CF^2*Mst2^2*nq*sthetat^2*Tf)/(3*ep*Mst1^2) - 
   (4*CF^2*Mfg^2*Mt^2*nq*sthetat^2*Tf)/(ep^3*Mst1^2*(Mst1^2 - Mst2^2)) + (4*CF^2*Mfg^2*Mt^2*nq*sthetat^2*Tf)/
    (ep^2*Mst1^2*(Mst1^2 - Mst2^2)) + (CA*CF*nt*sthetat^2*Tf)/(4*ep^3) + (CF^2*nt*sthetat^2*Tf)/(4*ep^3) - 
   (CA*CF*nt*sthetat^2*Tf)/(3*ep^2) - (CF^2*nt*sthetat^2*Tf)/(3*ep^2) - (CA*CF*nt*sthetat^2*Tf)/(24*ep) + (CF^2*nt*sthetat^2*Tf)/(3*ep) - 
   (CA*CF*Mst2^2*nt*sthetat^2*Tf)/(4*ep^3*Mst1^2) - (CF^2*Mst2^2*nt*sthetat^2*Tf)/(4*ep^3*Mst1^2) + 
   (CA*CF*Mst2^2*nt*sthetat^2*Tf)/(3*ep^2*Mst1^2) + (CF^2*Mst2^2*nt*sthetat^2*Tf)/(3*ep^2*Mst1^2) + 
   (CA*CF*Mst2^2*nt*sthetat^2*Tf)/(24*ep*Mst1^2) - (CF^2*Mst2^2*nt*sthetat^2*Tf)/(3*ep*Mst1^2) - 
   (4*CF^2*Mfg^2*Mt^2*nt*sthetat^2*Tf)/(ep^3*Mst1^2*(Mst1^2 - Mst2^2)) + (4*CF^2*Mfg^2*Mt^2*nt*sthetat^2*Tf)/
    (ep^2*Mst1^2*(Mst1^2 - Mst2^2)) + (3*CF^2*cthetat*Mfg*Mt*nq*sthetat^3*Tf)/(ep^3*Mst1^2) - 
   (3*CF^2*cthetat*Mfg*Mt*nq*sthetat^3*Tf)/(ep^2*Mst1^2) + (3*CF^2*cthetat*Mfg*Mt*nt*sthetat^3*Tf)/(ep^3*Mst1^2) - 
   (3*CF^2*cthetat*Mfg*Mt*nt*sthetat^3*Tf)/(ep^2*Mst1^2) - (CA*CF*nq*sthetat^4*Tf)/(4*ep^3) - (3*CF^2*nq*sthetat^4*Tf)/(4*ep^3) + 
   (CA*CF*nq*sthetat^4*Tf)/(3*ep^2) + (5*CF^2*nq*sthetat^4*Tf)/(6*ep^2) + (CA*CF*nq*sthetat^4*Tf)/(24*ep) - (CF^2*nq*sthetat^4*Tf)/(3*ep) + 
   (CA*CF*Mst2^2*nq*sthetat^4*Tf)/(4*ep^3*Mst1^2) + (3*CF^2*Mst2^2*nq*sthetat^4*Tf)/(4*ep^3*Mst1^2) - 
   (CA*CF*Mst2^2*nq*sthetat^4*Tf)/(3*ep^2*Mst1^2) - (5*CF^2*Mst2^2*nq*sthetat^4*Tf)/(6*ep^2*Mst1^2) - 
   (CA*CF*Mst2^2*nq*sthetat^4*Tf)/(24*ep*Mst1^2) + (CF^2*Mst2^2*nq*sthetat^4*Tf)/(3*ep*Mst1^2) + 
   (4*CF^2*Mfg^2*Mt^2*nq*sthetat^4*Tf)/(ep^3*Mst1^2*(Mst1^2 - Mst2^2)) - (4*CF^2*Mfg^2*Mt^2*nq*sthetat^4*Tf)/
    (ep^2*Mst1^2*(Mst1^2 - Mst2^2)) - (CA*CF*nt*sthetat^4*Tf)/(4*ep^3) - (3*CF^2*nt*sthetat^4*Tf)/(4*ep^3) + 
   (CA*CF*nt*sthetat^4*Tf)/(3*ep^2) + (5*CF^2*nt*sthetat^4*Tf)/(6*ep^2) + (CA*CF*nt*sthetat^4*Tf)/(24*ep) - (CF^2*nt*sthetat^4*Tf)/(3*ep) + 
   (CA*CF*Mst2^2*nt*sthetat^4*Tf)/(4*ep^3*Mst1^2) + (3*CF^2*Mst2^2*nt*sthetat^4*Tf)/(4*ep^3*Mst1^2) - 
   (CA*CF*Mst2^2*nt*sthetat^4*Tf)/(3*ep^2*Mst1^2) - (5*CF^2*Mst2^2*nt*sthetat^4*Tf)/(6*ep^2*Mst1^2) - 
   (CA*CF*Mst2^2*nt*sthetat^4*Tf)/(24*ep*Mst1^2) + (CF^2*Mst2^2*nt*sthetat^4*Tf)/(3*ep*Mst1^2) + 
   (4*CF^2*Mfg^2*Mt^2*nt*sthetat^4*Tf)/(ep^3*Mst1^2*(Mst1^2 - Mst2^2)) - (4*CF^2*Mfg^2*Mt^2*nt*sthetat^4*Tf)/
    (ep^2*Mst1^2*(Mst1^2 - Mst2^2)) - (3*CF^2*cthetat*Mfg*Mt*nq*sthetat^5*Tf)/(ep^3*Mst1^2) + 
   (3*CF^2*cthetat*Mfg*Mt*nq*sthetat^5*Tf)/(ep^2*Mst1^2) - (3*CF^2*cthetat*Mfg*Mt*nt*sthetat^5*Tf)/(ep^3*Mst1^2) + 
   (3*CF^2*cthetat*Mfg*Mt*nt*sthetat^5*Tf)/(ep^2*Mst1^2) + (CF^2*nq*sthetat^6*Tf)/ep^3 - (CF^2*nq*sthetat^6*Tf)/ep^2 - 
   (CF^2*Mst2^2*nq*sthetat^6*Tf)/(ep^3*Mst1^2) + (CF^2*Mst2^2*nq*sthetat^6*Tf)/(ep^2*Mst1^2) + (CF^2*nt*sthetat^6*Tf)/ep^3 - 
   (CF^2*nt*sthetat^6*Tf)/ep^2 - (CF^2*Mst2^2*nt*sthetat^6*Tf)/(ep^3*Mst1^2) + (CF^2*Mst2^2*nt*sthetat^6*Tf)/(ep^2*Mst1^2) - 
   (CF^2*nq*sthetat^8*Tf)/(2*ep^3) + (CF^2*nq*sthetat^8*Tf)/(2*ep^2) + (CF^2*Mst2^2*nq*sthetat^8*Tf)/(2*ep^3*Mst1^2) - 
   (CF^2*Mst2^2*nq*sthetat^8*Tf)/(2*ep^2*Mst1^2) - (CF^2*nt*sthetat^8*Tf)/(2*ep^3) + (CF^2*nt*sthetat^8*Tf)/(2*ep^2) + 
   (CF^2*Mst2^2*nt*sthetat^8*Tf)/(2*ep^3*Mst1^2) - (CF^2*Mst2^2*nt*sthetat^8*Tf)/(2*ep^2*Mst1^2) + (CF*Me^2*nq^2*Tf^2)/(8*ep^2*Mst1^2) - 
   (CF*Me^2*nq^2*Tf^2)/(8*ep*Mst1^2) - (CF*Mfg^2*nq^2*Tf^2)/(2*ep^3*Mst1^2) + (CF*Mfg^2*nq^2*Tf^2)/(2*ep^2*Mst1^2) + 
   (CF*Mfg^2*nq^2*Tf^2)/(2*ep*Mst1^2) + (CF*Msq^2*nq^2*Tf^2)/(4*ep^2*Mst1^2) - (CF*Msq^2*nq^2*Tf^2)/(4*ep*Mst1^2) - 
   (CF*Mt^2*nq^2*Tf^2)/(12*ep^3*Mst1^2) + (CF*Mt^2*nq^2*Tf^2)/(12*ep^2*Mst1^2) + (CF*Mt^2*nq^2*Tf^2)/(12*ep*Mst1^2) + 
   (CF*nq*nt*Tf^2)/(8*ep^2) - (CF*nq*nt*Tf^2)/(8*ep) + (CF*Me^2*nq*nt*Tf^2)/(4*ep^2*Mst1^2) - (CF*Me^2*nq*nt*Tf^2)/(4*ep*Mst1^2) - 
   (CF*Mfg^2*nq*nt*Tf^2)/(ep^3*Mst1^2) + (CF*Mfg^2*nq*nt*Tf^2)/(ep^2*Mst1^2) + (CF*Mfg^2*nq*nt*Tf^2)/(ep*Mst1^2) + 
   (CF*Msq^2*nq*nt*Tf^2)/(4*ep^2*Mst1^2) - (CF*Msq^2*nq*nt*Tf^2)/(4*ep*Mst1^2) + (CF*Mst2^2*nq*nt*Tf^2)/(8*ep^2*Mst1^2) - 
   (CF*Mst2^2*nq*nt*Tf^2)/(8*ep*Mst1^2) - (CF*Mt^2*nq*nt*Tf^2)/(6*ep^3*Mst1^2) - (CF*Mt^2*nq*nt*Tf^2)/(12*ep^2*Mst1^2) + 
   (5*CF*Mt^2*nq*nt*Tf^2)/(12*ep*Mst1^2) + (CF*nt^2*Tf^2)/(8*ep^2) - (CF*nt^2*Tf^2)/(8*ep) + (CF*Me^2*nt^2*Tf^2)/(8*ep^2*Mst1^2) - 
   (CF*Me^2*nt^2*Tf^2)/(8*ep*Mst1^2) - (CF*Mfg^2*nt^2*Tf^2)/(2*ep^3*Mst1^2) + (CF*Mfg^2*nt^2*Tf^2)/(2*ep^2*Mst1^2) + 
   (CF*Mfg^2*nt^2*Tf^2)/(2*ep*Mst1^2) + (CF*Mst2^2*nt^2*Tf^2)/(8*ep^2*Mst1^2) - (CF*Mst2^2*nt^2*Tf^2)/(8*ep*Mst1^2) - 
   (CF*Mt^2*nt^2*Tf^2)/(12*ep^3*Mst1^2) - (CF*Mt^2*nt^2*Tf^2)/(6*ep^2*Mst1^2) + (CF*Mt^2*nt^2*Tf^2)/(3*ep*Mst1^2) + 
   (CF*cthetat*Mfg*Mt*nq^2*sthetat*Tf^2)/(2*ep^3*Mst1^2) - (CF*cthetat*Mfg*Mt*nq^2*sthetat*Tf^2)/(2*ep^2*Mst1^2) - 
   (CF*cthetat*Mfg*Mt*nq^2*sthetat*Tf^2)/(2*ep*Mst1^2) + (CF*cthetat*Mfg*Mt*nq*nt*sthetat*Tf^2)/(ep^3*Mst1^2) - 
   (CF*cthetat*Mfg*Mt*nq*nt*sthetat*Tf^2)/(ep^2*Mst1^2) - (CF*cthetat*Mfg*Mt*nq*nt*sthetat*Tf^2)/(ep*Mst1^2) + 
   (CF*cthetat*Mfg*Mt*nt^2*sthetat*Tf^2)/(2*ep^3*Mst1^2) - (CF*cthetat*Mfg*Mt*nt^2*sthetat*Tf^2)/(2*ep^2*Mst1^2) - 
   (CF*cthetat*Mfg*Mt*nt^2*sthetat*Tf^2)/(2*ep*Mst1^2) - (CF*nq^2*sthetat^2*Tf^2)/(12*ep^3) + (CF*nq^2*sthetat^2*Tf^2)/(12*ep^2) + 
   (CF*nq^2*sthetat^2*Tf^2)/(12*ep) + (CF*Mst2^2*nq^2*sthetat^2*Tf^2)/(12*ep^3*Mst1^2) - (CF*Mst2^2*nq^2*sthetat^2*Tf^2)/(12*ep^2*Mst1^2) - 
   (CF*Mst2^2*nq^2*sthetat^2*Tf^2)/(12*ep*Mst1^2) - (CF*nq*nt*sthetat^2*Tf^2)/(6*ep^3) + (CF*nq*nt*sthetat^2*Tf^2)/(6*ep^2) + 
   (CF*nq*nt*sthetat^2*Tf^2)/(6*ep) + (CF*Mst2^2*nq*nt*sthetat^2*Tf^2)/(6*ep^3*Mst1^2) - (CF*Mst2^2*nq*nt*sthetat^2*Tf^2)/(6*ep^2*Mst1^2) - 
   (CF*Mst2^2*nq*nt*sthetat^2*Tf^2)/(6*ep*Mst1^2) - (CF*nt^2*sthetat^2*Tf^2)/(12*ep^3) + (CF*nt^2*sthetat^2*Tf^2)/(12*ep^2) + 
   (CF*nt^2*sthetat^2*Tf^2)/(12*ep) + (CF*Mst2^2*nt^2*sthetat^2*Tf^2)/(12*ep^3*Mst1^2) - (CF*Mst2^2*nt^2*sthetat^2*Tf^2)/(12*ep^2*Mst1^2) - 
   (CF*Mst2^2*nt^2*sthetat^2*Tf^2)/(12*ep*Mst1^2) + (CF*nq^2*sthetat^4*Tf^2)/(12*ep^3) - (CF*nq^2*sthetat^4*Tf^2)/(12*ep^2) - 
   (CF*nq^2*sthetat^4*Tf^2)/(12*ep) - (CF*Mst2^2*nq^2*sthetat^4*Tf^2)/(12*ep^3*Mst1^2) + (CF*Mst2^2*nq^2*sthetat^4*Tf^2)/(12*ep^2*Mst1^2) + 
   (CF*Mst2^2*nq^2*sthetat^4*Tf^2)/(12*ep*Mst1^2) + (CF*nq*nt*sthetat^4*Tf^2)/(6*ep^3) - (CF*nq*nt*sthetat^4*Tf^2)/(6*ep^2) - 
   (CF*nq*nt*sthetat^4*Tf^2)/(6*ep) - (CF*Mst2^2*nq*nt*sthetat^4*Tf^2)/(6*ep^3*Mst1^2) + (CF*Mst2^2*nq*nt*sthetat^4*Tf^2)/(6*ep^2*Mst1^2) + 
   (CF*Mst2^2*nq*nt*sthetat^4*Tf^2)/(6*ep*Mst1^2) + (CF*nt^2*sthetat^4*Tf^2)/(12*ep^3) - (CF*nt^2*sthetat^4*Tf^2)/(12*ep^2) - 
   (CF*nt^2*sthetat^4*Tf^2)/(12*ep) - (CF*Mst2^2*nt^2*sthetat^4*Tf^2)/(12*ep^3*Mst1^2) + (CF*Mst2^2*nt^2*sthetat^4*Tf^2)/(12*ep^2*Mst1^2) + 
   (CF*Mst2^2*nt^2*sthetat^4*Tf^2)/(12*ep*Mst1^2) + (3*CA*CF*Mfg^2*nq*Tf*z3)/(ep*Mst1^2) - (3*CF^2*Mfg^2*nq*Tf*z3)/(ep*Mst1^2) + 
   (CA*CF*Mt^2*nq*Tf*z3)/(2*ep*Mst1^2) - (CF^2*Mt^2*nq*Tf*z3)/(2*ep*Mst1^2) + (3*CA*CF*Mfg^2*nt*Tf*z3)/(ep*Mst1^2) - 
   (3*CF^2*Mfg^2*nt*Tf*z3)/(ep*Mst1^2) + (CA*CF*Mt^2*nt*Tf*z3)/(2*ep*Mst1^2) - (CF^2*Mt^2*nt*Tf*z3)/(2*ep*Mst1^2) - 
   (3*CA*CF*cthetat*Mfg*Mt*nq*sthetat*Tf*z3)/(ep*Mst1^2) + (3*CF^2*cthetat*Mfg*Mt*nq*sthetat*Tf*z3)/(ep*Mst1^2) - 
   (3*CA*CF*cthetat*Mfg*Mt*nt*sthetat*Tf*z3)/(ep*Mst1^2) + (3*CF^2*cthetat*Mfg*Mt*nt*sthetat*Tf*z3)/(ep*Mst1^2) + 
   (CA*CF*nq*sthetat^2*Tf*z3)/(2*ep) - (CF^2*nq*sthetat^2*Tf*z3)/(2*ep) - (CA*CF*Mst2^2*nq*sthetat^2*Tf*z3)/(2*ep*Mst1^2) + 
   (CF^2*Mst2^2*nq*sthetat^2*Tf*z3)/(2*ep*Mst1^2) + (CA*CF*nt*sthetat^2*Tf*z3)/(2*ep) - (CF^2*nt*sthetat^2*Tf*z3)/(2*ep) - 
   (CA*CF*Mst2^2*nt*sthetat^2*Tf*z3)/(2*ep*Mst1^2) + (CF^2*Mst2^2*nt*sthetat^2*Tf*z3)/(2*ep*Mst1^2) - (CA*CF*nq*sthetat^4*Tf*z3)/(2*ep) + 
   (CF^2*nq*sthetat^4*Tf*z3)/(2*ep) + (CA*CF*Mst2^2*nq*sthetat^4*Tf*z3)/(2*ep*Mst1^2) - (CF^2*Mst2^2*nq*sthetat^4*Tf*z3)/(2*ep*Mst1^2) - 
   (CA*CF*nt*sthetat^4*Tf*z3)/(2*ep) + (CF^2*nt*sthetat^4*Tf*z3)/(2*ep) + (CA*CF*Mst2^2*nt*sthetat^4*Tf*z3)/(2*ep*Mst1^2) - 
   (CF^2*Mst2^2*nt*sthetat^4*Tf*z3)/(2*ep*Mst1^2)));


Zmst2 = (1 + (as*(-(CF*Mfg^2) - CF*Mt^2 - 2*CF*cthetat*Mfg*Mt*sthetat + CF*Mst1^2*sthetat^2 - CF*Mst2^2*sthetat^2 - CF*Mst1^2*sthetat^4 + 
    CF*Mst2^2*sthetat^4))/(ep*Mst2^2) + as^2*((-3*CA*CF*Me^2)/(8*ep*Mst2^2) + (9*CA*CF*Mfg^2)/(8*ep^2*Mst2^2) - 
   (11*CA*CF*Mfg^2)/(8*ep*Mst2^2) + (3*CF^2*Mfg^2)/(4*ep*Mst2^2) + (3*CA*CF*Mt^2)/(8*ep^2*Mst2^2) + (CF^2*Mt^2)/(2*ep^2*Mst2^2) - 
   (3*CA*CF*Mt^2)/(8*ep*Mst2^2) + (CF^2*Mt^2)/(4*ep*Mst2^2) + (CF^2*Mfg^2*Mt^2)/(ep^2*Mst2^2*(-Mst1^2 + Mst2^2)) + 
   (3*CA*CF*cthetat*Mfg*Mt*sthetat)/(2*ep^2*Mst2^2) + (2*CF^2*cthetat*Mfg*Mt*sthetat)/(ep^2*Mst2^2) - 
   (3*CA*CF*cthetat*Mfg*Mt*sthetat)/(2*ep*Mst2^2) + (CF^2*cthetat*Mfg*Mt*sthetat)/(ep*Mst2^2) + (3*CA*CF*sthetat^2)/(8*ep^2) + 
   (CF^2*sthetat^2)/(2*ep^2) - (3*CA*CF*sthetat^2)/(8*ep) + (CF^2*sthetat^2)/(4*ep) - (3*CA*CF*Mst1^2*sthetat^2)/(8*ep^2*Mst2^2) - 
   (CF^2*Mst1^2*sthetat^2)/(2*ep^2*Mst2^2) + (3*CA*CF*Mst1^2*sthetat^2)/(8*ep*Mst2^2) - (CF^2*Mst1^2*sthetat^2)/(4*ep*Mst2^2) - 
   (4*CF^2*Mfg^2*Mt^2*sthetat^2)/(ep^2*Mst2^2*(-Mst1^2 + Mst2^2)) - (4*CF^2*cthetat*Mfg*Mt*sthetat^3)/(ep^2*Mst2^2) - 
   (3*CA*CF*sthetat^4)/(8*ep^2) - (3*CF^2*sthetat^4)/(2*ep^2) + (3*CA*CF*sthetat^4)/(8*ep) - (CF^2*sthetat^4)/(4*ep) + 
   (3*CA*CF*Mst1^2*sthetat^4)/(8*ep^2*Mst2^2) + (3*CF^2*Mst1^2*sthetat^4)/(2*ep^2*Mst2^2) - (3*CA*CF*Mst1^2*sthetat^4)/(8*ep*Mst2^2) + 
   (CF^2*Mst1^2*sthetat^4)/(4*ep*Mst2^2) + (4*CF^2*Mfg^2*Mt^2*sthetat^4)/(ep^2*Mst2^2*(-Mst1^2 + Mst2^2)) + 
   (4*CF^2*cthetat*Mfg*Mt*sthetat^5)/(ep^2*Mst2^2) + (2*CF^2*sthetat^6)/ep^2 - (2*CF^2*Mst1^2*sthetat^6)/(ep^2*Mst2^2) - 
   (CF^2*sthetat^8)/ep^2 + (CF^2*Mst1^2*sthetat^8)/(ep^2*Mst2^2) + (CF*Me^2*nq*Tf)/(4*ep*Mst2^2) - (3*CF*Mfg^2*nq*Tf)/(4*ep^2*Mst2^2) + 
   (3*CF*Mfg^2*nq*Tf)/(4*ep*Mst2^2) + (CF*Msq^2*nq*Tf)/(2*ep*Mst2^2) - (CF*Mt^2*nq*Tf)/(4*ep^2*Mst2^2) + (CF*Mt^2*nq*Tf)/(4*ep*Mst2^2) + 
   (CF*nt*Tf)/(4*ep) + (CF*Me^2*nt*Tf)/(4*ep*Mst2^2) - (3*CF*Mfg^2*nt*Tf)/(4*ep^2*Mst2^2) + (3*CF*Mfg^2*nt*Tf)/(4*ep*Mst2^2) + 
   (CF*Mst1^2*nt*Tf)/(4*ep*Mst2^2) - (CF*Mt^2*nt*Tf)/(4*ep^2*Mst2^2) - (CF*Mt^2*nt*Tf)/(4*ep*Mst2^2) - 
   (CF*cthetat*Mfg*Mt*nq*sthetat*Tf)/(ep^2*Mst2^2) + (CF*cthetat*Mfg*Mt*nq*sthetat*Tf)/(ep*Mst2^2) - 
   (CF*cthetat*Mfg*Mt*nt*sthetat*Tf)/(ep^2*Mst2^2) + (CF*cthetat*Mfg*Mt*nt*sthetat*Tf)/(ep*Mst2^2) - (CF*nq*sthetat^2*Tf)/(4*ep^2) + 
   (CF*nq*sthetat^2*Tf)/(4*ep) + (CF*Mst1^2*nq*sthetat^2*Tf)/(4*ep^2*Mst2^2) - (CF*Mst1^2*nq*sthetat^2*Tf)/(4*ep*Mst2^2) - 
   (CF*nt*sthetat^2*Tf)/(4*ep^2) + (CF*nt*sthetat^2*Tf)/(4*ep) + (CF*Mst1^2*nt*sthetat^2*Tf)/(4*ep^2*Mst2^2) - 
   (CF*Mst1^2*nt*sthetat^2*Tf)/(4*ep*Mst2^2) + (CF*nq*sthetat^4*Tf)/(4*ep^2) - (CF*nq*sthetat^4*Tf)/(4*ep) - 
   (CF*Mst1^2*nq*sthetat^4*Tf)/(4*ep^2*Mst2^2) + (CF*Mst1^2*nq*sthetat^4*Tf)/(4*ep*Mst2^2) + (CF*nt*sthetat^4*Tf)/(4*ep^2) - 
   (CF*nt*sthetat^4*Tf)/(4*ep) - (CF*Mst1^2*nt*sthetat^4*Tf)/(4*ep^2*Mst2^2) + (CF*Mst1^2*nt*sthetat^4*Tf)/(4*ep*Mst2^2)) + 
 as^3*((9*CA^2*CF*Me^2)/(32*ep^2*Mst2^2) - (15*CA^2*CF*Me^2)/(32*ep*Mst2^2) + (3*CA*CF^2*Me^2)/(16*ep*Mst2^2) - 
   (9*CA^2*CF*Mfg^2)/(8*ep^3*Mst2^2) + (33*CA^2*CF*Mfg^2)/(16*ep^2*Mst2^2) - (3*CA*CF^2*Mfg^2)/(4*ep^2*Mst2^2) - 
   (5*CA^2*CF*Mfg^2)/(4*ep*Mst2^2) + (7*CA*CF^2*Mfg^2)/(8*ep*Mst2^2) - (CF^3*Mfg^2)/(ep*Mst2^2) - (3*CA^2*CF*Mt^2)/(16*ep^3*Mst2^2) - 
   (3*CA*CF^2*Mt^2)/(8*ep^3*Mst2^2) - (CF^3*Mt^2)/(6*ep^3*Mst2^2) + (5*CA^2*CF*Mt^2)/(16*ep^2*Mst2^2) + (CA*CF^2*Mt^2)/(4*ep^2*Mst2^2) - 
   (CF^3*Mt^2)/(4*ep^2*Mst2^2) - (CA^2*CF*Mt^2)/(8*ep*Mst2^2) + (CA*CF^2*Mt^2)/(8*ep*Mst2^2) - (CF^3*Mt^2)/(6*ep*Mst2^2) - 
   (3*CA*CF^2*Mfg^2*Mt^2)/(2*ep^3*Mst2^2*(-Mst1^2 + Mst2^2)) - (CF^3*Mfg^2*Mt^2)/(ep^3*Mst2^2*(-Mst1^2 + Mst2^2)) + 
   (3*CA*CF^2*Mfg^2*Mt^2)/(2*ep^2*Mst2^2*(-Mst1^2 + Mst2^2)) - (CF^3*Mfg^2*Mt^2)/(ep^2*Mst2^2*(-Mst1^2 + Mst2^2)) - 
   (9*CA^2*CF*cthetat*Mfg*Mt*sthetat)/(8*ep^3*Mst2^2) - (9*CA*CF^2*cthetat*Mfg*Mt*sthetat)/(4*ep^3*Mst2^2) - 
   (CF^3*cthetat*Mfg*Mt*sthetat)/(ep^3*Mst2^2) + (15*CA^2*CF*cthetat*Mfg*Mt*sthetat)/(8*ep^2*Mst2^2) + 
   (3*CA*CF^2*cthetat*Mfg*Mt*sthetat)/(2*ep^2*Mst2^2) - (3*CF^3*cthetat*Mfg*Mt*sthetat)/(2*ep^2*Mst2^2) - 
   (3*CA^2*CF*cthetat*Mfg*Mt*sthetat)/(4*ep*Mst2^2) + (3*CA*CF^2*cthetat*Mfg*Mt*sthetat)/(4*ep*Mst2^2) - 
   (CF^3*cthetat*Mfg*Mt*sthetat)/(ep*Mst2^2) + (4*CF^3*cthetat*Mfg^3*Mt^3*sthetat)/(ep^3*Mst2^2*(-Mst1^2 + Mst2^2)^2) - 
   (3*CA^2*CF*sthetat^2)/(16*ep^3) - (3*CA*CF^2*sthetat^2)/(8*ep^3) - (CF^3*sthetat^2)/(6*ep^3) + (5*CA^2*CF*sthetat^2)/(16*ep^2) + 
   (CA*CF^2*sthetat^2)/(4*ep^2) - (CF^3*sthetat^2)/(4*ep^2) - (CA^2*CF*sthetat^2)/(8*ep) + (CA*CF^2*sthetat^2)/(8*ep) - 
   (CF^3*sthetat^2)/(6*ep) + (3*CA^2*CF*Mst1^2*sthetat^2)/(16*ep^3*Mst2^2) + (3*CA*CF^2*Mst1^2*sthetat^2)/(8*ep^3*Mst2^2) + 
   (CF^3*Mst1^2*sthetat^2)/(6*ep^3*Mst2^2) - (5*CA^2*CF*Mst1^2*sthetat^2)/(16*ep^2*Mst2^2) - (CA*CF^2*Mst1^2*sthetat^2)/(4*ep^2*Mst2^2) + 
   (CF^3*Mst1^2*sthetat^2)/(4*ep^2*Mst2^2) + (CA^2*CF*Mst1^2*sthetat^2)/(8*ep*Mst2^2) - (CA*CF^2*Mst1^2*sthetat^2)/(8*ep*Mst2^2) + 
   (CF^3*Mst1^2*sthetat^2)/(6*ep*Mst2^2) + (6*CA*CF^2*Mfg^2*Mt^2*sthetat^2)/(ep^3*Mst2^2*(-Mst1^2 + Mst2^2)) + 
   (10*CF^3*Mfg^2*Mt^2*sthetat^2)/(ep^3*Mst2^2*(-Mst1^2 + Mst2^2)) - (6*CA*CF^2*Mfg^2*Mt^2*sthetat^2)/(ep^2*Mst2^2*(-Mst1^2 + Mst2^2)) + 
   (4*CF^3*Mfg^2*Mt^2*sthetat^2)/(ep^2*Mst2^2*(-Mst1^2 + Mst2^2)) + (9*CA*CF^2*cthetat*Mfg*Mt*sthetat^3)/(2*ep^3*Mst2^2) + 
   (6*CF^3*cthetat*Mfg*Mt*sthetat^3)/(ep^3*Mst2^2) - (9*CA*CF^2*cthetat*Mfg*Mt*sthetat^3)/(2*ep^2*Mst2^2) + 
   (3*CF^3*cthetat*Mfg*Mt*sthetat^3)/(ep^2*Mst2^2) - (16*CF^3*cthetat*Mfg^3*Mt^3*sthetat^3)/(ep^3*Mst2^2*(-Mst1^2 + Mst2^2)^2) + 
   (3*CA^2*CF*sthetat^4)/(16*ep^3) + (9*CA*CF^2*sthetat^4)/(8*ep^3) + (7*CF^3*sthetat^4)/(6*ep^3) - (5*CA^2*CF*sthetat^4)/(16*ep^2) - 
   (CA*CF^2*sthetat^4)/ep^2 + (3*CF^3*sthetat^4)/(4*ep^2) + (CA^2*CF*sthetat^4)/(8*ep) - (CA*CF^2*sthetat^4)/(8*ep) + 
   (CF^3*sthetat^4)/(6*ep) - (3*CA^2*CF*Mst1^2*sthetat^4)/(16*ep^3*Mst2^2) - (9*CA*CF^2*Mst1^2*sthetat^4)/(8*ep^3*Mst2^2) - 
   (7*CF^3*Mst1^2*sthetat^4)/(6*ep^3*Mst2^2) + (5*CA^2*CF*Mst1^2*sthetat^4)/(16*ep^2*Mst2^2) + (CA*CF^2*Mst1^2*sthetat^4)/(ep^2*Mst2^2) - 
   (3*CF^3*Mst1^2*sthetat^4)/(4*ep^2*Mst2^2) - (CA^2*CF*Mst1^2*sthetat^4)/(8*ep*Mst2^2) + (CA*CF^2*Mst1^2*sthetat^4)/(8*ep*Mst2^2) - 
   (CF^3*Mst1^2*sthetat^4)/(6*ep*Mst2^2) - (6*CA*CF^2*Mfg^2*Mt^2*sthetat^4)/(ep^3*Mst2^2*(-Mst1^2 + Mst2^2)) - 
   (34*CF^3*Mfg^2*Mt^2*sthetat^4)/(ep^3*Mst2^2*(-Mst1^2 + Mst2^2)) + (6*CA*CF^2*Mfg^2*Mt^2*sthetat^4)/(ep^2*Mst2^2*(-Mst1^2 + Mst2^2)) - 
   (4*CF^3*Mfg^2*Mt^2*sthetat^4)/(ep^2*Mst2^2*(-Mst1^2 + Mst2^2)) - (9*CA*CF^2*cthetat*Mfg*Mt*sthetat^5)/(2*ep^3*Mst2^2) - 
   (18*CF^3*cthetat*Mfg*Mt*sthetat^5)/(ep^3*Mst2^2) + (9*CA*CF^2*cthetat*Mfg*Mt*sthetat^5)/(2*ep^2*Mst2^2) - 
   (3*CF^3*cthetat*Mfg*Mt*sthetat^5)/(ep^2*Mst2^2) + (16*CF^3*cthetat*Mfg^3*Mt^3*sthetat^5)/(ep^3*Mst2^2*(-Mst1^2 + Mst2^2)^2) - 
   (3*CA*CF^2*sthetat^6)/(2*ep^3) - (4*CF^3*sthetat^6)/ep^3 + (3*CA*CF^2*sthetat^6)/(2*ep^2) - (CF^3*sthetat^6)/ep^2 + 
   (3*CA*CF^2*Mst1^2*sthetat^6)/(2*ep^3*Mst2^2) + (4*CF^3*Mst1^2*sthetat^6)/(ep^3*Mst2^2) - (3*CA*CF^2*Mst1^2*sthetat^6)/(2*ep^2*Mst2^2) + 
   (CF^3*Mst1^2*sthetat^6)/(ep^2*Mst2^2) + (48*CF^3*Mfg^2*Mt^2*sthetat^6)/(ep^3*Mst2^2*(-Mst1^2 + Mst2^2)) + 
   (24*CF^3*cthetat*Mfg*Mt*sthetat^7)/(ep^3*Mst2^2) + (3*CA*CF^2*sthetat^8)/(4*ep^3) + (7*CF^3*sthetat^8)/ep^3 - 
   (3*CA*CF^2*sthetat^8)/(4*ep^2) + (CF^3*sthetat^8)/(2*ep^2) - (3*CA*CF^2*Mst1^2*sthetat^8)/(4*ep^3*Mst2^2) - 
   (7*CF^3*Mst1^2*sthetat^8)/(ep^3*Mst2^2) + (3*CA*CF^2*Mst1^2*sthetat^8)/(4*ep^2*Mst2^2) - (CF^3*Mst1^2*sthetat^8)/(2*ep^2*Mst2^2) - 
   (24*CF^3*Mfg^2*Mt^2*sthetat^8)/(ep^3*Mst2^2*(-Mst1^2 + Mst2^2)) - (12*CF^3*cthetat*Mfg*Mt*sthetat^9)/(ep^3*Mst2^2) - 
   (6*CF^3*sthetat^10)/ep^3 + (6*CF^3*Mst1^2*sthetat^10)/(ep^3*Mst2^2) + (2*CF^3*sthetat^12)/ep^3 - 
   (2*CF^3*Mst1^2*sthetat^12)/(ep^3*Mst2^2) - (3*CA*CF*Me^2*nq*Tf)/(8*ep^2*Mst2^2) + (CA*CF*Me^2*nq*Tf)/(2*ep*Mst2^2) - 
   (CF^2*Me^2*nq*Tf)/(8*ep*Mst2^2) + (3*CA*CF*Mfg^2*nq*Tf)/(2*ep^3*Mst2^2) - (17*CA*CF*Mfg^2*nq*Tf)/(8*ep^2*Mst2^2) - 
   (CF^2*Mfg^2*nq*Tf)/(2*ep^2*Mst2^2) + (9*CF^2*Mfg^2*nq*Tf)/(4*ep*Mst2^2) - (3*CA*CF*Msq^2*nq*Tf)/(8*ep^2*Mst2^2) + 
   (5*CA*CF*Msq^2*nq*Tf)/(8*ep*Mst2^2) - (CF^2*Msq^2*nq*Tf)/(4*ep*Mst2^2) + (CA*CF*Mt^2*nq*Tf)/(4*ep^3*Mst2^2) + 
   (CF^2*Mt^2*nq*Tf)/(4*ep^3*Mst2^2) - (CA*CF*Mt^2*nq*Tf)/(3*ep^2*Mst2^2) - (CF^2*Mt^2*nq*Tf)/(3*ep^2*Mst2^2) - 
   (CA*CF*Mt^2*nq*Tf)/(24*ep*Mst2^2) + (CF^2*Mt^2*nq*Tf)/(3*ep*Mst2^2) + (CF^2*Mfg^2*Mt^2*nq*Tf)/(ep^3*Mst2^2*(-Mst1^2 + Mst2^2)) - 
   (CF^2*Mfg^2*Mt^2*nq*Tf)/(ep^2*Mst2^2*(-Mst1^2 + Mst2^2)) - (3*CA*CF*nt*Tf)/(16*ep^2) + (5*CA*CF*nt*Tf)/(16*ep) - (CF^2*nt*Tf)/(8*ep) - 
   (3*CA*CF*Me^2*nt*Tf)/(8*ep^2*Mst2^2) + (CA*CF*Me^2*nt*Tf)/(2*ep*Mst2^2) - (CF^2*Me^2*nt*Tf)/(8*ep*Mst2^2) + 
   (3*CA*CF*Mfg^2*nt*Tf)/(2*ep^3*Mst2^2) - (17*CA*CF*Mfg^2*nt*Tf)/(8*ep^2*Mst2^2) - (CF^2*Mfg^2*nt*Tf)/(2*ep^2*Mst2^2) + 
   (9*CF^2*Mfg^2*nt*Tf)/(4*ep*Mst2^2) - (3*CA*CF*Mst1^2*nt*Tf)/(16*ep^2*Mst2^2) + (5*CA*CF*Mst1^2*nt*Tf)/(16*ep*Mst2^2) - 
   (CF^2*Mst1^2*nt*Tf)/(8*ep*Mst2^2) + (CA*CF*Mt^2*nt*Tf)/(4*ep^3*Mst2^2) + (CF^2*Mt^2*nt*Tf)/(4*ep^3*Mst2^2) + 
   (CA*CF*Mt^2*nt*Tf)/(24*ep^2*Mst2^2) - (CF^2*Mt^2*nt*Tf)/(3*ep^2*Mst2^2) - (2*CA*CF*Mt^2*nt*Tf)/(3*ep*Mst2^2) + 
   (7*CF^2*Mt^2*nt*Tf)/(12*ep*Mst2^2) + (CF^2*Mfg^2*Mt^2*nt*Tf)/(ep^3*Mst2^2*(-Mst1^2 + Mst2^2)) - 
   (CF^2*Mfg^2*Mt^2*nt*Tf)/(ep^2*Mst2^2*(-Mst1^2 + Mst2^2)) + (3*CA*CF*cthetat*Mfg*Mt*nq*sthetat*Tf)/(2*ep^3*Mst2^2) + 
   (3*CF^2*cthetat*Mfg*Mt*nq*sthetat*Tf)/(2*ep^3*Mst2^2) - (2*CA*CF*cthetat*Mfg*Mt*nq*sthetat*Tf)/(ep^2*Mst2^2) - 
   (2*CF^2*cthetat*Mfg*Mt*nq*sthetat*Tf)/(ep^2*Mst2^2) - (CA*CF*cthetat*Mfg*Mt*nq*sthetat*Tf)/(4*ep*Mst2^2) + 
   (2*CF^2*cthetat*Mfg*Mt*nq*sthetat*Tf)/(ep*Mst2^2) + (3*CA*CF*cthetat*Mfg*Mt*nt*sthetat*Tf)/(2*ep^3*Mst2^2) + 
   (3*CF^2*cthetat*Mfg*Mt*nt*sthetat*Tf)/(2*ep^3*Mst2^2) - (2*CA*CF*cthetat*Mfg*Mt*nt*sthetat*Tf)/(ep^2*Mst2^2) - 
   (2*CF^2*cthetat*Mfg*Mt*nt*sthetat*Tf)/(ep^2*Mst2^2) - (CA*CF*cthetat*Mfg*Mt*nt*sthetat*Tf)/(4*ep*Mst2^2) + 
   (2*CF^2*cthetat*Mfg*Mt*nt*sthetat*Tf)/(ep*Mst2^2) + (CA*CF*nq*sthetat^2*Tf)/(4*ep^3) + (CF^2*nq*sthetat^2*Tf)/(4*ep^3) - 
   (CA*CF*nq*sthetat^2*Tf)/(3*ep^2) - (CF^2*nq*sthetat^2*Tf)/(3*ep^2) - (CA*CF*nq*sthetat^2*Tf)/(24*ep) + (CF^2*nq*sthetat^2*Tf)/(3*ep) - 
   (CA*CF*Mst1^2*nq*sthetat^2*Tf)/(4*ep^3*Mst2^2) - (CF^2*Mst1^2*nq*sthetat^2*Tf)/(4*ep^3*Mst2^2) + 
   (CA*CF*Mst1^2*nq*sthetat^2*Tf)/(3*ep^2*Mst2^2) + (CF^2*Mst1^2*nq*sthetat^2*Tf)/(3*ep^2*Mst2^2) + 
   (CA*CF*Mst1^2*nq*sthetat^2*Tf)/(24*ep*Mst2^2) - (CF^2*Mst1^2*nq*sthetat^2*Tf)/(3*ep*Mst2^2) - 
   (4*CF^2*Mfg^2*Mt^2*nq*sthetat^2*Tf)/(ep^3*Mst2^2*(-Mst1^2 + Mst2^2)) + (4*CF^2*Mfg^2*Mt^2*nq*sthetat^2*Tf)/
    (ep^2*Mst2^2*(-Mst1^2 + Mst2^2)) + (CA*CF*nt*sthetat^2*Tf)/(4*ep^3) + (CF^2*nt*sthetat^2*Tf)/(4*ep^3) - 
   (CA*CF*nt*sthetat^2*Tf)/(3*ep^2) - (CF^2*nt*sthetat^2*Tf)/(3*ep^2) - (CA*CF*nt*sthetat^2*Tf)/(24*ep) + (CF^2*nt*sthetat^2*Tf)/(3*ep) - 
   (CA*CF*Mst1^2*nt*sthetat^2*Tf)/(4*ep^3*Mst2^2) - (CF^2*Mst1^2*nt*sthetat^2*Tf)/(4*ep^3*Mst2^2) + 
   (CA*CF*Mst1^2*nt*sthetat^2*Tf)/(3*ep^2*Mst2^2) + (CF^2*Mst1^2*nt*sthetat^2*Tf)/(3*ep^2*Mst2^2) + 
   (CA*CF*Mst1^2*nt*sthetat^2*Tf)/(24*ep*Mst2^2) - (CF^2*Mst1^2*nt*sthetat^2*Tf)/(3*ep*Mst2^2) - 
   (4*CF^2*Mfg^2*Mt^2*nt*sthetat^2*Tf)/(ep^3*Mst2^2*(-Mst1^2 + Mst2^2)) + (4*CF^2*Mfg^2*Mt^2*nt*sthetat^2*Tf)/
    (ep^2*Mst2^2*(-Mst1^2 + Mst2^2)) - (3*CF^2*cthetat*Mfg*Mt*nq*sthetat^3*Tf)/(ep^3*Mst2^2) + 
   (3*CF^2*cthetat*Mfg*Mt*nq*sthetat^3*Tf)/(ep^2*Mst2^2) - (3*CF^2*cthetat*Mfg*Mt*nt*sthetat^3*Tf)/(ep^3*Mst2^2) + 
   (3*CF^2*cthetat*Mfg*Mt*nt*sthetat^3*Tf)/(ep^2*Mst2^2) - (CA*CF*nq*sthetat^4*Tf)/(4*ep^3) - (3*CF^2*nq*sthetat^4*Tf)/(4*ep^3) + 
   (CA*CF*nq*sthetat^4*Tf)/(3*ep^2) + (5*CF^2*nq*sthetat^4*Tf)/(6*ep^2) + (CA*CF*nq*sthetat^4*Tf)/(24*ep) - (CF^2*nq*sthetat^4*Tf)/(3*ep) + 
   (CA*CF*Mst1^2*nq*sthetat^4*Tf)/(4*ep^3*Mst2^2) + (3*CF^2*Mst1^2*nq*sthetat^4*Tf)/(4*ep^3*Mst2^2) - 
   (CA*CF*Mst1^2*nq*sthetat^4*Tf)/(3*ep^2*Mst2^2) - (5*CF^2*Mst1^2*nq*sthetat^4*Tf)/(6*ep^2*Mst2^2) - 
   (CA*CF*Mst1^2*nq*sthetat^4*Tf)/(24*ep*Mst2^2) + (CF^2*Mst1^2*nq*sthetat^4*Tf)/(3*ep*Mst2^2) + 
   (4*CF^2*Mfg^2*Mt^2*nq*sthetat^4*Tf)/(ep^3*Mst2^2*(-Mst1^2 + Mst2^2)) - (4*CF^2*Mfg^2*Mt^2*nq*sthetat^4*Tf)/
    (ep^2*Mst2^2*(-Mst1^2 + Mst2^2)) - (CA*CF*nt*sthetat^4*Tf)/(4*ep^3) - (3*CF^2*nt*sthetat^4*Tf)/(4*ep^3) + 
   (CA*CF*nt*sthetat^4*Tf)/(3*ep^2) + (5*CF^2*nt*sthetat^4*Tf)/(6*ep^2) + (CA*CF*nt*sthetat^4*Tf)/(24*ep) - (CF^2*nt*sthetat^4*Tf)/(3*ep) + 
   (CA*CF*Mst1^2*nt*sthetat^4*Tf)/(4*ep^3*Mst2^2) + (3*CF^2*Mst1^2*nt*sthetat^4*Tf)/(4*ep^3*Mst2^2) - 
   (CA*CF*Mst1^2*nt*sthetat^4*Tf)/(3*ep^2*Mst2^2) - (5*CF^2*Mst1^2*nt*sthetat^4*Tf)/(6*ep^2*Mst2^2) - 
   (CA*CF*Mst1^2*nt*sthetat^4*Tf)/(24*ep*Mst2^2) + (CF^2*Mst1^2*nt*sthetat^4*Tf)/(3*ep*Mst2^2) + 
   (4*CF^2*Mfg^2*Mt^2*nt*sthetat^4*Tf)/(ep^3*Mst2^2*(-Mst1^2 + Mst2^2)) - (4*CF^2*Mfg^2*Mt^2*nt*sthetat^4*Tf)/
    (ep^2*Mst2^2*(-Mst1^2 + Mst2^2)) + (3*CF^2*cthetat*Mfg*Mt*nq*sthetat^5*Tf)/(ep^3*Mst2^2) - 
   (3*CF^2*cthetat*Mfg*Mt*nq*sthetat^5*Tf)/(ep^2*Mst2^2) + (3*CF^2*cthetat*Mfg*Mt*nt*sthetat^5*Tf)/(ep^3*Mst2^2) - 
   (3*CF^2*cthetat*Mfg*Mt*nt*sthetat^5*Tf)/(ep^2*Mst2^2) + (CF^2*nq*sthetat^6*Tf)/ep^3 - (CF^2*nq*sthetat^6*Tf)/ep^2 - 
   (CF^2*Mst1^2*nq*sthetat^6*Tf)/(ep^3*Mst2^2) + (CF^2*Mst1^2*nq*sthetat^6*Tf)/(ep^2*Mst2^2) + (CF^2*nt*sthetat^6*Tf)/ep^3 - 
   (CF^2*nt*sthetat^6*Tf)/ep^2 - (CF^2*Mst1^2*nt*sthetat^6*Tf)/(ep^3*Mst2^2) + (CF^2*Mst1^2*nt*sthetat^6*Tf)/(ep^2*Mst2^2) - 
   (CF^2*nq*sthetat^8*Tf)/(2*ep^3) + (CF^2*nq*sthetat^8*Tf)/(2*ep^2) + (CF^2*Mst1^2*nq*sthetat^8*Tf)/(2*ep^3*Mst2^2) - 
   (CF^2*Mst1^2*nq*sthetat^8*Tf)/(2*ep^2*Mst2^2) - (CF^2*nt*sthetat^8*Tf)/(2*ep^3) + (CF^2*nt*sthetat^8*Tf)/(2*ep^2) + 
   (CF^2*Mst1^2*nt*sthetat^8*Tf)/(2*ep^3*Mst2^2) - (CF^2*Mst1^2*nt*sthetat^8*Tf)/(2*ep^2*Mst2^2) + (CF*Me^2*nq^2*Tf^2)/(8*ep^2*Mst2^2) - 
   (CF*Me^2*nq^2*Tf^2)/(8*ep*Mst2^2) - (CF*Mfg^2*nq^2*Tf^2)/(2*ep^3*Mst2^2) + (CF*Mfg^2*nq^2*Tf^2)/(2*ep^2*Mst2^2) + 
   (CF*Mfg^2*nq^2*Tf^2)/(2*ep*Mst2^2) + (CF*Msq^2*nq^2*Tf^2)/(4*ep^2*Mst2^2) - (CF*Msq^2*nq^2*Tf^2)/(4*ep*Mst2^2) - 
   (CF*Mt^2*nq^2*Tf^2)/(12*ep^3*Mst2^2) + (CF*Mt^2*nq^2*Tf^2)/(12*ep^2*Mst2^2) + (CF*Mt^2*nq^2*Tf^2)/(12*ep*Mst2^2) + 
   (CF*nq*nt*Tf^2)/(8*ep^2) - (CF*nq*nt*Tf^2)/(8*ep) + (CF*Me^2*nq*nt*Tf^2)/(4*ep^2*Mst2^2) - (CF*Me^2*nq*nt*Tf^2)/(4*ep*Mst2^2) - 
   (CF*Mfg^2*nq*nt*Tf^2)/(ep^3*Mst2^2) + (CF*Mfg^2*nq*nt*Tf^2)/(ep^2*Mst2^2) + (CF*Mfg^2*nq*nt*Tf^2)/(ep*Mst2^2) + 
   (CF*Msq^2*nq*nt*Tf^2)/(4*ep^2*Mst2^2) - (CF*Msq^2*nq*nt*Tf^2)/(4*ep*Mst2^2) + (CF*Mst1^2*nq*nt*Tf^2)/(8*ep^2*Mst2^2) - 
   (CF*Mst1^2*nq*nt*Tf^2)/(8*ep*Mst2^2) - (CF*Mt^2*nq*nt*Tf^2)/(6*ep^3*Mst2^2) - (CF*Mt^2*nq*nt*Tf^2)/(12*ep^2*Mst2^2) + 
   (5*CF*Mt^2*nq*nt*Tf^2)/(12*ep*Mst2^2) + (CF*nt^2*Tf^2)/(8*ep^2) - (CF*nt^2*Tf^2)/(8*ep) + (CF*Me^2*nt^2*Tf^2)/(8*ep^2*Mst2^2) - 
   (CF*Me^2*nt^2*Tf^2)/(8*ep*Mst2^2) - (CF*Mfg^2*nt^2*Tf^2)/(2*ep^3*Mst2^2) + (CF*Mfg^2*nt^2*Tf^2)/(2*ep^2*Mst2^2) + 
   (CF*Mfg^2*nt^2*Tf^2)/(2*ep*Mst2^2) + (CF*Mst1^2*nt^2*Tf^2)/(8*ep^2*Mst2^2) - (CF*Mst1^2*nt^2*Tf^2)/(8*ep*Mst2^2) - 
   (CF*Mt^2*nt^2*Tf^2)/(12*ep^3*Mst2^2) - (CF*Mt^2*nt^2*Tf^2)/(6*ep^2*Mst2^2) + (CF*Mt^2*nt^2*Tf^2)/(3*ep*Mst2^2) - 
   (CF*cthetat*Mfg*Mt*nq^2*sthetat*Tf^2)/(2*ep^3*Mst2^2) + (CF*cthetat*Mfg*Mt*nq^2*sthetat*Tf^2)/(2*ep^2*Mst2^2) + 
   (CF*cthetat*Mfg*Mt*nq^2*sthetat*Tf^2)/(2*ep*Mst2^2) - (CF*cthetat*Mfg*Mt*nq*nt*sthetat*Tf^2)/(ep^3*Mst2^2) + 
   (CF*cthetat*Mfg*Mt*nq*nt*sthetat*Tf^2)/(ep^2*Mst2^2) + (CF*cthetat*Mfg*Mt*nq*nt*sthetat*Tf^2)/(ep*Mst2^2) - 
   (CF*cthetat*Mfg*Mt*nt^2*sthetat*Tf^2)/(2*ep^3*Mst2^2) + (CF*cthetat*Mfg*Mt*nt^2*sthetat*Tf^2)/(2*ep^2*Mst2^2) + 
   (CF*cthetat*Mfg*Mt*nt^2*sthetat*Tf^2)/(2*ep*Mst2^2) - (CF*nq^2*sthetat^2*Tf^2)/(12*ep^3) + (CF*nq^2*sthetat^2*Tf^2)/(12*ep^2) + 
   (CF*nq^2*sthetat^2*Tf^2)/(12*ep) + (CF*Mst1^2*nq^2*sthetat^2*Tf^2)/(12*ep^3*Mst2^2) - (CF*Mst1^2*nq^2*sthetat^2*Tf^2)/(12*ep^2*Mst2^2) - 
   (CF*Mst1^2*nq^2*sthetat^2*Tf^2)/(12*ep*Mst2^2) - (CF*nq*nt*sthetat^2*Tf^2)/(6*ep^3) + (CF*nq*nt*sthetat^2*Tf^2)/(6*ep^2) + 
   (CF*nq*nt*sthetat^2*Tf^2)/(6*ep) + (CF*Mst1^2*nq*nt*sthetat^2*Tf^2)/(6*ep^3*Mst2^2) - (CF*Mst1^2*nq*nt*sthetat^2*Tf^2)/(6*ep^2*Mst2^2) - 
   (CF*Mst1^2*nq*nt*sthetat^2*Tf^2)/(6*ep*Mst2^2) - (CF*nt^2*sthetat^2*Tf^2)/(12*ep^3) + (CF*nt^2*sthetat^2*Tf^2)/(12*ep^2) + 
   (CF*nt^2*sthetat^2*Tf^2)/(12*ep) + (CF*Mst1^2*nt^2*sthetat^2*Tf^2)/(12*ep^3*Mst2^2) - (CF*Mst1^2*nt^2*sthetat^2*Tf^2)/(12*ep^2*Mst2^2) - 
   (CF*Mst1^2*nt^2*sthetat^2*Tf^2)/(12*ep*Mst2^2) + (CF*nq^2*sthetat^4*Tf^2)/(12*ep^3) - (CF*nq^2*sthetat^4*Tf^2)/(12*ep^2) - 
   (CF*nq^2*sthetat^4*Tf^2)/(12*ep) - (CF*Mst1^2*nq^2*sthetat^4*Tf^2)/(12*ep^3*Mst2^2) + (CF*Mst1^2*nq^2*sthetat^4*Tf^2)/(12*ep^2*Mst2^2) + 
   (CF*Mst1^2*nq^2*sthetat^4*Tf^2)/(12*ep*Mst2^2) + (CF*nq*nt*sthetat^4*Tf^2)/(6*ep^3) - (CF*nq*nt*sthetat^4*Tf^2)/(6*ep^2) - 
   (CF*nq*nt*sthetat^4*Tf^2)/(6*ep) - (CF*Mst1^2*nq*nt*sthetat^4*Tf^2)/(6*ep^3*Mst2^2) + (CF*Mst1^2*nq*nt*sthetat^4*Tf^2)/(6*ep^2*Mst2^2) + 
   (CF*Mst1^2*nq*nt*sthetat^4*Tf^2)/(6*ep*Mst2^2) + (CF*nt^2*sthetat^4*Tf^2)/(12*ep^3) - (CF*nt^2*sthetat^4*Tf^2)/(12*ep^2) - 
   (CF*nt^2*sthetat^4*Tf^2)/(12*ep) - (CF*Mst1^2*nt^2*sthetat^4*Tf^2)/(12*ep^3*Mst2^2) + (CF*Mst1^2*nt^2*sthetat^4*Tf^2)/(12*ep^2*Mst2^2) + 
   (CF*Mst1^2*nt^2*sthetat^4*Tf^2)/(12*ep*Mst2^2) + (3*CA*CF*Mfg^2*nq*Tf*z3)/(ep*Mst2^2) - (3*CF^2*Mfg^2*nq*Tf*z3)/(ep*Mst2^2) + 
   (CA*CF*Mt^2*nq*Tf*z3)/(2*ep*Mst2^2) - (CF^2*Mt^2*nq*Tf*z3)/(2*ep*Mst2^2) + (3*CA*CF*Mfg^2*nt*Tf*z3)/(ep*Mst2^2) - 
   (3*CF^2*Mfg^2*nt*Tf*z3)/(ep*Mst2^2) + (CA*CF*Mt^2*nt*Tf*z3)/(2*ep*Mst2^2) - (CF^2*Mt^2*nt*Tf*z3)/(2*ep*Mst2^2) + 
   (3*CA*CF*cthetat*Mfg*Mt*nq*sthetat*Tf*z3)/(ep*Mst2^2) - (3*CF^2*cthetat*Mfg*Mt*nq*sthetat*Tf*z3)/(ep*Mst2^2) + 
   (3*CA*CF*cthetat*Mfg*Mt*nt*sthetat*Tf*z3)/(ep*Mst2^2) - (3*CF^2*cthetat*Mfg*Mt*nt*sthetat*Tf*z3)/(ep*Mst2^2) + 
   (CA*CF*nq*sthetat^2*Tf*z3)/(2*ep) - (CF^2*nq*sthetat^2*Tf*z3)/(2*ep) - (CA*CF*Mst1^2*nq*sthetat^2*Tf*z3)/(2*ep*Mst2^2) + 
   (CF^2*Mst1^2*nq*sthetat^2*Tf*z3)/(2*ep*Mst2^2) + (CA*CF*nt*sthetat^2*Tf*z3)/(2*ep) - (CF^2*nt*sthetat^2*Tf*z3)/(2*ep) - 
   (CA*CF*Mst1^2*nt*sthetat^2*Tf*z3)/(2*ep*Mst2^2) + (CF^2*Mst1^2*nt*sthetat^2*Tf*z3)/(2*ep*Mst2^2) - (CA*CF*nq*sthetat^4*Tf*z3)/(2*ep) + 
   (CF^2*nq*sthetat^4*Tf*z3)/(2*ep) + (CA*CF*Mst1^2*nq*sthetat^4*Tf*z3)/(2*ep*Mst2^2) - (CF^2*Mst1^2*nq*sthetat^4*Tf*z3)/(2*ep*Mst2^2) - 
   (CA*CF*nt*sthetat^4*Tf*z3)/(2*ep) + (CF^2*nt*sthetat^4*Tf*z3)/(2*ep) + (CA*CF*Mst1^2*nt*sthetat^4*Tf*z3)/(2*ep*Mst2^2) - 
   (CF^2*Mst1^2*nt*sthetat^4*Tf*z3)/(2*ep*Mst2^2)));

dThetat = ((as*(2*CF*Mfg*Mt - CF*cthetat*Mst1^2*sthetat + CF*cthetat*Mst2^2*sthetat - 4*CF*Mfg*Mt*sthetat^2 + 2*CF*cthetat*Mst1^2*sthetat^3 - 
    2*CF*cthetat*Mst2^2*sthetat^3))/(2*ep*(Mst1^2 - Mst2^2))
 + (as^2*(-12*CA*CF*Mfg*Mst1^2*Mt - 8*CF^2*Mfg*Mst1^2*Mt + 12*CA*CF*ep*Mfg*Mst1^2*Mt - 8*CF^2*ep*Mfg*Mst1^2*Mt + 12*CA*CF*Mfg*Mst2^2*Mt + 
    8*CF^2*Mfg*Mst2^2*Mt - 12*CA*CF*ep*Mfg*Mst2^2*Mt + 8*CF^2*ep*Mfg*Mst2^2*Mt + 3*CA*CF*cthetat*Mst1^4*sthetat + 
    2*CF^2*cthetat*Mst1^4*sthetat - 3*CA*CF*cthetat*ep*Mst1^4*sthetat + 2*CF^2*cthetat*ep*Mst1^4*sthetat - 
    6*CA*CF*cthetat*Mst1^2*Mst2^2*sthetat - 4*CF^2*cthetat*Mst1^2*Mst2^2*sthetat + 6*CA*CF*cthetat*ep*Mst1^2*Mst2^2*sthetat - 
    4*CF^2*cthetat*ep*Mst1^2*Mst2^2*sthetat + 3*CA*CF*cthetat*Mst2^4*sthetat + 2*CF^2*cthetat*Mst2^4*sthetat - 
    3*CA*CF*cthetat*ep*Mst2^4*sthetat + 2*CF^2*cthetat*ep*Mst2^4*sthetat - 64*CF^2*cthetat*Mfg^2*Mt^2*sthetat + 
    24*CA*CF*Mfg*Mst1^2*Mt*sthetat^2 + 80*CF^2*Mfg*Mst1^2*Mt*sthetat^2 - 24*CA*CF*ep*Mfg*Mst1^2*Mt*sthetat^2 + 
    16*CF^2*ep*Mfg*Mst1^2*Mt*sthetat^2 - 24*CA*CF*Mfg*Mst2^2*Mt*sthetat^2 - 80*CF^2*Mfg*Mst2^2*Mt*sthetat^2 + 
    24*CA*CF*ep*Mfg*Mst2^2*Mt*sthetat^2 - 16*CF^2*ep*Mfg*Mst2^2*Mt*sthetat^2 - 6*CA*CF*cthetat*Mst1^4*sthetat^3 - 
    20*CF^2*cthetat*Mst1^4*sthetat^3 + 6*CA*CF*cthetat*ep*Mst1^4*sthetat^3 - 4*CF^2*cthetat*ep*Mst1^4*sthetat^3 + 
    12*CA*CF*cthetat*Mst1^2*Mst2^2*sthetat^3 + 40*CF^2*cthetat*Mst1^2*Mst2^2*sthetat^3 - 12*CA*CF*cthetat*ep*Mst1^2*Mst2^2*sthetat^3 + 
    8*CF^2*cthetat*ep*Mst1^2*Mst2^2*sthetat^3 - 6*CA*CF*cthetat*Mst2^4*sthetat^3 - 20*CF^2*cthetat*Mst2^4*sthetat^3 + 
    6*CA*CF*cthetat*ep*Mst2^4*sthetat^3 - 4*CF^2*cthetat*ep*Mst2^4*sthetat^3 + 128*CF^2*cthetat*Mfg^2*Mt^2*sthetat^3 - 
    192*CF^2*Mfg*Mst1^2*Mt*sthetat^4 + 192*CF^2*Mfg*Mst2^2*Mt*sthetat^4 + 48*CF^2*cthetat*Mst1^4*sthetat^5 - 
    96*CF^2*cthetat*Mst1^2*Mst2^2*sthetat^5 + 48*CF^2*cthetat*Mst2^4*sthetat^5 + 128*CF^2*Mfg*Mst1^2*Mt*sthetat^6 - 
    128*CF^2*Mfg*Mst2^2*Mt*sthetat^6 - 32*CF^2*cthetat*Mst1^4*sthetat^7 + 64*CF^2*cthetat*Mst1^2*Mst2^2*sthetat^7 - 
    32*CF^2*cthetat*Mst2^4*sthetat^7 + 8*CF*Mfg*Mst1^2*Mt*nq*Tf - 8*CF*ep*Mfg*Mst1^2*Mt*nq*Tf - 8*CF*Mfg*Mst2^2*Mt*nq*Tf + 
    8*CF*ep*Mfg*Mst2^2*Mt*nq*Tf + 8*CF*Mfg*Mst1^2*Mt*nt*Tf - 8*CF*ep*Mfg*Mst1^2*Mt*nt*Tf - 8*CF*Mfg*Mst2^2*Mt*nt*Tf + 
    8*CF*ep*Mfg*Mst2^2*Mt*nt*Tf - 2*CF*cthetat*Mst1^4*nq*sthetat*Tf + 2*CF*cthetat*ep*Mst1^4*nq*sthetat*Tf + 
    4*CF*cthetat*Mst1^2*Mst2^2*nq*sthetat*Tf - 4*CF*cthetat*ep*Mst1^2*Mst2^2*nq*sthetat*Tf - 2*CF*cthetat*Mst2^4*nq*sthetat*Tf + 
    2*CF*cthetat*ep*Mst2^4*nq*sthetat*Tf - 2*CF*cthetat*Mst1^4*nt*sthetat*Tf + 2*CF*cthetat*ep*Mst1^4*nt*sthetat*Tf + 
    4*CF*cthetat*Mst1^2*Mst2^2*nt*sthetat*Tf - 4*CF*cthetat*ep*Mst1^2*Mst2^2*nt*sthetat*Tf - 2*CF*cthetat*Mst2^4*nt*sthetat*Tf + 
    2*CF*cthetat*ep*Mst2^4*nt*sthetat*Tf - 16*CF*Mfg*Mst1^2*Mt*nq*sthetat^2*Tf + 16*CF*ep*Mfg*Mst1^2*Mt*nq*sthetat^2*Tf + 
    16*CF*Mfg*Mst2^2*Mt*nq*sthetat^2*Tf - 16*CF*ep*Mfg*Mst2^2*Mt*nq*sthetat^2*Tf - 16*CF*Mfg*Mst1^2*Mt*nt*sthetat^2*Tf + 
    16*CF*ep*Mfg*Mst1^2*Mt*nt*sthetat^2*Tf + 16*CF*Mfg*Mst2^2*Mt*nt*sthetat^2*Tf - 16*CF*ep*Mfg*Mst2^2*Mt*nt*sthetat^2*Tf + 
    4*CF*cthetat*Mst1^4*nq*sthetat^3*Tf - 4*CF*cthetat*ep*Mst1^4*nq*sthetat^3*Tf - 8*CF*cthetat*Mst1^2*Mst2^2*nq*sthetat^3*Tf + 
    8*CF*cthetat*ep*Mst1^2*Mst2^2*nq*sthetat^3*Tf + 4*CF*cthetat*Mst2^4*nq*sthetat^3*Tf - 4*CF*cthetat*ep*Mst2^4*nq*sthetat^3*Tf + 
    4*CF*cthetat*Mst1^4*nt*sthetat^3*Tf - 4*CF*cthetat*ep*Mst1^4*nt*sthetat^3*Tf - 8*CF*cthetat*Mst1^2*Mst2^2*nt*sthetat^3*Tf + 
    8*CF*cthetat*ep*Mst1^2*Mst2^2*nt*sthetat^3*Tf + 4*CF*cthetat*Mst2^4*nt*sthetat^3*Tf - 4*CF*cthetat*ep*Mst2^4*nt*sthetat^3*Tf))/
  (16*ep^2*(Mst1^2 - Mst2^2)^2)
 + as^3*((9*CA^2*CF*Mfg*Mt)/(16*ep^3*(Mst1^2 - Mst2^2)) + (9*CA*CF^2*Mfg*Mt)/(16*ep^3*(Mst1^2 - Mst2^2)) + 
   (CF^3*Mfg*Mt)/(8*ep^3*(Mst1^2 - Mst2^2)) - (15*CA^2*CF*Mfg*Mt)/(16*ep^2*(Mst1^2 - Mst2^2)) - 
   (3*CA*CF^2*Mfg*Mt)/(16*ep^2*(Mst1^2 - Mst2^2)) + (3*CF^3*Mfg*Mt)/(8*ep^2*(Mst1^2 - Mst2^2)) + 
   (3*CA^2*CF*Mfg*Mt)/(8*ep*(Mst1^2 - Mst2^2)) - (3*CA*CF^2*Mfg*Mt)/(8*ep*(Mst1^2 - Mst2^2)) + (CF^3*Mfg*Mt)/(2*ep*(Mst1^2 - Mst2^2)) - 
   (4*CF^3*Mfg^3*Mt^3)/(3*ep^3*(Mst1^2 - Mst2^2)^3) - (3*CA^2*CF*cthetat*Mst1^2*sthetat)/(32*ep^3*(Mst1^2 - Mst2^2)) - 
   (3*CA*CF^2*cthetat*Mst1^2*sthetat)/(32*ep^3*(Mst1^2 - Mst2^2)) - (CF^3*cthetat*Mst1^2*sthetat)/(48*ep^3*(Mst1^2 - Mst2^2)) + 
   (5*CA^2*CF*cthetat*Mst1^2*sthetat)/(32*ep^2*(Mst1^2 - Mst2^2)) + (CA*CF^2*cthetat*Mst1^2*sthetat)/(32*ep^2*(Mst1^2 - Mst2^2)) - 
   (CF^3*cthetat*Mst1^2*sthetat)/(16*ep^2*(Mst1^2 - Mst2^2)) - (CA^2*CF*cthetat*Mst1^2*sthetat)/(16*ep*(Mst1^2 - Mst2^2)) + 
   (CA*CF^2*cthetat*Mst1^2*sthetat)/(16*ep*(Mst1^2 - Mst2^2)) - (CF^3*cthetat*Mst1^2*sthetat)/(12*ep*(Mst1^2 - Mst2^2)) + 
   (3*CA^2*CF*cthetat*Mst2^2*sthetat)/(32*ep^3*(Mst1^2 - Mst2^2)) + (3*CA*CF^2*cthetat*Mst2^2*sthetat)/(32*ep^3*(Mst1^2 - Mst2^2)) + 
   (CF^3*cthetat*Mst2^2*sthetat)/(48*ep^3*(Mst1^2 - Mst2^2)) - (5*CA^2*CF*cthetat*Mst2^2*sthetat)/(32*ep^2*(Mst1^2 - Mst2^2)) - 
   (CA*CF^2*cthetat*Mst2^2*sthetat)/(32*ep^2*(Mst1^2 - Mst2^2)) + (CF^3*cthetat*Mst2^2*sthetat)/(16*ep^2*(Mst1^2 - Mst2^2)) + 
   (CA^2*CF*cthetat*Mst2^2*sthetat)/(16*ep*(Mst1^2 - Mst2^2)) - (CA*CF^2*cthetat*Mst2^2*sthetat)/(16*ep*(Mst1^2 - Mst2^2)) + 
   (CF^3*cthetat*Mst2^2*sthetat)/(12*ep*(Mst1^2 - Mst2^2)) + (6*CA*CF^2*cthetat*Mfg^2*Mt^2*sthetat)/(ep^3*(Mst1^2 - Mst2^2)^2) + 
   (6*CF^3*cthetat*Mfg^2*Mt^2*sthetat)/(ep^3*(Mst1^2 - Mst2^2)^2) - (6*CA*CF^2*cthetat*Mfg^2*Mt^2*sthetat)/(ep^2*(Mst1^2 - Mst2^2)^2) + 
   (4*CF^3*cthetat*Mfg^2*Mt^2*sthetat)/(ep^2*(Mst1^2 - Mst2^2)^2) - (9*CA^2*CF*Mfg*Mt*sthetat^2)/(8*ep^3*(Mst1^2 - Mst2^2)) - 
   (45*CA*CF^2*Mfg*Mt*sthetat^2)/(8*ep^3*(Mst1^2 - Mst2^2)) - (17*CF^3*Mfg*Mt*sthetat^2)/(4*ep^3*(Mst1^2 - Mst2^2)) + 
   (15*CA^2*CF*Mfg*Mt*sthetat^2)/(8*ep^2*(Mst1^2 - Mst2^2)) + (39*CA*CF^2*Mfg*Mt*sthetat^2)/(8*ep^2*(Mst1^2 - Mst2^2)) - 
   (15*CF^3*Mfg*Mt*sthetat^2)/(4*ep^2*(Mst1^2 - Mst2^2)) - (3*CA^2*CF*Mfg*Mt*sthetat^2)/(4*ep*(Mst1^2 - Mst2^2)) + 
   (3*CA*CF^2*Mfg*Mt*sthetat^2)/(4*ep*(Mst1^2 - Mst2^2)) - (CF^3*Mfg*Mt*sthetat^2)/(ep*(Mst1^2 - Mst2^2)) + 
   (24*CF^3*Mfg^3*Mt^3*sthetat^2)/(ep^3*(Mst1^2 - Mst2^2)^3) + (3*CA^2*CF*cthetat*Mst1^2*sthetat^3)/(16*ep^3*(Mst1^2 - Mst2^2)) + 
   (15*CA*CF^2*cthetat*Mst1^2*sthetat^3)/(16*ep^3*(Mst1^2 - Mst2^2)) + (17*CF^3*cthetat*Mst1^2*sthetat^3)/(24*ep^3*(Mst1^2 - Mst2^2)) - 
   (5*CA^2*CF*cthetat*Mst1^2*sthetat^3)/(16*ep^2*(Mst1^2 - Mst2^2)) - (13*CA*CF^2*cthetat*Mst1^2*sthetat^3)/(16*ep^2*(Mst1^2 - Mst2^2)) + 
   (5*CF^3*cthetat*Mst1^2*sthetat^3)/(8*ep^2*(Mst1^2 - Mst2^2)) + (CA^2*CF*cthetat*Mst1^2*sthetat^3)/(8*ep*(Mst1^2 - Mst2^2)) - 
   (CA*CF^2*cthetat*Mst1^2*sthetat^3)/(8*ep*(Mst1^2 - Mst2^2)) + (CF^3*cthetat*Mst1^2*sthetat^3)/(6*ep*(Mst1^2 - Mst2^2)) - 
   (3*CA^2*CF*cthetat*Mst2^2*sthetat^3)/(16*ep^3*(Mst1^2 - Mst2^2)) - (15*CA*CF^2*cthetat*Mst2^2*sthetat^3)/(16*ep^3*(Mst1^2 - Mst2^2)) - 
   (17*CF^3*cthetat*Mst2^2*sthetat^3)/(24*ep^3*(Mst1^2 - Mst2^2)) + (5*CA^2*CF*cthetat*Mst2^2*sthetat^3)/(16*ep^2*(Mst1^2 - Mst2^2)) + 
   (13*CA*CF^2*cthetat*Mst2^2*sthetat^3)/(16*ep^2*(Mst1^2 - Mst2^2)) - (5*CF^3*cthetat*Mst2^2*sthetat^3)/(8*ep^2*(Mst1^2 - Mst2^2)) - 
   (CA^2*CF*cthetat*Mst2^2*sthetat^3)/(8*ep*(Mst1^2 - Mst2^2)) + (CA*CF^2*cthetat*Mst2^2*sthetat^3)/(8*ep*(Mst1^2 - Mst2^2)) - 
   (CF^3*cthetat*Mst2^2*sthetat^3)/(6*ep*(Mst1^2 - Mst2^2)) - (12*CA*CF^2*cthetat*Mfg^2*Mt^2*sthetat^3)/(ep^3*(Mst1^2 - Mst2^2)^2) - 
   (44*CF^3*cthetat*Mfg^2*Mt^2*sthetat^3)/(ep^3*(Mst1^2 - Mst2^2)^2) + (12*CA*CF^2*cthetat*Mfg^2*Mt^2*sthetat^3)/
    (ep^2*(Mst1^2 - Mst2^2)^2) - (8*CF^3*cthetat*Mfg^2*Mt^2*sthetat^3)/(ep^2*(Mst1^2 - Mst2^2)^2) + 
   (27*CA*CF^2*Mfg*Mt*sthetat^4)/(2*ep^3*(Mst1^2 - Mst2^2)) + (28*CF^3*Mfg*Mt*sthetat^4)/(ep^3*(Mst1^2 - Mst2^2)) - 
   (27*CA*CF^2*Mfg*Mt*sthetat^4)/(2*ep^2*(Mst1^2 - Mst2^2)) + (9*CF^3*Mfg*Mt*sthetat^4)/(ep^2*(Mst1^2 - Mst2^2)) - 
   (64*CF^3*Mfg^3*Mt^3*sthetat^4)/(ep^3*(Mst1^2 - Mst2^2)^3) - (9*CA*CF^2*cthetat*Mst1^2*sthetat^5)/(4*ep^3*(Mst1^2 - Mst2^2)) - 
   (14*CF^3*cthetat*Mst1^2*sthetat^5)/(3*ep^3*(Mst1^2 - Mst2^2)) + (9*CA*CF^2*cthetat*Mst1^2*sthetat^5)/(4*ep^2*(Mst1^2 - Mst2^2)) - 
   (3*CF^3*cthetat*Mst1^2*sthetat^5)/(2*ep^2*(Mst1^2 - Mst2^2)) + (9*CA*CF^2*cthetat*Mst2^2*sthetat^5)/(4*ep^3*(Mst1^2 - Mst2^2)) + 
   (14*CF^3*cthetat*Mst2^2*sthetat^5)/(3*ep^3*(Mst1^2 - Mst2^2)) - (9*CA*CF^2*cthetat*Mst2^2*sthetat^5)/(4*ep^2*(Mst1^2 - Mst2^2)) + 
   (3*CF^3*cthetat*Mst2^2*sthetat^5)/(2*ep^2*(Mst1^2 - Mst2^2)) + (96*CF^3*cthetat*Mfg^2*Mt^2*sthetat^5)/(ep^3*(Mst1^2 - Mst2^2)^2) - 
   (9*CA*CF^2*Mfg*Mt*sthetat^6)/(ep^3*(Mst1^2 - Mst2^2)) - (72*CF^3*Mfg*Mt*sthetat^6)/(ep^3*(Mst1^2 - Mst2^2)) + 
   (9*CA*CF^2*Mfg*Mt*sthetat^6)/(ep^2*(Mst1^2 - Mst2^2)) - (6*CF^3*Mfg*Mt*sthetat^6)/(ep^2*(Mst1^2 - Mst2^2)) + 
   (128*CF^3*Mfg^3*Mt^3*sthetat^6)/(3*ep^3*(Mst1^2 - Mst2^2)^3) + (3*CA*CF^2*cthetat*Mst1^2*sthetat^7)/(2*ep^3*(Mst1^2 - Mst2^2)) + 
   (12*CF^3*cthetat*Mst1^2*sthetat^7)/(ep^3*(Mst1^2 - Mst2^2)) - (3*CA*CF^2*cthetat*Mst1^2*sthetat^7)/(2*ep^2*(Mst1^2 - Mst2^2)) + 
   (CF^3*cthetat*Mst1^2*sthetat^7)/(ep^2*(Mst1^2 - Mst2^2)) - (3*CA*CF^2*cthetat*Mst2^2*sthetat^7)/(2*ep^3*(Mst1^2 - Mst2^2)) - 
   (12*CF^3*cthetat*Mst2^2*sthetat^7)/(ep^3*(Mst1^2 - Mst2^2)) + (3*CA*CF^2*cthetat*Mst2^2*sthetat^7)/(2*ep^2*(Mst1^2 - Mst2^2)) - 
   (CF^3*cthetat*Mst2^2*sthetat^7)/(ep^2*(Mst1^2 - Mst2^2)) - (64*CF^3*cthetat*Mfg^2*Mt^2*sthetat^7)/(ep^3*(Mst1^2 - Mst2^2)^2) + 
   (80*CF^3*Mfg*Mt*sthetat^8)/(ep^3*(Mst1^2 - Mst2^2)) - (40*CF^3*cthetat*Mst1^2*sthetat^9)/(3*ep^3*(Mst1^2 - Mst2^2)) + 
   (40*CF^3*cthetat*Mst2^2*sthetat^9)/(3*ep^3*(Mst1^2 - Mst2^2)) - (32*CF^3*Mfg*Mt*sthetat^10)/(ep^3*(Mst1^2 - Mst2^2)) + 
   (16*CF^3*cthetat*Mst1^2*sthetat^11)/(3*ep^3*(Mst1^2 - Mst2^2)) - (16*CF^3*cthetat*Mst2^2*sthetat^11)/(3*ep^3*(Mst1^2 - Mst2^2)) - 
   (3*CA*CF*Mfg*Mt*nq*Tf)/(4*ep^3*(Mst1^2 - Mst2^2)) - (3*CF^2*Mfg*Mt*nq*Tf)/(8*ep^3*(Mst1^2 - Mst2^2)) + 
   (CA*CF*Mfg*Mt*nq*Tf)/(ep^2*(Mst1^2 - Mst2^2)) + (5*CF^2*Mfg*Mt*nq*Tf)/(8*ep^2*(Mst1^2 - Mst2^2)) + 
   (CA*CF*Mfg*Mt*nq*Tf)/(8*ep*(Mst1^2 - Mst2^2)) - (CF^2*Mfg*Mt*nq*Tf)/(ep*(Mst1^2 - Mst2^2)) - 
   (3*CA*CF*Mfg*Mt*nt*Tf)/(4*ep^3*(Mst1^2 - Mst2^2)) - (3*CF^2*Mfg*Mt*nt*Tf)/(8*ep^3*(Mst1^2 - Mst2^2)) + 
   (CA*CF*Mfg*Mt*nt*Tf)/(ep^2*(Mst1^2 - Mst2^2)) + (5*CF^2*Mfg*Mt*nt*Tf)/(8*ep^2*(Mst1^2 - Mst2^2)) + 
   (CA*CF*Mfg*Mt*nt*Tf)/(8*ep*(Mst1^2 - Mst2^2)) - (CF^2*Mfg*Mt*nt*Tf)/(ep*(Mst1^2 - Mst2^2)) + 
   (CA*CF*cthetat*Mst1^2*nq*sthetat*Tf)/(8*ep^3*(Mst1^2 - Mst2^2)) + (CF^2*cthetat*Mst1^2*nq*sthetat*Tf)/(16*ep^3*(Mst1^2 - Mst2^2)) - 
   (CA*CF*cthetat*Mst1^2*nq*sthetat*Tf)/(6*ep^2*(Mst1^2 - Mst2^2)) - (5*CF^2*cthetat*Mst1^2*nq*sthetat*Tf)/(48*ep^2*(Mst1^2 - Mst2^2)) - 
   (CA*CF*cthetat*Mst1^2*nq*sthetat*Tf)/(48*ep*(Mst1^2 - Mst2^2)) + (CF^2*cthetat*Mst1^2*nq*sthetat*Tf)/(6*ep*(Mst1^2 - Mst2^2)) - 
   (CA*CF*cthetat*Mst2^2*nq*sthetat*Tf)/(8*ep^3*(Mst1^2 - Mst2^2)) - (CF^2*cthetat*Mst2^2*nq*sthetat*Tf)/(16*ep^3*(Mst1^2 - Mst2^2)) + 
   (CA*CF*cthetat*Mst2^2*nq*sthetat*Tf)/(6*ep^2*(Mst1^2 - Mst2^2)) + (5*CF^2*cthetat*Mst2^2*nq*sthetat*Tf)/(48*ep^2*(Mst1^2 - Mst2^2)) + 
   (CA*CF*cthetat*Mst2^2*nq*sthetat*Tf)/(48*ep*(Mst1^2 - Mst2^2)) - (CF^2*cthetat*Mst2^2*nq*sthetat*Tf)/(6*ep*(Mst1^2 - Mst2^2)) - 
   (4*CF^2*cthetat*Mfg^2*Mt^2*nq*sthetat*Tf)/(ep^3*(Mst1^2 - Mst2^2)^2) + (4*CF^2*cthetat*Mfg^2*Mt^2*nq*sthetat*Tf)/
    (ep^2*(Mst1^2 - Mst2^2)^2) + (CA*CF*cthetat*Mst1^2*nt*sthetat*Tf)/(8*ep^3*(Mst1^2 - Mst2^2)) + 
   (CF^2*cthetat*Mst1^2*nt*sthetat*Tf)/(16*ep^3*(Mst1^2 - Mst2^2)) - (CA*CF*cthetat*Mst1^2*nt*sthetat*Tf)/(6*ep^2*(Mst1^2 - Mst2^2)) - 
   (5*CF^2*cthetat*Mst1^2*nt*sthetat*Tf)/(48*ep^2*(Mst1^2 - Mst2^2)) - (CA*CF*cthetat*Mst1^2*nt*sthetat*Tf)/(48*ep*(Mst1^2 - Mst2^2)) + 
   (CF^2*cthetat*Mst1^2*nt*sthetat*Tf)/(6*ep*(Mst1^2 - Mst2^2)) - (CA*CF*cthetat*Mst2^2*nt*sthetat*Tf)/(8*ep^3*(Mst1^2 - Mst2^2)) - 
   (CF^2*cthetat*Mst2^2*nt*sthetat*Tf)/(16*ep^3*(Mst1^2 - Mst2^2)) + (CA*CF*cthetat*Mst2^2*nt*sthetat*Tf)/(6*ep^2*(Mst1^2 - Mst2^2)) + 
   (5*CF^2*cthetat*Mst2^2*nt*sthetat*Tf)/(48*ep^2*(Mst1^2 - Mst2^2)) + (CA*CF*cthetat*Mst2^2*nt*sthetat*Tf)/(48*ep*(Mst1^2 - Mst2^2)) - 
   (CF^2*cthetat*Mst2^2*nt*sthetat*Tf)/(6*ep*(Mst1^2 - Mst2^2)) - (4*CF^2*cthetat*Mfg^2*Mt^2*nt*sthetat*Tf)/(ep^3*(Mst1^2 - Mst2^2)^2) + 
   (4*CF^2*cthetat*Mfg^2*Mt^2*nt*sthetat*Tf)/(ep^2*(Mst1^2 - Mst2^2)^2) + (3*CA*CF*Mfg*Mt*nq*sthetat^2*Tf)/(2*ep^3*(Mst1^2 - Mst2^2)) + 
   (15*CF^2*Mfg*Mt*nq*sthetat^2*Tf)/(4*ep^3*(Mst1^2 - Mst2^2)) - (2*CA*CF*Mfg*Mt*nq*sthetat^2*Tf)/(ep^2*(Mst1^2 - Mst2^2)) - 
   (17*CF^2*Mfg*Mt*nq*sthetat^2*Tf)/(4*ep^2*(Mst1^2 - Mst2^2)) - (CA*CF*Mfg*Mt*nq*sthetat^2*Tf)/(4*ep*(Mst1^2 - Mst2^2)) + 
   (2*CF^2*Mfg*Mt*nq*sthetat^2*Tf)/(ep*(Mst1^2 - Mst2^2)) + (3*CA*CF*Mfg*Mt*nt*sthetat^2*Tf)/(2*ep^3*(Mst1^2 - Mst2^2)) + 
   (15*CF^2*Mfg*Mt*nt*sthetat^2*Tf)/(4*ep^3*(Mst1^2 - Mst2^2)) - (2*CA*CF*Mfg*Mt*nt*sthetat^2*Tf)/(ep^2*(Mst1^2 - Mst2^2)) - 
   (17*CF^2*Mfg*Mt*nt*sthetat^2*Tf)/(4*ep^2*(Mst1^2 - Mst2^2)) - (CA*CF*Mfg*Mt*nt*sthetat^2*Tf)/(4*ep*(Mst1^2 - Mst2^2)) + 
   (2*CF^2*Mfg*Mt*nt*sthetat^2*Tf)/(ep*(Mst1^2 - Mst2^2)) - (CA*CF*cthetat*Mst1^2*nq*sthetat^3*Tf)/(4*ep^3*(Mst1^2 - Mst2^2)) - 
   (5*CF^2*cthetat*Mst1^2*nq*sthetat^3*Tf)/(8*ep^3*(Mst1^2 - Mst2^2)) + (CA*CF*cthetat*Mst1^2*nq*sthetat^3*Tf)/(3*ep^2*(Mst1^2 - Mst2^2)) + 
   (17*CF^2*cthetat*Mst1^2*nq*sthetat^3*Tf)/(24*ep^2*(Mst1^2 - Mst2^2)) + (CA*CF*cthetat*Mst1^2*nq*sthetat^3*Tf)/(24*ep*(Mst1^2 - Mst2^2)) - 
   (CF^2*cthetat*Mst1^2*nq*sthetat^3*Tf)/(3*ep*(Mst1^2 - Mst2^2)) + (CA*CF*cthetat*Mst2^2*nq*sthetat^3*Tf)/(4*ep^3*(Mst1^2 - Mst2^2)) + 
   (5*CF^2*cthetat*Mst2^2*nq*sthetat^3*Tf)/(8*ep^3*(Mst1^2 - Mst2^2)) - (CA*CF*cthetat*Mst2^2*nq*sthetat^3*Tf)/(3*ep^2*(Mst1^2 - Mst2^2)) - 
   (17*CF^2*cthetat*Mst2^2*nq*sthetat^3*Tf)/(24*ep^2*(Mst1^2 - Mst2^2)) - (CA*CF*cthetat*Mst2^2*nq*sthetat^3*Tf)/(24*ep*(Mst1^2 - Mst2^2)) + 
   (CF^2*cthetat*Mst2^2*nq*sthetat^3*Tf)/(3*ep*(Mst1^2 - Mst2^2)) + (8*CF^2*cthetat*Mfg^2*Mt^2*nq*sthetat^3*Tf)/(ep^3*(Mst1^2 - Mst2^2)^2) - 
   (8*CF^2*cthetat*Mfg^2*Mt^2*nq*sthetat^3*Tf)/(ep^2*(Mst1^2 - Mst2^2)^2) - (CA*CF*cthetat*Mst1^2*nt*sthetat^3*Tf)/
    (4*ep^3*(Mst1^2 - Mst2^2)) - (5*CF^2*cthetat*Mst1^2*nt*sthetat^3*Tf)/(8*ep^3*(Mst1^2 - Mst2^2)) + 
   (CA*CF*cthetat*Mst1^2*nt*sthetat^3*Tf)/(3*ep^2*(Mst1^2 - Mst2^2)) + (17*CF^2*cthetat*Mst1^2*nt*sthetat^3*Tf)/
    (24*ep^2*(Mst1^2 - Mst2^2)) + (CA*CF*cthetat*Mst1^2*nt*sthetat^3*Tf)/(24*ep*(Mst1^2 - Mst2^2)) - 
   (CF^2*cthetat*Mst1^2*nt*sthetat^3*Tf)/(3*ep*(Mst1^2 - Mst2^2)) + (CA*CF*cthetat*Mst2^2*nt*sthetat^3*Tf)/(4*ep^3*(Mst1^2 - Mst2^2)) + 
   (5*CF^2*cthetat*Mst2^2*nt*sthetat^3*Tf)/(8*ep^3*(Mst1^2 - Mst2^2)) - (CA*CF*cthetat*Mst2^2*nt*sthetat^3*Tf)/(3*ep^2*(Mst1^2 - Mst2^2)) - 
   (17*CF^2*cthetat*Mst2^2*nt*sthetat^3*Tf)/(24*ep^2*(Mst1^2 - Mst2^2)) - (CA*CF*cthetat*Mst2^2*nt*sthetat^3*Tf)/(24*ep*(Mst1^2 - Mst2^2)) + 
   (CF^2*cthetat*Mst2^2*nt*sthetat^3*Tf)/(3*ep*(Mst1^2 - Mst2^2)) + (8*CF^2*cthetat*Mfg^2*Mt^2*nt*sthetat^3*Tf)/(ep^3*(Mst1^2 - Mst2^2)^2) - 
   (8*CF^2*cthetat*Mfg^2*Mt^2*nt*sthetat^3*Tf)/(ep^2*(Mst1^2 - Mst2^2)^2) - (9*CF^2*Mfg*Mt*nq*sthetat^4*Tf)/(ep^3*(Mst1^2 - Mst2^2)) + 
   (9*CF^2*Mfg*Mt*nq*sthetat^4*Tf)/(ep^2*(Mst1^2 - Mst2^2)) - (9*CF^2*Mfg*Mt*nt*sthetat^4*Tf)/(ep^3*(Mst1^2 - Mst2^2)) + 
   (9*CF^2*Mfg*Mt*nt*sthetat^4*Tf)/(ep^2*(Mst1^2 - Mst2^2)) + (3*CF^2*cthetat*Mst1^2*nq*sthetat^5*Tf)/(2*ep^3*(Mst1^2 - Mst2^2)) - 
   (3*CF^2*cthetat*Mst1^2*nq*sthetat^5*Tf)/(2*ep^2*(Mst1^2 - Mst2^2)) - (3*CF^2*cthetat*Mst2^2*nq*sthetat^5*Tf)/(2*ep^3*(Mst1^2 - Mst2^2)) + 
   (3*CF^2*cthetat*Mst2^2*nq*sthetat^5*Tf)/(2*ep^2*(Mst1^2 - Mst2^2)) + (3*CF^2*cthetat*Mst1^2*nt*sthetat^5*Tf)/(2*ep^3*(Mst1^2 - Mst2^2)) - 
   (3*CF^2*cthetat*Mst1^2*nt*sthetat^5*Tf)/(2*ep^2*(Mst1^2 - Mst2^2)) - (3*CF^2*cthetat*Mst2^2*nt*sthetat^5*Tf)/(2*ep^3*(Mst1^2 - Mst2^2)) + 
   (3*CF^2*cthetat*Mst2^2*nt*sthetat^5*Tf)/(2*ep^2*(Mst1^2 - Mst2^2)) + (6*CF^2*Mfg*Mt*nq*sthetat^6*Tf)/(ep^3*(Mst1^2 - Mst2^2)) - 
   (6*CF^2*Mfg*Mt*nq*sthetat^6*Tf)/(ep^2*(Mst1^2 - Mst2^2)) + (6*CF^2*Mfg*Mt*nt*sthetat^6*Tf)/(ep^3*(Mst1^2 - Mst2^2)) - 
   (6*CF^2*Mfg*Mt*nt*sthetat^6*Tf)/(ep^2*(Mst1^2 - Mst2^2)) - (CF^2*cthetat*Mst1^2*nq*sthetat^7*Tf)/(ep^3*(Mst1^2 - Mst2^2)) + 
   (CF^2*cthetat*Mst1^2*nq*sthetat^7*Tf)/(ep^2*(Mst1^2 - Mst2^2)) + (CF^2*cthetat*Mst2^2*nq*sthetat^7*Tf)/(ep^3*(Mst1^2 - Mst2^2)) - 
   (CF^2*cthetat*Mst2^2*nq*sthetat^7*Tf)/(ep^2*(Mst1^2 - Mst2^2)) - (CF^2*cthetat*Mst1^2*nt*sthetat^7*Tf)/(ep^3*(Mst1^2 - Mst2^2)) + 
   (CF^2*cthetat*Mst1^2*nt*sthetat^7*Tf)/(ep^2*(Mst1^2 - Mst2^2)) + (CF^2*cthetat*Mst2^2*nt*sthetat^7*Tf)/(ep^3*(Mst1^2 - Mst2^2)) - 
   (CF^2*cthetat*Mst2^2*nt*sthetat^7*Tf)/(ep^2*(Mst1^2 - Mst2^2)) + (CF*Mfg*Mt*nq^2*Tf^2)/(4*ep^3*(Mst1^2 - Mst2^2)) - 
   (CF*Mfg*Mt*nq^2*Tf^2)/(4*ep^2*(Mst1^2 - Mst2^2)) - (CF*Mfg*Mt*nq^2*Tf^2)/(4*ep*(Mst1^2 - Mst2^2)) + 
   (CF*Mfg*Mt*nq*nt*Tf^2)/(2*ep^3*(Mst1^2 - Mst2^2)) - (CF*Mfg*Mt*nq*nt*Tf^2)/(2*ep^2*(Mst1^2 - Mst2^2)) - 
   (CF*Mfg*Mt*nq*nt*Tf^2)/(2*ep*(Mst1^2 - Mst2^2)) + (CF*Mfg*Mt*nt^2*Tf^2)/(4*ep^3*(Mst1^2 - Mst2^2)) - 
   (CF*Mfg*Mt*nt^2*Tf^2)/(4*ep^2*(Mst1^2 - Mst2^2)) - (CF*Mfg*Mt*nt^2*Tf^2)/(4*ep*(Mst1^2 - Mst2^2)) - 
   (CF*cthetat*Mst1^2*nq^2*sthetat*Tf^2)/(24*ep^3*(Mst1^2 - Mst2^2)) + (CF*cthetat*Mst1^2*nq^2*sthetat*Tf^2)/(24*ep^2*(Mst1^2 - Mst2^2)) + 
   (CF*cthetat*Mst1^2*nq^2*sthetat*Tf^2)/(24*ep*(Mst1^2 - Mst2^2)) + (CF*cthetat*Mst2^2*nq^2*sthetat*Tf^2)/(24*ep^3*(Mst1^2 - Mst2^2)) - 
   (CF*cthetat*Mst2^2*nq^2*sthetat*Tf^2)/(24*ep^2*(Mst1^2 - Mst2^2)) - (CF*cthetat*Mst2^2*nq^2*sthetat*Tf^2)/(24*ep*(Mst1^2 - Mst2^2)) - 
   (CF*cthetat*Mst1^2*nq*nt*sthetat*Tf^2)/(12*ep^3*(Mst1^2 - Mst2^2)) + (CF*cthetat*Mst1^2*nq*nt*sthetat*Tf^2)/(12*ep^2*(Mst1^2 - Mst2^2)) + 
   (CF*cthetat*Mst1^2*nq*nt*sthetat*Tf^2)/(12*ep*(Mst1^2 - Mst2^2)) + (CF*cthetat*Mst2^2*nq*nt*sthetat*Tf^2)/(12*ep^3*(Mst1^2 - Mst2^2)) - 
   (CF*cthetat*Mst2^2*nq*nt*sthetat*Tf^2)/(12*ep^2*(Mst1^2 - Mst2^2)) - (CF*cthetat*Mst2^2*nq*nt*sthetat*Tf^2)/(12*ep*(Mst1^2 - Mst2^2)) - 
   (CF*cthetat*Mst1^2*nt^2*sthetat*Tf^2)/(24*ep^3*(Mst1^2 - Mst2^2)) + (CF*cthetat*Mst1^2*nt^2*sthetat*Tf^2)/(24*ep^2*(Mst1^2 - Mst2^2)) + 
   (CF*cthetat*Mst1^2*nt^2*sthetat*Tf^2)/(24*ep*(Mst1^2 - Mst2^2)) + (CF*cthetat*Mst2^2*nt^2*sthetat*Tf^2)/(24*ep^3*(Mst1^2 - Mst2^2)) - 
   (CF*cthetat*Mst2^2*nt^2*sthetat*Tf^2)/(24*ep^2*(Mst1^2 - Mst2^2)) - (CF*cthetat*Mst2^2*nt^2*sthetat*Tf^2)/(24*ep*(Mst1^2 - Mst2^2)) - 
   (CF*Mfg*Mt*nq^2*sthetat^2*Tf^2)/(2*ep^3*(Mst1^2 - Mst2^2)) + (CF*Mfg*Mt*nq^2*sthetat^2*Tf^2)/(2*ep^2*(Mst1^2 - Mst2^2)) + 
   (CF*Mfg*Mt*nq^2*sthetat^2*Tf^2)/(2*ep*(Mst1^2 - Mst2^2)) - (CF*Mfg*Mt*nq*nt*sthetat^2*Tf^2)/(ep^3*(Mst1^2 - Mst2^2)) + 
   (CF*Mfg*Mt*nq*nt*sthetat^2*Tf^2)/(ep^2*(Mst1^2 - Mst2^2)) + (CF*Mfg*Mt*nq*nt*sthetat^2*Tf^2)/(ep*(Mst1^2 - Mst2^2)) - 
   (CF*Mfg*Mt*nt^2*sthetat^2*Tf^2)/(2*ep^3*(Mst1^2 - Mst2^2)) + (CF*Mfg*Mt*nt^2*sthetat^2*Tf^2)/(2*ep^2*(Mst1^2 - Mst2^2)) + 
   (CF*Mfg*Mt*nt^2*sthetat^2*Tf^2)/(2*ep*(Mst1^2 - Mst2^2)) + (CF*cthetat*Mst1^2*nq^2*sthetat^3*Tf^2)/(12*ep^3*(Mst1^2 - Mst2^2)) - 
   (CF*cthetat*Mst1^2*nq^2*sthetat^3*Tf^2)/(12*ep^2*(Mst1^2 - Mst2^2)) - (CF*cthetat*Mst1^2*nq^2*sthetat^3*Tf^2)/(12*ep*(Mst1^2 - Mst2^2)) - 
   (CF*cthetat*Mst2^2*nq^2*sthetat^3*Tf^2)/(12*ep^3*(Mst1^2 - Mst2^2)) + (CF*cthetat*Mst2^2*nq^2*sthetat^3*Tf^2)/
    (12*ep^2*(Mst1^2 - Mst2^2)) + (CF*cthetat*Mst2^2*nq^2*sthetat^3*Tf^2)/(12*ep*(Mst1^2 - Mst2^2)) + 
   (CF*cthetat*Mst1^2*nq*nt*sthetat^3*Tf^2)/(6*ep^3*(Mst1^2 - Mst2^2)) - (CF*cthetat*Mst1^2*nq*nt*sthetat^3*Tf^2)/
    (6*ep^2*(Mst1^2 - Mst2^2)) - (CF*cthetat*Mst1^2*nq*nt*sthetat^3*Tf^2)/(6*ep*(Mst1^2 - Mst2^2)) - 
   (CF*cthetat*Mst2^2*nq*nt*sthetat^3*Tf^2)/(6*ep^3*(Mst1^2 - Mst2^2)) + (CF*cthetat*Mst2^2*nq*nt*sthetat^3*Tf^2)/
    (6*ep^2*(Mst1^2 - Mst2^2)) + (CF*cthetat*Mst2^2*nq*nt*sthetat^3*Tf^2)/(6*ep*(Mst1^2 - Mst2^2)) + 
   (CF*cthetat*Mst1^2*nt^2*sthetat^3*Tf^2)/(12*ep^3*(Mst1^2 - Mst2^2)) - (CF*cthetat*Mst1^2*nt^2*sthetat^3*Tf^2)/
    (12*ep^2*(Mst1^2 - Mst2^2)) - (CF*cthetat*Mst1^2*nt^2*sthetat^3*Tf^2)/(12*ep*(Mst1^2 - Mst2^2)) - 
   (CF*cthetat*Mst2^2*nt^2*sthetat^3*Tf^2)/(12*ep^3*(Mst1^2 - Mst2^2)) + (CF*cthetat*Mst2^2*nt^2*sthetat^3*Tf^2)/
    (12*ep^2*(Mst1^2 - Mst2^2)) + (CF*cthetat*Mst2^2*nt^2*sthetat^3*Tf^2)/(12*ep*(Mst1^2 - Mst2^2)) - 
   (3*CA*CF*Mfg*Mt*nq*Tf*z3)/(2*ep*(Mst1^2 - Mst2^2)) + (3*CF^2*Mfg*Mt*nq*Tf*z3)/(2*ep*(Mst1^2 - Mst2^2)) - 
   (3*CA*CF*Mfg*Mt*nt*Tf*z3)/(2*ep*(Mst1^2 - Mst2^2)) + (3*CF^2*Mfg*Mt*nt*Tf*z3)/(2*ep*(Mst1^2 - Mst2^2)) + 
   (CA*CF*cthetat*Mst1^2*nq*sthetat*Tf*z3)/(4*ep*(Mst1^2 - Mst2^2)) - (CF^2*cthetat*Mst1^2*nq*sthetat*Tf*z3)/(4*ep*(Mst1^2 - Mst2^2)) - 
   (CA*CF*cthetat*Mst2^2*nq*sthetat*Tf*z3)/(4*ep*(Mst1^2 - Mst2^2)) + (CF^2*cthetat*Mst2^2*nq*sthetat*Tf*z3)/(4*ep*(Mst1^2 - Mst2^2)) + 
   (CA*CF*cthetat*Mst1^2*nt*sthetat*Tf*z3)/(4*ep*(Mst1^2 - Mst2^2)) - (CF^2*cthetat*Mst1^2*nt*sthetat*Tf*z3)/(4*ep*(Mst1^2 - Mst2^2)) - 
   (CA*CF*cthetat*Mst2^2*nt*sthetat*Tf*z3)/(4*ep*(Mst1^2 - Mst2^2)) + (CF^2*cthetat*Mst2^2*nt*sthetat*Tf*z3)/(4*ep*(Mst1^2 - Mst2^2)) + 
   (3*CA*CF*Mfg*Mt*nq*sthetat^2*Tf*z3)/(ep*(Mst1^2 - Mst2^2)) - (3*CF^2*Mfg*Mt*nq*sthetat^2*Tf*z3)/(ep*(Mst1^2 - Mst2^2)) + 
   (3*CA*CF*Mfg*Mt*nt*sthetat^2*Tf*z3)/(ep*(Mst1^2 - Mst2^2)) - (3*CF^2*Mfg*Mt*nt*sthetat^2*Tf*z3)/(ep*(Mst1^2 - Mst2^2)) - 
   (CA*CF*cthetat*Mst1^2*nq*sthetat^3*Tf*z3)/(2*ep*(Mst1^2 - Mst2^2)) + (CF^2*cthetat*Mst1^2*nq*sthetat^3*Tf*z3)/(2*ep*(Mst1^2 - Mst2^2)) + 
   (CA*CF*cthetat*Mst2^2*nq*sthetat^3*Tf*z3)/(2*ep*(Mst1^2 - Mst2^2)) - (CF^2*cthetat*Mst2^2*nq*sthetat^3*Tf*z3)/(2*ep*(Mst1^2 - Mst2^2)) - 
   (CA*CF*cthetat*Mst1^2*nt*sthetat^3*Tf*z3)/(2*ep*(Mst1^2 - Mst2^2)) + (CF^2*cthetat*Mst1^2*nt*sthetat^3*Tf*z3)/(2*ep*(Mst1^2 - Mst2^2)) + 
   (CA*CF*cthetat*Mst2^2*nt*sthetat^3*Tf*z3)/(2*ep*(Mst1^2 - Mst2^2)) - (CF^2*cthetat*Mst2^2*nt*sthetat^3*Tf*z3)/(2*ep*(Mst1^2 - Mst2^2))));


Zmsq = (1 -((as*CF*Mfg^2)/(ep*Msq^2)) + as^2*((3*CF*Mfg^2*(3*CA - 2*(nq + nt)*Tf))/(8*ep^2*Msq^2) + 
   (CF*(-(CA*(3*Me^2 + 11*Mfg^2)) + 2*(3*CF*Mfg^2 + (2*Msq^2*nq + Mst1^2*nt + Mst2^2*nt - 2*Mt^2*nt + Me^2*(nq + nt) + 3*Mfg^2*(nq + nt))*
         Tf)))/(8*ep*Msq^2)) + as^3*(-(CF*Mfg^2*(3*CA - 2*(nq + nt)*Tf)^2)/(8*ep^3*Msq^2) + 
   (CF*(CA^2*(9*Me^2 + 66*Mfg^2) + 4*(nq + nt)*Tf*(-4*CF*Mfg^2 + (2*Msq^2*nq + Mst1^2*nt + Mst2^2*nt - 2*Mt^2*nt + Me^2*(nq + nt) + 
          4*Mfg^2*(nq + nt))*Tf) - 2*CA*(12*CF*Mfg^2 + (6*Me^2*(nq + nt) + 34*Mfg^2*(nq + nt) + 
          3*(2*Msq^2*nq + (Mst1^2 + Mst2^2 - 2*Mt^2)*nt))*Tf)))/(32*ep^2*Msq^2) - 
   (CF*(5*CA^2*(3*Me^2 + 8*Mfg^2) - 2*CA*(CF*(3*Me^2 + 14*Mfg^2) + Tf*(10*Msq^2*nq + 5*Mst1^2*nt + 5*Mst2^2*nt - 10*Mt^2*nt + 
          8*Me^2*(nq + nt) + 48*Mfg^2*nq*z3 + 48*Mfg^2*nt*z3)) + 
      4*(8*CF^2*Mfg^2 + (nq + nt)*(2*Msq^2*nq + Mst1^2*nt + Mst2^2*nt - 2*Mt^2*nt + Me^2*(nq + nt) - 4*Mfg^2*(nq + nt))*Tf^2 + 
        CF*Tf*(2*Msq^2*nq + Mst1^2*nt + Mst2^2*nt - 2*Mt^2*nt + Me^2*(nq + nt) + 6*Mfg^2*(nq + nt)*(-3 + 4*z3)))))/(32*ep*Msq^2)));


gamThetat = ((as*(2*CF*cthetat^2*Mfg*Mt - CF*cthetat^3*Mst1^2*sthetat + CF*cthetat^3*Mst2^2*sthetat - 2*CF*Mfg*Mt*sthetat^2 + 
    CF*cthetat*Mst1^2*sthetat^3 - CF*cthetat*Mst2^2*sthetat^3))/(2*(Mst1^2 - Mst2^2)*Thetat)
 + (as^2*(12*CA*CF*cthetat^2*Mfg*Mt - 8*CF^2*cthetat^2*Mfg*Mt - 8*CF*cthetat^2*Tf*Mfg*Mt*nq - 8*CF*cthetat^2*Tf*Mfg*Mt*nt - 
    3*CA*CF*cthetat^3*Mst1^2*sthetat + 2*CF^2*cthetat^3*Mst1^2*sthetat + 3*CA*CF*cthetat^3*Mst2^2*sthetat - 
    2*CF^2*cthetat^3*Mst2^2*sthetat + 2*CF*cthetat^3*Tf*Mst1^2*nq*sthetat - 2*CF*cthetat^3*Tf*Mst2^2*nq*sthetat + 
    2*CF*cthetat^3*Tf*Mst1^2*nt*sthetat - 2*CF*cthetat^3*Tf*Mst2^2*nt*sthetat - 12*CA*CF*Mfg*Mt*sthetat^2 + 8*CF^2*Mfg*Mt*sthetat^2 + 
    8*CF*Tf*Mfg*Mt*nq*sthetat^2 + 8*CF*Tf*Mfg*Mt*nt*sthetat^2 + 3*CA*CF*cthetat*Mst1^2*sthetat^3 - 2*CF^2*cthetat*Mst1^2*sthetat^3 - 
    3*CA*CF*cthetat*Mst2^2*sthetat^3 + 2*CF^2*cthetat*Mst2^2*sthetat^3 - 2*CF*cthetat*Tf*Mst1^2*nq*sthetat^3 + 
    2*CF*cthetat*Tf*Mst2^2*nq*sthetat^3 - 2*CF*cthetat*Tf*Mst1^2*nt*sthetat^3 + 2*CF*cthetat*Tf*Mst2^2*nt*sthetat^3))/
  (8*(Mst1^2 - Mst2^2)*Thetat)
 + (as^3*(18*CA^2*CF*cthetat^2*Mfg*Mt - 18*CA*CF^2*cthetat^2*Mfg*Mt + 24*CF^3*cthetat^2*Mfg*Mt + 
    6*CA*CF*cthetat^2*Tf*Mfg*Mt*nq - 48*CF^2*cthetat^2*Tf*Mfg*Mt*nq - 12*CF*cthetat^2*Tf^2*Mfg*Mt*nq^2 + 
    6*CA*CF*cthetat^2*Tf*Mfg*Mt*nt - 48*CF^2*cthetat^2*Tf*Mfg*Mt*nt - 24*CF*cthetat^2*Tf^2*Mfg*Mt*nq*nt - 
    12*CF*cthetat^2*Tf^2*Mfg*Mt*nt^2 - 3*CA^2*CF*cthetat^3*Mst1^2*sthetat + 3*CA*CF^2*cthetat^3*Mst1^2*sthetat - 
    4*CF^3*cthetat^3*Mst1^2*sthetat + 3*CA^2*CF*cthetat^3*Mst2^2*sthetat - 3*CA*CF^2*cthetat^3*Mst2^2*sthetat + 
    4*CF^3*cthetat^3*Mst2^2*sthetat - CA*CF*cthetat^3*Tf*Mst1^2*nq*sthetat + 8*CF^2*cthetat^3*Tf*Mst1^2*nq*sthetat + 
    CA*CF*cthetat^3*Tf*Mst2^2*nq*sthetat - 8*CF^2*cthetat^3*Tf*Mst2^2*nq*sthetat + 2*CF*cthetat^3*Tf^2*Mst1^2*nq^2*sthetat - 
    2*CF*cthetat^3*Tf^2*Mst2^2*nq^2*sthetat - CA*CF*cthetat^3*Tf*Mst1^2*nt*sthetat + 8*CF^2*cthetat^3*Tf*Mst1^2*nt*sthetat + 
    CA*CF*cthetat^3*Tf*Mst2^2*nt*sthetat - 8*CF^2*cthetat^3*Tf*Mst2^2*nt*sthetat + 4*CF*cthetat^3*Tf^2*Mst1^2*nq*nt*sthetat - 
    4*CF*cthetat^3*Tf^2*Mst2^2*nq*nt*sthetat + 2*CF*cthetat^3*Tf^2*Mst1^2*nt^2*sthetat - 2*CF*cthetat^3*Tf^2*Mst2^2*nt^2*sthetat - 
    18*CA^2*CF*Mfg*Mt*sthetat^2 + 18*CA*CF^2*Mfg*Mt*sthetat^2 - 24*CF^3*Mfg*Mt*sthetat^2 - 6*CA*CF*Tf*Mfg*Mt*nq*sthetat^2 + 
    48*CF^2*Tf*Mfg*Mt*nq*sthetat^2 + 12*CF*Tf^2*Mfg*Mt*nq^2*sthetat^2 - 6*CA*CF*Tf*Mfg*Mt*nt*sthetat^2 + 
    48*CF^2*Tf*Mfg*Mt*nt*sthetat^2 + 24*CF*Tf^2*Mfg*Mt*nq*nt*sthetat^2 + 12*CF*Tf^2*Mfg*Mt*nt^2*sthetat^2 + 
    3*CA^2*CF*cthetat*Mst1^2*sthetat^3 - 3*CA*CF^2*cthetat*Mst1^2*sthetat^3 + 4*CF^3*cthetat*Mst1^2*sthetat^3 - 
    3*CA^2*CF*cthetat*Mst2^2*sthetat^3 + 3*CA*CF^2*cthetat*Mst2^2*sthetat^3 - 4*CF^3*cthetat*Mst2^2*sthetat^3 + 
    CA*CF*cthetat*Tf*Mst1^2*nq*sthetat^3 - 8*CF^2*cthetat*Tf*Mst1^2*nq*sthetat^3 - CA*CF*cthetat*Tf*Mst2^2*nq*sthetat^3 + 
    8*CF^2*cthetat*Tf*Mst2^2*nq*sthetat^3 - 2*CF*cthetat*Tf^2*Mst1^2*nq^2*sthetat^3 + 2*CF*cthetat*Tf^2*Mst2^2*nq^2*sthetat^3 + 
    CA*CF*cthetat*Tf*Mst1^2*nt*sthetat^3 - 8*CF^2*cthetat*Tf*Mst1^2*nt*sthetat^3 - CA*CF*cthetat*Tf*Mst2^2*nt*sthetat^3 + 
    8*CF^2*cthetat*Tf*Mst2^2*nt*sthetat^3 - 4*CF*cthetat*Tf^2*Mst1^2*nq*nt*sthetat^3 + 4*CF*cthetat*Tf^2*Mst2^2*nq*nt*sthetat^3 - 
    2*CF*cthetat*Tf^2*Mst1^2*nt^2*sthetat^3 + 2*CF*cthetat*Tf^2*Mst2^2*nt^2*sthetat^3 - 72*CA*CF*cthetat^2*Tf*Mfg*Mt*nq*z3 + 
    72*CF^2*cthetat^2*Tf*Mfg*Mt*nq*z3 - 72*CA*CF*cthetat^2*Tf*Mfg*Mt*nt*z3 + 72*CF^2*cthetat^2*Tf*Mfg*Mt*nt*z3 + 
    12*CA*CF*cthetat^3*Tf*Mst1^2*nq*sthetat*z3 - 12*CF^2*cthetat^3*Tf*Mst1^2*nq*sthetat*z3 - 12*CA*CF*cthetat^3*Tf*Mst2^2*nq*sthetat*z3 + 
    12*CF^2*cthetat^3*Tf*Mst2^2*nq*sthetat*z3 + 12*CA*CF*cthetat^3*Tf*Mst1^2*nt*sthetat*z3 - 12*CF^2*cthetat^3*Tf*Mst1^2*nt*sthetat*z3 - 
    12*CA*CF*cthetat^3*Tf*Mst2^2*nt*sthetat*z3 + 12*CF^2*cthetat^3*Tf*Mst2^2*nt*sthetat*z3 + 72*CA*CF*Tf*Mfg*Mt*nq*sthetat^2*z3 - 
    72*CF^2*Tf*Mfg*Mt*nq*sthetat^2*z3 + 72*CA*CF*Tf*Mfg*Mt*nt*sthetat^2*z3 - 72*CF^2*Tf*Mfg*Mt*nt*sthetat^2*z3 - 
    12*CA*CF*cthetat*Tf*Mst1^2*nq*sthetat^3*z3 + 12*CF^2*cthetat*Tf*Mst1^2*nq*sthetat^3*z3 + 12*CA*CF*cthetat*Tf*Mst2^2*nq*sthetat^3*z3 - 
    12*CF^2*cthetat*Tf*Mst2^2*nq*sthetat^3*z3 - 12*CA*CF*cthetat*Tf*Mst1^2*nt*sthetat^3*z3 + 12*CF^2*cthetat*Tf*Mst1^2*nt*sthetat^3*z3 + 
    12*CA*CF*cthetat*Tf*Mst2^2*nt*sthetat^3*z3 - 12*CF^2*cthetat*Tf*Mst2^2*nt*sthetat^3*z3))/(16*(Mst1^2 - Mst2^2)*Thetat));


gamMst1 = ((as*(-8*CF*Mfg^2 - CF*Mst1^2 + CF*cthetat^4*Mst1^2 + CF*Mst2^2 - CF*cthetat^4*Mst2^2 - 8*CF*Mt^2 + 16*CF*cthetat*Mfg*Mt*sthetat - 
    6*CF*cthetat^2*Mst1^2*sthetat^2 + 6*CF*cthetat^2*Mst2^2*sthetat^2 + CF*Mst1^2*sthetat^4 - CF*Mst2^2*sthetat^4))/(8*Mst1^2) + 
 (as^2*(-24*CA*CF*Me^2 - 88*CA*CF*Mfg^2 + 48*CF^2*Mfg^2 - 3*CA*CF*Mst1^2 + 2*CF^2*Mst1^2 + 3*CA*CF*cthetat^4*Mst1^2 - 
    2*CF^2*cthetat^4*Mst1^2 + 3*CA*CF*Mst2^2 - 2*CF^2*Mst2^2 - 3*CA*CF*cthetat^4*Mst2^2 + 2*CF^2*cthetat^4*Mst2^2 - 24*CA*CF*Mt^2 + 
    16*CF^2*Mt^2 + 96*CA*CF*cthetat*Mfg*Mt*sthetat - 64*CF^2*cthetat*Mfg*Mt*sthetat - 18*CA*CF*cthetat^2*Mst1^2*sthetat^2 + 
    12*CF^2*cthetat^2*Mst1^2*sthetat^2 + 18*CA*CF*cthetat^2*Mst2^2*sthetat^2 - 12*CF^2*cthetat^2*Mst2^2*sthetat^2 + 
    3*CA*CF*Mst1^2*sthetat^4 - 2*CF^2*Mst1^2*sthetat^4 - 3*CA*CF*Mst2^2*sthetat^4 + 2*CF^2*Mst2^2*sthetat^4 + 16*CF*Me^2*nq*Tf + 
    48*CF*Mfg^2*nq*Tf + 32*CF*Msq^2*nq*Tf + 2*CF*Mst1^2*nq*Tf - 2*CF*cthetat^4*Mst1^2*nq*Tf - 2*CF*Mst2^2*nq*Tf + 
    2*CF*cthetat^4*Mst2^2*nq*Tf + 16*CF*Mt^2*nq*Tf + 16*CF*Me^2*nt*Tf + 48*CF*Mfg^2*nt*Tf + 18*CF*Mst1^2*nt*Tf - 
    2*CF*cthetat^4*Mst1^2*nt*Tf + 14*CF*Mst2^2*nt*Tf + 2*CF*cthetat^4*Mst2^2*nt*Tf - 16*CF*Mt^2*nt*Tf - 64*CF*cthetat*Mfg*Mt*nq*sthetat*Tf - 
    64*CF*cthetat*Mfg*Mt*nt*sthetat*Tf + 12*CF*cthetat^2*Mst1^2*nq*sthetat^2*Tf - 12*CF*cthetat^2*Mst2^2*nq*sthetat^2*Tf + 
    12*CF*cthetat^2*Mst1^2*nt*sthetat^2*Tf - 12*CF*cthetat^2*Mst2^2*nt*sthetat^2*Tf - 2*CF*Mst1^2*nq*sthetat^4*Tf + 
    2*CF*Mst2^2*nq*sthetat^4*Tf - 2*CF*Mst1^2*nt*sthetat^4*Tf + 2*CF*Mst2^2*nt*sthetat^4*Tf))/(32*Mst1^2) + 
 (as^3*(-90*CA^2*CF*Me^2 + 36*CA*CF^2*Me^2 - 240*CA^2*CF*Mfg^2 + 168*CA*CF^2*Mfg^2 - 192*CF^3*Mfg^2 - 3*CA^2*CF*Mst1^2 + 3*CA*CF^2*Mst1^2 - 
    4*CF^3*Mst1^2 + 3*CA^2*CF*cthetat^4*Mst1^2 - 3*CA*CF^2*cthetat^4*Mst1^2 + 4*CF^3*cthetat^4*Mst1^2 + 3*CA^2*CF*Mst2^2 - 
    3*CA*CF^2*Mst2^2 + 4*CF^3*Mst2^2 - 3*CA^2*CF*cthetat^4*Mst2^2 + 3*CA*CF^2*cthetat^4*Mst2^2 - 4*CF^3*cthetat^4*Mst2^2 - 24*CA^2*CF*Mt^2 + 
    24*CA*CF^2*Mt^2 - 32*CF^3*Mt^2 + 144*CA^2*CF*cthetat*Mfg*Mt*sthetat - 144*CA*CF^2*cthetat*Mfg*Mt*sthetat + 
    192*CF^3*cthetat*Mfg*Mt*sthetat - 18*CA^2*CF*cthetat^2*Mst1^2*sthetat^2 + 18*CA*CF^2*cthetat^2*Mst1^2*sthetat^2 - 
    24*CF^3*cthetat^2*Mst1^2*sthetat^2 + 18*CA^2*CF*cthetat^2*Mst2^2*sthetat^2 - 18*CA*CF^2*cthetat^2*Mst2^2*sthetat^2 + 
    24*CF^3*cthetat^2*Mst2^2*sthetat^2 + 3*CA^2*CF*Mst1^2*sthetat^4 - 3*CA*CF^2*Mst1^2*sthetat^4 + 4*CF^3*Mst1^2*sthetat^4 - 
    3*CA^2*CF*Mst2^2*sthetat^4 + 3*CA*CF^2*Mst2^2*sthetat^4 - 4*CF^3*Mst2^2*sthetat^4 + 96*CA*CF*Me^2*nq*Tf - 24*CF^2*Me^2*nq*Tf + 
    432*CF^2*Mfg^2*nq*Tf + 120*CA*CF*Msq^2*nq*Tf - 48*CF^2*Msq^2*nq*Tf - CA*CF*Mst1^2*nq*Tf + 8*CF^2*Mst1^2*nq*Tf + 
    CA*CF*cthetat^4*Mst1^2*nq*Tf - 8*CF^2*cthetat^4*Mst1^2*nq*Tf + CA*CF*Mst2^2*nq*Tf - 8*CF^2*Mst2^2*nq*Tf - CA*CF*cthetat^4*Mst2^2*nq*Tf + 
    8*CF^2*cthetat^4*Mst2^2*nq*Tf - 8*CA*CF*Mt^2*nq*Tf + 64*CF^2*Mt^2*nq*Tf + 96*CA*CF*Me^2*nt*Tf - 24*CF^2*Me^2*nt*Tf + 
    432*CF^2*Mfg^2*nt*Tf + 59*CA*CF*Mst1^2*nt*Tf - 16*CF^2*Mst1^2*nt*Tf + CA*CF*cthetat^4*Mst1^2*nt*Tf - 8*CF^2*cthetat^4*Mst1^2*nt*Tf + 
    61*CA*CF*Mst2^2*nt*Tf - 32*CF^2*Mst2^2*nt*Tf - CA*CF*cthetat^4*Mst2^2*nt*Tf + 8*CF^2*cthetat^4*Mst2^2*nt*Tf - 128*CA*CF*Mt^2*nt*Tf + 
    112*CF^2*Mt^2*nt*Tf + 48*CA*CF*cthetat*Mfg*Mt*nq*sthetat*Tf - 384*CF^2*cthetat*Mfg*Mt*nq*sthetat*Tf + 
    48*CA*CF*cthetat*Mfg*Mt*nt*sthetat*Tf - 384*CF^2*cthetat*Mfg*Mt*nt*sthetat*Tf - 6*CA*CF*cthetat^2*Mst1^2*nq*sthetat^2*Tf + 
    48*CF^2*cthetat^2*Mst1^2*nq*sthetat^2*Tf + 6*CA*CF*cthetat^2*Mst2^2*nq*sthetat^2*Tf - 48*CF^2*cthetat^2*Mst2^2*nq*sthetat^2*Tf - 
    6*CA*CF*cthetat^2*Mst1^2*nt*sthetat^2*Tf + 48*CF^2*cthetat^2*Mst1^2*nt*sthetat^2*Tf + 6*CA*CF*cthetat^2*Mst2^2*nt*sthetat^2*Tf - 
    48*CF^2*cthetat^2*Mst2^2*nt*sthetat^2*Tf + CA*CF*Mst1^2*nq*sthetat^4*Tf - 8*CF^2*Mst1^2*nq*sthetat^4*Tf - CA*CF*Mst2^2*nq*sthetat^4*Tf + 
    8*CF^2*Mst2^2*nq*sthetat^4*Tf + CA*CF*Mst1^2*nt*sthetat^4*Tf - 8*CF^2*Mst1^2*nt*sthetat^4*Tf - CA*CF*Mst2^2*nt*sthetat^4*Tf + 
    8*CF^2*Mst2^2*nt*sthetat^4*Tf - 24*CF*Me^2*nq^2*Tf^2 + 96*CF*Mfg^2*nq^2*Tf^2 - 48*CF*Msq^2*nq^2*Tf^2 + 2*CF*Mst1^2*nq^2*Tf^2 - 
    2*CF*cthetat^4*Mst1^2*nq^2*Tf^2 - 2*CF*Mst2^2*nq^2*Tf^2 + 2*CF*cthetat^4*Mst2^2*nq^2*Tf^2 + 16*CF*Mt^2*nq^2*Tf^2 - 
    48*CF*Me^2*nq*nt*Tf^2 + 192*CF*Mfg^2*nq*nt*Tf^2 - 48*CF*Msq^2*nq*nt*Tf^2 - 20*CF*Mst1^2*nq*nt*Tf^2 - 4*CF*cthetat^4*Mst1^2*nq*nt*Tf^2 - 
    28*CF*Mst2^2*nq*nt*Tf^2 + 4*CF*cthetat^4*Mst2^2*nq*nt*Tf^2 + 80*CF*Mt^2*nq*nt*Tf^2 - 24*CF*Me^2*nt^2*Tf^2 + 96*CF*Mfg^2*nt^2*Tf^2 - 
    22*CF*Mst1^2*nt^2*Tf^2 - 2*CF*cthetat^4*Mst1^2*nt^2*Tf^2 - 26*CF*Mst2^2*nt^2*Tf^2 + 2*CF*cthetat^4*Mst2^2*nt^2*Tf^2 + 
    64*CF*Mt^2*nt^2*Tf^2 - 96*CF*cthetat*Mfg*Mt*nq^2*sthetat*Tf^2 - 192*CF*cthetat*Mfg*Mt*nq*nt*sthetat*Tf^2 - 
    96*CF*cthetat*Mfg*Mt*nt^2*sthetat*Tf^2 + 12*CF*cthetat^2*Mst1^2*nq^2*sthetat^2*Tf^2 - 12*CF*cthetat^2*Mst2^2*nq^2*sthetat^2*Tf^2 + 
    24*CF*cthetat^2*Mst1^2*nq*nt*sthetat^2*Tf^2 - 24*CF*cthetat^2*Mst2^2*nq*nt*sthetat^2*Tf^2 + 12*CF*cthetat^2*Mst1^2*nt^2*sthetat^2*Tf^2 - 
    12*CF*cthetat^2*Mst2^2*nt^2*sthetat^2*Tf^2 - 2*CF*Mst1^2*nq^2*sthetat^4*Tf^2 + 2*CF*Mst2^2*nq^2*sthetat^4*Tf^2 - 
    4*CF*Mst1^2*nq*nt*sthetat^4*Tf^2 + 4*CF*Mst2^2*nq*nt*sthetat^4*Tf^2 - 2*CF*Mst1^2*nt^2*sthetat^4*Tf^2 + 
    2*CF*Mst2^2*nt^2*sthetat^4*Tf^2 + 576*CA*CF*Mfg^2*nq*Tf*z3 - 576*CF^2*Mfg^2*nq*Tf*z3 + 12*CA*CF*Mst1^2*nq*Tf*z3 - 
    12*CF^2*Mst1^2*nq*Tf*z3 - 12*CA*CF*cthetat^4*Mst1^2*nq*Tf*z3 + 12*CF^2*cthetat^4*Mst1^2*nq*Tf*z3 - 12*CA*CF*Mst2^2*nq*Tf*z3 + 
    12*CF^2*Mst2^2*nq*Tf*z3 + 12*CA*CF*cthetat^4*Mst2^2*nq*Tf*z3 - 12*CF^2*cthetat^4*Mst2^2*nq*Tf*z3 + 96*CA*CF*Mt^2*nq*Tf*z3 - 
    96*CF^2*Mt^2*nq*Tf*z3 + 576*CA*CF*Mfg^2*nt*Tf*z3 - 576*CF^2*Mfg^2*nt*Tf*z3 + 12*CA*CF*Mst1^2*nt*Tf*z3 - 12*CF^2*Mst1^2*nt*Tf*z3 - 
    12*CA*CF*cthetat^4*Mst1^2*nt*Tf*z3 + 12*CF^2*cthetat^4*Mst1^2*nt*Tf*z3 - 12*CA*CF*Mst2^2*nt*Tf*z3 + 12*CF^2*Mst2^2*nt*Tf*z3 + 
    12*CA*CF*cthetat^4*Mst2^2*nt*Tf*z3 - 12*CF^2*cthetat^4*Mst2^2*nt*Tf*z3 + 96*CA*CF*Mt^2*nt*Tf*z3 - 96*CF^2*Mt^2*nt*Tf*z3 - 
    576*CA*CF*cthetat*Mfg*Mt*nq*sthetat*Tf*z3 + 576*CF^2*cthetat*Mfg*Mt*nq*sthetat*Tf*z3 - 576*CA*CF*cthetat*Mfg*Mt*nt*sthetat*Tf*z3 + 
    576*CF^2*cthetat*Mfg*Mt*nt*sthetat*Tf*z3 + 72*CA*CF*cthetat^2*Mst1^2*nq*sthetat^2*Tf*z3 - 72*CF^2*cthetat^2*Mst1^2*nq*sthetat^2*Tf*z3 - 
    72*CA*CF*cthetat^2*Mst2^2*nq*sthetat^2*Tf*z3 + 72*CF^2*cthetat^2*Mst2^2*nq*sthetat^2*Tf*z3 + 
    72*CA*CF*cthetat^2*Mst1^2*nt*sthetat^2*Tf*z3 - 72*CF^2*cthetat^2*Mst1^2*nt*sthetat^2*Tf*z3 - 
    72*CA*CF*cthetat^2*Mst2^2*nt*sthetat^2*Tf*z3 + 72*CF^2*cthetat^2*Mst2^2*nt*sthetat^2*Tf*z3 - 12*CA*CF*Mst1^2*nq*sthetat^4*Tf*z3 + 
    12*CF^2*Mst1^2*nq*sthetat^4*Tf*z3 + 12*CA*CF*Mst2^2*nq*sthetat^4*Tf*z3 - 12*CF^2*Mst2^2*nq*sthetat^4*Tf*z3 - 
    12*CA*CF*Mst1^2*nt*sthetat^4*Tf*z3 + 12*CF^2*Mst1^2*nt*sthetat^4*Tf*z3 + 12*CA*CF*Mst2^2*nt*sthetat^4*Tf*z3 - 
    12*CF^2*Mst2^2*nt*sthetat^4*Tf*z3))/(64*Mst1^2));


gamMst2 = ((as*(-8*CF*Mfg^2 + CF*Mst1^2 - CF*cthetat^4*Mst1^2 - CF*Mst2^2 + CF*cthetat^4*Mst2^2 - 8*CF*Mt^2 - 16*CF*cthetat*Mfg*Mt*sthetat + 
    6*CF*cthetat^2*Mst1^2*sthetat^2 - 6*CF*cthetat^2*Mst2^2*sthetat^2 - CF*Mst1^2*sthetat^4 + CF*Mst2^2*sthetat^4))/(8*Mst2^2) + 
 (as^2*(-24*CA*CF*Me^2 - 88*CA*CF*Mfg^2 + 48*CF^2*Mfg^2 + 3*CA*CF*Mst1^2 - 2*CF^2*Mst1^2 - 3*CA*CF*cthetat^4*Mst1^2 + 
    2*CF^2*cthetat^4*Mst1^2 - 3*CA*CF*Mst2^2 + 2*CF^2*Mst2^2 + 3*CA*CF*cthetat^4*Mst2^2 - 2*CF^2*cthetat^4*Mst2^2 - 24*CA*CF*Mt^2 + 
    16*CF^2*Mt^2 - 96*CA*CF*cthetat*Mfg*Mt*sthetat + 64*CF^2*cthetat*Mfg*Mt*sthetat + 18*CA*CF*cthetat^2*Mst1^2*sthetat^2 - 
    12*CF^2*cthetat^2*Mst1^2*sthetat^2 - 18*CA*CF*cthetat^2*Mst2^2*sthetat^2 + 12*CF^2*cthetat^2*Mst2^2*sthetat^2 - 
    3*CA*CF*Mst1^2*sthetat^4 + 2*CF^2*Mst1^2*sthetat^4 + 3*CA*CF*Mst2^2*sthetat^4 - 2*CF^2*Mst2^2*sthetat^4 + 16*CF*Me^2*nq*Tf + 
    48*CF*Mfg^2*nq*Tf + 32*CF*Msq^2*nq*Tf - 2*CF*Mst1^2*nq*Tf + 2*CF*cthetat^4*Mst1^2*nq*Tf + 2*CF*Mst2^2*nq*Tf - 
    2*CF*cthetat^4*Mst2^2*nq*Tf + 16*CF*Mt^2*nq*Tf + 16*CF*Me^2*nt*Tf + 48*CF*Mfg^2*nt*Tf + 14*CF*Mst1^2*nt*Tf + 
    2*CF*cthetat^4*Mst1^2*nt*Tf + 18*CF*Mst2^2*nt*Tf - 2*CF*cthetat^4*Mst2^2*nt*Tf - 16*CF*Mt^2*nt*Tf + 64*CF*cthetat*Mfg*Mt*nq*sthetat*Tf + 
    64*CF*cthetat*Mfg*Mt*nt*sthetat*Tf - 12*CF*cthetat^2*Mst1^2*nq*sthetat^2*Tf + 12*CF*cthetat^2*Mst2^2*nq*sthetat^2*Tf - 
    12*CF*cthetat^2*Mst1^2*nt*sthetat^2*Tf + 12*CF*cthetat^2*Mst2^2*nt*sthetat^2*Tf + 2*CF*Mst1^2*nq*sthetat^4*Tf - 
    2*CF*Mst2^2*nq*sthetat^4*Tf + 2*CF*Mst1^2*nt*sthetat^4*Tf - 2*CF*Mst2^2*nt*sthetat^4*Tf))/(32*Mst2^2) + 
 (as^3*(-90*CA^2*CF*Me^2 + 36*CA*CF^2*Me^2 - 240*CA^2*CF*Mfg^2 + 168*CA*CF^2*Mfg^2 - 192*CF^3*Mfg^2 + 3*CA^2*CF*Mst1^2 - 3*CA*CF^2*Mst1^2 + 
    4*CF^3*Mst1^2 - 3*CA^2*CF*cthetat^4*Mst1^2 + 3*CA*CF^2*cthetat^4*Mst1^2 - 4*CF^3*cthetat^4*Mst1^2 - 3*CA^2*CF*Mst2^2 + 
    3*CA*CF^2*Mst2^2 - 4*CF^3*Mst2^2 + 3*CA^2*CF*cthetat^4*Mst2^2 - 3*CA*CF^2*cthetat^4*Mst2^2 + 4*CF^3*cthetat^4*Mst2^2 - 24*CA^2*CF*Mt^2 + 
    24*CA*CF^2*Mt^2 - 32*CF^3*Mt^2 - 144*CA^2*CF*cthetat*Mfg*Mt*sthetat + 144*CA*CF^2*cthetat*Mfg*Mt*sthetat - 
    192*CF^3*cthetat*Mfg*Mt*sthetat + 18*CA^2*CF*cthetat^2*Mst1^2*sthetat^2 - 18*CA*CF^2*cthetat^2*Mst1^2*sthetat^2 + 
    24*CF^3*cthetat^2*Mst1^2*sthetat^2 - 18*CA^2*CF*cthetat^2*Mst2^2*sthetat^2 + 18*CA*CF^2*cthetat^2*Mst2^2*sthetat^2 - 
    24*CF^3*cthetat^2*Mst2^2*sthetat^2 - 3*CA^2*CF*Mst1^2*sthetat^4 + 3*CA*CF^2*Mst1^2*sthetat^4 - 4*CF^3*Mst1^2*sthetat^4 + 
    3*CA^2*CF*Mst2^2*sthetat^4 - 3*CA*CF^2*Mst2^2*sthetat^4 + 4*CF^3*Mst2^2*sthetat^4 + 96*CA*CF*Me^2*nq*Tf - 24*CF^2*Me^2*nq*Tf + 
    432*CF^2*Mfg^2*nq*Tf + 120*CA*CF*Msq^2*nq*Tf - 48*CF^2*Msq^2*nq*Tf + CA*CF*Mst1^2*nq*Tf - 8*CF^2*Mst1^2*nq*Tf - 
    CA*CF*cthetat^4*Mst1^2*nq*Tf + 8*CF^2*cthetat^4*Mst1^2*nq*Tf - CA*CF*Mst2^2*nq*Tf + 8*CF^2*Mst2^2*nq*Tf + CA*CF*cthetat^4*Mst2^2*nq*Tf - 
    8*CF^2*cthetat^4*Mst2^2*nq*Tf - 8*CA*CF*Mt^2*nq*Tf + 64*CF^2*Mt^2*nq*Tf + 96*CA*CF*Me^2*nt*Tf - 24*CF^2*Me^2*nt*Tf + 
    432*CF^2*Mfg^2*nt*Tf + 61*CA*CF*Mst1^2*nt*Tf - 32*CF^2*Mst1^2*nt*Tf - CA*CF*cthetat^4*Mst1^2*nt*Tf + 8*CF^2*cthetat^4*Mst1^2*nt*Tf + 
    59*CA*CF*Mst2^2*nt*Tf - 16*CF^2*Mst2^2*nt*Tf + CA*CF*cthetat^4*Mst2^2*nt*Tf - 8*CF^2*cthetat^4*Mst2^2*nt*Tf - 128*CA*CF*Mt^2*nt*Tf + 
    112*CF^2*Mt^2*nt*Tf - 48*CA*CF*cthetat*Mfg*Mt*nq*sthetat*Tf + 384*CF^2*cthetat*Mfg*Mt*nq*sthetat*Tf - 
    48*CA*CF*cthetat*Mfg*Mt*nt*sthetat*Tf + 384*CF^2*cthetat*Mfg*Mt*nt*sthetat*Tf + 6*CA*CF*cthetat^2*Mst1^2*nq*sthetat^2*Tf - 
    48*CF^2*cthetat^2*Mst1^2*nq*sthetat^2*Tf - 6*CA*CF*cthetat^2*Mst2^2*nq*sthetat^2*Tf + 48*CF^2*cthetat^2*Mst2^2*nq*sthetat^2*Tf + 
    6*CA*CF*cthetat^2*Mst1^2*nt*sthetat^2*Tf - 48*CF^2*cthetat^2*Mst1^2*nt*sthetat^2*Tf - 6*CA*CF*cthetat^2*Mst2^2*nt*sthetat^2*Tf + 
    48*CF^2*cthetat^2*Mst2^2*nt*sthetat^2*Tf - CA*CF*Mst1^2*nq*sthetat^4*Tf + 8*CF^2*Mst1^2*nq*sthetat^4*Tf + CA*CF*Mst2^2*nq*sthetat^4*Tf - 
    8*CF^2*Mst2^2*nq*sthetat^4*Tf - CA*CF*Mst1^2*nt*sthetat^4*Tf + 8*CF^2*Mst1^2*nt*sthetat^4*Tf + CA*CF*Mst2^2*nt*sthetat^4*Tf - 
    8*CF^2*Mst2^2*nt*sthetat^4*Tf - 24*CF*Me^2*nq^2*Tf^2 + 96*CF*Mfg^2*nq^2*Tf^2 - 48*CF*Msq^2*nq^2*Tf^2 - 2*CF*Mst1^2*nq^2*Tf^2 + 
    2*CF*cthetat^4*Mst1^2*nq^2*Tf^2 + 2*CF*Mst2^2*nq^2*Tf^2 - 2*CF*cthetat^4*Mst2^2*nq^2*Tf^2 + 16*CF*Mt^2*nq^2*Tf^2 - 
    48*CF*Me^2*nq*nt*Tf^2 + 192*CF*Mfg^2*nq*nt*Tf^2 - 48*CF*Msq^2*nq*nt*Tf^2 - 28*CF*Mst1^2*nq*nt*Tf^2 + 4*CF*cthetat^4*Mst1^2*nq*nt*Tf^2 - 
    20*CF*Mst2^2*nq*nt*Tf^2 - 4*CF*cthetat^4*Mst2^2*nq*nt*Tf^2 + 80*CF*Mt^2*nq*nt*Tf^2 - 24*CF*Me^2*nt^2*Tf^2 + 96*CF*Mfg^2*nt^2*Tf^2 - 
    26*CF*Mst1^2*nt^2*Tf^2 + 2*CF*cthetat^4*Mst1^2*nt^2*Tf^2 - 22*CF*Mst2^2*nt^2*Tf^2 - 2*CF*cthetat^4*Mst2^2*nt^2*Tf^2 + 
    64*CF*Mt^2*nt^2*Tf^2 + 96*CF*cthetat*Mfg*Mt*nq^2*sthetat*Tf^2 + 192*CF*cthetat*Mfg*Mt*nq*nt*sthetat*Tf^2 + 
    96*CF*cthetat*Mfg*Mt*nt^2*sthetat*Tf^2 - 12*CF*cthetat^2*Mst1^2*nq^2*sthetat^2*Tf^2 + 12*CF*cthetat^2*Mst2^2*nq^2*sthetat^2*Tf^2 - 
    24*CF*cthetat^2*Mst1^2*nq*nt*sthetat^2*Tf^2 + 24*CF*cthetat^2*Mst2^2*nq*nt*sthetat^2*Tf^2 - 12*CF*cthetat^2*Mst1^2*nt^2*sthetat^2*Tf^2 + 
    12*CF*cthetat^2*Mst2^2*nt^2*sthetat^2*Tf^2 + 2*CF*Mst1^2*nq^2*sthetat^4*Tf^2 - 2*CF*Mst2^2*nq^2*sthetat^4*Tf^2 + 
    4*CF*Mst1^2*nq*nt*sthetat^4*Tf^2 - 4*CF*Mst2^2*nq*nt*sthetat^4*Tf^2 + 2*CF*Mst1^2*nt^2*sthetat^4*Tf^2 - 
    2*CF*Mst2^2*nt^2*sthetat^4*Tf^2 + 576*CA*CF*Mfg^2*nq*Tf*z3 - 576*CF^2*Mfg^2*nq*Tf*z3 - 12*CA*CF*Mst1^2*nq*Tf*z3 + 
    12*CF^2*Mst1^2*nq*Tf*z3 + 12*CA*CF*cthetat^4*Mst1^2*nq*Tf*z3 - 12*CF^2*cthetat^4*Mst1^2*nq*Tf*z3 + 12*CA*CF*Mst2^2*nq*Tf*z3 - 
    12*CF^2*Mst2^2*nq*Tf*z3 - 12*CA*CF*cthetat^4*Mst2^2*nq*Tf*z3 + 12*CF^2*cthetat^4*Mst2^2*nq*Tf*z3 + 96*CA*CF*Mt^2*nq*Tf*z3 - 
    96*CF^2*Mt^2*nq*Tf*z3 + 576*CA*CF*Mfg^2*nt*Tf*z3 - 576*CF^2*Mfg^2*nt*Tf*z3 - 12*CA*CF*Mst1^2*nt*Tf*z3 + 12*CF^2*Mst1^2*nt*Tf*z3 + 
    12*CA*CF*cthetat^4*Mst1^2*nt*Tf*z3 - 12*CF^2*cthetat^4*Mst1^2*nt*Tf*z3 + 12*CA*CF*Mst2^2*nt*Tf*z3 - 12*CF^2*Mst2^2*nt*Tf*z3 - 
    12*CA*CF*cthetat^4*Mst2^2*nt*Tf*z3 + 12*CF^2*cthetat^4*Mst2^2*nt*Tf*z3 + 96*CA*CF*Mt^2*nt*Tf*z3 - 96*CF^2*Mt^2*nt*Tf*z3 + 
    576*CA*CF*cthetat*Mfg*Mt*nq*sthetat*Tf*z3 - 576*CF^2*cthetat*Mfg*Mt*nq*sthetat*Tf*z3 + 576*CA*CF*cthetat*Mfg*Mt*nt*sthetat*Tf*z3 - 
    576*CF^2*cthetat*Mfg*Mt*nt*sthetat*Tf*z3 - 72*CA*CF*cthetat^2*Mst1^2*nq*sthetat^2*Tf*z3 + 72*CF^2*cthetat^2*Mst1^2*nq*sthetat^2*Tf*z3 + 
    72*CA*CF*cthetat^2*Mst2^2*nq*sthetat^2*Tf*z3 - 72*CF^2*cthetat^2*Mst2^2*nq*sthetat^2*Tf*z3 - 
    72*CA*CF*cthetat^2*Mst1^2*nt*sthetat^2*Tf*z3 + 72*CF^2*cthetat^2*Mst1^2*nt*sthetat^2*Tf*z3 + 
    72*CA*CF*cthetat^2*Mst2^2*nt*sthetat^2*Tf*z3 - 72*CF^2*cthetat^2*Mst2^2*nt*sthetat^2*Tf*z3 + 12*CA*CF*Mst1^2*nq*sthetat^4*Tf*z3 - 
    12*CF^2*Mst1^2*nq*sthetat^4*Tf*z3 - 12*CA*CF*Mst2^2*nq*sthetat^4*Tf*z3 + 12*CF^2*Mst2^2*nq*sthetat^4*Tf*z3 + 
    12*CA*CF*Mst1^2*nt*sthetat^4*Tf*z3 - 12*CF^2*Mst1^2*nt*sthetat^4*Tf*z3 - 12*CA*CF*Mst2^2*nt*sthetat^4*Tf*z3 + 
    12*CF^2*Mst2^2*nt*sthetat^4*Tf*z3))/(64*Mst2^2));


gamMsq = (-((as*CF*Mfg^2)/Msq^2) + (as^2*(-3*CA*CF*Me^2 - 11*CA*CF*Mfg^2 + 6*CF^2*Mfg^2 + 2*CF*Me^2*nq*Tf + 6*CF*Mfg^2*nq*Tf + 4*CF*Msq^2*nq*Tf + 
    2*CF*Me^2*nt*Tf + 6*CF*Mfg^2*nt*Tf + 2*CF*Mst1^2*nt*Tf + 2*CF*Mst2^2*nt*Tf - 4*CF*Mt^2*nt*Tf))/(4*Msq^2) - 
 (3*as^3*(15*CA^2*CF*Me^2 - 6*CA*CF^2*Me^2 + 40*CA^2*CF*Mfg^2 - 28*CA*CF^2*Mfg^2 + 32*CF^3*Mfg^2 - 16*CA*CF*Me^2*nq*Tf + 4*CF^2*Me^2*nq*Tf - 
    72*CF^2*Mfg^2*nq*Tf - 20*CA*CF*Msq^2*nq*Tf + 8*CF^2*Msq^2*nq*Tf - 16*CA*CF*Me^2*nt*Tf + 4*CF^2*Me^2*nt*Tf - 72*CF^2*Mfg^2*nt*Tf - 
    10*CA*CF*Mst1^2*nt*Tf + 4*CF^2*Mst1^2*nt*Tf - 10*CA*CF*Mst2^2*nt*Tf + 4*CF^2*Mst2^2*nt*Tf + 20*CA*CF*Mt^2*nt*Tf - 8*CF^2*Mt^2*nt*Tf + 
    4*CF*Me^2*nq^2*Tf^2 - 16*CF*Mfg^2*nq^2*Tf^2 + 8*CF*Msq^2*nq^2*Tf^2 + 8*CF*Me^2*nq*nt*Tf^2 - 32*CF*Mfg^2*nq*nt*Tf^2 + 
    8*CF*Msq^2*nq*nt*Tf^2 + 4*CF*Mst1^2*nq*nt*Tf^2 + 4*CF*Mst2^2*nq*nt*Tf^2 - 8*CF*Mt^2*nq*nt*Tf^2 + 4*CF*Me^2*nt^2*Tf^2 - 
    16*CF*Mfg^2*nt^2*Tf^2 + 4*CF*Mst1^2*nt^2*Tf^2 + 4*CF*Mst2^2*nt^2*Tf^2 - 8*CF*Mt^2*nt^2*Tf^2 - 96*CA*CF*Mfg^2*nq*Tf*z3 + 
    96*CF^2*Mfg^2*nq*Tf*z3 - 96*CA*CF*Mfg^2*nt*Tf*z3 + 96*CF^2*Mfg^2*nt*Tf*z3))/(32*Msq^2));


gamAt = 1/At*(as*CF*Mfg + (3*as^2*CA*CF*Mfg)/2 + (9*as^3*CA^2*CF*Mfg)/8 - as^2*CF^2*Mfg - (9*as^3*CA*CF^2*Mfg)/8 + (3*as^3*CF^3*Mfg)/2 - 
 as^2*CF*Mfg*(nq + nt)*Tf + (3*as^3*CA*CF*Mfg*(nq + nt)*Tf)/8 - 3*as^3*CF^2*Mfg*(nq + nt)*Tf - (3*as^3*CF*Mfg*(nq + nt)^2*Tf^2)/4 - (9*as^3*CA*CF*Mfg*(nq + nt)*Tf*z3)/2 + 
 (9*as^3*CF^2*Mfg*(nq + nt)*Tf*z3)/2);


gamMQ = (((-3*as^2*CA*CF*Me^2)/4 - (45*as^3*CA^2*CF*Me^2)/32 + (9*as^3*CA*CF^2*Me^2)/16 - as*CF*Mfg^2 - (11*as^2*CA*CF*Mfg^2)/4 - 
  (15*as^3*CA^2*CF*Mfg^2)/4 + (3*as^2*CF^2*Mfg^2)/2 + (21*as^3*CA*CF^2*Mfg^2)/8 - 3*as^3*CF^3*Mfg^2 + (as^2*CF*Me^2*nq*Tf)/2 + 
  (3*as^3*CA*CF*Me^2*nq*Tf)/2 - (3*as^3*CF^2*Me^2*nq*Tf)/8 + (3*as^2*CF*Mfg^2*nq*Tf)/2 + (27*as^3*CF^2*Mfg^2*nq*Tf)/4 + 
  (as^2*CF*(MD1^2 + MD2^2 + MD3^2 + MQ^2 + 2*MQ1^2 + 2*MQ2^2 + MU1^2 + MU2^2)*nq*Tf)/10 + 
  (3*as^3*CA*CF*(MD1^2 + MD2^2 + MD3^2 + MQ^2 + 2*MQ1^2 + 2*MQ2^2 + MU1^2 + MU2^2)*nq*Tf)/16 - 
  (3*as^3*CF^2*(MD1^2 + MD2^2 + MD3^2 + MQ^2 + 2*MQ1^2 + 2*MQ2^2 + MU1^2 + MU2^2)*nq*Tf)/40 + (as^2*CF*Me^2*nt*Tf)/2 + 
  (3*as^3*CA*CF*Me^2*nt*Tf)/2 - (3*as^3*CF^2*Me^2*nt*Tf)/8 + (3*as^2*CF*Mfg^2*nt*Tf)/2 + (27*as^3*CF^2*Mfg^2*nt*Tf)/4 + 
  (as^2*CF*MQ^2*nt*Tf)/2 + (15*as^3*CA*CF*MQ^2*nt*Tf)/16 - (3*as^3*CF^2*MQ^2*nt*Tf)/8 + (as^2*CF*MU^2*nt*Tf)/2 + 
  (15*as^3*CA*CF*MU^2*nt*Tf)/16 - (3*as^3*CF^2*MU^2*nt*Tf)/8 - (3*as^3*CF*Me^2*nq^2*Tf^2)/8 + (3*as^3*CF*Mfg^2*nq^2*Tf^2)/2 - 
  (3*as^3*CF*(MD1^2 + MD2^2 + MD3^2 + MQ^2 + 2*MQ1^2 + 2*MQ2^2 + MU1^2 + MU2^2)*nq^2*Tf^2)/40 - (3*as^3*CF*Me^2*nq*nt*Tf^2)/4 + 
  3*as^3*CF*Mfg^2*nq*nt*Tf^2 - (3*as^3*CF*MQ^2*nq*nt*Tf^2)/8 - (3*as^3*CF*MU^2*nq*nt*Tf^2)/8 - 
  (3*as^3*CF*(MD1^2 + MD2^2 + MD3^2 + MQ^2 + 2*MQ1^2 + 2*MQ2^2 + MU1^2 + MU2^2)*nq*nt*Tf^2)/40 - (3*as^3*CF*Me^2*nt^2*Tf^2)/8 + 
  (3*as^3*CF*Mfg^2*nt^2*Tf^2)/2 - (3*as^3*CF*MQ^2*nt^2*Tf^2)/8 - (3*as^3*CF*MU^2*nt^2*Tf^2)/8 + 9*as^3*CA*CF*Mfg^2*nq*Tf*z3 - 
  9*as^3*CF^2*Mfg^2*nq*Tf*z3 + 9*as^3*CA*CF*Mfg^2*nt*Tf*z3 - 9*as^3*CF^2*Mfg^2*nt*Tf*z3)/MQ^2);


gamMU = (((-3*as^2*CA*CF*Me^2)/4 - (45*as^3*CA^2*CF*Me^2)/32 + (9*as^3*CA*CF^2*Me^2)/16 - as*CF*Mfg^2 - (11*as^2*CA*CF*Mfg^2)/4 - 
  (15*as^3*CA^2*CF*Mfg^2)/4 + (3*as^2*CF^2*Mfg^2)/2 + (21*as^3*CA*CF^2*Mfg^2)/8 - 3*as^3*CF^3*Mfg^2 + (as^2*CF*Me^2*nq*Tf)/2 + 
  (3*as^3*CA*CF*Me^2*nq*Tf)/2 - (3*as^3*CF^2*Me^2*nq*Tf)/8 + (3*as^2*CF*Mfg^2*nq*Tf)/2 + (27*as^3*CF^2*Mfg^2*nq*Tf)/4 + 
  (as^2*CF*(MD1^2 + MD2^2 + MD3^2 + MQ^2 + 2*MQ1^2 + 2*MQ2^2 + MU1^2 + MU2^2)*nq*Tf)/10 + 
  (3*as^3*CA*CF*(MD1^2 + MD2^2 + MD3^2 + MQ^2 + 2*MQ1^2 + 2*MQ2^2 + MU1^2 + MU2^2)*nq*Tf)/16 - 
  (3*as^3*CF^2*(MD1^2 + MD2^2 + MD3^2 + MQ^2 + 2*MQ1^2 + 2*MQ2^2 + MU1^2 + MU2^2)*nq*Tf)/40 + (as^2*CF*Me^2*nt*Tf)/2 + 
  (3*as^3*CA*CF*Me^2*nt*Tf)/2 - (3*as^3*CF^2*Me^2*nt*Tf)/8 + (3*as^2*CF*Mfg^2*nt*Tf)/2 + (27*as^3*CF^2*Mfg^2*nt*Tf)/4 + 
  (as^2*CF*MQ^2*nt*Tf)/2 + (15*as^3*CA*CF*MQ^2*nt*Tf)/16 - (3*as^3*CF^2*MQ^2*nt*Tf)/8 + (as^2*CF*MU^2*nt*Tf)/2 + 
  (15*as^3*CA*CF*MU^2*nt*Tf)/16 - (3*as^3*CF^2*MU^2*nt*Tf)/8 - (3*as^3*CF*Me^2*nq^2*Tf^2)/8 + (3*as^3*CF*Mfg^2*nq^2*Tf^2)/2 - 
  (3*as^3*CF*(MD1^2 + MD2^2 + MD3^2 + MQ^2 + 2*MQ1^2 + 2*MQ2^2 + MU1^2 + MU2^2)*nq^2*Tf^2)/40 - (3*as^3*CF*Me^2*nq*nt*Tf^2)/4 + 
  3*as^3*CF*Mfg^2*nq*nt*Tf^2 - (3*as^3*CF*MQ^2*nq*nt*Tf^2)/8 - (3*as^3*CF*MU^2*nq*nt*Tf^2)/8 - 
  (3*as^3*CF*(MD1^2 + MD2^2 + MD3^2 + MQ^2 + 2*MQ1^2 + 2*MQ2^2 + MU1^2 + MU2^2)*nq*nt*Tf^2)/40 - (3*as^3*CF*Me^2*nt^2*Tf^2)/8 + 
  (3*as^3*CF*Mfg^2*nt^2*Tf^2)/2 - (3*as^3*CF*MQ^2*nt^2*Tf^2)/8 - (3*as^3*CF*MU^2*nt^2*Tf^2)/8 + 9*as^3*CA*CF*Mfg^2*nq*Tf*z3 - 
  9*as^3*CF^2*Mfg^2*nq*Tf*z3 + 9*as^3*CA*CF*Mfg^2*nt*Tf*z3 - 9*as^3*CF^2*Mfg^2*nt*Tf*z3)/MU^2);


gamMst1DRbarPrime = ((as*(-8*CF*Mfg^2 - CF*Mst1^2 + CF*cthetat^4*Mst1^2 + CF*Mst2^2 - CF*cthetat^4*Mst2^2 - 8*CF*Mt^2 + 16*CF*cthetat*Mfg*Mt*sthetat - 
    6*CF*cthetat^2*Mst1^2*sthetat^2 + 6*CF*cthetat^2*Mst2^2*sthetat^2 + CF*Mst1^2*sthetat^4 - CF*Mst2^2*sthetat^4))/(8*Mst1^2) + 
 (as^2*(-80*CA*CF*Mfg^2 + 48*CF^2*Mfg^2 - 3*CA*CF*Mst1^2 + 2*CF^2*Mst1^2 + 3*CA*CF*cthetat^4*Mst1^2 - 2*CF^2*cthetat^4*Mst1^2 + 
    3*CA*CF*Mst2^2 - 2*CF^2*Mst2^2 - 3*CA*CF*cthetat^4*Mst2^2 + 2*CF^2*cthetat^4*Mst2^2 - 24*CA*CF*Mt^2 + 16*CF^2*Mt^2 + 
    96*CA*CF*cthetat*Mfg*Mt*sthetat - 64*CF^2*cthetat*Mfg*Mt*sthetat - 18*CA*CF*cthetat^2*Mst1^2*sthetat^2 + 
    12*CF^2*cthetat^2*Mst1^2*sthetat^2 + 18*CA*CF*cthetat^2*Mst2^2*sthetat^2 - 12*CF^2*cthetat^2*Mst2^2*sthetat^2 + 
    3*CA*CF*Mst1^2*sthetat^4 - 2*CF^2*Mst1^2*sthetat^4 - 3*CA*CF*Mst2^2*sthetat^4 + 2*CF^2*Mst2^2*sthetat^4 + 48*CF*Mfg^2*nq*Tf + 
    16*CF*Msq^2*nq*Tf + 2*CF*Mst1^2*nq*Tf - 2*CF*cthetat^4*Mst1^2*nq*Tf - 2*CF*Mst2^2*nq*Tf + 2*CF*cthetat^4*Mst2^2*nq*Tf + 
    16*CF*Mt^2*nq*Tf + 48*CF*Mfg^2*nt*Tf + 10*CF*Mst1^2*nt*Tf - 2*CF*cthetat^4*Mst1^2*nt*Tf + 6*CF*Mst2^2*nt*Tf + 
    2*CF*cthetat^4*Mst2^2*nt*Tf - 64*CF*cthetat*Mfg*Mt*nq*sthetat*Tf - 64*CF*cthetat*Mfg*Mt*nt*sthetat*Tf + 
    12*CF*cthetat^2*Mst1^2*nq*sthetat^2*Tf - 12*CF*cthetat^2*Mst2^2*nq*sthetat^2*Tf + 12*CF*cthetat^2*Mst1^2*nt*sthetat^2*Tf - 
    12*CF*cthetat^2*Mst2^2*nt*sthetat^2*Tf - 2*CF*Mst1^2*nq*sthetat^4*Tf + 2*CF*Mst2^2*nq*sthetat^4*Tf - 2*CF*Mst1^2*nt*sthetat^4*Tf + 
    2*CF*Mst2^2*nt*sthetat^4*Tf))/(32*Mst1^2) + 
 (as^3*(-188*CA^2*CF*Mfg^2 + 160*CA*CF^2*Mfg^2 - 192*CF^3*Mfg^2 - 3*CA^2*CF*Mst1^2 + 3*CA*CF^2*Mst1^2 - 4*CF^3*Mst1^2 + 
    3*CA^2*CF*cthetat^4*Mst1^2 - 3*CA*CF^2*cthetat^4*Mst1^2 + 4*CF^3*cthetat^4*Mst1^2 + 3*CA^2*CF*Mst2^2 - 3*CA*CF^2*Mst2^2 + 
    4*CF^3*Mst2^2 - 3*CA^2*CF*cthetat^4*Mst2^2 + 3*CA*CF^2*cthetat^4*Mst2^2 - 4*CF^3*cthetat^4*Mst2^2 - 24*CA^2*CF*Mt^2 + 24*CA*CF^2*Mt^2 - 
    32*CF^3*Mt^2 + 144*CA^2*CF*cthetat*Mfg*Mt*sthetat - 144*CA*CF^2*cthetat*Mfg*Mt*sthetat + 192*CF^3*cthetat*Mfg*Mt*sthetat - 
    18*CA^2*CF*cthetat^2*Mst1^2*sthetat^2 + 18*CA*CF^2*cthetat^2*Mst1^2*sthetat^2 - 24*CF^3*cthetat^2*Mst1^2*sthetat^2 + 
    18*CA^2*CF*cthetat^2*Mst2^2*sthetat^2 - 18*CA*CF^2*cthetat^2*Mst2^2*sthetat^2 + 24*CF^3*cthetat^2*Mst2^2*sthetat^2 + 
    3*CA^2*CF*Mst1^2*sthetat^4 - 3*CA*CF^2*Mst1^2*sthetat^4 + 4*CF^3*Mst1^2*sthetat^4 - 3*CA^2*CF*Mst2^2*sthetat^4 + 
    3*CA*CF^2*Mst2^2*sthetat^4 - 4*CF^3*Mst2^2*sthetat^4 - 24*CA*CF*Mfg^2*nq*Tf + 400*CF^2*Mfg^2*nq*Tf + 64*CA*CF*Msq^2*nq*Tf - 
    32*CF^2*Msq^2*nq*Tf - CA*CF*Mst1^2*nq*Tf + 8*CF^2*Mst1^2*nq*Tf + CA*CF*cthetat^4*Mst1^2*nq*Tf - 8*CF^2*cthetat^4*Mst1^2*nq*Tf + 
    CA*CF*Mst2^2*nq*Tf - 8*CF^2*Mst2^2*nq*Tf - CA*CF*cthetat^4*Mst2^2*nq*Tf + 8*CF^2*cthetat^4*Mst2^2*nq*Tf - 8*CA*CF*Mt^2*nq*Tf + 
    64*CF^2*Mt^2*nq*Tf - 24*CA*CF*Mfg^2*nt*Tf + 400*CF^2*Mfg^2*nt*Tf + 31*CA*CF*Mst1^2*nt*Tf - 8*CF^2*Mst1^2*nt*Tf + 
    CA*CF*cthetat^4*Mst1^2*nt*Tf - 8*CF^2*cthetat^4*Mst1^2*nt*Tf + 33*CA*CF*Mst2^2*nt*Tf - 24*CF^2*Mst2^2*nt*Tf - 
    CA*CF*cthetat^4*Mst2^2*nt*Tf + 8*CF^2*cthetat^4*Mst2^2*nt*Tf - 72*CA*CF*Mt^2*nt*Tf + 96*CF^2*Mt^2*nt*Tf + 
    48*CA*CF*cthetat*Mfg*Mt*nq*sthetat*Tf - 384*CF^2*cthetat*Mfg*Mt*nq*sthetat*Tf + 48*CA*CF*cthetat*Mfg*Mt*nt*sthetat*Tf - 
    384*CF^2*cthetat*Mfg*Mt*nt*sthetat*Tf - 6*CA*CF*cthetat^2*Mst1^2*nq*sthetat^2*Tf + 48*CF^2*cthetat^2*Mst1^2*nq*sthetat^2*Tf + 
    6*CA*CF*cthetat^2*Mst2^2*nq*sthetat^2*Tf - 48*CF^2*cthetat^2*Mst2^2*nq*sthetat^2*Tf - 6*CA*CF*cthetat^2*Mst1^2*nt*sthetat^2*Tf + 
    48*CF^2*cthetat^2*Mst1^2*nt*sthetat^2*Tf + 6*CA*CF*cthetat^2*Mst2^2*nt*sthetat^2*Tf - 48*CF^2*cthetat^2*Mst2^2*nt*sthetat^2*Tf + 
    CA*CF*Mst1^2*nq*sthetat^4*Tf - 8*CF^2*Mst1^2*nq*sthetat^4*Tf - CA*CF*Mst2^2*nq*sthetat^4*Tf + 8*CF^2*Mst2^2*nq*sthetat^4*Tf + 
    CA*CF*Mst1^2*nt*sthetat^4*Tf - 8*CF^2*Mst1^2*nt*sthetat^4*Tf - CA*CF*Mst2^2*nt*sthetat^4*Tf + 8*CF^2*Mst2^2*nt*sthetat^4*Tf + 
    96*CF*Mfg^2*nq^2*Tf^2 - 32*CF*Msq^2*nq^2*Tf^2 + 2*CF*Mst1^2*nq^2*Tf^2 - 2*CF*cthetat^4*Mst1^2*nq^2*Tf^2 - 2*CF*Mst2^2*nq^2*Tf^2 + 
    2*CF*cthetat^4*Mst2^2*nq^2*Tf^2 + 16*CF*Mt^2*nq^2*Tf^2 + 192*CF*Mfg^2*nq*nt*Tf^2 - 32*CF*Msq^2*nq*nt*Tf^2 - 12*CF*Mst1^2*nq*nt*Tf^2 - 
    4*CF*cthetat^4*Mst1^2*nq*nt*Tf^2 - 20*CF*Mst2^2*nq*nt*Tf^2 + 4*CF*cthetat^4*Mst2^2*nq*nt*Tf^2 + 64*CF*Mt^2*nq*nt*Tf^2 + 
    96*CF*Mfg^2*nt^2*Tf^2 - 14*CF*Mst1^2*nt^2*Tf^2 - 2*CF*cthetat^4*Mst1^2*nt^2*Tf^2 - 18*CF*Mst2^2*nt^2*Tf^2 + 
    2*CF*cthetat^4*Mst2^2*nt^2*Tf^2 + 48*CF*Mt^2*nt^2*Tf^2 - 96*CF*cthetat*Mfg*Mt*nq^2*sthetat*Tf^2 - 
    192*CF*cthetat*Mfg*Mt*nq*nt*sthetat*Tf^2 - 96*CF*cthetat*Mfg*Mt*nt^2*sthetat*Tf^2 + 12*CF*cthetat^2*Mst1^2*nq^2*sthetat^2*Tf^2 - 
    12*CF*cthetat^2*Mst2^2*nq^2*sthetat^2*Tf^2 + 24*CF*cthetat^2*Mst1^2*nq*nt*sthetat^2*Tf^2 - 24*CF*cthetat^2*Mst2^2*nq*nt*sthetat^2*Tf^2 + 
    12*CF*cthetat^2*Mst1^2*nt^2*sthetat^2*Tf^2 - 12*CF*cthetat^2*Mst2^2*nt^2*sthetat^2*Tf^2 - 2*CF*Mst1^2*nq^2*sthetat^4*Tf^2 + 
    2*CF*Mst2^2*nq^2*sthetat^4*Tf^2 - 4*CF*Mst1^2*nq*nt*sthetat^4*Tf^2 + 4*CF*Mst2^2*nq*nt*sthetat^4*Tf^2 - 
    2*CF*Mst1^2*nt^2*sthetat^4*Tf^2 + 2*CF*Mst2^2*nt^2*sthetat^4*Tf^2 + 576*CA*CF*Mfg^2*nq*Tf*z3 - 576*CF^2*Mfg^2*nq*Tf*z3 + 
    12*CA*CF*Mst1^2*nq*Tf*z3 - 12*CF^2*Mst1^2*nq*Tf*z3 - 12*CA*CF*cthetat^4*Mst1^2*nq*Tf*z3 + 12*CF^2*cthetat^4*Mst1^2*nq*Tf*z3 - 
    12*CA*CF*Mst2^2*nq*Tf*z3 + 12*CF^2*Mst2^2*nq*Tf*z3 + 12*CA*CF*cthetat^4*Mst2^2*nq*Tf*z3 - 12*CF^2*cthetat^4*Mst2^2*nq*Tf*z3 + 
    96*CA*CF*Mt^2*nq*Tf*z3 - 96*CF^2*Mt^2*nq*Tf*z3 + 576*CA*CF*Mfg^2*nt*Tf*z3 - 576*CF^2*Mfg^2*nt*Tf*z3 + 12*CA*CF*Mst1^2*nt*Tf*z3 - 
    12*CF^2*Mst1^2*nt*Tf*z3 - 12*CA*CF*cthetat^4*Mst1^2*nt*Tf*z3 + 12*CF^2*cthetat^4*Mst1^2*nt*Tf*z3 - 12*CA*CF*Mst2^2*nt*Tf*z3 + 
    12*CF^2*Mst2^2*nt*Tf*z3 + 12*CA*CF*cthetat^4*Mst2^2*nt*Tf*z3 - 12*CF^2*cthetat^4*Mst2^2*nt*Tf*z3 + 96*CA*CF*Mt^2*nt*Tf*z3 - 
    96*CF^2*Mt^2*nt*Tf*z3 - 576*CA*CF*cthetat*Mfg*Mt*nq*sthetat*Tf*z3 + 576*CF^2*cthetat*Mfg*Mt*nq*sthetat*Tf*z3 - 
    576*CA*CF*cthetat*Mfg*Mt*nt*sthetat*Tf*z3 + 576*CF^2*cthetat*Mfg*Mt*nt*sthetat*Tf*z3 + 72*CA*CF*cthetat^2*Mst1^2*nq*sthetat^2*Tf*z3 - 
    72*CF^2*cthetat^2*Mst1^2*nq*sthetat^2*Tf*z3 - 72*CA*CF*cthetat^2*Mst2^2*nq*sthetat^2*Tf*z3 + 
    72*CF^2*cthetat^2*Mst2^2*nq*sthetat^2*Tf*z3 + 72*CA*CF*cthetat^2*Mst1^2*nt*sthetat^2*Tf*z3 - 
    72*CF^2*cthetat^2*Mst1^2*nt*sthetat^2*Tf*z3 - 72*CA*CF*cthetat^2*Mst2^2*nt*sthetat^2*Tf*z3 + 
    72*CF^2*cthetat^2*Mst2^2*nt*sthetat^2*Tf*z3 - 12*CA*CF*Mst1^2*nq*sthetat^4*Tf*z3 + 12*CF^2*Mst1^2*nq*sthetat^4*Tf*z3 + 
    12*CA*CF*Mst2^2*nq*sthetat^4*Tf*z3 - 12*CF^2*Mst2^2*nq*sthetat^4*Tf*z3 - 12*CA*CF*Mst1^2*nt*sthetat^4*Tf*z3 + 
    12*CF^2*Mst1^2*nt*sthetat^4*Tf*z3 + 12*CA*CF*Mst2^2*nt*sthetat^4*Tf*z3 - 12*CF^2*Mst2^2*nt*sthetat^4*Tf*z3))/(64*Mst1^2));


gamMst2DRbarPrime = ((as*(-8*CF*Mfg^2 + CF*Mst1^2 - CF*cthetat^4*Mst1^2 - CF*Mst2^2 + CF*cthetat^4*Mst2^2 - 8*CF*Mt^2 - 16*CF*cthetat*Mfg*Mt*sthetat + 
    6*CF*cthetat^2*Mst1^2*sthetat^2 - 6*CF*cthetat^2*Mst2^2*sthetat^2 - CF*Mst1^2*sthetat^4 + CF*Mst2^2*sthetat^4))/(8*Mst2^2) + 
 (as^2*(-80*CA*CF*Mfg^2 + 48*CF^2*Mfg^2 + 3*CA*CF*Mst1^2 - 2*CF^2*Mst1^2 - 3*CA*CF*cthetat^4*Mst1^2 + 2*CF^2*cthetat^4*Mst1^2 - 
    3*CA*CF*Mst2^2 + 2*CF^2*Mst2^2 + 3*CA*CF*cthetat^4*Mst2^2 - 2*CF^2*cthetat^4*Mst2^2 - 24*CA*CF*Mt^2 + 16*CF^2*Mt^2 - 
    96*CA*CF*cthetat*Mfg*Mt*sthetat + 64*CF^2*cthetat*Mfg*Mt*sthetat + 18*CA*CF*cthetat^2*Mst1^2*sthetat^2 - 
    12*CF^2*cthetat^2*Mst1^2*sthetat^2 - 18*CA*CF*cthetat^2*Mst2^2*sthetat^2 + 12*CF^2*cthetat^2*Mst2^2*sthetat^2 - 
    3*CA*CF*Mst1^2*sthetat^4 + 2*CF^2*Mst1^2*sthetat^4 + 3*CA*CF*Mst2^2*sthetat^4 - 2*CF^2*Mst2^2*sthetat^4 + 48*CF*Mfg^2*nq*Tf + 
    16*CF*Msq^2*nq*Tf - 2*CF*Mst1^2*nq*Tf + 2*CF*cthetat^4*Mst1^2*nq*Tf + 2*CF*Mst2^2*nq*Tf - 2*CF*cthetat^4*Mst2^2*nq*Tf + 
    16*CF*Mt^2*nq*Tf + 48*CF*Mfg^2*nt*Tf + 6*CF*Mst1^2*nt*Tf + 2*CF*cthetat^4*Mst1^2*nt*Tf + 10*CF*Mst2^2*nt*Tf - 
    2*CF*cthetat^4*Mst2^2*nt*Tf + 64*CF*cthetat*Mfg*Mt*nq*sthetat*Tf + 64*CF*cthetat*Mfg*Mt*nt*sthetat*Tf - 
    12*CF*cthetat^2*Mst1^2*nq*sthetat^2*Tf + 12*CF*cthetat^2*Mst2^2*nq*sthetat^2*Tf - 12*CF*cthetat^2*Mst1^2*nt*sthetat^2*Tf + 
    12*CF*cthetat^2*Mst2^2*nt*sthetat^2*Tf + 2*CF*Mst1^2*nq*sthetat^4*Tf - 2*CF*Mst2^2*nq*sthetat^4*Tf + 2*CF*Mst1^2*nt*sthetat^4*Tf - 
    2*CF*Mst2^2*nt*sthetat^4*Tf))/(32*Mst2^2) + 
 (as^3*(-188*CA^2*CF*Mfg^2 + 160*CA*CF^2*Mfg^2 - 192*CF^3*Mfg^2 + 3*CA^2*CF*Mst1^2 - 3*CA*CF^2*Mst1^2 + 4*CF^3*Mst1^2 - 
    3*CA^2*CF*cthetat^4*Mst1^2 + 3*CA*CF^2*cthetat^4*Mst1^2 - 4*CF^3*cthetat^4*Mst1^2 - 3*CA^2*CF*Mst2^2 + 3*CA*CF^2*Mst2^2 - 
    4*CF^3*Mst2^2 + 3*CA^2*CF*cthetat^4*Mst2^2 - 3*CA*CF^2*cthetat^4*Mst2^2 + 4*CF^3*cthetat^4*Mst2^2 - 24*CA^2*CF*Mt^2 + 24*CA*CF^2*Mt^2 - 
    32*CF^3*Mt^2 - 144*CA^2*CF*cthetat*Mfg*Mt*sthetat + 144*CA*CF^2*cthetat*Mfg*Mt*sthetat - 192*CF^3*cthetat*Mfg*Mt*sthetat + 
    18*CA^2*CF*cthetat^2*Mst1^2*sthetat^2 - 18*CA*CF^2*cthetat^2*Mst1^2*sthetat^2 + 24*CF^3*cthetat^2*Mst1^2*sthetat^2 - 
    18*CA^2*CF*cthetat^2*Mst2^2*sthetat^2 + 18*CA*CF^2*cthetat^2*Mst2^2*sthetat^2 - 24*CF^3*cthetat^2*Mst2^2*sthetat^2 - 
    3*CA^2*CF*Mst1^2*sthetat^4 + 3*CA*CF^2*Mst1^2*sthetat^4 - 4*CF^3*Mst1^2*sthetat^4 + 3*CA^2*CF*Mst2^2*sthetat^4 - 
    3*CA*CF^2*Mst2^2*sthetat^4 + 4*CF^3*Mst2^2*sthetat^4 - 24*CA*CF*Mfg^2*nq*Tf + 400*CF^2*Mfg^2*nq*Tf + 64*CA*CF*Msq^2*nq*Tf - 
    32*CF^2*Msq^2*nq*Tf + CA*CF*Mst1^2*nq*Tf - 8*CF^2*Mst1^2*nq*Tf - CA*CF*cthetat^4*Mst1^2*nq*Tf + 8*CF^2*cthetat^4*Mst1^2*nq*Tf - 
    CA*CF*Mst2^2*nq*Tf + 8*CF^2*Mst2^2*nq*Tf + CA*CF*cthetat^4*Mst2^2*nq*Tf - 8*CF^2*cthetat^4*Mst2^2*nq*Tf - 8*CA*CF*Mt^2*nq*Tf + 
    64*CF^2*Mt^2*nq*Tf - 24*CA*CF*Mfg^2*nt*Tf + 400*CF^2*Mfg^2*nt*Tf + 33*CA*CF*Mst1^2*nt*Tf - 24*CF^2*Mst1^2*nt*Tf - 
    CA*CF*cthetat^4*Mst1^2*nt*Tf + 8*CF^2*cthetat^4*Mst1^2*nt*Tf + 31*CA*CF*Mst2^2*nt*Tf - 8*CF^2*Mst2^2*nt*Tf + 
    CA*CF*cthetat^4*Mst2^2*nt*Tf - 8*CF^2*cthetat^4*Mst2^2*nt*Tf - 72*CA*CF*Mt^2*nt*Tf + 96*CF^2*Mt^2*nt*Tf - 
    48*CA*CF*cthetat*Mfg*Mt*nq*sthetat*Tf + 384*CF^2*cthetat*Mfg*Mt*nq*sthetat*Tf - 48*CA*CF*cthetat*Mfg*Mt*nt*sthetat*Tf + 
    384*CF^2*cthetat*Mfg*Mt*nt*sthetat*Tf + 6*CA*CF*cthetat^2*Mst1^2*nq*sthetat^2*Tf - 48*CF^2*cthetat^2*Mst1^2*nq*sthetat^2*Tf - 
    6*CA*CF*cthetat^2*Mst2^2*nq*sthetat^2*Tf + 48*CF^2*cthetat^2*Mst2^2*nq*sthetat^2*Tf + 6*CA*CF*cthetat^2*Mst1^2*nt*sthetat^2*Tf - 
    48*CF^2*cthetat^2*Mst1^2*nt*sthetat^2*Tf - 6*CA*CF*cthetat^2*Mst2^2*nt*sthetat^2*Tf + 48*CF^2*cthetat^2*Mst2^2*nt*sthetat^2*Tf - 
    CA*CF*Mst1^2*nq*sthetat^4*Tf + 8*CF^2*Mst1^2*nq*sthetat^4*Tf + CA*CF*Mst2^2*nq*sthetat^4*Tf - 8*CF^2*Mst2^2*nq*sthetat^4*Tf - 
    CA*CF*Mst1^2*nt*sthetat^4*Tf + 8*CF^2*Mst1^2*nt*sthetat^4*Tf + CA*CF*Mst2^2*nt*sthetat^4*Tf - 8*CF^2*Mst2^2*nt*sthetat^4*Tf + 
    96*CF*Mfg^2*nq^2*Tf^2 - 32*CF*Msq^2*nq^2*Tf^2 - 2*CF*Mst1^2*nq^2*Tf^2 + 2*CF*cthetat^4*Mst1^2*nq^2*Tf^2 + 2*CF*Mst2^2*nq^2*Tf^2 - 
    2*CF*cthetat^4*Mst2^2*nq^2*Tf^2 + 16*CF*Mt^2*nq^2*Tf^2 + 192*CF*Mfg^2*nq*nt*Tf^2 - 32*CF*Msq^2*nq*nt*Tf^2 - 20*CF*Mst1^2*nq*nt*Tf^2 + 
    4*CF*cthetat^4*Mst1^2*nq*nt*Tf^2 - 12*CF*Mst2^2*nq*nt*Tf^2 - 4*CF*cthetat^4*Mst2^2*nq*nt*Tf^2 + 64*CF*Mt^2*nq*nt*Tf^2 + 
    96*CF*Mfg^2*nt^2*Tf^2 - 18*CF*Mst1^2*nt^2*Tf^2 + 2*CF*cthetat^4*Mst1^2*nt^2*Tf^2 - 14*CF*Mst2^2*nt^2*Tf^2 - 
    2*CF*cthetat^4*Mst2^2*nt^2*Tf^2 + 48*CF*Mt^2*nt^2*Tf^2 + 96*CF*cthetat*Mfg*Mt*nq^2*sthetat*Tf^2 + 
    192*CF*cthetat*Mfg*Mt*nq*nt*sthetat*Tf^2 + 96*CF*cthetat*Mfg*Mt*nt^2*sthetat*Tf^2 - 12*CF*cthetat^2*Mst1^2*nq^2*sthetat^2*Tf^2 + 
    12*CF*cthetat^2*Mst2^2*nq^2*sthetat^2*Tf^2 - 24*CF*cthetat^2*Mst1^2*nq*nt*sthetat^2*Tf^2 + 24*CF*cthetat^2*Mst2^2*nq*nt*sthetat^2*Tf^2 - 
    12*CF*cthetat^2*Mst1^2*nt^2*sthetat^2*Tf^2 + 12*CF*cthetat^2*Mst2^2*nt^2*sthetat^2*Tf^2 + 2*CF*Mst1^2*nq^2*sthetat^4*Tf^2 - 
    2*CF*Mst2^2*nq^2*sthetat^4*Tf^2 + 4*CF*Mst1^2*nq*nt*sthetat^4*Tf^2 - 4*CF*Mst2^2*nq*nt*sthetat^4*Tf^2 + 
    2*CF*Mst1^2*nt^2*sthetat^4*Tf^2 - 2*CF*Mst2^2*nt^2*sthetat^4*Tf^2 + 576*CA*CF*Mfg^2*nq*Tf*z3 - 576*CF^2*Mfg^2*nq*Tf*z3 - 
    12*CA*CF*Mst1^2*nq*Tf*z3 + 12*CF^2*Mst1^2*nq*Tf*z3 + 12*CA*CF*cthetat^4*Mst1^2*nq*Tf*z3 - 12*CF^2*cthetat^4*Mst1^2*nq*Tf*z3 + 
    12*CA*CF*Mst2^2*nq*Tf*z3 - 12*CF^2*Mst2^2*nq*Tf*z3 - 12*CA*CF*cthetat^4*Mst2^2*nq*Tf*z3 + 12*CF^2*cthetat^4*Mst2^2*nq*Tf*z3 + 
    96*CA*CF*Mt^2*nq*Tf*z3 - 96*CF^2*Mt^2*nq*Tf*z3 + 576*CA*CF*Mfg^2*nt*Tf*z3 - 576*CF^2*Mfg^2*nt*Tf*z3 - 12*CA*CF*Mst1^2*nt*Tf*z3 + 
    12*CF^2*Mst1^2*nt*Tf*z3 + 12*CA*CF*cthetat^4*Mst1^2*nt*Tf*z3 - 12*CF^2*cthetat^4*Mst1^2*nt*Tf*z3 + 12*CA*CF*Mst2^2*nt*Tf*z3 - 
    12*CF^2*Mst2^2*nt*Tf*z3 - 12*CA*CF*cthetat^4*Mst2^2*nt*Tf*z3 + 12*CF^2*cthetat^4*Mst2^2*nt*Tf*z3 + 96*CA*CF*Mt^2*nt*Tf*z3 - 
    96*CF^2*Mt^2*nt*Tf*z3 + 576*CA*CF*cthetat*Mfg*Mt*nq*sthetat*Tf*z3 - 576*CF^2*cthetat*Mfg*Mt*nq*sthetat*Tf*z3 + 
    576*CA*CF*cthetat*Mfg*Mt*nt*sthetat*Tf*z3 - 576*CF^2*cthetat*Mfg*Mt*nt*sthetat*Tf*z3 - 72*CA*CF*cthetat^2*Mst1^2*nq*sthetat^2*Tf*z3 + 
    72*CF^2*cthetat^2*Mst1^2*nq*sthetat^2*Tf*z3 + 72*CA*CF*cthetat^2*Mst2^2*nq*sthetat^2*Tf*z3 - 
    72*CF^2*cthetat^2*Mst2^2*nq*sthetat^2*Tf*z3 - 72*CA*CF*cthetat^2*Mst1^2*nt*sthetat^2*Tf*z3 + 
    72*CF^2*cthetat^2*Mst1^2*nt*sthetat^2*Tf*z3 + 72*CA*CF*cthetat^2*Mst2^2*nt*sthetat^2*Tf*z3 - 
    72*CF^2*cthetat^2*Mst2^2*nt*sthetat^2*Tf*z3 + 12*CA*CF*Mst1^2*nq*sthetat^4*Tf*z3 - 12*CF^2*Mst1^2*nq*sthetat^4*Tf*z3 - 
    12*CA*CF*Mst2^2*nq*sthetat^4*Tf*z3 + 12*CF^2*Mst2^2*nq*sthetat^4*Tf*z3 + 12*CA*CF*Mst1^2*nt*sthetat^4*Tf*z3 - 
    12*CF^2*Mst1^2*nt*sthetat^4*Tf*z3 - 12*CA*CF*Mst2^2*nt*sthetat^4*Tf*z3 + 12*CF^2*Mst2^2*nt*sthetat^4*Tf*z3))/(64*Mst2^2));


gamMsqDRbarPrime = (-((as*CF*Mfg^2)/Msq^2) + (as^2*(-10*CA*CF*Mfg^2 + 6*CF^2*Mfg^2 + 6*CF*Mfg^2*nq*Tf + 2*CF*Msq^2*nq*Tf + 6*CF*Mfg^2*nt*Tf + CF*Mst1^2*nt*Tf + 
    CF*Mst2^2*nt*Tf - 2*CF*Mt^2*nt*Tf))/(4*Msq^2) + 
 (as^3*(-47*CA^2*CF*Mfg^2 + 40*CA*CF^2*Mfg^2 - 48*CF^3*Mfg^2 - 6*CA*CF*Mfg^2*nq*Tf + 100*CF^2*Mfg^2*nq*Tf + 16*CA*CF*Msq^2*nq*Tf - 
    8*CF^2*Msq^2*nq*Tf - 6*CA*CF*Mfg^2*nt*Tf + 100*CF^2*Mfg^2*nt*Tf + 8*CA*CF*Mst1^2*nt*Tf - 4*CF^2*Mst1^2*nt*Tf + 8*CA*CF*Mst2^2*nt*Tf - 
    4*CF^2*Mst2^2*nt*Tf - 16*CA*CF*Mt^2*nt*Tf + 8*CF^2*Mt^2*nt*Tf + 24*CF*Mfg^2*nq^2*Tf^2 - 8*CF*Msq^2*nq^2*Tf^2 + 48*CF*Mfg^2*nq*nt*Tf^2 - 
    8*CF*Msq^2*nq*nt*Tf^2 - 4*CF*Mst1^2*nq*nt*Tf^2 - 4*CF*Mst2^2*nq*nt*Tf^2 + 8*CF*Mt^2*nq*nt*Tf^2 + 24*CF*Mfg^2*nt^2*Tf^2 - 
    4*CF*Mst1^2*nt^2*Tf^2 - 4*CF*Mst2^2*nt^2*Tf^2 + 8*CF*Mt^2*nt^2*Tf^2 + 144*CA*CF*Mfg^2*nq*Tf*z3 - 144*CF^2*Mfg^2*nq*Tf*z3 + 
    144*CA*CF*Mfg^2*nt*Tf*z3 - 144*CF^2*Mfg^2*nt*Tf*z3))/(16*Msq^2));


gamMQDRbarPrime = ((-(as*CF*Mfg^2) + (as^2*(-10*CA*CF*Mfg^2 + 6*CF^2*Mfg^2 + 6*CF*Mfg^2*nq*Tf + 
     (CF*(MD1^2 + MD2^2 + MD3^2 + MQ^2 + 2*MQ1^2 + 2*MQ2^2 + MU1^2 + MU2^2)*nq*Tf)/5 + 6*CF*Mfg^2*nt*Tf + CF*MQ^2*nt*Tf + CF*MU^2*nt*Tf))/4 + 
  (as^3*(-47*CA^2*CF*Mfg^2 + 40*CA*CF^2*Mfg^2 - 48*CF^3*Mfg^2 - 6*CA*CF*Mfg^2*nq*Tf + 100*CF^2*Mfg^2*nq*Tf + 
     (8*CA*CF*(MD1^2 + MD2^2 + MD3^2 + MQ^2 + 2*MQ1^2 + 2*MQ2^2 + MU1^2 + MU2^2)*nq*Tf)/5 - 
     (4*CF^2*(MD1^2 + MD2^2 + MD3^2 + MQ^2 + 2*MQ1^2 + 2*MQ2^2 + MU1^2 + MU2^2)*nq*Tf)/5 - 6*CA*CF*Mfg^2*nt*Tf + 100*CF^2*Mfg^2*nt*Tf + 
     8*CA*CF*MQ^2*nt*Tf - 4*CF^2*MQ^2*nt*Tf + 8*CA*CF*MU^2*nt*Tf - 4*CF^2*MU^2*nt*Tf + 24*CF*Mfg^2*nq^2*Tf^2 - 
     (4*CF*(MD1^2 + MD2^2 + MD3^2 + MQ^2 + 2*MQ1^2 + 2*MQ2^2 + MU1^2 + MU2^2)*nq^2*Tf^2)/5 + 48*CF*Mfg^2*nq*nt*Tf^2 - 4*CF*MQ^2*nq*nt*Tf^2 - 
     4*CF*MU^2*nq*nt*Tf^2 - (4*CF*(MD1^2 + MD2^2 + MD3^2 + MQ^2 + 2*MQ1^2 + 2*MQ2^2 + MU1^2 + MU2^2)*nq*nt*Tf^2)/5 + 24*CF*Mfg^2*nt^2*Tf^2 - 
     4*CF*MQ^2*nt^2*Tf^2 - 4*CF*MU^2*nt^2*Tf^2 + 144*CA*CF*Mfg^2*nq*Tf*z3 - 144*CF^2*Mfg^2*nq*Tf*z3 + 144*CA*CF*Mfg^2*nt*Tf*z3 - 
     144*CF^2*Mfg^2*nt*Tf*z3))/16)/MQ^2);


gamMUDRbarPrime = ((-(as*CF*Mfg^2) + (as^2*(-10*CA*CF*Mfg^2 + 6*CF^2*Mfg^2 + 6*CF*Mfg^2*nq*Tf + 
     (CF*(MD1^2 + MD2^2 + MD3^2 + MQ^2 + 2*MQ1^2 + 2*MQ2^2 + MU1^2 + MU2^2)*nq*Tf)/5 + 6*CF*Mfg^2*nt*Tf + CF*MQ^2*nt*Tf + CF*MU^2*nt*Tf))/4 + 
  (as^3*(-47*CA^2*CF*Mfg^2 + 40*CA*CF^2*Mfg^2 - 48*CF^3*Mfg^2 - 6*CA*CF*Mfg^2*nq*Tf + 100*CF^2*Mfg^2*nq*Tf + 
     (8*CA*CF*(MD1^2 + MD2^2 + MD3^2 + MQ^2 + 2*MQ1^2 + 2*MQ2^2 + MU1^2 + MU2^2)*nq*Tf)/5 - 
     (4*CF^2*(MD1^2 + MD2^2 + MD3^2 + MQ^2 + 2*MQ1^2 + 2*MQ2^2 + MU1^2 + MU2^2)*nq*Tf)/5 - 6*CA*CF*Mfg^2*nt*Tf + 100*CF^2*Mfg^2*nt*Tf + 
     8*CA*CF*MQ^2*nt*Tf - 4*CF^2*MQ^2*nt*Tf + 8*CA*CF*MU^2*nt*Tf - 4*CF^2*MU^2*nt*Tf + 24*CF*Mfg^2*nq^2*Tf^2 - 
     (4*CF*(MD1^2 + MD2^2 + MD3^2 + MQ^2 + 2*MQ1^2 + 2*MQ2^2 + MU1^2 + MU2^2)*nq^2*Tf^2)/5 + 48*CF*Mfg^2*nq*nt*Tf^2 - 4*CF*MQ^2*nq*nt*Tf^2 - 
     4*CF*MU^2*nq*nt*Tf^2 - (4*CF*(MD1^2 + MD2^2 + MD3^2 + MQ^2 + 2*MQ1^2 + 2*MQ2^2 + MU1^2 + MU2^2)*nq*nt*Tf^2)/5 + 24*CF*Mfg^2*nt^2*Tf^2 - 
     4*CF*MQ^2*nt^2*Tf^2 - 4*CF*MU^2*nt^2*Tf^2 + 144*CA*CF*Mfg^2*nq*Tf*z3 - 144*CF^2*Mfg^2*nq*Tf*z3 + 144*CA*CF*Mfg^2*nt*Tf*z3 - 
     144*CF^2*Mfg^2*nt*Tf*z3))/16)/MU^2);


gamMst1withMeOS = ((as*(-8*CF*Mfg^2 - CF*Mst1^2 + CF*cthetat^4*Mst1^2 + CF*Mst2^2 - CF*cthetat^4*Mst2^2 - 8*CF*Mt^2 + 16*CF*cthetat*Mfg*Mt*sthetat - 
    6*CF*cthetat^2*Mst1^2*sthetat^2 + 6*CF*cthetat^2*Mst2^2*sthetat^2 + CF*Mst1^2*sthetat^4 - CF*Mst2^2*sthetat^4))/(8*Mst1^2) + 
 (as^2*(-24*CA*CF*MeOS^2 - 88*CA*CF*Mfg^2 + 48*CF^2*Mfg^2 - 3*CA*CF*Mst1^2 + 2*CF^2*Mst1^2 + 3*CA*CF*cthetat^4*Mst1^2 - 
    2*CF^2*cthetat^4*Mst1^2 + 3*CA*CF*Mst2^2 - 2*CF^2*Mst2^2 - 3*CA*CF*cthetat^4*Mst2^2 + 2*CF^2*cthetat^4*Mst2^2 - 24*CA*CF*Mt^2 + 
    16*CF^2*Mt^2 + 96*CA*CF*cthetat*Mfg*Mt*sthetat - 64*CF^2*cthetat*Mfg*Mt*sthetat - 18*CA*CF*cthetat^2*Mst1^2*sthetat^2 + 
    12*CF^2*cthetat^2*Mst1^2*sthetat^2 + 18*CA*CF*cthetat^2*Mst2^2*sthetat^2 - 12*CF^2*cthetat^2*Mst2^2*sthetat^2 + 
    3*CA*CF*Mst1^2*sthetat^4 - 2*CF^2*Mst1^2*sthetat^4 - 3*CA*CF*Mst2^2*sthetat^4 + 2*CF^2*Mst2^2*sthetat^4 + 16*CF*MeOS^2*nq*Tf + 
    48*CF*Mfg^2*nq*Tf + 32*CF*Msq^2*nq*Tf + 2*CF*Mst1^2*nq*Tf - 2*CF*cthetat^4*Mst1^2*nq*Tf - 2*CF*Mst2^2*nq*Tf + 
    2*CF*cthetat^4*Mst2^2*nq*Tf + 16*CF*Mt^2*nq*Tf + 16*CF*MeOS^2*nt*Tf + 48*CF*Mfg^2*nt*Tf + 18*CF*Mst1^2*nt*Tf - 
    2*CF*cthetat^4*Mst1^2*nt*Tf + 14*CF*Mst2^2*nt*Tf + 2*CF*cthetat^4*Mst2^2*nt*Tf - 16*CF*Mt^2*nt*Tf - 64*CF*cthetat*Mfg*Mt*nq*sthetat*Tf - 
    64*CF*cthetat*Mfg*Mt*nt*sthetat*Tf + 12*CF*cthetat^2*Mst1^2*nq*sthetat^2*Tf - 12*CF*cthetat^2*Mst2^2*nq*sthetat^2*Tf + 
    12*CF*cthetat^2*Mst1^2*nt*sthetat^2*Tf - 12*CF*cthetat^2*Mst2^2*nt*sthetat^2*Tf - 2*CF*Mst1^2*nq*sthetat^4*Tf + 
    2*CF*Mst2^2*nq*sthetat^4*Tf - 2*CF*Mst1^2*nt*sthetat^4*Tf + 2*CF*Mst2^2*nt*sthetat^4*Tf))/(32*Mst1^2) + 
 as^3*((-3*CA^2*CF)/64 + (3*CA*CF^2)/64 - CF^3/16 + (3*CA^2*CF*cthetat^4)/64 - (3*CA*CF^2*cthetat^4)/64 + (CF^3*cthetat^4)/16 + 
   (9*CA^2*CF*MeOS^2)/(32*Mst1^2) + (9*CA*CF^2*MeOS^2)/(16*Mst1^2) - (3*CA^2*CF*MeOS^4)/(64*Mfg^2*Mst1^2) - (51*CA^2*CF*Mfg^2)/(16*Mst1^2) + 
   (21*CA*CF^2*Mfg^2)/(8*Mst1^2) - (3*CF^3*Mfg^2)/Mst1^2 + (3*CA^2*CF*Mst2^2)/(64*Mst1^2) - (3*CA*CF^2*Mst2^2)/(64*Mst1^2) + 
   (CF^3*Mst2^2)/(16*Mst1^2) - (3*CA^2*CF*cthetat^4*Mst2^2)/(64*Mst1^2) + (3*CA*CF^2*cthetat^4*Mst2^2)/(64*Mst1^2) - 
   (CF^3*cthetat^4*Mst2^2)/(16*Mst1^2) - (3*CA^2*CF*Mt^2)/(8*Mst1^2) + (3*CA*CF^2*Mt^2)/(8*Mst1^2) - (CF^3*Mt^2)/(2*Mst1^2) + 
   (9*CA^2*CF*cthetat*Mfg*Mt*sthetat)/(4*Mst1^2) - (9*CA*CF^2*cthetat*Mfg*Mt*sthetat)/(4*Mst1^2) + (3*CF^3*cthetat*Mfg*Mt*sthetat)/Mst1^2 - 
   (9*CA^2*CF*cthetat^2*sthetat^2)/32 + (9*CA*CF^2*cthetat^2*sthetat^2)/32 - (3*CF^3*cthetat^2*sthetat^2)/8 + 
   (9*CA^2*CF*cthetat^2*Mst2^2*sthetat^2)/(32*Mst1^2) - (9*CA*CF^2*cthetat^2*Mst2^2*sthetat^2)/(32*Mst1^2) + 
   (3*CF^3*cthetat^2*Mst2^2*sthetat^2)/(8*Mst1^2) + (3*CA^2*CF*sthetat^4)/64 - (3*CA*CF^2*sthetat^4)/64 + (CF^3*sthetat^4)/16 - 
   (3*CA^2*CF*Mst2^2*sthetat^4)/(64*Mst1^2) + (3*CA*CF^2*Mst2^2*sthetat^4)/(64*Mst1^2) - (CF^3*Mst2^2*sthetat^4)/(16*Mst1^2) - 
   (CA*CF*nq*Tf)/64 + (CF^2*nq*Tf)/8 + (CA*CF*cthetat^4*nq*Tf)/64 - (CF^2*cthetat^4*nq*Tf)/8 - (3*CA*CF*MeOS^2*nq*Tf)/(4*Mst1^2) - 
   (3*CF^2*MeOS^2*nq*Tf)/(8*Mst1^2) + (CA*CF*MeOS^4*nq*Tf)/(32*Mfg^2*Mst1^2) - (3*CA*CF*Mfg^2*nq*Tf)/(8*Mst1^2) + 
   (27*CF^2*Mfg^2*nq*Tf)/(4*Mst1^2) + (3*CA*CF*Msq^2*nq*Tf)/(4*Mst1^2) - (3*CF^2*Msq^2*nq*Tf)/(4*Mst1^2) + 
   (CA*CF*Mst2^2*nq*Tf)/(64*Mst1^2) - (CF^2*Mst2^2*nq*Tf)/(8*Mst1^2) - (CA*CF*cthetat^4*Mst2^2*nq*Tf)/(64*Mst1^2) + 
   (CF^2*cthetat^4*Mst2^2*nq*Tf)/(8*Mst1^2) - (CA*CF*Mt^2*nq*Tf)/(8*Mst1^2) + (CF^2*Mt^2*nq*Tf)/Mst1^2 + (23*CA*CF*nt*Tf)/64 - 
   (CF^2*nt*Tf)/4 + (CA*CF*cthetat^4*nt*Tf)/64 - (CF^2*cthetat^4*nt*Tf)/8 + (3*CA*CF*MeOS^2*nt*Tf)/(8*Mst1^2) - 
   (3*CF^2*MeOS^2*nt*Tf)/(8*Mst1^2) + (CA*CF*MeOS^4*nt*Tf)/(32*Mfg^2*Mst1^2) - (3*CA*CF*Mfg^2*nt*Tf)/(8*Mst1^2) + 
   (27*CF^2*Mfg^2*nt*Tf)/(4*Mst1^2) + (25*CA*CF*Mst2^2*nt*Tf)/(64*Mst1^2) - (CF^2*Mst2^2*nt*Tf)/(2*Mst1^2) - 
   (CA*CF*cthetat^4*Mst2^2*nt*Tf)/(64*Mst1^2) + (CF^2*cthetat^4*Mst2^2*nt*Tf)/(8*Mst1^2) - (3*CA*CF*MeOS^4*nt*Tf)/(32*Mst1^2*Mt^2) - 
   (7*CA*CF*Mt^2*nt*Tf)/(8*Mst1^2) + (7*CF^2*Mt^2*nt*Tf)/(4*Mst1^2) + (3*CA*CF*cthetat*Mfg*Mt*nq*sthetat*Tf)/(4*Mst1^2) - 
   (6*CF^2*cthetat*Mfg*Mt*nq*sthetat*Tf)/Mst1^2 + (3*CA*CF*cthetat*Mfg*Mt*nt*sthetat*Tf)/(4*Mst1^2) - 
   (6*CF^2*cthetat*Mfg*Mt*nt*sthetat*Tf)/Mst1^2 - (3*CA*CF*cthetat^2*nq*sthetat^2*Tf)/32 + (3*CF^2*cthetat^2*nq*sthetat^2*Tf)/4 + 
   (3*CA*CF*cthetat^2*Mst2^2*nq*sthetat^2*Tf)/(32*Mst1^2) - (3*CF^2*cthetat^2*Mst2^2*nq*sthetat^2*Tf)/(4*Mst1^2) - 
   (3*CA*CF*cthetat^2*nt*sthetat^2*Tf)/32 + (3*CF^2*cthetat^2*nt*sthetat^2*Tf)/4 + (3*CA*CF*cthetat^2*Mst2^2*nt*sthetat^2*Tf)/(32*Mst1^2) - 
   (3*CF^2*cthetat^2*Mst2^2*nt*sthetat^2*Tf)/(4*Mst1^2) + (CA*CF*nq*sthetat^4*Tf)/64 - (CF^2*nq*sthetat^4*Tf)/8 - 
   (CA*CF*Mst2^2*nq*sthetat^4*Tf)/(64*Mst1^2) + (CF^2*Mst2^2*nq*sthetat^4*Tf)/(8*Mst1^2) + (CA*CF*nt*sthetat^4*Tf)/64 - 
   (CF^2*nt*sthetat^4*Tf)/8 - (CA*CF*Mst2^2*nt*sthetat^4*Tf)/(64*Mst1^2) + (CF^2*Mst2^2*nt*sthetat^4*Tf)/(8*Mst1^2) + (CF*nq^2*Tf^2)/32 - 
   (CF*cthetat^4*nq^2*Tf^2)/32 + (3*CF*MeOS^2*nq^2*Tf^2)/(8*Mst1^2) + (3*CF*Mfg^2*nq^2*Tf^2)/(2*Mst1^2) - 
   (CF*Mst2^2*nq^2*Tf^2)/(32*Mst1^2) + (CF*cthetat^4*Mst2^2*nq^2*Tf^2)/(32*Mst1^2) + (CF*Mt^2*nq^2*Tf^2)/(4*Mst1^2) + (CF*nq*nt*Tf^2)/16 - 
   (CF*cthetat^4*nq*nt*Tf^2)/16 + (3*CF*Mfg^2*nq*nt*Tf^2)/Mst1^2 - (CF*Mst2^2*nq*nt*Tf^2)/(16*Mst1^2) + 
   (CF*cthetat^4*Mst2^2*nq*nt*Tf^2)/(16*Mst1^2) + (CF*MeOS^4*nq*nt*Tf^2)/(16*Mst1^2*Mt^2) + (CF*Mt^2*nq*nt*Tf^2)/(2*Mst1^2) + 
   (CF*nt^2*Tf^2)/32 - (CF*cthetat^4*nt^2*Tf^2)/32 - (3*CF*MeOS^2*nt^2*Tf^2)/(8*Mst1^2) + (3*CF*Mfg^2*nt^2*Tf^2)/(2*Mst1^2) - 
   (CF*Mst2^2*nt^2*Tf^2)/(32*Mst1^2) + (CF*cthetat^4*Mst2^2*nt^2*Tf^2)/(32*Mst1^2) + (CF*MeOS^4*nt^2*Tf^2)/(16*Mst1^2*Mt^2) + 
   (CF*Mt^2*nt^2*Tf^2)/(4*Mst1^2) - (3*CF*cthetat*Mfg*Mt*nq^2*sthetat*Tf^2)/(2*Mst1^2) - (3*CF*cthetat*Mfg*Mt*nq*nt*sthetat*Tf^2)/Mst1^2 - 
   (3*CF*cthetat*Mfg*Mt*nt^2*sthetat*Tf^2)/(2*Mst1^2) + (3*CF*cthetat^2*nq^2*sthetat^2*Tf^2)/16 - 
   (3*CF*cthetat^2*Mst2^2*nq^2*sthetat^2*Tf^2)/(16*Mst1^2) + (3*CF*cthetat^2*nq*nt*sthetat^2*Tf^2)/8 - 
   (3*CF*cthetat^2*Mst2^2*nq*nt*sthetat^2*Tf^2)/(8*Mst1^2) + (3*CF*cthetat^2*nt^2*sthetat^2*Tf^2)/16 - 
   (3*CF*cthetat^2*Mst2^2*nt^2*sthetat^2*Tf^2)/(16*Mst1^2) - (CF*nq^2*sthetat^4*Tf^2)/32 + (CF*Mst2^2*nq^2*sthetat^4*Tf^2)/(32*Mst1^2) - 
   (CF*nq*nt*sthetat^4*Tf^2)/16 + (CF*Mst2^2*nq*nt*sthetat^4*Tf^2)/(16*Mst1^2) - (CF*nt^2*sthetat^4*Tf^2)/32 + 
   (CF*Mst2^2*nt^2*sthetat^4*Tf^2)/(32*Mst1^2) + (3*CA*CF*nq*Tf*z3)/16 - (3*CF^2*nq*Tf*z3)/16 - (3*CA*CF*cthetat^4*nq*Tf*z3)/16 + 
   (3*CF^2*cthetat^4*nq*Tf*z3)/16 + (9*CA*CF*Mfg^2*nq*Tf*z3)/Mst1^2 - (9*CF^2*Mfg^2*nq*Tf*z3)/Mst1^2 - 
   (3*CA*CF*Mst2^2*nq*Tf*z3)/(16*Mst1^2) + (3*CF^2*Mst2^2*nq*Tf*z3)/(16*Mst1^2) + (3*CA*CF*cthetat^4*Mst2^2*nq*Tf*z3)/(16*Mst1^2) - 
   (3*CF^2*cthetat^4*Mst2^2*nq*Tf*z3)/(16*Mst1^2) + (3*CA*CF*Mt^2*nq*Tf*z3)/(2*Mst1^2) - (3*CF^2*Mt^2*nq*Tf*z3)/(2*Mst1^2) + 
   (3*CA*CF*nt*Tf*z3)/16 - (3*CF^2*nt*Tf*z3)/16 - (3*CA*CF*cthetat^4*nt*Tf*z3)/16 + (3*CF^2*cthetat^4*nt*Tf*z3)/16 + 
   (9*CA*CF*Mfg^2*nt*Tf*z3)/Mst1^2 - (9*CF^2*Mfg^2*nt*Tf*z3)/Mst1^2 - (3*CA*CF*Mst2^2*nt*Tf*z3)/(16*Mst1^2) + 
   (3*CF^2*Mst2^2*nt*Tf*z3)/(16*Mst1^2) + (3*CA*CF*cthetat^4*Mst2^2*nt*Tf*z3)/(16*Mst1^2) - (3*CF^2*cthetat^4*Mst2^2*nt*Tf*z3)/(16*Mst1^2) + 
   (3*CA*CF*Mt^2*nt*Tf*z3)/(2*Mst1^2) - (3*CF^2*Mt^2*nt*Tf*z3)/(2*Mst1^2) - (9*CA*CF*cthetat*Mfg*Mt*nq*sthetat*Tf*z3)/Mst1^2 + 
   (9*CF^2*cthetat*Mfg*Mt*nq*sthetat*Tf*z3)/Mst1^2 - (9*CA*CF*cthetat*Mfg*Mt*nt*sthetat*Tf*z3)/Mst1^2 + 
   (9*CF^2*cthetat*Mfg*Mt*nt*sthetat*Tf*z3)/Mst1^2 + (9*CA*CF*cthetat^2*nq*sthetat^2*Tf*z3)/8 - (9*CF^2*cthetat^2*nq*sthetat^2*Tf*z3)/8 - 
   (9*CA*CF*cthetat^2*Mst2^2*nq*sthetat^2*Tf*z3)/(8*Mst1^2) + (9*CF^2*cthetat^2*Mst2^2*nq*sthetat^2*Tf*z3)/(8*Mst1^2) + 
   (9*CA*CF*cthetat^2*nt*sthetat^2*Tf*z3)/8 - (9*CF^2*cthetat^2*nt*sthetat^2*Tf*z3)/8 - (9*CA*CF*cthetat^2*Mst2^2*nt*sthetat^2*Tf*z3)/
    (8*Mst1^2) + (9*CF^2*cthetat^2*Mst2^2*nt*sthetat^2*Tf*z3)/(8*Mst1^2) - (3*CA*CF*nq*sthetat^4*Tf*z3)/16 + 
   (3*CF^2*nq*sthetat^4*Tf*z3)/16 + (3*CA*CF*Mst2^2*nq*sthetat^4*Tf*z3)/(16*Mst1^2) - (3*CF^2*Mst2^2*nq*sthetat^4*Tf*z3)/(16*Mst1^2) - 
   (3*CA*CF*nt*sthetat^4*Tf*z3)/16 + (3*CF^2*nt*sthetat^4*Tf*z3)/16 + (3*CA*CF*Mst2^2*nt*sthetat^4*Tf*z3)/(16*Mst1^2) - 
   (3*CF^2*Mst2^2*nt*sthetat^4*Tf*z3)/(16*Mst1^2) - (9*CA*CF*MeOS^2*nq*Tf*Log[-(mu^2/MeOS^2)])/(16*Mst1^2) + 
   (3*CF*MeOS^2*nq^2*Tf^2*Log[-(mu^2/MeOS^2)])/(8*Mst1^2) + (3*CF*MeOS^2*nq*nt*Tf^2*Log[-(mu^2/MeOS^2)])/(8*Mst1^2) + 
   (9*CA^2*CF*MeOS^2*Log[mu^2/MeOS^2])/(8*Mst1^2) - (3*CA*CF*MeOS^2*nq*Tf*Log[mu^2/MeOS^2])/(4*Mst1^2) - 
   (3*CA*CF*MeOS^2*nt*Tf*Log[mu^2/MeOS^2])/(4*Mst1^2) - (9*CA^2*CF*MeOS^2*Log[mu^2/Mfg^2])/(32*Mst1^2) + 
   (9*CA^2*CF*Mfg^2*Log[mu^2/Mfg^2])/(16*Mst1^2) + (3*CA*CF*MeOS^2*nq*Tf*Log[mu^2/Mfg^2])/(16*Mst1^2) - 
   (3*CA*CF*Mfg^2*nq*Tf*Log[mu^2/Mfg^2])/(8*Mst1^2) + (3*CA*CF*MeOS^2*nt*Tf*Log[mu^2/Mfg^2])/(16*Mst1^2) - 
   (3*CA*CF*Mfg^2*nt*Tf*Log[mu^2/Mfg^2])/(8*Mst1^2) - (9*CA*CF*Msq^2*nq*Tf*Log[mu^2/Msq^2])/(8*Mst1^2) + 
   (3*CF*Msq^2*nq^2*Tf^2*Log[mu^2/Msq^2])/(4*Mst1^2) + (3*CF*Msq^2*nq*nt*Tf^2*Log[mu^2/Msq^2])/(4*Mst1^2) - 
   (9*CA*CF*nt*Tf*Log[mu^2/Mst1^2])/16 + (3*CF*nq*nt*Tf^2*Log[mu^2/Mst1^2])/8 + (3*CF*nt^2*Tf^2*Log[mu^2/Mst1^2])/8 - 
   (9*CA*CF*Mst2^2*nt*Tf*Log[mu^2/Mst2^2])/(16*Mst1^2) + (3*CF*Mst2^2*nq*nt*Tf^2*Log[mu^2/Mst2^2])/(8*Mst1^2) + 
   (3*CF*Mst2^2*nt^2*Tf^2*Log[mu^2/Mst2^2])/(8*Mst1^2) - (9*CA*CF*MeOS^2*nt*Tf*Log[mu^2/Mt^2])/(16*Mst1^2) + 
   (9*CA*CF*Mt^2*nt*Tf*Log[mu^2/Mt^2])/(8*Mst1^2) + (3*CF*MeOS^2*nq*nt*Tf^2*Log[mu^2/Mt^2])/(8*Mst1^2) - 
   (3*CF*Mt^2*nq*nt*Tf^2*Log[mu^2/Mt^2])/(4*Mst1^2) + (3*CF*MeOS^2*nt^2*Tf^2*Log[mu^2/Mt^2])/(8*Mst1^2) - 
   (3*CF*Mt^2*nt^2*Tf^2*Log[mu^2/Mt^2])/(4*Mst1^2)));


gamMst2withMeOS = ((as*(-8*CF*Mfg^2 + CF*Mst1^2 - CF*cthetat^4*Mst1^2 - CF*Mst2^2 + CF*cthetat^4*Mst2^2 - 8*CF*Mt^2 - 16*CF*cthetat*Mfg*Mt*sthetat + 
    6*CF*cthetat^2*Mst1^2*sthetat^2 - 6*CF*cthetat^2*Mst2^2*sthetat^2 - CF*Mst1^2*sthetat^4 + CF*Mst2^2*sthetat^4))/(8*Mst2^2) + 
 (as^2*(-24*CA*CF*MeOS^2 - 88*CA*CF*Mfg^2 + 48*CF^2*Mfg^2 + 3*CA*CF*Mst1^2 - 2*CF^2*Mst1^2 - 3*CA*CF*cthetat^4*Mst1^2 + 
    2*CF^2*cthetat^4*Mst1^2 - 3*CA*CF*Mst2^2 + 2*CF^2*Mst2^2 + 3*CA*CF*cthetat^4*Mst2^2 - 2*CF^2*cthetat^4*Mst2^2 - 24*CA*CF*Mt^2 + 
    16*CF^2*Mt^2 - 96*CA*CF*cthetat*Mfg*Mt*sthetat + 64*CF^2*cthetat*Mfg*Mt*sthetat + 18*CA*CF*cthetat^2*Mst1^2*sthetat^2 - 
    12*CF^2*cthetat^2*Mst1^2*sthetat^2 - 18*CA*CF*cthetat^2*Mst2^2*sthetat^2 + 12*CF^2*cthetat^2*Mst2^2*sthetat^2 - 
    3*CA*CF*Mst1^2*sthetat^4 + 2*CF^2*Mst1^2*sthetat^4 + 3*CA*CF*Mst2^2*sthetat^4 - 2*CF^2*Mst2^2*sthetat^4 + 16*CF*MeOS^2*nq*Tf + 
    48*CF*Mfg^2*nq*Tf + 32*CF*Msq^2*nq*Tf - 2*CF*Mst1^2*nq*Tf + 2*CF*cthetat^4*Mst1^2*nq*Tf + 2*CF*Mst2^2*nq*Tf - 
    2*CF*cthetat^4*Mst2^2*nq*Tf + 16*CF*Mt^2*nq*Tf + 16*CF*MeOS^2*nt*Tf + 48*CF*Mfg^2*nt*Tf + 14*CF*Mst1^2*nt*Tf + 
    2*CF*cthetat^4*Mst1^2*nt*Tf + 18*CF*Mst2^2*nt*Tf - 2*CF*cthetat^4*Mst2^2*nt*Tf - 16*CF*Mt^2*nt*Tf + 64*CF*cthetat*Mfg*Mt*nq*sthetat*Tf + 
    64*CF*cthetat*Mfg*Mt*nt*sthetat*Tf - 12*CF*cthetat^2*Mst1^2*nq*sthetat^2*Tf + 12*CF*cthetat^2*Mst2^2*nq*sthetat^2*Tf - 
    12*CF*cthetat^2*Mst1^2*nt*sthetat^2*Tf + 12*CF*cthetat^2*Mst2^2*nt*sthetat^2*Tf + 2*CF*Mst1^2*nq*sthetat^4*Tf - 
    2*CF*Mst2^2*nq*sthetat^4*Tf + 2*CF*Mst1^2*nt*sthetat^4*Tf - 2*CF*Mst2^2*nt*sthetat^4*Tf))/(32*Mst2^2) + 
 as^3*((-3*CA^2*CF)/64 + (3*CA*CF^2)/64 - CF^3/16 + (3*CA^2*CF*cthetat^4)/64 - (3*CA*CF^2*cthetat^4)/64 + (CF^3*cthetat^4)/16 + 
   (9*CA^2*CF*MeOS^2)/(32*Mst2^2) + (9*CA*CF^2*MeOS^2)/(16*Mst2^2) - (3*CA^2*CF*MeOS^4)/(64*Mfg^2*Mst2^2) - (51*CA^2*CF*Mfg^2)/(16*Mst2^2) + 
   (21*CA*CF^2*Mfg^2)/(8*Mst2^2) - (3*CF^3*Mfg^2)/Mst2^2 + (3*CA^2*CF*Mst1^2)/(64*Mst2^2) - (3*CA*CF^2*Mst1^2)/(64*Mst2^2) + 
   (CF^3*Mst1^2)/(16*Mst2^2) - (3*CA^2*CF*cthetat^4*Mst1^2)/(64*Mst2^2) + (3*CA*CF^2*cthetat^4*Mst1^2)/(64*Mst2^2) - 
   (CF^3*cthetat^4*Mst1^2)/(16*Mst2^2) - (3*CA^2*CF*Mt^2)/(8*Mst2^2) + (3*CA*CF^2*Mt^2)/(8*Mst2^2) - (CF^3*Mt^2)/(2*Mst2^2) - 
   (9*CA^2*CF*cthetat*Mfg*Mt*sthetat)/(4*Mst2^2) + (9*CA*CF^2*cthetat*Mfg*Mt*sthetat)/(4*Mst2^2) - (3*CF^3*cthetat*Mfg*Mt*sthetat)/Mst2^2 - 
   (9*CA^2*CF*cthetat^2*sthetat^2)/32 + (9*CA*CF^2*cthetat^2*sthetat^2)/32 - (3*CF^3*cthetat^2*sthetat^2)/8 + 
   (9*CA^2*CF*cthetat^2*Mst1^2*sthetat^2)/(32*Mst2^2) - (9*CA*CF^2*cthetat^2*Mst1^2*sthetat^2)/(32*Mst2^2) + 
   (3*CF^3*cthetat^2*Mst1^2*sthetat^2)/(8*Mst2^2) + (3*CA^2*CF*sthetat^4)/64 - (3*CA*CF^2*sthetat^4)/64 + (CF^3*sthetat^4)/16 - 
   (3*CA^2*CF*Mst1^2*sthetat^4)/(64*Mst2^2) + (3*CA*CF^2*Mst1^2*sthetat^4)/(64*Mst2^2) - (CF^3*Mst1^2*sthetat^4)/(16*Mst2^2) - 
   (CA*CF*nq*Tf)/64 + (CF^2*nq*Tf)/8 + (CA*CF*cthetat^4*nq*Tf)/64 - (CF^2*cthetat^4*nq*Tf)/8 - (3*CA*CF*MeOS^2*nq*Tf)/(4*Mst2^2) - 
   (3*CF^2*MeOS^2*nq*Tf)/(8*Mst2^2) + (CA*CF*MeOS^4*nq*Tf)/(32*Mfg^2*Mst2^2) - (3*CA*CF*Mfg^2*nq*Tf)/(8*Mst2^2) + 
   (27*CF^2*Mfg^2*nq*Tf)/(4*Mst2^2) + (3*CA*CF*Msq^2*nq*Tf)/(4*Mst2^2) - (3*CF^2*Msq^2*nq*Tf)/(4*Mst2^2) + 
   (CA*CF*Mst1^2*nq*Tf)/(64*Mst2^2) - (CF^2*Mst1^2*nq*Tf)/(8*Mst2^2) - (CA*CF*cthetat^4*Mst1^2*nq*Tf)/(64*Mst2^2) + 
   (CF^2*cthetat^4*Mst1^2*nq*Tf)/(8*Mst2^2) - (CA*CF*Mt^2*nq*Tf)/(8*Mst2^2) + (CF^2*Mt^2*nq*Tf)/Mst2^2 + (23*CA*CF*nt*Tf)/64 - 
   (CF^2*nt*Tf)/4 + (CA*CF*cthetat^4*nt*Tf)/64 - (CF^2*cthetat^4*nt*Tf)/8 + (3*CA*CF*MeOS^2*nt*Tf)/(8*Mst2^2) - 
   (3*CF^2*MeOS^2*nt*Tf)/(8*Mst2^2) + (CA*CF*MeOS^4*nt*Tf)/(32*Mfg^2*Mst2^2) - (3*CA*CF*Mfg^2*nt*Tf)/(8*Mst2^2) + 
   (27*CF^2*Mfg^2*nt*Tf)/(4*Mst2^2) + (25*CA*CF*Mst1^2*nt*Tf)/(64*Mst2^2) - (CF^2*Mst1^2*nt*Tf)/(2*Mst2^2) - 
   (CA*CF*cthetat^4*Mst1^2*nt*Tf)/(64*Mst2^2) + (CF^2*cthetat^4*Mst1^2*nt*Tf)/(8*Mst2^2) - (3*CA*CF*MeOS^4*nt*Tf)/(32*Mst2^2*Mt^2) - 
   (7*CA*CF*Mt^2*nt*Tf)/(8*Mst2^2) + (7*CF^2*Mt^2*nt*Tf)/(4*Mst2^2) - (3*CA*CF*cthetat*Mfg*Mt*nq*sthetat*Tf)/(4*Mst2^2) + 
   (6*CF^2*cthetat*Mfg*Mt*nq*sthetat*Tf)/Mst2^2 - (3*CA*CF*cthetat*Mfg*Mt*nt*sthetat*Tf)/(4*Mst2^2) + 
   (6*CF^2*cthetat*Mfg*Mt*nt*sthetat*Tf)/Mst2^2 - (3*CA*CF*cthetat^2*nq*sthetat^2*Tf)/32 + (3*CF^2*cthetat^2*nq*sthetat^2*Tf)/4 + 
   (3*CA*CF*cthetat^2*Mst1^2*nq*sthetat^2*Tf)/(32*Mst2^2) - (3*CF^2*cthetat^2*Mst1^2*nq*sthetat^2*Tf)/(4*Mst2^2) - 
   (3*CA*CF*cthetat^2*nt*sthetat^2*Tf)/32 + (3*CF^2*cthetat^2*nt*sthetat^2*Tf)/4 + (3*CA*CF*cthetat^2*Mst1^2*nt*sthetat^2*Tf)/(32*Mst2^2) - 
   (3*CF^2*cthetat^2*Mst1^2*nt*sthetat^2*Tf)/(4*Mst2^2) + (CA*CF*nq*sthetat^4*Tf)/64 - (CF^2*nq*sthetat^4*Tf)/8 - 
   (CA*CF*Mst1^2*nq*sthetat^4*Tf)/(64*Mst2^2) + (CF^2*Mst1^2*nq*sthetat^4*Tf)/(8*Mst2^2) + (CA*CF*nt*sthetat^4*Tf)/64 - 
   (CF^2*nt*sthetat^4*Tf)/8 - (CA*CF*Mst1^2*nt*sthetat^4*Tf)/(64*Mst2^2) + (CF^2*Mst1^2*nt*sthetat^4*Tf)/(8*Mst2^2) + (CF*nq^2*Tf^2)/32 - 
   (CF*cthetat^4*nq^2*Tf^2)/32 + (3*CF*MeOS^2*nq^2*Tf^2)/(8*Mst2^2) + (3*CF*Mfg^2*nq^2*Tf^2)/(2*Mst2^2) - 
   (CF*Mst1^2*nq^2*Tf^2)/(32*Mst2^2) + (CF*cthetat^4*Mst1^2*nq^2*Tf^2)/(32*Mst2^2) + (CF*Mt^2*nq^2*Tf^2)/(4*Mst2^2) + (CF*nq*nt*Tf^2)/16 - 
   (CF*cthetat^4*nq*nt*Tf^2)/16 + (3*CF*Mfg^2*nq*nt*Tf^2)/Mst2^2 - (CF*Mst1^2*nq*nt*Tf^2)/(16*Mst2^2) + 
   (CF*cthetat^4*Mst1^2*nq*nt*Tf^2)/(16*Mst2^2) + (CF*MeOS^4*nq*nt*Tf^2)/(16*Mst2^2*Mt^2) + (CF*Mt^2*nq*nt*Tf^2)/(2*Mst2^2) + 
   (CF*nt^2*Tf^2)/32 - (CF*cthetat^4*nt^2*Tf^2)/32 - (3*CF*MeOS^2*nt^2*Tf^2)/(8*Mst2^2) + (3*CF*Mfg^2*nt^2*Tf^2)/(2*Mst2^2) - 
   (CF*Mst1^2*nt^2*Tf^2)/(32*Mst2^2) + (CF*cthetat^4*Mst1^2*nt^2*Tf^2)/(32*Mst2^2) + (CF*MeOS^4*nt^2*Tf^2)/(16*Mst2^2*Mt^2) + 
   (CF*Mt^2*nt^2*Tf^2)/(4*Mst2^2) + (3*CF*cthetat*Mfg*Mt*nq^2*sthetat*Tf^2)/(2*Mst2^2) + (3*CF*cthetat*Mfg*Mt*nq*nt*sthetat*Tf^2)/Mst2^2 + 
   (3*CF*cthetat*Mfg*Mt*nt^2*sthetat*Tf^2)/(2*Mst2^2) + (3*CF*cthetat^2*nq^2*sthetat^2*Tf^2)/16 - 
   (3*CF*cthetat^2*Mst1^2*nq^2*sthetat^2*Tf^2)/(16*Mst2^2) + (3*CF*cthetat^2*nq*nt*sthetat^2*Tf^2)/8 - 
   (3*CF*cthetat^2*Mst1^2*nq*nt*sthetat^2*Tf^2)/(8*Mst2^2) + (3*CF*cthetat^2*nt^2*sthetat^2*Tf^2)/16 - 
   (3*CF*cthetat^2*Mst1^2*nt^2*sthetat^2*Tf^2)/(16*Mst2^2) - (CF*nq^2*sthetat^4*Tf^2)/32 + (CF*Mst1^2*nq^2*sthetat^4*Tf^2)/(32*Mst2^2) - 
   (CF*nq*nt*sthetat^4*Tf^2)/16 + (CF*Mst1^2*nq*nt*sthetat^4*Tf^2)/(16*Mst2^2) - (CF*nt^2*sthetat^4*Tf^2)/32 + 
   (CF*Mst1^2*nt^2*sthetat^4*Tf^2)/(32*Mst2^2) + (3*CA*CF*nq*Tf*z3)/16 - (3*CF^2*nq*Tf*z3)/16 - (3*CA*CF*cthetat^4*nq*Tf*z3)/16 + 
   (3*CF^2*cthetat^4*nq*Tf*z3)/16 + (9*CA*CF*Mfg^2*nq*Tf*z3)/Mst2^2 - (9*CF^2*Mfg^2*nq*Tf*z3)/Mst2^2 - 
   (3*CA*CF*Mst1^2*nq*Tf*z3)/(16*Mst2^2) + (3*CF^2*Mst1^2*nq*Tf*z3)/(16*Mst2^2) + (3*CA*CF*cthetat^4*Mst1^2*nq*Tf*z3)/(16*Mst2^2) - 
   (3*CF^2*cthetat^4*Mst1^2*nq*Tf*z3)/(16*Mst2^2) + (3*CA*CF*Mt^2*nq*Tf*z3)/(2*Mst2^2) - (3*CF^2*Mt^2*nq*Tf*z3)/(2*Mst2^2) + 
   (3*CA*CF*nt*Tf*z3)/16 - (3*CF^2*nt*Tf*z3)/16 - (3*CA*CF*cthetat^4*nt*Tf*z3)/16 + (3*CF^2*cthetat^4*nt*Tf*z3)/16 + 
   (9*CA*CF*Mfg^2*nt*Tf*z3)/Mst2^2 - (9*CF^2*Mfg^2*nt*Tf*z3)/Mst2^2 - (3*CA*CF*Mst1^2*nt*Tf*z3)/(16*Mst2^2) + 
   (3*CF^2*Mst1^2*nt*Tf*z3)/(16*Mst2^2) + (3*CA*CF*cthetat^4*Mst1^2*nt*Tf*z3)/(16*Mst2^2) - (3*CF^2*cthetat^4*Mst1^2*nt*Tf*z3)/(16*Mst2^2) + 
   (3*CA*CF*Mt^2*nt*Tf*z3)/(2*Mst2^2) - (3*CF^2*Mt^2*nt*Tf*z3)/(2*Mst2^2) + (9*CA*CF*cthetat*Mfg*Mt*nq*sthetat*Tf*z3)/Mst2^2 - 
   (9*CF^2*cthetat*Mfg*Mt*nq*sthetat*Tf*z3)/Mst2^2 + (9*CA*CF*cthetat*Mfg*Mt*nt*sthetat*Tf*z3)/Mst2^2 - 
   (9*CF^2*cthetat*Mfg*Mt*nt*sthetat*Tf*z3)/Mst2^2 + (9*CA*CF*cthetat^2*nq*sthetat^2*Tf*z3)/8 - (9*CF^2*cthetat^2*nq*sthetat^2*Tf*z3)/8 - 
   (9*CA*CF*cthetat^2*Mst1^2*nq*sthetat^2*Tf*z3)/(8*Mst2^2) + (9*CF^2*cthetat^2*Mst1^2*nq*sthetat^2*Tf*z3)/(8*Mst2^2) + 
   (9*CA*CF*cthetat^2*nt*sthetat^2*Tf*z3)/8 - (9*CF^2*cthetat^2*nt*sthetat^2*Tf*z3)/8 - (9*CA*CF*cthetat^2*Mst1^2*nt*sthetat^2*Tf*z3)/
    (8*Mst2^2) + (9*CF^2*cthetat^2*Mst1^2*nt*sthetat^2*Tf*z3)/(8*Mst2^2) - (3*CA*CF*nq*sthetat^4*Tf*z3)/16 + 
   (3*CF^2*nq*sthetat^4*Tf*z3)/16 + (3*CA*CF*Mst1^2*nq*sthetat^4*Tf*z3)/(16*Mst2^2) - (3*CF^2*Mst1^2*nq*sthetat^4*Tf*z3)/(16*Mst2^2) - 
   (3*CA*CF*nt*sthetat^4*Tf*z3)/16 + (3*CF^2*nt*sthetat^4*Tf*z3)/16 + (3*CA*CF*Mst1^2*nt*sthetat^4*Tf*z3)/(16*Mst2^2) - 
   (3*CF^2*Mst1^2*nt*sthetat^4*Tf*z3)/(16*Mst2^2) - (9*CA*CF*MeOS^2*nq*Tf*Log[-(mu^2/MeOS^2)])/(16*Mst2^2) + 
   (3*CF*MeOS^2*nq^2*Tf^2*Log[-(mu^2/MeOS^2)])/(8*Mst2^2) + (3*CF*MeOS^2*nq*nt*Tf^2*Log[-(mu^2/MeOS^2)])/(8*Mst2^2) + 
   (9*CA^2*CF*MeOS^2*Log[mu^2/MeOS^2])/(8*Mst2^2) - (3*CA*CF*MeOS^2*nq*Tf*Log[mu^2/MeOS^2])/(4*Mst2^2) - 
   (3*CA*CF*MeOS^2*nt*Tf*Log[mu^2/MeOS^2])/(4*Mst2^2) - (9*CA^2*CF*MeOS^2*Log[mu^2/Mfg^2])/(32*Mst2^2) + 
   (9*CA^2*CF*Mfg^2*Log[mu^2/Mfg^2])/(16*Mst2^2) + (3*CA*CF*MeOS^2*nq*Tf*Log[mu^2/Mfg^2])/(16*Mst2^2) - 
   (3*CA*CF*Mfg^2*nq*Tf*Log[mu^2/Mfg^2])/(8*Mst2^2) + (3*CA*CF*MeOS^2*nt*Tf*Log[mu^2/Mfg^2])/(16*Mst2^2) - 
   (3*CA*CF*Mfg^2*nt*Tf*Log[mu^2/Mfg^2])/(8*Mst2^2) - (9*CA*CF*Msq^2*nq*Tf*Log[mu^2/Msq^2])/(8*Mst2^2) + 
   (3*CF*Msq^2*nq^2*Tf^2*Log[mu^2/Msq^2])/(4*Mst2^2) + (3*CF*Msq^2*nq*nt*Tf^2*Log[mu^2/Msq^2])/(4*Mst2^2) - 
   (9*CA*CF*Mst1^2*nt*Tf*Log[mu^2/Mst1^2])/(16*Mst2^2) + (3*CF*Mst1^2*nq*nt*Tf^2*Log[mu^2/Mst1^2])/(8*Mst2^2) + 
   (3*CF*Mst1^2*nt^2*Tf^2*Log[mu^2/Mst1^2])/(8*Mst2^2) - (9*CA*CF*nt*Tf*Log[mu^2/Mst2^2])/16 + (3*CF*nq*nt*Tf^2*Log[mu^2/Mst2^2])/8 + 
   (3*CF*nt^2*Tf^2*Log[mu^2/Mst2^2])/8 - (9*CA*CF*MeOS^2*nt*Tf*Log[mu^2/Mt^2])/(16*Mst2^2) + 
   (9*CA*CF*Mt^2*nt*Tf*Log[mu^2/Mt^2])/(8*Mst2^2) + (3*CF*MeOS^2*nq*nt*Tf^2*Log[mu^2/Mt^2])/(8*Mst2^2) - 
   (3*CF*Mt^2*nq*nt*Tf^2*Log[mu^2/Mt^2])/(4*Mst2^2) + (3*CF*MeOS^2*nt^2*Tf^2*Log[mu^2/Mt^2])/(8*Mst2^2) - 
   (3*CF*Mt^2*nt^2*Tf^2*Log[mu^2/Mt^2])/(4*Mst2^2)));


Zmes = (1 + (as*(-3*CA*Me^2 - 2*CA*Mfg^2 + 2*Me^2*nq*Tf - 4*Mq^2*nq*Tf + 4*Msq^2*nq*Tf + 2*Me^2*nt*Tf + 2*Mst1^2*nt*Tf + 2*Mst2^2*nt*Tf - 
    4*Mt^2*nt*Tf))/(4*ep*Me^2) + (as^2*(9*CA^2*Me^2 - 6*CA^2*ep*Me^2 + 12*CA^2*Mfg^2 - 10*CA^2*ep*Mfg^2 - 12*CA*Me^2*nq*Tf + 
    4*CA*ep*Me^2*nq*Tf - 8*CA*Mfg^2*nq*Tf - 8*CF*Mfg^2*nq*Tf + 4*CA*ep*Mfg^2*nq*Tf + 8*CF*ep*Mfg^2*nq*Tf - 12*CA*Msq^2*nq*Tf + 
    8*CA*ep*Msq^2*nq*Tf - 12*CA*Me^2*nt*Tf + 4*CA*ep*Me^2*nt*Tf - 8*CA*Mfg^2*nt*Tf - 8*CF*Mfg^2*nt*Tf + 4*CA*ep*Mfg^2*nt*Tf + 
    8*CF*ep*Mfg^2*nt*Tf - 6*CA*Mst1^2*nt*Tf + 4*CA*ep*Mst1^2*nt*Tf - 6*CA*Mst2^2*nt*Tf + 4*CA*ep*Mst2^2*nt*Tf + 12*CA*Mt^2*nt*Tf - 
    8*CA*ep*Mt^2*nt*Tf + 4*Me^2*nq^2*Tf^2 + 8*Msq^2*nq^2*Tf^2 + 8*Me^2*nq*nt*Tf^2 + 8*Msq^2*nq*nt*Tf^2 + 4*Mst1^2*nq*nt*Tf^2 + 
    4*Mst2^2*nq*nt*Tf^2 - 8*Mt^2*nq*nt*Tf^2 + 4*Me^2*nt^2*Tf^2 + 4*Mst1^2*nt^2*Tf^2 + 4*Mst2^2*nt^2*Tf^2 - 8*Mt^2*nt^2*Tf^2))/(16*ep^2*Me^2));


gamMsqwithMeOS = (-((as*CF*Mfg^2)/Msq^2) + (as^2*(-3*CA*CF*MeOS^2 - 11*CA*CF*Mfg^2 + 6*CF^2*Mfg^2 + 2*CF*MeOS^2*nq*Tf + 6*CF*Mfg^2*nq*Tf + 4*CF*Msq^2*nq*Tf + 
    2*CF*MeOS^2*nt*Tf + 6*CF*Mfg^2*nt*Tf + 2*CF*Mst1^2*nt*Tf + 2*CF*Mst2^2*nt*Tf - 4*CF*Mt^2*nt*Tf))/(4*Msq^2) + 
 as^3*((9*CA^2*CF*MeOS^2)/(32*Msq^2) + (9*CA*CF^2*MeOS^2)/(16*Msq^2) - (3*CA^2*CF*MeOS^4)/(64*Mfg^2*Msq^2) - (51*CA^2*CF*Mfg^2)/(16*Msq^2) + 
   (21*CA*CF^2*Mfg^2)/(8*Msq^2) - (3*CF^3*Mfg^2)/Msq^2 + (3*CA*CF*nq*Tf)/4 - (3*CF^2*nq*Tf)/4 - (3*CA*CF*MeOS^2*nq*Tf)/(4*Msq^2) - 
   (3*CF^2*MeOS^2*nq*Tf)/(8*Msq^2) + (CA*CF*MeOS^4*nq*Tf)/(32*Mfg^2*Msq^2) - (3*CA*CF*Mfg^2*nq*Tf)/(8*Msq^2) + 
   (27*CF^2*Mfg^2*nq*Tf)/(4*Msq^2) + (3*CA*CF*MeOS^2*nt*Tf)/(8*Msq^2) - (3*CF^2*MeOS^2*nt*Tf)/(8*Msq^2) + 
   (CA*CF*MeOS^4*nt*Tf)/(32*Mfg^2*Msq^2) - (3*CA*CF*Mfg^2*nt*Tf)/(8*Msq^2) + (27*CF^2*Mfg^2*nt*Tf)/(4*Msq^2) + 
   (3*CA*CF*Mst1^2*nt*Tf)/(8*Msq^2) - (3*CF^2*Mst1^2*nt*Tf)/(8*Msq^2) + (3*CA*CF*Mst2^2*nt*Tf)/(8*Msq^2) - (3*CF^2*Mst2^2*nt*Tf)/(8*Msq^2) - 
   (3*CA*CF*MeOS^4*nt*Tf)/(32*Msq^2*Mt^2) - (3*CA*CF*Mt^2*nt*Tf)/(4*Msq^2) + (3*CF^2*Mt^2*nt*Tf)/(4*Msq^2) + 
   (3*CF*MeOS^2*nq^2*Tf^2)/(8*Msq^2) + (3*CF*Mfg^2*nq^2*Tf^2)/(2*Msq^2) + (3*CF*Mfg^2*nq*nt*Tf^2)/Msq^2 + 
   (CF*MeOS^4*nq*nt*Tf^2)/(16*Msq^2*Mt^2) - (3*CF*MeOS^2*nt^2*Tf^2)/(8*Msq^2) + (3*CF*Mfg^2*nt^2*Tf^2)/(2*Msq^2) + 
   (CF*MeOS^4*nt^2*Tf^2)/(16*Msq^2*Mt^2) + (9*CA*CF*Mfg^2*nq*Tf*z3)/Msq^2 - (9*CF^2*Mfg^2*nq*Tf*z3)/Msq^2 + (9*CA*CF*Mfg^2*nt*Tf*z3)/Msq^2 - 
   (9*CF^2*Mfg^2*nt*Tf*z3)/Msq^2 - (9*CA*CF*MeOS^2*nq*Tf*Log[-(mu^2/MeOS^2)])/(16*Msq^2) + 
   (3*CF*MeOS^2*nq^2*Tf^2*Log[-(mu^2/MeOS^2)])/(8*Msq^2) + (3*CF*MeOS^2*nq*nt*Tf^2*Log[-(mu^2/MeOS^2)])/(8*Msq^2) + 
   (9*CA^2*CF*MeOS^2*Log[mu^2/MeOS^2])/(8*Msq^2) - (3*CA*CF*MeOS^2*nq*Tf*Log[mu^2/MeOS^2])/(4*Msq^2) - 
   (3*CA*CF*MeOS^2*nt*Tf*Log[mu^2/MeOS^2])/(4*Msq^2) - (9*CA^2*CF*MeOS^2*Log[mu^2/Mfg^2])/(32*Msq^2) + 
   (9*CA^2*CF*Mfg^2*Log[mu^2/Mfg^2])/(16*Msq^2) + (3*CA*CF*MeOS^2*nq*Tf*Log[mu^2/Mfg^2])/(16*Msq^2) - 
   (3*CA*CF*Mfg^2*nq*Tf*Log[mu^2/Mfg^2])/(8*Msq^2) + (3*CA*CF*MeOS^2*nt*Tf*Log[mu^2/Mfg^2])/(16*Msq^2) - 
   (3*CA*CF*Mfg^2*nt*Tf*Log[mu^2/Mfg^2])/(8*Msq^2) - (9*CA*CF*nq*Tf*Log[mu^2/Msq^2])/8 + (3*CF*nq^2*Tf^2*Log[mu^2/Msq^2])/4 + 
   (3*CF*nq*nt*Tf^2*Log[mu^2/Msq^2])/4 - (9*CA*CF*Mst1^2*nt*Tf*Log[mu^2/Mst1^2])/(16*Msq^2) + 
   (3*CF*Mst1^2*nq*nt*Tf^2*Log[mu^2/Mst1^2])/(8*Msq^2) + (3*CF*Mst1^2*nt^2*Tf^2*Log[mu^2/Mst1^2])/(8*Msq^2) - 
   (9*CA*CF*Mst2^2*nt*Tf*Log[mu^2/Mst2^2])/(16*Msq^2) + (3*CF*Mst2^2*nq*nt*Tf^2*Log[mu^2/Mst2^2])/(8*Msq^2) + 
   (3*CF*Mst2^2*nt^2*Tf^2*Log[mu^2/Mst2^2])/(8*Msq^2) - (9*CA*CF*MeOS^2*nt*Tf*Log[mu^2/Mt^2])/(16*Msq^2) + 
   (9*CA*CF*Mt^2*nt*Tf*Log[mu^2/Mt^2])/(8*Msq^2) + (3*CF*MeOS^2*nq*nt*Tf^2*Log[mu^2/Mt^2])/(8*Msq^2) - 
   (3*CF*Mt^2*nq*nt*Tf^2*Log[mu^2/Mt^2])/(4*Msq^2) + (3*CF*MeOS^2*nt^2*Tf^2*Log[mu^2/Mt^2])/(8*Msq^2) - 
   (3*CF*Mt^2*nt^2*Tf^2*Log[mu^2/Mt^2])/(4*Msq^2)));
