/[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.10 by mocchiut, Thu Nov 29 14:20:30 2007 UTC revision 1.14 by mocchiut, Thu Jan 16 15:29:09 2014 UTC
# Line 73  C Line 73  C
73        CALL VZERO(YY,NPLAV)              CALL VZERO(YY,NPLAV)      
74        VARFIT(1) = 0.        VARFIT(1) = 0.
75        VARFIT(2) = 0.        VARFIT(2) = 0.
76        NPFIT(1) = 0.        NPFIT(1) = 0
77        NPFIT(2) = 0.        NPFIT(2) = 0
78        VARCFIT(1) = 0.        VARCFIT(1) = 0.
79        VARCFIT(2) = 0.        VARCFIT(2) = 0.
80        NPCFIT(1) = 0.        NPCFIT(1) = 0
81        NPCFIT(2) = 0.        NPCFIT(2) = 0
82        ELEN = 0.        ELEN = 0.
83        SELEN = 0.        SELEN = 0.
84        NSTRIP = 0.        NSTRIP = 0.
# Line 202  C                  NN = 1 Line 202  C                  NN = 1
202  C                  IF (M.EQ.2) NN = 1  C                  IF (M.EQ.2) NN = 1
203                    IF (MOD(I,2).EQ.NN) THEN                    IF (MOD(I,2).EQ.NN) THEN
204                       IF (REVERSE.EQ.0) THEN                       IF (REVERSE.EQ.0) THEN
                         SHIFT = +0.5  
                      ELSE  
205                          SHIFT = -0.5                          SHIFT = -0.5
206                         ELSE
207                            SHIFT = +0.5
208                       ENDIF                       ENDIF
209                    ELSE                    ELSE
210                       IF (REVERSE.EQ.0) THEN                       IF (REVERSE.EQ.0) THEN
                         SHIFT = -0.5  
                      ELSE  
211                          SHIFT = +0.5                          SHIFT = +0.5
212                         ELSE
213                            SHIFT = -0.5
214                       ENDIF                       ENDIF
215                    ENDIF                    ENDIF
216  C      C    
# Line 247  C Line 247  C
247  C     CALCULATE ELEN AND SELEN  C     CALCULATE ELEN AND SELEN
248  C  C
249  c      ELEN =  40.82 * QTOT * 1.06 * 1E-4    c      ELEN =  40.82 * QTOT * 1.06 * 1E-4  
250        ELEN =  260. * QTOT    C      ELEN =  260. * QTOT  
251          ELEN =  QTOT / 260.  
252        IF (ELEN.GT.0.) THEN        IF (ELEN.GT.0.) THEN
253           SELEN = ELEN * (0.01183 + 0.121 / SQRT(ELEN))           SELEN = ELEN * (0.01183 + 0.121 / SQRT(ELEN))
254        ELSE        ELSE
# Line 256  c      ELEN =  40.82 * QTOT * 1.06 * 1E- Line 257  c      ELEN =  40.82 * QTOT * 1.06 * 1E-
257  C  C
258   9696 CONTINUE   9696 CONTINUE
259  C  C
260   45   continue  c 45   continue
261    
262   50   continue  c 50   continue
263    
264        return        return
265        END        END
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.23