/[PAMELA software]/DarthVader/CalorimeterLevel2/src/calol2cm.for
ViewVC logotype

Diff of /DarthVader/CalorimeterLevel2/src/calol2cm.for

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

revision 1.6 by mocchiut, Mon Mar 26 14:02:07 2007 UTC revision 1.7 by mocchiut, Fri Mar 30 11:17:16 2007 UTC
# Line 106  C Line 106  C
106           DO J = 1,96           DO J = 1,96
107              IF ( MOD(I,2).NE.0 ) THEN              IF ( MOD(I,2).NE.0 ) THEN
108                 IF ( ESTRIP(2,I,J).GT.EMIN ) THEN                 IF ( ESTRIP(2,I,J).GT.EMIN ) THEN
                   ENPLANEY = ENPLANEY + ESTRIP(2,I,J)  
109                    DEXY(2,I,J) = ESTRIP(2,I,J)                    DEXY(2,I,J) = ESTRIP(2,I,J)
110                      ENPLANEY = ENPLANEY + DEXY(2,I,J)
111                    NSTRIP = NSTRIP + 1.                    NSTRIP = NSTRIP + 1.
112                    QTOT = QTOT + ESTRIP(2,I,J)                    QTOT = QTOT + DEXY(2,I,J)
113                    IF (I.LT.11) QQ(1) = QQ(1) + ESTRIP(2,I,J)                    IF (I.LT.11) QQ(1) = QQ(1) + DEXY(2,I,J)
114                 ENDIF                 ENDIF
115                 IF ( ESTRIP(1,I,J).GT.EMIN ) THEN                 IF ( ESTRIP(1,I,J).GT.EMIN ) THEN
116                    ENPLANEX = ENPLANEX + ESTRIP(1,I,J)  C               IF ( ESTRIP(1,I+1,J).GT.EMIN ) THEN
117    C                  
118    c                  print *,' i ',i,' i + 1 ',I+1
119    C                  DEXY(1,I,J) = ESTRIP(1,I+1,J)
120                    DEXY(1,I,J) = ESTRIP(1,I,J)                    DEXY(1,I,J) = ESTRIP(1,I,J)
121    C
122                      ENPLANEX = ENPLANEX + DEXY(1,I,J)
123                    NSTRIP = NSTRIP + 1.                    NSTRIP = NSTRIP + 1.
124                    QTOT = QTOT + ESTRIP(1,I,J)                    QTOT = QTOT + DEXY(1,I,J)
125                    if (i.lt.11) QQ(2) = QQ(2) + ESTRIP(1,I,J)                    if (i.lt.11) QQ(2) = QQ(2) + DEXY(1,I,J)
126    c                  IF (I.EQ.21) THEN
127    c                     NX22 = NX22 + 1.
128    c                     QX22 = QX22 + DEXY(1,I,J)
129    c                  ENDIF
130                 ENDIF                 ENDIF
131              ENDIF              ENDIF
132              IF ( MOD(I,2).EQ.0 ) THEN              IF ( MOD(I,2).EQ.0 ) THEN
133                 IF (ESTRIP(2,I,J).GT.EMIN) THEN                 IF (ESTRIP(2,I,J).GT.EMIN) THEN
                   ENPLANEY = ENPLANEY + ESTRIP(2,I,J)  
134                    DEXY(2,I,J) = ESTRIP(2,I,J)                    DEXY(2,I,J) = ESTRIP(2,I,J)
135                      ENPLANEY = ENPLANEY + DEXY(2,I,J)
136                    NSTRIP = NSTRIP + 1.                    NSTRIP = NSTRIP + 1.
137                    QTOT = QTOT + ESTRIP(2,I,J)                    QTOT = QTOT + DEXY(2,I,J)
138                    if (i.lt.11) QQ(3) = QQ(3) + ESTRIP(2,I,J)                    if (i.lt.11) QQ(3) = QQ(3) + DEXY(2,I,J)
139                 ENDIF                 ENDIF
140                 IF (ESTRIP(1,I,J).GT.EMIN) THEN                 IF (ESTRIP(1,I,J).GT.EMIN) THEN
141                    ENPLANEX = ENPLANEX + ESTRIP(1,I,J)  C               IF (ESTRIP(1,I-1,J).GT.EMIN) THEN
142    C                  
143    c                  print *,' i ',i,' i - 1 ',I-1
144    C                  DEXY(1,I,J) = ESTRIP(1,I-1,J)
145                    DEXY(1,I,J) = ESTRIP(1,I,J)                    DEXY(1,I,J) = ESTRIP(1,I,J)
146    C
147                      ENPLANEX = ENPLANEX + DEXY(1,I,J)
148                    NSTRIP = NSTRIP + 1.                    NSTRIP = NSTRIP + 1.
149                    QTOT = QTOT + ESTRIP(1,I,J)                    QTOT = QTOT + DEXY(1,I,J)
150                    IF (I.EQ.22) THEN                    IF (I.EQ.22) THEN
151                       NX22 = NX22 + 1.                       NX22 = NX22 + 1.
152                       QX22 = QX22 + ESTRIP(1,I,J)                       QX22 = QX22 + DEXY(1,I,J)
153                    ENDIF                    ENDIF
154                    IF (I.LT.11) QQ(4) = QQ(4) + ESTRIP(1,I,J)                    IF (I.LT.11) QQ(4) = QQ(4) + DEXY(1,I,J)
155                 ENDIF                 ENDIF
156              ENDIF              ENDIF
157  C  C
# Line 176  C Line 190  C
190                 ENDIF                 ENDIF
191  C                C              
192                 DO I = 1,NPLA                     DO I = 1,NPLA    
193                    NN = 0                    NN = 1
194                    IF (M.EQ.2) NN = 1  C                  NN = 0
195    C                  IF (M.EQ.2) NN = 1
196                    IF (MOD(I,2).EQ.NN) THEN                    IF (MOD(I,2).EQ.NN) THEN
197                       SHIFT = +0.5                       SHIFT = +0.5
198                    ELSE                    ELSE
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.23