]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | DOUBLE PRECISION FUNCTION SSHWW2(XX) | |
3 | C----------------------------------------------------------------------- | |
4 | C SSHWW: hi -> z + z* | |
5 | C Bisset's FUNZZ | |
6 | C----------------------------------------------------------------------- | |
7 | #if defined(CERNLIB_IMPNONE) | |
8 | IMPLICIT NONE | |
9 | #endif | |
10 | #include "isajet/sspar.inc" | |
11 | #include "isajet/sssm.inc" | |
12 | #include "isajet/sstmp.inc" | |
13 | C | |
14 | DOUBLE PRECISION XX,EPLN,PROP,TEMP,FN,AAA,MZ,DELTAZ,MH | |
15 | C | |
16 | C | |
17 | MZ=AMZ | |
18 | DELTAZ=GAMZ | |
19 | MH=TMP(1) | |
20 | C | |
21 | EPLN=MZ/MH | |
22 | PROP=(1.D0-XX)**2 | |
23 | PROP=PROP+(EPLN**2)*DELTAZ**2/MH**2 | |
24 | TEMP=XX**2-12*XX*EPLN**2 | |
25 | TEMP=TEMP+8*EPLN**2+12*EPLN**4 | |
26 | AAA=XX**2-4*EPLN**2 | |
27 | IF(AAA.LT.0) THEN | |
28 | AAA=0 | |
29 | ENDIF | |
30 | FN=TEMP*SQRT(AAA)/PROP | |
31 | SSHWW2=FN | |
32 | RETURN | |
33 | END |