| 0795afa3 |
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 |
git://git.uio.no / u / mrichter / AliRoot.git / blame