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

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

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.2 - (show annotations) (download)
Fri Dec 12 16:18:18 2008 UTC (15 years, 11 months ago) by mocchiut
Branch: MAIN
CVS Tags: v10RED, v9r00, v9r01, v10REDr01, v6r01, v6r00, HEAD
Changes since 1.1: +18 -2 lines
Calo optimization + bug fixed, ToF errors -315/6/7 fixed

1 C
2 C------------------------------------------------------
3 c SUBROUTINE BARIC(Q,INF,ISUP,BAR,EPLA4,INDEX)
4 INTEGER FUNCTION BARIC(Q,INF,ISUP,BAR,EPLA4,INDEX)
5 C------------------------------------------------------
6 INCLUDE 'INTEST.TXT'
7 REAL Q(NCHA)
8 REAL BAR
9 REAL CM, EPLA4, ENN
10 INTEGER INF,ISUP,INDEX, LL
11 COMMON / SHIFT / SHIFT
12 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.
21 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
30 ENN = Q(LL)
31 IF (ENN.GT.EMIN) THEN
32 C
33 CALL MILLIM(LL,RMM)
34 C
35 CM = CM + RMM * ENN
36 EPLA4 = EPLA4 + ENN
37 ENDIF
38 ENDDO
39 IF (EPLA4.EQ.0) GO TO 300
40 BAR = CM / EPLA4
41 GOTO 400
42 C
43 300 INDEX = 0
44 GOTO 500
45 C
46 400 INDEX = 1
47 C
48 500 CONTINUE
49 RETURN
50 END
51 C

  ViewVC Help
Powered by ViewVC 1.1.23