| fe4da5cc |
1 | * |
| 2 | * $Id$ |
| 3 | * |
| 4 | * $Log$ |
| 5 | * Revision 1.1.1.1 1995/10/24 10:20:01 cernlib |
| 6 | * Geant |
| 7 | * |
| 8 | * |
| 9 | #include "geant321/pilot.h" |
| 10 | *CMZ : 3.21/02 29/03/94 15.41.44 by S.Giani |
| 11 | *-- Author : |
| 12 | *$ CREATE KPOIS.FOR |
| 13 | *COPY KPOIS |
| 14 | * |
| 15 | *=== kpois ============================================================* |
| 16 | * |
| 17 | FUNCTION KPOIS(EXPLAM) |
| 18 | |
| 19 | #include "geant321/dblprc.inc" |
| 20 | #include "geant321/dimpar.inc" |
| 21 | #include "geant321/iounit.inc" |
| 22 | REAL RNDM(1) |
| 23 | C*** EXPLAM IS EXP(-LAMBDA) |
| 24 | C*** SQMPLING FROM POISSON DISTRIBUTION, AVERAGE LAMBDA |
| 25 | C*** PK = EXP(-LAMBDA)*LAMBDA**K/(K-FACTORIAL) |
| 26 | Y=1.D0 |
| 27 | KPOISS=0 |
| 28 | 5 CONTINUE |
| 29 | CALL GRNDM(RNDM,1) |
| 30 | Y=Y*RNDM(1) |
| 31 | IF (Y.LT.EXPLAM) GO TO 10 |
| 32 | KPOISS=KPOISS+1 |
| 33 | GO TO 5 |
| 34 | 10 CONTINUE |
| 35 | KPOIS=KPOISS |
| 36 | RETURN |
| 37 | END |
git://git.uio.no / u / mrichter / AliRoot.git / blame