GEANT321/fluka/altra.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:19:58  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/02 29/03/94  15.41.44  by  S.Giani
  11 *-- Author :
  12 *$ CREATE ALTRA.FOR
  13 *COPY ALTRA
  14 *
  15 *=== altra ============================================================*
  16 *
  17       SUBROUTINE ALTRA(GA,BGX,BGY,BGZ,PCX,PCY,PCZ,EC,P,PX,PY,PZ,E)
  18
  19 #include "geant321/dblprc.inc"
  20 #include "geant321/dimpar.inc"
  21 #include "geant321/iounit.inc"
  22 C*** ARBITRARY LORENTZ TRANSFORM
  23       EP = PCX * BGX + PCY * BGY + PCZ * BGZ
  24       E  = GA * EC + EP
  25       PE = EP / (GA + 1.D0) + EC
  26       PX = PCX + BGX * PE
  27       PY = PCY + BGY * PE
  28       PZ = PCZ + BGZ * PE
  29       TMP30 = 1.D30
  30       PX = MIN(TMP30, MAX(-TMP30,PX))
  31       PY = MIN(TMP30, MAX(-TMP30,PY))
  32       PZ = MIN(TMP30, MAX(-TMP30,PZ))
  33       P  = SQRT (PX * PX + PY * PY + PZ * PZ)
  34       RETURN
  35       END