GEANT321/gkine/glund.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 GLUND
  16 C.
  17 C.    ******************************************************************
  18 C.    *                                                                *
  19 C.    *        GEANT INTERFACE TO THE LUND MONTE-CARLO                 *
  20 C.    *                                                                *
  21 C.    *  FOR TECHNICAL DETAILS OF THE LUND MONTE CARLO,SEE             *
  22 C.    *  1. THE LUND MONTE CARLO FOR JET FRAGMENTATION - T.SJOSTRAND   *
  23 C.    *     LU TP 82-3  and COMPUTER PHYSICS COMMUNICATION 27(1982)243 *
  24 C.    *  2. THE LUND MONTE CARLO FOR E+E- JET PHYSICS  - T.SJOSTRAND   *
  25 C.    *     LU TP 82-7 and COMPUTER PHYSICS COMMUNICATION 28(1983)229  *
  26 C.    *                                                                *
  27 C.    *    The parameters to use the LUND interface are set in GLUNDI  *
  28 C.    *    They can be changed by data cards (see GLUNDI)              *
  29 C.    *                                                                *
  30 C.    *    ==>Called by : <USER>, GUKINE                               *
  31 C.    *       Author    R.Hemingway  *********                         *
  32 C.    *       modified by R.Brun  18/11/85 **                          *
  33 C.    *                                                                *
  34 C.    ******************************************************************
  35 C.
  36       DIMENSION VERT(3)
  37 #include "geant321/gcbank.inc"
  38 #include "geant321/gcflag.inc"
  39 #include "geant321/gcpush.inc"
  40 #include "geant321/lujets.inc"
  41 #include "geant321/gclund.inc"
  42       PARAMETER (NSEL=44)
  43       DIMENSION IPLUND(NSEL)
  44       REAL      PGEANT(3),BUF(1)
  45       EQUIVALENCE (BUF,K)
  46       DATA IPLUND/
  47      +    22,   -11,    11,    12,   -13,    13,   111,   211,
  48      +  -211,   130,   321,  -321,  2112,  2212, -2212,   310,
  49      +   221,  3122,  3222,  3212,  3112,  3322,  3312,  3334,
  50      + -2112, -3122, -3112, -3212, -3222, -3322, -3312, -3334,
  51      +   -15,    15,   411,  -411,   421,  -421,   431,  -431,
  52      +  4122,    24,   -24,    23/
  53 *.
  54 *.    ------------------------------------------------------------------
  55 *.
  56       IF (IFINIT(3).EQ.0)  CALL GLUNDI
  57 *
  58       CALL PYEEVT(IFLUND,ECLUND)
  59 *
  60       NCVERT = 10
  61       NCKINE = N+10
  62       NCJXYZ = NCKINE
  63       NPVERT =  5
  64       NPKINE = 10
  65       NPJXYZ = 10
  66 *
  67       VERT(1)   = 0.
  68       VERT(2)   = 0.
  69       VERT(3)   = 0.
  70       NWBUF=2+10*N
  71       CALL GSVERT(V,0,0,BUF,NWBUF,NVTX)
  72       IF (NVTX.GT.0) THEN
  73          DO 20 I = 1,N
  74             IF(K(I,1).GE.1.AND.K(I,1).LE.10) THEN
  75                DO 10 IGE=1, NSEL
  76                   IF(IPLUND(IGE).EQ.K(I,2)) THEN
  77                      PGEANT(1)=P(I,1)
  78                      PGEANT(2)=P(I,2)
  79                      PGEANT(3)=P(I,3)
  80                      CALL GSKINE(PGEANT,IGE,NVTX,BUF,0,NT)
  81                      GOTO 20
  82                   ENDIF
  83    10          CONTINUE
  84             ENDIF
  85    20    CONTINUE
  86 *
  87          JV = LQ(JVERTX-1)
  88          JUV = LQ(JV-1)
  89          Q(JUV+1)= 527.0
  90          Q(JUV+2)= N
  91 *
  92 *     Copy the LUND record into a bank of real variables.
  93 *
  94          DO 30 NP=1, N
  95             Q(JUV+2    +NP) = K(NP,1)
  96             Q(JUV+2+  N+NP) = K(NP,2)
  97             Q(JUV+2+2*N+NP) = K(NP,3)
  98             Q(JUV+2+3*N+NP) = K(NP,4)
  99             Q(JUV+2+4*N+NP) = K(NP,5)
 100             Q(JUV+2+5*N+NP) = P(NP,1)
 101             Q(JUV+2+6*N+NP) = P(NP,2)
 102             Q(JUV+2+7*N+NP) = P(NP,3)
 103             Q(JUV+2+8*N+NP) = P(NP,4)
 104             Q(JUV+2+9*N+NP) = P(NP,5)
 105    30    CONTINUE
 106 *
 107       ENDIF
 108       END