]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | REAL FUNCTION SSZZF1(X) | |
3 | C----------------------------------------------------------------------- | |
4 | C SSWZBF: ziss -> zjss f fbar | |
5 | C Baer's TFUNC | |
6 | C----------------------------------------------------------------------- | |
7 | #if defined(CERNLIB_IMPNONE) | |
8 | IMPLICIT NONE | |
9 | #endif | |
10 | #include "isajet/sslun.inc" | |
11 | #include "isajet/sssm.inc" | |
12 | #include "isajet/sspar.inc" | |
13 | #include "isajet/sstmp.inc" | |
14 | C | |
15 | REAL X | |
16 | DOUBLE PRECISION N,S,Q,D,SQBRKT,TERM1,QS | |
17 | DOUBLE PRECISION MZ2,MZ1,M,DSN | |
18 | C | |
19 | MZ2=TMP(1) | |
20 | MZ1=TMP(2) | |
21 | M=TMP(3) | |
22 | DSN=TMP(4) | |
23 | C | |
24 | N=MZ1**2/MZ2**2 | |
25 | S=M**2/MZ2**2 | |
26 | Q=X*(1.D0-N)/2. | |
27 | QS=Q**2 | |
28 | D=(S-2*S*Q-N)/((1.D0-2*Q)*(S-2*Q-N)) | |
29 | IF(D.LE.0.) THEN | |
30 | WRITE(LOUT,*) 'ERROR IN SSZZF1: D,S,Q,N=',D,S,Q,N | |
31 | SSZZF1=0 | |
32 | RETURN | |
33 | END IF | |
34 | SQBRKT=-Q*(1.D0-2*Q-N)/(1.D0-2*Q)-(2*Q-S+N)/2.D0*DLOG(D) | |
35 | TERM1=QS*(1.D0-2*Q-N)**2/(1.D0-2*Q-S)**2/(1.D0-2*Q) | |
36 | $+DSQRT(N)/2.D0/(1.D0-2*Q-S)*SQBRKT*DSN | |
37 | SSZZF1=(1.D0-N)/2.D0*TERM1 | |
38 | RETURN | |
39 | END |