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