/[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.7 by mocchiut, Thu Nov 22 13:53:28 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
81                 X(I) = PIANO(I) - 5.81                 X(I) = PIANO(I)
82    C               X(I) = PIANO(I) - 5.81
83              ELSE              ELSE
84                 X(I) = PIANO(I)                                   X(I) = PIANO(I) - 5.81
85    C               X(I) = PIANO(I)                  
86              ENDIF              ENDIF
87              Y(I) = 0.              Y(I) = 0.
88  C  C
# Line 104  C Line 107  C
107                 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
108                    ICONT = ICONT + 1                    ICONT = ICONT + 1
109                    Y(I) = CLUS(M,I,1)                    Y(I) = CLUS(M,I,1)
110  c                  W(I) = ((CLUS(M,I,1+NCHA/2)**0.29))**2.                    W(I) = ((CLUS(M,I,1+NCHA/2)**0.29))**2.
111                    W(I) = 1.  c                  W(I) = 1.
112                 ENDIF                 ENDIF
113  c               IF (W(I).GT.2000.) THEN  c               IF (W(I).GT.2000.) THEN
114  cx                  W(I) = 2000.  cx                  W(I) = 2000.
# Line 120  c Line 123  c
123                 IF (NCL(M,I).EQ.1) THEN                 IF (NCL(M,I).EQ.1) THEN
124                    ICONT = ICONT + 1                    ICONT = ICONT + 1
125                    Y(I) = CLUS(M,I,1)                    Y(I) = CLUS(M,I,1)
126  c                  W(I) = ((CLUS(M,I,1+NCHA/2)**0.79))**2.                    W(I) = ((CLUS(M,I,1+NCHA/2)**0.79))**2.
127                    W(I) = 1.  c                  W(I) = 1.
128                 ENDIF                 ENDIF
129              ENDIF              ENDIF
130  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.7
  ViewVC Help
Powered by ViewVC 1.1.23