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