gpamela/gptof/gptofv.F

*
* $Id: gptofv.F,v 3.6 2006/11/16 18:45:29 pam-ba Exp $
*
* $Log: gptofv.F,v $
* Revision 3.6  2006/11/16 18:45:29  pam-ba
* Simulated an aluminum container for S4
*
* Revision 3.5  2006/05/18 10:52:32  pam-ba
* TOF geometry completed and a new material, the polystyrene (density 35 g/l), added
*
* Revision 3.4  2006/05/02 10:40:13  pam-ba
* TOF geometry and position updated and a new material, the mylar, added
*
* Revision 3.3  2006/04/09 23:29:22  cafagna
* Several new things, among this: ND and CARD
*
* Revision 3.2  2005/06/21 02:42:29  cafagna
* Major modification to the geometry and to the random number chain
*
* Revision 3.1.1.1  2002/07/11 16:02:01  cafagna
* First GPAMELA release on CVS
*
*
*CMZ :  2.01/00 05/04/2000  14.35.18  by  Marialuigia Ambriola
*CMZ :  2.00/00 03/03/2000  15.39.05  by  Francesco Cafagna
*CMZ :  1.02/00 21/10/99  16.30.25  by  Francesco Cafagna
*CMZ :  1.00/02 18/03/96  17.18.03  by  Francesco Cafagna
*-- Author :    Francesco Cafagna   05/12/95
      SUBROUTINE GPTOFV
************************************************************************
*                                                                      *
* Volumes definition for the scintillators counters                    *
* Called by: GPGEO                                                     *
* Author: Francesco Cafagna, 05/12/95 13.39.41                         *
*                                                                      *
************************************************************************
#include "gpgeo.inc"
#include "gpmed.inc"
*
      INTEGER IROT,IVOLU,N,NMED,NUM
      REAL X,Y,Z
*
      NMED=MAL
*
*  Define the S1     volume
*
      CALL GSVOLU('S1  ','BOX ',NMED,S1, 3,IVOLU)
*
*  Define the S2     volume
*
      CALL GSVOLU('S2  ','BOX ',NMED,S2, 3,IVOLU)
*
      NMED=MN2
*
*  Define the S1A     volume
*
      CALL GSVOLU('S1A ','BOX ',NMED,S1A, 3,IVOLU)
*
*  Define the S2A     volume
*
      CALL GSVOLU('S2A ','BOX ',NMED,S2A, 3,IVOLU)
*
*  Define the S3     volume
*
      CALL GSVOLU('S3  ','BOX ',NMED,S3, 3,IVOLU)
*
      NMED=MMYL
*
*  Define the S11M     volume
*
      CALL GSVOLU('S11M','BOX ',NMED,S11M, 3,IVOLU)
*
*  Define the S12M     volume
*
      CALL GSVOLU('S12M','BOX ',NMED,S12M, 3,IVOLU)
*
*  Define the S21M     volume
*
      CALL GSVOLU('S21M','BOX ',NMED,S21M, 3,IVOLU)
*
*  Define the S22M     volume
*
      CALL GSVOLU('S22M','BOX ',NMED,S22M, 3,IVOLU)
*
*  Define the S31M     volume
*
      CALL GSVOLU('S31M','BOX ',NMED,S31M, 3,IVOLU)
*
*  Define the S32M     volume
*
      CALL GSVOLU('S32M','BOX ',NMED,S32M, 3,IVOLU)
*
*  Define the POL1     volume
*
      NMED=MPOL
      CALL GSVOLU('POL1','BOX ',NMED,POL1, 3,IVOLU)
*
*  Define the POL2     volume
*
      CALL GSVOLU('POL2','BOX ',NMED,POL2, 3,IVOLU)
*
*  Define the POLY     volume
*
      CALL GSVOLU('POLY','BOX ',NMED,POLY, 3,IVOLU)
*
*
*  S4 container volume.
*
      NMED=MAL
      CALL GSVOLU('S4AL','BOX ',NMED,S4AL, 3,IVOLU)
*
*
*  Change medium just for S4.
*
      NMED=MSCIN
*
*  Define the S4     volume
*
      CALL GSVOLU('S4  ','BOX ',NMED,S4, 3,IVOLU)
*
*  Start subvolume definition.
*
*
*  Define the S11Y     volume
*
      CALL GSVOLU('S11Y','BOX ',NMED,S11Y, 3,IVOLU)
