/[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.2 by mocchiut, Tue May 30 14:57:59 2006 UTC revision 1.6 by mocchiut, Thu Aug 10 06:32:02 2006 UTC
# Line 142  C Line 142  C
142              call DOTRACK(NPLA,ZIN,XOUT,YOUT,AL_P,IFAIL)              call DOTRACK(NPLA,ZIN,XOUT,YOUT,AL_P,IFAIL)
143              if(IFAIL.ne.0)then              if(IFAIL.ne.0)then
144                 GOOD2 = 0                 GOOD2 = 0
145                 print *,' CALORIMETER - WARNING F77: tracking failed '  c               print *,' CALORIMETER - WARNING F77: tracking failed '
146                 goto 969                 goto 969
147              endif              endif
148              TX = TAN(ASIN(AL_P(3))) * COS(AL_P(4))              TX = TAN(ASIN(AL_P(3))) * COS(AL_P(4))
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 159  C     Line 159  C    
159  C     CHECK IF XOUT OR YOUT ARE NaN  C     CHECK IF XOUT OR YOUT ARE NaN
160  C      C    
161                 IF (XOUT(I).NE.XOUT(I).OR.YOUT(I).NE.YOUT(I)) THEN                 IF (XOUT(I).NE.XOUT(I).OR.YOUT(I).NE.YOUT(I)) THEN
162                    print *,  c                  print *,
163       &         ' CALORIMETER - WARNING F77: tracking error (NaN values)'                    c     &         ' CALORIMETER - WARNING F77: tracking error (NaN values)'                  
164                    GOOD2 = 0                    GOOD2 = 0
165                    GOTO 969                    GOTO 969
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
# Line 306  C IF THE TRACK IS OUTSIDE THE CALORIMETE Line 306  C IF THE TRACK IS OUTSIDE THE CALORIMETE
306  C  C
307        IF (CHTRACK.EQ.44) THEN        IF (CHTRACK.EQ.44) THEN
308           GOOD2 = 0           GOOD2 = 0
309           PRINT *,' CALORIMETER - WARNING F77: track outside calorimeter'  c         PRINT *,' CALORIMETER - WARNING F77: track outside calorimeter'
310           GOTO 50           GOTO 50
311        ELSE        ELSE
312           IF ( IWPL(1).LE.IWPL(2) ) THEN           IF ( IWPL(1).LE.IWPL(2) ) THEN
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.23