[u/mrichter/AliRoot.git] / GEANT321 / gcons / gsdk.F

*
* $Id$
*
* $Log$
* Revision 1.1.1.1  1995/10/24 10:20:16  cernlib
* Geant
*
*
#include "geant321/pilot.h"
*CMZ :  3.21/03 02/08/94  18.13.20  by  S.Ravndal
*-- Author :
      SUBROUTINE GSDK(IPART,BRATIO,MODE)
C.
C.    ******************************************************************
C.    *                                                                *
C.    *  Defines branching ratios and decay modes for standard         *
C.    *  GEANT particles.                                              *
C.    *                                                                *
C.    *    ==>Called by : <USER>, GPART                                *
C.    *       Author    G.Patrick  *********                           *
C.    *                                                                *
C.    ******************************************************************
C.
#include "geant321/gcbank.inc"
#include "geant321/gcphys.inc"
#include "geant321/gcnum.inc"
#include "geant321/gcunit.inc"
      DIMENSION BRATIO(6)
      DIMENSION MODE(6)
C.
C.    ------------------------------------------------------------------
C.
      IF (IPART.LE.0)                              GO TO 99
      IF (IPART.GT.NPART)                          GO TO 99
C
C     Particle pointer.
C
      JPA = LQ(JPART-IPART)
      IF (JPA.LE.0)                                GO TO 99
C
C     Book decay links and data banks.
C
      JDK1=LQ(JPA-1)
      JDK2=LQ(JPA-2)
      IF(JDK1+JDK2.NE.0) THEN
         WRITE(CHMAIL, 10000)
         CALL GMAIL(1,0)
         CALL GPPART(IPART)
         CALL MZDROP(IXCONS,LQ(JPA-1),' ')
         CALL MZDROP(IXCONS,LQ(JPA-2),' ')
      ENDIF
      CALL MZBOOK(IXCONS,JDK1,JPA,-1,'PABR',0,0,6,3,0)
      JPA=LQ(JPART-IPART)
      CALL MZBOOK(IXCONS,JDK2,JPA,-2,'PAMO',0,0,6,2,0)
      JPA=LQ(JPART-IPART)
      JDK1=LQ(JPA-1)
      IQ(JDK1-5)=IPART
      IQ(JDK2-5)=IPART
C
C     Store branching ratios & decay modes.
C
      DO 20 I=1,6
         Q(JDK1+I) = BRATIO(I)
         IQ(JDK2+I) = MODE(I)
   20 CONTINUE
C
  99  RETURN
10000 FORMAT(' *** GSDK ***: Warning, redefinition of decay ',
     +       'for particle:')
      END
Commit	Line	Data
fe4da5cc	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