| 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 |
git://git.uio.no / u / mrichter / AliRoot.git / blame