]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | REAL FUNCTION SSWWF1(EE) | |
3 | C----------------------------------------------------------------------- | |
4 | C SSWZBF: wiss -> zjss f fbar | |
5 | C Baer's WIWFUN | |
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 | DOUBLE PRECISION E,MW2,MW1,MZ,SN,XX,YY,T1,T2,T3,T4 | |
15 | REAL EE | |
16 | C | |
17 | E=EE | |
18 | MW2=TMP(1) | |
19 | MW1=TMP(2) | |
20 | MZ=TMP(3) | |
21 | SN=TMP(4) | |
22 | XX=TMP(5) | |
23 | YY=TMP(6) | |
24 | C | |
25 | T1=DSQRT(E**2-MW1**2)/(MW2**2+MW1**2-MZ**2-2*MW2*E)**2 | |
26 | T2=3*E*MW2*(MW2**2+MW1**2)-2*MW2**2*MW1**2-4*MW2**2*E*E | |
27 | T3=2*E*MW2-MW2**2-MW1**2 | |
28 | T4=T1*((XX**2+YY**2)*T2-3*SN*(XX**2-YY**2)*MW2*MW1*T3) | |
29 | SSWWF1=T4 | |
30 | RETURN | |
31 | END |