/[PAMELA software]/DarthVader/CalorimeterLevel2/src/direction.for

Diff of /DarthVader/CalorimeterLevel2/src/direction.for

Parent Directory | Revision Log | View Patch Patch

-revision 1.6 by mocchiut,
Fri Jul 20 08:24:53 2007 UTC
+revision 1.9 by mocchiut,
Mon Jul 27 13:58:40 2009 UTC
 Line 38 
 c      PARAMETER (FLIMIT=500)
        INTEGER ICONT, M, I, SNPFIT,SSNPFIT
        INTEGER NPFIT(2), CHTRACK,IWPL(2)
+       INTEGER FORCEMODE
        INTEGER IBAR(2,NPLAV)
        COMMON/ANGOLO/BAR,IBAR
-Line 61 
 C
+Line 62 
 C
        VARFIT(1) = 0.
        VARFIT(2) = 0.
  C
+       FORCEMODE = -1
+       IF (FMODE(1).GE.1000.AND.FMODE(2).GE.1000) THEN
+          FORCEMODE = FMODE(1)-1000
+       ENDIF
+ C
  C THE FIT IS PERFORMED WHEN THERE ARE AT LEAST TWO
  C LAYERS WITH SIGNAL . FOR THE FIT THE CLUSTER WITH THE
  C HIGHEST DETECTED ENERGY IS USED .
-Line 78 
 C
+Line 84 
 C
           DO I = 1,NPLA
  C
              IF (MOD(M,2).EQ.0) THEN
-                X(I) = PIANO(I) - 5.81
+                X(I) = PIANO(I)
+ C               X(I) = PIANO(I) - 5.81
              ELSE
-                X(I) = PIANO(I)
+                X(I) = PIANO(I) - 5.81
+ C               X(I) = PIANO(I)
              ENDIF
              Y(I) = 0.
  C
  C
  C
-             IF (FMODE(M).EQ.0) THEN
+             IF ((FMODE(M).EQ.0.AND.FORCEMODE.EQ.-1)
+      +           .OR.FORCEMODE.EQ.0) THEN
                 IF (NCL(M,I).GE.1) THEN
                    ICONT = ICONT + 1
                    Y(I) = CLUS(M,I,1)
-Line 101 
 C                     ICONT = ICONT -1
+Line 110 
 C                     ICONT = ICONT -1
  C
  C No weigth and less cluster per plane
  C
-             IF(FMODE(M).EQ.1) THEN
+             IF((FMODE(M).EQ.1.AND.FORCEMODE.EQ.-1).OR.
+      +           FORCEMODE.EQ.1) THEN
                 IF (NCL(M,I).GE.1.AND.NCL(M,I).LE.2) THEN
                    ICONT = ICONT + 1
                    Y(I) = CLUS(M,I,1)
-Line 117 
 c               ENDIF
+Line 127 
 c               ENDIF
  c
  c One cluster per plane
  c
-             IF(FMODE(M).EQ.2) THEN
+             IF((FMODE(M).EQ.2.AND.FORCEMODE.EQ.-1).OR.
+      +           FORCEMODE.EQ.2) THEN
                 IF (NCL(M,I).EQ.1) THEN
                    ICONT = ICONT + 1
                    Y(I) = CLUS(M,I,1)
-Line 140 
 c         print *,' icont = ',icont
+Line 151 
 c         print *,' icont = ',icont
              NPFIT(M) = 0
           ENDIF
  C
- C TRY A BETTER FIT IF VAR IS BIG ON THE FIRST TRY
+          IF (FORCEMODE.EQ.-1) THEN
  C