*
*  Define the S12X     volume
*
      CALL GSVOLU('S12X','BOX ',NMED,S12X, 3,IVOLU)
*
*  Define the S21X     volume
*
      CALL GSVOLU('S21X','BOX ',NMED,S21X, 3,IVOLU)
*
*  Define the S22Y     volume
*
      CALL GSVOLU('S22Y','BOX ',NMED,S22Y, 3,IVOLU)
*
*  Define the S31Y     volume
*
      CALL GSVOLU('S31Y','BOX ',NMED,S31Y, 3,IVOLU)
*
*  Define the S32X     volume
*
      CALL GSVOLU('S32X','BOX ',NMED,S32X, 3,IVOLU)

c$$$<<<<<<< gptofv.F
c$$$=======
c$$$c ml: 21/12/04:
c$$$      NMED=MAL
c$$$*
c$$$*  Define the SC11     volume
c$$$*
c$$$      CALL GSVOLU('SC11','BOX ',NMED,SC11, 3,IVOLU)
c$$$*
c$$$*  Define the SC12     volume
c$$$*
c$$$      CALL GSVOLU('SC12','BOX ',NMED,SC12, 3,IVOLU)
c$$$*
c$$$*  Define the SC21     volume
c$$$*
c$$$      CALL GSVOLU('SC21','BOX ',NMED,SC21, 3,IVOLU)
c$$$*
c$$$*  Define the SC22     volume
c$$$*
c$$$      CALL GSVOLU('SC22','BOX ',NMED,SC22, 3,IVOLU)
c$$$
c$$$>>>>>>> 3.3
*
* Positioning volumes
*
      X=0.
      Y=0.
      Z=0.
      CALL GSPOS('S1A ',N,'S1  ',X,Y,Z,0,'ONLY')
*
      N=1
      X=0.
      Y=0.
      Z=-S1A(3)+HGLUTOF+HPVCMIN/2.
      CALL GSPOS('POL2',N,'S1A ',X,Y,Z,0,'ONLY')
*
      X=0.
      Z=-S1A(3)+HGLUTOF+HPVCMIN+S12M(3)
      DO I=1,NS12X
         N= I
c ab mod         Y= -S1A(2) + HMYLTOF + (2*I-1)*S12M(2)
         Y= -S1A(2) + (2*I-1)*S12M(2)
         CALL GSPOS('S12M',N,'S1A ',X,Y,Z,0,'ONLY')
      ENDDO
*
      Y= 0.
      Z=-S1A(3)+HGLUTOF+HPVCMIN+2*S12M(3)+S11M(3)
      DO I=1,NS11Y
         N= I
c ab mod         X= -S1A(1) + HMYLTOF + (2*I-1)*S11M(1)
         X= -S1A(1) + (2*I-1)*S11M(1)
         CALL GSPOS('S11M',N,'S1A ',X,Y,Z,0,'ONLY')
      ENDDO
*
      N=1
      X=0.
      Y=0.
      Z=-S1A(3)+HGLUTOF+HPVCMIN+2*S12M(3)+2.*S11M(3)+HPVCMAX/2.
      CALL GSPOS('POL1',N,'S1A ',X,Y,Z,0,'ONLY')
C      Z=-S1A(3)+HGLUTOF+HPVCMAX+S12M(3)
C      CALL GSPOS('S12M',N,'S1A ',X,Y,Z,0,'ONLY')
C      Z=-S1A(3)+HGLUTOF+0.4+2*S12M(3)+S11M(3)
C      CALL GSPOS('S11M',N,'S1A ',X,Y,Z,0,'ONLY')
*
      X=0.
      Y=0.
      Z=0.
      CALL GSPOS('S2A ',N,'S2  ',X,Y,Z,0,'ONLY')
*
      Y= 0.
      Z=-S2A(3)+HGLUTOF+S22M(3)
      DO I=1,NS22Y
         N= I
c ab mod         X= -S2A(1) + HMYLTOF + (2*I-1)*S22M(1)
         X= -S2A(1) + (2*I-1)*S22M(1)
         CALL GSPOS('S22M',N,'S2A ',X,Y,Z,0,'ONLY')
      ENDDO
*
      X= 0.
      Z=-S2A(3)+HGLUTOF+2*S22M(3)+S21M(3)
      DO I=1,NS21X
         N= I
