/[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.6 by mocchiut, Fri Jul 20 08:24:53 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(NPLAV),Y(NPLAV),W(NPLAV)
31        REAL DEVIA(2),TG(2)        REAL DEVIA(2),TG(2)
32        REAL BAR(2,NPLA)        REAL BAR(2,NPLAV)
33        REAL A, B, VAR        REAL A, B, VAR
34        REAL SA, SB, SVAR        REAL SA, SB, SVAR
35        REAL SSA, SSB, SSVAR        REAL SSA, SSB, SSVAR
36        REAL PIANO(22)        REAL PIANO(NPLAV)
37        REAL VARFIT(2)        REAL VARFIT(2)
38    
39        INTEGER ICONT, M, I, SNPFIT,SSNPFIT        INTEGER ICONT, M, I, SNPFIT,SSNPFIT
40        INTEGER NPFIT(2), CHTRACK,IWPL(2)        INTEGER NPFIT(2), CHTRACK,IWPL(2)
41                
42        INTEGER IBAR(2,NPLA)              INTEGER IBAR(2,NPLAV)      
43        COMMON/ANGOLO/BAR,IBAR        COMMON/ANGOLO/BAR,IBAR
44        SAVE /ANGOLO/        SAVE /ANGOLO/
45    
# Line 71  C Line 72  C
72           FMODE(M) = 0           FMODE(M) = 0
73   10      CONTINUE           10      CONTINUE        
74           ICONT = 0           ICONT = 0
75           CALL VZERO(X,NPLA)           CALL VZERO(X,NPLAV)
76           CALL VZERO(Y,NPLA)           CALL VZERO(Y,NPLAV)
77           CALL VZERO(W,NPLA)           CALL VZERO(W,NPLAV)
78           DO I = 1,NPLA           DO I = 1,NPLA
79  C      C    
80              IF (MOD(M,2).EQ.0) THEN              IF (MOD(M,2).EQ.0) THEN
# 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.6
  ViewVC Help
Powered by ViewVC 1.1.23