]>
Commit | Line | Data |
---|---|---|
3820ca8e | 1 | |
2 | CDECK ID>, HWURAP. | |
3 | ||
4 | *CMZ :- -26/04/91 11.11.56 by Bryan Webber | |
5 | ||
6 | *-- Author : Bryan Webber | |
7 | ||
8 | C----------------------------------------------------------------------- | |
9 | ||
10 | FUNCTION HWURAP(P) | |
11 | ||
12 | C----------------------------------------------------------------------- | |
13 | ||
14 | C LONGITUDINAL RAPIDITY (SET TO +/-1000 IF TOO LARGE) | |
15 | ||
16 | C----------------------------------------------------------------------- | |
17 | ||
18 | DOUBLE PRECISION HWURAP,EMT2,P(5),ZERO | |
19 | ||
20 | PARAMETER (ZERO=0.D0) | |
21 | ||
22 | EMT2=P(1)**2+P(2)**2+P(5)**2 | |
23 | ||
24 | IF (P(3).GT.ZERO) THEN | |
25 | ||
26 | IF (EMT2.EQ.ZERO) THEN | |
27 | ||
28 | HWURAP=1000. | |
29 | ||
30 | ELSE | |
31 | ||
32 | HWURAP= 0.5*LOG((P(3)+P(4))**2/EMT2) | |
33 | ||
34 | ENDIF | |
35 | ||
36 | ELSEIF (P(3).LT.ZERO) THEN | |
37 | ||
38 | IF (EMT2.EQ.ZERO) THEN | |
39 | ||
40 | HWURAP=-1000. | |
41 | ||
42 | ELSE | |
43 | ||
44 | HWURAP=-0.5*LOG((P(3)-P(4))**2/EMT2) | |
45 | ||
46 | ENDIF | |
47 | ||
48 | ELSE | |
49 | ||
50 | HWURAP=0. | |
51 | ||
52 | ENDIF | |
53 | ||
54 | END |