[u/mrichter/AliRoot.git] / GEANT321 / gheisha / pmltpc.F

*
* $Id$
*
* $Log$
* Revision 1.1.1.1  1995/10/24 10:21:00  cernlib
* Geant
*
*
#include "geant321/pilot.h"
*CMZ :  3.21/02 29/03/94  15.41.38  by  S.Giani
*-- Author :
      FUNCTION PMLTPC(NP,NM,NZ,N,B,C)
C
C *** NVE 03-MAR-1988 CERN GENEVA ***
C
C ORIGIN : H.FESEFELDT (14-SEP-1987)
C
#include "geant321/limits.inc"
      RLNNPF=0.
      IF(NP.LE.1) GOTO 2
      DO 1 I=2,NP
    1 RLNNPF=RLNNPF+LOG(I*1.)
    2 RLNNMF=0.
      IF(NM.LE.1) GOTO 4
      DO 3 I=2,NM
    3 RLNNMF=RLNNMF+LOG(I*1.)
    4 RLNNZF=0.
      IF(NZ.LE.1) GOTO 6
      DO 5 I=2,NZ
    5 RLNNZF=RLNNZF+LOG(I*1.)
    6 PMLTPC=-(NP-NM+NZ+B)**2/(2*(C*N)**2)-RLNNPF-RLNNMF-RLNNZF
      IF(PMLTPC.LT.EXPXL) PMLTPC=EXPXL
      PMLTPC=EXP(PMLTPC)
      RETURN
      END
Commit	Line	Data
fe4da5cc	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(CN)*2)-RLNNPF-RLNNMF-RLNNZF
	32	IF(PMLTPC.LT.EXPXL) PMLTPC=EXPXL
	33	PMLTPC=EXP(PMLTPC)
	34	RETURN
	35	END