ISAJET/isasusy/sswzbf.F

   1 #include "isajet/pilot.h"
   2         SUBROUTINE SSWZBF
   3 C-----------------------------------------------------------------------
   4 C       This subroutine gives chargino (wi) and neutralino (zi)
   5 C       decays according to Baer, et al.
   6 C       Valid for all scalar masses (functions in double precision)
   7 C       Includes Higgs sector radiative corrections (Aug. 31)
   8 C
   9 C       Auxiliary functions are called SSWxyi, SSZxyi, where normally
  10 C       x indicates the SUSY particle, y the SM particle(s), and i is
  11 C       a counter.
  12 C
  13 C       Ver. 7.07: For w_i -> z_j or z_i -> w_j decays, require that
  14 C                  decay be allowed by a factor FUDGE = 1.0
  15 C       Ver. 7.28: Calculate full large tan(beta) decays.
  16 C                  Calculate tau polarizations.
  17 C       Ver. 7.33: Add GMSB modes.
  18 C       Ver. 7.41: Add decay matrix elements.
  19 C                  Split real work into sszibf and sswfbf.
  20 C
  21 C       Baer's GAUGBF
  22 C
  23 C-----------------------------------------------------------------------
  24 #if defined(CERNLIB_IMPNONE)
  25       IMPLICIT NONE
  26 #endif
  27 #include "isajet/sslun.inc"
  28 #include "isajet/ssmode.inc"
  29 #include "isajet/sssm.inc"
  30 #include "isajet/sspar.inc"
  31 #include "isajet/sstype.inc"
  32 #include "isajet/sstmp.inc"
  33 #include "isajet/sspols.inc"
  34 C
  35       INTEGER IZ
  36 C
  37 C          Initialize polarizations to zero
  38 C
  39       DO 150 IZ=1,4
  40         PTAU1(IZ)=0
  41         PTAU2(IZ)=0
  42 150   CONTINUE
  43       PTAUZZ=0
  44       PTAUWZ=0
  45 C-----------------------------------------------------------------------
  46 C          Generate Neutralino zi Branching Fractions
  47 C-----------------------------------------------------------------------
  48       CALL SSZIBF
  49 C-----------------------------------------------------------------------
  50 C          Generate Chargino Branching Fractions
  51 C-----------------------------------------------------------------------
  52       CALL SSWIBF
  53 C
  54 C          Set more neutralino polarizations
  55 C
  56       IF (ABS(AMZISS(2)).GT.(AMTAU+AML1SS)) PTAUZ2(1)=PTAU1(2)
  57       IF (ABS(AMZISS(2)).GT.(AMTAU+AML2SS)) PTAUZ2(2)=PTAU2(2)
  58       IF (ABS(AMZISS(3)).GT.(AMTAU+AML1SS)) PTAUZ3(1)=PTAU1(3)
  59       IF (ABS(AMZISS(3)).GT.(AMTAU+AML2SS)) PTAUZ3(2)=PTAU2(3)
  60       IF (ABS(AMZISS(4)).GT.(AMTAU+AML1SS)) PTAUZ4(1)=PTAU1(4)
  61       IF (ABS(AMZISS(4)).GT.(AMTAU+AML2SS)) PTAUZ4(2)=PTAU2(4)
  62 C
  63       RETURN
  64       END