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

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

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

revision 1.1 by mocchiut, Fri May 19 13:15:50 2006 UTC revision 1.2 by mocchiut, Fri Dec 12 16:18:18 2008 UTC
# Line 1  Line 1 
1  C  C
2  C------------------------------------------------------  C------------------------------------------------------
3        SUBROUTINE BARIC(Q,INF,ISUP,BAR,EPLA4,INDEX)  c      SUBROUTINE BARIC(Q,INF,ISUP,BAR,EPLA4,INDEX)
4          INTEGER FUNCTION BARIC(Q,INF,ISUP,BAR,EPLA4,INDEX)
5  C------------------------------------------------------  C------------------------------------------------------
6        INCLUDE 'INTEST.TXT'        INCLUDE 'INTEST.TXT'
7        REAL Q(NCHA)        REAL Q(NCHA)
8        REAL BAR        REAL BAR
9          REAL CM, EPLA4, ENN
10          INTEGER INF,ISUP,INDEX, LL
11        COMMON / SHIFT / SHIFT        COMMON / SHIFT / SHIFT
12        SAVE / SHIFT /        SAVE / SHIFT /
13          DATA CM / 0. /
14    c      DATA EPLA4 / 0. /
15          DATA ENN / 0. /
16    
17          BARIC = 0
18    
19          LL = 0
20        CM = 0.        CM = 0.
21        EPLA4 = 0.        EPLA4 = 0.
22          if ( index*Q(1) .eq. -2000.323311 ) then
23    c         index = 10
24             print *,' BEWARE!'
25             print *,' this is the GHOST of the compilation '
26             print *,' you have awakened me too soon '
27             print *,' I will launch an anatema against you!'
28          endif
29        DO LL = INF,ISUP        DO LL = INF,ISUP
30          ENN = Q(LL)          ENN = Q(LL)
31          IF (ENN.GT.EMIN) THEN          IF (ENN.GT.EMIN) THEN
32  C  C
33            CALL MILLIM(LL,RMM)            CALL MILLIM(LL,RMM)
34    C
35            CM = CM + RMM * ENN            CM = CM + RMM * ENN
36            EPLA4 = EPLA4 + ENN            EPLA4 = EPLA4 + ENN
37          ENDIF          ENDIF
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.23