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

Diff of /gpamela/gpcal/gpcalv.F

Parent Directory | Revision Log | View Patch Patch

-revision 3.1 by cafagna,
Thu Jul 11 16:02:14 2002 UTC
+revision 3.6 by cafagna,
Fri Mar 30 15:24:54 2007 UTC
 Line 1
  *
- * $Id$
+ * $Id: gpcalv.F,v 3.5 2006/05/09 12:01:15 cafagna Exp $
+ *
+ * $Log: gpcalv.F,v $
+ * Revision 3.5  2006/05/09 12:01:15  cafagna
+ * Calo geometry updated
+ *
+ * Revision 3.4  2005/07/25 11:53:20  cafagna
+ * Several updates. See history for details
+ *
+ * Revision 3.3  2003/12/17 11:32:48  pamela
+ * CALO SIMULATION COMPLETED: geometry and special tracking parameters updated and simulation checked by a comparison with the Trieste's standalone Monte Carlo simulation
+ *
+ * Revision 3.2  2002/12/05 10:17:41  pamela
+ * Update CAS and CALO geometries and positions. Makefile updated as well
+ *
+ * Revision 3.1.1.1  2002/07/11 16:02:14  cafagna
+ * First GPAMELA release on CVS
  *
- * $Log$
  *
  *CMZ :  3.00/00 14/11/2000  10.17.50  by  Emiliano Mocchiutti
  *CMZ :  2.02/00 12/10/2000  18.19.42  by  Francesco Cafagna
-Line 24
+Line 39
  *
        INTEGER IROT,IVOLU,N,NMED
        REAL X,Y,Z
+ * MA
+       REAL X1, Y1
+       INTEGER m
+ * END MA
  *
  *  Define the CALB     volume
  *
-Line 39
+Line 58
  *
        NMED= MN2
        CALL GSVOLU('CAPL','BOX ',NMED,CAPL, 3,IVOLU)
+ * F.V.&ML:
+ *
+ *  Define the CANS     volume
+ *
+       NMED= MN2
+       CALL GSVOLU('CANS','BOX ',NMED,CANS, 3,IVOLU)
+ * END  F.V.&ML.
  *
  *  Define the CAPD     volume
  *
-Line 52
+Line 78
  *
  *  Division in CASI
  *
-       CALL GSDVN ('CAST','CASI',NCASTR,1)
+ * ML:
+ C      CALL GSDVN ('CAST','CASI',NCASTR,1)
+       CALL GSDVN ('CAST','CASI',NCASTR,2)
+ * END ML:
  *
  *  Define the CG10     volume
  *
-Line 93
+Line 122
        X = 0.
        Y = 0.
  *EM:
-       Z = -CAPL(3) + CAKA(3)
+ *F.V.&ML:
-       CALL GSPOS('CAKA',N,'CAPL',X,Y,Z,0,'ONLY')
+ c      CALL GPVOLU(-2)
-       CALL GSPOS('CAKA',N,'CAPD',X,Y,Z,0,'ONLY')
+ C      Z = -CAPL(3) + CAKA(3)
+ C      CALL GSPOS('CAKA',N,'CAPL',X,Y,Z,0,'ONLY')
+       Z = -CANS(3) + CAKA(3)
+       CALL GSPOS('CAKA',N,'CANS',X,Y,Z,0,'ONLY')
+ C      CALL GSPOS('CAKA',N,'CAPD',X,Y,Z,0,'ONLY')
        Z = Z + CAKA(3) + C10C(3)
-       CALL GSPOS('C10C',N,'CAPL',X,Y,Z,0,'ONLY')
+ C      CALL GSPOS('C10C',N,'CAPL',X,Y,Z,0,'ONLY')
-       CALL GSPOS('C10C',N,'CAPD',X,Y,Z,0,'ONLY')
+       CALL GSPOS('C10C',N,'CANS',X,Y,Z,0,'ONLY')
+ C      CALL GSPOS('C10C',N,'CAPD',X,Y,Z,0,'ONLY')
        Z = Z + C10C(3) + CAGL(3)
-       CALL GSPOS('CAGL',N,'CAPL',X,Y,Z,0,'ONLY')
+ C      CALL GSPOS('CAGL',N,'CAPL',X,Y,Z,0,'ONLY')
+       CALL GSPOS('CAGL',N,'CANS',X,Y,Z,0,'ONLY')
        Z = Z + CAGL(3) + CAKP(3)
