/[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.7 by mocchiut, Tue Nov 14 14:08:50 2006 UTC revision 1.9 by mocchiut, Mon Jan 22 09:44:49 2007 UTC
# Line 55  C Line 55  C
55        COMMON / CH / CHECK        COMMON / CH / CHECK
56        SAVE / CH /        SAVE / CH /
57  C  C
58        COMMON / CALOFIT / VARFIT, NPFIT        COMMON / CALOFIT / VARFIT, NPFIT, IWPL,CHTRACK
59        SAVE / CALOFIT /        SAVE / CALOFIT /
60  C  C
61        COMMON / pawcd / hmemor        COMMON / pawcd / hmemor
# Line 89  C Line 89  C
89        NLAST = 0.        NLAST = 0.
90        PLANETOT = 0.        PLANETOT = 0.
91        QMEAN = 0.        QMEAN = 0.
92        SELFTRIGGER = 0  C      SELFTRIGGER = 0
93  C  C
94  C     BEGIN WITH THE FISRT TRACK IF WE HAVE A TRACK FROM TRACKER  C     BEGIN WITH THE FISRT TRACK IF WE HAVE A TRACK FROM TRACKER
95  C  C
# Line 232  C Line 232  C
232  C  C
233  C     CALL SELFTRIGGER SUBROUTINE  C     CALL SELFTRIGGER SUBROUTINE
234  C  C
235             CALL VZERO(IWPL,2)
236           CALL VZERO(VARCFIT,2)           CALL VZERO(VARCFIT,2)
237           CALL VZERO(NPCFIT,2)           CALL VZERO(NPCFIT,2)
238           CALL VZERO(TBAR,2*NPLA)           CALL VZERO(TBAR,2*NPLA)
# Line 276  C     Line 277  C    
277                       DISTX = PIANO(I) - 5.81                       DISTX = PIANO(I) - 5.81
278                       Y(I) = DISTX * TANX + CX                       Y(I) = DISTX * TANX + CX
279                       BAR(M,I) = Y(I)                       BAR(M,I) = Y(I)
280                       CBAR(M,I) = Y(I)          c                     CBAR(M,I) = Y(I)
281                         CBAR(M,I) = (Y(I) - XALIG)/10.        
282                       IF (I.EQ.22) MX=ABS(Y(1)-Y(22))/ABS(ZIN(1)-ZIN(22))                       IF (I.EQ.22) MX=ABS(Y(1)-Y(22))/ABS(ZIN(1)-ZIN(22))
283  C      C    
284                    ELSE                    ELSE
285                       DISTY = PIANO(I)                                       DISTY = PIANO(I)                
286                       YY(I) = DISTY * TANY + CY                       YY(I) = DISTY * TANY + CY
287                       BAR(M,I) = YY(I)                       BAR(M,I) = YY(I)
288                       CBAR(M,I) = YY(I)  c                     CBAR(M,I) = YY(I)
289                         CBAR(M,I) = (-YALIG + YY(I))/10.
290                       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))
291  C      C    
292                    ENDIF                    ENDIF
293                    CALL LASTRISCIA(BAR(M,I),IBAR(M,I))                    CALL LASTRISCIA(BAR(M,I),IBAR(M,I))
294                    cibar(M,I) = ibar(m,i)                    cibar(M,I) = ibar(m,i)
295                      IF (ibar(m,i).EQ.-1) THEN
296                         CHTRACK = CHTRACK + 1
297                      ELSE
298                         IWPL(M) = IWPL(M) + 1
299                      ENDIF
300                 ENDDO                             ENDDO            
301              ENDIF              ENDIF
302  C  C
303           ENDDO           ENDDO
304  C  C
305        ENDIF        ENDIF
306  C      IF (TRIGTY.GE.2.AND.HZN.NE.0) THEN  C
307  C         IF (GOOD2.EQ.0) THEN        IF (TRIGTY.GE.2.AND.HZN.NE.0) THEN
308  C            PRINT *,' CALORIMETER - WARNING F77: unknown request'           IF (GOOD2.EQ.1) THEN
309  C            GOOD2 = 1              PRINT *,' CALORIMETER - WARNING F77: unknown request'
310  C            GOTO 50              GOOD2 = 1
311  C         ENDIF              GOTO 50
312  C      ENDIF           ENDIF
313             IF ( NPCFIT(1).EQ.0.OR.NPCFIT(2).EQ.0 ) THEN
314                GOOD2 = 1
315                GOTO 50
316             ENDIF
317          ENDIF
318  C  C
319   6996 CONTINUE   6996 CONTINUE
320  C  C
# Line 344  C Line 357  C
357        IF (TRIGTY.GE.2.AND.HZN.EQ.0.AND.GOOD2.EQ.0) THEN        IF (TRIGTY.GE.2.AND.HZN.EQ.0.AND.GOOD2.EQ.0) THEN
358           RIG = ELEN ! SELFTRIGGER RIGIDITY           RIG = ELEN ! SELFTRIGGER RIGIDITY
359           IF ( RIG.EQ.0. ) THEN           IF ( RIG.EQ.0. ) THEN
360              GOOD2 = 0              GOOD2 = 1
361              PRINT *,' CALORIMETER - WARNING F77: ST track with R = 0'              PRINT *,' CALORIMETER - WARNING F77: ST track with R = 0'
362              GOTO 50              GOTO 50
363           ENDIF           ENDIF
# Line 602  c            ISUPY = NNY + 8 Line 615  c            ISUPY = NNY + 8
615        ENDDO        ENDDO
616  C  C
617        EINF = EMIN        EINF = EMIN
618        ESUP = 50.        ESUP = 150.
619    C
620          IF (TRIGTY.GE.2.AND.HZN.NE.0) THEN
621             EINF = 50.
622             ESUP = 1500.
623          ENDIF
624  C  C
625  C     CALCULATE PLANETOT AND QMEAN  C     CALCULATE PLANETOT AND QMEAN
626  C  C
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.23