/[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.13 by mocchiut, Sat Jan 27 06:22:19 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 252  C Line 253  C
253  C      C    
254           DO M = 1,2           DO M = 1,2
255  C  C
256    c            print *,' ax ',ax,' ay ',ay
257    c            print *,' bx ',bx,' by ',by
258              IF (NPCFIT(M).GE.2) THEN              IF (NPCFIT(M).GE.2) THEN
259                 IF (M.EQ.1) THEN                 IF (M.EQ.1) THEN
260                    VARCFIT(1) = CHI2X                    VARCFIT(1) = CHI2X
261                    IMPX = 10. * ( AX + 12.1 )                    IMPX = AX ! PAMELA REF
262                    TANX = BX                    TANX = BX
263                 ELSE                 ELSE
264                    VARCFIT(2) = CHI2Y                    VARCFIT(2) = CHI2Y
265                    IMPY = 10. * ( AY + 12.2 )                    IMPY = AY ! PAMELA REF
266                    TANY = BY                    TANY = BY
267                 ENDIF                 ENDIF
268  C  C
# Line 274  C?????                  IF (M.EQ.2) NN = Line 277  C?????                  IF (M.EQ.2) NN =
277  C      C    
278                    IF (M.EQ.1) THEN                    IF (M.EQ.1) THEN
279                       DISTX = PIANO(I) - 5.81                       DISTX = PIANO(I) - 5.81
280                       Y(I) = DISTX * TANX + CX                       Y(I) = (DISTX * TANX) +  AX - XALIG
281    c                     CBAR(M,I) = Y(I)
282                       BAR(M,I) = Y(I)                       BAR(M,I) = Y(I)
283                       CBAR(M,I) = Y(I)                               CBAR(M,I) = (Y(I) + XALIG)/10.
284                       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))
285  C      C    
286                    ELSE                    ELSE
287                       DISTY = PIANO(I)                                       DISTY = PIANO(I)                
288                       YY(I) = DISTY * TANY + CY                       YY(I) = (DISTY * TANY) + AY - YALIG
289    c                     CBAR(M,I) = YY(I)
290                       BAR(M,I) = YY(I)                       BAR(M,I) = YY(I)
291                       CBAR(M,I) = YY(I)                       CBAR(M,I) = (YY(I) + YALIG)/10.
292                       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))
293  C      C    
294                    ENDIF                    ENDIF
295                    CALL LASTRISCIA(BAR(M,I),IBAR(M,I))                    CALL LASTRISCIA(BAR(M,I),IBAR(M,I))
296                    cibar(M,I) = ibar(m,i)                    cibar(M,I) = ibar(m,i)
297                      IF (ibar(m,i).EQ.-1) THEN
298                         CHTRACK = CHTRACK + 1
299                      ELSE
300                         IWPL(M) = IWPL(M) + 1
301                      ENDIF
302                 ENDDO                             ENDDO            
303              ENDIF              ENDIF
304  C  C
305           ENDDO           ENDDO
306  C  C
307        ENDIF        ENDIF
308  C      IF (TRIGTY.GE.2.AND.HZN.NE.0) THEN  C
309  C         IF (GOOD2.EQ.0) THEN        IF (TRIGTY.GE.2.AND.HZN.NE.0) THEN
310  C            PRINT *,' CALORIMETER - WARNING F77: unknown request'           IF (GOOD2.EQ.1) THEN
311  C            GOOD2 = 1              PRINT *,' CALORIMETER - WARNING F77: unknown request'
312  C            GOTO 50              GOOD2 = 1
313  C         ENDIF              GOTO 50
314  C      ENDIF           ENDIF
315             IF ( NPCFIT(1).EQ.0.OR.NPCFIT(2).EQ.0 ) THEN
316                GOOD2 = 1
317                GOTO 50
318             ENDIF
319          ENDIF
320  C  C
321   6996 CONTINUE   6996 CONTINUE
322  C  C
# Line 344  C Line 359  C
359        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
360           RIG = ELEN ! SELFTRIGGER RIGIDITY           RIG = ELEN ! SELFTRIGGER RIGIDITY
361           IF ( RIG.EQ.0. ) THEN           IF ( RIG.EQ.0. ) THEN
362              GOOD2 = 0              GOOD2 = 1
363              PRINT *,' CALORIMETER - WARNING F77: ST track with R = 0'              PRINT *,' CALORIMETER - WARNING F77: ST track with R = 0'
364              GOTO 50              GOTO 50
365           ENDIF           ENDIF
# Line 601  c            ISUPY = NNY + 8 Line 616  c            ISUPY = NNY + 8
616           ENDIF           ENDIF
617        ENDDO        ENDDO
618  C  C
       EINF = EMIN  
       ESUP = 50.  
619  C  C
620  C     CALCULATE PLANETOT AND QMEAN  C     CALCULATE PLANETOT AND QMEAN
621  C  C
# Line 613  C Line 626  C
626        NPIANI = 5        NPIANI = 5
627        QMEAN = 0.        QMEAN = 0.
628        INDEX = 0        INDEX = 0
629        CALL ELIO(RPIANO,NPIANI,QMEAN,NTOT,INDEX)  C
630        PLANETOT = RPIANO(1) + RPIANO(2)          IF (TRIGTY.GE.2.AND.HZN.NE.0) THEN
631             EINF = 50.
632             ESUP = 15000.
633             CALL NUCLEI(RPIANO,NPIANI,QMEAN,NTOT,INDEX)
634             PLANETOT = RPIANO(1) + RPIANO(2)  
635          ELSE
636             EINF = EMIN
637             ESUP = 15000.
638             CALL ELIO(RPIANO,NPIANI,QMEAN,NTOT,INDEX)
639             PLANETOT = RPIANO(1) + RPIANO(2)  
640          ENDIF
641  C  C
642   50   CONTINUE   50   CONTINUE
643  C  C
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.23