c======================================================================= c If you copy this code, please send a message to c RIDOLFI@VXCERN.CERN.CH c or to c MLM@VXCERN.CERN.CH c so that we can keep you informed of upgrades or modifications. c======================================================================= implicit real*8 (a-h,o-z) parameter(sin2w=.232d0,xmz=91.188d0,gz=2.5d0,xmw=80.2) dimension y(4,4),zr(4,4),work(4),axm(4) real*4 s_y(4,4),s_xm(4),s_zr(4,4),s_work(4) common/prcpar/xmrsel,xmlsel,ine,jne common/beams/s,sthr common/susy/xn(4,4),xm(4),tb common/bveg1/xl(10),xu(10),acc,ndim,ncall,itmx,nprn external radsig c Initialize parameters for Vegas do j=1,10 xl(j) = 0 xu(j) = 1 enddo acc = -1 ndim = 1 ncall = 10000 itmx = 4 1 write(*,*)'Enter sqrt(s), M, mu, tan(beta), m(selR), m(selL),', # 'i, j (i <= j)' read (*,*)sr,rm,rmu,tb,xmrsel,xmlsel,ix,jx if(sr.lt.0)stop if(ix.gt.jx)then write(*,*)' i > j forbidden' goto 1 endif s=sr**2 cos2w=1-sin2w cw=sqrt(cos2w) sw=sqrt(sin2w) cb=1.d0/sqrt(1.d0+tb**2) sb=tb/sqrt(1.d0+tb**2) c2b=cb**2-sb**2 s2b=2*cb*sb c Y is the neutralino mass matrix as in Nucl. Phys. B278(1986)1 rmp=rm*5.d0/3.d0*sin2w/cos2w y(1,1)=rmp*cos2w+rm*sin2w y(1,2)=(rm-rmp)*sw*cw y(1,3)=0 y(1,4)=0 y(2,1)=y(1,2) y(2,2)=rmp*sin2w+rm*cos2w y(2,3)=xmz y(2,4)=0 y(3,1)=y(1,3) y(3,2)=y(2,3) y(3,3)=rmu*s2b y(3,4)=-rmu*c2b y(4,1)=y(1,4) y(4,2)=y(2,4) y(4,3)=y(3,4) y(4,4)=-rmu*s2b c Diagonalization of Y; the eigenvalues (WITH THEIR SIGN!) are c contained in the vector XM. c c Transfer first to sngl precision, as eisrs1 works with real*4 do i=1,4 do j=1,4 s_y(i,j)=sngl(y(i,j)) enddo enddo c call eisrs1 (4,4,s_y,s_xm,s_zr,ierr,s_work) c c Transfer back to real*8 do i=1,4 do j=1,4 zr(i,j)=dble(s_zr(i,j)) enddo xm(i)=dble(s_xm(i)) work(i)=dble(s_work(i)) enddo c do i=1,4 do j=1,4 xn(i,j)=zr(j,i) enddo enddo c Now the matrix XN is the same as N in Nucl. Phys. B278(1986)1 c c ordinamento dei neutralini [axm(l1)