GEANT321/gphys/grangi1.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:21:33  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 #if defined(CERNLIB_NEVER)
  11 *CMZ :  3.21/02 29/03/94  15.41.23  by  S.Giani
  12 *-- Author :
  13       SUBROUTINE GRANGI
  14 C.
  15 C.    ******************************************************************
  16 C.    *                                                                *
  17 C.    *  Calculates the stopping range , in cm.                        *
  18 C.    *     (i.e. the maximum step due to the continuous energy loss)  *
  19 C.    *      The stopping range is the integral of the inverse of the  *
  20 C.    *      DE/DX table.                                              *
  21 C.    *     There are 4 tables : electron,positron,muon,proton         *
  22 C.    *                                                                *
  23 C.    *    ==>Called by : GPHYSI                                       *
  24 C.    *       Author      M.Maire   *********                          *
  25 C.    *                                                                *
  26 C.    ******************************************************************
  27 C.
  28 #include "geant321/gcbank.inc"
  29 #include "geant321/gctrak.inc"
  30 #include "geant321/gcjloc.inc"
  31 #include "geant321/gcmulo.inc"
  32 #include "geant321/gconsp.inc"
  33 #include "geant321/gcmate.inc"
  34 *
  35       DO 50 IEKBIN=1,NEK1
  36 *
  37 * *** Electrons
  38 *
  39          JLOSS = LQ(JMA-1)
  40          JRANG = LQ(JMA-15)
  41          IF (Z.LT.1.) THEN
  42             Q(JRANG+IEKBIN) = BIG
  43          ELSE IF ((IEKBIN.EQ.1).OR.(Q(JLOSS+IEKBIN-1).LE.0.)) THEN
  44             Q(JRANG+IEKBIN) = 0.
  45          ELSE
  46             DEBIN=ELOW(IEKBIN) - ELOW(IEKBIN-1)
  47             FRACT= Q(JLOSS+IEKBIN)/Q(JLOSS+IEKBIN-1) - 1.
  48             IF (ABS(FRACT).GE.0.1) THEN
  49                DRANG=DEBIN/(Q(JLOSS+IEKBIN)-Q(JLOSS+IEKBIN-1))
  50                DRANG=DRANG*LOG(1.+FRACT)
  51             ELSE
  52                DRANG=DEBIN/Q(JLOSS+IEKBIN-1)
  53                DRANG= DRANG*(1. - 0.5*FRACT + 0.3333*FRACT*FRACT)
  54             ENDIF
  55             Q(JRANG+IEKBIN)=Q(JRANG+IEKBIN-1) + DRANG
  56          ENDIF
  57 *
  58 * *** Positons
  59 *
  60          JLOSS = LQ(JMA-1) + NEK1
  61          JRANG = LQ(JMA-15) + NEK1
  62          IF (Z.LT.1.) THEN
  63             Q(JRANG+IEKBIN) = BIG
  64          ELSE IF ((IEKBIN.EQ.1).OR.(Q(JLOSS+IEKBIN-1).LE.0.)) THEN
  65             Q(JRANG+IEKBIN) = 0.
  66          ELSE
  67             DEBIN=ELOW(IEKBIN) - ELOW(IEKBIN-1)
  68             FRACT= Q(JLOSS+IEKBIN)/Q(JLOSS+IEKBIN-1) - 1.
  69             IF (ABS(FRACT).GE.0.1) THEN
  70                DRANG=DEBIN/(Q(JLOSS+IEKBIN)-Q(JLOSS+IEKBIN-1))
  71                DRANG=DRANG*LOG(1.+FRACT)
  72             ELSE
  73                DRANG=DEBIN/Q(JLOSS+IEKBIN-1)
  74                DRANG= DRANG*(1. - 0.5*FRACT + 0.3333*FRACT*FRACT)
  75             ENDIF
  76             Q(JRANG+IEKBIN)=Q(JRANG+IEKBIN-1) + DRANG
  77          ENDIF
  78 *
  79 * *** Muons
  80 *
  81          JLOSS = LQ(JMA-2)
  82          JRANG = LQ(JMA-16)
  83          IF (Z.LT.1.) THEN
  84             Q(JRANG+IEKBIN) = BIG
  85          ELSE IF ((IEKBIN.EQ.1).OR.(Q(JLOSS+IEKBIN-1).LE.0.)) THEN
  86             Q(JRANG+IEKBIN) = 0.
  87          ELSE
  88             DEBIN=ELOW(IEKBIN) - ELOW(IEKBIN-1)
  89             FRACT= Q(JLOSS+IEKBIN)/Q(JLOSS+IEKBIN-1) - 1.
  90             IF (ABS(FRACT).GE.0.1) THEN
  91                DRANG=DEBIN/(Q(JLOSS+IEKBIN)-Q(JLOSS+IEKBIN-1))
  92                DRANG=DRANG*LOG(1.+FRACT)
  93             ELSE
  94                DRANG=DEBIN/Q(JLOSS+IEKBIN-1)
  95                DRANG= DRANG*(1. - 0.5*FRACT + 0.3333*FRACT*FRACT)
  96             ENDIF
  97             Q(JRANG+IEKBIN)=Q(JRANG+IEKBIN-1) + DRANG
  98          ENDIF
  99 *
 100 * *** Protons
 101 *
 102          JLOSS = LQ(JMA-3)
 103          JRANG = LQ(JMA-16) + NEK1
 104          IF (Z.LT.1.) THEN
 105             Q(JRANG+IEKBIN) = BIG
 106          ELSE IF ((IEKBIN.EQ.1).OR.(Q(JLOSS+IEKBIN-1).LE.0.)) THEN
 107             Q(JRANG+IEKBIN) = 0.
 108          ELSE
 109             DEBIN=ELOW(IEKBIN) - ELOW(IEKBIN-1)
 110             FRACT= Q(JLOSS+IEKBIN)/Q(JLOSS+IEKBIN-1) - 1.
 111             IF (ABS(FRACT).GE.0.1) THEN
 112                DRANG=DEBIN/(Q(JLOSS+IEKBIN)-Q(JLOSS+IEKBIN-1))
 113                DRANG=DRANG*LOG(1.+FRACT)
 114             ELSE
 115                DRANG=DEBIN/Q(JLOSS+IEKBIN-1)
 116                DRANG= DRANG*(1. - 0.5*FRACT + 0.3333*FRACT*FRACT)
 117             ENDIF
 118             Q(JRANG+IEKBIN)=Q(JRANG+IEKBIN-1) + DRANG
 119          ENDIF
 120    50 CONTINUE
 121 *
 122       END
 123 #endif