1 #include "isajet/pilot.h"
2 SUBROUTINE DBLVEC(P,DP)
4 C Calculate double precision vector DP for 5-vector P.
5 C Exact components are 1,2,5 and larger of +,-
6 C Ver 6.44: Always use this, even if IF=SINGLE.
8 #if defined(CERNLIB_IMPNONE)
12 DOUBLE PRECISION DP(5),DPPL,DPMN
17 IF(DP(4)+ABS(DP(3)).EQ.0.) RETURN
20 DPMN=(DP(1)**2+DP(2)**2+DP(5)**2)/DPPL
23 DPPL=(DP(1)**2+DP(2)**2+DP(5)**2)/DPMN
25 DP(3)=0.5D0*(DPPL-DPMN)
26 DP(4)=0.5D0*(DPPL+DPMN)