ISAJET/isasusy/sshcc.F

   1 #include "isajet/pilot.h"
   2       SUBROUTINE SSHCC
   3 C-----------------------------------------------------------------------
   4 C     Calculates the decay widths of all neutral Higgses into all
   5 C     possible pairs of charginos.
   6 C
   7 C     Bisset's CHGINO
   8 C-----------------------------------------------------------------------
   9 #if defined(CERNLIB_IMPNONE)
  10       IMPLICIT NONE
  11 #endif
  12 #include "isajet/sspar.inc"
  13 #include "isajet/sssm.inc"
  14 #include "isajet/sstype.inc"
  15 C
  16       DOUBLE PRECISION PI,SR2,XM,THETX,YM,THETY,SGL,CGL,SGR,CGR
  17      $,MW1,MW2,THETM,THETP,G2,GP2,BETA,ALPHA,T1,MH,M1,M2,LAMB
  18      $,DWID,TEMP,TEMPXY
  19       DOUBLE PRECISION MHI(3),IDHI(3),SHP(3),SHM(3),SH(3),PH(3)
  20       DOUBLE PRECISION SSDLAM
  21       REAL WID
  22       INTEGER NUMH,IDHHA
  23 C
  24       PI=4.*ATAN(1.D0)
  25       SR2=SQRT(2.D0)
  26       XM=1./TAN(GAMMAL)
  27       THETX=SIGN(1.D0,XM)
  28       YM=1./TAN(GAMMAR)
  29       THETY=SIGN(1.D0,YM)
  30       SGL=1/(DSQRT(1+XM**2))
  31       CGL=SGL*XM
  32       SGR=1/(DSQRT(1+YM**2))
  33       CGR=SGR*YM
  34       MW1=DBLE(ABS(AMW1SS))
  35       MW2=DBLE(ABS(AMW2SS))
  36       THETM=SIGN(1.,AMW1SS)
  37       THETP=SIGN(1.,AMW2SS)
  38       G2=4*PI*ALFAEM/SN2THW
  39       GP2=4*PI*ALFAEM/(1-SN2THW)
  40       BETA=ATAN(1.0/RV2V1)
  41       ALPHA=ALFAH
  42 C
  43 C          Arrays for loops
  44 C
  45       MHI(1)=AMHL
  46       MHI(2)=AMHH
  47       MHI(3)=AMHA
  48       IDHI(1)=ISHL
  49       IDHI(2)=ISHH
  50       IDHI(3)=ISHA
  51 C          The following came from Bisset's MASZIN, but with L,H,P
  52 C          replaced by a generic H and a subscript.
  53       TEMPXY=0.5*THETX*THETY*(-THETP)
  54       SHP(1)=SIN(ALPHA)*CGR*SGL+COS(ALPHA)*CGL*SGR
  55       SHP(1)=SHP(1)*TEMPXY
  56       SHM(1)=SIN(ALPHA)*SGR*CGL+COS(ALPHA)*SGL*CGR
  57       SHM(1)=SHM(1)*0.5*THETM
  58       SH(1)=-THETX*SGR*SGL*SIN(ALPHA)*THETM
  59       PH(1)=-SH(1)
  60       T1=THETX*CGL*CGR*COS(ALPHA)*THETM
  61       SH(1)=SH(1)+T1
  62       PH(1)=PH(1)-T1
  63       T1=THETY*SGL*SGR*COS(ALPHA)*THETP
  64       SH(1)=SH(1)-T1
  65       PH(1)=PH(1)-T1
  66       T1=THETY*CGL*CGR*SIN(ALPHA)*THETP
  67       SH(1)=SH(1)+T1
  68       PH(1)=PH(1)+T1
  69       SH(1)=0.5*SH(1)
  70       PH(1)=0.5*PH(1)
  71       SHP(2)=COS(ALPHA)*CGR*SGL-SIN(ALPHA)*CGL*SGR
  72       SHP(2)=SHP(2)*TEMPXY
  73       SHM(2)=COS(ALPHA)*SGR*CGL-SIN(ALPHA)*SGL*CGR
  74       SHM(2)=SHM(2)*0.5*THETM
  75       SH(2)=-THETX*SGR*SGL*COS(ALPHA)*THETM
  76       PH(2)=-SH(2)
  77       T1=THETX*CGL*CGR*SIN(ALPHA)*THETM
  78       SH(2)=SH(2)-T1
  79       PH(2)=PH(2)+T1
  80       T1=THETY*SGL*SGR*SIN(ALPHA)*THETP
  81       SH(2)=SH(2)+T1
  82       PH(2)=PH(2)+T1
  83       T1=THETY*CGL*CGR*COS(ALPHA)*THETP
  84       SH(2)=SH(2)+T1
  85       PH(2)=PH(2)+T1
  86       SH(2)=0.5*SH(2)
  87       PH(2)=0.5*PH(2)
  88       SHP(3)=SIN(BETA)*CGR*SGL+COS(BETA)*CGL*SGR
  89       SHP(3)=SHP(3)*0.5*THETX*THETY*(-THETP)
  90       SHM(3)=SIN(BETA)*SGR*CGL+COS(BETA)*SGL*CGR
  91       SHM(3)=SHM(3)*0.5*THETM
  92       SH(3)=-THETX*SGR*SGL*SIN(BETA)*THETM
  93       PH(3)=-SH(3)
  94       T1=THETX*CGL*CGR*COS(BETA)*THETM
  95       SH(3)=SH(3)+T1
  96       PH(3)=PH(3)-T1
  97       T1=THETY*SGL*SGR*COS(BETA)*THETP
  98       SH(3)=SH(3)+T1
  99       PH(3)=PH(3)+T1
 100       T1=THETY*CGL*CGR*SIN(BETA)*THETP
 101       SH(3)=SH(3)-T1
 102       PH(3)=PH(3)-T1
 103       SH(3)=0.5*SH(3)
 104       PH(3)=0.5*PH(3)
 105 C
 106 C          Loop over neutral Higgs
 107 C
 108       DO 100 NUMH=1,3
 109         MH=MHI(NUMH)
 110         IDHHA=IDHI(NUMH)
 111 C          w1 + w1
 112         M1=ABS(AMW1SS)
 113         M2=M1
 114         IF(MH.GT.M1+M2) THEN
 115           LAMB=SSDLAM(MH**2,M1**2,M2**2)
 116           TEMP=1-4*M1**2/MH**2
 117           DWID=G2*MH*SHM(NUMH)**2/(4.0*PI)
 118           DWID=DWID*SQRT(TEMP**3)
 119           WID=DWID
 120           CALL SSSAVE(IDHHA,WID,ISW1,-ISW1,0,0,0)
 121         ENDIF
 122 C          w2 + w2
 123         M1=ABS(AMW2SS)
 124         M2=M1
 125         IF(MH.GT.M1+M2) THEN
 126           TEMP=1-4*M1**2/MH**2
 127           DWID=G2*MH*SHP(NUMH)**2/(4*PI)
 128           DWID=DWID*SQRT(TEMP**3)
 129           WID=DWID
 130           CALL SSSAVE(IDHHA,WID,ISW2,-ISW2,0,0,0)
 131         ENDIF
 132 C          w1 + w2
 133         M1=ABS(AMW1SS)
 134         M2=ABS(AMW2SS)
 135         IF(MH.GT.M1+M2) THEN
 136           LAMB=SSDLAM(MH**2,M1**2,M2**2)
 137           DWID=PH(NUMH)**2*(MH**2-(M1-M2)**2)
 138           DWID=DWID+SH(NUMH)**2*(MH**2-(M1+M2)**2)
 139           DWID=DWID*G2*SQRT(LAMB)/(16.0*PI*(MH**3))
 140           WID=DWID
 141           CALL SSSAVE(IDHHA,WID,ISW1,-ISW2,0,0,0)
 142           CALL SSSAVE(IDHHA,WID,-ISW1,ISW2,0,0,0)
 143         ENDIF
 144 100   CONTINUE
 145 C
 146       RETURN
 147       END