GEANT321/fluka/trafo.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:20:01  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 TRAFO.FOR
  13 *COPY TRAFO
  14 *
  15 *=== trafo ============================================================*
  16 *
  17       SUBROUTINE TRAFO(GAM,BGAM,CX,CY,CZ,COD,COF,SIF,P,ECM,
  18      1PL,CXL,CYL,CZL,EL)
  19
  20 #include "geant321/dblprc.inc"
  21 #include "geant321/dimpar.inc"
  22 #include "geant321/iounit.inc"
  23 C     LORENTZ TRANSFORMATION INTO THE LAB - SYSTEM
  24       SID=SQRT(1.D0-COD*COD)
  25       PLX=P*SID*COF
  26       PLY=P*SID*SIF
  27       PCMZ=P*COD
  28       PLZ=GAM*PCMZ+BGAM*ECM
  29       PL=SQRT(PLX*PLX+PLY*PLY+PLZ*PLZ)
  30       EL=GAM*ECM+BGAM*PCMZ
  31 C     ROTATION INTO THE ORIGINAL DIRECTION
  32       COZ=PLZ/PL
  33       SIZ=SQRT(1.D0-COZ**2)
  34       CALL TTRANS(CX,CY,CZ,COZ,SIZ,SIF,COF,CXL,CYL,CZL)
  35       RETURN
  36       END