/[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.1 by mocchiut, Fri May 19 13:15:50 2006 UTC revision 1.3 by mocchiut, Wed May 31 09:31:11 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 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))
# 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
# 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
# Line 559  C Line 559  C
559        DO J = NPLA-4,NPLA        DO J = NPLA-4,NPLA
560           NNX = IBAR(1,J)           NNX = IBAR(1,J)
561           IF (NNX.NE.-1) THEN           IF (NNX.NE.-1) THEN
562              IF (NNX.LT.9) NNX = 9              IF (NNX.LT.5) NNX = 5
563              IF (NNX.GT.88) NNX = 88              IF (NNX.GT.92) NNX = 92
564              INFX = NNX - 8  c            IF (NNX.LT.9) NNX = 9
565              ISUPX = NNX + 8  c            IF (NNX.GT.88) NNX = 88
566                INFX = NNX - 4
567                ISUPX = NNX + 4
568    c            INFX = NNX - 8
569    c            ISUPX = NNX + 8
570              DO I = INFX,ISUPX              DO I = INFX,ISUPX
571                 IF (DEXY(1,J,I).GE.EMIN) THEN                 IF (DEXY(1,J,I).GE.EMIN) THEN
572                    NLAST = NLAST + 1                    NLAST = NLAST + 1
# Line 573  C Line 577  C
577  C  C
578           NNY = IBAR(2,J)           NNY = IBAR(2,J)
579           IF (NNY.NE.-1) THEN           IF (NNY.NE.-1) THEN
580              IF (NNY.LT.9) NNY = 9              IF (NNY.LT.5) NNY = 5
581              IF (NNY.GT.88) NNY = 88              IF (NNY.GT.92) NNY = 92
582              INFY = NNY - 8  c            IF (NNY.LT.9) NNY = 9
583              ISUPY = NNY + 8  c            IF (NNY.GT.88) NNY = 88
584                INFY = NNY - 4
585                ISUPY = NNY + 4
586    c            INFY = NNY - 8
587    c            ISUPY = NNY + 8
588              DO I=INFY,ISUPY              DO I=INFY,ISUPY
589                 IF (DEXY(2,J,I).GE.EMIN) THEN                 IF (DEXY(2,J,I).GE.EMIN) THEN
590                    NLAST = NLAST + 1                    NLAST = NLAST + 1
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.23