This commit was generated by cvs2svn to compensate for changes in r2,
[u/mrichter/AliRoot.git] / GEANT321 / erdecks / ermcsc.F
1 *
2 * $Id$
3 *
4 * $Log$
5 * Revision 1.1.1.1  1996/03/06 15:37:35  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 ERMCSC (E)
13 *
14 C.    ******************************************************************
15 C.    *                                                                *
16 C.    *       Calculates Multiple Scattering Matrix in a step          *
17 C.    *                                                                *
18 C.    *           E        - Error Matrix (Triangular)                 *
19 C.    *           STEP     - Step length in cm                         *
20 C.    *           RADL     - Radiation length in cm                    *
21 C.    *           VECT(7)  - Momentum in GeV                           *
22 C.    *           CLA      - Cosine of the deep angle                  *
23 C.    *                                                                *
24 C.    *    ==>Called by : ERPROP                                       *
25 C.    *       Author    E.Nagy  *********                              *
26 C.    *                                                                *
27 C.    ******************************************************************
28 *
29 #include "geant321/gcmate.inc"
30 #include "geant321/gckine.inc"
31 #include "geant321/gctrak.inc"
32 *
33       DOUBLE PRECISION E(15)
34 *
35       PBETA = VECT(7)*VECT(7)/GETOT
36       RI    = STEP/RADL
37 *
38       CRLOG = 1.
39 CC      CRLOG = 1. + 0.201*LOG(RI)
40 CC      IF (CRLOG.LE.0.) CRLOG = 0.
41       D     = 1.8496E-4*RI*((CHARGE*CRLOG/PBETA)**2)
42       S1    = D*STEP**2/3.
43       S2    = D
44       S3    = D*STEP/2.
45       CLA = SQRT (VECT(4)**2+VECT(5)**2)
46       E(6)  = E(6)+S2
47       E(9)  = E(9)-S3
48       E(10) = E(10)+S2/CLA**2
49       E(11) = E(11)+S3/CLA
50       E(13) = E(13)+S1
51       E(15) = E(15)+S1
52 *
53       END