/[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.4 by mocchiut, Tue Jan 23 11:52:26 2007 UTC revision 1.5 by mocchiut, Sat Jan 27 06:22:20 2007 UTC
# Line 24  C--------------------------------------- Line 24  C---------------------------------------
24        PARAMETER (NVAR=1,NPAR=2)                  PARAMETER (NVAR=1,NPAR=2)          
25  c      PARAMETER (FLIMIT=10000)  c      PARAMETER (FLIMIT=10000)
26  c      PARAMETER (FLIMIT2=2000)  c      PARAMETER (FLIMIT2=2000)
27        PARAMETER (FLIMIT=500)  c      PARAMETER (FLIMIT=500)
28          PARAMETER (FLIMIT=2000)
29        PARAMETER (FLIMIT2=500)        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)
# Line 104  C Line 105  C
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  c                  W(I) = ((CLUS(M,I,1+NCHA/2)**0.29))**2.                    W(I) = ((CLUS(M,I,1+NCHA/2)**0.29))**2.
109                    W(I) = 1.  c                  W(I) = 1.
110                 ENDIF                 ENDIF
111  c               IF (W(I).GT.2000.) THEN  c               IF (W(I).GT.2000.) THEN
112  cx                  W(I) = 2000.  cx                  W(I) = 2000.
# Line 120  c Line 121  c
121                 IF (NCL(M,I).EQ.1) THEN                 IF (NCL(M,I).EQ.1) THEN
122                    ICONT = ICONT + 1                    ICONT = ICONT + 1
123                    Y(I) = CLUS(M,I,1)                    Y(I) = CLUS(M,I,1)
124  c                  W(I) = ((CLUS(M,I,1+NCHA/2)**0.79))**2.                    W(I) = ((CLUS(M,I,1+NCHA/2)**0.79))**2.
125                    W(I) = 1.  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
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.23