]>
Commit | Line | Data |
---|---|---|
3820ca8e | 1 | |
2 | CDECK ID>, HWUPCM. | |
3 | ||
4 | *CMZ :- -26/04/91 11.11.56 by Bryan Webber | |
5 | ||
6 | *-- Author : Bryan Webber | |
7 | ||
8 | C----------------------------------------------------------------------- | |
9 | ||
10 | FUNCTION HWUPCM(EM0,EM1,EM2) | |
11 | ||
12 | C----------------------------------------------------------------------- | |
13 | ||
14 | C C.M. MOMENTUM FOR DECAY MASSES EM0 -> EM1 + EM2 | |
15 | ||
16 | C SET TO -1 BELOW THRESHOLD | |
17 | ||
18 | C----------------------------------------------------------------------- | |
19 | ||
20 | DOUBLE PRECISION HWUPCM,EM0,EM1,EM2,EMS,EMD | |
21 | ||
22 | EMS=ABS(EM1+EM2) | |
23 | ||
24 | EMD=ABS(EM1-EM2) | |
25 | ||
26 | IF (EM0.LT.EMS.OR.EM0.LT.EMD) THEN | |
27 | ||
28 | HWUPCM=-1. | |
29 | ||
30 | ELSEIF (EM0.EQ.EMS.OR.EM0.EQ.EMD) THEN | |
31 | ||
32 | HWUPCM=0. | |
33 | ||
34 | ELSE | |
35 | ||
36 | HWUPCM=SQRT((EM0+EMD)*(EM0-EMD)* | |
37 | ||
38 | & (EM0+EMS)*(EM0-EMS))*.5/EM0 | |
39 | ||
40 | ENDIF | |
41 | ||
42 | END |