Bugfix in AliPoints2Memory
[u/mrichter/AliRoot.git] / GEANT321 / erdecks / ermcsc.F
CommitLineData
fe4da5cc 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*
14C. ******************************************************************
15C. * *
16C. * Calculates Multiple Scattering Matrix in a step *
17C. * *
18C. * E - Error Matrix (Triangular) *
19C. * STEP - Step length in cm *
20C. * RADL - Radiation length in cm *
21C. * VECT(7) - Momentum in GeV *
22C. * CLA - Cosine of the deep angle *
23C. * *
24C. * ==>Called by : ERPROP *
25C. * Author E.Nagy ********* *
26C. * *
27C. ******************************************************************
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.
39CC CRLOG = 1. + 0.201*LOG(RI)
40CC 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