]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | REAL FUNCTION SSWZF3(QQ) | |
3 | C----------------------------------------------------------------------- | |
4 | C SSWZBF: wiss -> zjss f fbar | |
5 | C Baer's PHIFUN | |
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 QQ | |
15 | DOUBLE PRECISION M1,M2,M3,T,B,XLOG,SQBKT,Q,PI | |
16 | DATA PI/3.14159265D0/ | |
17 | C | |
18 | Q=QQ | |
19 | M1=TMP(1) | |
20 | M2=TMP(2) | |
21 | M3=TMP(3) | |
22 | C | |
23 | T=M2**2*(M1-2*Q)-M1*M3**2 | |
24 | B=(M1-2*Q)*(M2**2-2*M1*Q-M3**2) | |
25 | XLOG=DLOG(T/B) | |
26 | SQBKT=-Q*(M1**2-M3**2-2*M1*Q)/M1/(M1-2*Q)- | |
27 | $(2*M1*Q-M2**2+M3**2)*XLOG/2./M1 | |
28 | SSWZF3=.5*PI**2*M1*M3*SQBKT/(M1**2-M2**2-2*M1*Q) | |
29 | RETURN | |
30 | END |