]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | REAL FUNCTION SSLRT1(SS) | |
3 | C----------------------------------------------------------------------- | |
4 | C SSLRT1: l_R -> l+tau+stau_1 | |
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 | REAL SS | |
13 | DOUBLE PRECISION ETMX,ETMN,S,MT1,MT,MLR,BEZI,BEZJ,TM, | |
14 | ,AI,AJ,BI,BJ,MZI,MZJ,SNZI,SNZJ,XL,BK1,BK2,BK3,BK,WID,SSDLAM | |
15 | S=SS | |
16 | MT1=AML1SS | |
17 | MT=AMTAU | |
18 | MLR=TMP(1) | |
19 | BEZI=TMP(2) | |
20 | BEZJ=TMP(3) | |
21 | AI=TMP(4) | |
22 | AJ=TMP(5) | |
23 | BI=TMP(6) | |
24 | BJ=TMP(7) | |
25 | MZI=ABS(TMP(8)) | |
26 | MZJ=ABS(TMP(9)) | |
27 | SNZI=SIGN(1.0,TMP(8)) | |
28 | SNZJ=SIGN(1.0,TMP(9)) | |
29 | TM=SSDLAM(S,MT**2,MT1**2) | |
30 | XL=DSQRT(MAX(0.D0,TM)) | |
31 | ETMN=(S+MT**2-MT1**2-XL*(MLR**2-S)/(MLR**2+S))*(MLR**2+S)/ | |
32 | , (2*S)/(2*MLR) | |
33 | ETMX=(S+MT**2-MT1**2+XL*(MLR**2-S)/(MLR**2+S))*(MLR**2+S)/ | |
34 | , (2*S)/(2*MLR) | |
35 | BK1=-(ETMX-ETMN)*((ETMX+ETMN)*MLR*S- | |
36 | , (S+MT**2-MT1**2)*MLR**2)/2.D0 | |
37 | BK2=(ETMX-ETMN)*((ETMX+ETMN)*MLR-S-MT**2+MT1**2)/2.D0 | |
38 | BK3=SNZJ*BI*AJ*MZJ+SNZI*BJ*AI*MZI | |
39 | BK=BI*BJ*BK1+AI*AJ*MZI*MZJ*SNZI*SNZJ*BK2+BK3*MT*(MLR**2-S)* | |
40 | , (ETMX-ETMN)/2.D0 | |
41 | WID=BEZI*BEZJ*BK/(S-MZI**2)/(S-MZJ**2) | |
42 | SSLRT1=WID | |
43 | RETURN | |
44 | END |