GEANT321/gphys/gloren.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:21:26  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/02 29/03/94  15.41.22  by  S.Giani
  11 *-- Author :
  12       SUBROUTINE GLOREN(BETA,PA,PB)
  13 C.
  14 C.    ******************************************************************
  15 C.    *                                                                *
  16 C     *       Routine to transform momentum and energy from the        *
  17 C     *       Lorentz frame A to the Lorentz frame B                   *
  18 C     *                                                                *
  19 C     *       PA(1)                                                    *
  20 C     *       PA(2)     Momentum components in frame A                 *
  21 C     *       PA(3)                                                    *
  22 C     *       PA(4)     Energy                                         *
  23 C     *       PB(..)   same quantities in frame B                      *
  24 C     *                                                                *
  25 C     *       BETA(1)    Components of velocity of frame B             *
  26 C     *       BETA(2)        as seen from frame A                      *
  27 C     *       BETA(3)                                                  *
  28 C     *       BETA(4)    1./SQRT(1.-BETA**2)                           *
  29 C.    *                                                                *
  30 C.    *    ==>Called by : GDECAY,GDECA3                                *
  31 C.    *       Author    M.Hansroul  *********                          *
  32 C.    *                                                                *
  33 C.    ******************************************************************
  34 C.
  35       DIMENSION BETA(4),PA(4),PB(4)
  36 C.
  37 C.    ------------------------------------------------------------------
  38 C.
  39       BETPA  = BETA(1)*PA(1) + BETA(2)*PA(2) + BETA(3)*PA(3)
  40       BPGAM  = (BETPA * BETA(4)/(BETA(4) + 1.) - PA(4)) * BETA(4)
  41       PB(1) = PA(1) + BPGAM  * BETA(1)
  42       PB(2) = PA(2) + BPGAM  * BETA(2)
  43       PB(3) = PA(3) + BPGAM  * BETA(3)
  44       PB(4) =(PA(4) - BETPA) * BETA(4)
  45       END
  46