| fe4da5cc |
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 |
git://git.uio.no / u / mrichter / AliRoot.git / blame