Bugfix in AliPoints2Memory
[u/mrichter/AliRoot.git] / GEANT321 / gheisha / dpdnuc.F
1 *
2 * $Id$
3 *
4 * $Log$
5 * Revision 1.1.1.1  1995/10/24 10:21:06  cernlib
6 * Geant
7 *
8 *
9 #include "geant321/pilot.h"
10 *CMZ :  3.21/02 29/03/94  15.41.40  by  S.Giani
11 *-- Author :
12       FUNCTION DPDNUC(A,B,C)
13 C
14 C CALLED BY : PHPNUC
15 C ORIGIN    : H.FESEFELDT
16 C
17 #if !defined(CERNLIB_SINGLE)
18       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
19 #endif
20 C
21 C     DPDK = SQRT(A*A+(B*B-C*C)**2/(A*A) - 2.0*(B*B+C*C))/2.0
22       A2 = A*A
23       B2 = B*B
24       C2 = C*C
25       IF(A2) 21,21,61
26    61 CONTINUE
27       ARG=A2+(B2-C2)**2/A2-2.0*(B2+C2)
28       IF (ARG) 21,21,31
29    21 DPDNUC=0.0
30       GOTO 41
31    31 CONTINUE
32       DPDNUC = 0.5*SQRT(A2 + (B2-C2)**2/A2 - 2.0*(B2+C2))
33    41 CONTINUE
34       RETURN
35       END