pxuvec.f

Go to the documentation of this file.

CDECK  ID>, PXUVEC.
*
******............................................................******
       SUBROUTINE pxuvec(NTRAK,PP,PU,IERR)
******............................................................******
*
*** Routine to calculate unit vectors PU of all particles PP
      IMPLICIT NONE
      INTEGER mxtrak
      parameter(mxtrak=200)
      INTEGER ntrak, ierr
      REAL pp(4,mxtrak)
      REAL pu(3,mxtrak)
      INTEGER n,mu
      REAL mag
       DO 100 n=1,ntrak
          mag=0.0
          DO 110 mu=1,3
             mag=mag+pp(mu,n)**2
110       CONTINUE
          mag=sqrt(mag)
          IF (mag.EQ.0.0) THEN
             WRITE(*,*)' PXCONE: An input particle has zero mod(p)'
             ierr=-1
             RETURN
          ENDIF
          DO 120 mu=1,3
           pu(mu,n)=pp(mu,n)/mag
120       CONTINUE
100    CONTINUE
       RETURN
       END