]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 1 | * |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
5 | * Revision 1.1.1.1 1995/10/24 10:20:03 cernlib | |
6 | * Geant | |
7 | * | |
8 | * | |
9 | #include "geant321/pilot.h" | |
10 | #if defined(CERNLIB_OLDNAME) | |
11 | *CMZ : 3.21/02 29/03/94 15.41.44 by S.Giani | |
12 | *-- Author : | |
13 | *=== vertex ===========================================================* | |
14 | * | |
15 | SUBROUTINE VERTEX(IT,LT,LL,KFA,E0,IV,RE,KFR1,KFR2,AME,IOPT) | |
16 | ||
17 | #include "geant321/dblprc.inc" | |
18 | #include "geant321/dimpar.inc" | |
19 | #include "geant321/iounit.inc" | |
20 | REAL RNDM(2) | |
21 | C CHOICE OF THE VERTEX | |
22 | DIMENSION RE(*),KFR1(*),KFR2(*),IV(*) | |
23 | C WAEHLT PASSENDEN VERTEX IV AUS | |
24 | IF(LT.EQ.1) WRITE(LUNOUT,288)IT,LT,LL,KFA,E0,AME,IOPT | |
25 | 288 FORMAT (4I5,2E12.4,I5,' VERTEX IT,LT,LL,KFA,E0,AME') | |
26 | PDB=0.8D0 | |
27 | IVA=1 | |
28 | I=IT | |
29 | J=IT-1 | |
30 | IVX=IVA | |
31 | IF(IOPT.EQ.2) IVX=2 | |
32 | IF(IOPT.EQ.4.AND.KFA.GT.6.AND.LL.EQ.0) IVX=2 | |
33 | IF(IOPT.EQ.4.AND.KFA.LE.6.AND.LL.EQ.1) IVX=2 | |
34 | IF(IOPT.EQ.5) IVX=2 | |
35 | IF(I.GT.1)IVX=IV(J) | |
36 | IF(LL.EQ.1.AND.I.GT.1)IVX=IV(J)-5 | |
37 | CALL GRNDM(RNDM,2) | |
38 | Z=RNDM(1) | |
39 | Y=RNDM(2) | |
40 | GO TO (1,2,1,2,2),IVX | |
41 | 1 IVX=1 | |
42 | IF(Z.GE.AME) IVX=2 | |
43 | GO TO 20 | |
44 | 2 IVX=3 | |
45 | IF(Z.GT.PDB) IVX=4 | |
46 | IF(Z.GT.PDB.AND.Y.LE.0.5D0) IVX=5 | |
47 | GO TO 20 | |
48 | 20 CONTINUE | |
49 | IF(LL.EQ.1)IVX=IVX+5 | |
50 | IV(I)=IVX | |
51 | IF(LT.EQ.0)GO TO 60 | |
52 | C AUSDRUCKEN | |
53 | IF(I.EQ.1)GO TO 30 | |
54 | WRITE(LUNOUT,50)IV(J),IV(I),LL | |
55 | GO TO 40 | |
56 | 30 WRITE(LUNOUT,50)IVA,IV(I),LL | |
57 | 40 CONTINUE | |
58 | 50 FORMAT(1H0,15HIV(J),IV(I),LL=,3I5) | |
59 | 60 CONTINUE | |
60 | RETURN | |
61 | END | |
62 | #endif |