]>
Commit | Line | Data |
---|---|---|
3820ca8e | 1 | |
2 | CDECK ID>, HWDRM4. | |
3 | ||
4 | *CMZ :- -20/07/99 10:56:12 by Peter Richardson | |
5 | ||
6 | *-- Author : Peter Richardson | |
7 | ||
8 | C----------------------------------------------------------------------- | |
9 | ||
10 | FUNCTION HWDRM4(X,Y,MA,MB,MC,MD,MR1,MR2,GAM1,GAM2,A1,A2,B1,B2) | |
11 | ||
12 | C----------------------------------------------------------------------- | |
13 | ||
14 | C Function to compute the interference part of a 3-body | |
15 | ||
16 | C R-parity decay | |
17 | ||
18 | C----------------------------------------------------------------------- | |
19 | ||
20 | IMPLICIT NONE | |
21 | ||
22 | DOUBLE PRECISION X,Y,MA,MB,MC,MD,A1,A2,B1,B2,HWDRM4,MR1,MR2,GAM1 | |
23 | ||
24 | & ,GAM2 | |
25 | ||
26 | C | |
27 | ||
28 | HWDRM4 = 2*((GAM1*GAM2*MR1*MR2 + (X - MR1**2)*(Y - MR2**2))* | |
29 | ||
30 | & (A2*B1*MC*MD*(X - MA**2 - MB**2) + | |
31 | ||
32 | & A1*A2*MA*MC*(X + Y - MA**2 - MC**2) + | |
33 | ||
34 | & A1*B2*MA*MD*(Y - MB**2 - MC**2) + | |
35 | ||
36 | & B1*B2*(X*Y - MA**2*MC**2 - MB**2*MD**2)))/ | |
37 | ||
38 | & (((X-MR1**2)**2+GAM1**2*MR1**2)*((Y-MR2**2)**2+GAM2**2*MR2**2)) | |
39 | ||
40 | END |