]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | REAL FUNCTION SSWZF7(SS) | |
3 | C----------------------------------------------------------------------- | |
4 | C SSWZBF: wiss -> zjss f fbar | |
5 | C Drees' function for charged Higgs/sfermion interference | |
6 | C----------------------------------------------------------------------- | |
7 | #if defined(CERNLIB_IMPNONE) | |
8 | IMPLICIT NONE | |
9 | #endif | |
10 | #include "isajet/sssm.inc" | |
11 | #include "isajet/sspar.inc" | |
12 | #include "isajet/sstmp.inc" | |
13 | C | |
14 | REAL SS | |
15 | DOUBLE PRECISION S,MS,MW,MZ,MH,AL,BE,SN,RES, | |
16 | $EQ,Q,MUS,XL,TERM | |
17 | C | |
18 | S=SS | |
19 | MW=TMP(1) | |
20 | MZ=TMP(2) | |
21 | MH=TMP(3) | |
22 | MS=TMP(4) | |
23 | AL=TMP(5) | |
24 | BE=TMP(6) | |
25 | SN=TMP(7) | |
26 | C | |
27 | EQ=(S+MW**2-MZ**2)/2.D0/MW | |
28 | Q=SQRT(MAX(0.D0,EQ**2-S)) | |
29 | MUS=S+MS**2-MZ**2 | |
30 | XL=LOG((MW*(EQ+Q)-MUS)/(MW*(EQ-Q)-MUS)) | |
31 | TERM=BE*S*MS**2+SN*AL*MW*MZ*S | |
32 | RES=(S*Q*BE/2.D0+TERM*XL/4.D0/MW)/(S-MH**2) | |
33 | SSWZF7=RES | |
34 | RETURN | |
35 | END |