c ab mod         Y= -S2A(2) + HMYLTOF+ (2*I-1)*S21M(2)
         Y= -S2A(2) + (2*I-1)*S21M(2)
         CALL GSPOS('S21M',N,'S2A ',X,Y,Z,0,'ONLY')
      ENDDO
*
      N=1
      X=0.
      Y=0.
      Z=-S2A(3)+HGLUTOF+2*S22M(3)+2.*S21M(3)+HPVCMIN/2.
      CALL GSPOS('POLY',N,'S2A ',X,Y,Z,0,'ONLY')
C      Z=-S2A(3)+HGLUTOF+S22M(3)
C      CALL GSPOS('S22M',N,'S2A ',X,Y,Z,0,'ONLY')
C      Z=-S2A(3)+HGLUTOF+2*S22M(3)+S21M(3)
C      CALL GSPOS('S21M',N,'S2A ',X,Y,Z,0,'ONLY')
*
      X= 0.
ccC      Z=-S3(3)+HAIRS3+S32M(3)
      Z=-S3(3)+S32M(3)
      DO I=1,NS32X
         N= I
         Y= -S3(2)+ (2*I-1)*S32M(2)
         CALL GSPOS('S32M',N,'S3  ',X,Y,Z,0,'ONLY')
      ENDDO
*
      Y= 0.
ccC      Z= -S3(3)+HAIRS3+2*S32M(3)+HAIRS3+S31M(3)
      Z=-S3(3)+2*S32M(3)+HAIRS3+S31M(3)
      DO I=1,NS31Y
         N= I
         X= -S3(1)+(2*I-1)*S31M(1)
         CALL GSPOS('S31M',N,'S3  ',X,Y,Z,0,'ONLY')
      ENDDO

*
C      Z=-S3(3)+HAIRS3+S32M(3)
C      CALL GSPOS('S32M',N,'S3  ',X,Y,Z,0,'ONLY')
C      Z=-S32(3)+HAIRS3+2*S32M(3)+HAIRS3+S31M(3)
C      CALL GSPOS('S31M',N,'S3  ',X,Y,Z,0,'ONLY')
****************************************
c ml: 21/12/04:
      X=0.
      Y= 0.
      Z= 0.
      N=1
c ml: 27/10/04:
c      DO I=1,NS11Y
c         N= I
c         X= -S11M(1) + HMYLTOF + (2*I-1)*S11Y(1)
      CALL GSPOS('S11Y',N,'S11M',X,Y,Z,0,'ONLY')
c      ENDDO
c end ml.
*
      X= 0.
      Y=0.
      Z= 0.
      N=1
c ml: 27/10/04
C      DO I=1,NS12X
C         N= I
C         Y= -S11M(2) + HMYLTOF + (2*I-1)*S12X(2)
      CALL GSPOS('S12X',N,'S12M',X,Y,Z,0,'ONLY')
C      ENDDO
c end ml.
*
      X=0.
      y=0.
      Z= 0.
      N=1
C      DO I=1,NS21X
C         N= I
c      Y= -S21M(2) + HMYLTOF+ (2*I-1)*S21X(2)
      CALL GSPOS('S21X',N,'S21M',X,Y,Z,0,'ONLY')
C      ENDDO
*
      X=0.
      Y= 0.
      Z= 0.
      N=1
C      DO I=1,NS22Y
C         N= I
C      X= -S22M(1) + HMYLTOF + (2*I-1)*S22Y(1)
      CALL GSPOS('S22Y',N,'S22M',X,Y,Z,0,'ONLY')
C      ENDDO
      Y= 0.
      Z= -S31M(3)+HMYLTOF+S31Y(3)
      N=1
C      DO I=1,NS31Y
C         N= I
C         X= -S31M(1) + HMYLTOF+(2*I-1)*S31Y(1)
      X= -S31M(1) + 2*(HMYLTOF/2)+S31Y(1)
      CALL GSPOS('S31Y',N,'S31M',X,Y,Z,0,'ONLY')
C      ENDDO
*
      X= 0.
      Z= -S32M(3)+HMYLTOF+S32X(3)
      N=1
C      DO I=1,NS32X
C         N= I
C      Y= -S32M(2) + HMYLTOF+ (2*I-1)*S32X(2)
      Y= -S32M(2) + 2*(HMYLTOF/2)+ S32X(2)
      CALL GSPOS('S32X',N,'S32M',X,Y,Z,0,'ONLY')
