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