* MFLWI.F by G. Passarino, November 1999 * i mu,al,be,alp,bep,i,j; * v pp,qp,k,kb,Qp,Qm,BQ,dv; * s A0; s RA0; s one; s [pm]; s Deta; s Qu,Qd; s p,m,W; s [Y-x2]; s M2,aux; s B0s,B0t; s Pfs,Pft; s tau,upsi; s Gs,Gt,Gw; s DSgs,DSgt; s M2_s,M2_t; s [1+Gw*Sgw]; s Sgs,Sgt,Sgw; s eta_s,eta_t; s [1-Gw*sigw]; s Neta_s,Deta_s; s Neta_t,Deta_t; s sigs,sigt,sigw; s aux1,aux2,aux3; s sw,fws,fwt,fww; s sigps,sigpt,sigpw; s Rm_s,Rm_t,[1+Rm_s]; s sela,selb,selc,seld; s [QpS-Spt],[BQS-Sps]; s B0s_MM,B0s_mm,B0s_mM; s B0t_MM,B0t_mm,B0t_mM; s B0w_MM,B0w_mm,B0w_mM; s [BQS+M2_s],[QpS+M2_t]; s up,tp,zp,zm,x1,T,U,sp; s B21s_MM,B21s_mm,B21s_Mm; s B21t_MM,B21t_mm,B21t_Mm; s B21w_MM,B21w_mm,B21w_Mm; s swvs1,swvs2,swsv1,swsv2,swss; s [em^2-y*s],[um^2-t],[dm^2-u]; s g,st,srt,y,s,x2,X,km,kp,z,Y,selm,set; s [Gw/g^2*(y*s)/sw*x2-1],[Gw/g^2*(y*s)/sw*Y-1]; s wvv0,wvv1,wvv2,wvv3,wvv4,wvs1,wvs2,wsv1,wsv2,wss; s [(1+Gw*Sgw)*(y*s)/sw*x2-1],[(1+Gw*Sgw)*(y*s)/sw*Y-1]; s [-x2*y*s+M2],[-Y*y*s+M2],[x2-M2/(y*s)],[Y-M2/(y*s)]; s QmS,Sws,Swt,[BQS+M2],[QpS+M2],[(qp+BQ)^2+em^2],[(Qm-k)^2+um^2], [(Qm-kb)^2+dm^2],[(qp+BQ)^2],[(Qm-k)^2],[(Qm-kb)^2]; s [(1+Gw*Sgw)*x2*y*s/sw+Gw*sigw-1],[(1+Gw*Sgw)*Y*y*s/sw+Gw*sigw-1]; s em,um,dm,tm,tmr,M,QpS,QpQ,BQS,BQQ,D_s,D_t,Sps,Spt,N_s,N_t,Ns_s,Ns_t, Pwt,Pws,Pps,Ppt; s z1,z2,z3,z4,z5,z6,z7,z8,z9,z10,z11,z12,z13,z14,z15,z16,z17,z18,z19,z20, z21,z22,z23; * f gd,gd5,gdp,gdm; f vv1,vs1,sv1,ss1; f dbt_vv,dbt_vs,dbt_sv,dbt_ss; f v,vb,u,ub,slv,sls,sqv,sqs,vv0,vs0,sv0,ss0,dbs_vv,dbs_vs,dbs_sv,dbs_ss; * cf ST; cf ST_l,ST_q; cf d1,d2,d3,d4,d1v; * .global * #procedure genesisnv(dummy) * g WI_dd= Qm(mu)*(d1(mu)+d2(mu)+d3(mu)+d4(mu))*selm; g WI_extra= Qm(mu)*(d1(mu)+d2(mu)+d3(mu)+d4(mu))*(1-selm)*set*BQS^2*QpS^2; * #call diagnv{dummy} #call tlvertices{dummy} #call prop{dummy} * id selm*D_s^-1= 0; id selm*D_t^-1= 0; id selm= 1; * id set= Pws*Pwt*Pps*Ppt*D_s*D_t; id srt^-1= 1/2; * id Ns_s= BQS+M^2-M^2/BQS*Sps; al Ns_t= QpS+M^2-M^2/QpS*Spt; id N_s= M^2*Sps-BQS*Sws+M^2*Sps^2/Pps; al N_t= M^2*Spt-QpS*Swt+M^2*Spt^2/Ppt; * id D_s*Pps= [BQS+M2]^2*(1-Sps/BQS); id D_t*Ppt= [QpS+M2]^2*(1-Spt/QpS); * id Pws= [BQS+M2]-Sws; al Pwt= [QpS+M2]-Swt; al Pps= [BQS+M2]-Sps; al Ppt= [QpS+M2]-Spt; * id Sws= BQS+M^2/BQS*Sps-M^2/M2_s*(BQS-Sps); al Swt= QpS+M^2/QpS*Spt-M^2/M2_t*(QpS-Spt); * #call dirac{dummy} * b i_,g,st,y,s,em,um,dm,u,ub,v,vb,M,Pws,Pwt,Pps,Ppt,D_s,D_t,Sws,Swt,Sps,Spt, [BQS+M2],[QpS+M2],[(qp+BQ)^2+em^2],[(Qm-k)^2+um^2],[(Qm-kb)^2+dm^2], Qu,Qd,M2_s,M2_t,[BQS-Sps],[QpS-Spt],ST_l,ST_q; print; * .store * save outputwi.sav WI_dd,WI_extra; * #endprocedure * #procedure genesisv(dummy) * g WIv_dd= d1v*selm; g WIv_extra= d1v*(1-selm)*set*BQS^2*QpS^2; * #call diagv{dummy} #call tlvertices{dummy} #call prop{dummy} * id selm*D_s^-1= 0; id selm*D_t^-1= 0; id selm= 1; * id set= Pws*Pwt*Pps*Ppt*D_s*D_t; id srt^-1= 1/2; * id Ns_s= BQS+M^2-M^2/BQS*Sps; al Ns_t= QpS+M^2-M^2/QpS*Spt; id N_s= M^2*Sps-BQS*Sws+M^2*Sps^2/Pps; al N_t= M^2*Spt-QpS*Swt+M^2*Spt^2/Ppt; * id D_s*Pps= [BQS+M2]^2*(1-Sps/BQS); id D_t*Ppt= [QpS+M2]^2*(1-Spt/QpS); * id Pws= [BQS+M2]-Sws; al Pwt= [QpS+M2]-Swt; al Pps= [BQS+M2]-Sps; al Ppt= [QpS+M2]-Spt; * id Sws= BQS+M^2/BQS*Sps-M^2/M2_s*(BQS-Sps); al Swt= QpS+M^2/QpS*Spt-M^2/M2_t*(QpS-Spt); * #call dirac{dummy} * b i_,g,st,y,s,em,um,dm,u,ub,v,vb,M,Pws,Pwt,Pps,Ppt,D_s,D_t,Sws,Swt,Sps,Spt, [BQS+M2],[QpS+M2],[(qp+BQ)^2+em^2],[(Qm-k)^2+um^2],[(Qm-kb)^2+dm^2], Qu,Qd,M2_s,M2_t,[BQS-Sps],[QpS-Spt],ST_l,ST_q, wvv0,wvv1,wvv2,wvv3,wvv4,wvs1,wvs2,wsv1,wsv2,wss; print; * .store * save outputwiv.sav WIv_dd,WIv_extra; * #endprocedure * #procedure diag_nv(dummy) * id d1(mu?)= (i_*g/2/srt)^2*g*st*( vb(i,pp)*slv(i,1,be)*v(i,qp)*ub(j,k)*sqv(j,1,al)*v(j,kb)* vv0(mu,alp,bep)*dbs_vv(al,alp)*dbt_vv(be,bep)+ vb(i,pp)*slv(i,1,be)*v(i,qp)*ub(j,k)*sqs(j,1)*v(j,kb)* vv0(mu,alp,bep)*dbs_vs(alp)*dbt_vv(be,bep)+ vb(i,pp)*slv(i,1,be)*v(i,qp)*ub(j,k)*sqv(j,1,al)*v(j,kb)* vs0(mu,bep)*dbs_sv(al)*dbt_vv(be,bep)+ vb(i,pp)*slv(i,1,be)*v(i,qp)*ub(j,k)*sqs(j,1)*v(j,kb)* vs0(mu,bep)*dbs_ss*dbt_vv(be,bep)+ * vb(i,pp)*slv(i,1,be)*v(i,qp)*ub(j,k)*sqv(j,1,al)*v(j,kb)* sv0(mu,alp)*dbs_vv(al,alp)*dbt_vs(be)+ vb(i,pp)*slv(i,1,be)*v(i,qp)*ub(j,k)*sqs(j,1)*v(j,kb)* sv0(mu,alp)*dbs_vs(alp)*dbt_vs(be)+ vb(i,pp)*slv(i,1,be)*v(i,qp)*ub(j,k)*sqv(j,1,al)*v(j,kb)* ss0(mu)*dbs_sv(al)*dbt_vs(be)+ vb(i,pp)*slv(i,1,be)*v(i,qp)*ub(j,k)*sqs(j,1)*v(j,kb)* ss0(mu)*dbs_ss*dbt_vs(be)+ * vb(i,pp)*sls(i,1)*v(i,qp)*ub(j,k)*sqv(j,1,al)*v(j,kb)* vv0(mu,alp,bep)*dbs_vv(al,alp)*dbt_sv(bep)+ vb(i,pp)*sls(i,1)*v(i,qp)*ub(j,k)*sqs(j,1)*v(j,kb)* vv0(mu,alp,bep)*dbs_vs(alp)*dbt_sv(bep)+ vb(i,pp)*sls(i,1)*v(i,qp)*ub(j,k)*sqv(j,1,al)*v(j,kb)* vs0(mu,bep)*dbs_sv(al)*dbt_sv(bep)+ vb(i,pp)*sls(i,1)*v(i,qp)*ub(j,k)*sqs(j,1)*v(j,kb)* vs0(mu,bep)*dbs_ss*dbt_sv(bep)+ * vb(i,pp)*sls(i,1)*v(i,qp)*ub(j,k)*sqv(j,1,al)*v(j,kb)* sv0(mu,alp)*dbs_vv(al,alp)*dbt_ss+ vb(i,pp)*sls(i,1)*v(i,qp)*ub(j,k)*sqs(j,1)*v(j,kb)* sv0(mu,alp)*dbs_vs(alp)*dbt_ss+ vb(i,pp)*sls(i,1)*v(i,qp)*ub(j,k)*sqv(j,1,al)*v(j,kb)* ss0(mu)*dbs_sv(al)*dbt_ss+ vb(i,pp)*sls(i,1)*v(i,qp)*ub(j,k)*sqs(j,1)*v(j,kb)* ss0(mu)*dbs_ss*dbt_ss); * id d2(mu?)= (i_*g/2/srt)^2*(-i_*g*st)*( vb(i,pp)*slv(i,2,mu,be)*v(i,qp)*ub(j,k)*sqv(j,2,al)*v(j,kb)* dbs_vv(be,al)+ vb(i,pp)*slv(i,2,mu,be)*v(i,qp)*ub(j,k)*sqs(j,2)*v(j,kb)* dbs_vs(be)+ vb(i,pp)*sls(i,2,mu)*v(i,qp)*ub(j,k)*sqv(j,2,al)*v(j,kb)* dbs_sv(al)+ vb(i,pp)*sls(i,2,mu)*v(i,qp)*ub(j,k)*sqs(j,2)*v(j,kb)* dbs_ss); * id d3(mu?)= (i_*g/2/srt)^2*(i_*g*Qu*st)*( vb(i,pp)*slv(i,3,be)*v(i,qp)*ub(j,k)*sqv(j,3,mu,al)*v(j,kb)* dbt_vv(be,al)+ vb(i,pp)*slv(i,3,be)*v(i,qp)*ub(j,k)*sqs(j,3,mu)*v(j,kb)* dbt_vs(be)+ vb(i,pp)*sls(i,3)*v(i,qp)*ub(j,k)*sqv(j,3,mu,al)*v(j,kb)* dbt_sv(al)+ vb(i,pp)*sls(i,3)*v(i,qp)*ub(j,k)*sqs(j,3,mu)*v(j,kb)* dbt_ss); * id d4(mu?)= (i_*g/2/srt)^2*(i_*g*Qd*st)*( vb(i,pp)*slv(i,4,be)*v(i,qp)*ub(j,k)*sqv(j,4,mu,al)*v(j,kb)* dbt_vv(be,al)+ vb(i,pp)*slv(i,4,be)*v(i,qp)*ub(j,k)*sqs(j,4,mu)*v(j,kb)* dbt_vs(be)+ vb(i,pp)*sls(i,4)*v(i,qp)*ub(j,k)*sqv(j,4,mu,al)*v(j,kb)* dbt_sv(al)+ vb(i,pp)*sls(i,4)*v(i,qp)*ub(j,k)*sqs(j,4,mu)*v(j,kb)* dbt_ss); * sum al,alp,be,bep; * #endprocedure * #procedure diagv(dummy) * id d1v= (i_*g/2/srt)^2*g^3*st*( vb(i,pp)*slv(i,1,be)*v(i,qp)*ub(j,k)*sqv(j,1,al)*v(j,kb)* vv1(alp,bep)*dbs_vv(al,alp)*dbt_vv(be,bep)+ vb(i,pp)*slv(i,1,be)*v(i,qp)*ub(j,k)*sqs(j,1)*v(j,kb)* vv1(alp,bep)*dbs_vs(alp)*dbt_vv(be,bep)+ vb(i,pp)*slv(i,1,be)*v(i,qp)*ub(j,k)*sqv(j,1,al)*v(j,kb)* vs1(bep)*dbs_sv(al)*dbt_vv(be,bep)+ vb(i,pp)*slv(i,1,be)*v(i,qp)*ub(j,k)*sqs(j,1)*v(j,kb)* vs1(bep)*dbs_ss*dbt_vv(be,bep)+ * vb(i,pp)*slv(i,1,be)*v(i,qp)*ub(j,k)*sqv(j,1,al)*v(j,kb)* sv1(alp)*dbs_vv(al,alp)*dbt_vs(be)+ vb(i,pp)*slv(i,1,be)*v(i,qp)*ub(j,k)*sqs(j,1)*v(j,kb)* sv1(alp)*dbs_vs(alp)*dbt_vs(be)+ vb(i,pp)*slv(i,1,be)*v(i,qp)*ub(j,k)*sqv(j,1,al)*v(j,kb)* ss1*dbs_sv(al)*dbt_vs(be)+ vb(i,pp)*slv(i,1,be)*v(i,qp)*ub(j,k)*sqs(j,1)*v(j,kb)* ss1*dbs_ss*dbt_vs(be)+ * vb(i,pp)*sls(i,1)*v(i,qp)*ub(j,k)*sqv(j,1,al)*v(j,kb)* vv1(alp,bep)*dbs_vv(al,alp)*dbt_sv(bep)+ vb(i,pp)*sls(i,1)*v(i,qp)*ub(j,k)*sqs(j,1)*v(j,kb)* vv1(alp,bep)*dbs_vs(alp)*dbt_sv(bep)+ vb(i,pp)*sls(i,1)*v(i,qp)*ub(j,k)*sqv(j,1,al)*v(j,kb)* vs1(bep)*dbs_sv(al)*dbt_sv(bep)+ vb(i,pp)*sls(i,1)*v(i,qp)*ub(j,k)*sqs(j,1)*v(j,kb)* vs1(bep)*dbs_ss*dbt_sv(bep)+ * vb(i,pp)*sls(i,1)*v(i,qp)*ub(j,k)*sqv(j,1,al)*v(j,kb)* sv1(alp)*dbs_vv(al,alp)*dbt_ss+ vb(i,pp)*sls(i,1)*v(i,qp)*ub(j,k)*sqs(j,1)*v(j,kb)* sv1(alp)*dbs_vs(alp)*dbt_ss+ vb(i,pp)*sls(i,1)*v(i,qp)*ub(j,k)*sqv(j,1,al)*v(j,kb)* ss1*dbs_sv(al)*dbt_ss+ vb(i,pp)*sls(i,1)*v(i,qp)*ub(j,k)*sqs(j,1)*v(j,kb)* ss1*dbs_ss*dbt_ss); * sum al,alp,be,bep; * #endprocedure * #procedure tlvertices(dummy) * id slv(i?,3,mu?)= gd(i,mu)*gdp(i); al sls(i?,3)= em/M*gdp(i); id sqv(i?,3,mu?,al?)= gd(i,mu)*(i_*(gd(i,Qm)-gd(i,k))+um)/[(Qm-k)^2+um^2]* gd(i,al)*gdp(i); al sqs(i?,3,mu?)= gd(i,mu)*(i_*(gd(i,Qm)-gd(i,k))+um)/[(Qm-k)^2+um^2]* (um/M*gdp(i)-dm/M*gdm(i)); * id slv(i?,4,mu?)= gd(i,mu)*gdp(i); al sls(i?,4)= em/M*gdp(i); id sqv(i?,4,mu?,al?)= gd(i,al)*gdp(i)*(-i_*(gd(i,Qm)-gd(i,kb))+dm)/ [(Qm-kb)^2+dm^2]*gd(i,mu); al sqs(i?,4,mu?)= (um/M*gdp(i)-dm/M*gdm(i))*(-i_*(gd(i,Qm)-gd(i,kb))+dm)/ [(Qm-kb)^2+dm^2]*gd(i,mu); * id slv(i?,1,mu?)= gd(i,mu)*gdp(i); id sls(i?,1)= em/M*gdp(i); * id sqv(i?,1,mu?)= gd(i,mu)*gdp(i); al sqs(i?,1)= 1/2*((um-dm)/M*(gdp(i)+gdm(i))+(um+dm)/M*(gdp(i)-gdm(i))); * id vv0(mu?,al?,be?)= d_(mu,be)*(Qm(al)-Qp(al))+ d_(al,be)*(2*Qp(mu)+Qm(mu))- d_(mu,al)*(2*Qm(be)+Qp(be)); al vs0(mu?,al?)= i_*M*d_(mu,al); al sv0(mu?,al?)= -i_*M*d_(mu,al); al ss0(mu?)= Qp(mu)+BQ(mu); * id slv(i?,2,mu?,be?)= gd(i,mu)*(i_*(gd(i,qp)+gd(i,BQ))+em)/[(qp+BQ)^2+em^2]* gd(i,be)*gdp(i); al sls(i?,2,mu?)= gd(i,mu)*(i_*(gd(i,qp)+gd(i,BQ))+em)/[(qp+BQ)^2+em^2]* em/M*gdp(i); al sqv(i?,2,mu?)= gd(i,mu)*gdp(i); al sqs(i?,2)= 1/2*((um-dm)/M*(gdp(i)+gdm(i))+(um+dm)/M*(gdp(i)-gdm(i))); * id vv1(al?,be?)= wvv0*d_(al,be)+wvv1*Qp(al)*Qm(be)+wvv2*Qp(be)*Qm(al)+ wvv3*Qm(al)*Qm(be)+wvv4*Qp(al)*Qp(be); al vs1(be?)= wvs1*Qp(be)+wvs2*Qm(be); al sv1(al?)= wsv1*Qp(al)+wsv2*Qm(al); al ss1= wss; * #endprocedure * #procedure prop(dummy) * id dbt_vv(al?,be?)= 1/Pwt*(d_(al,be)+N_t/(QpQ*D_t)*Qp(al)*Qp(be)); al dbt_vs(al?)= i_*M*Spt/QpS*Qp(al)/Ppt/D_t; al dbt_sv(al?)= -i_*M*Spt/QpS*Qp(al)/Ppt/D_t; al dbt_ss= Ns_t/D_t/Ppt; * id dbs_vv(al?,be?)= 1/Pws*(d_(al,be)+N_s/(BQQ*D_s)*BQ(al)*BQ(be)); al dbs_vs(al?)= i_*M*Sps/BQS*BQ(al)/Pps/D_s; al dbs_sv(al?)= -i_*M*Sps/BQS*BQ(al)/Pps/D_s; al dbs_ss= Ns_s/D_s/Pps; * #endprocedure * #procedure dirac(dummy) * id gd(i?,Qm)*gd(i?,Qm)= QmS; al gd(i,Qp)*gd(i,BQ)= QpS+gd(i,Qp)*gd(i,Qm); * id ub(j,k)*gd(j,Qm)*gd(j,k)= -i_*um*ub(j,k)*gd(j,Qm)+2*k.Qm*ub(j,k); id gd(j,kb)*gd(j,Qm)*v(j,kb)= i_*dm*gd(j,Qm)*v(j,kb)+2*kb.Qm*v(j,kb); id ub(j,k)*gd(j,k)= i_*um*ub(j,k); * id gd(i,BQ)*gd(i,BQ)= BQS; id gd(i,Qm)*gd(i,BQ)= BQS-gd(i,Qp)*gd(i,BQ); id vb(i,pp)*gd(i,Qp)*gd(i,BQ)*gdp(i)*v(i,qp)= -i_*em*vb(i,pp)*gd(i,BQ)*gdp(i)*v(i,qp)- 2*qp.BQ*vb(i,pp)*gdp(i)*v(i,qp); * id vb(i,pp)*gd(i,Qp)*gdp(i)*v(i,qp)= -i_*em*vb(i,pp)*gdp(i)*v(i,qp); id ub(j,k)*gd(j,BQ)*gdp(j)*v(j,kb)= i_*ub(j,k)*(um*gdp(j)-dm*gdm(j))*v(j,kb); id ub(j,k)*gd(j,Qm)*gdp(j)*v(j,kb)= -ub(j,k)*gd(j,Qp)*gdp(j)*v(j,kb)+ i_*ub(j,k)*(um*gdp(j)-dm*gdm(j))*v(j,kb); id vb(i,pp)*gd(i,BQ)*gdp(i)*v(i,qp)= vb(i,pp)*gd(i,Qm)*gdp(i)*v(i,qp)- i_*em*vb(i,pp)*gdp(i)*v(i,qp); repeat; id gd(i,qp)*gdp(i)= gdm(i)*gd(i,qp); id gd(i,qp)*gd(i,mu?)= 2*qp(mu)-gd(i,mu)*gd(i,qp); endrepeat; id gd(i,qp)*v(i,qp)= 0; * id ub(j,k)*gd(j,Qp)*gdp(j)*gd(j,Qm)*v(j,kb)= -QmS*ub(j,k)*gdm(j)*v(j,kb)+ i_*um*ub(j,k)*gd(j,Qm)*gdm(j)*v(j,kb)+ i_*dm*ub(j,k)*gd(j,Qm)*gdp(j)*v(j,kb)+ 2*kb.Qm*ub(j,k)*gdm(j)*v(j,kb); id ub(j,k)*gd(j,Qp)*gdm(j)*gd(j,Qm)*v(j,kb)= -QmS*ub(j,k)*gdp(j)*v(j,kb)+ i_*um*ub(j,k)*gd(j,Qm)*gdp(j)*v(j,kb)+ i_*dm*ub(j,k)*gd(j,Qm)*gdm(j)*v(j,kb)+ 2*kb.Qm*ub(j,k)*gdp(j)*v(j,kb); * id ub(j,k)*gd(j,Qm)*gd(j,Qp)*gdp(j)*v(j,kb)= -QmS*ub(j,k)*gdp(j)*v(j,kb)- i_*um*ub(j,k)*gd(j,Qm)*gdp(j)*v(j,kb)+ 2*k.Qm*ub(j,k)*gdp(j)*v(j,kb)- i_*dm*ub(j,k)*gd(j,Qm)*gdm(j)*v(j,kb); id ub(j,k)*gd(j,Qm)*gd(j,Qp)*gdm(j)*v(j,kb)= -QmS*ub(j,k)*gdm(j)*v(j,kb)- i_*um*ub(j,k)*gd(j,Qm)*gdm(j)*v(j,kb)+ 2*k.Qm*ub(j,k)*gdm(j)*v(j,kb)- i_*dm*ub(j,k)*gd(j,Qm)*gdp(j)*v(j,kb); * id ub(j,k)*gd(j,Qp)*gdp(j)*v(j,kb)= -ub(j,k)*gd(j,Qm)*gdp(j)*v(j,kb)+ i_*um*ub(j,k)*gdp(j)*v(j,kb)- i_*dm*ub(j,k)*gdm(j)*v(j,kb); id ub(j,k)*gd(j,Qp)*gdm(j)*v(j,kb)= -ub(j,k)*gd(j,Qm)*gdm(j)*v(j,kb)+ i_*um*ub(j,k)*gdm(j)*v(j,kb)- i_*dm*ub(j,k)*gdp(j)*v(j,kb); * id vb(i,pp)*gd(i,Qp)*gd(i,Qm)*gd(i,mu?)*gdp(i)*v(i,qp)= -i_*em*vb(i,pp)*gd(i,Qm)*gd(i,mu)*gdp(i)*v(i,qp)- 2*qp.Qm*vb(i,pp)*gd(i,mu)*gdp(i)*v(i,qp)+ 2*qp(mu)*vb(i,pp)*gd(i,Qm)*gdp(i)*v(i,qp); * id gd(i?,Qm)*gd(i?,Qm)= QmS; al gd(i,Qp)*gd(i,BQ)= QpS+gd(i,Qp)*gd(i,Qm); * id qp.BQ= 1/2*([(qp+BQ)^2+em^2]-BQS-em^2); id Qm.Qm= QmS; al Qp.Qp= QpS; * id k.kb= 1/2*(BQS+um^2+dm^2); id k.Qp= k.BQ-k.Qm; id k.BQ= 1/2*(BQS+dm^2-um^2); id k.Qm= -1/2*([(Qm-k)^2+um^2]-QmS); id kb.Qm= 1/2*(QmS-[(Qm-kb)^2+dm^2]); * id Qm.BQ= BQS-Qp.BQ; al Qm.Qp= BQ.Qp-QpS; al QmS= BQS+QpS-2*BQ.Qp; * id BQQ= BQS^2; al BQQ^-1= 1/BQS^2; id QpQ= QpS^2; al QpQ^-1= 1/QpS^2; * id BQS= [BQS+M2]-M^2; id QpS= [QpS+M2]-M^2; * id Sps= -[BQS-Sps]+BQS; id Spt= -[QpS-Spt]+QpS; * id BQS= [BQS+M2]-M^2; id QpS= [QpS+M2]-M^2; * id gdp(i?)*gd(i?,Qm)= gd(i,Qm)*gdm(i); id gdm(i?)*gd(i?,Qm)= gd(i,Qm)*gdp(i); * id ub(j,k)*gd(j,mu?)*gd(j,Qm)= -ub(j,k)*gd(j,Qm)*gd(j,mu)+2*Qm(mu)*ub(j,k); * id vb(i,pp)*gd(i,Qp)*gd(i,Qm)*gd(i,mu?)*gdp(i)*v(i,qp)= -i_*em*vb(i,pp)*gd(i,Qm)*gd(i,mu)*gdp(i)*v(i,qp)- 2*qp.Qm*vb(i,pp)*gd(i,mu)*gdp(i)*v(i,qp)+ 2*qp(mu)*vb(i,pp)*gd(i,Qm)*gdp(i)*v(i,qp); * id vb(i,pp)*gd(i,mu?)*gdp(i)*v(i,qp)= ST_l(mu,p); id vb(i,pp)*gd(i,mu?)*gdm(i)*v(i,qp)= ST_l(mu,m); id ub(j,k)*gd(j,mu?)*gdp(j)*v(j,kb)= ST_q(mu,p); id ub(j,k)*gd(j,mu?)*gdm(j)*v(j,kb)= ST_q(mu,m); * id vb(i,pp)*gd(i,Qm)*gdp(i)*v(i,qp)= ST_l(Qm,p); id vb(i,pp)*gd(i,Qm)*gdm(i)*v(i,qp)= ST_l(Qm,m); id vb(i,pp)*gdp(i)*v(i,qp)= ST_l(1,p); id vb(i,pp)*gdm(i)*v(i,qp)= ST_l(1,m); id ub(j,k)*gd(j,Qm)*gdp(j)*v(j,kb)= ST_q(Qm,p); id ub(j,k)*gd(j,Qm)*gdm(j)*v(j,kb)= ST_q(Qm,m); id ub(j,k)*gdp(j)*v(j,kb)= ST_q(1,p); id ub(j,k)*gdm(j)*v(j,kb)= ST_q(1,m); * #endprocedure * #procedure massless(dummy) * load outputwi.sav; * g WI_dd_m= WI_dd; g Wi_extra_m= WI_extra; * id em= 0; al um= 0; al dm= 0; * id [(qp+BQ)^2+em^2]^-1= [(qp+BQ)^2]^-1; al [(Qm-k)^2+um^2]^-1= [(Qm-k)^2]^-1; al [(Qm-kb)^2+dm^2]^-1= [(Qm-kb)^2]^-1; * b i_,g,st,y,s,em,um,dm,u,ub,v,vb,M,Pws,Pwt,Pps,Ppt,D_s,D_t,Sws,Swt,Sps,Spt, [BQS+M2],[QpS+M2],[(qp+BQ)^2],[(Qm-k)^2],[(Qm-kb)^2], Qu,Qd,M2_s,M2_t,[BQS-Sps],[QpS-Spt],ST_l,ST_q; print; * .store * #endprocedure * #procedure mfixedw(dummy) * g FW_m_WI= WI_dd_m*Pws*Pwt; * id Pws= [BQS+M2]-i_*M*W; al Pwt= [QpS+M2]-i_*M*W; * id [BQS+M2]= BQS+M^2; al [QpS+M2]= QpS+M^2; * id Qp.BQ= QpS+Qp.Qm; al Qm.k= 1/2*(QmS+tau*y*s); id BQS= QpS+QmS+2*Qp.Qm; id QpS= -2*pp.qp; id Qp.Qm= pp.Qm-qp.Qm; * id pp.Qm= -1/2*(QmS+y*s); id Qm.qp= 1/2*(QmS+km); al pp.qp= 1/2*kp; * id km= z*y*s; id kp= y*s*(-1-X+x2-z); id QmS= X*y*s; * id [(qp+BQ)^2]^-1= -1/(y*s); al [(Qm-k)^2]^-1= -1/(tau*y*s); al [(Qm-kb)^2]^-1= -1/(upsi*y*s); * id Qu= Qd+1; * b i_,g,st,y,s,em,um,dm,u,ub,v,vb,M,Pws,Pwt,Pps,Ppt,D_s,D_t,Sws,Swt,Sps,Spt, [BQS+M2],[QpS+M2],[(qp+BQ)^2],[(Qm-k)^2],[(Qm-kb)^2], Qu,Qd,M2_s,M2_t,[BQS-Sps],[QpS-Spt],ST_l,ST_q,tau,upsi; print; * .store * #endprocedure * #procedure evolution(dummy) * load outputwi.sav; * g NWI= (WI_dd+WI_extra/BQS/QpS/[BQS+M2]^2/[QpS+M2]^2/[BQS-Sps]/[QpS-Spt]/ Pws/Pwt)/g/st; * id [BQS-Sps]= BQS*M2_s/M^2*Pws/[BQS+M2_s]; id [QpS-Spt]= QpS*M2_t/M^2*Pwt/[QpS+M2_t]; * id [BQS+M2]= BQS+M^2; al [QpS+M2]= QpS+M^2; * ratio,BQS,[BQS+M2],M2; ratio,QpS,[QpS+M2],M2; id M2= M^2; al M2^-1= 1/M^2; * ratio,BQS,[BQS+M2_s],M2_s; ratio,QpS,[QpS+M2_t],M2_t; * id BQS= QpS+QmS+2*Qp.Qm; id BQ.dv?= Qp.dv+Qm.dv; id Qp.Qm= pp.Qm-qp.Qm; id Qp.Qp= QpS; id QpS= -em^2-2*pp.qp; * id pp.Qm= 1/2*(em^2-QmS-y*s); al k.kb= 1/2*(-sp+um^2+dm^2); al Qm.k= 1/2*(QmS-um^2+T); al pp.kb= 1/2*(-em^2-dm^2+tp); al Qm.kb= 1/2*(QmS-dm^2+U); al pp.k= 1/2*(-em^2-um^2+up); al pp.qp= 1/2*(kp-em^2); al Qm.qp= 1/2*(QmS+km); al kb.qp= 1/2*(-dm^2+zp); al k.qp= 1/2*(um^2+zm); * id T= tau*y*s; al sp= x2*y*s; al km= z*y*s; al zm= (x2-x1)*y*s; id QmS= X*y*s; * id kp= em^2 * ( 1 ) + y*s * ( - 1 - X + x2 - z ); al U= em^2 * ( 1 ) + um^2 * ( 1 ) + dm^2 * ( 1 ) + y*s * ( - 1 - tau - 2*X - z ); al up= em^2 * ( 1 ) + um^2 * ( 2 ) + dm^2 * ( 1 ) + y*s * ( - tau - X - x1 ); al tp= um^2 * ( - 1 ) + y*s * ( tau + X + x1 - x2 + z ); al zp= um^2 * ( - 1 ) + dm^2 * ( 1 ) + y*s * ( - 1 + x1 ); * id [(qp+BQ)^2+em^2]^-1= 1/[em^2-y*s]; al [(Qm-k)^2+um^2]^-1= 1/[um^2-t]; al [(Qm-kb)^2+dm^2]^-1= 1/[dm^2-u]; id BQS^-1= -1/(x2*y*s); *al QpS^-1= 1/(Y*y*s); al QpS^-1= -1/(Y*y*s); * id Qu= Qd+1; * id em^2*[em^2-y*s]^-1= y*s*[em^2-y*s]^-1+1; * id [BQS+M2]^-1= [-x2*y*s+M2]^-1; al [QpS+M2]^-1= [-Y*y*s+M2]^-1; id z= 1/(y*s)*(em^2+(-1-X+x2-Y)*y*s); * id [-x2*y*s+M2]^-1= -1/(y*s)/[x2-M2/(y*s)]; al [-Y*y*s+M2]^-1= -1/(y*s)/[Y-M2/(y*s)]; * ratio,x2,[x2-M2/(y*s)],aux; id aux= -M^2/(y*s); al aux^-1= -y*s/M^2; ratio,Y,[Y-M2/(y*s)],aux; id aux= -M^2/(y*s); al aux^-1= -y*s/M^2; * id ST_l(Qm,p)*ST_q(1,p)*um= ST_l(Qm,p)*ST_q(1,m)*dm+ST(Qm,p,1,[pm]); id ST_l(1,p)*ST_q(1,p)*um= ST_l(1,p)*ST_q(1,m)*dm+ST(1,p,1,[pm]); * id M2_s^-1= (1+Rm_s)/M^2; al M2_t^-1= (1+Rm_t)/M^2; * repeat; id M^2*[x2-M2/(y*s)]^-1= -y*s*(1-x2/[x2-M2/(y*s)]); al M^2*[Y-M2/(y*s)]^-1= -y*s*(1-Y/[Y-M2/(y*s)]); endrepeat; * b i_,g,st,y,s,em,um,dm,M,u,ub,v,vb,Pws,Pwt,Pps,Ppt,D_s,D_t,Sws,Swt,Sps,Spt, [BQS+M2],[QpS+M2],[(qp+BQ)^2+em^2],[(Qm-k)^2+um^2],[(Qm-kb)^2+dm^2], Qu,Qd,[BQS-Sps],[QpS-Spt],ST_l,ST_q,ST,[em^2-y*s],[um^2-t],[dm^2-u], [BQS+M2_s],[QpS+M2_t],[x2-M2/(y*s)],[Y-M2/(y*s)],Rm_s,Rm_t,[1+Rm_s]; print; * .store * #endprocedure * #procedure evolutionv(dummy) * load outputwiv.sav; * g NWIv= (WIv_dd+WIv_extra/BQS/QpS/[BQS+M2]^2/[QpS+M2]^2/[BQS-Sps]/[QpS-Spt]/ Pws/Pwt)/g/st; * id [BQS-Sps]= BQS*M2_s/M^2*Pws/[BQS+M2_s]; id [QpS-Spt]= QpS*M2_t/M^2*Pwt/[QpS+M2_t]; * id [BQS+M2]= BQS+M^2; al [QpS+M2]= QpS+M^2; * ratio,BQS,[BQS+M2],M2; ratio,QpS,[QpS+M2],M2; id M2= M^2; al M2^-1= 1/M^2; * ratio,BQS,[BQS+M2_s],M2_s; ratio,QpS,[QpS+M2_t],M2_t; * id BQS= QpS+QmS+2*Qp.Qm; id BQ.dv?= Qp.dv+Qm.dv; id Qp.Qm= pp.Qm-qp.Qm; id Qp.Qp= QpS; id QpS= -em^2-2*pp.qp; * id pp.Qm= 1/2*(em^2-QmS-y*s); al k.kb= 1/2*(-sp+um^2+dm^2); al Qm.k= 1/2*(QmS-um^2+T); al pp.kb= 1/2*(-em^2-dm^2+tp); al Qm.kb= 1/2*(QmS-dm^2+U); al pp.k= 1/2*(-em^2-um^2+up); al pp.qp= 1/2*(kp-em^2); al Qm.qp= 1/2*(QmS+km); al kb.qp= 1/2*(-dm^2+zp); al k.qp= 1/2*(um^2+zm); * id T= tau*y*s; al sp= x2*y*s; al km= z*y*s; al zm= (x2-x1)*y*s; id QmS= X*y*s; * id kp= em^2 * ( 1 ) + y*s * ( - 1 - X + x2 - z ); al U= em^2 * ( 1 ) + um^2 * ( 1 ) + dm^2 * ( 1 ) + y*s * ( - 1 - tau - 2*X - z ); al up= em^2 * ( 1 ) + um^2 * ( 2 ) + dm^2 * ( 1 ) + y*s * ( - tau - X - x1 ); al tp= um^2 * ( - 1 ) + y*s * ( tau + X + x1 - x2 + z ); al zp= um^2 * ( - 1 ) + dm^2 * ( 1 ) + y*s * ( - 1 + x1 ); * id [(qp+BQ)^2+em^2]^-1= 1/[em^2-y*s]; al [(Qm-k)^2+um^2]^-1= 1/[um^2-t]; al [(Qm-kb)^2+dm^2]^-1= 1/[dm^2-u]; id BQS^-1= -1/(x2*y*s); *al QpS^-1= 1/(Y*y*s); al QpS^-1= -1/(Y*y*s); * id Qu= Qd+1; * id em^2*[em^2-y*s]^-1= y*s*[em^2-y*s]^-1+1; * id [BQS+M2]^-1= [-x2*y*s+M2]^-1; al [QpS+M2]^-1= [-Y*y*s+M2]^-1; id z= 1/(y*s)*(em^2+(-1-X+x2-Y)*y*s); * id [-x2*y*s+M2]^-1= -1/(y*s)/[x2-M2/(y*s)]; al [-Y*y*s+M2]^-1= -1/(y*s)/[Y-M2/(y*s)]; * ratio,x2,[x2-M2/(y*s)],aux; id aux= -M^2/(y*s); al aux^-1= -y*s/M^2; ratio,Y,[Y-M2/(y*s)],aux; id aux= -M^2/(y*s); al aux^-1= -y*s/M^2; *-- id ST_l(Qm,p)*ST_q(1,p)*um= ST_l(Qm,p)*ST_q(1,m)*dm+ST(Qm,p,1,[pm]); id ST_l(1,p)*ST_q(1,p)*um= ST_l(1,p)*ST_q(1,m)*dm+ST(1,p,1,[pm]); * id M2_s^-1= (1+Rm_s)/M^2; al M2_t^-1= (1+Rm_t)/M^2; * repeat; id M^2*[x2-M2/(y*s)]^-1= -y*s*(1-x2/[x2-M2/(y*s)]); al M^2*[Y-M2/(y*s)]^-1= -y*s*(1-Y/[Y-M2/(y*s)]); endrepeat; * id wvs1= 0; al wvs2= 0; al wsv1= 0; al wsv2= 0; al wss= 0; * id wvv2= wvv1; al wvv3= wvv1; * b i_,g,st,y,s,em,um,dm,M,u,ub,v,vb,Pws,Pwt,Pps,Ppt,D_s,D_t,Sws,Swt,Sps,Spt, [BQS+M2],[QpS+M2],[(qp+BQ)^2+em^2],[(Qm-k)^2+um^2],[(Qm-kb)^2+dm^2], Qu,Qd,[BQS-Sps],[QpS-Spt],ST_l,ST_q,ST,[em^2-y*s],[um^2-t],[dm^2-u], [BQS+M2_s],[QpS+M2_t],[x2-M2/(y*s)],[Y-M2/(y*s)],Rm_s,Rm_t,[1+Rm_s], wvv0,wvv1,wvv2,wvv3,wvv4,wvs1,wvs2,wsv1,wsv2,wss; print; * .store * #endprocedure * #procedure Mfixedw(dummy) * g NNWI= NWI*Pws*Pwt; * id Pws= -x2*y*s+M^2-i_*M*W; al Pwt= -Y*y*s+M^2-i_*M*W; * ratio,x2,[x2-M2/(y*s)],aux; id aux= -M^2/(y*s); al aux^-1= -y*s/M^2; ratio,Y,[Y-M2/(y*s)],aux; id aux= -M^2/(y*s); al aux^-1= -y*s/M^2; * id M2_s= M^2; id M2_t= M^2; id M2_s^-1= 1/M^2; id M2_t^-1= 1/M^2; id Rm_s= 0; al Rm_t= 0; id [1+Rm_s]^-1= 1; * id Qu= Qd+1; * b i_,g,st,y,s,em,um,dm,u,ub,v,vb,Pws,Pwt,Pps,Ppt,D_s,D_t,Sws,Swt,Sps,Spt, [BQS+M2],[QpS+M2],[(qp+BQ)^2+em^2],[(Qm-k)^2+um^2],[(Qm-kb)^2+dm^2], Qu,Qd,[BQS-Sps],[QpS-Spt],ST_l,ST_q,ST,[em^2-y*s],[um^2-t],[dm^2-u], [BQS+M2_s],[QpS+M2_t],[x2-M2/(y*s)],[Y-M2/(y*s)],Rm_s,Rm_t; print; * .store * #endprocedure * #procedure renorm(dummy) * g RWI= NWI*Pfs*Pft/aux; * id Pws^-1= -Gs/g^2/Pfs; al Pwt^-1= -Gt/g^2/Pft; * id Pfs= x2*y*s-sw+(1-Gs/Gw)*sw+Gs*fws; id Pft= Y*y*s-sw+(1-Gt/Gw)*sw+Gt*fwt; * id fws= 0; al fwt= 0; * id Rm_s= g^2/Gs-1; al Rm_t= g^2/Gt-1; id M^2= g^2/Gw*sw; al M^-2= Gw/g^2/sw; id [x2-M2/(y*s)]^-1= Gw/g^2*y*s/sw/[Gw/g^2*(y*s)/sw*x2-1]; id [Y-M2/(y*s)]^-1= Gw/g^2*y*s/sw/[Gw/g^2*(y*s)/sw*Y-1]; * id g^-2= 1/2*(1/Gt+Sgt+1/Gs+Sgs); * id Gs*Gt*Gw= aux1; id Gs*Gt= aux2; id Gt= Gs+Gs*Gt*(Sgt-Sgs); id aux1= Gs*Gt*Gw; id aux2= Gs*Gt; * id Gs*Gt^-1= 1+Gs*(Sgs-Sgt); * id Gs*Gt*Gw= aux1; id Gt*Gw= Gs*Gw+Gs*Gt*Gw*(Sgt-Sgs); id aux1= Gs*Gt*Gw; * id [Gw/g^2*(y*s)/sw*x2-1]^-1= 1/[(1+Gw*Sgw)*(y*s)/sw*x2-1]; id [Gw/g^2*(y*s)/sw*Y-1]^-1= 1/[(1+Gw*Sgw)*(y*s)/sw*Y-1]; * id aux^-1= 1/(Gs*Gt); id Gw*Gt^-1= 1+Gw*(Sgw-Sgt); * id Sgs= 4/9-4/3*B0s; al Sgt= 4/9-4/3*B0t; * id Gw*Sgw*[(1+Gw*Sgw)*(y*s)/sw*Y-1]^-1= -[(1+Gw*Sgw)*(y*s)/sw*Y-1]^-1+sw/(y*s*Y)*(1+[(1+Gw*Sgw)*(y*s)/sw*Y-1]^-1); al Gw*Sgw*[(1+Gw*Sgw)*(y*s)/sw*x2-1]^-1= -[(1+Gw*Sgw)*(y*s)/sw*x2-1]^-1+sw/(y*s*x2)*(1+[(1+Gw*Sgw)*(y*s)/sw*x2-1]^-1); * b i_,g,st,y,s,em,um,dm,u,ub,v,vb,Pws,Pwt,Pps,Ppt,D_s,D_t,Sws,Swt,Sps,Spt, [BQS+M2],[QpS+M2],[(qp+BQ)^2+em^2],[(Qm-k)^2+um^2],[(Qm-kb)^2+dm^2], Qu,Qd,[BQS-Sps],[QpS-Spt],ST_l,ST_q,ST,[em^2-y*s],[um^2-t],[dm^2-u], [BQS+M2_s],[QpS+M2_t],[x2-M2/(y*s)],[Y-M2/(y*s)],Rm_s,Rm_t,Gs,Gt,Gw,fws,fwt, fww,Sgs,Sgt,Sgw,[(1+Gw*Sgw)*(y*s)/sw*x2-1],[(1+Gw*Sgw)*(y*s)/sw*Y-1], B0s,B0t; print; * .store * #endprocedure * #procedure renormv(dummy) * g RWIv= NWIv*Pfs*Pft/aux; * id Pws^-1= -Gs/g^2/Pfs; al Pwt^-1= -Gt/g^2/Pft; * id Pfs= x2*y*s-sw+(1-Gs/Gw)*sw+Gs*fws; id Pft= Y*y*s-sw+(1-Gt/Gw)*sw+Gt*fwt; * id fws= 0; al fwt= 0; * id wvs1= 0; al wvs2= 0; al wsv1= 0; al wsv2= 0; al wss= 0; * id wvv2= wvv1; al wvv3= wvv1; * id aux^-1= 1/(Gs*Gt); * .store * g NRWIv= RWIv; * id Rm_s= g^2/Gs-1; al Rm_t= g^2/Gt-1; id M^2= g^2/Gw*sw; al M^-2= Gw/g^2/sw; * id [x2-M2/(y*s)]^-1*[Y-M2/(y*s)]^-1= ([x2-M2/(y*s)]^-1-[Y-M2/(y*s)]^-1)/ [Y-x2]; * id [x2-M2/(y*s)]^-1= Gw/g^2*y*s/sw/[Gw/g^2*(y*s)/sw*x2-1]; id [Y-M2/(y*s)]^-1= Gw/g^2*y*s/sw/[Gw/g^2*(y*s)/sw*Y-1]; * id g^-2= 1/2*(1/Gt+Sgt+1/Gs+Sgs); * id [Gw/g^2*(y*s)/sw*x2-1]^-1= 1/[(1+Gw*Sgw)*(y*s)/sw*x2-1]; id [Gw/g^2*(y*s)/sw*Y-1]^-1= 1/[(1+Gw*Sgw)*(y*s)/sw*Y-1]; * id wvv0= y*s*(4/9*(x2-Y)+4/3*(Y*B0t-x2*B0s)); al wvv1= 4/9-4/3*B0s; al wvv2= 4/9-4/3*B0s; al wvv3= 4/9-4/3*B0s; al wvv4= 4/3*(B0t-B0s); * ratio,[Y-x2],Y,x2; * id Gw*Gs^-1= 1+Gw*(Sgw-Sgs); id Gw*Gt^-1= 1+Gw*(Sgw-Sgt); * id Sgs= 4/9-4/3*B0s; al Sgt= 4/9-4/3*B0t; * b i_,g,st,y,s,em,um,dm,u,ub,v,vb,Pws,Pwt,Pps,Ppt,D_s,D_t,Sws,Swt,Sps,Spt, [BQS+M2],[QpS+M2],[(qp+BQ)^2+em^2],[(Qm-k)^2+um^2],[(Qm-kb)^2+dm^2], Qu,Qd,[BQS-Sps],[QpS-Spt],ST_l,ST_q,ST,[em^2-y*s],[um^2-t],[dm^2-u], [BQS+M2_s],[QpS+M2_t],[x2-M2/(y*s)],[Y-M2/(y*s)],Rm_s,Rm_t,Gs,Gt,Gw,fws,fwt, fww,Sgs,Sgt,Sgw,[(1+Gw*Sgw)*(y*s)/sw*x2-1],[(1+Gw*Sgw)*(y*s)/sw*Y-1], B0s,B0t; print; * .store * #endprocedure * #procedure renormM(dummy) * g RWI= NWI*Pfs*Pft/aux; * id Pws^-1= -Gs/g^2/Pfs; al Pwt^-1= -Gt/g^2/Pft; * id Pfs= x2*y*s-sw+(1-Gs/Gw)*sw+Gs*(fws-fww); id Pft= Y*y*s-sw+(1-Gt/Gw)*sw+Gt*(fwt-fww); * id Rm_s= g^2/Gs*eta_s*(1-Gw*sigw)-1; al Rm_t= g^2/Gt*eta_t*(1-Gw*sigw)-1; * id M^2= g^2*sw*(1/Gw-sigw); al M^-2= Gw/g^2/sw/[1-Gw*sigw]; * id [x2-M2/(y*s)]^-1= Gw/g^2*y*s/sw/[(1+Gw*Sgw)*x2*y*s/sw+Gw*sigw-1]; id [Y-M2/(y*s)]^-1= Gw/g^2*y*s/sw/[(1+Gw*Sgw)*Y*y*s/sw+Gw*sigw-1]; * id g^-2= 1/2*(1/Gt+Sgt+1/Gs+Sgs); * id Gs*Gt*Gw= aux1; id Gs*Gt= aux2; id Gt= Gs+Gs*Gt*(Sgt-Sgs); id aux1= Gs*Gt*Gw; id aux2= Gs*Gt; * id Gs*Gt^-1= 1+Gs*(Sgs-Sgt); * id Gs*Gt*Gw= aux1; id Gt*Gw= Gs*Gw+Gs*Gt*Gw*(Sgt-Sgs); id aux1= Gs*Gt*Gw; * id aux^-1= 1/(Gs*Gt); id Gw*Gt^-1= 1+Gw*(Sgw-Sgt); id Gw*Gs^-1= 1+Gw*(Sgw-Sgs); id Gt^-1= Gs^-1+Sgs-Sgt; * id Gw*Sgw*[(1+Gw*Sgw)*x2*y*s/sw+Gw*sigw-1]^-1= -1/[(1+Gw*Sgw)*x2*y*s/sw+Gw*sigw-1]+sw/(x2*y*s)*(1- (Gw*sigw-1)/[(1+Gw*Sgw)*x2*y*s/sw+Gw*sigw-1]); * al Gw*Sgw*[(1+Gw*Sgw)*Y*y*s/sw+Gw*sigw-1]^-1= -1/[(1+Gw*Sgw)*Y*y*s/sw+Gw*sigw-1]+sw/(Y*y*s)*(1- (Gw*sigw-1)/[(1+Gw*Sgw)*Y*y*s/sw+Gw*sigw-1]); * id sigw= -fww/sw; * b i_,g,st,y,s,em,um,dm,u,ub,v,vb,Pws,Pwt,Pps,Ppt,D_s,D_t,Sws,Swt,Sps,Spt, [BQS+M2],[QpS+M2],[(qp+BQ)^2+em^2],[(Qm-k)^2+um^2],[(Qm-kb)^2+dm^2], Qu,Qd,[BQS-Sps],[QpS-Spt],ST_l,ST_q,ST,[em^2-y*s],[um^2-t],[dm^2-u], [BQS+M2_s],[QpS+M2_t],[x2-M2/(y*s)],[Y-M2/(y*s)],Rm_s,Rm_t,Gs,Gt,Gw,fws,fwt, fww,Sgs,Sgt,Sgw,[(1+Gw*Sgw)*x2*y*s/sw+Gw*sigw-1], [(1+Gw*Sgw)*Y*y*s/sw+Gw*sigw-1],tm,eta_s,eta_t,sw,[1-Gw*sigw]; print; * .store * #endprocedure * #procedure renormMv(dummy) * g RWIv= NWIv*Pfs*Pft/aux; * id Pws^-1= -Gs/g^2/Pfs; al Pwt^-1= -Gt/g^2/Pft; * id Pfs= x2*y*s-sw+(1-Gs/Gw)*sw+Gs*(fws-fww); id Pft= Y*y*s-sw+(1-Gt/Gw)*sw+Gt*(fwt-fww); * id wvs1= i_/M*swvs1; al wvs2= i_/M*swvs2; id wsv1= i_/M*swsv1; al wsv2= i_/M*swsv2; id wss= 1/M^2*swss; * id Rm_s= g^2/Gs*eta_s*(1-Gw*sigw)-1; al Rm_t= g^2/Gt*eta_t*(1-Gw*sigw)-1; id M^2= g^2*sw*(1/Gw-sigw); al M^-2= Gw/g^2/sw/[1-Gw*sigw]; * id sigw= -fww/sw; * id aux^-1= 1/(Gs*Gt); * .store * g NRWIv= RWIv; * id [x2-M2/(y*s)]^-1*[Y-M2/(y*s)]^-1= ([x2-M2/(y*s)]^-1-[Y-M2/(y*s)]^-1)/ [Y-x2]; * id [x2-M2/(y*s)]^-1= Gw/g^2*y*s/sw/[(1+Gw*Sgw)*x2*y*s/sw+Gw*sigw-1]; id [Y-M2/(y*s)]^-1= Gw/g^2*y*s/sw/[(1+Gw*Sgw)*Y*y*s/sw+Gw*sigw-1]; * id g^-2= 1/2*(1/Gt+Sgt+1/Gs+Sgs); * ratio,[Y-x2],Y,x2; * id Gw*Gs^-1= 1+Gw*(Sgw-Sgs); id Gw*Gt^-1= 1+Gw*(Sgw-Sgt); * id Gw*Sgw*[(1+Gw*Sgw)*x2*y*s/sw+Gw*sigw-1]^-1= -1/[(1+Gw*Sgw)*x2*y*s/sw+Gw*sigw-1]+sw/(x2*y*s)*(1- (Gw*sigw-1)/[(1+Gw*Sgw)*x2*y*s/sw+Gw*sigw-1]); al Gw*Sgw*[(1+Gw*Sgw)*Y*y*s/sw+Gw*sigw-1]^-1= -1/[(1+Gw*Sgw)*Y*y*s/sw+Gw*sigw-1]+sw/(Y*y*s)*(1- (Gw*sigw-1)/[(1+Gw*Sgw)*Y*y*s/sw+Gw*sigw-1]); ratio,[Y-x2],Y,x2; * id sigw= -fww/sw; * b i_,g,st,y,s,em,um,dm,u,ub,v,vb,Pws,Pwt,Pps,Ppt,D_s,D_t,Sws,Swt,Sps,Spt, [BQS+M2],[QpS+M2],[(qp+BQ)^2+em^2],[(Qm-k)^2+um^2],[(Qm-kb)^2+dm^2], Qu,Qd,[BQS-Sps],[QpS-Spt],ST_l,ST_q,ST,[em^2-y*s],[um^2-t],[dm^2-u], [BQS+M2_s],[QpS+M2_t],[x2-M2/(y*s)],[Y-M2/(y*s)],Rm_s,Rm_t,Gs,Gt,Gw,fws,fwt, fww,Sgs,Sgt,Sgw,[(1+Gw*Sgw)*x2*y*s/sw+Gw*sigw-1], [(1+Gw*Sgw)*Y*y*s/sw+Gw*sigw-1],tm,sw,[1-Gw*sigw]; print; * .store * #endprocedure * #procedure totalm(dummy) * g TOTa= (RWI+NRWIv)*sela; g TOTb= (RWI+NRWIv)*selb; g TOTc= (RWI+NRWIv)*selc; g TOTd= (RWI+NRWIv)*seld; * id Gw*Sgw*[(1+Gw*Sgw)*(y*s)/sw*Y-1]^-1= -[(1+Gw*Sgw)*(y*s)/sw*Y-1]^-1+sw/(y*s*Y)*(1+[(1+Gw*Sgw)*(y*s)/sw*Y-1]^-1); al Gw*Sgw*[(1+Gw*Sgw)*(y*s)/sw*x2-1]^-1= -[(1+Gw*Sgw)*(y*s)/sw*x2-1]^-1+sw/(y*s*x2)*(1+[(1+Gw*Sgw)*(y*s)/sw*x2-1]^-1); * id sela*ST(Qm,p,1,[pm])= 1; id selb*ST(1,p,1,[pm])= 1; id selc*ST_l(1,p)*ST_q(Qm,p)= 1; id seld*ST_l(mu?,p)*ST_q(mu?,p)= 1; id sela= 0; id selb= 0; id selc= 0; id seld= 0; * b i_,em,ST_l,ST_q,ST,Gw,y,s,sw,x2,Y,X, [(1+Gw*Sgw)*(y*s)/sw*x2-1],[(1+Gw*Sgw)*(y*s)/sw*Y-1]; print; * #endprocedure * #procedure totalM(dummy) * g TOTa= (RWI+NRWIv)*sela*Deta_s; g TOTc= (RWI+NRWIv)*selc*Deta_t; g TOTd= (RWI+NRWIv)*seld; * id sela*ST(Qm,p,1,[pm])= 1; id selb*ST(1,p,1,[pm])= 1; id selc*ST_l(1,p)*ST_q(Qm,p)= 1; id seld*ST_l(mu?,p)*ST_q(mu?,p)= 1; id sela= 0; id selb= 0; id selc= 0; id seld= 0; * id eta_s= Neta_s/Deta_s; id Neta_s= 1-Gs*(-fws+tm^2*sigps)/(x2*y*s); al Deta_s= 1-Gw*(-fww+tm^2*sigps)/sw; * id eta_t= Neta_t/Deta_t; id Neta_t= 1-Gt*(-fwt+tm^2*sigpt)/(Y*y*s); al Deta_t= 1-Gw*(-fww+tm^2*sigpt)/sw; * id Gw*Gs^-1= 1+Gw*(Sgw-Sgs); * id Gt*Gs^-1= 1+Gt*(Sgt-Sgs); * id Gt= Gw+Gw*Gt*(Sgt-Sgw); id Gs= Gw+Gw*Gs*(Sgs-Sgw); * id sigw= -fww/sw; * ratio,[Y-x2],Y,x2; * .store * g NTOTa= TOTa*[1-Gw*sigw]; g NTOTc= TOTc*[1-Gw*sigw]; g NTOTd= TOTd; * id [1-Gw*sigw]= 1-Gw*sigw; * id sigw= -fww/sw; * repeat; id Gs*Gw*Sgs= Gs*Gw*Sgw+Gs-Gw; id Gt*Gw*Sgt= Gt*Gw*Sgw+Gt-Gw; endrepeat; * id Gw*Sgw*[(1+Gw*Sgw)*x2*y*s/sw+Gw*sigw-1]^-1= -1/[(1+Gw*Sgw)*x2*y*s/sw+Gw*sigw-1]+sw/(x2*y*s)*(1- (Gw*sigw-1)/[(1+Gw*Sgw)*x2*y*s/sw+Gw*sigw-1]); * al Gw*Sgw*[(1+Gw*Sgw)*Y*y*s/sw+Gw*sigw-1]^-1= -1/[(1+Gw*Sgw)*Y*y*s/sw+Gw*sigw-1]+sw/(Y*y*s)*(1- (Gw*sigw-1)/[(1+Gw*Sgw)*Y*y*s/sw+Gw*sigw-1]); * id sigw= -fww/sw; ratio,[Y-x2],Y,x2; * .store * g NNTOTa= NTOTa*[(1+Gw*Sgw)*x2*y*s/sw+Gw*sigw-1]; g NNTOTc= NTOTc*[(1+Gw*Sgw)*Y*y*s/sw+Gw*sigw-1]; g NNTOTd= NTOTd; * id [(1+Gw*Sgw)*x2*y*s/sw+Gw*sigw-1]= (1+Gw*Sgw)*x2*y*s/sw+Gw*sigw-1; id [(1+Gw*Sgw)*Y*y*s/sw+Gw*sigw-1]= (1+Gw*Sgw)*Y*y*s/sw+Gw*sigw-1; * repeat; id Gs*Gw*Sgs= Gs*Gw*Sgw+Gs-Gw; id Gt*Gw*Sgt= Gt*Gw*Sgw+Gt-Gw; endrepeat; * id sigw= -fww/sw; ratio,[Y-x2],Y,x2; * *------------------------------------------ * id fws= -x2*y*s*sigs; al fwt= -Y*y*s*sigt; al fww= -sw*sigw; * id wvv2= wvv1; id wvv3= wvv1; * id swsv1= -swvs1; id swsv2= -swvs2; * *------------------------------------------ * .store * g NNNTOTd= NNTOTd; * id wvv0= 3*y*s*(4/9*(x2-Y)+4/3*(Y*B0t_mm-x2*B0s_mm)) + s*y * ( 1/3*x2 - 1/3*Y ) + B0t_mM*s*y^-1*tmr^4 * ( - 1/2*Y^-1 ) + B0t_mM*s*y * ( Y ) + B0t_mM*s*tmr^2 * ( - 1/2 ) + B0s_mM*s*y^-1*tmr^4 * ( 1/2*x2^-1 ) + B0s_mM*s*y * ( - x2 ) + B0s_mM*s*tmr^2 * ( 1/2 ) + A0*y^-1*tmr^2 * ( 1/2*x2^-1 - 1/2*Y^-1 ) - s*y * ( 1/3*x2 - 1/3*Y ) - B0t_mm*s*y * ( Y ) - B0s_mm*s*y * ( - x2 ); * id Sgs= 12*(2*B21s_mm-B0s_mm)+2*(2*(B21s_MM-B21s_mm)-(B0s_MM-B0s_mm)); al Sgt= 12*(2*B21t_mm-B0t_mm)+2*(2*(B21t_MM-B21t_mm)-(B0t_MM-B0t_mm)); id sigs= -fws/(x2*y*s); id sigt= -fwt/(Y*y*s); id fws= -1/3*tm^2 -4/3*B0s_MM*tm^2 -2/3*x2*y*s*B0s_MM -1/2*B0s_mM*tm^2 -1/2/(x2*y*s)*B0s_mM*tm^4 +x2*y*s*B0s_mM -1/3*x2*y*s*B0s_mm -1/2/(x2*y*s)*A0*tm^2 -1/3*A0; al fwt= -1/3*tm^2 -4/3*B0t_MM*tm^2 -2/3*Y*y*s*B0t_MM -1/2*B0t_mM*tm^2 -1/2/(Y*y*s)*B0t_mM*tm^4 +Y*y*s*B0t_mM -1/3*Y*y*s*B0t_mm -1/2/(Y*y*s)*A0*tm^2 -1/3*A0; id B21s_MM= -1/3/(x2*y*s)*tm^2+1/3*(1-tm^2/(x2*y*s))*B0s_MM- 1/3/(x2*y*s)*A0+1/18; al B21t_MM= -1/3/(Y*y*s)*tm^2+1/3*(1-tm^2/(Y*y*s))*B0t_MM-1/3/(Y*y*s)*A0+1/18; * id B21s_mm= 1/18+1/3*B0s_mm; id B21t_mm= 1/18+1/3*B0t_mm; id A0= tm^2*RA0; id tmr^2= tm^2/s; * b i_,g,x2,Y,y,s,em,tm,Gs,Gt,Gw,sw; print; * .store * g NNNTOTa= NNTOTa; g NNNTOTc= NNTOTc; * id sigs= -fws/(x2*y*s); al sigt= -fwt/(Y*y*s); al sigw= -fww/sw; * id fws= -1/3*tm^2 -4/3*B0s_MM*tm^2 -2/3*x2*y*s*B0s_MM -1/2*B0s_mM*tm^2 -1/2/(x2*y*s)*B0s_mM*tm^4 +x2*y*s*B0s_mM -1/3*x2*y*s*B0s_mm -1/2/(x2*y*s)*A0*tm^2 -1/3*A0; al fwt= -1/3*tm^2 -4/3*B0t_MM*tm^2 -2/3*Y*y*s*B0t_MM -1/2*B0t_mM*tm^2 -1/2/(Y*y*s)*B0t_mM*tm^4 +Y*y*s*B0t_mM -1/3*Y*y*s*B0t_mm -1/2/(Y*y*s)*A0*tm^2 -1/3*A0; id fww= -1/3*tm^2 -4/3*B0w_MM*tm^2 -2/3*sw*B0w_MM -1/2*B0w_mM*tm^2 -1/2/sw*B0w_mM*tm^4 +sw*B0w_mM -1/3*sw*B0w_mm -1/2/sw*A0*tm^2 -1/3*A0; * id sigps= 3/2/(x2*y*s)*(A0+(tm^2-x2*y*s)*B0s_mM); al sigpt= 3/2/(Y*y*s)*(A0+(tm^2-Y*y*s)*B0t_mM); * id Sgs= 12*(2*B21s_mm-B0s_mm)+2*(2*(B21s_MM-B21s_mm)-(B0s_MM-B0s_mm)); al Sgt= 12*(2*B21t_mm-B0t_mm)+2*(2*(B21t_MM-B21t_mm)-(B0t_MM-B0t_mm)); al Sgw= 12*(2*B21w_mm-B0w_mm)+2*(2*(B21w_MM-B21w_mm)-(B0w_MM-B0w_mm)); * id wvv0= 3*y*s*(4/9*(x2-Y)+4/3*(Y*B0t_mm-x2*B0s_mm)) + s*y * ( 1/3*x2 - 1/3*Y ) + B0t_mM*s*y^-1*tmr^4 * ( - 1/2*Y^-1 ) + B0t_mM*s*y * ( Y ) + B0t_mM*s*tmr^2 * ( - 1/2 ) + B0s_mM*s*y^-1*tmr^4 * ( 1/2*x2^-1 ) + B0s_mM*s*y * ( - x2 ) + B0s_mM*s*tmr^2 * ( 1/2 ) + A0*y^-1*tmr^2 * ( 1/2*x2^-1 - 1/2*Y^-1 ) - s*y * ( 1/3*x2 - 1/3*Y ) - B0t_mm*s*y * ( Y ) - B0s_mm*s*y * ( - x2 ); al wvv1= 3*(4/9-4/3*B0s_mm) + y^-1*tmr^2 * ( - x2^-1 ) + 1/3 + B0s_mM*y^-2*tmr^4 * ( 2*x2^-2 ) + B0s_mM*y^-1*tmr^2 * ( - x2^-1 ) - B0s_mM + A0*s^-1*y^-2*tmr^2 * ( 2*x2^-2 ) + A0*s^-1*y^-1 * ( - x2^-1 ) - 1/3 + B0s_mm; id wvv4= 3*4/3*(B0t_mm-B0s_mm) + y^-1*tmr^2 * ( - x2^-1 + Y^-1 ) + B0t_mM*y^-2*tmr^4 * ( - 2*Y^-2 ) + B0t_mM*y^-1*tmr^2 * ( Y^-1 ) + B0t_mM + B0s_mM*y^-2*tmr^4 * ( 2*x2^-2 ) + B0s_mM*y^-1*tmr^2 * ( - x2^-1 ) - B0s_mM + A0*s^-1*y^-2*tmr^2 * ( 2*x2^-2 - 2*Y^-2 ) + A0*s^-1*y^-1 * ( - x2^-1 + Y^-1 ) - B0t_mm + B0s_mm; * id swvs1= B0t_mM*s*y^-1*tmr^4 * ( 3/2*Y^-1 ) + B0t_mM*s*tmr^2 * ( - 3/2 ) + B0s_mM*s*y^-1*tmr^4 * ( - 3/2*x2^-1 ) + B0s_mM*s*tmr^2 * ( 3/2 ) + A0*y^-1*tmr^2 * ( - 3/2*x2^-1 + 3/2*Y^-1 ); al swvs2= B0s_mM*s*y^-1*tmr^4 * ( - 3/2*x2^-1 ) + B0s_mM*s*tmr^2 * ( 3/2 ) + A0*y^-1*tmr^2 * ( - 3/2*x2^-1 ); * *id swsv1= B0t_mM*s*y^-1*tmr^4 * ( - 3/2*Y^-1 ) * + B0t_mM*s*tmr^2 * ( 3/2 ) * + B0s_mM*s*y^-1*tmr^4 * ( 3/2*x2^-1 ) * + B0s_mM*s*tmr^2 * ( - 3/2 ) * + A0*y^-1*tmr^2 * ( 3/2*x2^-1 - 3/2*Y^-1 ); *al swsv2= B0s_mM*s*y^-1*tmr^4 * ( 3/2*x2^-1 ) * + B0s_mM*s*tmr^2 * ( - 3/2 ) * + A0*y^-1*tmr^2 * ( 3/2*x2^-1 ); * id swss= B0t_mM*s^2*y*tmr^2 * ( 3/2*Y ) + B0t_mM*s^2*tmr^4 * ( - 3/2 ) + B0s_mM*s^2*y*tmr^2 * ( - 3/2*x2 ) + B0s_mM*s^2*tmr^4 * ( 3/2 ); * id tmr^2= tm^2/s; * id B21s_MM= -1/3/(x2*y*s)*tm^2+1/3*(1-tm^2/(x2*y*s))*B0s_MM- 1/3/(x2*y*s)*A0+1/18; al B21t_MM= -1/3/(Y*y*s)*tm^2+1/3*(1-tm^2/(Y*y*s))*B0t_MM-1/3/(Y*y*s)*A0+1/18; id B21w_MM= -1/3/sw*tm^2+1/3*(1-tm^2/sw)*B0w_MM-1/3/sw*A0+1/18; * id B21s_mm= 1/18+1/3*B0s_mm; id B21t_mm= 1/18+1/3*B0t_mm; id B21w_mm= 1/18+1/3*B0w_mm; id A0= tm^2*RA0; * b i_,g,x2,Y,y,s,em,tm,Gs,Gt,Gw,sw; print; * .store * #endprocedure * #procedure totalMB(dummy) * g ATOTb= (RWI+NRWIv)*selb*Deta_s*Deta_t; * id sela*ST(Qm,p,1,[pm])= 1; id selb*ST(1,p,1,[pm])= 1; id selc*ST_l(1,p)*ST_q(Qm,p)= 1; id seld*ST_l(mu?,p)*ST_q(mu?,p)= 1; id sela= 0; id selb= 0; id selc= 0; id seld= 0; * id eta_s= Neta_s/Deta_s; id Neta_s= 1-Gs*(-fws+tm^2*sigps)/(x2*y*s); al Deta_s= 1-Gw*(-fww+tm^2*sigps)/sw; * id eta_t= Neta_t/Deta_t; id Neta_t= 1-Gt*(-fwt+tm^2*sigpt)/(Y*y*s); al Deta_t= 1-Gw*(-fww+tm^2*sigpt)/sw; * id Gw*Gs^-1= 1+Gw*(Sgw-Sgs); * id Gt*Gs^-1= 1+Gt*(Sgt-Sgs); * id Gt= Gw+Gw*Gt*(Sgt-Sgw); id Gs= Gw+Gw*Gs*(Sgs-Sgw); * id sigw= -fww/sw; * ratio,[Y-x2],Y,x2; * .store * g BTOTb= ATOTb*[1-Gw*sigw]^2; * id [1-Gw*sigw]= 1-Gw*sigw; * id sigw= -fww/sw; * id Gw*Sgw*[(1+Gw*Sgw)*x2*y*s/sw+Gw*sigw-1]^-1= -1/[(1+Gw*Sgw)*x2*y*s/sw+Gw*sigw-1]+sw/(x2*y*s)*(1- (Gw*sigw-1)/[(1+Gw*Sgw)*x2*y*s/sw+Gw*sigw-1]); * al Gw*Sgw*[(1+Gw*Sgw)*Y*y*s/sw+Gw*sigw-1]^-1= -1/[(1+Gw*Sgw)*Y*y*s/sw+Gw*sigw-1]+sw/(Y*y*s)*(1- (Gw*sigw-1)/[(1+Gw*Sgw)*Y*y*s/sw+Gw*sigw-1]); * repeat; id Gs*Gw*Sgs= Gs*Gw*Sgw+Gs-Gw; id Gt*Gw*Sgt= Gt*Gw*Sgw+Gt-Gw; endrepeat; * id sigw= -fww/sw; ratio,[Y-x2],Y,x2; * .store * g CTOTb= BTOTb; * id Gw*Gt^-1= 1+Gw*(Sgw-Sgt); id Gs*Gt^-1= 1+Gs*(Sgs-Sgt); * id Gs*Gt*Sgt= Gs*Gt*Sgs-Gs+Gt; id Gs*Gw*Sgs= Gs*Gw*Sgw+Gs-Gw; id Gt*Gw*Sgt= Gt*Gw*Sgw+Gt-Gw; * .store * g DTOTb= CTOTb*[1+Gw*Sgw]; * id g^2= Gw/[1+Gw*Sgw]; id [1+Gw*Sgw]= 1+Gw*Sgw; id Gw*Gt^-1= 1+Gw*(Sgw-Sgt); * id Gw*Sgw*[(1+Gw*Sgw)*x2*y*s/sw+Gw*sigw-1]^-1= -1/[(1+Gw*Sgw)*x2*y*s/sw+Gw*sigw-1]+sw/(x2*y*s)*(1- (Gw*sigw-1)/[(1+Gw*Sgw)*x2*y*s/sw+Gw*sigw-1]); * al Gw*Sgw*[(1+Gw*Sgw)*Y*y*s/sw+Gw*sigw-1]^-1= -1/[(1+Gw*Sgw)*Y*y*s/sw+Gw*sigw-1]+sw/(Y*y*s)*(1- (Gw*sigw-1)/[(1+Gw*Sgw)*Y*y*s/sw+Gw*sigw-1]); * .store * g ETOTb= DTOTb*[(1+Gw*Sgw)*x2*y*s/sw+Gw*sigw-1]* [(1+Gw*Sgw)*Y*y*s/sw+Gw*sigw-1]; * id [(1+Gw*Sgw)*x2*y*s/sw+Gw*sigw-1]= (1+Gw*Sgw)*x2*y*s/sw+Gw*sigw-1; id [(1+Gw*Sgw)*Y*y*s/sw+Gw*sigw-1]= (1+Gw*Sgw)*Y*y*s/sw+Gw*sigw-1; * repeat; id Gs*Gw*Sgs= Gs*Gw*Sgw+Gs-Gw; id Gt*Gw*Sgt= Gt*Gw*Sgw+Gt-Gw; endrepeat; * ratio,[Y-x2],Y,x2; * *------------------------------------------ * id fws= -x2*y*s*sigs; al fwt= -Y*y*s*sigt; al fww= -sw*sigw; * id wvv2= wvv1; id wvv3= wvv1; * id swsv1= -swvs1; id swsv2= -swvs2; * *------------------------------------------ * .store * save outputB.sav ETOTb; * #endprocedure * #procedure totalMBcnt(dummy) * load outputB.sav; * g FTOTb= ETOTb; * id sigw= (1-[1-Gw*sigw])/Gw; id Sgw= ([1+Gw*Sgw]-1)/Gw; * .store * g GTOTb= FTOTb; * id sigs= -fws/(x2*y*s); al sigt= -fwt/(Y*y*s); al sigw= -fww/sw; * id fws= -1/3*tm^2 -4/3*B0s_MM*tm^2 -2/3*x2*y*s*B0s_MM -1/2*B0s_mM*tm^2 -1/2/(x2*y*s)*B0s_mM*tm^4 +x2*y*s*B0s_mM -1/3*x2*y*s*B0s_mm -1/2/(x2*y*s)*A0*tm^2 -1/3*A0; al fwt= -1/3*tm^2 -4/3*B0t_MM*tm^2 -2/3*Y*y*s*B0t_MM -1/2*B0t_mM*tm^2 -1/2/(Y*y*s)*B0t_mM*tm^4 +Y*y*s*B0t_mM -1/3*Y*y*s*B0t_mm -1/2/(Y*y*s)*A0*tm^2 -1/3*A0; id fww= -1/3*tm^2 -4/3*B0w_MM*tm^2 -2/3*sw*B0w_MM -1/2*B0w_mM*tm^2 -1/2/sw*B0w_mM*tm^4 +sw*B0w_mM -1/3*sw*B0w_mm -1/2/sw*A0*tm^2 -1/3*A0; * id sigps= 3/2/(x2*y*s)*(A0+(tm^2-x2*y*s)*B0s_mM); al sigpt= 3/2/(Y*y*s)*(A0+(tm^2-Y*y*s)*B0t_mM); * id Sgs= 12*(2*B21s_mm-B0s_mm)+2*(2*(B21s_MM-B21s_mm)-(B0s_MM-B0s_mm)); al Sgt= 12*(2*B21t_mm-B0t_mm)+2*(2*(B21t_MM-B21t_mm)-(B0t_MM-B0t_mm)); al Sgw= 12*(2*B21w_mm-B0w_mm)+2*(2*(B21w_MM-B21w_mm)-(B0w_MM-B0w_mm)); * id wvv0= 3*y*s*(4/9*(x2-Y)+4/3*(Y*B0t_mm-x2*B0s_mm)) + s*y * ( 1/3*x2 - 1/3*Y ) + B0t_mM*s*y^-1*tmr^4 * ( - 1/2*Y^-1 ) + B0t_mM*s*y * ( Y ) + B0t_mM*s*tmr^2 * ( - 1/2 ) + B0s_mM*s*y^-1*tmr^4 * ( 1/2*x2^-1 ) + B0s_mM*s*y * ( - x2 ) + B0s_mM*s*tmr^2 * ( 1/2 ) + A0*y^-1*tmr^2 * ( 1/2*x2^-1 - 1/2*Y^-1 ) - s*y * ( 1/3*x2 - 1/3*Y ) - B0t_mm*s*y * ( Y ) - B0s_mm*s*y * ( - x2 ); al wvv1= 3*(4/9-4/3*B0s_mm) + y^-1*tmr^2 * ( - x2^-1 ) + 1/3 + B0s_mM*y^-2*tmr^4 * ( 2*x2^-2 ) + B0s_mM*y^-1*tmr^2 * ( - x2^-1 ) - B0s_mM + A0*s^-1*y^-2*tmr^2 * ( 2*x2^-2 ) + A0*s^-1*y^-1 * ( - x2^-1 ) - 1/3 + B0s_mm; id wvv4= 3*4/3*(B0t_mm-B0s_mm) + y^-1*tmr^2 * ( - x2^-1 + Y^-1 ) + B0t_mM*y^-2*tmr^4 * ( - 2*Y^-2 ) + B0t_mM*y^-1*tmr^2 * ( Y^-1 ) + B0t_mM + B0s_mM*y^-2*tmr^4 * ( 2*x2^-2 ) + B0s_mM*y^-1*tmr^2 * ( - x2^-1 ) - B0s_mM + A0*s^-1*y^-2*tmr^2 * ( 2*x2^-2 - 2*Y^-2 ) + A0*s^-1*y^-1 * ( - x2^-1 + Y^-1 ) - B0t_mm + B0s_mm; * id swvs1= B0t_mM*s*y^-1*tmr^4 * ( 3/2*Y^-1 ) + B0t_mM*s*tmr^2 * ( - 3/2 ) + B0s_mM*s*y^-1*tmr^4 * ( - 3/2*x2^-1 ) + B0s_mM*s*tmr^2 * ( 3/2 ) + A0*y^-1*tmr^2 * ( - 3/2*x2^-1 + 3/2*Y^-1 ); al swvs2= B0s_mM*s*y^-1*tmr^4 * ( - 3/2*x2^-1 ) + B0s_mM*s*tmr^2 * ( 3/2 ) + A0*y^-1*tmr^2 * ( - 3/2*x2^-1 ); * *id swsv1= B0t_mM*s*y^-1*tmr^4 * ( - 3/2*Y^-1 ) * + B0t_mM*s*tmr^2 * ( 3/2 ) * + B0s_mM*s*y^-1*tmr^4 * ( 3/2*x2^-1 ) * + B0s_mM*s*tmr^2 * ( - 3/2 ) * + A0*y^-1*tmr^2 * ( 3/2*x2^-1 - 3/2*Y^-1 ); *al swsv2= B0s_mM*s*y^-1*tmr^4 * ( 3/2*x2^-1 ) * + B0s_mM*s*tmr^2 * ( - 3/2 ) * + A0*y^-1*tmr^2 * ( 3/2*x2^-1 ); * id swss= B0t_mM*s^2*y*tmr^2 * ( 3/2*Y ) + B0t_mM*s^2*tmr^4 * ( - 3/2 ) + B0s_mM*s^2*y*tmr^2 * ( - 3/2*x2 ) + B0s_mM*s^2*tmr^4 * ( 3/2 ); * id tmr^2= tm^2/s; * id B21s_MM= -1/3/(x2*y*s)*tm^2+1/3*(1-tm^2/(x2*y*s))*B0s_MM- 1/3/(x2*y*s)*A0+1/18; al B21t_MM= -1/3/(Y*y*s)*tm^2+1/3*(1-tm^2/(Y*y*s))*B0t_MM-1/3/(Y*y*s)*A0+1/18; id B21w_MM= -1/3/sw*tm^2+1/3*(1-tm^2/sw)*B0w_MM-1/3/sw*A0+1/18; * id B21s_mm= 1/18+1/3*B0s_mm; id B21t_mm= 1/18+1/3*B0t_mm; id B21w_mm= 1/18+1/3*B0w_mm; id A0= tm^2*RA0; * b i_,g,x2,Y,X,y,s,em,tm,Gs,Gt,Gw,sw,[1-Gw*sigw],[1+Gw*Sgw]; print; * #endprocedure * #procedure connect(dummy) * g eq1= wvv0+z1*Sgs+z2*Sgt+z3*sigs+z4*sigt; g eq2= wvv1+z5*Sgs+z6*sigs+z7*sigps; g eq3= wvv4+z8*Sgs+z9*Sgt+z10*sigs+z11*sigt+z12*sigps+z13*sigpt; g eq4= swvs1+z14*sigs+z15*sigt+z16*sigps+z17*sigpt; g eq5= swvs2+z18*sigs+z19*sigps; g eq6= swss+z20*sigs+z21*sigt+z22*sigps+z23*sigpt; * id wvv0= 3*y*s*(4/9*(x2-Y)+4/3*(Y*B0t_mm-x2*B0s_mm)) + s*y * ( 1/3*x2 - 1/3*Y ) + B0t_mM*s*y^-1*tmr^4 * ( - 1/2*Y^-1 ) + B0t_mM*s*y * ( Y ) + B0t_mM*s*tmr^2 * ( - 1/2 ) + B0s_mM*s*y^-1*tmr^4 * ( 1/2*x2^-1 ) + B0s_mM*s*y * ( - x2 ) + B0s_mM*s*tmr^2 * ( 1/2 ) + A0*y^-1*tmr^2 * ( 1/2*x2^-1 - 1/2*Y^-1 ) - s*y * ( 1/3*x2 - 1/3*Y ) - B0t_mm*s*y * ( Y ) - B0s_mm*s*y * ( - x2 ); al wvv1= 3*(4/9-4/3*B0s_mm) + y^-1*tmr^2 * ( - x2^-1 ) + 1/3 + B0s_mM*y^-2*tmr^4 * ( 2*x2^-2 ) + B0s_mM*y^-1*tmr^2 * ( - x2^-1 ) - B0s_mM + A0*s^-1*y^-2*tmr^2 * ( 2*x2^-2 ) + A0*s^-1*y^-1 * ( - x2^-1 ) - 1/3 + B0s_mm; id wvv4= 3*4/3*(B0t_mm-B0s_mm) + y^-1*tmr^2 * ( - x2^-1 + Y^-1 ) + B0t_mM*y^-2*tmr^4 * ( - 2*Y^-2 ) + B0t_mM*y^-1*tmr^2 * ( Y^-1 ) + B0t_mM + B0s_mM*y^-2*tmr^4 * ( 2*x2^-2 ) + B0s_mM*y^-1*tmr^2 * ( - x2^-1 ) - B0s_mM + A0*s^-1*y^-2*tmr^2 * ( 2*x2^-2 - 2*Y^-2 ) + A0*s^-1*y^-1 * ( - x2^-1 + Y^-1 ) - B0t_mm + B0s_mm; id swvs1= B0t_mM*s*y^-1*tmr^4 * ( 3/2*Y^-1 ) + B0t_mM*s*tmr^2 * ( - 3/2 ) + B0s_mM*s*y^-1*tmr^4 * ( - 3/2*x2^-1 ) + B0s_mM*s*tmr^2 * ( 3/2 ) + A0*y^-1*tmr^2 * ( - 3/2*x2^-1 + 3/2*Y^-1 ); al swvs2= B0s_mM*s*y^-1*tmr^4 * ( - 3/2*x2^-1 ) + B0s_mM*s*tmr^2 * ( 3/2 ) + A0*y^-1*tmr^2 * ( - 3/2*x2^-1 ); id swss= B0t_mM*s^2*y*tmr^2 * ( 3/2*Y ) + B0t_mM*s^2*tmr^4 * ( - 3/2 ) + B0s_mM*s^2*y*tmr^2 * ( - 3/2*x2 ) + B0s_mM*s^2*tmr^4 * ( 3/2 ); * id Sgs= 12*(2*B21s_mm-B0s_mm)+2*(2*(B21s_MM-B21s_mm)-(B0s_MM-B0s_mm)); al Sgt= 12*(2*B21t_mm-B0t_mm)+2*(2*(B21t_MM-B21t_mm)-(B0t_MM-B0t_mm)); * id sigs= -fws/(x2*y*s); al sigt= -fwt/(Y*y*s); * id fws= -1/3*tm^2 -4/3*B0s_MM*tm^2 -2/3*x2*y*s*B0s_MM -1/2*B0s_mM*tm^2 -1/2/(x2*y*s)*B0s_mM*tm^4 +x2*y*s*B0s_mM -1/3*x2*y*s*B0s_mm -1/2/(x2*y*s)*A0*tm^2 -1/3*A0; al fwt= -1/3*tm^2 -4/3*B0t_MM*tm^2 -2/3*Y*y*s*B0t_MM -1/2*B0t_mM*tm^2 -1/2/(Y*y*s)*B0t_mM*tm^4 +Y*y*s*B0t_mM -1/3*Y*y*s*B0t_mm -1/2/(Y*y*s)*A0*tm^2 -1/3*A0; * id sigps= 3/2/(x2*y*s)*(A0+(tm^2-x2*y*s)*B0s_mM); al sigpt= 3/2/(Y*y*s)*(A0+(tm^2-Y*y*s)*B0t_mM); * id B21s_MM= -1/3/(x2*y*s)*tm^2+1/3*(1-tm^2/(x2*y*s))*B0s_MM- 1/3/(x2*y*s)*A0+1/18; al B21t_MM= -1/3/(Y*y*s)*tm^2+1/3*(1-tm^2/(Y*y*s))*B0t_MM-1/3/(Y*y*s)*A0+1/18; * id B21s_mm= 1/18+1/3*B0s_mm; id B21t_mm= 1/18+1/3*B0t_mm; id A0= tm^2*RA0; id tmr^2= tm^2/s; * id z4= z2; al z3= z1; * id z1= -x2*y*s; al z2= Y*y*s; * id z6= z5; id z5= -1; id z7= -tm^2/(x2*y*s); * id z10= z8; al z11= z9; id z9= 1; al z8= -1; id z12= -tm^2/(x2*y*s); al z13= tm^2/(Y*y*s); * id z14= 0; al z15= 0; id z16= tm^2; al z17= -tm^2; * id z18= 0; al z19= tm^2; * id z20= 0; al z21= 0; id z22= -tm^2*x2*y*s; al z23= tm^2*Y*y*s; * b tm,x2,Y,y,s,RA0,B0s_MM,B0s_mm,B0s_mM,B0t_MM,B0t_mm,B0t_mM; print; * #endprocedure * */////////////////////////////////////////////////////////////////////// * * Here goes the MAIN * .end