]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | REAL FUNCTION SSZZF4(E) | |
3 | C----------------------------------------------------------------------- | |
4 | C Z_I -> Z_J +B +BBAR VIA HIGGS | |
5 | C----------------------------------------------------------------------- | |
6 | #if defined(CERNLIB_IMPNONE) | |
7 | IMPLICIT NONE | |
8 | #endif | |
9 | #include "isajet/sssm.inc" | |
10 | #include "isajet/sspar.inc" | |
11 | #include "isajet/sstmp.inc" | |
12 | C | |
13 | REAL E | |
14 | DOUBLE PRECISION MZ1,MZ2,SN,EE,MH1,MH2,P,MF,BF,TM | |
15 | C | |
16 | MZ2=TMP(1) | |
17 | MZ1=TMP(2) | |
18 | MH1=TMP(3) | |
19 | SN=TMP(4) | |
20 | MH2=TMP(5) | |
21 | MF=TMP(6) | |
22 | EE=E | |
23 | C | |
24 | P=SQRT(MAX(0.D0,EE**2-MZ1**2)) | |
25 | TM=1.D0-4*MF**2/(MZ2**2+MZ1**2-2*E*MZ2) | |
26 | BF=DSQRT(MAX(0.D0,TM)) | |
27 | SSZZF4=P*BF*(MZ2**2+MZ1**2-2*MZ2*EE-2*MF**2)* | |
28 | $(MZ2*EE+SN*MZ2*MZ1)/(MZ2**2+MZ1**2-2*MZ2*EE-MH1**2)/ | |
29 | $(MZ2**2+MZ1**2-2*MZ2*EE-MH2**2) | |
30 | RETURN | |
31 | END |