HERWIG/src/hwdrm3.f

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