GEANT321/gheisha/pmltpc.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:21:00  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/02 29/03/94  15.41.38  by  S.Giani
  11 *-- Author :
  12       FUNCTION PMLTPC(NP,NM,NZ,N,B,C)
  13 C
  14 C *** NVE 03-MAR-1988 CERN GENEVA ***
  15 C
  16 C ORIGIN : H.FESEFELDT (14-SEP-1987)
  17 C
  18 #include "geant321/limits.inc"
  19       RLNNPF=0.
  20       IF(NP.LE.1) GOTO 2
  21       DO 1 I=2,NP
  22     1 RLNNPF=RLNNPF+LOG(I*1.)
  23     2 RLNNMF=0.
  24       IF(NM.LE.1) GOTO 4
  25       DO 3 I=2,NM
  26     3 RLNNMF=RLNNMF+LOG(I*1.)
  27     4 RLNNZF=0.
  28       IF(NZ.LE.1) GOTO 6
  29       DO 5 I=2,NZ
  30     5 RLNNZF=RLNNZF+LOG(I*1.)
  31     6 PMLTPC=-(NP-NM+NZ+B)**2/(2*(C*N)**2)-RLNNPF-RLNNMF-RLNNZF
  32       IF(PMLTPC.LT.EXPXL) PMLTPC=EXPXL
  33       PMLTPC=EXP(PMLTPC)
  34       RETURN
  35       END