Bugfix in AliPoints2Memory
[u/mrichter/AliRoot.git] / GEANT321 / erpremc / trptsc.F
CommitLineData
fe4da5cc 1*
2* $Id$
3*
4* $Log$
5* Revision 1.1.1.1 1996/03/06 15:37:36 mclareni
6* Add geane321 source directories
7*
8*
9#include "geant321/pilot.h"
10*CMZ : 3.21/02 29/03/94 15.41.49 by S.Giani
11*-- Author :
12 SUBROUTINE TRPTSC(PC,RC,PD,RD,IERR)
13C
14C *** TRANSFORMS ERROR MATRIX
15C FROM SC VARIABLES (1/Pt,LAMBDA,PHI,YT,ZT)
16C FROM SC VARIABLES (1/P,LAMBDA,PHI,YT,ZT)
17
18
19#if !defined(CERNLIB_SINGLE)
20 IMPLICIT DOUBLE PRECISION (A-H,O-Z)
21 REAL PC,PD,RC,RD
22#endif
23#include "geant321/trcom3.inc"
24 DIMENSION PC(3),PD(3),RC(15),RD(15)
25*
26*______________________________________________________________________
27*
28 IERR = 0
29 COSL = COS(PC(2))
30 IF (ABS(COSL).EQ.0) GO TO 901
31 SINL = SIN(PC(2))
32*
33 PD(1) = PC(1)*COSL
34 PD(2) = PC(2)
35 PD(3) = PC(3)
36*
37 J=0
38*
39 DO 10 I=1,5
40 DO 5 K=I,5
41 J=J+1
42 A(I,K)=0.
43 A(K,I)=0.
44 S(J)=RC(J)
45 5 CONTINUE
46 10 CONTINUE
47*
48 A(1,1) = COSL
49 A(2,2) = 1.D0
50 A(3,3) = 1.D0
51 A(4,4) = 1.D0
52 A(5,5) = 1.D0
53*
54 A(1,2) = -PC(1)*SINL
55C
56 CALL SSMT5T(A,S,S)
57C
58 DO 25 J=1,15
59 RD(J)=S(J)
60 25 CONTINUE
61C
62 RETURN
63C
64C *** ERROR EXITS
65C
66 901 IERR=1
67 910 CONTINUE
68C
69 END