GEANT321/gkine/glundi.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1999/05/18 15:55:20  fca
   6 * AliRoot sources
   7 *
   8 * Revision 1.1.1.1  1995/10/24 10:21:18  cernlib
   9 * Geant
  10 *
  11 *
  12 #include "geant321/pilot.h"
  13 *CMZ :  3.21/02 29/03/94  15.41.21  by  S.Giani
  14 *-- Author :
  15       SUBROUTINE GLUNDI
  16 C.
  17 C.    ******************************************************************
  18 C.    *                                                                *
  19 C.    *        GEANT INTERFACE TO THE LUND MONTE-CARLO                 *
  20 C.    *            (Initialization routine)                            *
  21 C.    *                                                                *
  22 C.    *  FOR TECHNICAL DETAILS OF THE LUND MONTE CARLO,SEE             *
  23 C.    *  1. THE LUND MONTE CARLO FOR JET FRAGMENTATION - T.SJOSTRAND   *
  24 C.    *     LU TP 82-3  and COMPUTER PHYSICS COMMUNICATION 27(1982)243 *
  25 C.    *  2. THE LUND MONTE CARLO FOR E+E- JET PHYSICS  - T.SJOSTRAND   *
  26 C.    *     LU TP 82-7 and COMPUTER PHYSICS COMMUNICATION 28(1983)229  *
  27 C.    *                                                                *
  28 C.    *  THE PARTICLES FORWARDED TO GEANT ARE                          *
  29 C.    *  GAMMA,ELECTRON,NEUTRINO(E,MU,TAU),MU-,PI+,K+,K0S,K0L,         *
  30 C.    *  PROTON,NEUTRON,SIG(+,-),XI(0,-),LAM,OM-                       *
  31 C.    *  ....  AND THEIR ANTIPARTICLES                                 *
  32 C.    *                                                                *
  33 C.    *                                                                *
  34 C.    *  THE USER MUST SUPPLY TWO PARAMETERS FOR EVENT GENERATION      *
  35 C.    *  VIA     CALL LUEEVT(IFLUND,ECLUND)     WHERE                  *
  36 C.    *  IFLUND=THE FLAVOUR CODE (1=U,2=D,3=S,4=C,5=B,6=T,0=MIX)       *
  37 C.    *  ECLUND=THE TOTAL CMS ENERGY                                   *
  38 C.    *                                                                *
  39 C.    *     GLUNDI is automatically called by GLUND.                   *
  40 C.    *     GLUNDI can be called by the user in the initialization     *
  41 C.    *     routine UGINIT just before reading data cards (GFFGO)      *
  42 C.    *     In this case the default parameters can be changed by      *
  43 C.    *     the data cards KTYP,PMAS,PWID,MSTE,LUND and IDB            *
  44 C.    *                                                                *
  45 C.    *    ==>Called by : <USER>, GUKINE                               *
  46 C.    *       Author    R.Hemingway , D.Ward *********                 *
  47 C.    *       modified by R.Brun 18/11/85 ******                       *
  48 C.    *                                                                *
  49 C.    ******************************************************************
  50 C.
  51 #include "geant321/gcflag.inc"
  52 #include "geant321/ludat1.inc"
  53 #include "geant321/ludat3.inc"
  54 #include "geant321/gclund.inc"
  55 C.
  56 C.    ------------------------------------------------------------------
  57 C.
  58       IF (IFINIT(3).EQ.0)THEN
  59 *
  60 *              Set defaults
  61 *
  62          IFINIT(3) = 1
  63          IFLUND = 0
  64          ECLUND = 92.25
  65 *
  66 *--  Initial state radiative effects     MSTJ(107) =  1
  67 *
  68          MSTJ(107)= 1
  69 *
  70 *--  Set some particles to be stable in LUND so that GEANT can handle
  71 *--  their tracking and decays. (K0s, sigma+-, xi0-, lambda0, omega-)
  72 *
  73          MDCY(PYCOMP(310) ,1)=0
  74          MDCY(PYCOMP(3222),1)=0
  75          MDCY(PYCOMP(3112),1)=0
  76          MDCY(PYCOMP(3322),1)=0
  77          MDCY(PYCOMP(3312),1)=0
  78          MDCY(PYCOMP(3122),1)=0
  79          MDCY(PYCOMP(3334),1)=0
  80 *
  81          CALL FFKEY('LUND',IFLUND,2,'MIXED')
  82       ENDIF
  83 C
  84       END