]>
Commit | Line | Data |
---|---|---|
3820ca8e | 1 | |
2 | CDECK ID>, HWDTWO. | |
3 | ||
4 | *CMZ :- -27/01/94 17.38.49 by Mike Seymour | |
5 | ||
6 | *-- Author : Bryan Webber & Mike Seymour | |
7 | ||
8 | C----------------------------------------------------------------------- | |
9 | ||
10 | SUBROUTINE HWDTWO(P0,P1,P2,PCM,COSTH,ZAXIS) | |
11 | ||
12 | C----------------------------------------------------------------------- | |
13 | ||
14 | C GENERATES DECAY 0 -> 1+2 | |
15 | ||
16 | C | |
17 | ||
18 | C PCM IS CM MOMENTUM | |
19 | ||
20 | C | |
21 | ||
22 | C COSTH = COS THETA IN P0 REST FRAME (>1 FOR ISOTROPIC) | |
23 | ||
24 | C IF ZAXIS=.TRUE., COS THETA IS MEASURED FROM THE ZAXIS | |
25 | ||
26 | C IF .FALSE., IT IS MEASURED FROM P0'S DIRECTION | |
27 | ||
28 | C----------------------------------------------------------------------- | |
29 | ||
30 | DOUBLE PRECISION HWRUNI,ONE,ZERO,PCM,COSTH,C,S,P0(5),P1(5),P2(5), | |
31 | ||
32 | & PP(5),R(9) | |
33 | ||
34 | LOGICAL ZAXIS | |
35 | ||
36 | EXTERNAL HWRUNI | |
37 | ||
38 | PARAMETER (ZERO=0.D0, ONE=1.D0) | |
39 | ||
40 | C--CHOOSE C.M. ANGLES | |
41 | ||
42 | C=COSTH | |
43 | ||
44 | IF (C.GT.ONE) C=HWRUNI(0,-ONE,ONE) | |
45 | ||
46 | S=SQRT(ONE-C*C) | |
47 | ||
48 | CALL HWRAZM(PCM*S,PP(1),PP(2)) | |
49 | ||
50 | C--PP IS MOMENTUM OF 2 IN C.M. | |
51 | ||
52 | PP(3)=-PCM*C | |
53 | ||
54 | PP(4)=SQRT(P2(5)**2+PCM**2) | |
55 | ||
56 | PP(5)=P2(5) | |
57 | ||
58 | C--ROTATE IF NECESSARY | |
59 | ||
60 | IF (COSTH.LE.ONE.AND..NOT.ZAXIS) THEN | |
61 | ||
62 | CALL HWUROT(P0,ONE,ZERO,R) | |
63 | ||
64 | CALL HWUROB(R,PP,PP) | |
65 | ||
66 | ENDIF | |
67 | ||
68 | C--BOOST FROM C.M. TO LAB FRAME | |
69 | ||
70 | CALL HWULOB(P0,PP,P2) | |
71 | ||
72 | CALL HWVDIF(4,P0,P2,P1) | |
73 | ||
74 | END |