ISAJET/isasusy/ssstbf.F

   1 #include "isajet/pilot.h"
   2         SUBROUTINE SSSTBF
   3 C-----------------------------------------------------------------------
   4 C
   5 C        This program gives stop squark branching fractions to gauginos
   6 C        according to Baer and Tata.
   7 C        If no other modes are allowed, stop -> c z_i through loops is
   8 C        used as the default.
   9 C
  10 C-----------------------------------------------------------------------
  11 #if defined(CERNLIB_IMPNONE)
  12       IMPLICIT NONE
  13 #endif
  14 #include "isajet/sslun.inc"
  15 #include "isajet/ssmode.inc"
  16 #include "isajet/sssm.inc"
  17 #include "isajet/sspar.inc"
  18 #include "isajet/sstype.inc"
  19 C
  20       COMPLEX ZI,ZONE,ZA,ZB,ZPP,ZPM,ZAUIZ,ZBUIZ
  21       DOUBLE PRECISION SSALFS,SSMQCD
  22       REAL SSXLAM
  23       REAL WID,AWD(2),BW(2),FB,FT,XM,YM,THX,THY,AU1,MZ1,WIDC1
  24       REAL PI,SR2,G,GP,TANB,COTB,MPL,MMI,AH
  25       REAL AUIZ,MZIZ,SINT,COST,AS,BS,SNZI,THIZ
  26       INTEGER IZ,ISTOP,IDSTOP
  27       REAL AMSTOP,BWP(2),A
  28       REAL MW1,MW2,SNW1,SNW2,CS2THW,BETA,TN2THW,SINB,COSB
  29       REAL ASMB,MBMB,MBQ,ASMT,MTMT,MTQ,SUALFS
  30       INTEGER ISZIZ(4)
  31       DATA ZONE/(1.,0.)/,ZI/(0.,1.)/
  32 C
  33 C          Partly duplicated from SSMASS
  34 C
  35       CS2THW=1.-SN2THW
  36       TN2THW=SN2THW/CS2THW
  37       PI=4.*ATAN(1.)
  38       SR2=SQRT(2.)
  39       G=SQRT(4*PI*ALFAEM/SN2THW)
  40       GP=G*SQRT(SN2THW/(1.-SN2THW))
  41       TANB=1./RV2V1
  42       COTB=RV2V1
  43       BETA=ATAN(TANB)
  44 C          Reconstruct masses from SSMASS
  45       ASMB=SUALFS(AMBT**2,.36,AMTP,3)
  46       MBMB=AMBT*(1.-4*ASMB/3./PI)
  47       MBQ=SSMQCD(DBLE(MBMB),DBLE(AMT1SS))
  48       ASMT=SUALFS(AMTP**2,.36,AMTP,3)
  49       MTMT=AMTP/(1.+4*ASMT/3./PI+(16.11-1.04*(5.-6.63/AMTP))*
  50      $(ASMT/PI)**2)
  51       MTQ=SSMQCD(DBLE(MTMT),DBLE(AMT1SS))
  52       FB=G*MBQ/SR2/AMW/COS(BETA)
  53       FT=G*MTQ/SR2/AMW/SIN(BETA)
  54       MW1=ABS(AMW1SS)
  55       MW2=ABS(AMW2SS)
  56       SNW1=SIGN(1.,AMW1SS)
  57       SNW2=SIGN(1.,AMW2SS)
  58       XM=1./TAN(GAMMAL)
  59       YM=1./TAN(GAMMAR)
  60       THX=SIGN(1.,XM)
  61       THY=SIGN(1.,YM)
  62 C
  63       AWD(1)=-G*SNW1*SIN(GAMMAR)
  64       AWD(2)=-G*SNW2*THY*COS(GAMMAR)
  65       BW(1)=-FT*SNW1*COS(GAMMAR)
  66       BW(2)=FT*SNW2*THY*SIN(GAMMAR)
  67       BWP(1)=-FB*COS(GAMMAL)
  68       BWP(2)=FB*THX*SIN(GAMMAL)
  69       MMI=AMW1SS
  70       MPL=AMW2SS
  71       COST=COS(THETAT)
  72       SINT=SIN(THETAT)
  73       COSB=COS(THETAB)
  74       SINB=SIN(THETAB)
  75 C
  76 C          Compute stop_i branching fractions to charm + zi if no other
  77 C          modes are allowed. WIDC1 is an unknown amplitude arbitrarily
  78 C          set equal to 1.0.
  79 C
  80       ISZIZ(1)=ISZ1
  81       ISZIZ(2)=ISZ2
  82       ISZIZ(3)=ISZ3
  83       ISZIZ(4)=ISZ4
  84       AU1=-G/SR2*ZMIXSS(3,1)-GP/3./SR2*ZMIXSS(4,1)
  85       MZ1=ABS(AMZISS(1))
  86       DO 100 ISTOP=1,2
  87         IF(ISTOP.EQ.1) THEN
  88           AMSTOP=AMT1SS
  89           IDSTOP=ISTP1
  90           WIDC1=1.0E-6
  91         ELSE
  92           AMSTOP=AMT2SS
  93           IDSTOP=ISTP2
  94           WIDC1=1.0E-6
  95         ENDIF
  96         IF(AMSTOP.LT.(MW1+AMBT).AND.AMSTOP.GT.(AMCH+MZ1)) THEN
  97           DO 110 IZ=1,4
  98             MZIZ=ABS(AMZISS(IZ))
  99             AUIZ=-G/SR2*ZMIXSS(3,IZ)-GP/3./SR2*ZMIXSS(4,IZ)
 100             IF (AMT1SS.GT.(AMCH+MZIZ)) THEN
 101               WID=AUIZ**2*(AMSTOP**2-MZIZ**2)/AU1**2
 102      $        /(AMSTOP**2-MZ1**2)*WIDC1
 103               CALL SSSAVE(IDSTOP,WID,ISZIZ(IZ),IDCH,0,0,0)
 104             END IF
 105 110       CONTINUE
 106         ELSEIF(AMSTOP.LT.(MW1+AMBT).AND.AMSTOP.LE.(AMCH+MZ1)) THEN
 107           WRITE(LOUT,1000) ISTOP
 108 1000      FORMAT(' ERROR IN SSSTBF: NO ALLOWED MODE FOR STOP',I2)
 109         END IF
 110 100   CONTINUE
 111 C
 112 C          stop_i -> gluino + top
 113 C
 114       IF (AMT1SS.GT.(AMGLSS+AMTP)) THEN
 115         WID=2*SSALFS(DBLE(AMT1SS**2))*AMT1SS*((1.-AMGLSS**2/AMT1SS**2-
 116      $  AMTP**2/AMT1SS**2)-2*SIN(2*THETAT)*AMTP*AMGLSS/AMT1SS**2)
 117      $  *SQRT(SSXLAM(1.,AMGLSS**2/AMT1SS**2,AMTP**2/AMT1SS**2))/3.
 118         CALL SSSAVE(ISTP1,WID,ISGL,IDTP,0,0,0)
 119       END IF
 120 C
 121       IF (AMT2SS.GT.(AMGLSS+AMTP)) THEN
 122         WID=2*SSALFS(DBLE(AMT2SS**2))*AMT2SS*((1.-AMGLSS**2/AMT2SS**2-
 123      $  AMTP**2/AMT2SS**2)+2*SIN(2*THETAT)*AMTP*AMGLSS/AMT2SS**2)
 124      $  *SQRT(SSXLAM(1.,AMGLSS**2/AMT2SS**2,AMTP**2/AMT2SS**2))/3.
 125         CALL SSSAVE(ISTP2,WID,ISGL,IDTP,0,0,0)
 126       END IF
 127 C
 128 C          stop_1 -> top + zino_i
 129 C
 130       DO 200 IZ=1,4
 131         MZIZ=ABS(AMZISS(IZ))
 132         SNZI=SIGN(1.,AMZISS(IZ))
 133         IF (SNZI.EQ.1.) THEN
 134            THIZ=0.
 135         ELSE
 136            THIZ=1.
 137         END IF
 138         ZAUIZ=ZI**(THIZ-1.)*SNZI
 139      $  *(-G/SR2*ZMIXSS(3,IZ)-GP/3./SR2*ZMIXSS(4,IZ))
 140         ZBUIZ=ZI**(THIZ-1.)*4*GP*ZMIXSS(4,IZ)/3./SR2
 141         ZPP=ZI**THIZ
 142         ZPM=(-ZI)**THIZ
 143         ZA=((ZI*ZAUIZ-ZPP*FT*ZMIXSS(1,IZ))*COST
 144      $  -(ZI*ZBUIZ-ZPM*FT*ZMIXSS(1,IZ))*SINT)/2.
 145         ZB=((-ZI*ZAUIZ-ZPP*FT*ZMIXSS(1,IZ))*COST
 146      $  -(ZI*ZBUIZ+ZPM*FT*ZMIXSS(1,IZ))*SINT)/2.
 147         AS=ZA*CONJG(ZA)
 148         BS=ZB*CONJG(ZB)
 149         IF (AMT1SS.GT.(AMTP+MZIZ)) THEN
 150           WID=(AS*(AMT1SS**2-(AMTP+MZIZ)**2)+BS
 151      $    *(AMT1SS**2-(AMTP-MZIZ)**2))/8./PI/AMT1SS
 152      $    *SQRT(SSXLAM(1.,AMTP**2/AMT1SS**2,MZIZ**2/AMT1SS**2))
 153           CALL SSSAVE(ISTP1,WID,ISZIZ(IZ),IDTP,0,0,0)
 154         END IF
 155 200   CONTINUE
 156 C
 157 C          Wino decays
 158 C
 159       IF (AMT1SS.GT.(AMBT+MW1)) THEN
 160         A=AWD(1)*COST-BW(1)*SINT
 161         AS=A*A
 162         WID=AMT1SS*((AS+BWP(1)**2*COST**2)*(1.-MW1**2/AMT1SS**2-
 163      $   AMBT**2/AMT1SS**2)-4*MW1*AMBT*BWP(1)*COST*A/AMT1SS**2)
 164      $   *SQRT(SSXLAM(1.,MW1**2/AMT1SS**2,AMBT**2/AMT1SS**2))/16./PI
 165         CALL SSSAVE(ISTP1,WID,ISW1,IDBT,0,0,0)
 166       END IF
 167       IF (AMT1SS.GT.(AMBT+MW2)) THEN
 168         A=AWD(2)*COST-BW(2)*SINT
 169         AS=A*A
 170         WID=AMT1SS*((AS+BWP(2)**2*COST**2)*(1.-MW2**2/AMT1SS**2-
 171      $   AMBT**2/AMT1SS**2)-4*MW2*AMBT*BWP(2)*COST*A/AMT1SS**2)
 172      $   *SQRT(SSXLAM(1.,MW2**2/AMT1SS**2,AMBT**2/AMT1SS**2))/16./PI
 173         CALL SSSAVE(ISTP1,WID,ISW2,IDBT,0,0,0)
 174       END IF
 175 C
 176       IF (AMT2SS.GT.(AMBT+MW1)) THEN
 177         A=AWD(1)*SINT+BW(1)*COST
 178         AS=A*A
 179         WID=AMT2SS*((AS+BWP(1)**2*SINT**2)*(1.-MW1**2/AMT2SS**2-
 180      $   AMBT**2/AMT2SS**2)-4*MW1*AMBT*BWP(1)*SINT*A/AMT2SS**2)
 181      $   *SQRT(SSXLAM(1.,MW1**2/AMT2SS**2,AMBT**2/AMT2SS**2))/16./PI
 182         CALL SSSAVE(ISTP2,WID,ISW1,IDBT,0,0,0)
 183       END IF
 184       IF (AMT2SS.GT.(AMBT+MW2)) THEN
 185         A=AWD(2)*SINT+BW(2)*COST
 186         AS=A*A
 187         WID=AMT2SS*((AS+BWP(2)**2*SINT**2)*(1.-MW2**2/AMT2SS**2-
 188      $   AMBT**2/AMT2SS**2)-4*MW2*AMBT*BWP(2)*SINT*A/AMT2SS**2)
 189      $  *SQRT(SSXLAM(1.,MW2**2/AMT2SS**2,AMBT**2/AMT2SS**2))/16./PI
 190         CALL SSSAVE(ISTP2,WID,ISW2,IDBT,0,0,0)
 191       END IF
 192 C
 193 C          stop_2 -> stop_1 + X modes
 194 C
 195       IF (AMT2SS.GT.(AMT1SS+AMZ)) THEN
 196         WID=G**2*COST**2*SINT**2
 197      $  *(SQRT(SSXLAM(AMT2SS**2,AMZ**2,AMT1SS**2)))**3
 198      $  /64./PI/CS2THW/AMT2SS**3/AMZ**2
 199         CALL SSSAVE(ISTP2,WID,IDZ,ISTP1,0,0,0)
 200       END IF
 201 C
 202       IF (AMT2SS.GT.(AMT1SS+AMHL)) THEN
 203         AH=G*AMW*SIN(BETA-ALFAH)*(1.-5.*TN2THW/3.)*SINT*COST/2.
 204      $  +G*AMTP*COS(2.*THETAT)*(TWOM1*SIN(ALFAH)+AAT*COS(ALFAH))/2.
 205      $  /AMW/SIN(BETA)
 206         WID=AH**2/16./PI/AMT2SS**3
 207      $  *SQRT(SSXLAM(AMT2SS**2,AMHL**2,AMT1SS**2))
 208         CALL SSSAVE(ISTP2,WID,ISHL,ISTP1,0,0,0)
 209       END IF
 210 C
 211       IF (AMT2SS.GT.(AMT1SS+AMHH)) THEN
 212         AH=-G*AMW*COS(BETA-ALFAH)*(1.-5.*TN2THW/3.)*SINT*COST/2.
 213      $  +G*AMTP*COS(2.*THETAT)*(TWOM1*COS(ALFAH)-AAT*SIN(ALFAH))/2.
 214      $  /AMW/SIN(BETA)
 215         WID=AH**2/16./PI/AMT2SS**3
 216      $  *SQRT(SSXLAM(AMT2SS**2,AMHH**2,AMT1SS**2))
 217         CALL SSSAVE(ISTP2,WID,ISHH,ISTP1,0,0,0)
 218       END IF
 219 C
 220       IF (AMT2SS.GT.(AMT1SS+AMHA)) THEN
 221         AH=G*AMTP*(TWOM1-AAT/TANB)/2./AMW
 222         WID=AH**2/16./PI/AMT2SS**3
 223      $  *SQRT(SSXLAM(AMT2SS**2,AMHA**2,AMT1SS**2))
 224         CALL SSSAVE(ISTP2,WID,ISHA,ISTP1,0,0,0)
 225       END IF
 226 C
 227 C          t_i --> b_i + W decays
 228 C
 229       IF (AMT1SS.GT.(AMB1SS+AMW)) THEN
 230         WID=G**2*COST**2*COSB**2*(SSXLAM(AMT1SS**2,AMB1SS**2,
 231      $AMW**2))**1.5/32./PI/AMT1SS**3/AMW**2
 232         CALL SSSAVE(ISTP1,WID,IDW,ISBT1,0,0,0)
 233       END IF
 234 C
 235       IF (AMT1SS.GT.(AMB2SS+AMW)) THEN
 236         WID=G**2*COST**2*SINB**2*(SSXLAM(AMT1SS**2,AMB2SS**2,
 237      $AMW**2))**1.5/32./PI/AMT1SS**3/AMW**2
 238         CALL SSSAVE(ISTP1,WID,IDW,ISBT2,0,0,0)
 239       END IF
 240 C
 241       IF (AMT2SS.GT.(AMB1SS+AMW)) THEN
 242         WID=G**2*SINT**2*COSB**2*(SSXLAM(AMT2SS**2,AMB1SS**2,
 243      $AMW**2))**1.5/32./PI/AMT2SS**3/AMW**2
 244         CALL SSSAVE(ISTP2,WID,IDW,ISBT1,0,0,0)
 245       END IF
 246 C
 247       IF (AMT2SS.GT.(AMB2SS+AMW)) THEN
 248         WID=G**2*SINT**2*SINB**2*(SSXLAM(AMT2SS**2,AMB2SS**2,
 249      $AMW**2))**1.5/32./PI/AMT2SS**3/AMW**2
 250         CALL SSSAVE(ISTP2,WID,IDW,ISBT2,0,0,0)
 251       END IF
 252 C
 253 C          t_i --> b_i + H+ decays
 254 C
 255       IF (AMT1SS.GT.(AMB1SS+AMHC)) THEN
 256         A=G/SR2/AMW*(AMTP*AMBT*(COTB+TANB)*SINT*SINB+
 257      $(AMBT**2*TANB+AMTP**2*COTB-AMW**2*SIN(2*BETA))*
 258      $COST*COSB-AMTP*(TWOM1-AAT*COTB)*SINT*COSB-AMBT*
 259      $(TWOM1-AAB*TANB)*SINB*COST)
 260         AS=A*A
 261         WID=AS*SQRT(SSXLAM(AMT1SS**2,AMB1SS**2,AMHC**2))/
 262      $      16./PI/AMT1SS**3
 263         CALL SSSAVE(ISTP1,WID,ISHC,ISBT1,0,0,0)
 264       END IF
 265 C
 266       IF (AMT1SS.GT.(AMB2SS+AMHC)) THEN
 267         A=G/SR2/AMW*(-AMTP*AMBT*(COTB+TANB)*SINT*COSB+
 268      $(AMBT**2*TANB+AMTP**2*COTB-AMW**2*SIN(2*BETA))*
 269      $COST*SINB-AMTP*(TWOM1-AAT*COTB)*SINT*SINB+AMBT*
 270      $(TWOM1-AAB*TANB)*COST*COSB)
 271         AS=A*A
 272         WID=AS*SQRT(SSXLAM(AMT1SS**2,AMB2SS**2,AMHC**2))/
 273      $      16./PI/AMT1SS**3
 274         CALL SSSAVE(ISTP1,WID,ISHC,ISBT2,0,0,0)
 275       END IF
 276 C
 277       IF (AMT2SS.GT.(AMB1SS+AMHC)) THEN
 278         A=G/SR2/AMW*(-AMTP*AMBT*(COTB+TANB)*COST*SINT+
 279      $(AMBT**2*TANB+AMTP**2*COTB-AMW**2*SIN(2*BETA))*
 280      $SINT*COSB+AMTP*(TWOM1-AAT*COTB)*COST*COSB-AMBT*
 281      $(TWOM1-AAB*TANB)*SINT*SINB)
 282         AS=A*A
 283         WID=AS*SQRT(SSXLAM(AMT2SS**2,AMB1SS**2,AMHC**2))/
 284      $      16./PI/AMT2SS**3
 285         CALL SSSAVE(ISTP2,WID,ISHC,ISBT1,0,0,0)
 286       END IF
 287 C
 288       IF (AMT2SS.GT.(AMB2SS+AMHC)) THEN
 289         A=G/SR2/AMW*(AMTP*AMBT*(COTB+TANB)*COST*COSB+
 290      $(AMBT**2*TANB+AMTP**2*COTB-AMW**2*SIN(2*BETA))*
 291      $SINT*SINB+AMTP*(TWOM1-AAT*COTB)*SINB*COST+AMBT*
 292      $(TWOM1-AAB*TANB)*COSB*SINT)
 293         AS=A*A
 294         WID=AS*SQRT(SSXLAM(AMT2SS**2,AMB2SS**2,AMHC**2))/
 295      $      16./PI/AMT2SS**3
 296         CALL SSSAVE(ISTP2,WID,ISHC,ISBT2,0,0,0)
 297       END IF
 298 C
 299 C
 300 C          stop_2 -> top + zino_i
 301 C
 302       DO 500 IZ=1,4
 303         MZIZ=ABS(AMZISS(IZ))
 304         SNZI=SIGN(1.,AMZISS(IZ))
 305         IF (SNZI.EQ.1.) THEN
 306            THIZ=0.
 307         ELSE
 308            THIZ=1.
 309         END IF
 310         ZAUIZ=ZI**(THIZ-1.)*SNZI
 311      $  *(-G/SR2*ZMIXSS(3,IZ)-GP/3./SR2*ZMIXSS(4,IZ))
 312         ZBUIZ=ZI**(THIZ-1.)*4*GP*ZMIXSS(4,IZ)/3./SR2
 313         ZPP=ZI**THIZ
 314         ZPM=(-ZI)**THIZ
 315         ZA=((ZI*ZAUIZ-ZPP*FT*ZMIXSS(1,IZ))*SINT
 316      $  +(ZI*ZBUIZ-ZPM*FT*ZMIXSS(1,IZ))*COST)/2.
 317         ZB=((-ZI*ZAUIZ-ZPP*FT*ZMIXSS(1,IZ))*SINT
 318      $  +(ZI*ZBUIZ+ZPM*FT*ZMIXSS(1,IZ))*COST)/2.
 319         AS=ZA*CONJG(ZA)
 320         BS=ZB*CONJG(ZB)
 321         IF (AMT2SS.GT.(AMTP+MZIZ)) THEN
 322           WID=(AS*(AMT2SS**2-(AMTP+MZIZ)**2)+BS
 323      $    *(AMT2SS**2-(AMTP-MZIZ)**2))/8./PI/AMT2SS
 324      $    *SQRT(SSXLAM(1.,AMTP**2/AMT2SS**2,MZIZ**2/AMT2SS**2))
 325           CALL SSSAVE(ISTP2,WID,ISZIZ(IZ),IDTP,0,0,0)
 326         END IF
 327 500   CONTINUE
 328 C
 329 C          Normalize branching ratios
 330 C
 331        CALL SSNORM(ISTP1)
 332        CALL SSNORM(ISTP2)
 333 C
 334        RETURN
 335        END