pyapps.f

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file pyapps.f

 
C*********************************************************************
 
C...PYAPPS
C...Uses approximate analytical formulae to determine the full set of
C...MSSM parameters from SUGRA input.
C...See M. Drees and S.P. Martin, hep-ph/9504124
 
      SUBROUTINE pyapps
 
C...Double precision and integer declarations.
      IMPLICIT DOUBLE PRECISION(a-h, o-z)
      IMPLICIT INTEGER(i-n)
      INTEGER pyk,pychge,pycomp
C...Parameter statement to help give large particle numbers.
      parameter(ksusy1=1000000,ksusy2=2000000,ktechn=3000000,
     &kexcit=4000000,kdimen=5000000)
C...Commonblocks.
      common/pydat1/mstu(200),paru(200),mstj(200),parj(200)
      common/pydat2/kchg(500,4),pmas(500,4),parf(2000),vckm(4,4)
      common/pymssm/imss(0:99),rmss(0:99)
      SAVE /pydat1/,/pydat2/,/pymssm/

      WRITE(mstu(11),*) '(PYAPPS:) approximate mSUGRA relations'//
     &' not intended for serious physics studies'
      imss(5)=0
      imss(8)=0
      xmt=pmas(6,1)
      xmz2=pmas(23,1)**2
      xmw2=pmas(24,1)**2
      tanb=rmss(5)
      beta=atan(tanb)
      xw=paru(102)
      xmg=rmss(1)
      xmg2=xmg*xmg
      xm0=rmss(8)
      xm02=xm0*xm0
C...Temporary sign change for AT. Others unchanged.
      at=-rmss(16)
      rmss(15)=rmss(16)
      rmss(17)=rmss(16)
      sinb=tanb/sqrt(tanb**2+1d0)
      cosb=sinb/tanb
 
      dterm=xmz2*cos(2d0*beta)
      xmer=sqrt(xm02+0.15d0*xmg2-xw*dterm)
      xmel=sqrt(xm02+0.52d0*xmg2-(0.5d0-xw)*dterm)
      rmss(6)=xmel
      rmss(7)=xmer
      xmur=sqrt(pyrnmq(2,2d0/3d0*xw*dterm))
      xmdr=sqrt(pyrnmq(3,-1d0/3d0*xw*dterm))
      xmul=sqrt(pyrnmq(1,(0.5d0-2d0/3d0*xw)*dterm))
      xmdl=sqrt(pyrnmq(1,-(0.5d0-1d0/3d0*xw)*dterm))
      DO 100 i=1,5,2
        pmas(pycomp(ksusy1+i),1)=xmdl
        pmas(pycomp(ksusy2+i),1)=xmdr
        pmas(pycomp(ksusy1+i+1),1)=xmul
        pmas(pycomp(ksusy2+i+1),1)=xmur
  100 CONTINUE
      xarg=xmel**2-xmw2*abs(cos(2d0*beta))
      IF(xarg.LT.0d0) THEN
        WRITE(mstu(11),*) ' SNEUTRINO MASS IS NEGATIVE'//
     &  ' FROM THE SUM RULE. '
        WRITE(mstu(11),*) '  TRY A SMALLER VALUE OF TAN(BETA). '
        RETURN
      ELSE
        xarg=sqrt(xarg)
      ENDIF
      DO 110 i=11,15,2
        pmas(pycomp(ksusy1+i),1)=xmel
        pmas(pycomp(ksusy2+i),1)=xmer
        pmas(pycomp(ksusy1+i+1),1)=xarg
        pmas(pycomp(ksusy2+i+1),1)=9999d0
  110 CONTINUE
      rmt=pymrun(6,pmas(6,1)**2)
      xtop=(rmt/150d0/sinb)**2*(.9d0*xm02+2.1d0*xmg2+
     &(1d0-(rmt/190d0/sinb)**3)*(.24d0*at**2+at*xmg))
      rmb=pymrun(5,pmas(6,1)**2)
      xbot=(rmb/150d0/cosb)**2*(.9d0*xm02+2.1d0*xmg2+
     &(1d0-(rmb/190d0/cosb)**3)*(.24d0*at**2+at*xmg))
      xtau=1d-4/cosb**2*(xm02+0.15d0*xmg2+at**2/3d0)
      atp=at*(1d0-(rmt/190d0/sinb)**2)+xmg*(3.47d0-1.9d0*(rmt/190d0/
     &sinb)**2)
      rmss(16)=-atp
      xmu2=-.5d0*xmz2+(sinb**2*(xm02+.52d0*xmg2-xtop)-
     &cosb**2*(xm02+.52d0*xmg2-xbot-xtau/3d0))/(cosb**2-sinb**2)
      xma2=2d0*(xm02+.52d0*xmg2+xmu2)-xtop-xbot-xtau/3d0
      xmu=sign(sqrt(xmu2),rmss(4))
      rmss(4)=xmu
      IF(xma2.GT.0d0) THEN
        rmss(19)=sqrt(xma2)
      ELSE
        WRITE(mstu(11),*) ' PYAPPS:: PSEUDOSCALAR MASS**2 < 0 '
        CALL pystop(102)
      ENDIF
      arg=xm02+0.15d0*xmg2-2d0*xtau/3d0-xw*dterm
      IF(arg.GT.0d0) THEN
        rmss(14)=sqrt(arg)
      ELSE
        WRITE(mstu(11),*) ' PYAPPS:: RIGHT STAU MASS**2 < 0 '
        CALL pystop(102)
      ENDIF
      arg=xm02+0.52d0*xmg2-xtau/3d0-(0.5d0-xw)*dterm
      IF(arg.GT.0d0) THEN
        rmss(13)=sqrt(arg)
      ELSE
        WRITE(mstu(11),*) ' PYAPPS::  LEFT STAU MASS**2 < 0 '
        CALL pystop(102)
      ENDIF
      arg=pyrnmq(1,-(xbot+xtop)/3d0)
      IF(arg.GT.0d0) THEN
        rmss(10)=sqrt(arg)
      ELSE
        rmss(10)=-sqrt(-arg)
      ENDIF
      arg=pyrnmq(2,-2d0*xtop/3d0)
      IF(arg.GT.0d0) THEN
        rmss(12)=sqrt(arg)
      ELSE
        rmss(12)=-sqrt(-arg)
      ENDIF
      arg=pyrnmq(3,-2d0*xbot/3d0)
      IF(arg.GT.0d0) THEN
        rmss(11)=sqrt(arg)
      ELSE
        rmss(11)=-sqrt(-arg)
      ENDIF
 
      RETURN
      END