/[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.16 by mocchiut, Fri Jul 20 08:24:53 2007 UTC revision 1.19 by mocchiut, Thu Nov 29 16:11:53 2007 UTC
# Line 92  c      print *,' sono qui' Line 92  c      print *,' sono qui'
92        QMEAN = 0.        QMEAN = 0.
93  C      SELFTRIGGER = 0  C      SELFTRIGGER = 0
94  C  C
95  C     BEGIN WITH THE FISRT TRACK IF WE HAVE A TRACK FROM TRACKER  C     BEGIN WITH THE FIRST TRACK IF WE HAVE A TRACK FROM TRACKER
96  C  C
97        T = 1        T = 1
98  C  C
# Line 126  c            print *,' al_p(',m,') = ',a Line 126  c            print *,' al_p(',m,') = ',a
126                 YOUT(I) = 0.                 YOUT(I) = 0.
127                 IF (MOD(M,2).EQ.0) THEN                 IF (MOD(M,2).EQ.0) THEN
128                    DISTX = PIANO(I) + ZALIG                    DISTX = PIANO(I) + ZALIG
129    c                  print *,'T Y PLANE I= ',I,' Z = ',DISTX
130                 ELSE                 ELSE
131                    DISTX = PIANO(I) - 5.81 + ZALIG                    DISTX = PIANO(I) - 5.81 + ZALIG
132    c                  print *,'T X PLANE I= ',I,' Z = ',DISTX
133                 ENDIF                               ENDIF              
134                 ZIN(I) = distx / 10.                 ZIN(I) = distx / 10.
135  c               print *,' ZIN(',I,') = ',ZIN(I)  c               print *,' ZIN(',I,') = ',ZIN(I)
# Line 145  c               print *,' CALORIMETER - Line 147  c               print *,' CALORIMETER -
147              TY = TAN(ASIN(AL_P(3))) * SIN(AL_P(4))              TY = TAN(ASIN(AL_P(3))) * SIN(AL_P(4))
148              DO I = 1, NPLA              DO I = 1, NPLA
149                 NN = 0                 NN = 0
150                 IF (M.EQ.2) NN = 1  c               IF (M.EQ.2) NN = 1
151                 IF (MOD(I,2).EQ.NN) THEN                 IF (MOD(I,2).EQ.NN) THEN
152                    IF (REVERSE.EQ.0) THEN                    IF (REVERSE.EQ.0) THEN
                      SHIFT = +0.5  
                   ELSE  
153                       SHIFT = -0.5                       SHIFT = -0.5
154                      ELSE
155                         SHIFT = +0.5
156                    ENDIF                    ENDIF
157                 ELSE                 ELSE
158                    IF (REVERSE.EQ.0) THEN                    IF (REVERSE.EQ.0) THEN
                      SHIFT = -0.5  
                   ELSE  
159                       SHIFT = +0.5                       SHIFT = +0.5
160                      ELSE
161                         SHIFT = -0.5
162                    ENDIF                    ENDIF
163                 ENDIF                 ENDIF
164  C      C    
# Line 281  c            print *,' bx ',bx,' by ',by Line 283  c            print *,' bx ',bx,' by ',by
283  C  C
284                 DO I = 1,NPLA                     DO I = 1,NPLA    
285                    NN = 0                    NN = 0
286                    IF (M.EQ.2) NN = 1  c                  IF (M.EQ.2) NN = 1
287                    IF (MOD(I,2).EQ.NN) THEN                    IF (MOD(I,2).EQ.NN) THEN
288                       IF (REVERSE.EQ.0) THEN                       IF (REVERSE.EQ.0) THEN
                         SHIFT = +0.5  
                      ELSE  
289                          SHIFT = -0.5                          SHIFT = -0.5
290                         ELSE
291                            SHIFT = +0.5
292                       ENDIF                       ENDIF
293                    ELSE                    ELSE
294                       IF (REVERSE.EQ.0) THEN                       IF (REVERSE.EQ.0) THEN
                         SHIFT = -0.5  
                      ELSE  
295                          SHIFT = +0.5                          SHIFT = +0.5
296                         ELSE
297                            SHIFT = -0.5
298                       ENDIF                       ENDIF
299                    ENDIF                    ENDIF
300  C      C    
Legend:
Removed from v.1.16  
changed lines
  Added in v.1.19
  ViewVC Help
Powered by ViewVC 1.1.23