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

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

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

revision 1.2 by mocchiut, Tue May 23 08:35:45 2006 UTC revision 1.7 by mocchiut, Fri Mar 30 11:17:16 2007 UTC
# Line 17  c Line 17  c
17        INTEGER gtr        INTEGER gtr
18        INTEGER j, m, nn, i        INTEGER j, m, nn, i
19        INTEGER IBAR(2,NPLA), NPFIT(2)        INTEGER IBAR(2,NPLA), NPFIT(2)
20        INTEGER Iquest(100)            INTEGER Iquest(100), IWPL(2),CHTRACK
21  c  c
22        COMMON / TAGLIOEN / EINF, ESUP, ENER(2)        COMMON / TAGLIOEN / EINF, ESUP, ENER(2)
23        SAVE / TAGLIOEN /        SAVE / TAGLIOEN /
# Line 37  C Line 37  C
37        COMMON / CH / CHECK        COMMON / CH / CHECK
38        SAVE / CH /        SAVE / CH /
39  C  C
40        COMMON / CALOFIT / VARFIT, NPFIT        COMMON / CALOFIT / VARFIT, NPFIT, IWPL,CHTRACK
41        SAVE / CALOFIT /        SAVE / CALOFIT /
42  C  C
43        COMMON / pawcd / hmemor        COMMON / pawcd / hmemor
# Line 86  C Line 86  C
86        QX22 = 0.        QX22 = 0.
87        QMAX = 0.        QMAX = 0.
88        GTR = 0        GTR = 0
89          CHTRACK = 0
90  C      C    
91        DISTX = 0.        DISTX = 0.
92        DISTY = 0.        DISTY = 0.
# Line 105  C Line 106  C
106           DO J = 1,96           DO J = 1,96
107              IF ( MOD(I,2).NE.0 ) THEN              IF ( MOD(I,2).NE.0 ) THEN
108                 IF ( ESTRIP(2,I,J).GT.EMIN ) THEN                 IF ( ESTRIP(2,I,J).GT.EMIN ) THEN
                   ENPLANEY = ENPLANEY + ESTRIP(2,I,J)  
109                    DEXY(2,I,J) = ESTRIP(2,I,J)                    DEXY(2,I,J) = ESTRIP(2,I,J)
110                      ENPLANEY = ENPLANEY + DEXY(2,I,J)
111                    NSTRIP = NSTRIP + 1.                    NSTRIP = NSTRIP + 1.
112                    QTOT = QTOT + ESTRIP(2,I,J)                    QTOT = QTOT + DEXY(2,I,J)
113                    IF (I.LT.11) QQ(1) = QQ(1) + ESTRIP(2,I,J)                    IF (I.LT.11) QQ(1) = QQ(1) + DEXY(2,I,J)
114                 ENDIF                 ENDIF
115                 IF ( ESTRIP(1,I,J).GT.EMIN ) THEN                 IF ( ESTRIP(1,I,J).GT.EMIN ) THEN
116                    ENPLANEX = ENPLANEX + ESTRIP(1,I,J)  C               IF ( ESTRIP(1,I+1,J).GT.EMIN ) THEN
117    C                  
118    c                  print *,' i ',i,' i + 1 ',I+1
119    C                  DEXY(1,I,J) = ESTRIP(1,I+1,J)
120                    DEXY(1,I,J) = ESTRIP(1,I,J)                    DEXY(1,I,J) = ESTRIP(1,I,J)
121    C
122                      ENPLANEX = ENPLANEX + DEXY(1,I,J)
123                    NSTRIP = NSTRIP + 1.                    NSTRIP = NSTRIP + 1.
124                    QTOT = QTOT + ESTRIP(1,I,J)                    QTOT = QTOT + DEXY(1,I,J)
125                    if (i.lt.11) QQ(2) = QQ(2) + ESTRIP(1,I,J)                    if (i.lt.11) QQ(2) = QQ(2) + DEXY(1,I,J)
126    c                  IF (I.EQ.21) THEN
127    c                     NX22 = NX22 + 1.
128    c                     QX22 = QX22 + DEXY(1,I,J)
129    c                  ENDIF
130                 ENDIF                 ENDIF
131              ENDIF              ENDIF
132              IF ( MOD(I,2).EQ.0 ) THEN              IF ( MOD(I,2).EQ.0 ) THEN
133                 IF (ESTRIP(2,I,J).GT.EMIN) THEN                 IF (ESTRIP(2,I,J).GT.EMIN) THEN
                   ENPLANEY = ENPLANEY + ESTRIP(2,I,J)  
