]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | REAL FUNCTION SSPOLE(MGMS,MUSQ,AS) | |
3 | C********************************************************************* | |
4 | C* Computes the on-shell (pole) gluino mass for given running (MSbar)* | |
5 | C* gluino mass, defined at scale MUSQ, and given alpha_s (AS). The * | |
6 | C* squark masses are stored in the SQUARK COMMON block. * | |
7 | C* This function needs the complex functions B0 and B1. * | |
8 | C* Contributed by M. Drees; modified by H. Baer * | |
9 | C * | |
10 | C Version 7.30: Cast COMPLEX*16 to REAL*8 in standard way. :-( * | |
11 | C********************************************************************* | |
12 | #if defined(CERNLIB_IMPNONE) | |
13 | IMPLICIT NONE | |
14 | #endif | |
15 | #include "isajet/sssm.inc" | |
16 | #include "isajet/sspar.inc" | |
17 | #include "isajet/sstmp.inc" | |
18 | #include "isajet/ssinf.inc" | |
19 | REAL MGMS,MUSQ,AS,MGSQ,FAC | |
20 | DOUBLE PRECISIONDMUSQ,DFAC | |
21 | COMPLEX*16 SSB1 | |
22 | DMUSQ=MUSQ | |
23 | XLAM = LOG(DMUSQ) | |
24 | MGSQ = MGMS*MGMS | |
25 | C | |
26 | C Cast COMPLEX*16 to REAL*8: | |
27 | C | |
28 | DFAC = -8*SSB1(MGSQ,0.,AMULSS) - (SSB1(MGSQ,AMTP,AMT1SS)+ | |
29 | $SSB1(MGSQ,AMTP,AMT2SS)+SSB1(MGSQ,4.0,AMB1SS)+ | |
30 | $SSB1(MGSQ,4.0,AMB2SS)) | |
31 | DFAC = DFAC + 12.D0 + 9.D0*LOG(DMUSQ/MGSQ) | |
32 | FAC=DFAC | |
33 | SSPOLE = MGMS*(1.0 + .0796*AS*FAC ) | |
34 | RETURN | |
35 | END |