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

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

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

revision 1.3 by mocchiut, Wed May 31 09:31:11 2006 UTC revision 1.6 by mocchiut, Thu Aug 10 06:32:02 2006 UTC
# Line 121  C Line 121  C
121              al_p(m) = al_pp(t,m)              al_p(m) = al_pp(t,m)
122           enddo           enddo
123           if (al_p(5).eq.0.) THEN           if (al_p(5).eq.0.) THEN
124              PRINT *,' CALORIMETER - WARNING F77: track with R = 0, discarded'         PRINT *,' CALORIMETER - WARNING F77: track with R = 0, discarded'
125              GOOD2 = 0              GOOD2 = 0
126              GOTO 969              GOTO 969
127           ENDIF           ENDIF
# Line 149  c               print *,' CALORIMETER - Line 149  c               print *,' CALORIMETER -
149              TY = TAN(ASIN(AL_P(3))) * SIN(AL_P(4))              TY = TAN(ASIN(AL_P(3))) * SIN(AL_P(4))
150              DO I = 1, NPLA              DO I = 1, NPLA
151                 NN = 0                 NN = 0
152                 IF (M.EQ.2) NN = 1  C?????               IF (M.EQ.2) NN = 1
153                 IF (MOD(I,2).EQ.NN) THEN                 IF (MOD(I,2).EQ.NN) THEN
154                    SHIFT = +0.5                    SHIFT = +0.5
155                 ELSE                 ELSE
# Line 166  c     &         ' CALORIMETER - WARNING Line 166  c     &         ' CALORIMETER - WARNING
166                 ENDIF                 ENDIF
167  C  C
168                 CX = XOUT(I)*10. + XALIG                 CX = XOUT(I)*10. + XALIG
169                 CY = -YOUT(I)*10. + YALIG                 CY = YOUT(I)*10. + YALIG
170  C      C    
171                 IF (I.EQ.1) THEN                 IF (I.EQ.1) THEN
172                    TIMPX = CX                    TIMPX = CX
# Line 180  C     Line 180  C    
180                 ELSE                 ELSE
181                    YY(I) = CY                    YY(I) = CY
182                    BAR(M,I) = YY(I)                                      BAR(M,I) = YY(I)                  
183                    TBAR(M,I) = (YALIG - YY(I))/10.                        TBAR(M,I) = (-YALIG + YY(I))/10.    
184                    IF (I.EQ.22) MY=ABS(Y(1)-Y(22))/ABS(ZIN(1)-ZIN(22))                    IF (I.EQ.22) MY=ABS(Y(1)-Y(22))/ABS(ZIN(1)-ZIN(22))
185                 ENDIF                 ENDIF
186                 CALL LASTRISCIA(BAR(M,I),IBAR(M,I))                 CALL LASTRISCIA(BAR(M,I),IBAR(M,I))
# Line 261  C Line 261  C
261  C  C
262                 DO I = 1,NPLA                     DO I = 1,NPLA    
263                    NN = 0                    NN = 0
264                    IF (M.EQ.2) NN = 1  C?????                  IF (M.EQ.2) NN = 1
265                    IF (MOD(I,2).EQ.NN) THEN                    IF (MOD(I,2).EQ.NN) THEN
266                       SHIFT = +0.5                       SHIFT = +0.5
267                    ELSE                    ELSE
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.23