/[PAMELA software]/gpamela/gpcal/gpcalv.F
ViewVC logotype

Diff of /gpamela/gpcal/gpcalv.F

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 3.2 by pamela, Thu Dec 5 10:17:41 2002 UTC revision 3.5 by cafagna, Tue May 9 12:01:15 2006 UTC
# Line 1  Line 1 
1  *  *
2  * $Id: gpcalv.F,v 3.1.1.1 2002/07/11 16:02:14 cafagna Exp $  * $Id: gpcalv.F,v 3.4 2005/07/25 11:53:20 cafagna Exp $
3  *  *
4  * $Log: gpcalv.F,v $  * $Log: gpcalv.F,v $
5    * Revision 3.4  2005/07/25 11:53:20  cafagna
6    * Several updates. See history for details
7    *
8    * Revision 3.3  2003/12/17 11:32:48  pamela
9    * CALO SIMULATION COMPLETED: geometry and special tracking parameters updated and simulation checked by a comparison with the Trieste's standalone Monte Carlo simulation
10    *
11    * Revision 3.2  2002/12/05 10:17:41  pamela
12    * Update CAS and CALO geometries and positions. Makefile updated as well
13    *
14  * Revision 3.1.1.1  2002/07/11 16:02:14  cafagna  * Revision 3.1.1.1  2002/07/11 16:02:14  cafagna
15  * First GPAMELA release on CVS  * First GPAMELA release on CVS
16  *  *
# Line 27  Line 36 
36  *  *
37        INTEGER IROT,IVOLU,N,NMED        INTEGER IROT,IVOLU,N,NMED
38        REAL X,Y,Z        REAL X,Y,Z
39    * MA
40          REAL X1, Y1
41          INTEGER m
42    * END MA
43  *  *
44  *  Define the CALB     volume  *  Define the CALB     volume
45  *  *
# Line 107  C      CALL GSDVN ('CAST','CASI',NCASTR, Line 120  C      CALL GSDVN ('CAST','CASI',NCASTR,
120        Y = 0.        Y = 0.
121  *EM:  *EM:
122  *F.V.&ML:  *F.V.&ML:
123        CALL GPVOLU(-2)  c      CALL GPVOLU(-2)
124  C      Z = -CAPL(3) + CAKA(3)  C      Z = -CAPL(3) + CAKA(3)
125  C      CALL GSPOS('CAKA',N,'CAPL',X,Y,Z,0,'ONLY')  C      CALL GSPOS('CAKA',N,'CAPL',X,Y,Z,0,'ONLY')
126        Z = -CANS(3) + CAKA(3)        Z = -CANS(3) + CAKA(3)
# Line 145  C      Z= Z + CAKP(3) + CASI(3) Line 158  C      Z= Z + CAKP(3) + CASI(3)
158           ENDDO           ENDDO
159        ENDDO        ENDDO
160  *END EM.  *END EM.
161  *F.V.&ML:  * MA
162  C commenting the original code with a 'C' per raw:  * 11 MODULI DI 2 W
163  C*  c     Z = CALB(3) - CAPL(3) - DCASIZ/2
164  C* Positioning volume CAPL and CAAB into the mother CALB        Z = CALB(3) - CAPL(3) - DCASIZ
165  C*        N = 0
166  C      X= 0.        DO M = 0,20,2
167  C      Y= 0.  *     Positioning CAPL CANS CAAB in CALB to form a module
168  CC      Z= CALB(3) + CAPL(3)           X = .10
169  C*EM:           y = .05
170  C      Z= CALB(3) + CAPL(3) - DCASIZ/2.           DO I = 1,2
171  C*END EM. (we have 0.15 cm from the top to the first Si)  *     Positioning CAPL
172  C      N= 0              N = N + 1
173  C      DO I=1,NCAPL  C     if(N.EQ.1) print *,'z ini = ',Z
174  C*              IF(.NOT.NOCAPL(N))
175  C* X view, absorber and Y view on the other       +           CALL GSPOS('CAPL',N,'CALB',X,(Y-0.1),Z,0,'ONLY')
176  C*  *     Positioning CANS
177  C         Z= Z - 2*CAPL(3)              Z = Z - CAPL(3) - CANS(3)
178  C         N= N + 1              CALL GSPOS('CANS',N,'CALB',0,-0.1,Z,0,'ONLY')
179  C*  *     Positioning CAAB
180  C* Check if this plane is excluded              Z = Z - CANS(3) - CAAB(3)
181  C*              IF(.NOT.NOCAAB(M+I))
182  C         IF(.NOT.NOCAPL(N))       +           CALL GSPOS('CAAB',(M+I),'CALB',0,-0.1,Z,0,'ONLY')
183  C     +        CALL GSPOS('CAPL',N,'CALB',X,Y,Z,0,'ONLY')  *     Positioning CANS (ruotato)
184  C         Z= Z - CAPL(3) - CAAB(3)              N = N + 1
185  C*              Z = Z - CAAB(3) - CANS(3)
186  C* Check if this absorber is excluded              CALL GSPOS('CANS',N,'CALB',0,-0.1,Z,3,'ONLY')
187  C*  *     Positioning CAPL (ruotato e shiftato)
188  C         IF(.NOT.NOCAAB(I))              X1 = X
189  C     +        CALL GSPOS('CAAB',(I-1),'CALB',X,Y,Z,0,'ONLY')              Y1 = Y
190  C         Z= Z - CAAB(3) - CAPL(3)              X = -Y1
191  C         N= N + 1              Y = X1
192  C*              Z = Z - CANS(3) - CAPL(3)
193  C* Check if this detector is excluded  C     if(N.EQ.44) print *,'z fin = ',Z
194  C*              IF(.NOT.NOCAPL(N))
195  C         IF(.NOT.NOCAPL(N))       +           CALL GSPOS('CAPL',N,'CALB',X,(Y-0.1),Z,3,'ONLY')
196  C     +        CALL GSPOS('CAPL',N,'CALB',X,Y,Z,3,'ONLY')  
197  C      ENDDO              IF(I.EQ.1) THEN
198  C end comments.  c              Z = Z - 2*CAPL(3) - DZMO ! SPESSORE TRA I SILICI CHE SI GUARDANO
199  *                 Z = Z - 2*CAPL(3) - DZM0 ! corretto
200  * Positioning volume CANS ,CAPL and CAAB into the mother CALB                 X1 = X                   ! SHIFT PER IL CAPL PROSSIMO CICLO
201  *                 Y1 = Y
202        X= 0.                 X = -Y1
203        Y= 0.                 Y = X1
204        Z= CALB(3) + CAPL(3) - DCASIZ/2.              ENDIF
205  *END EM. (we have 0.15 cm from the top to the first Si)           ENDDO
206        N= 0           Z = Z - DZST - 2*CAPL(3)
 *  
 *       I DISPARI  
       DO I=1,NCAPL,2  
 *  
 * X view, absorber and Y view on the other  
 *  
          X=0.10  
          Y=0.05  
          Z= Z - 2*CAPL(3)  
          N= N + 1  
 *  
 * Check if this plane is excluded  
 *  
 *  
          IF(.NOT.NOCAPL(N))  
      +        CALL GSPOS('CAPL',N,'CALB',X,Y,Z,0,'ONLY')  
 *  
          X = 0.  
          Y = 0.  
          Z= Z - CAPL(3) - CANS(3)  
 *  
          CALL GSPOS('CANS',N,'CALB',X,Y,Z,0,'ONLY')  
 *  
          Z= Z - CANS(3) - CAAB(3)  
 *  
 * Check if this absorber is excluded  
 *  
 *  
          IF(.NOT.NOCAAB(I))  
      +        CALL GSPOS('CAAB',(I-1),'CALB',X,Y,Z,0,'ONLY')  
 *  
          Z = Z - CAAB(3)- CANS(3)  
          N = N + 1  
 *      
          CALL GSPOS('CANS',N,'CALB',X,Y,Z,3,'ONLY')  
          Z= Z - CANS(3) - CAPL(3)  
          X=0.05  
          Y=-0.10  
 *  
 * Check if this detector is excluded  
 *  
          IF(.NOT.NOCAPL(N))  
      +        CALL GSPOS('CAPL',N,'CALB',X,Y,Z,3,'ONLY')  
          N=N+2  
          Z = Z -2*(2*CAPL(3)+2*CANS(3)+CAAB(3))  
       ENDDO  
 *  
 *       I PARI  
 *  
       X= 0.  
       Y= 0.  
       Z= CALB(3) + CAPL(3) - DCASIZ/2.  
 *END EM. (we have 0.15 cm from the top to the first Si)  
       N= 2  
 *  
 *  
       DO I=2,NCAPL,2  
 *  
 * X view, absorber and Y view on the other  
 *  
          Z= Z - 2*CAPL(3)-2*(2*CAPL(3)+2*CANS(3)+CAAB(3))  
          N= N + 1  
 *  
 * Check if this plane is excluded  
          X = -0.10  
          Y = -0.05  
 *  
          IF(.NOT.NOCAPL(N))  
      +        CALL GSPOS('CAPL',N,'CALB',X,Y,Z,0,'ONLY')  
          X = 0.  
          Y = 0.  
          Z = Z - CAPL(3) - CANS(3)  
 *  
          CALL GSPOS('CANS',N,'CALB',X,Y,Z,0,'ONLY')  
 *  
 * Check if this absorber is excluded  
 *  
          N = N + 1  
          Z = Z - CANS(3) - CAAB(3)  
 *      
          IF(.NOT.NOCAAB(I))  
      +        CALL GSPOS('CAAB',(I-1),'CALB',X,Y,Z,0,'ONLY')  
 *  
          Z= Z - CAAB(3) - CANS(3)  
 *  
          CALL GSPOS('CANS',N,'CALB',X,Y,Z,3,'ONLY')  
          Z= Z - CANS(3) - CAPL(3)  
          X = -0.05  
          Y = -0.10  
 *  
 *  
 * Check if this detector is excluded  
 *  
          IF(.NOT.NOCAPL(N))  
      +        CALL GSPOS('CAPL',N,'CALB',X,Y,Z,3,'ONLY')  
          N = N + 2  
 *  
