/[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.19 by mocchiut, Thu Nov 29 16:11:53 2007 UTC revision 1.23 by mocchiut, Wed May 23 12:28:33 2012 UTC
# Line 21  C Line 21  C
21        REAL ax,bx,eax,ebx,chi2x        REAL ax,bx,eax,ebx,chi2x
22        REAL ay,by,eay,eby,chi2y        REAL ay,by,eay,eby,chi2y
23        REAL parzen3, TMISD        REAL parzen3, TMISD
24        INTEGER Nfitx,Nfity        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        INTEGER j, m, i, IWPL(2), timpx, timpy, T, nn
# Line 96  C     BEGIN WITH THE FIRST TRACK IF WE H Line 96  C     BEGIN WITH THE FIRST TRACK IF WE H
96  C  C
97        T = 1        T = 1
98  C  C
99   10   CONTINUE  c 10   CONTINUE
100          CONTINUE
101  C  C
102        IF (GOOD2.EQ.1) THEN        IF (GOOD2.EQ.1) THEN
103  C      C    
# Line 116  C Line 117  C
117  c            print *,' al_p(',m,') = ',al_p(m)  c            print *,' al_p(',m,') = ',al_p(m)
118           enddo           enddo
119           if (al_p(5).eq.0.) THEN           if (al_p(5).eq.0.) THEN
120         PRINT *,' CALORIMETER - WARNING F77: track with R = 0, discarded'  c       PRINT *,' CALORIMETER - WARNING F77: track with R = 0, discarded'
121              GOOD2 = 0              GOOD2 = 0
122              GOTO 969              GOTO 969
123           ENDIF           ENDIF
# Line 273  c            print *,' bx ',bx,' by ',by Line 274  c            print *,' bx ',bx,' by ',by
274              IF (NPCFIT(M).GE.2) THEN              IF (NPCFIT(M).GE.2) THEN
275                 IF (M.EQ.1) THEN                 IF (M.EQ.1) THEN
276                    VARCFIT(1) = CHI2X                    VARCFIT(1) = CHI2X
277                    IMPX = AX ! PAMELA REF                    IMPX = AX + BX * (ZALIG/10.) ! PAMELA REF
278                    TANX = BX                    TANX = BX
279                 ELSE                 ELSE
280                    VARCFIT(2) = CHI2Y                    VARCFIT(2) = CHI2Y
281                    IMPY = AY ! PAMELA REF                    IMPY = AY + BY * (ZALIG/10.) ! PAMELA REF
282                    TANY = BY                    TANY = BY
283                 ENDIF                 ENDIF
284  C  C
# Line 300  c                  IF (M.EQ.2) NN = 1 Line 301  c                  IF (M.EQ.2) NN = 1
301  C      C    
302                    IF (M.EQ.1) THEN                    IF (M.EQ.1) THEN
303                       DISTX = PIANO(I) - 5.81                       DISTX = PIANO(I) - 5.81
304                       Y(I) = (DISTX * TANX) +  AX - XALIG                       Y(I) = (DISTX * TANX) +  IMPX*10. !- XALIG
305  c                     CBAR(M,I) = Y(I)  c                     CBAR(M,I) = Y(I)
306                       BAR(M,I) = Y(I)                       BAR(M,I) =  Y(I) + XALIG
307                       CBAR(M,I) = (Y(I) + XALIG)/10.                       CBAR(M,I) = Y(I) / 10.
308                       IF (I.EQ.NPLA) MX=ABS(Y(1)-Y(NPLA))/                       IF (I.EQ.NPLA) MX=ABS(Y(1)-Y(NPLA))/
309       &                    ABS(ZIN(1)-ZIN(NPLA))       &                    ABS(ZIN(1)-ZIN(NPLA))
310  C      C    
311                    ELSE                    ELSE
312                       DISTY = PIANO(I)                                       DISTY = PIANO(I)                
313                       YY(I) = (DISTY * TANY) + AY - YALIG                       YY(I) = (DISTY * TANY) + IMPY*10. !- YALIG
314    c                     print *,' I ',i,' YY ',YY(I)
315  c                     CBAR(M,I) = YY(I)  c                     CBAR(M,I) = YY(I)
316                       BAR(M,I) = YY(I)                       BAR(M,I) =  YY(I) + YALIG
317                       CBAR(M,I) = (YY(I) + YALIG)/10.                       CBAR(M,I) = YY(I) / 10.
318                       IF (I.EQ.NPLA) MY=ABS(Y(1)-Y(NPLA))/                       IF (I.EQ.NPLA) MY=ABS(Y(1)-Y(NPLA))/
319       &                    ABS(ZIN(1)-ZIN(NPLA))       &                    ABS(ZIN(1)-ZIN(NPLA))
320  C      C    
# Line 377  C Line 379  C
379              RIG = 1./(AL_PP(GTR,5))              RIG = 1./(AL_PP(GTR,5))
380           ELSE           ELSE
381              GOOD2 = 0              GOOD2 = 0
382              PRINT *,' CALORIMETER - WARNING F77: track with R = 0'  c            PRINT *,' CALORIMETER - WARNING F77: track with R = 0'
383              GOTO 50              GOTO 50
384           ENDIF           ENDIF
385        ENDIF        ENDIF
# Line 385  C Line 387  C
387           RIG = ELEN ! SELFTRIGGER RIGIDITY           RIG = ELEN ! SELFTRIGGER RIGIDITY
388           IF ( RIG.EQ.0. ) THEN           IF ( RIG.EQ.0. ) THEN
389              GOOD2 = 1              GOOD2 = 1
390              PRINT *,' CALORIMETER - WARNING F77: ST track with R = 0'  c            PRINT *,' CALORIMETER - WARNING F77: ST track with R = 0'
391              GOTO 50              GOTO 50
392           ENDIF           ENDIF
393        ENDIF        ENDIF
# Line 479  C     Line 481  C    
481              INFX = NNX - 8              INFX = NNX - 8
482              ISUPX = NNX + 8              ISUPX = NNX + 8
483              DO I = INFX,ISUPX              DO I = INFX,ISUPX
484                 IF (DEXY(1,J,I).LT.EMIN) GO TO 710                 IF (DEXY(1,J,I).GE.EMIN) THEN
485                 NCYL = NCYL + 1                    NCYL = NCYL + 1
486                 QCYL = QCYL + DEXY(1,J,I)                    QCYL = QCYL + DEXY(1,J,I)
487   710        ENDDO                 ENDIF
488                ENDDO
489           ENDIF           ENDIF
490           NNY = IBAR(2,J)           NNY = IBAR(2,J)
491           IF (NNY.NE.-1) THEN           IF (NNY.NE.-1) THEN
# Line 491  C     Line 494  C    
494              INFY = NNY - 8              INFY = NNY - 8
495              ISUPY = NNY + 8              ISUPY = NNY + 8
496              DO I=INFY,ISUPY              DO I=INFY,ISUPY
497                 IF (DEXY(2,J,I).LT.EMIN) GO TO 810                 IF (DEXY(2,J,I).GE.EMIN) THEN
498                 NCYL = NCYL + 1                    NCYL = NCYL + 1
499                 QCYL = QCYL + DEXY(2,J,I)                    QCYL = QCYL + DEXY(2,J,I)
500   810        ENDDO                 ENDIF
501                ENDDO
502           ENDIF           ENDIF
503  C      C    
504  C     QTR = DETECTED ENERGY AND NTR = NUMBER OF HIT STRIPS IN A CYLINDER oF  C     QTR = DETECTED ENERGY AND NTR = NUMBER OF HIT STRIPS IN A CYLINDER oF
# Line 606  C Line 610  C
610  C      C    
611  C     CALCULATE NLAST AND QLAST  C     CALCULATE NLAST AND QLAST
612  C  C
613        DO J = NPLA-4,NPLA        MNPLA = NPLA -4
614          IF ( MNPLA .LT. 1 ) MNPLA = 1
615          DO J = MNPLA,NPLA
616           NNX = IBAR(1,J)           NNX = IBAR(1,J)
617           IF (NNX.NE.-1) THEN           IF (NNX.NE.-1) THEN
618              IF (NNX.LT.5) NNX = 5              IF (NNX.LT.5) NNX = 5
Legend:
Removed from v.1.19  
changed lines
  Added in v.1.23
  ViewVC Help
Powered by ViewVC 1.1.23