/[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.2 by mocchiut, Mon Jan 22 09:17:01 2007 UTC revision 1.5 by mocchiut, Sat Jan 27 06:22:20 2007 UTC
# Line 19  C--------------------------------------- Line 19  C---------------------------------------
19    
20        INCLUDE 'INTEST.TXT'        INCLUDE 'INTEST.TXT'
21    
22        REAL FLIMIT        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    c      PARAMETER (FLIMIT2=2000)
27    c      PARAMETER (FLIMIT=500)
28          PARAMETER (FLIMIT=2000)
29          PARAMETER (FLIMIT2=500)
30        REAL X(NPLA),Y(NPLA),W(NPLA)        REAL X(NPLA),Y(NPLA),W(NPLA)
31        REAL DEVIA(2),TG(2)        REAL DEVIA(2),TG(2)
32        REAL BAR(2,NPLA)        REAL BAR(2,NPLA)
# Line 95  C                     ICONT = ICONT -1 Line 99  C                     ICONT = ICONT -1
99                 ENDIF                 ENDIF
100              ENDIF              ENDIF
101  C  C
102  C More importance to weigth and less cluster per plane  C No weigth and less cluster per plane
103  C  C
104              IF(FMODE(M).EQ.1) THEN              IF(FMODE(M).EQ.1) THEN
105                 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
106                    ICONT = ICONT + 1                    ICONT = ICONT + 1
107                    Y(I) = CLUS(M,I,1)                    Y(I) = CLUS(M,I,1)
108                    W(I) = ((CLUS(M,I,1+NCHA/2)**0.29))**2.                    W(I) = ((CLUS(M,I,1+NCHA/2)**0.29))**2.
109    c                  W(I) = 1.
110                 ENDIF                 ENDIF
111                 IF (W(I).GT.2000.) THEN  c               IF (W(I).GT.2000.) THEN
112                    W(I) = 2000.  cx                  W(I) = 2000.
113  c                  Y(I) = 0.  cc                  Y(I) = 0.
114  c                  ICONT = ICONT -1  cc                  ICONT = ICONT -1
115                 ENDIF  c               ENDIF
116              ENDIF              ENDIF
117  c  c
118  c One cluster per plane  c One cluster per plane
# Line 117  c Line 122  c
122                    ICONT = ICONT + 1                    ICONT = ICONT + 1
123                    Y(I) = CLUS(M,I,1)                    Y(I) = CLUS(M,I,1)
124                    W(I) = ((CLUS(M,I,1+NCHA/2)**0.79))**2.                    W(I) = ((CLUS(M,I,1+NCHA/2)**0.79))**2.
125    c                  W(I) = 1.
126                 ENDIF                 ENDIF
127              ENDIF              ENDIF
128  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 146  C Line 152  C
152              GOTO 10              GOTO 10
153           ENDIF           ENDIF
154  C  C
155  C TRY AGAIN IF STILL BAD AND WE HAVE AT LEAST 5 POINT TO MAKE THE FIT  C TRY AGAIN IF STILL BAD AND WE HAVE AT LEAST 10 POINT TO MAKE THE FIT
156  C  C
157           IF ((ABS(VAR).GT.FLIMIT.OR.VAR.EQ.0.).AND.FMODE(M).EQ.1           IF ((ABS(VAR).GT.FLIMIT2.OR.VAR.EQ.0.).AND.FMODE(M).EQ.1
158       &        .AND.NPFIT(M).GT.4) THEN       &        .AND.SNPFIT.GT.5) THEN
159              SSVAR = VAR              SSVAR = VAR
160              SSB = B              SSB = B
161              SSA = A              SSA = A
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.23