]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 1 | * |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
57b2d457 | 5 | * Revision 1.1.1.1 1999/05/18 15:55:20 fca |
6 | * AliRoot sources | |
7 | * | |
fe4da5cc | 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 | * | |
57b2d457 | 58 | CALL PYEEVT(IFLUND,ECLUND) |
fe4da5cc | 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 |