/[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.8 by mocchiut, Mon Jan 22 09:17:01 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 289  C     Line 290  C    
290                    ENDIF                    ENDIF
291                    CALL LASTRISCIA(BAR(M,I),IBAR(M,I))                    CALL LASTRISCIA(BAR(M,I),IBAR(M,I))
292                    cibar(M,I) = ibar(m,i)                    cibar(M,I) = ibar(m,i)
293                      IF (ibar(m,i).EQ.-1) THEN
294                         CHTRACK = CHTRACK + 1
295                      ELSE
296                         IWPL(M) = IWPL(M) + 1
297                      ENDIF
298                 ENDDO                             ENDDO            
299              ENDIF              ENDIF
300  C  C
301           ENDDO           ENDDO
302  C  C
303        ENDIF        ENDIF
304  C      IF (TRIGTY.GE.2.AND.HZN.NE.0) THEN  C
305  C         IF (GOOD2.EQ.0) THEN        IF (TRIGTY.GE.2.AND.HZN.NE.0) THEN
306  C            PRINT *,' CALORIMETER - WARNING F77: unknown request'           IF (GOOD2.EQ.1) THEN
307  C            GOOD2 = 1              PRINT *,' CALORIMETER - WARNING F77: unknown request'
308  C            GOTO 50              GOOD2 = 1
309  C         ENDIF              GOTO 50
310  C      ENDIF           ENDIF
311             IF ( NPCFIT(1).EQ.0.OR.NPCFIT(2).EQ.0 ) THEN
312                GOOD2 = 1
313                GOTO 50
314             ENDIF
315          ENDIF
316  C  C
317   6996 CONTINUE   6996 CONTINUE
318  C  C
# Line 344  C Line 355  C
355        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
356           RIG = ELEN ! SELFTRIGGER RIGIDITY           RIG = ELEN ! SELFTRIGGER RIGIDITY
357           IF ( RIG.EQ.0. ) THEN           IF ( RIG.EQ.0. ) THEN
358              GOOD2 = 0              GOOD2 = 1
359              PRINT *,' CALORIMETER - WARNING F77: ST track with R = 0'              PRINT *,' CALORIMETER - WARNING F77: ST track with R = 0'
360              GOTO 50              GOTO 50
361           ENDIF           ENDIF
# Line 602  c            ISUPY = NNY + 8 Line 613  c            ISUPY = NNY + 8
613        ENDDO        ENDDO
614  C  C
615        EINF = EMIN        EINF = EMIN
616        ESUP = 50.        ESUP = 150.
617    C
618          IF (TRIGTY.GE.2.AND.HZN.NE.0) THEN
619             EINF = 50.
620             ESUP = 1500.
621          ENDIF
622  C  C
623  C     CALCULATE PLANETOT AND QMEAN  C     CALCULATE PLANETOT AND QMEAN
624  C  C
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.23