/[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.23 by mocchiut, Wed May 23 12:28:33 2012 UTC revision 1.24 by mocchiut, Thu Jan 16 15:29:11 2014 UTC
# Line 24  C Line 24  C
24        INTEGER Nfitx,Nfity, MNPLA        INTEGER Nfitx,Nfity, MNPLA
25  C  C
26        INTEGER INDEX, NTOT(2), NPIANI, GTR        INTEGER INDEX, NTOT(2), NPIANI, GTR
27        INTEGER j, m, i, IWPL(2), timpx, timpy, T, nn  c      INTEGER j, m, i, IWPL(2), timpx, timpy, T, nn
28          INTEGER j, m, i, IWPL(2), T, nn
29        INTEGER IPLANE, NNX, NNY, INFX, INFY, ISUPX, ISUPY        INTEGER IPLANE, NNX, NNY, INFX, INFY, ISUPX, ISUPY
30        INTEGER IBAR(2,NPLAV), NPFIT(2), CHTRACK,IWPLU        INTEGER IBAR(2,NPLAV), NPFIT(2), CHTRACK,IWPLU
31        INTEGER Iquest(100), ICONTROL5, nin, IFAIL        INTEGER Iquest(100), ICONTROL5, nin, IFAIL
# Line 144  c               print *,' ZIN(',I,') = ' Line 145  c               print *,' ZIN(',I,') = '
145  c               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 = REAL(DTAN(DASIN(AL_P(3))) * DCOS(AL_P(4)))
149              TY = TAN(ASIN(AL_P(3))) * SIN(AL_P(4))              TY = REAL(DTAN(DASIN(AL_P(3))) * DSIN(AL_P(4)))
150              DO I = 1, NPLA              DO I = 1, NPLA
151                 NN = 0                 NN = 0
152  c               IF (M.EQ.2) NN = 1  c               IF (M.EQ.2) NN = 1
# Line 172  c     &         ' CALORIMETER - WARNING Line 173  c     &         ' CALORIMETER - WARNING
173                    GOTO 969                    GOTO 969
174                 ENDIF                 ENDIF
175  C  C
176                 CX = XOUT(I)*10. + XALIG                 CX = REAL(XOUT(I))*10. + XALIG
177                 CY = YOUT(I)*10. + YALIG                 CY = REAL(YOUT(I))*10. + YALIG
178  C      C    
179                 IF (I.EQ.1) THEN  c               IF (I.EQ.1) THEN   !EM GCC4.7  TIMPX/Y are not used in che code...
180                    TIMPX = CX  c                  TIMPX = NINT(CX)
181                    TIMPY = CY  c                  TIMPY = NINT(CY)
182                 ENDIF  c               ENDIF
183                 IF (M.EQ.1) THEN                 IF (M.EQ.1) THEN
184                    Y(I) = CX                    Y(I) = CX
185                    BAR(M,I) = Y(I)                        BAR(M,I) = Y(I)    
186                    TBAR(M,I) = (Y(I) - XALIG)/10.                    TBAR(M,I) = (Y(I) - XALIG)/10.
187                    IF (I.EQ.NPLA) MX=ABS(Y(1)-Y(NPLA))/                    IF (I.EQ.NPLA) MX=ABS(Y(1)-Y(NPLA))/
188       &                 ABS(ZIN(1)-ZIN(NPLA))       &                 ABS(REAL(ZIN(1)-ZIN(NPLA)))
189                 ELSE                 ELSE
190                    YY(I) = CY                    YY(I) = CY
191                    BAR(M,I) = YY(I)                                      BAR(M,I) = YY(I)                  
192                    TBAR(M,I) = (-YALIG + YY(I))/10.                        TBAR(M,I) = (-YALIG + YY(I))/10.    
193                    IF (I.EQ.NPLA) MY=ABS(Y(1)-Y(NPLA))/                    IF (I.EQ.NPLA) MY=ABS(Y(1)-Y(NPLA))/
194       &                 ABS(ZIN(1)-ZIN(NPLA))       &                 ABS(REAL(ZIN(1)-ZIN(NPLA)))
195                 ENDIF                 ENDIF
196                 CALL LASTRISCIA(BAR(M,I),IBAR(M,I))                 CALL LASTRISCIA(BAR(M,I),IBAR(M,I))
197                 tibar(M,I) = ibar(m,i)                 tibar(M,I) = ibar(m,i)
# Line 306  c                     CBAR(M,I) = Y(I) Line 307  c                     CBAR(M,I) = Y(I)
307                       BAR(M,I) =  Y(I) + XALIG                       BAR(M,I) =  Y(I) + XALIG
308                       CBAR(M,I) = Y(I) / 10.                       CBAR(M,I) = Y(I) / 10.
309                       IF (I.EQ.NPLA) MX=ABS(Y(1)-Y(NPLA))/                       IF (I.EQ.NPLA) MX=ABS(Y(1)-Y(NPLA))/
310       &                    ABS(ZIN(1)-ZIN(NPLA))       &                    ABS(REAL(ZIN(1)-ZIN(NPLA)))
311  C      C    
312                    ELSE                    ELSE
313                       DISTY = PIANO(I)                                       DISTY = PIANO(I)                
# Line 316  c                     CBAR(M,I) = YY(I) Line 317  c                     CBAR(M,I) = YY(I)
317                       BAR(M,I) =  YY(I) + YALIG                       BAR(M,I) =  YY(I) + YALIG
318                       CBAR(M,I) = YY(I) / 10.                       CBAR(M,I) = YY(I) / 10.
319                       IF (I.EQ.NPLA) MY=ABS(Y(1)-Y(NPLA))/                       IF (I.EQ.NPLA) MY=ABS(Y(1)-Y(NPLA))/
320       &                    ABS(ZIN(1)-ZIN(NPLA))       &                    ABS(REAL(ZIN(1)-ZIN(NPLA)))
321  C      C    
322                    ENDIF                    ENDIF
323                    CALL LASTRISCIA(BAR(M,I),IBAR(M,I))                    CALL LASTRISCIA(BAR(M,I),IBAR(M,I))
# Line 376  C Line 377  C
377           GTR = 1           GTR = 1
378           IF (TRKCHI2.LT.0) GTR = 2           IF (TRKCHI2.LT.0) GTR = 2
379           IF ( AL_PP(GTR,5).NE.0. ) THEN           IF ( AL_PP(GTR,5).NE.0. ) THEN
380              RIG = 1./(AL_PP(GTR,5))              RIG = REAL(1./(AL_PP(GTR,5)))
381           ELSE           ELSE
382              GOOD2 = 0              GOOD2 = 0
383  c            PRINT *,' CALORIMETER - WARNING F77: track with R = 0'  c            PRINT *,' CALORIMETER - WARNING F77: track with R = 0'
Legend:
Removed from v.1.23  
changed lines
  Added in v.1.24
  ViewVC Help
Powered by ViewVC 1.1.23