* * General stuff * * A T1 diagram is made of a Left 3-Vertex * Up 3-Vertex * Down 3-Vertex * 1/2/3 Propagator * * A T2 diagram is made of a Left 3-Vertex * Right 4-Vertex * * A T3 diagram is made of a Left 4-Vertex * Up 3-Vertex * * A T4 diagram is made of a Left 4-Vertex * Down 3-Vertex * * * A Form diagram is defined by a momentum-arrow (in short, m-arrow) * * Everything is oriented counter-clockwise (this defines the loop m-arrow) * For external lines the m-arrow is inwards * Upper W^+ has outgoing arrow * Lower W^- has ingoing arrow * All fields but gamma, Z, H, phi0 have a flow-arrow (in short f-arrow) * **************************** * i must be a neutral line * **************************** * * 0 = vertex * p+p1+p2= 0 * den1 q * den2 q+p1 * den3 p+p1+p2 = q-p * * TOPOLOGIES * * * / l,al,p2 * 0 * mup,ku/|nup,qu * q-p / | * sl,kl / | * i,mu,p ----0 | q+p1 T1 * nl,ql\ | * q \ | * mdo,pd\|sdo,kd * 0 * \ j,be,p1 * * * mup,ku * _____ l,al,p2 * / \ / * sl,kl/ \ / * i,mu,p ----0 0 T2 * nl,ql\ / \ * \_____/ \ * j,be,p1 * mdo,pd * * * _____ * / \ mup,ku * sl,kl/ \ * i,mu,p ----0 0---- l,al,p2 T3 * |\ /nup,qu * | \_____/ * | sdo,kd * |______________ j,be,p1 * * * ______________ l,al,p2 * | nup,qu * | _____ * | / \ sdo,kd * |/ \ * i,mu,p ----0 0---- j,be,p1 T4 * \ /mdo,pd * \_____/ * nl,ql * * * List of fields * * 1 = gamma * 2 = z * 3 = W f-arrow along m-arrow * 4 = W f-arrow against m-arrow * 5 = H * 6 = phi0 * 7 = phi f-arrow along m-arrow * 8 = phi f-arrow against m-arrow * 9 = xp f-arrow along m-arrow * 10 = xp f-arrow against m-arrow * 11 = xm f-arrow along m-arrow * 12 = xm f-arrow against m-arrow * 13 = yz f-arrow along m-arrow * 14 = yz f-arrow against m-arrow * 15 = ya f-arrow along m-arrow * 16 = ya f-arrow against m-arrow * 17 = up f-arrow along m-arrow * 18 = up f-arrow against m-arrow * 19 = do f-arrow along m-arrow * 20 = do f-arrow against m-arrow * s n; * i ii; i m1=n,m2=n,m3=n,m4=n,m5=n,m6=n; i de1=n,de2=n,de3=n,de4=n,de5=n,de6=n,de7=n,de8=n,de9=n; i mu=n,nu=n,al=n,be=n,nl=n,sl=n,mup=n,nup=n,sdo=n,mdo=n; i la1=n,la2=n,la3=n,la4=n,la5=n,la6=n,la7=n,la8=n,la9=n,la10=n,la11=n,la12=n; * cf dmf; cf frum; cf w,ws,wss,wb,wh,gram; cf b21_0,b1_0; cf b22,b21,b1,b0,a0,b0p,ra0,db0; cf den1,den2,den3,den2p,den3p; cf qf,pf,p1f,p2f,kl,ku,qu,kd,pd,ql; * f sym211,sym122,pd1,pd2; f vertl,vertu,vertd,verttq,vertqdt,vertqut; f diagvs,diagg,diagf,diagtq,diagqdt,diagqut; f f1,f2,R1,R2,R3,R4,R5,R6,R8.R9,R11,R12,R13,R15; f c0,c11,c12,c21,c22,c23,c24,c31,c32,c33,c34,c35,c36; f cp_21,cp_22,cp_23,cp_24,cp_31,cp_32,cp_33,cp_34,cp_35,cp_36; f cm_21,cm_22,cm_23,cm_24,cm_31,cm_32,cm_33,cm_34,cm_35,cm_36; f prop,em,vertl,vertu,vertd,vert4r,vert4u,vert4d,int,gd,gp,gm,btxx; f ac0,ac11,ac12,ac21,ac22,ac23,ac24,ac31,ac32,ac33,ac34,ac35,ac36; * v q,p,p1,p2; v Qp,Qm,pp,qp; * s wi; s sing; s frac; s zero; s s,ss,S,x2,X,Y,y; s xt,xh; s am,bm; s gramd; s lp,ls,lt; s m,mp,ms,eps(:1); s lepton,quark; s ctws,ctwsm1,one; s dlt,Lnw,Lns,Lnmu; s z,km,QmS,Ys,Xpo,aux; s xm11,xm12,xm21,xm22; s rwm,wm,wms,rhm,rum,rtm,tmr; s fcsx,fcq,fcs,stws,stwc; s pq,p1q,p2q,vpu,vmu,vpd,vmd; s s1,s2,p11,p12,p22,ps,M1,M2,M3; s sf13,sf14,sf15,sf16,sf17,sf18,sf19; s wm,zm,pm,hm,um,dm,tm,g,stw,ctw,Qu,Qd,sr2,q2,trace; s sel1,sel2,sel3,sel4,sel5,sel6,sel7,sel8,sel9,sel10,selM,selm; s sel11,sel12,sel13,sel14,sel15,sel16,sel17,sel18,sel19,sel20,sel21,sel22, sel23,sel24; s sf1,sf2,sf3,sf4,sf5,sf6,sf7,sf8,sf9,sf10,sf11,sf12; s mu2,betaww,betapww,betazz,betapzz,betazg,betapzg,betagg,lkhz,quat,quatc,ctws; cf aux2,aux4; * .global * *------------PROCEDURES------------------------------------------- * * Total fermionic for one iso-doublet * #procedure totf(i,j,l,total) * g 'total'= sum_(lp,17,20,sum_(ls,17,20,sum_(lt,17,20,diagf(lp,ls,lt)))); * #do k1=17,20 #do k2=17,20 #do k3=17,20 id diagf('k1','k2','k3')= vertf'i''j''l''k1''k2''k3'; #enddo #enddo #enddo * #endprocedure * * * Full scalarization * #procedure sclrzng(j,l,dummy) * #call scalarizing3b{'j'|'l'|dummy} #call scalarizing3a{'j'|'l'|dummy} #call scalarizing2{'j'|'l'|dummy} #call scalarizing1{'j'|'l'|dummy} #call scalarb{dummy} * #call mass{dummy} #call sincos{dummy} #call offs{dummy} #call sincos{dummy} * Symmetrize c0 1,3; * #endprocedure * #procedure nsclrzng(j,l,dummy) * #call mass{dummy} #call sincos{dummy} #call offs{dummy} * Symmetrize c0 1,3; * #endprocedure * * T1 - Vector and scalar internal lines * #procedure vertvs(i,j,l) * #do k1=1,8 #do k2=1,8 #do k3=1,8 * g vertvs'i''j''l''k1''k2''k3' = vertl('i','k3','k1')*prop('k3',3)* vertu('l','k2','k3')*prop('k2',2)* vertd('j','k1','k2')*prop('k1',1)*int; nwrite statistics; #enddo #enddo #enddo * #endprocedure * * T1 - Fermion internal lines -- trace is taken against the arrow * #procedure vertf(i,j,l) * #do k1=17,20 #do k2=17,20 g vertf'i''j''l''k1''k2'17 = -vertl('i',17,'k1')*prop(17,3)* vertu('l','k2',17)*prop('k2',2)* vertd('j','k1','k2')*prop('k1',1)*int*trace; g vertf'i''j''l''k1''k2'18 = -vertl('i',18,'k1')*prop('k1',1)* vertd('j','k1','k2')*prop('k2',2)* vertu('l','k2',18)*prop(18,3)*int*trace; g vertf'i''j''l''k1''k2'19 = -vertl('i',19,'k1')*prop(19,3)* vertu('l','k2',19)*prop('k2',2)* vertd('j','k1','k2')*prop('k1',1)*int*trace; g vertf'i''j''l''k1''k2'20 = -vertl('i',20,'k1')*prop('k1',1)* vertd('j','k1','k2')*prop('k2',2)* vertu('l','k2',20)*prop(20,3)*int*trace; nwrite statistics; #enddo #enddo * id qf(mu?)= q(mu); * * 1 \pm gamma^5 to the left * repeat; id gd(ii?,mu?)*gp(ii?)= gm(ii)*gd(ii,mu); id gd(ii?,mu?)*gm(ii?)= gp(ii)*gd(ii,mu); id gp(ii?)*gm(ii?)= 0; al gm(ii?)*gp(ii?)= 0; al gp(ii?)*gp(ii?)= 2*gp(ii); al gm(ii?)*gm(ii?)= 2*gm(ii); endrepeat; * * Reduction of gd(i,q) ... gd(i,q) * repeat; id gd(ii,q)*gd(ii,q)= q2; * id gd(ii,q)*gd(ii,mu?)*gd(ii,nu?)*gd(ii,al?)*gd(ii,q)= -gd(ii,q)*gd(ii,mu)*gd(ii,nu)*gd(ii,q)*gd(ii,al)+ 2*q(al)*gd(ii,q)*gd(ii,mu)*gd(ii,nu); * id gd(ii,q)*gd(ii,mu?)*gd(ii,nu?)*gd(ii,q)= q2*gd(ii,mu)*gd(ii,nu)- 2*q(mu)*gd(ii,q)*gd(ii,nu)+2*q(nu)*gd(ii,q)*gd(ii,mu); * id gd(ii,q)*gd(ii,mu?)*gd(ii,q)= -q2*gd(ii,mu)+2*q(mu)*gd(ii,q); endrepeat; * * Traces * id gp(ii)*gd(ii,m1?)*gd(ii,m2?)*gd(ii,m3?)*gd(ii,m4?)* gd(ii,m5?)*gd(ii,m6?)*trace= 4*d_(m1,m2)*d_(m3,m4)*d_(m5,m6)-4*d_(m1,m2)*d_(m3,m5)*d_(m4,m6)+4* d_(m1,m2)*d_(m3,m6)*d_(m4,m5)-4*d_(m1,m3)*d_(m2,m4)*d_(m5,m6)+4* d_(m1,m3)*d_(m2,m5)*d_(m4,m6)-4*d_(m1,m3)*d_(m2,m6)*d_(m4,m5)+4* d_(m1,m4)*d_(m2,m3)*d_(m5,m6)-4*d_(m1,m4)*d_(m2,m5)*d_(m3,m6)+4* d_(m1,m4)*d_(m2,m6)*d_(m3,m5)-4*d_(m1,m5)*d_(m2,m3)*d_(m4,m6)+4* d_(m1,m5)*d_(m2,m4)*d_(m3,m6)-4*d_(m1,m5)*d_(m2,m6)*d_(m3,m4)+4* d_(m1,m6)*d_(m2,m3)*d_(m4,m5)-4*d_(m1,m6)*d_(m2,m4)*d_(m3,m5)+4* d_(m1,m6)*d_(m2,m5)*d_(m3,m4)+4*e_(m1,m2,m3,m4)*d_(m5,m6)-4*e_(m1,m2 ,m3,m5)*d_(m4,m6)+4*e_(m1,m2,m3,m6)*d_(m4,m5)+4*e_(m1,m4,m5,m6)* d_(m2,m3)-4*e_(m2,m4,m5,m6)*d_(m1,m3)+4*e_(m3,m4,m5,m6)*d_(m1,m2); id gm(ii)*gd(ii,m1?)*gd(ii,m2?)*gd(ii,m3?)*gd(ii,m4?)* gd(ii,m5?)*gd(ii,m6?)*trace= 4*d_(m1,m2)*d_(m3,m4)*d_(m5,m6)-4*d_(m1,m2)*d_(m3,m5)*d_(m4,m6)+4* d_(m1,m2)*d_(m3,m6)*d_(m4,m5)-4*d_(m1,m3)*d_(m2,m4)*d_(m5,m6)+4* d_(m1,m3)*d_(m2,m5)*d_(m4,m6)-4*d_(m1,m3)*d_(m2,m6)*d_(m4,m5)+4* d_(m1,m4)*d_(m2,m3)*d_(m5,m6)-4*d_(m1,m4)*d_(m2,m5)*d_(m3,m6)+4* d_(m1,m4)*d_(m2,m6)*d_(m3,m5)-4*d_(m1,m5)*d_(m2,m3)*d_(m4,m6)+4* d_(m1,m5)*d_(m2,m4)*d_(m3,m6)-4*d_(m1,m5)*d_(m2,m6)*d_(m3,m4)+4* d_(m1,m6)*d_(m2,m3)*d_(m4,m5)-4*d_(m1,m6)*d_(m2,m4)*d_(m3,m5)+4* d_(m1,m6)*d_(m2,m5)*d_(m3,m4)-4*e_(m1,m2,m3,m4)*d_(m5,m6)+4*e_(m1,m2 ,m3,m5)*d_(m4,m6)-4*e_(m1,m2,m3,m6)*d_(m4,m5)-4*e_(m1,m4,m5,m6)* d_(m2,m3)+4*e_(m2,m4,m5,m6)*d_(m1,m3)-4*e_(m3,m4,m5,m6)*d_(m1,m2); id gd(ii,m1?)*gd(ii,m2?)*gd(ii,m3?)*gd(ii,m4?)*gd(ii,m5?)= 0; id gp(ii)*gd(ii,m1?)*gd(ii,m2?)*gd(ii,m3?)*gd(ii,m4?)*trace= 4*d_(m1,m2)*d_(m3,m4)-4*d_(m1,m3)*d_(m2,m4)+4*d_(m1,m4)*d_(m2,m3)+ 4*e_(m1,m2,m3,m4); id gm(ii)*gd(ii,m1?)*gd(ii,m2?)*gd(ii,m3?)*gd(ii,m4?)*trace= 4*d_(m1,m2)*d_(m3,m4)-4*d_(m1,m3)*d_(m2,m4)+4*d_(m1,m4)*d_(m2,m3)- 4*e_(m1,m2,m3,m4); id gd(ii,m1?)*gd(ii,m2?)*gd(ii,m3?)= 0; id gp(ii)*gd(ii,m1?)*gd(ii,m2?)*trace= 4*d_(m1,m2); id gm(ii)*gd(ii,m1?)*gd(ii,m2?)*trace= 4*d_(m1,m2); id gd(ii,m1?)= 0; * id q.q= q2; id q(mu?)= qf(mu); id pf(q)= pq; id p1f(q)= p1q; id p2f(q)= p2q; id pq^2= pf(de8)*pf(de9)*qf(de8)*qf(de9); id pq= pf(de8)*qf(de8); id p1q^2= p1f(de1)*p1f(de2)*qf(de1)*qf(de2); id p1q*q2= p1f(de1)*qf(de2)*qf(de1)*qf(de2); id p1q= p1f(de1)*qf(de1); id p2q^2= p2f(de3)*p2f(de4)*qf(de3)*qf(de4); id p1q*q2= p2f(de3)*qf(de4)*qf(de3)*qf(de4); id p2q= p2f(de3)*qf(de3); al gd(ii,q)= gd(ii,de5)*qf(de5); al e_(q,mu?,al?,be?)= e_(de6,mu,al,be)*qf(de6); al q2^2= qf(de7)*qf(de7)*qf(de8)*qf(de8); al q2= qf(de7)*qf(de7); * *sum de1,de2,de3,de4,de5,de6,de7,de8,de9; * #endprocedure * * PFD of sin/cos * #procedure sincos(dummy) * id stw^2= -ctwsm1; al stw^-2= -1/ctwsm1; al ctw^-2= 1/ctws; ratio,ctwsm1,ctws,one; id one= 1; al one^-1= 1; id ctws= ctw^2; al ctws^-1= 1/ctw^2; al ctwsm1= -stw^2; al ctwsm1^-1= -1/stw^2; * #endprocedure * * Off shell j,l * #procedure offs(dummy) * id sr2^-2= 1/2; id p1.p1= p11; al p2.p2= p22; al p.p1= -p11-p1.p2; al p.p2= -p22-p1.p2; al p1.p2= 1/2*(-s-p11-p22); al p.p= ps; id ps= -s; al ps^-1= -1/s; *id zm= wm/ctw; *id wm^2= rwm^2*s; *al hm^2= rhm^2*s; *al wm^-2= 1/rwm^2/s; *al wm= rwm*ss; *al wm^-1= 1/rwm/ss; al um^2= rum^2*s; * id ss^2= s; al ss^-2= 1/s; al s^-1*ss= 1/ss; al s^2*ss^-1= s*ss; * #endprocedure * * Total fermionic(third family) * #procedure tfms(tf,totf) * *g 'totf'= 'tf'*(lepton+3*quark); g 'totf'= 'tf'*(3*quark); * id rum^m?= frum(m); id frum(0)= 1; id frum(m?)*lepton= 0; id frum(m?)*quark= rtm^m*quark; id a0(um)*lepton= 0; al a0(um)*quark= a0(tm)*quark; * id vpu= 1-2*Qu*stw^2; al vmu= -2*Qu*stw^2; al vpd= -1-2*Qd*stw^2; al vmd= -2*Qd*stw^2; * id Qu^3*lepton= 0; al Qu^2*lepton= 0; al Qu*lepton= 0; id Qu^3*quark= 8/27; al Qu^2*quark= 4/9; al Qu*quark= 2/3; * id Qd^3*lepton= -1; al Qd^2*lepton= 1; al Qd*lepton= -1; id Qd^3*quark= -1/27; al Qd^2*quark= 1/9; al Qd*quark= -1/3; * id lepton= 1; al quark= 1; * Symmetrize b0 2,3; id a0(0)= 0; id a0(zero)= 0; id ctw^2= 1-stw^2; * #endprocedure * * Total fermionic(one massless family) * #procedure tfml(tf,totf) * g 'totf'= 'tf'*(lepton+3*quark); * id rum= 0; argument; id um= zero; endargument; * id vpu= 1-2*Qu*stw^2; al vmu= -2*Qu*stw^2; al vpd= -1-2*Qd*stw^2; al vmd= -2*Qd*stw^2; * id Qu^3*lepton= 0; al Qu^2*lepton= 0; al Qu*lepton= 0; id Qu^3*quark= 8/27; al Qu^2*quark= 4/9; al Qu*quark= 2/3; * id Qd^3*lepton= -1; al Qd^2*lepton= 1; al Qd*lepton= -1; id Qd^3*quark= -1/27; al Qd^2*quark= 1/9; al Qd*quark= -1/3; * id lepton= 1; al quark= 1; * Symmetrize b0 2,3; id a0(0)= 0; id a0(zero)= 0; id ctw^2= 1-stw^2; * #endprocedure * * Masses * #procedure mass(dummy) * id dm= 0; al pm= 0; argument; id dm= zero; endargument; id zm= wm/ctw; al zm^-1= ctw/wm; * #endprocedure * * External masses * #procedure exmasses(am,bm) * id p11= -'am'^2; al p22= -'bm'^2; al p12= 1/2*(-s+'am'^2+'bm'^2); id p11^-1= -1/'am'^2; al p22^-1= -1/'bm'^2; * #endprocedure * * eps X Singular parts * #procedure sing(dummy) * id n^2= 16+8*eps; id n= 4+eps; id eps*c36(M1?,M2?,M3?)= 1/6; id eps*c35(M1?,M2?,M3?)= 1/3; id eps*c24(M1?,M2?,M3?)= -1/2; id eps*b22(ps?,M1?,M2?)= 1/2*(M1^2+M2^2+ps/3); id eps*b21(ps?,M1?,M2?)= -2/3; id eps*b1(ps?,M1?,M2?)= 1; id eps*b0(ps?,M1?,M2?)= -2; id eps= 0; * #endprocedure * * Massless fermion B0 functions * #procedure massless(dummy); * id b0(wms,1,1,0,0)= dlt+2-Lnw; id b0(p11,1,1,0,0)= dlt+2-Lnw; al b0(p22,1,1,0,0)= dlt+2-Lnw; al b0(ps,1,1,0,0)= dlt+2-Lns; * id Lnw= Lnmu+Lns; * #endprocedure * * Integration * #procedure reduction(dummy) * id kl(mu?)= qf(mu)+p1f(mu)+p2f(mu); al ku(mu?)= -qf(mu)-p1f(mu)-p2f(mu); al qu(mu?)= qf(mu)+p1f(mu); al kd(mu?)= -qf(mu)-p1f(mu); al pd(mu?)= qf(mu); al ql(mu?)= -qf(mu); * id qf(mu?)*qf(nu?)*qf(al?)*den1(m?)*den2(mp?)*den3(ms?)*int= c31(m,mp,ms)*p1f(mu)*p1f(nu)*p1f(al)+ c32(m,mp,ms)*p2f(mu)*p2f(nu)*p2f(al)+ c33(m,mp,ms)*sym211(mu,nu,al)+ c34(m,mp,ms)*sym122(mu,nu,al)+ c35(m,mp,ms)*pd1(mu,nu,al)+ c36(m,mp,ms)*pd2(mu,nu,al); id sym211(mu?,nu?,al?)= p2f(mu)*p1f(nu)*p1f(al)+p1f(mu)*p2f(nu)*p1f(al)+ p1f(mu)*p1f(nu)*p2f(al); id sym122(mu?,nu?,al?)= p1f(mu)*p2f(nu)*p2f(al)+p2f(mu)*p1f(nu)*p2f(al)+ p2f(mu)*p2f(nu)*p1f(al); id pd1(mu?,nu?,al?)= p1f(mu)*d_(nu,al)+p1f(nu)*d_(mu,al)+p1f(al)*d_(mu,nu); id pd2(mu?,nu?,al?)= p2f(mu)*d_(nu,al)+p2f(nu)*d_(mu,al)+p2f(al)*d_(mu,nu); * id qf(mu?)*qf(nu?)*den1(m?)*den2(mp?)*den3(ms?)*int= c21(m,mp,ms)*p1f(mu)*p1f(nu)+ c22(m,mp,ms)*p2f(mu)*p2f(nu)+ c23(m,mp,ms)*(p1f(mu)*p2f(nu)+p1f(nu)*p2f(mu))+ c24(m,mp,ms)*d_(mu,nu); * id qf(mu?)*den1(m?)*den2(mp?)*den3(ms?)*int= c11(m,mp,ms)*p1f(mu)+ c12(m,mp,ms)*p2f(mu); * id den1(m?)*den2(mp?)*den3(ms?)*int= c0(m,mp,ms); * id qf(mu?)*qf(nu?)*den1(m?)*den3(ms?)*int= b21(ps,m,ms)*(p1f(mu)+p2f(mu))*(p1f(nu)+p2f(nu))+ b22(ps,m,ms)*d_(mu,nu); id qf(mu?)*den1(m?)*den3(ms?)*int= b1(ps,m,ms)*(p1f(mu)+p2f(mu)); id den1(m?)*den3(ms?)*int= b0(ps,m,ms); * id qf(mu?)*qf(nu?)*den1(m?)*den2(mp?)*int= b21(p11,m,mp)*p1f(mu)*p1f(nu)+ b22(p11,m,mp)*d_(mu,nu); id qf(mu?)*den1(m?)*den2(mp?)*int= b1(p11,m,mp)*p1f(mu); id den1(m?)*den2(mp?)*int= b0(p11,m,mp); * * q -> q-p1 * id qf(mu?)*qf(nu?)*den2(mp?)*den3(ms?)= (qf(mu)-p1f(mu))*(qf(nu)-p1f(nu))*den2p(mp)*den3p(ms); id qf(mu?)*den2(mp?)*den3(ms?)= (qf(mu)-p1f(mu))*den2p(mp)*den3p(ms); id den2(mp?)*den3(ms?)= den2p(mp)*den3p(ms); * id qf(mu?)*qf(nu?)*den2p(mp?)*den3p(ms?)*int= b21(p22,mp,ms)*p2f(mu)*p2f(nu)+ b22(p22,mp,ms)*d_(mu,nu); id qf(mu?)*den2p(mp?)*den3p(ms?)*int= b1(p22,mp,ms)*p2f(mu); id den2p(mp?)*den3p(ms?)*int= b0(p22,mp,ms); * id pf(mu?)= p(mu); al p1f(mu?)= p1(mu); al p2f(mu?)= p2(mu); * #endprocedure * * 7 Independent form factors for physical sources * #procedure formfact(dummy) * id p1(be)= 0; al p(be)= -p2(be); id p2(al)= 0; al p1(al)= -p(al); id p(mu)= 0; al p2(mu)= -p1(mu); * id p(al)*p1(mu)*p2(be)= 1/2*w(4)*s*ss; id d_(al,be)*p1(mu)= 1/4*ss*w(1)-d_(mu,be)*p(al)-d_(mu,al)*p2(be); * id p(be)= -p2(be); al p2(mu)= -p1(mu); id d_(mu,be)*p(al)= 1/4*(w(2)+w(3))*ss; id d_(mu,al)*p2(be)= 1/4*(w(2)-w(3))*ss; * id e_(p2,mu,al,be)= 1/2*(w(5)+w(6))*ss; id e_(p1,mu,al,be)= 1/2*(w(5)-w(6))*ss; id e_(p1,p2,al,be)*p1(mu)= -w(7)*s*ss; id e_(p1,p2,al,be)*p2(mu)= w(7)*s*ss; * * After Schouten identity * id e_(p1,p2,mu,al)*p2(be)= -1/2*w(5)*ss*(p12-p22)- 1/2*w(6)*ss*(p12+p22)-w(7)*s*ss; id e_(p1,p2,mu,be)*p(al)= -1/2*w(5)*ss*(p11-p12)- 1/2*w(6)*ss*(p11+p12)+w(7)*s*ss; * id p(al)*p2(be)= ws(1)*s; al d_(al,be)= ws(2); al e_(p1,p2,al,be)= ws(3)*s; * id p1(mu)*p2(be)= wh(1)*s; al d_(mu,be)= wh(2); al e_(p1,p2,mu,be)= wh(3)*s; * #endprocedure * * 10 Form factors for unphysical i -- p(mu).not.0 * #procedure formfactwin(dummy) * id p1(be)= 0; al p(be)= -p2(be); id p2(al)= 0; al p1(al)= -p(al); al p2(mu)= -p(mu)-p1(mu); * id p(al)*p1(mu)*p2(be)= -1/2*w(4)*(X*y*S)*ss; id p(al)*p(mu)*p2(be)= -1/2*w(8)*(X*y*S)*ss; id d_(al,be)*p(mu)= 1/2*w(9)*ss; id d_(al,be)*p1(mu)= 1/4*ss*w(1)-d_(mu,be)*p(al)-d_(mu,al)*p2(be); * id p(be)= -p2(be); al p2(mu)= -p(mu)-p1(mu); id d_(mu,be)*p(al)= 1/4*(w(2)+w(3))*ss; id d_(mu,al)*p2(be)= 1/4*(w(2)-w(3))*ss; * id p(mu)= -p1(mu)-p2(mu); id e_(p2,mu,al,be)= 1/2*(w(5)+w(6))*ss; id e_(p1,mu,al,be)= 1/2*(w(5)-w(6))*ss; id e_(p1,p2,al,be)*p1(mu)= -w(10)*(X*y*S)*ss; id e_(p1,p2,al,be)*p2(mu)= -w(7)*(X*y*S)*ss; * * After Schouten identity * id e_(p1,p2,mu,al)*p2(be)= -1/2*w(5)*ss*(p12-p22)- 1/2*w(6)*ss*(p12+p22)+w(7)*(X*y*S)*ss; id e_(p1,p2,mu,be)*p(al)= -1/2*w(5)*ss*(p11-p12)- 1/2*w(6)*ss*(p11+p12)+w(10)*(X*y*S)*ss; * id p(al)*p2(be)= ws(1)*s; al d_(al,be)= ws(2); al e_(p1,p2,al,be)= ws(3)*s; * id p2(mu)= -p(mu)-p1(mu); id p1(mu)*p2(be)= wh(1)*s; al d_(mu,be)= wh(2); al p(mu)*p2(be)= wh(3)*s; * #endprocedure * #procedure formfactwinM(dummy) * id p1(mu)*p1(al)*p1(be)= -1/2*w(11)*(X*y*S)*ss; id p1(mu)*p1(al)*p2(be)= -1/2*w(12)*(X*y*S)*ss; id p1(mu)*p1(be)*p2(al)= -1/2*w(13)*(X*y*S)*ss; id p1(mu)*p2(al)*p2(be)= -1/2*w(14)*(X*y*S)*ss; id p1(al)*p1(be)*p2(mu)= -1/2*w(15)*(X*y*S)*ss; id p1(al)*p2(mu)*p2(be)= -1/2*w(16)*(X*y*S)*ss; id p1(be)*p2(mu)*p2(al)= -1/2*w(17)*(X*y*S)*ss; id p2(mu)*p2(al)*p2(be)= -1/2*w(18)*(X*y*S)*ss; id d_(mu,al)*p1(be)= 1/2*w(19)*ss; id d_(mu,be)*p1(al)= 1/2*w(20)*ss; id d_(mu,be)*p2(al)= 1/2*w(21)*ss; id d_(al,be)*p2(mu)= 1/2*w(22)*ss; id e_(p1,p2,mu,al)*p1(be)= e_(p1,p2,mu,be)*p1(al)-w(23)*(X*y*S)*ss; * id p(al)*p1(mu)*p2(be)= -1/2*w(4)*(X*y*S)*ss; id p(al)*p(mu)*p2(be)= -1/2*w(8)*(X*y*S)*ss; id d_(al,be)*p(mu)= 1/2*w(9)*ss; id d_(al,be)*p1(mu)= 1/4*ss*w(1)-d_(mu,be)*p(al)-d_(mu,al)*p2(be); * al p2(mu)= -p(mu)-p1(mu); id d_(mu,be)*p(al)= 1/4*(w(2)+w(3))*ss; id d_(mu,al)*p2(be)= 1/4*(w(2)-w(3))*ss; * id p(mu)= -p1(mu)-p2(mu); id e_(p2,mu,al,be)= 1/2*(w(5)+w(6))*ss; id e_(p1,mu,al,be)= 1/2*(w(5)-w(6))*ss; id e_(p1,p2,al,be)*p1(mu)= -w(10)*(X*y*S)*ss; id e_(p1,p2,al,be)*p2(mu)= -w(7)*(X*y*S)*ss; * id e_(p1,p2,mu,be)*p1(al)= -e_(p1,p2,mu,be)*p2(al)+e_(p1,p2,mu,be)*p(al); id e_(p1,p2,mu,be)*p2(al)= -w(24)*(X*y*S)*ss; * * After Schouten identity * id e_(p1,p2,mu,al)*p2(be)= -1/2*w(5)*ss*(p12-p22)- 1/2*w(6)*ss*(p12+p22)+w(7)*(X*y*S)*ss; id e_(p1,p2,mu,be)*p(al)= -1/2*w(5)*ss*(p11-p12)- 1/2*w(6)*ss*(p11+p12)+w(10)*(X*y*S)*ss; * id d_(mu,be)= ws(1); id p1(mu)*p2(be)= -ws(2)*(X*y*S); id e_(p1,p2,mu,be)= -ws(3)*(X*y*S); id p1(mu)*p1(be)= -ws(4)*(X*y*S); id p1(be)*p2(mu)= -ws(5)*(X*y*S); id p2(mu)*p2(be)= -ws(6)*(X*y*S); * id d_(mu,al)= wss(1); id p1(mu)*p2(al)= -wss(2)*(X*y*S); id e_(p1,p2,mu,al)= -wss(3)*(X*y*S); id p1(mu)*p1(al)= -wss(4)*(X*y*S); id p1(al)*p2(mu)= -wss(5)*(X*y*S); id p2(mu)*p2(al)= -wss(6)*(X*y*S); * id p1(mu)= wb(1)*ss; id p2(mu)= wb(2)*ss; * #endprocedure * #procedure formfactwic(dummy) * id p1(be)= 0; al p(be)= -p2(be); al p1(al)= -p2(al)-p(al); id p(mu)= 0; al p2(mu)= -p1(mu); * id p(al)*p1(mu)*p2(be)= 1/2*w(4)*s*ss; id p2(al)*p1(mu)*p2(be)= 1/2*w(8)*s*ss; id d_(mu,be)*p2(al)= 1/2*w(9)*ss; id d_(al,be)*p1(mu)= 1/4*ss*w(1)-d_(mu,be)*p(al)-d_(mu,al)*p2(be); * id p(be)= -p2(be); al p2(mu)= -p1(mu); id d_(mu,be)*p(al)= 1/4*(w(2)+w(3))*ss; id d_(mu,al)*p2(be)= 1/4*(w(2)-w(3))*ss; * id e_(p2,mu,al,be)= 1/2*(w(5)+w(6))*ss; id e_(p1,mu,al,be)= 1/2*(w(5)-w(6))*ss; id e_(p1,p2,al,be)*p2(mu)= w(7)*s*ss; id e_(p1,p2,mu,be)*p2(al)= w(10)*s*ss; * * After Schouten identity * id e_(p1,p2,mu,al)*p2(be)= -1/2*ss*w(5)*(p12-p22)- 1/2*ss*w(6)*(p12+p22)-s*ss*(w(7)-w(10)); id e_(p1,p2,mu,be)*p(al)= -1/2*ss*w(5)*(p11-p12)- 1/2*ss*w(6)*(p11+p12)+s*ss*(w(7)-w(10)); * id p1(mu)*p2(be)= ws(1)*s; al d_(mu,be)= ws(2); al e_(p1,p2,mu,be)= ws(3)*s; * #endprocedure * * Reduction to scalar integrals of c1 is performed * #procedure scalarizing1(j,l,dummy) * id c11(M1?,M2?,M3?)= xm11*R1(M1,M2,M3)+ xm12*R2(M1,M2,M3); al c12(M1?,M2?,M3?)= xm21*R1(M1,M2,M3)+ xm22*R2(M1,M2,M3); * id R1(M1?,M2?,M3?)= 1/2*(f1(M1,M2,M3)* c0(M1,M2,M3)+ b0(ps,M1,M3)-b0(p22,M2,M3)); al R2(M1?,M2?,M3?)= 1/2*(f2(M1,M2,M3)* c0(M1,M2,M3)+ b0(p11,M1,M2)-b0(ps,M1,M3)); * id f1(M1?,M2?,M3?)= M1^2-M2^2-p11; al f2(M1?,M2?,M3?)= M2^2-M3^2-ps+p11; id xm11= p22*gram('j','l'); al xm12= -p12*gram('j','l'); al xm21= -p12*gram('j','l'); al xm22= p11*gram('j','l'); id pm= 0; * #endprocedure * * Reduction to scalar integrals of c3 is performed * #procedure scalarizing3a(j,l,dummy) * id c35(M1?,M2?,M3?)= xm11*R11(M1,M2,M3)+ xm12*R15(M1,M2,M3); id c36(M1?,M2?,M3?)= xm21*R11(M1,M2,M3)+ xm22*R15(M1,M2,M3); * id R11(M1?,M2?,M3?)= 1/2*(f1(M1,M2,M3)* c24(M1,M2,M3)+ b22(ps,M1,M3)-b22(p22,M2,M3)); id R15(M1?,M2?,M3?)= 1/2*(f2(M1,M2,M3)* c24(M1,M2,M3)+ b22(p11,M1,M2)-b22(ps,M1,M3)); * id c24(M1?,M2?,M3?)= 1/4-1/2*M1^2*c0(M1,M2,M3)+ 1/4*(b0(p22,M2,M3)- f1(M1,M2,M3)*c11(M1,M2,M3)- f2(M1,M2,M3)*c12(M1,M2,M3)); * id f1(M1?,M2?,M3?)= M1^2-M2^2-p11; al f2(M1?,M2?,M3?)= M2^2-M3^2-ps+p11; id xm11= p22*gram('j','l'); al xm12= -p12*gram('j','l'); al xm21= -p12*gram('j','l'); al xm22= p11*gram('j','l'); id pm= 0; * #endprocedure * * Reduction to scalar integrals of c3 is performed * #procedure scalarizing3b(j,l,dummy) * id c31(M1?,M2?,M3?)= xm11*R8(M1,M2,M3)+ xm12*R12(M1,M2,M3); id c33(M1?,M2?,M3?)= xm21*R8(M1,M2,M3)+ xm22*R12(M1,M2,M3); id c34(M1?,M2?,M3?)= xm11*R9(M1,M2,M3)+ xm12*R13(M1,M2,M3); id c32(M1?,M2?,M3?)= xm21*R9(M1,M2,M3)+ xm22*R13(M1,M2,M3); * id R8(M1?,M2?,M3?)= 1/2*(f1(M1,M2,M3)* c21(M1,M2,M3)+b21(ps,M1,M3)-b0(p22,M2,M3))- 2*c35(M1,M2,M3); al R9(M1?,M2?,M3?)= 1/2*(f1(M1,M2,M3)* c22(M1,M2,M3)+b21(ps,M1,M3)-b21(p22,M2,M3)); id R12(M1?,M2?,M3?)= 1/2*(f2(M1,M2,M3)* c21(M1,M2,M3)+b21(p11,M1,M2)-b21(ps,M1,M3)); id R13(M1?,M2?,M3?)= 1/2*(f2(M1,M2,M3)* c22(M1,M2,M3)-b21(ps,M1,M3))- 2*c36(M1,M2,M3); * id f1(M1?,M2?,M3?)= M1^2-M2^2-p11; al f2(M1?,M2?,M3?)= M2^2-M3^2-ps+p11; id xm11= p22*gram('j','l'); al xm12= -p12*gram('j','l'); al xm21= -p12*gram('j','l'); al xm22= p11*gram('j','l'); id pm= 0; * #endprocedure * * Reduction to scalar integrals of c2 is performed * #procedure scalarizing2(j,l,dummy) * id c21(M1?,M2?,M3?)= xm11*R3(M1,M2,M3)+ xm12*R5(M1,M2,M3); id c23(M1?,M2?,M3?)= xm21*R3(M1,M2,M3)+ xm22*R5(M1,M2,M3); id c22(M1?,M2?,M3?)= xm21*R4(M1,M2,M3)+ xm22*R6(M1,M2,M3); * id R3(M1?,M2?,M3?)= 1/2*(f1(M1,M2,M3)* c11(M1,M2,M3)+ b1(ps,M1,M3)+b0(p22,M2,M3))- c24(M1,M2,M3); al R4(M1?,M2?,M3?)= 1/2*(f1(M1,M2,M3)* c12(M1,M2,M3)+ b1(ps,M1,M3)-b1(p22,M2,M3)); al R5(M1?,M2?,M3?)= 1/2*(f2(M1,M2,M3)* c11(M1,M2,M3)+ b1(p11,M1,M2)-b1(ps,M1,M3)); al R6(M1?,M2?,M3?)= 1/2*(f2(M1,M2,M3)* c12(M1,M2,M3)- b1(ps,M1,M3))-c24(M1,M2,M3); * id c24(M1?,M2?,M3?)= 1/4-1/2*M1^2*c0(M1,M2,M3)+ 1/4*(b0(p22,M2,M3)- f1(M1,M2,M3)*c11(M1,M2,M3)- f2(M1,M2,M3)*c12(M1,M2,M3)); * id f1(M1?,M2?,M3?)= M1^2-M2^2-p11; al f2(M1?,M2?,M3?)= M2^2-M3^2-ps+p11; id xm11= p22*gram('j','l'); al xm12= -p12*gram('j','l'); al xm21= -p12*gram('j','l'); al xm22= p11*gram('j','l'); id pm= 0; * #endprocedure * * Reduction to scalar integrals of b2i,b1i * #procedure scalarb(dummy) * id b22(ps?,M1?,M2?)= -ps*b21(ps,M1,M2)+1/2*(a0(M2)+ (M1^2-M2^2-ps)*b1(ps,M1,M2)); id b1(ps?,M1?,M1?) = -1/2*b0(ps,M1,M1); id b21(ps?,M1?,M2?)= ((3*(M1^2+M2^2)+ps)/18+ 1/3/ps*(M1^2-M2^2-ps)*a0(M1)+ 1/3/ps*(M2^2-M1^2+2*ps)*a0(M2)+ 1/3*(ps+2*M2^2-M1^2+ (M1^2-M2^2)^2/ps)*b0(ps,M1,M2))/ps; id b1(ps?,M1?,M2?)= 1/2*(a0(M1)-a0(M2)-(ps+M2^2-M1^2)* b0(ps,M1,M2))/ps; id pm= 0; * #endprocedure * *----------TABLES------------------------------ * Table relax vert4r(1:8,1:8,1:8,1:8); * #do k1=1,8 #do k2=1,8 #do k3=1,8 #do k4=1,8 * Fill vert4r('k1','k2','k3','k4')= 0; * #enddo #enddo #enddo #enddo * * XXAA * Fill vert4r(3,1,1,3)= -g^2*stw^2*(2*d_(al,be)*d_(mup,mdo)- d_(al,mdo)*d_(be,mup)- d_(al,mup)*d_(be,mdo)); Fill vert4r(7,1,1,7)= -2*g^2*stw^2*d_(al,be); * * XXZZ * Fill vert4r(3,2,2,3)= -g^2*ctw^2*(2*d_(al,be)*d_(mup,mdo)- d_(al,mdo)*d_(be,mup)- d_(al,mup)*d_(be,mdo)); Fill vert4r(7,2,2,7)= -g^2*(2*ctw^2-1)^2/2/ctw^2*d_(al,be); * * XXZA * Fill vert4r(3,1,2,3)= -g^2*stw*ctw*(2*d_(al,be)*d_(mup,mdo)- d_(al,mdo)*d_(be,mup)- d_(al,mup)*d_(be,mdo)); Fill vert4r(3,2,1,3)= -g^2*stw*ctw*(2*d_(al,be)*d_(mup,mdo)- d_(al,mdo)*d_(be,mup)- d_(al,mup)*d_(be,mdo)); Fill vert4r(7,1,2,7)= -g^2*stw/ctw*(2*ctw^2-1)*d_(al,be); Fill vert4r(7,2,1,7)= -g^2*stw/ctw*(2*ctw^2-1)*d_(al,be); * * XXWW * Fill vert4r(3,3,4,3)= g^2*(2*d_(mdo,be)*d_(al,mup)- d_(mdo,al)*d_(be,mup)- d_(mdo,mup)*d_(al,be)); Fill vert4r(7,3,4,7)= -1/2*g^2*d_(al,be); * * XXWphi * Fill vert4r(2,3,8,5)= i_*g^2/2*stw^2/ctw*d_(mdo,be); Fill vert4r(4,3,8,8)= -g^2/2*d_(mdo,be); * * XXZH * Fill vert4r(5,2,5,2)= -g^2/2/ctw^2*d_(mup,be); * Fill vert4r(7,2,5,3)= i_*g^2*stw^2/ctw*d_(mup,be); * Fill vert4r(8,2,5,4)= -i_*g^2*stw^2/ctw*d_(mup,be); * Table relax vert4u(1:8,1:8,1:8,1:8); * #do k1=1,8 #do k2=1,8 #do k3=1,8 #do k4=1,8 * Fill vert4u('k1','k2','k3','k4')= 0; * #enddo #enddo #enddo #enddo * * XWAX * Fill vert4u(3,3,1,1)= -g^2*stw^2*(2*d_(mu,sl)*d_(be,sdo)- d_(mu,sdo)*d_(be,sl)- d_(mu,be)*d_(sl,sdo)); Fill vert4u(3,3,1,2)= -g^2*stw*ctw*(2*d_(mu,sl)*d_(be,sdo)- d_(mu,sdo)*d_(be,sl)- d_(mu,be)*d_(sl,sdo)); Fill vert4u(7,3,1,6)= g^2/2*stw*d_(mu,be); Fill vert4u(7,3,1,5)= i_/2*g^2*stw*d_(mu,be); * * XWZX * Fill vert4u(3,3,2,2)= -g^2*ctw^2*(2*d_(mu,sl)*d_(be,sdo)- d_(mu,sdo)*d_(be,sl)- d_(mu,be)*d_(sl,sdo)); Fill vert4u(3,3,2,1)= -g^2*stw*ctw*(2*d_(mu,sl)*d_(be,sdo)- d_(mu,sdo)*d_(be,sl)- d_(mu,be)*d_(sl,sdo)); Fill vert4u(7,3,2,6)= -g^2/2*stw^2/ctw*d_(mu,be); Fill vert4u(7,3,2,5)= -i_/2*g^2*stw^2/ctw*d_(mu,be); * * XWphi0X * Fill vert4u(1,3,6,8)= g^2/2*stw*d_(be,sdo); Fill vert4u(2,3,6,8)= -g^2/2*stw^2/ctw*d_(be,sdo); * * XAAX * Fill vert4u(3,1,1,3)= -g^2*stw^2*(d_(mu,be)*d_(sl,sdo)- d_(mu,sl)*d_(be,sdo)-d_(mu,sdo)*d_(be,sl)); Fill vert4u(7,1,1,7)= -2*g^2*stw^2*d_(mu,be); * * XAZX * Fill vert4u(3,1,2,3)= -g^2*stw*ctw*(d_(mu,be)*d_(sl,sdo)- d_(mu,sl)*d_(be,sdo)-d_(mu,sdo)*d_(be,sl)); Fill vert4u(7,1,2,7)= -g^2*stw/ctw*(2*ctw^2-1)*d_(mu,be); * * XZAX * Fill vert4u(3,2,1,3)= -g^2*stw*ctw*(d_(mu,be)*d_(sl,sdo)- d_(mu,sl)*d_(be,sdo)-d_(mu,sdo)*d_(be,sl)); Fill vert4u(7,2,1,7)= -g^2*stw/ctw*(2*ctw^2-1)*d_(mu,be); * * XZZX * Fill vert4u(3,2,2,3)= -g^2*ctw^2*(d_(mu,be)*d_(sl,sdo)- d_(mu,sl)*d_(be,sdo)-d_(mu,sdo)*d_(be,sl)); Fill vert4u(7,2,2,7)= -g^2/2/ctw^2*(2*ctw^2-1)^2*d_(mu,be); Fill vert4u(6,2,2,6)= -g^2/2/ctw^2*d_(mu,be); * * XAphi0X * Fill vert4u(8,1,6,4)= g^2/2*stw*d_(be,sl); Fill vert4u(7,1,6,3)= g^2/2*stw*d_(be,sl); * * XAHX * Fill vert4u(8,1,5,4)= -i_*g^2/2*stw*d_(be,sl); Fill vert4u(7,1,5,3)= i_*g^2/2*stw*d_(be,sl); * * XZphi0X * Fill vert4u(8,2,6,4)= -g^2/2*stw^2/ctw*d_(be,sl); Fill vert4u(7,2,6,3)= -g^2/2*stw^2/ctw*d_(be,sl); Fill vert4u(2,2,6,6)= -g^2/2/ctw^2*d_(be,sl); * Table relax vert4d(1:8,1:8,1:8,1:8); * #do k1=1,8 #do k2=1,8 #do k3=1,8 #do k4=1,8 * Fill vert4d('k1','k2','k3','k4')= 0; * #enddo #enddo #enddo #enddo * * XAWX * Fill vert4d(4,1,4,1)= -g^2*stw^2*(2*d_(mu,nup)*d_(al,nl)- d_(al,nup)*d_(mu,nl)- d_(nl,nup)*d_(mu,al)); Fill vert4d(4,1,4,2)= -g^2*stw*ctw*(2*d_(mu,nup)*d_(al,nl)- d_(al,nup)*d_(mu,nl)- d_(nl,nup)*d_(mu,al)); Fill vert4d(8,1,4,6)= g^2/2*stw*d_(mu,al); Fill vert4d(8,1,4,5)= i_/2*g^2*stw*d_(mu,al); * * XZWX * Fill vert4d(4,2,4,2)= -g^2*ctw^2*(2*d_(mu,nup)*d_(al,nl)- d_(al,nup)*d_(mu,nl)- d_(nl,nup)*d_(mu,al)); Fill vert4d(4,2,4,1)= -g^2*stw*ctw*(2*d_(mu,nup)*d_(al,nl)- d_(al,nup)*d_(mu,nl)- d_(nl,nup)*d_(mu,al)); Fill vert4d(8,2,4,6)= -g^2/2*stw^2/ctw*d_(mu,al); Fill vert4d(8,2,4,5)= -i_/2*g^2*stw^2/ctw*d_(mu,al); * * XAAX * Fill vert4d(7,1,1,7)= -2*g^2*stw^2*d_(mu,al); Fill vert4d(3,1,1,3)= -g^2*stw^2*(2*d_(mu,al)*d_(nl,nup)- d_(mu,nl)*d_(al,nup)-d_(mu,nup)*d_(al,nl)); * * XAZX * Fill vert4d(7,1,2,7)= -g^2*stw/ctw*(2*ctw^2-1)*d_(mu,al); Fill vert4d(3,1,2,3)= -g^2*stw*ctw*(2*d_(mu,al)*d_(nl,nup)- d_(mu,nl)*d_(al,nup)-d_(mu,nup)*d_(al,nl)); * * XZAX * Fill vert4d(7,2,1,7)= -g^2*stw/ctw*(2*ctw^2-1)*d_(mu,al); Fill vert4d(3,2,1,3)= -g^2*stw*ctw*(2*d_(mu,al)*d_(nl,nup)- d_(mu,nl)*d_(al,nup)-d_(mu,nup)*d_(al,nl)); * * XZZX * Fill vert4d(7,2,2,7)= -g^2/2/ctw^2*(2*ctw^2-1)^2*d_(mu,al); Fill vert4d(3,2,2,3)= -g^2*ctw^2*(2*d_(mu,al)*d_(nl,nup)- d_(mu,nl)*d_(al,nup)-d_(mu,nup)*d_(al,nl)); Fill vert4d(6,2,2,6)= -g^2/2/ctw^2*d_(mu,al); * * Xphi0WX * Fill vert4d(1,6,4,7)= g^2/2*stw*d_(al,nl); Fill vert4d(2,6,4,7)= -g^2/2*stw^2/ctw*d_(al,nl); * * Xphi0AX * Fill vert4d(8,6,1,4)= g^2/2*stw*d_(al,nup); Fill vert4d(7,6,1,3)= g^2/2*stw*d_(al,nup); * * XHAX * Fill vert4d(8,5,1,4)= -i_*g^2/2*stw*d_(al,nup); Fill vert4d(7,5,1,3)= i_*g^2/2*stw*d_(al,nup); * * Xphi0ZX * Fill vert4d(8,6,2,4)= -g^2/2*stw^2/ctw*d_(al,nup); Fill vert4d(7,6,2,3)= -g^2/2*stw^2/ctw*d_(al,nup); Fill vert4d(2,6,2,6)= -g^2/2/ctw^2*d_(al,nup); * * XAphiX * Fill vert4d(8,1,8,1)= -2*g^2*stw^2*d_(mu,nup); Fill vert4d(8,1,8,2)= -g^2*stw/ctw*(2*ctw^2-1)*d_(mu,nup); Fill vert4d(6,1,8,4)= g^2/2*stw*d_(mu,nup); Fill vert4d(5,1,8,4)= i_*g^2/2*stw*d_(mu,nup); * * XZphi0X * Fill vert4d(8,2,8,1)= -g^2*stw/ctw*(2*ctw^2-1)*d_(mu,nup); Fill vert4d(8,2,8,2)= -g^2/2/ctw^2*(2*ctw^2-1)^2*d_(mu,nup); Fill vert4d(6,2,8,4)= -g^2/2*stw^2/ctw*d_(mu,nup); Fill vert4d(5,2,8,4)= -i_*g^2/2*stw^2/ctw*d_(mu,nup); * * XAHX * Fill vert4d(4,1,5,8)= i_*g^2/2*stw*d_(mu,nl); Fill vert4d(3,1,5,7)= -i_*g^2/2*stw*d_(mu,nl); * * XZHX * Fill vert4d(4,2,5,8)= -i_*g^2/2*stw^2/ctw*d_(mu,nl); Fill vert4d(3,2,5,7)= i_*g^2/2*stw^2/ctw*d_(mu,nl); Fill vert4d(2,2,5,5)= -g^2/2/ctw^2*d_(mu,nl); * Table relax vertl(1:8,1:20,1:20); * #do k1=1,8 #do k2=1,20 #do k3=1,20 * Fill vertl('k1','k2','k3')= 0; * #enddo #enddo #enddo * * AXX * Fill vertl(1,3,3)= g*stw*(d_(mu,sl)*(pf(nl)-kl(nl))+ d_(nl,sl)*(kl(mu)-ql(mu))+ d_(mu,nl)*(ql(sl)-pf(sl))); Fill vertl(1,4,4)= -g*stw*(d_(mu,sl)*(pf(nl)-kl(nl))+ d_(nl,sl)*(kl(mu)-ql(mu))+ d_(mu,nl)*(ql(sl)-pf(sl))); * Fill vertl(1,7,3)= -i_*g*stw*wm*d_(mu,nl); Fill vertl(1,3,7)= i_*g*stw*wm*d_(mu,sl); * Fill vertl(1,8,4)= i_*g*stw*wm*d_(mu,nl); Fill vertl(1,4,8)= -i_*g*stw*wm*d_(mu,sl); * Fill vertl(1,7,7)= g*stw*(kl(mu)-ql(mu)); Fill vertl(1,8,8)= g*stw*(ql(mu)-kl(mu)); * Fill vertl(1,9,9)= -g*stw*ql(mu); Fill vertl(1,10,10)= -g*stw*kl(mu); Fill vertl(1,11,11)= g*stw*ql(mu); Fill vertl(1,12,12)= g*stw*kl(mu); * Fill vertl(1,17,17)= i_*g*Qu*stw*gd(ii,mu); Fill vertl(1,18,18)= i_*g*Qu*stw*gd(ii,mu); Fill vertl(1,19,19)= i_*g*Qd*stw*gd(ii,mu); Fill vertl(1,20,20)= i_*g*Qd*stw*gd(ii,mu); * * ZXX * Fill vertl(2,3,3)= g*ctw*(d_(mu,sl)*(pf(nl)-kl(nl))+ d_(nl,sl)*(kl(mu)-ql(mu))+ d_(mu,nl)*(ql(sl)-pf(sl))); Fill vertl(2,4,4)= -g*ctw*(d_(mu,sl)*(pf(nl)-kl(nl))+ d_(nl,sl)*(kl(mu)-ql(mu))+ d_(mu,nl)*(ql(sl)-pf(sl))); * Fill vertl(2,7,3)= i_*g*stw^2/ctw*wm*d_(mu,nl); Fill vertl(2,3,7)= -i_*g*stw^2/ctw*wm*d_(mu,sl); * Fill vertl(2,8,4)= -i_*g*stw^2/ctw*wm*d_(mu,nl); Fill vertl(2,4,8)= i_*g*stw^2/ctw*wm*d_(mu,sl); * Fill vertl(2,7,7)= g*(2*ctw^2-1)/2/ctw*(kl(mu)-ql(mu)); Fill vertl(2,8,8)= g*(2*ctw^2-1)/2/ctw*(ql(mu)-kl(mu)); * Fill vertl(2,2,5)= -g*wm/ctw^2*d_(mu,sl); Fill vertl(2,5,2)= -g*wm/ctw^2*d_(mu,nl); * Fill vertl(2,6,5)= i_/2*g/ctw*(kl(mu)-ql(mu)); Fill vertl(2,5,6)= i_/2*g/ctw*(ql(mu)-kl(mu)); * Fill vertl(2,9,9)= -g*ctw*ql(mu); Fill vertl(2,10,10)= -g*ctw*kl(mu); Fill vertl(2,11,11)= g*ctw*ql(mu); Fill vertl(2,12,12)= g*ctw*kl(mu); * Fill vertl(2,17,17)= i_*g/4/ctw*gd(ii,mu)*(vpu*gp(ii)+vmu*gm(ii)); Fill vertl(2,18,18)= i_*g/4/ctw*gd(ii,mu)*(vpu*gp(ii)+vmu*gm(ii)); Fill vertl(2,19,19)= i_*g/4/ctw*gd(ii,mu)*(vpd*gp(ii)+vmd*gm(ii)); Fill vertl(2,20,20)= i_*g/4/ctw*gd(ii,mu)*(vpd*gp(ii)+vmd*gm(ii)); * * phi0XX * Fill vertl(6,4,8)= g/2*(p(sl)-ql(sl)); Fill vertl(6,8,4)= g/2*(kl(nl)-p(nl)); * Fill vertl(6,7,3)= g/2*(p(nl)-kl(nl)); Fill vertl(6,3,7)= g/2*(ql(sl)-p(sl)); * Fill vertl(6,2,5)= i_/2*g/ctw*(p(sl)-ql(sl)); Fill vertl(6,5,2)= i_/2*g/ctw*(p(nl)-kl(nl)); * Fill vertl(6,5,6)= -g/2*hm^2/wm; Fill vertl(6,6,5)= -g/2*hm^2/wm; * Fill vertl(6,9,9)= i_/2*g*wm; Fill vertl(6,10,10)= i_/2*g*wm; Fill vertl(6,11,11)= -i_/2*g*wm; Fill vertl(6,12,12)= -i_/2*g*wm; * Fill vertl(6,17,17)= i_/4*g*um/wm*(gp(ii)-gm(ii)); Fill vertl(6,18,18)= i_/4*g*um/wm*(gp(ii)-gm(ii)); Fill vertl(6,19,19)= -i_/4*g*dm/wm*(gp(ii)-gm(ii)); Fill vertl(6,20,20)= -i_/4*g*dm/wm*(gp(ii)-gm(ii)); * * HXX * Fill vertl(5,2,2)= -g*wm/ctw^2*d_(nl,sl); Fill vertl(5,3,3)= -g*wm*d_(nl,sl); Fill vertl(5,4,4)= -g*wm*d_(nl,sl); Fill vertl(5,2,6)= i_*g/2/ctw*(ql(sl)-p(sl)); Fill vertl(5,6,2)= i_*g/2/ctw*(kl(nl)-p(nl)); Fill vertl(5,4,8)= i_*g/2*(ql(sl)-p(sl)); Fill vertl(5,7,3)= i_*g/2*(kl(nl)-p(nl)); Fill vertl(5,3,7)= i_*g/2*(ql(sl)-p(sl)); Fill vertl(5,8,4)= i_*g/2*(kl(nl)-p(nl)); Fill vertl(5,6,6)= -g/2*hm^2/wm; Fill vertl(5,7,7)= -g/2*hm^2/wm; Fill vertl(5,8,8)= -g/2*hm^2/wm; Fill vertl(5,9,9)= -g/2*wm; Fill vertl(5,10,10)= -g/2*wm; Fill vertl(5,11,11)= -g/2*wm; Fill vertl(5,12,12)= -g/2*wm; Fill vertl(5,13,13)= -g/2/ctw^2*wm; Fill vertl(5,14,14)= -g/2/ctw^2*wm; Fill vertl(5,17,17)= -g/4*um/wm*(gp(ii)+gm(ii)); Fill vertl(5,18,18)= -g/4*um/wm*(gp(ii)+gm(ii)); Fill vertl(5,19,19)= -g/4*dm/wm*(gp(ii)+gm(ii)); Fill vertl(5,20,20)= -g/4*dm/wm*(gp(ii)+gm(ii)); * Table relax vertu(1:8,1:20,1:20); * #do k1=1,8 #do k2=1,20 #do k3=1,20 * Fill vertu('k1','k2','k3')= 0; * #enddo #enddo #enddo * * WXX * Fill vertu(4,3,1)= g*stw*(d_(mup,nup)*(ku(al)-qu(al))+ d_(nup,al)*(qu(mup)-p2f(mup))+ d_(mup,al)*(p2f(nup)-ku(nup))); Fill vertu(4,3,2)= g*ctw*(d_(mup,nup)*(ku(al)-qu(al))+ d_(nup,al)*(qu(mup)-p2f(mup))+ d_(mup,al)*(p2f(nup)-ku(nup))); * Fill vertu(4,1,4)= g*stw*(d_(mup,nup)*(qu(al)-ku(al))+ d_(mup,al)*(ku(nup)-p2f(nup))+ d_(nup,al)*(p2f(mup)-qu(mup))); Fill vertu(4,2,4)= g*ctw*(d_(mup,nup)*(qu(al)-ku(al))+ d_(mup,al)*(ku(nup)-p2f(nup))+ d_(nup,al)*(p2f(mup)-qu(mup))); * Fill vertu(4,7,1)= -i_*g*stw*wm*d_(mup,al); Fill vertu(4,1,8)= -i_*g*stw*wm*d_(al,nup); * Fill vertu(4,7,2)= i_*g*stw^2/ctw*wm*d_(mup,al); Fill vertu(4,2,8)= i_*g*stw^2/ctw*wm*d_(al,nup); * Fill vertu(4,3,5)= -g*wm*d_(nup,al); Fill vertu(4,5,4)= -g*wm*d_(mup,al); * Fill vertu(4,7,6)= 1/2*g*(ku(al)-qu(al)); Fill vertu(4,6,8)= 1/2*g*(qu(al)-ku(al)); * Fill vertu(4,7,5)= i_/2*g*(qu(al)-ku(al)); Fill vertu(4,5,8)= i_/2*g*(ku(al)-qu(al)); * Fill vertu(4,12,14)= -g*ctw*qu(al); Fill vertu(4,13,11)= -g*ctw*ku(al); * Fill vertu(4,12,16)= -g*stw*qu(al); Fill vertu(4,15,11)= -g*stw*ku(al); * Fill vertu(4,14,10)= g*ctw*qu(al); Fill vertu(4,9,13)= g*ctw*ku(al); * Fill vertu(4,16,10)= g*stw*qu(al); Fill vertu(4,9,15)= g*stw*ku(al); * Fill vertu(4,20,18)= i_*g/2/sr2*gd(ii,al)*gp(ii); Fill vertu(4,17,19)= i_*g/2/sr2*gd(ii,al)*gp(ii); * * AXX * Fill vertu(1,3,3)= g*stw*(d_(al,nup)*(p2f(mup)-qu(mup))+ d_(mup,nup)*(qu(al)-ku(al))+ d_(al,mup)*(ku(nup)-p2f(nup))); Fill vertu(1,4,4)= -g*stw*(d_(al,nup)*(p2f(mup)-qu(mup))+ d_(mup,nup)*(qu(al)-ku(al))+ d_(al,mup)*(ku(nup)-p2f(nup))); *Fill vertu(1,4,4)= g*stw*(d_(al,mup)*(p2f(nup)-ku(nup))+ * d_(mup,nup)*(ku(al)-qu(al))+ * d_(al,nup)*(qu(mup)-p2f(mup))); * Fill vertu(1,7,3)= -i_*g*stw*wm*d_(al,mup); Fill vertu(1,3,7)= i_*g*stw*wm*d_(al,nup); * Fill vertu(1,8,4)= i_*g*stw*wm*d_(al,mup); Fill vertu(1,4,8)= -i_*g*stw*wm*d_(al,nup); * Fill vertu(1,7,7)= g*stw*(qu(al)-ku(al)); Fill vertu(1,8,8)= g*stw*(ku(al)-qu(al)); * Fill vertu(1,9,9)= -g*stw*ku(al); Fill vertu(1,10,10)= -g*stw*qu(al); Fill vertu(1,11,11)= g*stw*ku(al); Fill vertu(1,12,12)= g*stw*qu(al); * Fill vertu(1,17,17)= i_*g*Qu*stw*gd(ii,al); Fill vertu(1,18,18)= i_*g*Qu*stw*gd(ii,al); Fill vertu(1,19,19)= i_*g*Qd*stw*gd(ii,al); Fill vertu(1,20,20)= i_*g*Qd*stw*gd(ii,al); * * ZXX * Fill vertu(2,3,3)= g*ctw*(d_(al,nup)*(p2f(mup)-qu(mup))+ d_(mup,nup)*(qu(al)-ku(al))+ d_(al,mup)*(ku(nup)-p2f(nup))); Fill vertu(2,4,4)= -g*ctw*(d_(al,nup)*(p2f(mup)-qu(mup))+ d_(mup,nup)*(qu(al)-ku(al))+ d_(al,mup)*(ku(nup)-p2f(nup))); *Fill vertu(2,4,4)= g*ctw*(d_(al,mup)*(p2f(nup)-ku(nup))+ * d_(mup,nup)*(ku(al)-qu(al))+ * d_(al,nup)*(qu(mup)-p2f(mup))); * Fill vertu(2,7,3)= i_*g*stw^2/ctw*wm*d_(al,mup); Fill vertu(2,3,7)= -i_*g*stw^2/ctw*wm*d_(al,nup); * Fill vertu(2,8,4)= -i_*g*stw^2/ctw*wm*d_(al,mup); Fill vertu(2,4,8)= i_*g*stw^2/ctw*wm*d_(al,nup); * Fill vertu(2,7,7)= g*(2*ctw^2-1)/2/ctw*(qu(al)-ku(al)); Fill vertu(2,8,8)= g*(2*ctw^2-1)/2/ctw*(ku(al)-qu(al)); * Fill vertu(2,2,5)= -g*wm/ctw^2*d_(al,nup); Fill vertu(2,5,2)= -g*wm/ctw^2*d_(al,mup); * Fill vertu(2,6,5)= i_/2*g/ctw*(ku(al)-qu(al)); Fill vertu(2,5,6)= -i_/2*g/ctw*(ku(al)-qu(al)); * Fill vertu(2,9,9)= -g*ctw*ku(al); Fill vertu(2,10,10)= -g*ctw*qu(al); Fill vertu(2,11,11)= g*ctw*ku(al); Fill vertu(2,12,12)= g*ctw*qu(al); * Fill vertu(2,17,17)= i_*g/4/ctw*gd(ii,al)*(vpu*gp(ii)+vmu*gm(ii)); Fill vertu(2,18,18)= i_*g/4/ctw*gd(ii,al)*(vpu*gp(ii)+vmu*gm(ii)); Fill vertu(2,19,19)= i_*g/4/ctw*gd(ii,al)*(vpd*gp(ii)+vmd*gm(ii)); Fill vertu(2,20,20)= i_*g/4/ctw*gd(ii,al)*(vpd*gp(ii)+vmd*gm(ii)); * * phiXX * Fill vertu(8,3,1)= i_*g*stw*wm*d_(mup,nup); Fill vertu(8,1,4)= i_*g*stw*wm*d_(mup,nup); * Fill vertu(8,3,2)= -i_*g*stw^2/ctw*wm*d_(mup,nup); Fill vertu(8,2,4)= -i_*g*stw^2/ctw*wm*d_(mup,nup); * Fill vertu(8,7,1)= g*stw*(qu(mup)-p2f(mup)); Fill vertu(8,1,8)= g*stw*(ku(nup)-p2f(nup)); * Fill vertu(8,7,2)= g*(2*ctw^2-1)/2/ctw*(qu(mup)-p2f(mup)); Fill vertu(8,2,8)= g*(2*ctw^2-1)/2/ctw*(ku(nup)-p2f(nup)); * Fill vertu(8,6,4)= g/2*(p2f(mup)-qu(mup)); Fill vertu(8,3,6)= g/2*(p2f(nup)-ku(nup)); * Fill vertu(8,5,4)= i_*g/2*(p2f(mup)-qu(mup)); Fill vertu(8,3,5)= i_*g/2*(p2f(nup)-ku(nup)); * Fill vertu(8,5,8)= -g/2*hm^2/wm; Fill vertu(8,7,5)= -g/2*hm^2/wm; * Fill vertu(8,12,14)= i_*g*(2*ctw^2-1)/2/ctw*wm; Fill vertu(8,13,11)= i_*g*(2*ctw^2-1)/2/ctw*wm; * Fill vertu(8,12,16)= i_*g*stw*wm; Fill vertu(8,15,11)= i_*g*stw*wm; * Fill vertu(8,14,10)= -i_*g/2/ctw*wm; Fill vertu(8,9,13)= -i_*g/2/ctw*wm; * Fill vertu(8,20,18)= i_*g/4/sr2*(-(um-dm)/wm*(gp(ii)+gm(ii))+ (um+dm)/wm*(gp(ii)-gm(ii))); Fill vertu(8,17,19)= i_*g/4/sr2*(-(um-dm)/wm*(gp(ii)+gm(ii))+ (um+dm)/wm*(gp(ii)-gm(ii))); * * HXX * Fill vertu(5,2,2)= -g*wm/ctw^2*d_(mup,nup); Fill vertu(5,3,3)= -g*wm*d_(mup,nup); Fill vertu(5,4,4)= -g*wm*d_(mup,nup); Fill vertu(5,5,5)= -3*g/2*hm^2/wm; Fill vertu(5,6,6)= -g/2*hm^2/wm; Fill vertu(5,7,7)= -g/2*hm^2/wm; Fill vertu(5,8,8)= -g/2*hm^2/wm; * Fill vertu(5,6,2)= i_*g/2/ctw*(qu(mup)-p2f(mup)); Fill vertu(5,2,6)= i_*g/2/ctw*(ku(nup)-p2f(nup)); * Fill vertu(5,7,3)= i_*g/2*(qu(mup)-p2f(mup)); Fill vertu(5,3,7)= i_*g/2*(ku(nup)-p2f(nup)); * Fill vertu(5,4,8)= i_*g/2*(ku(nup)-p2f(nup)); Fill vertu(5,8,4)= i_*g/2*(qu(mup)-p2f(mup)); * Fill vertu(5,9,9)= -g/2*wm; Fill vertu(5,10,10)= -g/2*wm; Fill vertu(5,11,11)= -g/2*wm; Fill vertu(5,12,12)= -g/2*wm; Fill vertu(5,13,13)= -g/2*wm/ctw^2; Fill vertu(5,14,14)= -g/2*wm/ctw^2; * Fill vertu(5,17,17)= -g/2*um/wm; Fill vertu(5,18,18)= -g/2*um/wm; Fill vertu(5,19,19)= -g/2*dm/wm; Fill vertu(5,20,20)= -g/2*dm/wm; * Table relax vertd(1:8,1:20,1:20); * #do k1=1,8 #do k2=1,20 #do k3=1,20 * Fill vertd('k1','k2','k3')= 0; * #enddo #enddo #enddo * * WXX * Fill vertd(3,1,3)= g*stw*(d_(mdo,be)*(pd(sdo)-p1f(sdo))+ d_(be,sdo)*(p1f(mdo)-kd(mdo))+ d_(mdo,sdo)*(kd(be)-pd(be))); Fill vertd(3,2,3)= g*ctw*(d_(mdo,be)*(pd(sdo)-p1f(sdo))+ d_(be,sdo)*(p1f(mdo)-kd(mdo))+ d_(mdo,sdo)*(kd(be)-pd(be))); * Fill vertd(3,4,1)= g*stw*(d_(sdo,be)*(kd(mdo)-p1f(mdo))+ d_(be,mdo)*(p1f(sdo)-pd(sdo))+ d_(mdo,sdo)*(pd(be)-kd(be))); Fill vertd(3,4,2)= g*ctw*(d_(sdo,be)*(kd(mdo)-p1f(mdo))+ d_(be,mdo)*(p1f(sdo)-pd(sdo))+ d_(mdo,sdo)*(pd(be)-kd(be))); * Fill vertd(3,1,7)= i_*g*stw*wm*d_(mdo,be); Fill vertd(3,8,1)= i_*g*stw*wm*d_(sdo,be); Fill vertd(3,2,7)= -i_*g*stw^2/ctw*wm*d_(mdo,be); Fill vertd(3,8,2)= -i_*g*stw^2/ctw*wm*d_(sdo,be); * Fill vertd(3,4,5)= -g*wm*d_(mdo,be); Fill vertd(3,5,3)= -g*wm*d_(sdo,be); * Fill vertd(3,6,7)= 1/2*g*(kd(be)-pd(be)); Fill vertd(3,8,6)= 1/2*g*(pd(be)-kd(be)); * Fill vertd(3,5,7)= i_/2*g*(kd(be)-pd(be)); Fill vertd(3,8,5)= i_/2*g*(pd(be)-kd(be)); * Fill vertd(3,11,13)= -g*ctw*kd(be); Fill vertd(3,14,12)= -g*ctw*pd(be); * Fill vertd(3,11,15)= -g*stw*kd(be); Fill vertd(3,16,12)= -g*stw*pd(be); * Fill vertd(3,13,9)= g*ctw*kd(be); Fill vertd(3,10,14)= g*ctw*pd(be); * Fill vertd(3,15,9)= g*stw*kd(be); Fill vertd(3,10,16)= g*stw*pd(be); * Fill vertd(3,19,17)= i_*g/2/sr2*gd(ii,be)*gp(ii); Fill vertd(3,18,20)= i_*g/2/sr2*gd(ii,be)*gp(ii); * * AXX * Fill vertd(1,3,3)= g*stw*(d_(be,mdo)*(p1f(sdo)-pd(sdo))+ d_(sdo,mdo)*(pd(be)-kd(be))+ d_(be,sdo)*(kd(mdo)-p1f(mdo))); Fill vertd(1,4,4)= -g*stw*(d_(be,mdo)*(p1f(sdo)-pd(sdo))+ d_(sdo,mdo)*(pd(be)-kd(be))+ d_(be,sdo)*(kd(mdo)-p1f(mdo))); *Fill vertd(1,4,4)= g*stw*(d_(be,sdo)*(p1f(mdo)-kd(mdo))+ * d_(sdo,mdo)*(kd(be)-pd(be))+ * d_(be,mdo)*(pd(sdo)-p1f(sdo))); * Fill vertd(1,4,8)= -i_*g*stw*wm*d_(be,mdo); Fill vertd(1,8,4)= i_*g*stw*wm*d_(be,sdo); * Fill vertd(1,7,3)= -i_*g*stw*wm*d_(be,sdo); Fill vertd(1,3,7)= i_*g*stw*wm*d_(be,mdo); * Fill vertd(1,7,7)= g*stw*(pd(be)-kd(be)); Fill vertd(1,8,8)= g*stw*(kd(be)-pd(be)); * Fill vertd(1,9,9)= -g*stw*kd(be); Fill vertd(1,10,10)= -g*stw*pd(be); Fill vertd(1,11,11)= g*stw*kd(be); Fill vertd(1,12,12)= g*stw*pd(be); * Fill vertd(1,17,17)= i_*g*Qu*stw*gd(ii,be); Fill vertd(1,18,18)= i_*g*Qu*stw*gd(ii,be); Fill vertd(1,19,19)= i_*g*Qd*stw*gd(ii,be); Fill vertd(1,20,20)= i_*g*Qd*stw*gd(ii,be); * * ZXX * Fill vertd(2,3,3)= g*ctw*(d_(be,mdo)*(p1f(sdo)-pd(sdo))+ d_(sdo,mdo)*(pd(be)-kd(be))+ d_(be,sdo)*(kd(mdo)-p1f(mdo))); Fill vertd(2,4,4)= -g*ctw*(d_(be,mdo)*(p1f(sdo)-pd(sdo))+ d_(sdo,mdo)*(pd(be)-kd(be))+ d_(be,sdo)*(kd(mdo)-p1f(mdo))); *Fill vertd(2,4,4)= g*ctw*(d_(be,sdo)*(p1f(mdo)-kd(mdo))+ * d_(sdo,mdo)*(kd(be)-pd(be))+ * d_(be,mdo)*(pd(sdo)-p1f(sdo))); Fill vertd(2,7,3)= i_*g*stw^2/ctw*wm*d_(be,sdo); Fill vertd(2,3,7)= -i_*g*stw^2/ctw*wm*d_(be,mdo); * Fill vertd(2,8,4)= -i_*g*stw^2/ctw*wm*d_(be,sdo); Fill vertd(2,4,8)= i_*g*stw^2/ctw*wm*d_(be,mdo); * Fill vertd(2,7,7)= g*(2*ctw^2-1)/2/ctw*(pd(be)-kd(be)); Fill vertd(2,8,8)= g*(2*ctw^2-1)/2/ctw*(kd(be)-pd(be)); * Fill vertd(2,2,5)= -g*wm/ctw^2*d_(be,mdo); Fill vertd(2,5,2)= -g*wm/ctw^2*d_(be,sdo); * Fill vertd(2,6,5)= i_/2*g/ctw*(kd(be)-pd(be)); Fill vertd(2,5,6)= -i_/2*g/ctw*(kd(be)-pd(be)); * Fill vertd(2,9,9)= -g*ctw*kd(be); Fill vertd(2,10,10)= -g*ctw*pd(be); Fill vertd(2,11,11)= g*ctw*kd(be); Fill vertd(2,12,12)= g*ctw*pd(be); * Fill vertd(2,17,17)= i_*g/4/ctw*gd(ii,be)*(vpu*gp(ii)+vmu*gm(ii)); Fill vertd(2,18,18)= i_*g/4/ctw*gd(ii,be)*(vpu*gp(ii)+vmu*gm(ii)); Fill vertd(2,19,19)= i_*g/4/ctw*gd(ii,be)*(vpd*gp(ii)+vmd*gm(ii)); Fill vertd(2,20,20)= i_*g/4/ctw*gd(ii,be)*(vpd*gp(ii)+vmd*gm(ii)); * * phiXX * Fill vertd(7,19,17)= i_*g/4/sr2*((um-dm)/wm*(gp(ii)+gm(ii))+ (um+dm)/wm*(gp(ii)-gm(ii))); Fill vertd(7,18,20)= i_*g/4/sr2*((um-dm)/wm*(gp(ii)+gm(ii))+ (um+dm)/wm*(gp(ii)-gm(ii))); * Table relax em(1:8); Fill em(1)= 0; Fill em(2)= zm; Fill em(3)= wm; Fill em(4)= wm; Fill em(5)= hm; Fill em(6)= zm; Fill em(7)= wm; Fill em(8)= wm; * Table relax btxx(1:5,1:5); * Fill btxx(1,1)= betagg; Fill btxx(1,2)= betazg; Fill btxx(2,1)= betazg; Fill btxx(2,2)= betazz; Fill btxx(3,4)= betaww; *Fill btxx(2,5)= lkhz; Fill btxx(2,5)= 1; * Table relax prop(1:20,1:3); * Fill prop(1,1)= den1(pm)*d_(nl,mdo); Fill prop(1,2)= den2(pm)*d_(sdo,nup); Fill prop(1,3)= den3(pm)*d_(mup,sl); * Fill prop(2,1)= den1(zm)*d_(nl,mdo); Fill prop(2,2)= den2(zm)*d_(sdo,nup); Fill prop(2,3)= den3(zm)*d_(mup,sl); * Fill prop(3,1)= den1(wm)*d_(nl,mdo); Fill prop(3,2)= den2(wm)*d_(sdo,nup); Fill prop(3,3)= den3(wm)*d_(mup,sl); * #do k=1,3 Fill prop(4,'k')= prop(3,'k'); #enddo * Fill prop(5,1)= den1(hm); Fill prop(5,2)= den2(hm); Fill prop(5,3)= den3(hm); Fill prop(6,1)= den1(zm); Fill prop(6,2)= den2(zm); Fill prop(6,3)= den3(zm); Fill prop(7,1)= den1(wm); Fill prop(7,2)= den2(wm); Fill prop(7,3)= den3(wm); Fill prop(8,1)= den1(wm); Fill prop(8,2)= den2(wm); Fill prop(8,3)= den3(wm); Fill prop(9,1)= den1(wm); Fill prop(9,2)= den2(wm); Fill prop(9,3)= den3(wm); Fill prop(10,1)= den1(wm); Fill prop(10,2)= den2(wm); Fill prop(10,3)= den3(wm); Fill prop(11,1)= den1(wm); Fill prop(11,2)= den2(wm); Fill prop(11,3)= den3(wm); Fill prop(12,1)= den1(wm); Fill prop(12,2)= den2(wm); Fill prop(12,3)= den3(wm); Fill prop(13,1)= den1(zm); Fill prop(13,2)= den2(zm); Fill prop(13,3)= den3(zm); Fill prop(14,1)= den1(zm); Fill prop(14,2)= den2(zm); Fill prop(14,3)= den3(zm); Fill prop(15,1)= den1(pm); Fill prop(15,2)= den2(pm); Fill prop(15,3)= den3(pm); Fill prop(16,1)= den1(pm); Fill prop(16,2)= den2(pm); Fill prop(16,3)= den3(pm); Fill prop(17,1)= (-i_*gd(ii,la1)*qf(la1)+um)*den1(um); Fill prop(17,2)= (-i_*gd(ii,la2)*(qf(la2)+p1f(la2))+um)*den2(um); Fill prop(17,3)= (-i_*gd(ii,la3)*(qf(la3)+p1f(la3)+p2f(la3))+um)*den3(um); Fill prop(18,1)= (i_*gd(ii,la4)*qf(la4)+um)*den1(um); Fill prop(18,2)= (i_*gd(ii,la5)*(qf(la5)+p1f(la5))+um)*den2(um); Fill prop(18,3)= (i_*gd(ii,la6)*(qf(la6)+p1f(la6)+p2f(la6))+um)*den3(um); Fill prop(19,1)= (-i_*gd(ii,la7)*qf(la7)+dm)*den1(dm); Fill prop(19,2)= (-i_*gd(ii,la8)*(qf(la8)+p1f(la8))+dm)*den2(dm); Fill prop(19,3)= (-i_*gd(ii,la9)*(qf(la9)+p1f(la9)+p2f(la9))+dm)*den3(dm); Fill prop(20,1)= (i_*gd(ii,la10)*qf(la10)+dm)*den1(dm); Fill prop(20,2)= (i_*gd(ii,la11)*(qf(la11)+p1f(la11))+dm)*den2(dm); Fill prop(20,3)= (i_*gd(ii,la12)*(qf(la12)+p1f(la12)+p2f(la12))+dm)*den3(dm); * *sum la1,la2,la3,la4,la5,la6,la7,la8,la9,la10,la11,la12; * .global * *-----MAIN----------------------------------------------------------- * #procedure genesismM(i,j,l) * #call vertf{'i'|'j'|'l'} * #call reduction{dummy} #call sing{dummy} #call formfactwinM{dummy} #call offs{dummy} * .store * #call totf{'i'|'j'|'l'|tf} #call nsclrzng{'j'|'l'|dummy} #call sclrzng{'j'|'l'|dummy} * #call tfml{tf|tfml} * argument; id zero= 0; endargument; id zero= 0; al a0(0)= 0; * al p12= 1/2*(-s-p11-p22); id p11= -s1; al p22= -s2; id s= -X*y*S; al s1= Y*y*S; al s2= x2*y*S; id gram(3,4)= -4/gramd/(y*S)^2; id gram(3,8)= -4/gramd/(y*S)^2; id gram(7,4)= -4/gramd/(y*S)^2; id gram(7,8)= -4/gramd/(y*S)^2; * b g,stw,y,S,ss,w,ws,c0,b0; print; * #endprocedure * #procedure genesisMM(i,j,l) * #call vertf{'i'|'j'|'l'} * #call reduction{dummy} #call sing{dummy} #call formfactwinM{dummy} #call offs{dummy} * .store * #call totf{'i'|'j'|'l'|tf} #call nsclrzng{'j'|'l'|dummy} *#call sclrzng{'j'|'l'|dummy} nwrite statistics; * #call tfms{tf|tfms} nwrite statistics; * argument; id zero= 0; endargument; id zero= 0; al a0(0)= 0; * al p12= 1/2*(-s-p11-p22); id p11= -s1; al p22= -s2; id s= -X*y*S; al s1= Y*y*S; al s2= x2*y*S; id gram(3,4)= -4/gramd/(y*S)^2; id gram(3,8)= -4/gramd/(y*S)^2; id gram(7,4)= -4/gramd/(y*S)^2; id gram(7,8)= -4/gramd/(y*S)^2; * b g,stw,y,S,ss,w,ws,c0,b0,p1,p2,p,d_,e_; print; * #endprocedure * #procedure genesism(i,j,l) * #call vertf{'i'|'j'|'l'} * #call reduction{dummy} #call sing{dummy} *#call formfactwin{dummy} #call formfactwinM{dummy} #call offs{dummy} * .store * #call totf{'i'|'j'|'l'|tf} #call nsclrzng{'j'|'l'|dummy} *#call sclrzng{'j'|'l'|dummy} * #call tfml{tf|tfml} * print; * #endprocedure * #procedure genesisM(i,j,l) * #call vertf{'i'|'j'|'l'} * #call reduction{dummy} #call sing{dummy} #call formfactwin{dummy} #call offs{dummy} * .store * #call totf{'i'|'j'|'l'|tf} #call nsclrzng{'j'|'l'|dummy} *#call sclrzng{'j'|'l'|dummy} * #call tfms{tf|tfms} * nwrite statistics; print; * #endprocedure * #procedure singpart(dummy) * id c36(m?,mp?,ms?)= -1/12*dlt; al c35(m?,mp?,ms?)= -1/6*dlt; al c34(m?,mp?,ms?)= 0; al c33(m?,mp?,ms?)= 0; al c32(m?,mp?,ms?)= 0; al c31(m?,mp?,ms?)= 0; al c24(m?,mp?,ms?)= 1/4*dlt; al c23(m?,mp?,ms?)= 0; al c22(m?,mp?,ms?)= 0; al c21(m?,mp?,ms?)= 0; al c12(m?,mp?,ms?)= 0; al c11(m?,mp?,ms?)= 0; al c0(m?,mp?,ms?)= 0; al b1(ps?,m?,mp?)= -1/2*dlt; al b0(ps?,m?,mp?)= dlt; al a0(m?)= -m^2*dlt; * #endprocedure * #procedure coefficientm(dummy) * #do k1=1,10 g Wfm'k1'= tfml*sel'k1'/stw; #enddo nwrite statistics; * #do k1=1,10 id w('k1')*sel'k1'= 1; #enddo * #do k1=1,10 id sel'k1'= 0; #enddo * argument; id zero= 0; id um= tm; endargument; id zero= 0; * id zm= wm/ctw; id wm= rwm*ss; al wm^-1= 1/rwm/ss; id tm= rtm*ss; id wm= rwm*ss; al hm= rhm*ss; * id g^3= 1; al p12= 1/2*(-s-p11-p22); id gram(3,4)= -4/gramd/(y*S)^2; * id rtm^2= xt*rwm^2; id rhm^2= xh*rwm^2; * * id ctw^2= 1-stw^2; *id stw^2= 1-ctw^2; * .store * * g Cf1= Wfm1/ss; g Cf2= Wfm2/ss; g Cf3= Wfm3/ss; g Cf4= Wfm4/ss; g Cf5= Wfm5/ss; g Cf6= Wfm6/ss; g Cf7= Wfm7/ss; g Cf8= Wfm8/ss; g Cf9= Wfm9/ss; g Cf10= Wfm10/ss; * id p11= -s1; al p22= -s2; id s= -X*y*S; al s1= Y*y*S; al s2= x2*y*S; * b stw,ctw,x2,X,Y,y,S,c0,b0,a0; print; * #endprocedure * #procedure coefficientmoff(dummy) * #do k1=1,24 g Wfm'k1'= tfml*sel'k1'/stw; #enddo nwrite statistics; * #do k1=1,24 id w('k1')*sel'k1'= 1; #enddo * #do k1=1,24 id sel'k1'= 0; #enddo * argument; id zero= 0; id um= tm; endargument; id zero= 0; * id zm= wm/ctw; id wm= rwm*ss; al wm^-1= 1/rwm/ss; id tm= rtm*ss; id wm= rwm*ss; al hm= rhm*ss; * id g^3= 1; al p12= 1/2*(-s-p11-p22); id gram(3,4)= -4/gramd/(y*S)^2; * id rtm^2= xt*rwm^2; id rhm^2= xh*rwm^2; * * id ctw^2= 1-stw^2; *id stw^2= 1-ctw^2; * .store * g Cf1= Wfm1/ss; g Cf2= Wfm2/ss; g Cf3= Wfm3/ss; g Cf4= Wfm4/ss; g Cf5= Wfm5/ss; g Cf6= Wfm6/ss; g Cf7= Wfm7/ss; g Cf8= Wfm8/ss; g Cf9= Wfm9/ss; g Cf10= Wfm10/ss; g Cf11= Wfm11/ss; g Cf12= Wfm12/ss; g Cf13= Wfm13/ss; g Cf14= Wfm14/ss; g Cf15= Wfm15/ss; g Cf16= Wfm16/ss; g Cf17= Wfm17/ss; g Cf18= Wfm18/ss; g Cf19= Wfm19/ss; g Cf20= Wfm20/ss; g Cf21= Wfm21/ss; g Cf22= Wfm22/ss; g Cf23= Wfm23/ss; g Cf24= Wfm24/ss; * id p11= -s1; al p22= -s2; id s= -X*y*S; al s1= Y*y*S; al s2= x2*y*S; * b stw,ctw,x2,X,Y,y,S,c0,b0,a0; print; * #endprocedure * #procedure coefficientmoffs(dummy) * #do k1=1,6 g Wfm'k1'= tfms*sel'k1'/stw; #enddo nwrite statistics; * #do k1=1,6 id ws('k1')*sel'k1'= 1; id wss('k1')*sel'k1'= 1; #enddo * #do k1=1,6 id sel'k1'= 0; #enddo * argument; id zero= 0; id um= tm; endargument; id zero= 0; * id zm= wm/ctw; *id wm= rwm*ss; *al wm^-1= 1/rwm/ss; *id tm= rtm*ss; *id wm= rwm*ss; al hm= rhm*ss; id rtm= tm/s; * id g^3= 1; al p12= 1/2*(-s-p11-p22); id gram(3,4)= -4/gramd/(y*S)^2; * *id rtm^2= xt*rwm^2; id rhm^2= xh*rwm^2; * * id ctw^2= 1-stw^2; *id stw^2= 1-ctw^2; * .store * g Cf1= Wfm1/ss; g Cf2= Wfm2/ss; g Cf3= Wfm3/ss; g Cf4= Wfm4/ss; g Cf5= Wfm5/ss; g Cf6= Wfm6/ss; * id p11= -s1; al p22= -s2; id s= -X*y*S; al s^-1= -1/(X*y*S); al s1= Y*y*S; al s2= x2*y*S; id ss^-2= -1/(X*y*S); * b stw,ctw,x2,X,Y,y,S,c0,b0,a0; print; * #endprocedure * #procedure coefficientM(dummy) * #do k1=1,10 g WfM'k1'= tfms*sel'k1'/stw; #enddo nwrite statistics; * #do k1=1,10 id w('k1')*sel'k1'= 1; #enddo * #do k1=1,10 id sel'k1'= 0; #enddo * argument; id zero= 0; id um= tm; endargument; id zero= 0; * id tm= rtm*ss; * id g^3= 1; al p12= 1/2*(-s-p11-p22); id gram(3,4)= -4/gramd/(y*S)^2; * * id ctw^2= 1-stw^2; *id stw^2= 1-ctw^2; * .store * *g Cf1= WfM1*(1-selM)/ss; *g Cf2= WfM2*(1-selM)/ss; *g Cf3= WfM3*(1-selM)/ss; *g Cf4= WfM4*(1-selM)/ss; *g Cf5= WfM5*(1-selM)/ss; *g Cf6= WfM6*(1-selM)/ss; *g Cf7= WfM7*(1-selM)/ss; *g Cf8= WfM8*(1-selM)/ss; *g Cf9= WfM9*(1-selM)/ss; *g Cf10= WfM10*(1-selM)/ss; * g Cf1= WfM1/ss; g Cf2= WfM2/ss; g Cf3= WfM3/ss; g Cf4= WfM4/ss; g Cf5= WfM5/ss; g Cf6= WfM6/ss; g Cf7= WfM7/ss; g Cf8= WfM8/ss; g Cf9= WfM9/ss; g Cf10= WfM10/ss; * id p11= -s1; al p22= -s2; id s= -X*y*S; al s1= Y*y*S; al s2= x2*y*S; * id rtm^2= -tmr^2/(X*y); id ss^-2= -1/(X*y*S); id p11^-1= -1/(Y*y*S); al p22^-1= -1/(x2*y*S); * id tmr^4*selM= 0; id tmr^2*selM= 0; id a0(tm)*selM= 0; id b0(ps?,0,tm)*selM= b0(ps,0,0); id b0(ps?,tm,tm)*selM= b0(ps,0,0); id c0(tm,0,tm)*selM= c0(0,0,0); id c0(0,tm,0)*selM= c0(0,0,0); id selM= 1; * *id b0(ps?,0,tm)*tmr^4= aux4(ps); *id b0(ps?,0,tm)*tmr^2= aux2(ps); *id b0(ps?,0,tm)= b0(ps,0,0)+db0(ps,0,tm); *id aux4(ps?)= b0(ps,0,tm)*tmr^4; *al aux2(ps?)= b0(ps,0,tm)*tmr^2; * *id b0(ps?,tm,tm)*tmr^4= aux4(ps); *id b0(ps?,tm,tm)*tmr^2= aux2(ps); *id b0(ps?,tm,tm)= b0(ps,0,0)+db0(ps,tm,tm); *id aux4(ps?)= b0(ps,tm,tm)*tmr^4; *al aux2(ps?)= b0(ps,tm,tm)*tmr^2; * *id c21(tm,0,tm)= 1/2*(cp_21+cm_21); *al c21(0,tm,0)= cp_21-cm_21; *id c22(tm,0,tm)= 1/2*(cp_22+cm_22); *al c22(0,tm,0)= cp_22-cm_22; *id c23(tm,0,tm)= 1/2*(cp_23+cm_23); *al c23(0,tm,0)= cp_23-cm_23; *id c24(tm,0,tm)= 1/2*(cp_24+cm_24); *al c24(0,tm,0)= cp_24-cm_24; *id c31(tm,0,tm)= 1/2*(cp_31+cm_31); *al c31(0,tm,0)= cp_31-cm_31; *id c32(tm,0,tm)= 1/2*(cp_32+cm_32); *al c32(0,tm,0)= cp_32-cm_32; *id c33(tm,0,tm)= 1/2*(cp_33+cm_33); *al c33(0,tm,0)= cp_33-cm_33; *id c34(tm,0,tm)= 1/2*(cp_34+cm_34); *al c34(0,tm,0)= cp_34-cm_34; *id c35(tm,0,tm)= 1/2*(cp_35+cm_35); *al c35(0,tm,0)= cp_35-cm_35; *id c36(tm,0,tm)= 1/2*(cp_36+cm_36); *al c36(0,tm,0)= cp_36-cm_36; * b stw,ctw,X,Y,x2,y,S,gramd,tmr,c0,b0,db0,a0; print; * #endprocedure * #procedure coefficientMoff(dummy) * #do k1=1,24 g WfM'k1'= tfms*sel'k1'/stw; #enddo nwrite statistics; * #do k1=1,24 id w('k1')*sel'k1'= 1; #enddo * #do k1=1,24 id sel'k1'= 0; #enddo * argument; id zero= 0; id um= tm; endargument; id zero= 0; * id tm= rtm*ss; * id g^3= 1; al p12= 1/2*(-s-p11-p22); id gram(3,4)= -4/gramd/(y*S)^2; * id ctw^2= 1-stw^2; *id stw^2= 1-ctw^2; * .store * g Cf1= WfM1/ss; g Cf2= WfM2/ss; g Cf3= WfM3/ss; g Cf4= WfM4/ss; g Cf5= WfM5/ss; g Cf6= WfM6/ss; g Cf7= WfM7/ss; g Cf8= WfM8/ss; g Cf9= WfM9/ss; g Cf10= WfM10/ss; g Cf11= WfM11/ss; g Cf12= WfM12/ss; g Cf13= WfM13/ss; g Cf14= WfM14/ss; g Cf15= WfM15/ss; g Cf16= WfM16/ss; g Cf17= WfM17/ss; g Cf18= WfM18/ss; g Cf19= WfM19/ss; g Cf20= WfM20/ss; g Cf21= WfM21/ss; g Cf22= WfM22/ss; g Cf23= WfM23/ss; g Cf24= WfM24/ss; * id p11= -s1; al p22= -s2; id s= -X*y*S; al s1= Y*y*S; al s2= x2*y*S; * id rtm^2= -tmr^2/(X*y); id ss^-2= -1/(X*y*S); id p11^-1= -1/(Y*y*S); al p22^-1= -1/(x2*y*S); * b stw,ctw,X,Y,x2,y,S,gramd,tmr,c0,b0,db0,a0; print; * #endprocedure * #procedure flm(dummy) * g Wid= tfml/g^3/stw; * id w(1)= 4*p(mu)*(d_(mu,be)*p(al)+d_(be,al)*p1(mu)+d_(mu,al)*p2(be))/ss; id w(2)= 2*p(mu)*(d_(mu,be)*p(al)+d_(mu,al)*p2(be))/ss; id w(3)= 2*p(mu)*(d_(mu,be)*p(al)-d_(mu,al)*p2(be))/ss; id w(4)= -2*p(al)*p1.p*p2(be)/(X*y*S)/ss; id w(5)= 0; id w(6)= 2*e_(p1,p2,al,be)/ss; id w(7)= -e_(p1,p2,al,be)*p2.p/(X*y*S)/ss; id w(8)= -2*p(al)*ps*p2(be)/(X*y*S)/ss; id w(9)= 2*d_(al,be)*ps/ss; id w(10)= -e_(p1,p2,al,be)*p1.p/(X*y*S)/ss; * id e_(p1,p2,al?,be?)= -e_(Qp,Qm,al,be); id p.p= -s; al p.p1= -p11-p12; al p.p2= -p22-p12; id ps= -s; id p1(mu?)= Qp(mu); al p2(mu?)= -Qp(mu)-Qm(mu); al p(mu?)= Qm(mu); id p1(mu?)= Qp(mu); al p2(mu?)= -Qp(mu)-Qm(mu); al p(mu?)= Qm(mu); * id Qp(be)= 0; al Qp(al)= -Qm(al); id Qp.Qm= pp.Qm-qp.Qm; al Qm.Qm= QmS; * id pp.Qm= -1/2*(QmS+y*S); al Qm.qp= 1/2*(QmS+km); * id km= z*y*S; id QmS= X*y*S; * id z= -Y-1-X+x2; * id ps= -s; argument; id zero= 0; endargument; id zero= 0; * id p12= 1/2*(-s-p11-p22); id gram(3,4)= -4/gramd/(y*S)^2; id p11= -s1; al p22= -s2; id s= -X*y*S; al s1= Y*y*S; al s2= x2*y*S; * .store * g NWid= Wid*gramd^3; * id gramd= Y^2+X^2+x2^2-2*(-Y*X+Y*x2-X*x2); * repeat; id x2^2= gramd-X^2-Y^2+2*(-Y*X+Y*x2-X*x2); endrepeat; * b y,S,ss,c0,b0,a0,gramd,Qp,Qm,d_,e_; print; * #endprocedure * #procedure flmM(dummy) * g Wid= tfml/g^3/stw; * id w(1)= 4*p(mu)*(d_(mu,be)*p(al)+d_(be,al)*p1(mu)+d_(mu,al)*p2(be))/ss; id w(2)= 2*p(mu)*(d_(mu,be)*p(al)+d_(mu,al)*p2(be))/ss; id w(3)= 2*p(mu)*(d_(mu,be)*p(al)-d_(mu,al)*p2(be))/ss; id w(4)= -2*p(al)*p1.p*p2(be)/(X*y*S)/ss; id w(5)= 0; id w(6)= 2*e_(p1,p2,al,be)/ss; id w(7)= -e_(p1,p2,al,be)*p2.p/(X*y*S)/ss; id w(8)= -2*p(al)*ps*p2(be)/(X*y*S)/ss; id w(9)= 2*d_(al,be)*ps/ss; id w(10)= -e_(p1,p2,al,be)*p1.p/(X*y*S)/ss; id w(11)= -2/ss/(X*y*S)*p1.p*p1(al)*p1(be); id w(12)= -2/ss/(X*y*S)*p1.p*p1(al)*p2(be); id w(13)= -2/ss/(X*y*S)*p1.p*p1(be)*p2(al); id w(14)= -2/ss/(X*y*S)*p1.p*p2(al)*p2(be); id w(15)= -2/ss/(X*y*S)*p1(al)*p1(be)*p2.p; id w(16)= -2/ss/(X*y*S)*p1(al)*p2.p*p2(be); id w(17)= -2/ss/(X*y*S)*p1(be)*p2.p*p2(al); id w(18)= -2/ss/(X*y*S)*p2.p*p2(al)*p2(be); id w(19)= 2/ss*d_(p,al)*p1(be); id w(20)= 2/ss*d_(p,be)*p1(al); id w(21)= 2/ss*d_(p,be)*p2(al); id w(22)= 2/ss*d_(al,be)*p2.p; id w(23)= 0; * id e_(p1,p2,al?,be?)= -e_(Qp,Qm,al,be); id p.p= -s; al p.p1= -p11-p12; al p.p2= -p22-p12; id ps= -s; id p1(mu?)= Qp(mu); al p2(mu?)= -Qp(mu)-Qm(mu); al p(mu?)= Qm(mu); * id Qp.Qm= pp.Qm-qp.Qm; al Qm.Qm= QmS; * id pp.Qm= -1/2*(QmS+y*S); al Qm.qp= 1/2*(QmS+km); * id km= z*y*S; id QmS= X*y*S; * id z= -Y-1-X+x2; * id ps= -s; argument; id zero= 0; endargument; id zero= 0; * id p12= 1/2*(-s-p11-p22); id gram(3,4)= -4/gramd/(y*S)^2; id p11= -s1; al p22= -s2; id s= -X*y*S; al s1= Y*y*S; al s2= x2*y*S; * .store * g NWid= Wid*gramd^3; * id gramd= Y^2+X^2+x2^2-2*(-Y*X+Y*x2-X*x2); * repeat; id x2^2= gramd-X^2-Y^2+2*(-Y*X+Y*x2-X*x2); endrepeat; * b y,S,ss,c0,b0,a0,gramd,Qp,Qm,d_,e_; print; * #endprocedure * #procedure flMM(dummy) * g Wid= tfms/g^3/stw; nwrite statistics; * id w(1)= 4*p(mu)*(d_(mu,be)*p(al)+d_(be,al)*p1(mu)+d_(mu,al)*p2(be))/ss; id w(2)= 2*p(mu)*(d_(mu,be)*p(al)+d_(mu,al)*p2(be))/ss; id w(3)= 2*p(mu)*(d_(mu,be)*p(al)-d_(mu,al)*p2(be))/ss; id w(4)= -2*p(al)*p1.p*p2(be)/(X*y*S)/ss; id w(5)= 0; id w(6)= 2*e_(p1,p2,al,be)/ss; id w(7)= -e_(p1,p2,al,be)*p2.p/(X*y*S)/ss; id w(8)= -2*p(al)*ps*p2(be)/(X*y*S)/ss; id w(9)= 2*d_(al,be)*ps/ss; id w(10)= -e_(p1,p2,al,be)*p1.p/(X*y*S)/ss; id w(11)= -2/ss/(X*y*S)*p1.p*p1(al)*p1(be); id w(12)= -2/ss/(X*y*S)*p1.p*p1(al)*p2(be); id w(13)= -2/ss/(X*y*S)*p1.p*p1(be)*p2(al); id w(14)= -2/ss/(X*y*S)*p1.p*p2(al)*p2(be); id w(15)= -2/ss/(X*y*S)*p1(al)*p1(be)*p2.p; id w(16)= -2/ss/(X*y*S)*p1(al)*p2.p*p2(be); id w(17)= -2/ss/(X*y*S)*p1(be)*p2.p*p2(al); id w(18)= -2/ss/(X*y*S)*p2.p*p2(al)*p2(be); id w(19)= 2/ss*d_(p,al)*p1(be); id w(20)= 2/ss*d_(p,be)*p1(al); id w(21)= 2/ss*d_(p,be)*p2(al); id w(22)= 2/ss*d_(al,be)*p2.p; id w(23)= 0; id w(24)= 0; * id ws(1)= p(be); id ws(2)= -p1.p*p2(be)/(X*y*S); id ws(3)= 0; id ws(4)= -p1.p*p1(be)/(X*y*S); id ws(5)= -p2.p*p1(be)/(X*y*S); id ws(6)= -p2.p*p2(be)/(X*y*S); * id wss(1)= p(al); id wss(2)= -p1.p*p2(al)/(X*y*S); id wss(3)= 0; id wss(4)= -p1.p*p1(al)/(X*y*S); id wss(5)= -p2.p*p1(al)/(X*y*S); id wss(6)= -p2.p*p2(al)/(X*y*S); * id wb(1)= p1.p/ss; id wb(2)= p2.p/ss; * id e_(p1,p2,al?,be?)= -e_(Qp,Qm,al,be); id p.p= -s; al p.p1= -p11-p12; al p.p2= -p22-p12; id ps= -s; id p1(mu?)= Qp(mu); al p2(mu?)= -Qp(mu)-Qm(mu); al p(mu?)= Qm(mu); * id Qp.Qm= pp.Qm-qp.Qm; al Qm.Qm= QmS; * id pp.Qm= -1/2*(QmS+y*S); al Qm.qp= 1/2*(QmS+km); * id km= z*y*S; id QmS= X*y*S; * id z= -Y-1-X+x2; * id ps= -s; argument; id zero= 0; endargument; id zero= 0; * id p12= 1/2*(-s-p11-p22); id gram(3,4)= -4/gramd/(y*S)^2; id p11= -s1; al p22= -s2; id s= -X*y*S; al s1= Y*y*S; al s2= x2*y*S; * .store * g NWid= Wid*gramd^3; nwrite statistics; * id gramd= Y^2+X^2+x2^2-2*(-Y*X+Y*x2-X*x2); * id p11^-1= -1/(Y*y*S); al p22^-1= -1/(x2*y*S); id rtm^2= -tmr^2/(X*y); id ss^-2= -1/(X*y*S); al ss^2= -(X*y*S); * repeat; id x2^2= gramd-X^2-Y^2+2*(-Y*X+Y*x2-X*x2); endrepeat; * id gramd^3*x2= aux; id x2= (gramd-X^2-Y^2+2*(-Y*X+Y*x2-X*x2))/x2; id aux= x2*gramd^3; * .store * g NNWid= NWid; nwrite statistics; * id gramd^3= aux; id gramd= Y^2+X^2+x2^2-2*(-Y*X+Y*x2-X*x2); * repeat; id Y^2= gramd-X^2-x2^2+2*(-Y*X+Y*x2-X*x2); endrepeat; * id aux= gramd^3; * b i_,y,S,ss,c0,b0,a0,gramd,Qp,Qm,d_,e_,tmr,tm,wm; print; * #endprocedure * #procedure flM(dummy) * g Wid= tfms/g^3/stw; * id w(1)= 4*p(mu)*(d_(mu,be)*p(al)+d_(be,al)*p1(mu)+d_(mu,al)*p2(be))/ss; id w(2)= 2*p(mu)*(d_(mu,be)*p(al)+d_(mu,al)*p2(be))/ss; id w(3)= 2*p(mu)*(d_(mu,be)*p(al)-d_(mu,al)*p2(be))/ss; id w(4)= -2*p(al)*p1.p*p2(be)/(X*y*S)/ss; id w(5)= 0; id w(6)= 2*e_(p1,p2,al,be)/ss; id w(7)= -e_(p1,p2,al,be)*p2.p/(X*y*S)/ss; id w(8)= -2*p(al)*ps*p2(be)/(X*y*S)/ss; id w(9)= 2*d_(al,be)*ps/ss; id w(10)= -e_(p1,p2,al,be)*p1.p/(X*y*S)/ss; * id e_(p1,p2,al?,be?)= -e_(Qp,Qm,al,be); id p.p= -s; al p.p1= -p11-p12; al p.p2= -p22-p12; id ps= -s; id p1(mu?)= Qp(mu); al p2(mu?)= -Qp(mu)-Qm(mu); al p(mu?)= Qm(mu); * id Qp(be)= 0; al Qp(al)= -Qm(al); id Qp.Qm= pp.Qm-qp.Qm; al Qm.Qm= QmS; * id pp.Qm= -1/2*(QmS+y*S); al Qm.qp= 1/2*(QmS+km); * id km= z*y*S; id QmS= X*y*S; * id z= -Y-1-X+x2; * id ps= -s; argument; id zero= 0; endargument; id zero= 0; * id p12= 1/2*(-s-p11-p22); id gram(3,4)= -4/gramd/(y*S)^2; id p11= -s1; al p22= -s2; id s= -X*y*S; al s1= Y*y*S; al s2= x2*y*S; * .store * g NWid= Wid*gramd^3; * id gramd= Y^2+X^2+x2^2-2*(-Y*X+Y*x2-X*x2); * id p11^-1= -1/(Y*y*S); al p22^-1= -1/(x2*y*S); id rtm^2= -tmr^2/(X*y); id ss^-2= -1/(X*y*S); * repeat; id x2^2= gramd-X^2-Y^2+2*(-Y*X+Y*x2-X*x2); endrepeat; * id gramd^3*x2= aux; id x2= (gramd-X^2-Y^2+2*(-Y*X+Y*x2-X*x2))/x2; id aux= x2*gramd^3; * b y,S,ss,c0,b0,a0,gramd,Qp,Qm,d_,e_,tmr; print; * #endprocedure * */////////////////////////////////////////////////////////////// * *#call genesism{1|3|4} *#call genesisM{1|3|4} *#call genesismM{1|3|4} *#call genesismM{1|3|8} *#call genesisMM{1|3|4} *#call genesisMM{1|3|8} #call genesisMM{1|7|4} *#call genesisMM{1|7|8} *#call coefficientm{dummy} *#call coefficientM{dummy} *#call coefficientMoff{dummy} *#call coefficientmoff{dummy} #call coefficientmoffs{dummy} *#call flm{dummy} *#call flmM{dummy} *#call flMM{dummy} *#call flM{dummy} #call singpart{dummy} * */////////////////////////////////////////////////////////////// * .end