GEANT321/gcons/gsdk.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:20:16  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/03 02/08/94  18.13.20  by  S.Ravndal
  11 *-- Author :
  12       SUBROUTINE GSDK(IPART,BRATIO,MODE)
  13 C.
  14 C.    ******************************************************************
  15 C.    *                                                                *
  16 C.    *  Defines branching ratios and decay modes for standard         *
  17 C.    *  GEANT particles.                                              *
  18 C.    *                                                                *
  19 C.    *    ==>Called by : <USER>, GPART                                *
  20 C.    *       Author    G.Patrick  *********                           *
  21 C.    *                                                                *
  22 C.    ******************************************************************
  23 C.
  24 #include "geant321/gcbank.inc"
  25 #include "geant321/gcphys.inc"
  26 #include "geant321/gcnum.inc"
  27 #include "geant321/gcunit.inc"
  28       DIMENSION BRATIO(6)
  29       DIMENSION MODE(6)
  30 C.
  31 C.    ------------------------------------------------------------------
  32 C.
  33       IF (IPART.LE.0)                              GO TO 99
  34       IF (IPART.GT.NPART)                          GO TO 99
  35 C
  36 C     Particle pointer.
  37 C
  38       JPA = LQ(JPART-IPART)
  39       IF (JPA.LE.0)                                GO TO 99
  40 C
  41 C     Book decay links and data banks.
  42 C
  43       JDK1=LQ(JPA-1)
  44       JDK2=LQ(JPA-2)
  45       IF(JDK1+JDK2.NE.0) THEN
  46          WRITE(CHMAIL, 10000)
  47          CALL GMAIL(1,0)
  48          CALL GPPART(IPART)
  49          CALL MZDROP(IXCONS,LQ(JPA-1),' ')
  50          CALL MZDROP(IXCONS,LQ(JPA-2),' ')
  51       ENDIF
  52       CALL MZBOOK(IXCONS,JDK1,JPA,-1,'PABR',0,0,6,3,0)
  53       JPA=LQ(JPART-IPART)
  54       CALL MZBOOK(IXCONS,JDK2,JPA,-2,'PAMO',0,0,6,2,0)
  55       JPA=LQ(JPART-IPART)
  56       JDK1=LQ(JPA-1)
  57       IQ(JDK1-5)=IPART
  58       IQ(JDK2-5)=IPART
  59 C
  60 C     Store branching ratios & decay modes.
  61 C
  62       DO 20 I=1,6
  63          Q(JDK1+I) = BRATIO(I)
  64          IQ(JDK2+I) = MODE(I)
  65    20 CONTINUE
  66 C
  67   99  RETURN
  68 10000 FORMAT(' *** GSDK ***: Warning, redefinition of decay ',
  69      +       'for particle:')
  70       END