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

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

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.3 by mocchiut, Mon Jan 22 10:46:00 2007 UTC revision 1.4 by mocchiut, Tue Jan 23 11:52:26 2007 UTC
# Line 22  C--------------------------------------- Line 22  C---------------------------------------
22        REAL FLIMIT, FLIMIT2        REAL FLIMIT, FLIMIT2
23        INTEGER NVAR,NPAR        INTEGER NVAR,NPAR
24        PARAMETER (NVAR=1,NPAR=2)                  PARAMETER (NVAR=1,NPAR=2)          
25        PARAMETER (FLIMIT=10000)  c      PARAMETER (FLIMIT=10000)
26        PARAMETER (FLIMIT2=2000)  c      PARAMETER (FLIMIT2=2000)
27          PARAMETER (FLIMIT=500)
28          PARAMETER (FLIMIT2=500)
29        REAL X(NPLA),Y(NPLA),W(NPLA)        REAL X(NPLA),Y(NPLA),W(NPLA)
30        REAL DEVIA(2),TG(2)        REAL DEVIA(2),TG(2)
31        REAL BAR(2,NPLA)        REAL BAR(2,NPLA)
# Line 96  C                     ICONT = ICONT -1 Line 98  C                     ICONT = ICONT -1
98                 ENDIF                 ENDIF
99              ENDIF              ENDIF
100  C  C
101  C More importance to weigth and less cluster per plane  C No weigth and less cluster per plane
102  C  C
103              IF(FMODE(M).EQ.1) THEN              IF(FMODE(M).EQ.1) THEN
104                 IF (NCL(M,I).GE.1.AND.NCL(M,I).LE.2) THEN                 IF (NCL(M,I).GE.1.AND.NCL(M,I).LE.2) THEN
105                    ICONT = ICONT + 1                    ICONT = ICONT + 1
106                    Y(I) = CLUS(M,I,1)                    Y(I) = CLUS(M,I,1)
107                    W(I) = ((CLUS(M,I,1+NCHA/2)**0.29))**2.  c                  W(I) = ((CLUS(M,I,1+NCHA/2)**0.29))**2.
108                 ENDIF                    W(I) = 1.
                IF (W(I).GT.2000.) THEN  
                   W(I) = 2000.  
 c                  Y(I) = 0.  
 c                  ICONT = ICONT -1  
109                 ENDIF                 ENDIF
110    c               IF (W(I).GT.2000.) THEN
111    cx                  W(I) = 2000.
112    cc                  Y(I) = 0.
113    cc                  ICONT = ICONT -1
114    c               ENDIF
115              ENDIF              ENDIF
116  c  c
117  c One cluster per plane  c One cluster per plane
# Line 117  c Line 120  c
120                 IF (NCL(M,I).EQ.1) THEN                 IF (NCL(M,I).EQ.1) THEN
121                    ICONT = ICONT + 1                    ICONT = ICONT + 1
122                    Y(I) = CLUS(M,I,1)                    Y(I) = CLUS(M,I,1)
123                    W(I) = ((CLUS(M,I,1+NCHA/2)**0.79))**2.  c                  W(I) = ((CLUS(M,I,1+NCHA/2)**0.79))**2.
124                      W(I) = 1.
125                 ENDIF                 ENDIF
126              ENDIF              ENDIF
127  c            if (Y(I).eq.0.or.W(I).eq.2000.) NPFIT(M) = NPFIT(M) - 1  c            if (Y(I).eq.0.or.W(I).eq.2000.) NPFIT(M) = NPFIT(M) - 1
# Line 147  C Line 151  C
151              GOTO 10              GOTO 10
152           ENDIF           ENDIF
153  C  C
154  C TRY AGAIN IF STILL BAD AND WE HAVE AT LEAST 15 POINT TO MAKE THE FIT  C TRY AGAIN IF STILL BAD AND WE HAVE AT LEAST 10 POINT TO MAKE THE FIT
155  C  C
156           IF ((ABS(VAR).GT.FLIMIT2.OR.VAR.EQ.0.).AND.FMODE(M).EQ.1           IF ((ABS(VAR).GT.FLIMIT2.OR.VAR.EQ.0.).AND.FMODE(M).EQ.1
157       &        .AND.NPFIT(M).GT.14) THEN       &        .AND.SNPFIT.GT.5) THEN
158              SSVAR = VAR              SSVAR = VAR
159              SSB = B              SSB = B
160              SSA = A              SSA = A
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.23