1 #include "isajet/pilot.h"
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)
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
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.
23 C-----------------------------------------------------------------------
24 #if defined(CERNLIB_IMPNONE)
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"
37 C Initialize polarizations to zero
45 C-----------------------------------------------------------------------
46 C Generate Neutralino zi Branching Fractions
47 C-----------------------------------------------------------------------
49 C-----------------------------------------------------------------------
50 C Generate Chargino Branching Fractions
51 C-----------------------------------------------------------------------
54 C Set more neutralino polarizations
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)