| fe4da5cc |
1 | * |
| 2 | * $Id$ |
| 3 | * |
| 4 | * $Log$ |
| 5 | * Revision 1.1.1.1 1995/10/24 10:19:59 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 AMGA.FOR |
| 13 | *COPY AMGA |
| 14 | * |
| 15 | *=== amga =============================================================* |
| 16 | * |
| 17 | FUNCTION AMGA(IT) |
| 18 | |
| 19 | #include "geant321/dblprc.inc" |
| 20 | #include "geant321/dimpar.inc" |
| 21 | #include "geant321/iounit.inc" |
| 22 | C*** RANDOM SELECTION OF MASSES OF DECAYING PARTICLES |
| 23 | COMMON/FKABLT/AMM(110),GA(110),TAU(110),ICH(110) |
| 24 | *,IBAR(110),K1(110),K2(110) |
| 25 | COMMON /FKGAMR/ REDU,AMO,AM (15 ) |
| 26 | DIMENSION GASUNI(14) |
| 27 | REAL RNDM(1) |
| 28 | SAVE GASUNI, GAUNO, GAUNON, IO, NSTAB |
| 29 | DATA GASUNI/ |
| 30 | *-1.D0,-.98D0,-.95D0,-.87D0,-.72D0,-.48D0,-.17D0,.17D0,.48D0, |
| 31 | *.72D0,.87D0,.95D0,.98D0,1.D0/ |
| 32 | DATA GAUNO/2.352D0/ |
| 33 | DATA GAUNON/2.4D0/ |
| 34 | DATA IO/14/ |
| 35 | DATA NSTAB/23/ |
| 36 | I=1 |
| 37 | IF (IT.LE.0) GO TO 3 |
| 38 | IF (IT.LE.NSTAB) GO TO 2 |
| 39 | DGAUNI=GAUNO*GAUNON/(IO-1.D0) |
| 40 | CALL GRNDM(RNDM,1) |
| 41 | VV = RNDM(1) |
| 42 | VV = VV*2.D0-1.D0+1.D-16 |
| 43 | 1 CONTINUE |
| 44 | VO=GASUNI(I) |
| 45 | I=I+1 |
| 46 | V1=GASUNI(I) |
| 47 | IF (VV.GT.V1) GO TO 1 |
| 48 | UNIGA=DGAUNI*(I-2.D0+(VV-VO+1.D-16)/(V1-VO)-(IO-1.D0)*.5D0) |
| 49 | DAM=GA(IT)*UNIGA/GAUNO |
| 50 | AAM=AMM(IT)+DAM |
| 51 | AMGA=AAM |
| 52 | RETURN |
| 53 | 2 CONTINUE |
| 54 | AMGA=AMM(IT) |
| 55 | RETURN |
| 56 | 3 CONTINUE |
| 57 | AMGA=0.D0 |
| 58 | RETURN |
| 59 | END |
git://git.uio.no / u / mrichter / AliRoot.git / blame