-       CALL GSPOS('CAKP',N,'CAPL',X,Y,Z,0,'ONLY')
+ C      CALL GSPOS('CAKP',N,'CAPL',X,Y,Z,0,'ONLY')
+       CALL GSPOS('CAKP',N,'CANS',X,Y,Z,0,'ONLY')
+ * END F.V.&ML
  *
  * Positioning volume CASI into the mother CAPL
  *
        N = 0
-       Z = Z + CAKP(3) + CASI(3)
+ * F.V.&ML:
+ C      Z= Z + CAKP(3) + CASI(3)
+       Z = -CAPL(3)+CASI(3)
+ * END F.V.&ML
        DO I=-1,1
           Y = I*(2.*CASI(2)+ DCASIY)
           DO II=-1,1
-Line 121
+Line 161
           ENDDO
        ENDDO
  *END EM.
- *
+ * MA
- * Positioning volume CAPL and CAAB into the mother CALB
+ * 11 MODULI DI 2 W
- *
+ c     Z = CALB(3) - CAPL(3) - DCASIZ/2
-       X= 0.
+       Z = CALB(3) - CAPL(3) - DCASIZ
-       Y= 0.
+       N = 0
- C      Z= CALB(3) + CAPL(3)
+       DO M = 0,20,2
- *EM:
+ *     Positioning CAPL CANS CAAB in CALB to form a module
-       Z= CALB(3) + CAPL(3) - DCASIZ/2.
+          X = .10
- *END EM. (we have 0.15 cm from the top to the first Si)
+          y = .05
-       N= 0
+          DO I = 1,2
-       DO I=1,NCAPL
+ *     Positioning CAPL
- *
+             N = N + 1
- * X view, absorber and Y view on the other
+ C     if(N.EQ.1) print *,'z ini = ',Z
- *
+             IF(.NOT.NOCAPL(N))
-          Z= Z - 2*CAPL(3)
+      +           CALL GSPOS('CAPL',N,'CALB',X,(Y-0.1),Z,0,'ONLY')
-          N= N + 1
+ *     Positioning CANS
- *
+             Z = Z - CAPL(3) - CANS(3)
- * Check if this plane is excluded
+             CALL GSPOS('CANS',N,'CALB',0,-0.1,Z,0,'ONLY')
- *
+ *     Positioning CAAB
-          IF(.NOT.NOCAPL(N))
+             Z = Z - CANS(3) - CAAB(3)
-      +        CALL GSPOS('CAPL',N,'CALB',X,Y,Z,0,'ONLY')
+             IF(.NOT.NOCAAB(M+I))
-          Z= Z - CAPL(3) - CAAB(3)
+      +           CALL GSPOS('CAAB',(M+I),'CALB',0,-0.1,Z,0,'ONLY')
- *
+ *     Positioning CANS (ruotato)
- * Check if this absorber is excluded
+             N = N + 1
- *
+             Z = Z - CAAB(3) - CANS(3)
-          IF(.NOT.NOCAAB(I))
+             CALL GSPOS('CANS',N,'CALB',0,-0.1,Z,3,'ONLY')
-      +        CALL GSPOS('CAAB',(I-1),'CALB',X,Y,Z,0,'ONLY')
+ *     Positioning CAPL (ruotato e shiftato)
-          Z= Z - CAAB(3) - CAPL(3)
+             X1 = X
-          N= N + 1
+             Y1 = Y
- *
+             X = -Y1
- * Check if this detector is excluded
+             Y = X1
- *
+             Z = Z - CANS(3) - CAPL(3)
-          IF(.NOT.NOCAPL(N))
+ c Mod by Caf on 30th Marc 2007, the X position of the X-view planes should be swapped
-      +        CALL GSPOS('CAPL',N,'CALB',X,Y,Z,3,'ONLY')
+ c$$$            IF(.NOT.NOCAPL(N))
-       ENDDO
+ c$$$     +           CALL GSPOS('CAPL',N,'CALB',X,(Y-0.1),Z,3,'ONLY')
- *
+             IF(.NOT.NOCAPL(N))
- * Positioning volume CAPD and CAAD into the mother CALB
+      +           CALL GSPOS('CAPL',N,'CALB',-X,(Y-0.1),Z,3,'ONLY')
- *
-       N= 0
+             IF(I.EQ.1) THEN
-       DO I=1,NCAPLD
+ c              Z = Z - 2*CAPL(3) - DZMO ! SPESSORE TRA I SILICI CHE SI GUARDANO
-          IF(I.EQ.1) THEN
+                Z = Z - 2*CAPL(3) - DZM0 ! corretto
-             Z= Z - CAPD(3) - CAPL(3)
+                X1 = X                   ! SHIFT PER IL CAPL PROSSIMO CICLO
-          ELSE
+                Y1 = Y
-             Z= Z - 2*CAPD(3)
+                X = -Y1
-          ENDIF
+                Y = X1
-          N= N + 1
+             ENDIF
-          CALL GSPOS('CAPD',N,'CALB',X,Y,Z,0,'ONLY')
+          ENDDO
-          Z= Z - CAPD(3) - CAAD(3)
+          Z = Z - DZST - 2*CAPL(3)
-          CALL GSPOS('CAAD',(I-1),'CALB',X,Y,Z,0,'ONLY')
-          Z= Z - CAAD(3) - CAPD(3)
-          N= N + 1
-          CALL GSPOS('CAPD',N,'CALB',X,Y,Z,3,'ONLY')
        ENDDO
