GEANT321/fluka/kinpar.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 KINPAR.FOR
  13 *COPY KINPAR
  14 *
  15 *=== kinpar ===========================================================*
  16 *
  17       SUBROUTINE KINPAR(PXYZ,E,RCPMV)
  18
  19 #include "geant321/dblprc.inc"
  20 #include "geant321/dimpar.inc"
  21 #include "geant321/iounit.inc"
  22       COMMON /FKKINP/AMTA,AMIT,PXTA,PYTA,PZTA,PSX,PSY,PSZ,PGES2,PGES,
  23      *PGESM,PTA2,PTAM,UMO2,CCX,CCY,CCZ,CXTA,CYTA,CZTA,GA,BGA
  24 C*** CALCULATION OF KINEMATIC PARAMETERS
  25       DIMENSION PXYZ(10),E(2)
  26       PXTA=PXYZ(  2)
  27       PYTA=PXYZ(  6)
  28       PZTA=PXYZ( 10)
  29       PSX=PXYZ( 1)+P XTA
  30       PSY=PXYZ( 5)+P YTA
  31       PSZ=PXYZ( 9)+P ZTA
  32       PGES2=PSX*PSX+PSY*PSY+PSZ*PSZ
  33       PGES=SQRT(PGES2)
  34       PGESM=1.D0/PGES
  35       PTA2=PXTA*PXTA+PYTA*PYTA+PZTA*PZTA
  36       PTAM=1.D0/SQRT(PTA2)
  37       UMO2=((E(1)+E(2))**2-PGES2)/(RCPMV*RCPMV)
  38       UMO2=UMO2*0.000001D0
  39       CCX=PS X*PGESM
  40       CCY=PS Y*PGESM
  41       CCZ=PS Z*PGESM
  42       CXTA=P XTA*PTAM
  43       CYTA=P YTA*PTAM
  44       CZTA=P ZTA*PTAM
  45       GA=(E(1)+E(2))*0.001D0/(RCPMV*SQRT(UMO2))
  46       BGA=SQRT(GA*GA-1.D0)
  47       RETURN
  48       END