207        ENDDO        ENDDO
208  *END F.V.&ML.  *     END MA
209  * F.V.&ML.:  
210  * again we have to comment with a 'C' per raw:  C # *F.V.&ML:
211  C*  C # C commenting the original code with a 'C' per raw:
212  C* Positioning volume CAPD and CAAD into the mother CALB  C # C*
213  C*  C # C* Positioning volume CAPL and CAAB into the mother CALB
214  C      N= 0  C # C*
215  C      DO I=1,NCAPLD  C # C      X= 0.
216  C         IF(I.EQ.1) THEN  C # C      Y= 0.
217  C            Z= Z - CAPD(3) - CAPL(3)  C # CC      Z= CALB(3) + CAPL(3)
218  C         ELSE  C # C*EM:
219  C            Z= Z - 2*CAPD(3)  C # C      Z= CALB(3) + CAPL(3) - DCASIZ/2.
220  C         ENDIF  C # C*END EM. (we have 0.15 cm from the top to the first Si)
221  C         N= N + 1  C # C      N= 0
222  C         CALL GSPOS('CAPD',N,'CALB',X,Y,Z,0,'ONLY')  C # C      DO I=1,NCAPL
223  C         Z= Z - CAPD(3) - CAAD(3)  C # C*
224  C         CALL GSPOS('CAAD',(I-1),'CALB',X,Y,Z,0,'ONLY')  C # C* X view, absorber and Y view on the other
225  C         Z= Z - CAAD(3) - CAPD(3)  C # C*
226  C         N= N + 1  C # C         Z= Z - 2*CAPL(3)
227  C         CALL GSPOS('CAPD',N,'CALB',X,Y,Z,3,'ONLY')  C # C         N= N + 1
228  C      ENDDO  C # C*
229  * end comments.  C # C* Check if this plane is excluded
230  *END F.V.&ML.  C # C*
231    C # C         IF(.NOT.NOCAPL(N))
232    C # C     +        CALL GSPOS('CAPL',N,'CALB',X,Y,Z,0,'ONLY')
233    C # C         Z= Z - CAPL(3) - CAAB(3)
234    C # C*
235    C # C* Check if this absorber is excluded
236    C # C*
237    C # C         IF(.NOT.NOCAAB(I))
238    C # C     +        CALL GSPOS('CAAB',(I-1),'CALB',X,Y,Z,0,'ONLY')
239    C # C         Z= Z - CAAB(3) - CAPL(3)
240    C # C         N= N + 1
241    C # C*
242    C # C* Check if this detector is excluded
243    C # C*
244    C # C         IF(.NOT.NOCAPL(N))
245    C # C     +        CALL GSPOS('CAPL',N,'CALB',X,Y,Z,3,'ONLY')
246    C # C      ENDDO
247    C # C end comments.
248    C # *
249    C # * Positioning volume CANS ,CAPL and CAAB into the mother CALB
250    C # *
251    C #       X= 0.
252    C #       Y= 0.
253    C #       Z= CALB(3) + CAPL(3) - DCASIZ/2.
254    C # *END EM. (we have 0.15 cm from the top to the first Si)
255    C #       N= 0
256    C # *
257    C # *   I DISPARI
258    C #       DO I=1,NCAPL,2
259    C # *
260    C # * X view, absorber and Y view on the other
261    C # *
262    C #      X=0.10
263    C #      Y=0.05
264    C #          Z= Z - 2*CAPL(3)
265    C #          N= N + 1
266    C # *
267    C # * Check if this plane is excluded
268    C # *
269    C # *
270    C #          IF(.NOT.NOCAPL(N))
271    C #      +        CALL GSPOS('CAPL',N,'CALB',X,Y,Z,0,'ONLY')
272    C # *
273    C #          X = 0.
274    C #          Y = 0.
275    C #          Z= Z - CAPL(3) - CANS(3)
276    C # *
277    C #          CALL GSPOS('CANS',N,'CALB',X,Y,Z,0,'ONLY')
278    C # *
279    C #          Z= Z - CANS(3) - CAAB(3)
280    C # *
281    C # * Check if this absorber is excluded
282    C # *
283    C # *
284    C #          IF(.NOT.NOCAAB(I))
285    C #      +        CALL GSPOS('CAAB',(I-1),'CALB',X,Y,Z,0,'ONLY')
286    C # *
287    C #          Z = Z - CAAB(3)- CANS(3)
288    C #          N = N + 1
289    C # *    
290    C #          CALL GSPOS('CANS',N,'CALB',X,Y,Z,3,'ONLY')
291    C #          Z= Z - CANS(3) - CAPL(3)
292    C #          X=-0.05
293    C #          Y=0.10
294    C # *
295    C # * Check if this detector is excluded
296    C # *
297    C #          IF(.NOT.NOCAPL(N))
298    C #      +        CALL GSPOS('CAPL',N,'CALB',X,Y,Z,3,'ONLY')
299    C #      N=N+2
300    C #      Z = Z -2*(2*CAPL(3)+2*CANS(3)+CAAB(3))
301    C #       ENDDO
302    C # *
303    C # *   I PARI
304    C # *
305    C #       X= 0.
306    C #       Y= 0.
307    C #       Z= CALB(3) + CAPL(3) - DCASIZ/2.
308    C # *END EM. (we have 0.15 cm from the top to the first Si)
309    C #       N= 2
310    C # *
311    C # *
312    C #       DO I=2,NCAPL,2
313    C # *
314    C # * X view, absorber and Y view on the other
315    C # *
316    C #          Z= Z - 2*CAPL(3)-2*(2*CAPL(3)+2*CANS(3)+CAAB(3))
317    C #          N= N + 1
318    C # *
319    C # * Check if this plane is excluded
320    C #          X= -0.10
321    C #          Y= -0.05
322    C # *
323    C #          IF(.NOT.NOCAPL(N))
324    C #      +        CALL GSPOS('CAPL',N,'CALB',X,Y,Z,0,'ONLY')
325    C #      X = 0.
326    C #      Y = 0.
327    C #          Z = Z - CAPL(3) - CANS(3)
328    C # *
329    C #              CALL GSPOS('CANS',N,'CALB',X,Y,Z,0,'ONLY')
330    C # *
331    C # * Check if this absorber is excluded
332    C # *
333    C #          N = N + 1
334    C #          Z = Z - CANS(3) - CAAB(3)
335    C # *    
336    C #          IF(.NOT.NOCAAB(I))
337    C #      +        CALL GSPOS('CAAB',(I-1),'CALB',X,Y,Z,0,'ONLY')
338    C # *
339    C #      Z= Z - CAAB(3) - CANS(3)
340    C # *
341    C #          CALL GSPOS('CANS',N,'CALB',X,Y,Z,3,'ONLY')
342    C #          Z= Z - CANS(3) - CAPL(3)
343    C #          X= 0.05
344    C #          Y= -0.10
345    C # *
346    C # *
347    C # * Check if this detector is excluded
348    C # *
349    C #          IF(.NOT.NOCAPL(N))
350    C #      +        CALL GSPOS('CAPL',N,'CALB',X,Y,Z,3,'ONLY')
351    C #          N = N + 2
352    C # *
353    C #       ENDDO
354    C # *END F.V.&ML.
355    C # * F.V.&ML.:
356    C # * again we have to comment with a 'C' per raw:
357    C # C*
358    C # C* Positioning volume CAPD and CAAD into the mother CALB
359    C # C*
360    C # C      N= 0
361    C # C      DO I=1,NCAPLD
362    C # C         IF(I.EQ.1) THEN
363    C # C            Z= Z - CAPD(3) - CAPL(3)
364    C # C         ELSE
365    C # C            Z= Z - 2*CAPD(3)
366    C # C         ENDIF
367    C # C         N= N + 1
368    C # C         CALL GSPOS('CAPD',N,'CALB',X,Y,Z,0,'ONLY')
369    C # C         Z= Z - CAPD(3) - CAAD(3)
370    C # C         CALL GSPOS('CAAD',(I-1),'CALB',X,Y,Z,0,'ONLY')
371    C # C         Z= Z - CAAD(3) - CAPD(3)
372    C # C         N= N + 1
373    C # C         CALL GSPOS('CAPD',N,'CALB',X,Y,Z,3,'ONLY')
374    C # C      ENDDO
375    C # * end comments.
376    C # *END F.V.&ML.
377    C # end MA comments
378  *  *
379  * Put the shell around it  * Put the shell around it
380  *  *
Legend:
Removed from v.3.2  
changed lines
  Added in v.3.5
  ViewVC Help
Powered by ViewVC 1.1.23