134                    DEXY(2,I,J) = ESTRIP(2,I,J)                    DEXY(2,I,J) = ESTRIP(2,I,J)
135                      ENPLANEY = ENPLANEY + DEXY(2,I,J)
136                    NSTRIP = NSTRIP + 1.                    NSTRIP = NSTRIP + 1.
137                    QTOT = QTOT + ESTRIP(2,I,J)                    QTOT = QTOT + DEXY(2,I,J)
138                    if (i.lt.11) QQ(3) = QQ(3) + ESTRIP(2,I,J)                    if (i.lt.11) QQ(3) = QQ(3) + DEXY(2,I,J)
139                 ENDIF                 ENDIF
140                 IF (ESTRIP(1,I,J).GT.EMIN) THEN                 IF (ESTRIP(1,I,J).GT.EMIN) THEN
141                    ENPLANEX = ENPLANEX + ESTRIP(1,I,J)  C               IF (ESTRIP(1,I-1,J).GT.EMIN) THEN
142    C                  
143    c                  print *,' i ',i,' i - 1 ',I-1
144    C                  DEXY(1,I,J) = ESTRIP(1,I-1,J)
145                    DEXY(1,I,J) = ESTRIP(1,I,J)                    DEXY(1,I,J) = ESTRIP(1,I,J)
146    C
147                      ENPLANEX = ENPLANEX + DEXY(1,I,J)
148                    NSTRIP = NSTRIP + 1.                    NSTRIP = NSTRIP + 1.
149                    QTOT = QTOT + ESTRIP(1,I,J)                    QTOT = QTOT + DEXY(1,I,J)
150                    IF (I.EQ.22) THEN                    IF (I.EQ.22) THEN
151                       NX22 = NX22 + 1.                       NX22 = NX22 + 1.
152                       QX22 = QX22 + ESTRIP(1,I,J)                       QX22 = QX22 + DEXY(1,I,J)
153                    ENDIF                    ENDIF
154                    IF (I.LT.11) QQ(4) = QQ(4) + ESTRIP(1,I,J)                    IF (I.LT.11) QQ(4) = QQ(4) + DEXY(1,I,J)
155                 ENDIF                 ENDIF
156              ENDIF              ENDIF
157  C  C
# Line 170  C Line 185  C
185                    IMPX = (CX - XALIG)/10.                    IMPX = (CX - XALIG)/10.
186                    TANX = TG(1)                    TANX = TG(1)
187                 ELSE                 ELSE
188                    IMPY =( YALIG - CY)/10.                    IMPY =( CY - YALIG)/10.
189                    TANY = TG(2)                    TANY = TG(2)
190                 ENDIF                 ENDIF
191  C                C              
192                 DO I = 1,NPLA                     DO I = 1,NPLA    
193                    NN = 0                    NN = 1
194                    IF (M.EQ.2) NN = 1  C                  NN = 0
195    C                  IF (M.EQ.2) NN = 1
196                    IF (MOD(I,2).EQ.NN) THEN                    IF (MOD(I,2).EQ.NN) THEN
197                       SHIFT = +0.5                       SHIFT = +0.5
198                    ELSE                    ELSE
# Line 193  C     Line 209  C    
209                       DISTY = PIANO(I)                                       DISTY = PIANO(I)                
210                       YY(I) = DISTY * TG(2) + CY                       YY(I) = DISTY * TG(2) + CY
211                       BAR(M,I) = YY(I)                       BAR(M,I) = YY(I)
212                       CBAR(M,I) = (YALIG - YY(I))/10.                       CBAR(M,I) = (-YALIG + YY(I))/10.
213  C      C    
214                    ENDIF                    ENDIF
215                    CALL LASTRISCIA(BAR(M,I),IBAR(M,I))                    CALL LASTRISCIA(BAR(M,I),IBAR(M,I))
216                    CIBAR(M,I) = IBAR(M,I)                    CIBAR(M,I) = IBAR(M,I)
217                      IF (ibar(m,i).EQ.-1) THEN
218                         CHTRACK = CHTRACK + 1
219                      ELSE
220                         IWPL(M) = IWPL(M) + 1
221                      ENDIF
222                 ENDDO                             ENDDO            
223              ENDIF              ENDIF
224  C  C
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.23