Bugfix in AliPoints2Memory
[u/mrichter/AliRoot.git] / GEANT321 / gheisha / lor.F
CommitLineData
fe4da5cc 1*
2* $Id$
3*
4* $Log$
5* Revision 1.1.1.1 1995/10/24 10:21:07 cernlib
6* Geant
7*
8*
9#include "geant321/pilot.h"
10*CMZ : 3.21/02 29/03/94 15.41.38 by S.Giani
11*-- Author :
12C
13C *** VARIOUS VECTOR OPERATIONS ***
14C
15C COPIED FROM F14BLO.PAMLIB 23.4.82
16C AUTHOR: V.BLOBEL (UNIVERSITY OF HAMBURG)
17C DESYLIB
18C
19C *** BLANK COMMON REPLACED BY /VECUTY/ TO MATCH GEANT/GHEISHA CODE ***
20C *** NOTE THAT P(10,100) HAS BECOME PV(10,200) DUE TO THIS ***
21C
22C UN-USED ENTRIES REMOVED :
23C "PCOP" "PEXC" "PZER" "PWRT" "DOT4" "IMPU" "IMPULI" "ADD3"
24C "SUB3" "CROSS" "DOT" "SMUL" "NORZ" "PARPER" "PUNIT" "TRAP"
25C
26C *** ALL ENTRIES RE-WRITTEN AS SUBROUTINES USING ONLY NECESSARY ***
27C *** "DOUBLE PRECISION" STMTS. AND ALL SPECIFIC FUNCTIONS HAVE ***
28C *** BEEN CHANGED TO THEIR GENERIC EQUIVALENCES ***
29C *** NVE 29-MAR-1988 CERN GENEVA ***
30C
31C ORIGIN : H.FESEFELDT (22-JUNE-1984)
32C
33 SUBROUTINE LOR(K,L,M)
34C
35C --- PV-ARRAY ---
36#include "geant321/mxgkgh.inc"
37#include "geant321/s_blankp.inc"
38C
39 DOUBLE PRECISION A,B,C
40C
41 A=0.0
42 DO 6 I=1,3
43 A=A+PV(I,K)*PV(I,L)
44 6 CONTINUE
45 A=(A/(PV(4,L)+PV(5,L))-PV(4,K))/PV(5,L)
46 B=PV(5,K)*PV(5,K)
47 DO 8 I=1,3
48 C=PV(I,K)+A*PV(I,L)
49 B=B+C*C
50 PV(I,M)=C
51 8 CONTINUE
52 PV(4,M)=SQRT(B)
53 PV(5,M)=PV(5,K)
54 RETURN
55 END