-          IF ((ABS(VAR).GT.FLIMIT.OR.VAR.EQ.0.).AND.FMODE(M).EQ.0) THEN
+ C     TRY A BETTER FIT IF VAR IS BIG ON THE FIRST TRY
-             SVAR = VAR
+ C
-             SB = B
+             IF ((ABS(VAR).GT.FLIMIT.OR.VAR.EQ.0.)
-             SA = A
+      &           .AND.FMODE(M).EQ.0) THEN
-             SNPFIT = NPFIT(M)
+                SVAR = VAR
-             NPFIT(M) = NPLA
+                SB = B
-             FMODE(M) = 1
+                SA = A
-             GOTO 10
+                SNPFIT = NPFIT(M)
-          ENDIF
+                NPFIT(M) = NPLA
- C
- C TRY AGAIN IF STILL BAD AND WE HAVE AT LEAST 10 POINT TO MAKE THE FIT
- C
-          IF ((ABS(VAR).GT.FLIMIT2.OR.VAR.EQ.0.).AND.FMODE(M).EQ.1
-      &        .AND.SNPFIT.GT.5) THEN
-             SSVAR = VAR
-             SSB = B
-             SSA = A
-             SSNPFIT = NPFIT(M)
-             NPFIT(M) = NPLA
-             FMODE(M) = 2
-             GOTO 10
-          ENDIF
- C
- C IF WE ARE AT THE THIRD TRY AND THE FIT IS WORST THAN THE SECOND TRY TAKE THE FIRST TRY
- C
-          IF (FMODE(M).EQ.2) THEN
-             IF (ABS(VAR).GT.ABS(SSVAR).OR.ABS(VAR).EQ.0.
-      &           .OR.NPFIT(M).LE.4) THEN
                 FMODE(M) = 1
-                VAR=SSVAR
+                GOTO 10
-                B = SSB
-                A = SSA
-                NPFIT(M) = SSNPFIT
              ENDIF
-          ENDIF
+ C
- C
+ C     TRY AGAIN IF STILL BAD AND WE HAVE AT LEAST 10 POINT TO MAKE THE FIT
- C IF THE SECOND TRY AND THE FIT IS WORST THAN THE FIRST TRY TAKE THE FIRST TRY
+ C
- C
+             IF ((ABS(VAR).GT.FLIMIT2.OR.VAR.EQ.0.).AND.FMODE(M).EQ.1
-          IF (FMODE(M).EQ.1) THEN
+      &           .AND.SNPFIT.GT.5) THEN
-             IF (ABS(VAR).GT.ABS(SVAR).OR.ABS(VAR).EQ.0.) THEN
+                SSVAR = VAR
-                FMODE(M) = 0
+                SSB = B
-                VAR=SVAR
+                SSA = A
-                B = SB
+                SSNPFIT = NPFIT(M)
-                A = SA
+                NPFIT(M) = NPLA
-                NPFIT(M) = SNPFIT
+                FMODE(M) = 2
+                GOTO 10
+             ENDIF
+ C
+ C     IF WE ARE AT THE THIRD TRY AND THE FIT IS WORST THAN THE SECOND TRY TAKE THE FIRST TRY
+ C
+             IF (FMODE(M).EQ.2) THEN
+                IF (ABS(VAR).GT.ABS(SSVAR).OR.ABS(VAR).EQ.0.
+      &              .OR.NPFIT(M).LE.4) THEN
+                   FMODE(M) = 1
+                   VAR=SSVAR
+                   B = SSB
+                   A = SSA
+                   NPFIT(M) = SSNPFIT
+                ENDIF
+             ENDIF
+ C
+ C     IF THE SECOND TRY AND THE FIT IS WORST THAN THE FIRST TRY TAKE THE FIRST TRY
+ C
+             IF (FMODE(M).EQ.1) THEN
+                IF (ABS(VAR).GT.ABS(SVAR).OR.ABS(VAR).EQ.0.) THEN
+                   FMODE(M) = 0
+                   VAR=SVAR
+                   B = SB
+                   A = SA
+                   NPFIT(M) = SNPFIT
+                ENDIF
              ENDIF
           ENDIF
  C

 Legend:



Removed from v.1.6
 


changed lines


 
Added in v.1.9
 Legend:



Removed from v.1.6
 


changed lines


 
Added in v.1.9
-Removed from v.1.6
+Added in v.1.9

	ViewVC Help
Powered by ViewVC 1.1.23