PYTHIA/pythia/rkzpr.F

   1
   2 *==================================================================
   3
   4       FUNCTION RKZPR(L,Q1,Q2)
   5       IMPLICIT DOUBLE PRECISION(A-H,O-Z)
   6       COMPLEX RKZPR
   7       DIMENSION Q1(0:4),Q2(0:4)
   8       IF(IABS(L).NE.1) THEN
   9         WRITE(6,*) ' RKZPR: ERROR   L=',L
  10         STOP
  11       ENDIF
  12 C...Introduce cutoff to check that R1 and R2 not zero.
  13       R1=DSQRT(MAX(1D-10,Q1(0)-Q1(1)))
  14       R2=DSQRT(MAX(1D-10,Q2(0)-Q2(1)))
  15       RKZPR=CMPLX(SNGL(Q1(2)),SNGL(Q1(3)))*R2/R1
  16      .     -CMPLX(SNGL(Q2(2)),SNGL(Q2(3)))*R1/R2
  17       IF(L.EQ.-1) RKZPR=-CONJG(RKZPR)
  18       END