]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 1 | * |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
af0a46b3 | 5 | * Revision 1.1.1.1 1999/05/18 15:55:21 fca |
6 | * AliRoot sources | |
7 | * | |
fe4da5cc | 8 | * Revision 1.1.1.1 1995/10/24 10:21:43 cernlib |
9 | * Geant | |
10 | * | |
11 | * | |
12 | #include "geant321/pilot.h" | |
13 | *CMZ : 3.21/02 29/03/94 15.41.24 by S.Giani | |
14 | *-- Author : | |
15 | SUBROUTINE GSXYZ | |
16 | C. | |
17 | C. ****************************************************************** | |
18 | C. * * | |
19 | C. * Store space point VECT in banks JXYZ * | |
20 | C. * * | |
21 | C. * ==>Called by : <USER> * | |
22 | C. * Authors R.Brun, F.Carena ********* * | |
23 | C. * * | |
24 | C. ****************************************************************** | |
25 | C. | |
26 | #include "geant321/gcbank.inc" | |
27 | #include "geant321/gcnum.inc" | |
28 | #include "geant321/gctrak.inc" | |
29 | #include "geant321/gcpush.inc" | |
30 | #include "geant321/gckine.inc" | |
31 | C. | |
32 | C. ------------------------------------------------------------------ | |
33 | C. | |
34 | IF(JXYZ.EQ.0)THEN | |
35 | C | |
36 | C Create mother bank JXYZ | |
37 | C | |
38 | CALL MZBOOK(IXDIV,JXYZ,JXYZ,1,'JXYZ',NCKINE,NCKINE,0,3,0) | |
39 | IQ(JXYZ-5)=0 | |
40 | ENDIF | |
41 | C | |
42 | NXYZ=IQ(JXYZ-2) | |
43 | IF(ITRA.GT.NXYZ)THEN | |
44 | NPUSH=MAX(10,ITRA-NXYZ) | |
45 | CALL MZPUSH(IXDIV,JXYZ,NPUSH,0,'I') | |
46 | ENDIF | |
47 | C | |
48 | C Create bank to store points for primary track ITRA | |
49 | C | |
50 | JX=LQ(JXYZ-ITRA) | |
51 | IF(JX.EQ.0)THEN | |
52 | CALL MZBOOK(IXDIV,JX,JXYZ,-ITRA,'JXYZ',0,0,NCJXYZ,3,0) | |
53 | Q(JX+1)=NCJXYZ-2. | |
54 | Q(JX+2)=3. | |
55 | IQ(JX)=NTMULT | |
56 | C | |
57 | C Check if there is enough space to store current point | |
58 | C | |
59 | ELSE | |
60 | NFREE=Q(JX+1) | |
61 | IF(NFREE.LE.6)THEN | |
62 | CALL MZPUSH(IXDIV,JX,0,NPJXYZ,'I') | |
63 | JX=LQ(JXYZ-ITRA) | |
64 | Q(JX+1)=Q(JX+1)+NPJXYZ | |
65 | ENDIF | |
66 | C | |
67 | IF(NTMULT.NE.IQ(JX))THEN | |
68 | IQ(JX)=NTMULT | |
69 | IPOINT=JX+Q(JX+2) | |
70 | IF(Q(IPOINT).NE.0.) THEN | |
71 | Q(JX+2)=Q(JX+2)+3.*Q(IPOINT)+2. | |
72 | Q(JX+1)=Q(JX+1)-2. | |
af0a46b3 | 73 | IPOINT=JX+Q(JX+2) |
74 | Q(IPOINT)=0 | |
fe4da5cc | 75 | ENDIF |
76 | ENDIF | |
77 | ENDIF | |
78 | C | |
79 | IPOINT=JX+Q(JX+2) | |
80 | JP=IPOINT+3.*Q(IPOINT)+1. | |
81 | Q(IPOINT)=Q(IPOINT)+1. | |
82 | Q(IPOINT+1)=IPART | |
83 | DO 10 I=1,3 | |
84 | 10 Q(JP+I)=VECT(I) | |
85 | Q(JX+1)=Q(JX+1)-3. | |
86 | C | |
87 | END |