+ *     END MA
+ C # *F.V.&ML:
+ C # C commenting the original code with a 'C' per raw:
+ C # C*
+ C # C* Positioning volume CAPL and CAAB into the mother CALB
+ C # C*
+ C # C      X= 0.
+ C # C      Y= 0.
+ C # CC      Z= CALB(3) + CAPL(3)
+ C # C*EM:
+ C # C      Z= CALB(3) + CAPL(3) - DCASIZ/2.
+ C # C*END EM. (we have 0.15 cm from the top to the first Si)
+ C # C      N= 0
+ C # C      DO I=1,NCAPL
+ C # C*
+ C # C* X view, absorber and Y view on the other
+ C # C*
+ C # C         Z= Z - 2*CAPL(3)
+ C # C         N= N + 1
+ C # C*
+ C # C* Check if this plane is excluded
+ C # C*
+ C # C         IF(.NOT.NOCAPL(N))
+ C # C     +        CALL GSPOS('CAPL',N,'CALB',X,Y,Z,0,'ONLY')
+ C # C         Z= Z - CAPL(3) - CAAB(3)
+ C # C*
+ C # C* Check if this absorber is excluded
+ C # C*
+ C # C         IF(.NOT.NOCAAB(I))
+ C # C     +        CALL GSPOS('CAAB',(I-1),'CALB',X,Y,Z,0,'ONLY')
+ C # C         Z= Z - CAAB(3) - CAPL(3)
+ C # C         N= N + 1
+ C # C*
+ C # C* Check if this detector is excluded
+ C # C*
+ C # C         IF(.NOT.NOCAPL(N))
+ C # C     +        CALL GSPOS('CAPL',N,'CALB',X,Y,Z,3,'ONLY')
+ C # C      ENDDO
+ C # C end comments.
+ C # *
+ C # * Positioning volume CANS ,CAPL and CAAB into the mother CALB
+ C # *
+ C #       X= 0.
+ C #       Y= 0.
+ C #       Z= CALB(3) + CAPL(3) - DCASIZ/2.
+ C # *END EM. (we have 0.15 cm from the top to the first Si)
+ C #       N= 0
+ C # *
+ C # *   I DISPARI
+ C #       DO I=1,NCAPL,2
+ C # *
+ C # * X view, absorber and Y view on the other
+ C # *
+ C #      X=0.10
+ C #      Y=0.05
+ C #          Z= Z - 2*CAPL(3)
+ C #          N= N + 1
+ C # *
+ C # * Check if this plane is excluded
+ C # *
+ C # *
+ C #          IF(.NOT.NOCAPL(N))
+ C #      +        CALL GSPOS('CAPL',N,'CALB',X,Y,Z,0,'ONLY')
+ C # *
+ C #          X = 0.
+ C #          Y = 0.
+ C #          Z= Z - CAPL(3) - CANS(3)
+ C # *
+ C #          CALL GSPOS('CANS',N,'CALB',X,Y,Z,0,'ONLY')
+ C # *
+ C #          Z= Z - CANS(3) - CAAB(3)
+ C # *
+ C # * Check if this absorber is excluded
+ C # *
+ C # *
+ C #          IF(.NOT.NOCAAB(I))
+ C #      +        CALL GSPOS('CAAB',(I-1),'CALB',X,Y,Z,0,'ONLY')
+ C # *
+ C #          Z = Z - CAAB(3)- CANS(3)
+ C #          N = N + 1
+ C # *
+ C #          CALL GSPOS('CANS',N,'CALB',X,Y,Z,3,'ONLY')
+ C #          Z= Z - CANS(3) - CAPL(3)
+ C #          X=-0.05
+ C #          Y=0.10
+ C # *
+ C # * Check if this detector is excluded
+ C # *
+ C #          IF(.NOT.NOCAPL(N))
+ C #      +        CALL GSPOS('CAPL',N,'CALB',X,Y,Z,3,'ONLY')
+ C #      N=N+2
+ C #      Z = Z -2*(2*CAPL(3)+2*CANS(3)+CAAB(3))
+ C #       ENDDO
+ C # *
+ C # *   I PARI
+ C # *
+ C #       X= 0.
+ C #       Y= 0.
+ C #       Z= CALB(3) + CAPL(3) - DCASIZ/2.
+ C # *END EM. (we have 0.15 cm from the top to the first Si)
+ C #       N= 2
+ C # *
+ C # *
+ C #       DO I=2,NCAPL,2
+ C # *
+ C # * X view, absorber and Y view on the other
+ C # *
+ C #          Z= Z - 2*CAPL(3)-2*(2*CAPL(3)+2*CANS(3)+CAAB(3))
+ C #          N= N + 1
+ C # *
+ C # * Check if this plane is excluded
+ C #          X= -0.10
+ C #          Y= -0.05
+ C # *
+ C #          IF(.NOT.NOCAPL(N))
+ C #      +        CALL GSPOS('CAPL',N,'CALB',X,Y,Z,0,'ONLY')
+ C #      X = 0.
+ C #      Y = 0.
+ C #          Z = Z - CAPL(3) - CANS(3)
+ C # *
+ C #              CALL GSPOS('CANS',N,'CALB',X,Y,Z,0,'ONLY')
+ C # *
+ C # * Check if this absorber is excluded
+ C # *
+ C #          N = N + 1
+ C #          Z = Z - CANS(3) - CAAB(3)
+ C # *
+ C #          IF(.NOT.NOCAAB(I))
+ C #      +        CALL GSPOS('CAAB',(I-1),'CALB',X,Y,Z,0,'ONLY')
+ C # *
+ C #      Z= Z - CAAB(3) - CANS(3)
+ C # *
+ C #          CALL GSPOS('CANS',N,'CALB',X,Y,Z,3,'ONLY')
+ C #          Z= Z - CANS(3) - CAPL(3)
+ C #          X= 0.05
+ C #          Y= -0.10
+ C # *
+ C # *
+ C # * Check if this detector is excluded
+ C # *
+ C #          IF(.NOT.NOCAPL(N))
+ C #      +        CALL GSPOS('CAPL',N,'CALB',X,Y,Z,3,'ONLY')
+ C #          N = N + 2
+ C # *
+ C #       ENDDO
+ C # *END F.V.&ML.
+ C # * F.V.&ML.:
+ C # * again we have to comment with a 'C' per raw:
+ C # C*
+ C # C* Positioning volume CAPD and CAAD into the mother CALB
+ C # C*
+ C # C      N= 0
+ C # C      DO I=1,NCAPLD
+ C # C         IF(I.EQ.1) THEN
+ C # C            Z= Z - CAPD(3) - CAPL(3)
+ C # C         ELSE
+ C # C            Z= Z - 2*CAPD(3)
+ C # C         ENDIF
+ C # C         N= N + 1
+ C # C         CALL GSPOS('CAPD',N,'CALB',X,Y,Z,0,'ONLY')
+ C # C         Z= Z - CAPD(3) - CAAD(3)
+ C # C         CALL GSPOS('CAAD',(I-1),'CALB',X,Y,Z,0,'ONLY')
+ C # C         Z= Z - CAAD(3) - CAPD(3)
+ C # C         N= N + 1
+ C # C         CALL GSPOS('CAPD',N,'CALB',X,Y,Z,3,'ONLY')
+ C # C      ENDDO
+ C # * end comments.
+ C # *END F.V.&ML.
+ C # end MA comments
  *
  * Put the shell around it
  *

 Legend:



Removed from v.3.1
 


changed lines


 
Added in v.3.6
 Legend:



Removed from v.3.1
 


changed lines


 
Added in v.3.6
-Removed from v.3.1
+Added in v.3.6

	ViewVC Help
Powered by ViewVC 1.1.23