GEANT321/fluka/vertex.F

   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