GEANT321/gkine/gludky.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:21:18  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/04 09/01/95  16.20.45  by  S.Ravndal
  11 *-- Author :
  12       SUBROUTINE GLUDKY
  13 C.
  14 C.    ******************************************************************
  15 C.    *                                                                *
  16 C.    *        GEANT INTERFACE TO THE LUND MONTE-CARLO                 *
  17 C.    *        FOR TAU DECAY GENERATION ONLY                           *
  18 C.    *                                                                *
  19 C.    *  THIS INTERFACE TO LUND                                        *
  20 C.    *  A) INHIBITS DECAYS OF K0S,SIG+,SIG-,XI0,XI-,LAM,OM-           *
  21 C.    *  B) ACCEPTS STABLE PARTICLES ONLY                              *
  22 C.    *  C) CONVERTS THE PARTICLE CODES TO GEANT CONVENTIONS           *
  23 C.    *                                                                *
  24 C.    *  THUS THE PARTICLES FORWARDED TO GEANT ARE                     *
  25 C.    *  GAMMA,ELECTRON,NEUTRINO(E,MU,TAU),MU-,PI+,K+,K0S,K0L,         *
  26 C.    *  PROTON,NEUTRON,SIG(+,-),XI(0,-),LAM,OM-                       *
  27 C.    *  ....  AND THEIR ANTIPARTICLES                                 *
  28 C.    *                                                                *
  29 C.    *       That routine can be called from the user routine GUDCAY  *
  30 C.    *                                                                *
  31 C.    *       Author    R.Hemingway  26 SEPT 1985                      *
  32 C.    *                                                                *
  33 C.    ******************************************************************
  34 C.
  35 C
  36 #include "geant321/gcflag.inc"
  37 #include "geant321/gckine.inc"
  38 #include "geant321/gcking.inc"
  39 #include "geant321/gctrak.inc"
  40 #include "geant321/gconsp.inc"
  41 #include "geant321/lujets.inc"
  42 #include "geant321/ludat3.inc"
  43       PARAMETER (NSEL=44)
  44       DIMENSION IPLUND(NSEL)
  45       DATA IPLUND/
  46      +    22,   -11,    11,    12,   -13,    13,   111,   211,
  47      +  -211,   130,   321,  -321,  2112,  2212, -2212,   310,
  48      +   221,  3122,  3222,  3212,  3112,  3322,  3312,  3334,
  49      + -2112, -3122, -3112, -3212, -3222, -3322, -3312, -3334,
  50      +   -15,    15,   411,  -411,   421,  -421,   431,  -431,
  51      +  4122,    24,   -24,    23/
  52 *.
  53 *.    ------------------------------------------------------------------
  54 *.
  55       IF (IFINIT(6).EQ.0)THEN
  56 *
  57 *--  Set some particles to be stable in LUND so that GEANT can handle
  58 *--  their tracking and decays. (K0s, sigma+-, xi0-, lambda0, omega-)
  59 *
  60          MDCY(LUCOMP(310) ,1)=0
  61          MDCY(LUCOMP(3222),1)=0
  62          MDCY(LUCOMP(3112),1)=0
  63          MDCY(LUCOMP(3322),1)=0
  64          MDCY(LUCOMP(3312),1)=0
  65          MDCY(LUCOMP(3122),1)=0
  66          MDCY(LUCOMP(3334),1)=0
  67          IFINIT(6) = 1
  68       ENDIF
  69 C
  70       IF(IPART.EQ.33.OR.IPART.EQ.34) THEN
  71          N=0
  72          EDKY=GETOT
  73          THDKY=ACOS(VECT(6))
  74          PHDKY=ATAN2(VECT(5),VECT(4))
  75 *
  76 *            for TAU+ TAU-
  77 *            could be extended to all charm particles
  78 *
  79          CALL LU1ENT(1,IPLUND(IPART),EDKY,THDKY,PHDKY)
  80          CALL LUEXEC
  81          NGKINE=0
  82          DO 20 I = 1,N
  83             IF(K(I,1).GE.1.AND.K(I,1).LE.10) THEN
  84                DO 10 IGE=1, NSEL
  85                   IF(IPLUND(IGE).EQ.K(I,2)) THEN
  86                      NGKINE=NGKINE+1
  87                      GKIN(1,NGKINE)=P(I,1)
  88                      GKIN(2,NGKINE)=P(I,2)
  89                      GKIN(3,NGKINE)=P(I,3)
  90                      GKIN(4,NGKINE)=P(I,4)
  91                      GKIN(5,NGKINE)=IGE
  92                      GPOS(1,NGKINE)=V(I,1)
  93                      GPOS(2,NGKINE)=V(I,2)
  94                      GPOS(3,NGKINE)=V(I,3)
  95                      TOFD(NGKINE)  =V(I,4)/(10.*CLIGHT)
  96                   ENDIF
  97    10          CONTINUE
  98             ENDIF
  99    20    CONTINUE
 100       ENDIF
 101 *
 102       END