/[PAMELA software]/gpamela/gpcard/gpcardv.F
ViewVC logotype

Diff of /gpamela/gpcard/gpcardv.F

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

revision 1.3 by cafagna, Sun Dec 17 14:53:05 2006 UTC revision 1.4 by pamela, Tue Jan 29 18:25:16 2008 UTC
# Line 24  Line 24 
24        H2=CAR(1)-H1        H2=CAR(1)-H1
25  *  *
26  * Greater Scintillator CAR1 (or type B)  * Greater Scintillator CAR1 (or type B)
27  *            *
28        B1=CARB(2)+CAR(4)+CAR(5)+H1/TAN(DEGRAD*CAR(7))        B1=CARB(2)+CAR(4)+CAR(5)+H1/TAN(DEGRAD*CAR(7))
29        C1=H2/TAN(DEGRAD*CAR(7))        C1=H2/TAN(DEGRAD*CAR(7))
30        A1=CARB(1)-C1        A1=CARB(1)-C1
31        D1=CARB(1)-CARB(2)-CAR(4)-CAR(5)        D1=CARB(1)-CARB(2)-CAR(4)-CAR(5)
32        F1=C1+H1/TAN(DEGRAD*CAR(7))        F1=C1+H1/TAN(DEGRAD*CAR(7))
33        ANG1C=ATAN((B1/2.+D1-A1/2.)/CAR(1))        ANG1C=ATAN((B1/2.+D1-A1/2.)/CAR(1))
34          
35        C1D1(1) = CAR(1)/2.        C1D1(1) = CAR(1)/2.
36        C1D1(2) = -ANG1C*RADDEG        C1D1(2) = -ANG1C*RADDEG
37        C1D1(3) = 0.        C1D1(3) = 0.
# Line 46  Line 46 
46    
47        NMED=MSCIN        NMED=MSCIN
48        CALL GSVOLU('C1D1','TRAP',NMED,C1D1,11,IVOLU)        CALL GSVOLU('C1D1','TRAP',NMED,C1D1,11,IVOLU)
49  *        *
50  * Smaller Scintillator CAR2 (or type A)  * Smaller Scintillator CAR2 (or type A)
51  *            *
52        B2=CARA(2)+CAR(4)+CAR(5)+H1/TAN(DEGRAD*CAR(7))        B2=CARA(2)+CAR(4)+CAR(5)+H1/TAN(DEGRAD*CAR(7))
53        C2=H2/TAN(DEGRAD*CAR(7))        C2=H2/TAN(DEGRAD*CAR(7))
54        A2=CARA(1)-C2        A2=CARA(1)-C2
# Line 95  Line 95 
95        NMED=MAL        NMED=MAL
96        CALL GSVOLU('CAR1','TRAP',NMED,CAR1,11,IVOLU)        CALL GSVOLU('CAR1','TRAP',NMED,CAR1,11,IVOLU)
97  *  *
98  * CR1P PLASTIC BOX TO CREATE RIM  * CR1P PLASTIC BOX TO CREATE RIM
99  *  *
100        BB=C1D1(9)+PT*(1./COS(ATAN(D1/CAR(1)))-        BB=C1D1(9)+PT*(1./COS(ATAN(D1/CAR(1)))-
101       +              D1/CAR(1)+1./TAN(DEGRAD*CAR(7))+       +              D1/CAR(1)+1./TAN(DEGRAD*CAR(7))+
# Line 139  Line 139 
139        CAR2(9) = BB        CAR2(9) = BB
140        CAR2(10) = CAR2(9)        CAR2(10) = CAR2(9)
141        CAR2(11) = 0.        CAR2(11) = 0.
142                
143        NMED=MAL        NMED=MAL
144        CALL GSVOLU('CAR2','TRAP',NMED,CAR2,11,IVOLU)        CALL GSVOLU('CAR2','TRAP',NMED,CAR2,11,IVOLU)
145  *  *
146  * CR1P PLASTIC BOX TO CREATE RIM  * CR1P PLASTIC BOX TO CREATE RIM
147  *  *
148        BB=C2D1(9)+PT*(1./COS(ATAN(D2/CAR(1)))-        BB=C2D1(9)+PT*(1./COS(ATAN(D2/CAR(1)))-
149       +              D2/CAR(1)+1./TAN(DEGRAD*CAR(7))+       +              D2/CAR(1)+1./TAN(DEGRAD*CAR(7))+
# Line 164  Line 164 
164        CR2P(10) = CR2P(9)        CR2P(10) = CR2P(9)
165        CR2P(11) = 0.        CR2P(11) = 0.
166    
167        NMED=MPLAS              NMED=MPLAS
168        CALL GSVOLU('CR2P','TRAP',NMED,CR2P,11,IVOLU)        CALL GSVOLU('CR2P','TRAP',NMED,CR2P,11,IVOLU)
169  *  *
170  *  Volumes to be subtracted from C1D1 and C2D1  *  Volumes to be subtracted from C1D1 and C2D1
171  *  *
172        AA=CAR(2)+CAR(3)        AA=CAR(2)+CAR(3)
173        BB=AA/TAN(DEGRAD*CAR(7))+AA*CAR(4)/CAR(2)        BB=AA/TAN(DEGRAD*CAR(7))+AA*CAR(4)/CAR(2)
174        
175        VSN1(1) = 0.5*AA        VSN1(1) = 0.5*AA
176        VSN1(2) = RADDEG*ATAN((BB/2.-AA/TAN(DEGRAD*CAR(7)))/AA)        VSN1(2) = RADDEG*ATAN((BB/2.-AA/TAN(DEGRAD*CAR(7)))/AA)
177        VSN1(3) = 0.        VSN1(3) = 0.
# Line 184  Line 184 
184        VSN1(10) = VSN1(9)        VSN1(10) = VSN1(9)
185        VSN1(11) = 0.        VSN1(11) = 0.
186    
187        AA=CAR(3)        AA=CAR(3)
188        BB=AA*(CAR(5)/CAR(3)-CAR(4)/CAR(2))        BB=AA*(CAR(5)/CAR(3)-CAR(4)/CAR(2))
189        
190        VSN2(1) = 0.5*AA        VSN2(1) = 0.5*AA
191        VSN2(2) = RADDEG*ATAN((BB/2.+(AA)*CAR(4)/CAR(2))/AA)        VSN2(2) = RADDEG*ATAN((BB/2.+(AA)*CAR(4)/CAR(2))/AA)
192        VSN2(3) = 0.        VSN2(3) = 0.
# Line 209  Line 209 
209        AA=(PT/SIN(DEGRAD*CAR(7)))/(1./TAN(DEGRAD*CAR(7))+        AA=(PT/SIN(DEGRAD*CAR(7)))/(1./TAN(DEGRAD*CAR(7))+
210       +    CAR(4)/CAR(2))+CAR(2)+CAR(3)+PT       +    CAR(4)/CAR(2))+CAR(2)+CAR(3)+PT
211        BB=AA/TAN(DEGRAD*CAR(7))+AA*CAR(4)/CAR(2)        BB=AA/TAN(DEGRAD*CAR(7))+AA*CAR(4)/CAR(2)
212      
213        VPN1(1) = 0.5*AA        VPN1(1) = 0.5*AA
214        VPN1(2) = VSN1(2)        VPN1(2) = VSN1(2)
215        VPN1(3) = 0.        VPN1(3) = 0.
# Line 221  Line 221 
221        VPN1(9) = 0.5*BB        VPN1(9) = 0.5*BB
222        VPN1(10) = VPN1(9)        VPN1(10) = VPN1(9)
223        VPN1(11) = 0.        VPN1(11) = 0.
224          
225        AA=CAR(3) + PT -        AA=CAR(3) + PT -
226       +    COS(ATAN(CAR(4)/CAR(2)))*PT/SIN(ATAN(CAR(5)/CAR(3))-       +    COS(ATAN(CAR(4)/CAR(2)))*PT/SIN(ATAN(CAR(5)/CAR(3))-
227       +    ATAN(CAR(4)/CAR(2)))       +    ATAN(CAR(4)/CAR(2)))
228        BB=AA*(CAR(5)/CAR(3)-CAR(4)/CAR(2))        BB=AA*(CAR(5)/CAR(3)-CAR(4)/CAR(2))
229        
230        VPN2(1) = 0.5*AA        VPN2(1) = 0.5*AA
231        VPN2(2) = RADDEG*ATAN((BB/2.+AA*CAR(4)/CAR(2))/AA)        VPN2(2) = RADDEG*ATAN((BB/2.+AA*CAR(4)/CAR(2))/AA)
232        VPN2(3) = 0.        VPN2(3) = 0.
# Line 249  Line 249 
249        AA=(ATZ/SIN(DEGRAD*CAR(7)))/(1./TAN(DEGRAD*CAR(7))+        AA=(ATZ/SIN(DEGRAD*CAR(7)))/(1./TAN(DEGRAD*CAR(7))+
250       +    CAR(4)/CAR(2))+CAR(2)+CAR(3)+ATZ       +    CAR(4)/CAR(2))+CAR(2)+CAR(3)+ATZ
251        BB=AA/TAN(DEGRAD*CAR(7))+AA*CAR(4)/CAR(2)        BB=AA/TAN(DEGRAD*CAR(7))+AA*CAR(4)/CAR(2)
252        
253        VAN1(1) = 0.5*AA        VAN1(1) = 0.5*AA
254        VAN1(2) = VSN1(2)        VAN1(2) = VSN1(2)
255        VAN1(3) = 0.        VAN1(3) = 0.
# Line 262  Line 262 
262        VAN1(10) = VAN1(9)        VAN1(10) = VAN1(9)
263        VAN1(11) = 0.        VAN1(11) = 0.
264    
265        AA=CAR(3) + ATZ -        AA=CAR(3) + ATZ -
266       +    COS(ATAN(CAR(4)/CAR(2)))*ATZ/SIN(ATAN(CAR(5)/CAR(3))-       +    COS(ATAN(CAR(4)/CAR(2)))*ATZ/SIN(ATAN(CAR(5)/CAR(3))-
267       +    ATAN(CAR(4)/CAR(2)))       +    ATAN(CAR(4)/CAR(2)))
268        BB=AA*(CAR(5)/CAR(3)-CAR(4)/CAR(2))        BB=AA*(CAR(5)/CAR(3)-CAR(4)/CAR(2))
269        
270        VAN2(1) = 0.5*AA        VAN2(1) = 0.5*AA
271        VAN2(2) = RADDEG*ATAN((BB/2.+(AA)*CAR(4)/CAR(2))/AA)        VAN2(2) = RADDEG*ATAN((BB/2.+(AA)*CAR(4)/CAR(2))/AA)
272        VAN2(3) = 0.        VAN2(3) = 0.
# Line 279  Line 279 
279        VAN2(10) = VAN2(9)        VAN2(10) = VAN2(9)
280        VAN2(11) = 0.        VAN2(11) = 0.
281    
282        NMED=MN2            NMED=MN2
283        CALL GSVOLU('VAN1','TRAP',NMED,VAN1,11,IVOLU)        CALL GSVOLU('VAN1','TRAP',NMED,VAN1,11,IVOLU)
284        CALL GSVOLU('VAN2','TRAP',NMED,VAN2,11,IVOLU)        CALL GSVOLU('VAN2','TRAP',NMED,VAN2,11,IVOLU)
285  *  *
# Line 354  Line 354 
354        X=-(CAR1(9)-CAR1(1)*TAN(CAR1(2)*DEGRAD)-        X=-(CAR1(9)-CAR1(1)*TAN(CAR1(2)*DEGRAD)-
355       +     (2.*VAN1(9)+VAN2(9)-VAN2(1)*TAN(VAN2(2)*DEGRAD)))       +     (2.*VAN1(9)+VAN2(9)-VAN2(1)*TAN(VAN2(2)*DEGRAD)))
356        CALL GSPOS('VAN2',N,'CAR1',X,Y,Z,0,'ONLY')        CALL GSPOS('VAN2',N,'CAR1',X,Y,Z,0,'ONLY')
357  *    *
358        N=1        N=1
359        X=PT*(1./COS(ATAN(D1/CAR(1)))-        X=PT*(1./COS(ATAN(D1/CAR(1)))-
360       +    (1./SIN(DEGRAD*CAR(7))))/2.       +    (1./SIN(DEGRAD*CAR(7))))/2.
# Line 412  Line 412 
412        X=-(CAR2(9)-CAR2(1)*TAN(CAR2(2)*DEGRAD)-        X=-(CAR2(9)-CAR2(1)*TAN(CAR2(2)*DEGRAD)-
413       +   (2.*VAN1(9)+VAN2(9)-VAN2(1)*TAN(VAN2(2)*DEGRAD)))       +   (2.*VAN1(9)+VAN2(9)-VAN2(1)*TAN(VAN2(2)*DEGRAD)))
414        CALL GSPOS('VAN2',N,'CAR2',X,Y,Z,0,'ONLY')        CALL GSPOS('VAN2',N,'CAR2',X,Y,Z,0,'ONLY')
415  *    *
416        N=1        N=1
417        X=PT*(1./COS(ATAN(D2/CAR(1)))-1./SIN(DEGRAD*CAR(7)))/2.        X=PT*(1./COS(ATAN(D2/CAR(1)))-1./SIN(DEGRAD*CAR(7)))/2.
418        Y=0        Y=0
# Line 434  Line 434 
434  *           CAR1B  *           CAR1B
435  *  *
436        ZH=49.229-ZTPLA-TPLA(3)        ZH=49.229-ZTPLA-TPLA(3)
437  *--- CAR1B        *--- CAR1B
438        N=1        N=1
439        Y=-CARB(5)+CAR(8)/2.*SIN(ANGLCAR1)        Y=-CARB(5)+CAR(8)/2.*SIN(ANGLCAR1)
440        Z=CARB(6)+ZH-CARDB(3)-CAR(8)*COS(ANGLCAR1)/2.        Z=CARB(6)+ZH-CARDB(3)-CAR(8)*COS(ANGLCAR1)/2.
441        X=((CARB(7)/2.)-TAN(ANG1C)*CAR(8)/2.)-CARB(4)        X=((CARB(7)/2.)-TAN(ANG1C)*CAR(8)/2.)-CARB(4)
442        CALL GSPOS('CAR1',N,'CARD',X,Y,Z,IRZCAR1B,'ONLY')        CALL GSPOS('CAR1',N,'CARD',X,Y,Z,IRZCAR1B,'ONLY')
443  *--- CAR1A        *--- CAR1A
444        N=2        N=2
445        X=-X        X=-X
446        Y=-Y        Y=-Y
# Line 451  Line 451 
451        Z=CARA(6)+ZH-CARDB(3)-CAR(8)*COS(ANGLCAR2)/2.        Z=CARA(6)+ZH-CARDB(3)-CAR(8)*COS(ANGLCAR2)/2.
452        Y=(CARA(7)/2.-TAN(ANG2C)*CAR(8)/2.)-CARA(5)        Y=(CARA(7)/2.-TAN(ANG2C)*CAR(8)/2.)-CARA(5)
453        CALL GSPOS('CAR2',N,'CARD',X,Y,Z,IRZCAR2A,'ONLY')        CALL GSPOS('CAR2',N,'CARD',X,Y,Z,IRZCAR2A,'ONLY')
454  *--- CAR2B        *--- CAR2B
455        N=2        N=2
456        X=-X        X=-X
457        Y=-Y        Y=-Y
# Line 463  c--------------------------------------- Line 463  c---------------------------------------
463        NMED=MAL        NMED=MAL
464        CALL GSVOLU('LEGB','BOX ',NMED,LEGB, 3,IVOLU)        CALL GSVOLU('LEGB','BOX ',NMED,LEGB, 3,IVOLU)
465  *  *
466  *     empty volumes to be subtracted from LEGB        *     empty volumes to be subtracted from LEGB
467  *  *
468        NMED=MN2        NMED=MN2
469        CALL GSVOLU('LEG1','TRAP',NMED,LEG1, 11,IVOLU)                    CALL GSVOLU('LEG1','TRAP',NMED,LEG1, 11,IVOLU)
470        X=-LEGB(1)+LEG1(1)        X=-LEGB(1)+LEG1(1)
471        Y=-LEGB(2)+LEG1(4)        Y=-LEGB(2)+LEG1(4)
472        Z=0.        Z=0.
473        CALL GSPOS('LEG1',1,'LEGB',X,Y,Z,RY90,'ONLY')              CALL GSPOS('LEG1',1,'LEGB',X,Y,Z,RY90,'ONLY')
474        Y=-LEGB(2)+2.*LEGP(1)-LEG1(4)        Y=-LEGB(2)+2.*LEGP(1)-LEG1(4)
475        CALL GSPOS('LEG1',2,'LEGB',X,Y,Z,RY90,'ONLY')              CALL GSPOS('LEG1',2,'LEGB',X,Y,Z,RY90,'ONLY')
476  *        *
477        CALL GSVOLU('LEG2','BOX ',NMED,LEG2, 3,IVOLU)              CALL GSVOLU('LEG2','BOX ',NMED,LEG2, 3,IVOLU)
478        X=-LEGB(1)+LEG2(1)        X=-LEGB(1)+LEG2(1)
479        Y=-LEGB(2)+LEGP(1)        Y=-LEGB(2)+LEGP(1)
480        Z=0.        Z=0.
481        CALL GSPOS('LEG2',1,'LEGB',X,Y,Z,0,'ONLY')        CALL GSPOS('LEG2',1,'LEGB',X,Y,Z,0,'ONLY')
482  *  *
483        CALL GSVOLU('LEG3','BOX ',NMED,LEG3, 3,IVOLU)              CALL GSVOLU('LEG3','BOX ',NMED,LEG3, 3,IVOLU)
484        X=-LEGB(1)+LEGP(2)*2.+LEG3(1)+0.2        X=-LEGB(1)+LEGP(2)*2.+LEG3(1)+0.2
485        Y=LEGB(2)-LEG3(2)        Y=LEGB(2)-LEG3(2)
486        Z=LEGB(3)-LEG3(3)        Z=LEGB(3)-LEG3(3)
487        CALL GSPOS('LEG3',1,'LEGB',X,Y,Z,0,'ONLY')        CALL GSPOS('LEG3',1,'LEGB',X,Y,Z,0,'ONLY')
488        Z=-LEGB(3)+LEG3(3)        Z=-LEGB(3)+LEG3(3)
489        CALL GSPOS('LEG3',2,'LEGB',X,Y,Z,0,'ONLY')              CALL GSPOS('LEG3',2,'LEGB',X,Y,Z,0,'ONLY')
490  *  *
491        CALL GSVOLU('LEG4','BOX ',NMED,LEG4, 3,IVOLU)              CALL GSVOLU('LEG4','BOX ',NMED,LEG4, 3,IVOLU)
492        X=-LEGB(1)+LEG4(1)        X=-LEGB(1)+LEG4(1)
493        Y=LEGB(2)-LEG4(2)        Y=LEGB(2)-LEG4(2)
494        Z=0.        Z=0.
495        CALL GSPOS('LEG4',1,'LEGB',X,Y,Z,0,'ONLY')        CALL GSPOS('LEG4',1,'LEGB',X,Y,Z,0,'ONLY')
496  *        *
497        CALL GSVOLU('LEG5','BOX ',NMED,LEG5, 3,IVOLU)              CALL GSVOLU('LEG5','BOX ',NMED,LEG5, 3,IVOLU)
498        X=LEGB(1)-LEG5(1)        X=LEGB(1)-LEG5(1)
499        Y=LEGB(2)-LEG5(2)        Y=LEGB(2)-LEG5(2)
500        Z=0.        Z=0.
501        CALL GSPOS('LEG5',1,'LEGB',X,Y,Z,0,'ONLY')        CALL GSPOS('LEG5',1,'LEGB',X,Y,Z,0,'ONLY')
502  *  *
503        CALL GSVOLU('LEG6','BOX ',NMED,LEG6, 3,IVOLU)              CALL GSVOLU('LEG6','BOX ',NMED,LEG6, 3,IVOLU)
504        X=-LEGB(1)+LEGP(2)*2.+LEG3(1)        X=-LEGB(1)+LEGP(2)*2.+LEG3(1)
505        Y=LEGB(2)-LEG6(2)        Y=LEGB(2)-LEG6(2)
506        Z=0.        Z=0.
507        CALL GSPOS('LEG6',1,'LEGB',X,Y,Z,0,'ONLY')              CALL GSPOS('LEG6',1,'LEGB',X,Y,Z,0,'ONLY')
508  *  *
509        CALL GSVOLU('LEG7','BOX ',NMED,LEG7, 3,IVOLU)              CALL GSVOLU('LEG7','BOX ',NMED,LEG7, 3,IVOLU)
510        X=LEGB(1)-LEGP(3)-LEGP(4)        X=LEGB(1)-LEGP(3)-LEGP(4)
511        Y=-LEGB(2)+LEG7(2)        Y=-LEGB(2)+LEG7(2)
512        Z=LEGP(5)-LEGP(8)-LEG7(3)        Z=LEGP(5)-LEGP(8)-LEG7(3)
513        CALL GSPOS('LEG7',1,'LEGB',X,Y,Z,0,'ONLY')              CALL GSPOS('LEG7',1,'LEGB',X,Y,Z,0,'ONLY')
514  *  *
515        CALL GSVOLU('LEG8','BOX ',NMED,LEG8, 3,IVOLU)              CALL GSVOLU('LEG8','BOX ',NMED,LEG8, 3,IVOLU)
516        Z=0.        Z=0.
517        CALL GSPOS('LEG8',1,'LEGB',X,Y,Z,0,'ONLY')              CALL GSPOS('LEG8',1,'LEGB',X,Y,Z,0,'ONLY')
518  *  *
519        CALL GSVOLU('LEG9','BOX ',NMED,LEG9, 3,IVOLU)              CALL GSVOLU('LEG9','BOX ',NMED,LEG9, 3,IVOLU)
520        Z=-(LEGP(5)-LEGP(8)-LEG7(3))        Z=-(LEGP(5)-LEGP(8)-LEG7(3))
521        CALL GSPOS('LEG9',1,'LEGB',X,Y,Z,0,'ONLY')              CALL GSPOS('LEG9',1,'LEGB',X,Y,Z,0,'ONLY')
522  *  *
523  *     Positioning LEG BOXES in CARD  *     Positioning LEG BOXES in CARD
524  *  *
# Line 558  c--------------------------------------- Line 558  c---------------------------------------
558  *  *
559        CALL GSVOLU('TH11','BOX ',NMED,TH11, 3,IVOLU)        CALL GSVOLU('TH11','BOX ',NMED,TH11, 3,IVOLU)
560        CALL GSVOLU('TH12','BOX ',NMED,TH12, 3,IVOLU)        CALL GSVOLU('TH12','BOX ',NMED,TH12, 3,IVOLU)
561        CALL GSVOLU('TH21','BOX ',NMED,TH21, 3,IVOLU)              CALL GSVOLU('TH21','BOX ',NMED,TH21, 3,IVOLU)
562        CALL GSVOLU('TH22','BOX ',NMED,TH22, 3,IVOLU)                              CALL GSVOLU('TH22','BOX ',NMED,TH22, 3,IVOLU)
563  *  *
564        CALL GSVOLU('TOPT','TRAP',NMED,TOPT,11,IVOLU)        CALL GSVOLU('TOPT','TRAP',NMED,TOPT,11,IVOLU)
565  *  *
# Line 573  c--------------------------------------- Line 573  c---------------------------------------
573        N=1        N=1
574        X=0.        X=0.
575        Y=0.        Y=0.
576        Z=0.            Z=0.
577        CALL GSPOS('TOPC',N,'TPLT',X,Y,Z,0,'ONLY')              CALL GSPOS('TOPC',N,'TPLT',X,Y,Z,0,'ONLY')
578  *  *
579  *  *
580  *     Positioning H11 in TPLT  *     Positioning H11 in TPLT
# Line 583  c--------------------------------------- Line 583  c---------------------------------------
583        X=-TOPP(1)        X=-TOPP(1)
584        Y=TOPP(2)        Y=TOPP(2)
585        Z=0.        Z=0.
586        CALL GSPOS('TH11',N,'TPLT',X,Y,Z,0,'ONLY')                        CALL GSPOS('TH11',N,'TPLT',X,Y,Z,0,'ONLY')
587        N=2        N=2
588        X=-TOPP(1)        X=-TOPP(1)
589        Y=-TOPP(2)        Y=-TOPP(2)
590        CALL GSPOS('TH11',N,'TPLT',X,Y,Z,0,'ONLY')                        CALL GSPOS('TH11',N,'TPLT',X,Y,Z,0,'ONLY')
591        N=3        N=3
592        X=TOPP(1)        X=TOPP(1)
593        Y=TOPP(2)        Y=TOPP(2)
594        CALL GSPOS('TH11',N,'TPLT',X,Y,Z,0,'ONLY')                        CALL GSPOS('TH11',N,'TPLT',X,Y,Z,0,'ONLY')
595        N=4        N=4
596        X=TOPP(1)        X=TOPP(1)
597        Y=-TOPP(2)        Y=-TOPP(2)
598        CALL GSPOS('TH11',N,'TPLT',X,Y,Z,0,'ONLY')                        CALL GSPOS('TH11',N,'TPLT',X,Y,Z,0,'ONLY')
599  *  *
600  *     Positioning H12 in TPLT  *     Positioning H12 in TPLT
601  *  *
# Line 603  c--------------------------------------- Line 603  c---------------------------------------
603        X=-TOPP(1)        X=-TOPP(1)
604        Y=TOPP(3)        Y=TOPP(3)
605        Z=0.        Z=0.
606        CALL GSPOS('TH12',N,'TPLT',X,Y,Z,0,'ONLY')                        CALL GSPOS('TH12',N,'TPLT',X,Y,Z,0,'ONLY')
607        N=2        N=2
608        X=-TOPP(1)        X=-TOPP(1)
609        Y=-TOPP(3)        Y=-TOPP(3)
610        CALL GSPOS('TH12',N,'TPLT',X,Y,Z,0,'ONLY')                        CALL GSPOS('TH12',N,'TPLT',X,Y,Z,0,'ONLY')
611        N=3        N=3
612        X=TOPP(1)        X=TOPP(1)
613        Y=TOPP(3)        Y=TOPP(3)
614        CALL GSPOS('TH12',N,'TPLT',X,Y,Z,0,'ONLY')                        CALL GSPOS('TH12',N,'TPLT',X,Y,Z,0,'ONLY')
615        N=4        N=4
616        X=TOPP(1)        X=TOPP(1)
617        Y=-TOPP(3)        Y=-TOPP(3)
618        CALL GSPOS('TH12',N,'TPLT',X,Y,Z,0,'ONLY')                        CALL GSPOS('TH12',N,'TPLT',X,Y,Z,0,'ONLY')
619  *  *
620  *     Positioning H21 in TPLT  *     Positioning H21 in TPLT
621  *  *
# Line 623  c--------------------------------------- Line 623  c---------------------------------------
623        X=-TOPP(4)        X=-TOPP(4)
624        Y=TOPP(5)        Y=TOPP(5)
625        Z=0.        Z=0.
626        CALL GSPOS('TH21',N,'TPLT',X,Y,Z,0,'ONLY')                        CALL GSPOS('TH21',N,'TPLT',X,Y,Z,0,'ONLY')
627        N=2        N=2
628        X=-TOPP(4)        X=-TOPP(4)
629        Y=-TOPP(5)        Y=-TOPP(5)
630        CALL GSPOS('TH21',N,'TPLT',X,Y,Z,0,'ONLY')                        CALL GSPOS('TH21',N,'TPLT',X,Y,Z,0,'ONLY')
631        N=3        N=3
632        X=TOPP(4)        X=TOPP(4)
633        Y=TOPP(5)        Y=TOPP(5)
634        CALL GSPOS('TH21',N,'TPLT',X,Y,Z,0,'ONLY')                        CALL GSPOS('TH21',N,'TPLT',X,Y,Z,0,'ONLY')
635        N=4        N=4
636        X=TOPP(4)        X=TOPP(4)
637        Y=-TOPP(5)        Y=-TOPP(5)
638        CALL GSPOS('TH21',N,'TPLT',X,Y,Z,0,'ONLY')                        CALL GSPOS('TH21',N,'TPLT',X,Y,Z,0,'ONLY')
639  *  *
640  *     Positioning H22 in TPLT  *     Positioning H22 in TPLT
641  *  *
# Line 643  c--------------------------------------- Line 643  c---------------------------------------
643        X=-TOPP(6)        X=-TOPP(6)
644        Y=TOPP(5)        Y=TOPP(5)
645        Z=0.        Z=0.
646        CALL GSPOS('TH22',N,'TPLT',X,Y,Z,0,'ONLY')                        CALL GSPOS('TH22',N,'TPLT',X,Y,Z,0,'ONLY')
647        N=2        N=2
648        X=-TOPP(6)        X=-TOPP(6)
649        Y=-TOPP(5)        Y=-TOPP(5)
650        CALL GSPOS('TH22',N,'TPLT',X,Y,Z,0,'ONLY')                        CALL GSPOS('TH22',N,'TPLT',X,Y,Z,0,'ONLY')
651        N=3        N=3
652        X=TOPP(6)        X=TOPP(6)
653        Y=TOPP(5)        Y=TOPP(5)
654        CALL GSPOS('TH22',N,'TPLT',X,Y,Z,0,'ONLY')                        CALL GSPOS('TH22',N,'TPLT',X,Y,Z,0,'ONLY')
655        N=4        N=4
656        X=TOPP(6)        X=TOPP(6)
657        Y=-TOPP(5)        Y=-TOPP(5)
658        CALL GSPOS('TH22',N,'TPLT',X,Y,Z,0,'ONLY')                        CALL GSPOS('TH22',N,'TPLT',X,Y,Z,0,'ONLY')
659        N=5        N=5
660        X=0.        X=0.
661        Y=TOPP(5)        Y=TOPP(5)
662        CALL GSPOS('TH22',N,'TPLT',X,Y,Z,0,'ONLY')                        CALL GSPOS('TH22',N,'TPLT',X,Y,Z,0,'ONLY')
663        N=6        N=6
664        X=0.        X=0.
665        Y=-TOPP(5)        Y=-TOPP(5)
666        CALL GSPOS('TH22',N,'TPLT',X,Y,Z,0,'ONLY')                        CALL GSPOS('TH22',N,'TPLT',X,Y,Z,0,'ONLY')
667  *  *
668  *     Positioning TOPT in TPLT  *     Positioning TOPT in TPLT
669  *  *
# Line 705  c--------------------------------------- Line 705  c---------------------------------------
705        CALL GSPOS('TOHO',N,'TPLT',X,Y,Z,0,'ONLY')        CALL GSPOS('TOHO',N,'TPLT',X,Y,Z,0,'ONLY')
706  *  *
707  *     Positioning TPLT in CARD volume  *     Positioning TPLT in CARD volume
708  *  *
709        N=1        N=1
710        X=0.        X=0.
711        Y=0.        Y=0.
712        Z=CARDB(3)-TPLT(3)        Z=CARDB(3)-TPLT(3)
713        CALL GSPOS('TPLT',N,'CARD',X,Y,Z,0,'ONLY')            CALL GSPOS('TPLT',N,'CARD',X,Y,Z,0,'ONLY')
714  C---------------------------------------------------------------        C---------------------------------------------------------------
715  *  *
716  *     MAGNETIC SCREEN  *     MAGNETIC SCREEN
717  *  *
718        NMED=10        NMED=MIRON
719        CALL GSVOLU('MGSC','BOX ',NMED,MGSC, 3,IVOLU)              CALL GSVOLU('MGSC','BOX ',NMED,MGSC, 3,IVOLU)
720  *  *
721  *  Empty volume to be subtracted from MGSC  *  Empty volume to be subtracted from MGSC
722  *  *
723        NMED=MN2        NMED=MN2
724        CALL GSVOLU('MGSH','BOX ',NMED,MGSH, 3,IVOLU)              CALL GSVOLU('MGSH','BOX ',NMED,MGSH, 3,IVOLU)
725    
726  *  *
727  *     Positioning in CARD...  *     Positioning in CARD...
728  *      *
729        N=1        N=1
730        X=0.        X=0.
731        Y=0.        Y=0.
732        Z=0.        Z=0.
733        CALL GSPOS('MGSH',N,'MGSC',X,Y,Z,0,'ONLY')        CALL GSPOS('MGSH',N,'MGSC',X,Y,Z,0,'ONLY')
734        Z=-CARDB(3)+MGSC(3)        Z=-CARDB(3)+MGSC(3)
735        CALL GSPOS('MGSC',N,'CARD',X,Y,Z,0,'ONLY')        CALL GSPOS('MGSC',N,'CARD',X,Y,Z,0,'ONLY')
736  c---------------------------------------------------------------        c---------------------------------------------------------------
737        RETURN        RETURN
738        END        END

Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4

  ViewVC Help
Powered by ViewVC 1.1.23