C      ENDDO
      X=0.
      Y=0.
      Z=0.
      CALL GSPOS('S4  ',N,'S4AL',X,Y,Z,0,'ONLY')
*
      RETURN
      END
1	cafagna	3.1	*
2	pamela	3.7	* $Id: gptofv.F,v 3.6 2006/11/16 18:45:29 pam-ba Exp $
3	cafagna	3.2	*
4			* $Log: gptofv.F,v $
5	pamela	3.7	* Revision 3.6 2006/11/16 18:45:29 pam-ba
6			* Simulated an aluminum container for S4
7			*
8	pam-ba	3.6	* Revision 3.5 2006/05/18 10:52:32 pam-ba
9			* TOF geometry completed and a new material, the polystyrene (density 35 g/l), added
10			*
11	pam-ba	3.5	* Revision 3.4 2006/05/02 10:40:13 pam-ba
12			* TOF geometry and position updated and a new material, the mylar, added
13			*
14	pam-ba	3.4	* Revision 3.3 2006/04/09 23:29:22 cafagna
15			* Several new things, among this: ND and CARD
16			*
17	cafagna	3.3	* Revision 3.2 2005/06/21 02:42:29 cafagna
18			* Major modification to the geometry and to the random number chain
19			*
20	cafagna	3.2	* Revision 3.1.1.1 2002/07/11 16:02:01 cafagna
21			* First GPAMELA release on CVS
22	cafagna	3.1	*
23			*
24			*CMZ : 2.01/00 05/04/2000 14.35.18 by Marialuigia Ambriola
25			*CMZ : 2.00/00 03/03/2000 15.39.05 by Francesco Cafagna
26			*CMZ : 1.02/00 21/10/99 16.30.25 by Francesco Cafagna
27			*CMZ : 1.00/02 18/03/96 17.18.03 by Francesco Cafagna
28			*-- Author : Francesco Cafagna 05/12/95
29			SUBROUTINE GPTOFV
30			************************************************************************
31			* *
32			* Volumes definition for the scintillators counters *
33			* Called by: GPGEO *
34			* Author: Francesco Cafagna, 05/12/95 13.39.41 *
35			* *
36			************************************************************************
37			#include "gpgeo.inc"
38			#include "gpmed.inc"
39			*
40			INTEGER IROT,IVOLU,N,NMED,NUM
41			REAL X,Y,Z
42			*
43	pam-ba	3.4	NMED=MAL
44			*
45			* Define the S1 volume
46			*
47			CALL GSVOLU('S1 ','BOX ',NMED,S1, 3,IVOLU)
48			*
49			* Define the S2 volume
50			*
51			CALL GSVOLU('S2 ','BOX ',NMED,S2, 3,IVOLU)
52			*
53	cafagna	3.1	NMED=MN2
54			*
55	pam-ba	3.4	* Define the S1A volume
56			*
57			CALL GSVOLU('S1A ','BOX ',NMED,S1A, 3,IVOLU)
58			*
59			* Define the S2A volume
60			*
61			CALL GSVOLU('S2A ','BOX ',NMED,S2A, 3,IVOLU)
62			*
63	pam-ba	3.5	* Define the S3 volume
64	cafagna	3.1	*
65	pam-ba	3.5	CALL GSVOLU('S3 ','BOX ',NMED,S3, 3,IVOLU)
66	cafagna	3.1	*
67	pam-ba	3.4	NMED=MMYL
68	cafagna	3.1	*
69	pam-ba	3.4	* Define the S11M volume
70	cafagna	3.1	*
71	pam-ba	3.4	CALL GSVOLU('S11M','BOX ',NMED,S11M, 3,IVOLU)
72	cafagna	3.1	*
73	pam-ba	3.4	* Define the S12M volume
74	cafagna	3.1	*
75	pam-ba	3.4	CALL GSVOLU('S12M','BOX ',NMED,S12M, 3,IVOLU)
76	cafagna	3.1	*
77	pam-ba	3.4	* Define the S21M volume
78	cafagna	3.1	*
79	pam-ba	3.4	CALL GSVOLU('S21M','BOX ',NMED,S21M, 3,IVOLU)
80	cafagna	3.1	*
81	pam-ba	3.4	* Define the S22M volume
82	cafagna	3.1	*
83	pam-ba	3.4	CALL GSVOLU('S22M','BOX ',NMED,S22M, 3,IVOLU)
84			*
85			* Define the S31M volume
86			*
87			CALL GSVOLU('S31M','BOX ',NMED,S31M, 3,IVOLU)
88			*
89			* Define the S32M volume
90			*
91			CALL GSVOLU('S32M','BOX ',NMED,S32M, 3,IVOLU)
92	cafagna	3.1	*
93	pam-ba	3.5	* Define the POL1 volume
94			*
95			NMED=MPOL
96			CALL GSVOLU('POL1','BOX ',NMED,POL1, 3,IVOLU)
97			*
98			* Define the POL2 volume
99			*
100			CALL GSVOLU('POL2','BOX ',NMED,POL2, 3,IVOLU)
101			*
102			* Define the POLY volume
103			*
104			CALL GSVOLU('POLY','BOX ',NMED,POLY, 3,IVOLU)
105			*
106	cafagna	3.1	*
107	pam-ba	3.6	* S4 container volume.
108			*
109			NMED=MAL
110			CALL GSVOLU('S4AL','BOX ',NMED,S4AL, 3,IVOLU)
111			*
112			*
113	cafagna	3.1	* Change medium just for S4.
114			*
115			NMED=MSCIN
116			*
117			* Define the S4 volume
118			*
119			CALL GSVOLU('S4 ','BOX ',NMED,S4, 3,IVOLU)
120			*
121			* Start subvolume definition.
122			*
123			*
124	cafagna	3.2	* Define the S11Y volume
125	cafagna	3.1	*
126	cafagna	3.2	CALL GSVOLU('S11Y','BOX ',NMED,S11Y, 3,IVOLU)
127	cafagna	3.1	*
128	cafagna	3.2	* Define the S12X volume
129	cafagna	3.1	*
130	cafagna	3.2	CALL GSVOLU('S12X','BOX ',NMED,S12X, 3,IVOLU)
131	cafagna	3.1	*
132			* Define the S21X volume
133			*
134			CALL GSVOLU('S21X','BOX ',NMED,S21X, 3,IVOLU)
135			*
136			* Define the S22Y volume
137			*
138			CALL GSVOLU('S22Y','BOX ',NMED,S22Y, 3,IVOLU)
139			*
140			* Define the S31Y volume
141			*
142			CALL GSVOLU('S31Y','BOX ',NMED,S31Y, 3,IVOLU)
143			*
144			* Define the S32X volume
145			*
146			CALL GSVOLU('S32X','BOX ',NMED,S32X, 3,IVOLU)
147
148	pam-ba	3.4	c$$$<<<<<<< gptofv.F
149			c$$$=======
150			c$$$c ml: 21/12/04:
151			c$$$ NMED=MAL
152			c$$$*
153			c$$$* Define the SC11 volume
154			c$$$*
155			c$$$ CALL GSVOLU('SC11','BOX ',NMED,SC11, 3,IVOLU)
156			c$$$*
157			c$$$* Define the SC12 volume
158			c$$$*
159			c$$$ CALL GSVOLU('SC12','BOX ',NMED,SC12, 3,IVOLU)
160			c$$$*
161			c$$$* Define the SC21 volume
162			c$$$*
163			c$$$ CALL GSVOLU('SC21','BOX ',NMED,SC21, 3,IVOLU)
164			c$$$*
165			c$$$* Define the SC22 volume
166			c$$$*
167			c$$$ CALL GSVOLU('SC22','BOX ',NMED,SC22, 3,IVOLU)
168			c$$$
169			c$$$>>>>>>> 3.3
170	cafagna	3.1	*
171	pam-ba	3.4	* Positioning volumes
172	cafagna	3.1	*
173	pam-ba	3.4	X=0.
174			Y=0.
175			Z=0.
176			CALL GSPOS('S1A ',N,'S1 ',X,Y,Z,0,'ONLY')
177	cafagna	3.1	*
178	pam-ba	3.5	N=1
179			X=0.
180			Y=0.
181			Z=-S1A(3)+HGLUTOF+HPVCMIN/2.
182			CALL GSPOS('POL2',N,'S1A ',X,Y,Z,0,'ONLY')
183			*
184	pam-ba	3.4	X=0.
185			Z=-S1A(3)+HGLUTOF+HPVCMIN+S12M(3)
186			DO I=1,NS12X
187	cafagna	3.2	N= I
188	pamela	3.7	c ab mod Y= -S1A(2) + HMYLTOF + (2I-1)S12M(2)
189			Y= -S1A(2) + (2I-1)S12M(2)
190	pam-ba	3.4	CALL GSPOS('S12M',N,'S1A ',X,Y,Z,0,'ONLY')
191	cafagna	3.1	ENDDO
192	pam-ba	3.5	*
193	cafagna	3.1	Y= 0.
194	pam-ba	3.4	Z=-S1A(3)+HGLUTOF+HPVCMIN+2*S12M(3)+S11M(3)
195			DO I=1,NS11Y
196	cafagna	3.1	N= I
197	pamela	3.7	c ab mod X= -S1A(1) + HMYLTOF + (2I-1)S11M(1)
198			X= -S1A(1) + (2I-1)S11M(1)
199	pam-ba	3.4	CALL GSPOS('S11M',N,'S1A ',X,Y,Z,0,'ONLY')
200	cafagna	3.1	ENDDO
201	pam-ba	3.5	*
202			N=1
203			X=0.
204			Y=0.
205			Z=-S1A(3)+HGLUTOF+HPVCMIN+2S12M(3)+2.S11M(3)+HPVCMAX/2.
206			CALL GSPOS('POL1',N,'S1A ',X,Y,Z,0,'ONLY')
207	pam-ba	3.4	C Z=-S1A(3)+HGLUTOF+HPVCMAX+S12M(3)
208	pamela	3.7	C CALL GSPOS('S12M',N,'S1A ',X,Y,Z,0,'ONLY')
209	pam-ba	3.4	C Z=-S1A(3)+HGLUTOF+0.4+2*S12M(3)+S11M(3)
210	pamela	3.7	C CALL GSPOS('S11M',N,'S1A ',X,Y,Z,0,'ONLY')
211			*
212	pam-ba	3.4	X=0.
213			Y=0.
214			Z=0.
215	pamela	3.7	CALL GSPOS('S2A ',N,'S2 ',X,Y,Z,0,'ONLY')
216	cafagna	3.1	*
217			Y= 0.
218	pam-ba	3.4	Z=-S2A(3)+HGLUTOF+S22M(3)
219			DO I=1,NS22Y
220	cafagna	3.1	N= I
221	pamela	3.7	c ab mod X= -S2A(1) + HMYLTOF + (2I-1)S22M(1)
222			X= -S2A(1) + (2I-1)S22M(1)
223	pam-ba	3.4	CALL GSPOS('S22M',N,'S2A ',X,Y,Z,0,'ONLY')
224	cafagna	3.1	ENDDO
225			*
226	pam-ba	3.4	X= 0.
227			Z=-S2A(3)+HGLUTOF+2*S22M(3)+S21M(3)
228	cafagna	3.1	DO I=1,NS21X
229			N= I
230	pamela	3.7	c ab mod Y= -S2A(2) + HMYLTOF+ (2I-1)S21M(2)
231			Y= -S2A(2) + (2I-1)S21M(2)
232	pam-ba	3.4	CALL GSPOS('S21M',N,'S2A ',X,Y,Z,0,'ONLY')
233	cafagna	3.1	ENDDO
234	pam-ba	3.5	*
235			N=1
236			X=0.
237			Y=0.
238			Z=-S2A(3)+HGLUTOF+2S22M(3)+2.S21M(3)+HPVCMIN/2.
239			CALL GSPOS('POLY',N,'S2A ',X,Y,Z,0,'ONLY')
240	pam-ba	3.4	C Z=-S2A(3)+HGLUTOF+S22M(3)
241	pamela	3.7	C CALL GSPOS('S22M',N,'S2A ',X,Y,Z,0,'ONLY')
242	pam-ba	3.4	C Z=-S2A(3)+HGLUTOF+2*S22M(3)+S21M(3)
243	pamela	3.7	C CALL GSPOS('S21M',N,'S2A ',X,Y,Z,0,'ONLY')
244	pam-ba	3.4	*
245	cafagna	3.1	X= 0.
246	pam-ba	3.5	ccC Z=-S3(3)+HAIRS3+S32M(3)
247			Z=-S3(3)+S32M(3)
248	pam-ba	3.4	DO I=1,NS32X
249	cafagna	3.1	N= I
250	pam-ba	3.5	Y= -S3(2)+ (2I-1)S32M(2)
251			CALL GSPOS('S32M',N,'S3 ',X,Y,Z,0,'ONLY')
252	cafagna	3.1	ENDDO
253			*
254	pam-ba	3.4	Y= 0.
255	pam-ba	3.5	ccC Z= -S3(3)+HAIRS3+2*S32M(3)+HAIRS3+S31M(3)
256			Z=-S3(3)+2*S32M(3)+HAIRS3+S31M(3)
257	pam-ba	3.4	DO I=1,NS31Y
258	cafagna	3.1	N= I
259	pam-ba	3.5	X= -S3(1)+(2I-1)S31M(1)
260			CALL GSPOS('S31M',N,'S3 ',X,Y,Z,0,'ONLY')
261	cafagna	3.1	ENDDO
262	pam-ba	3.5
263	pam-ba	3.4	*
264	pam-ba	3.5	C Z=-S3(3)+HAIRS3+S32M(3)
265	pamela	3.7	C CALL GSPOS('S32M',N,'S3 ',X,Y,Z,0,'ONLY')
266	pam-ba	3.4	C Z=-S32(3)+HAIRS3+2*S32M(3)+HAIRS3+S31M(3)
267	pamela	3.7	C CALL GSPOS('S31M',N,'S3 ',X,Y,Z,0,'ONLY')
268	pam-ba	3.4	****************************************
269			c ml: 21/12/04:
270			X=0.
271			Y= 0.
272			Z= 0.
273			N=1
274			c ml: 27/10/04:
275			c DO I=1,NS11Y
276			c N= I
277			c X= -S11M(1) + HMYLTOF + (2I-1)S11Y(1)
278			CALL GSPOS('S11Y',N,'S11M',X,Y,Z,0,'ONLY')
279			c ENDDO
280			c end ml.
281			*
282	cafagna	3.1	X= 0.
283	pam-ba	3.4	Y=0.
284			Z= 0.
285			N=1
286			c ml: 27/10/04
287			C DO I=1,NS12X
288			C N= I
289			C Y= -S11M(2) + HMYLTOF + (2I-1)S12X(2)
290			CALL GSPOS('S12X',N,'S12M',X,Y,Z,0,'ONLY')
291			C ENDDO
292			c end ml.
293			*
294			X=0.
295			y=0.
296			Z= 0.
297			N=1
298			C DO I=1,NS21X
299			C N= I
300			c Y= -S21M(2) + HMYLTOF+ (2I-1)S21X(2)
301			CALL GSPOS('S21X',N,'S21M',X,Y,Z,0,'ONLY')
302			C ENDDO
303			*
304			X=0.
305	cafagna	3.1	Y= 0.
306	pam-ba	3.4	Z= 0.
307			N=1
308			C DO I=1,NS22Y
309			C N= I
310	pamela	3.7	C X= -S22M(1) + HMYLTOF + (2I-1)S22Y(1)
311	pam-ba	3.4	CALL GSPOS('S22Y',N,'S22M',X,Y,Z,0,'ONLY')
312			C ENDDO
313	cafagna	3.1	Y= 0.
314	pam-ba	3.4	Z= -S31M(3)+HMYLTOF+S31Y(3)
315			N=1
316			C DO I=1,NS31Y
317			C N= I
318			C X= -S31M(1) + HMYLTOF+(2I-1)S31Y(1)
319			X= -S31M(1) + 2*(HMYLTOF/2)+S31Y(1)
320			CALL GSPOS('S31Y',N,'S31M',X,Y,Z,0,'ONLY')
321			C ENDDO
322	cafagna	3.1	*
323			X= 0.
324	pam-ba	3.4	Z= -S32M(3)+HMYLTOF+S32X(3)
325			N=1
326			C DO I=1,NS32X
327			C N= I
328			C Y= -S32M(2) + HMYLTOF+ (2I-1)S32X(2)
329			Y= -S32M(2) + 2*(HMYLTOF/2)+ S32X(2)
330			CALL GSPOS('S32X',N,'S32M',X,Y,Z,0,'ONLY')
331			C ENDDO
332	pam-ba	3.6	X=0.
333			Y=0.
334			Z=0.
335	pamela	3.7	CALL GSPOS('S4 ',N,'S4AL',X,Y,Z,0,'ONLY')
336	cafagna	3.1	*
337			RETURN
338			END