input_case nil$ in "kpack.red","pol.red"$ % Nonleptonic hyperon decays % -------------------------- mass p=1,p1=m1,p2=m2; mshell p,p1,p2; let p.p1=(1+m1^2-m2^2)/2,p.p2=(1-m1^2+m2^2)/2, p1.p2=(1-m1^2-m2^2)/2; MM:=4*(1+b*g(f,a))*(g(f,p)+1)*(1-bc*g(f,a))*(g(f,p1)+m1)$ on gcd; Gamma:=MM/2*ph2(1,m1,m2); procedure NonRel(x); begin scalar u,v,w; let abs(d)=d; w:=sub(m1=1-d,m2=d*m2,b=2*b/d/sqrt(1-m2^2), bc=2*bc/d/sqrt(1-m2^2),x); u:=num(w); v:=den(w); while (w:=sub(d=0,u))=0 do u:=u/d; u:=w; while (w:=sub(d=0,v))=0 do v:=v/d; return u/w end$ NonRel(Gamma); % polarized hyperon vector ee; let p.ee=0,p1.ee=-p2.ee,p2.ee=-PP(1,m1,m2)*c; MM:=(1+b*g(f,a))*(g(f,p)+1)*(1-g(f,a,ee)) *(1-bc*g(f,a))*(g(f,p1)+m1)$ df(MM,c)/sub(c=0,MM); NonRel(ws); % final baryon polarization let abs(m1)=m1; l1:=-second(ort({p1.l,p.l}))$ index l; pe:=p.l*l1$ remind l; let p1.ee=0,p2.ee=p.ee,p.ee=z*pe; MM:=(1+b*g(f,a))*(g(f,p)+1) *(1-bc*g(f,a))*(g(f,p1)+m1)*(1-g(f,a,ee))$ df(MM,z)/sub(z=0,MM); NonRel(ws); end;