GEANT321/gcons/gfcoul.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:20:13  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/02 29/03/94  15.41.19  by  S.Giani
  11 *-- Author :
  12       FUNCTION GFCOUL(Z)
  13 C.
  14 C.    ******************************************************************
  15 C.    *                                                                *
  16 C.    *  Compute Coulomb correction for pair production and Brem       *
  17 C.    *  REFERENCE : EGS MANUAL SLAC 210 - UC32 - JUNE 78              *
  18 C.    *                        FORMULA 2.7.17                          *
  19 C.    *                                                                *
  20 C.    *    ==>Called by : GSMIXT                                       *
  21 C.    *       Author    M.Maire *********                              *
  22 C.    *                                                                *
  23 C.    ******************************************************************
  24 C.
  25       PARAMETER (ALPHA = 7.29927E-03)
  26 C
  27       AZ2   = (ALPHA*Z)**2
  28       AZ4   =   AZ2 * AZ2
  29       FP    = ( 0.0083*AZ4 + 0.20206 + 1./(1.+AZ2) ) * AZ2
  30       FM    = ( 0.0020*AZ4 + 0.0369  ) * AZ4
  31       GFCOUL = FP - FM
  32       END