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

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

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.3 - (show annotations) (download)
Tue Aug 4 14:01:21 2009 UTC (15 years, 3 months ago) by mocchiut
Branch: MAIN
CVS Tags: v10RED, v9r00, v9r01, v10REDr01, HEAD
Changes since 1.2: +2 -1 lines
Changed to work with GCC 4.x (gfortran) + ROOT >= 5.24

1 SUBROUTINE NOINTER(IGOOD)
2 C-
3 C- Created 20-OCT-1994 MIRKO BOEZIO
4 C-
5 INCLUDE 'INTEST.TXT'
6 C
7 REAL BAR(2,NPLAV)
8 REAL RIL(NCHA/2)
9 INTEGER IBAR(2,NPLAV)
10 INTEGER NPIANI
11 COMMON/ ANGOLO / BAR, IBAR
12 SAVE / ANGOLO /
13 COMMON / GENERAL / RIG, RMASS
14 SAVE / GENERAL /
15 C
16 NPIANI = 0
17 C
18 EMAX = 2. * (RIG * RIG + RMASS * RMASS) / (RIG * RIG)
19 C
20 IF (RIG.NE.0) NPIANI = INT( ANINT( 9. + LOG(ABS(RIG)) ) )
21 C
22 IF ( NPIANI.gt.NPLA ) NPIANI = NPLA
23 IF ( NPIANI.lt.1 ) NPIANI = 1
24 C
25 DO K = 1,NCHA/2
26 RIL(K) = 1000.
27 ENDDO
28 C
29 NGIR = 0
30 IBAD = 0
31 IGOOD = 0
32 DO M = 1,2
33 DO J = 1,NPIANI
34 DO K = 1,NGIR
35 RIL(K) = 1000.
36 ENDDO
37 C
38 NNN = 0
39 BB = BAR(M,J)
40 C
41 C EM: IF TRACK EXIT THE CALORIMETER DO NOT CONSIDER THAT PLANE
42 C
43 IF (BB.LT.-6..OR.BB.GT.246.) GOTO 5
44 C
45 C END EM
46 C
47 NGIR = NCL(M,J)
48 IF (NGIR.EQ.0) GOTO 5
49 IF (NGIR.GT.1) IBAD = IBAD + 1
50 IF (NGIR.GT.1) GOTO 5
51 DO L = 1,NGIR
52 RIL(L) = ABS(CLUS(M,J,L) - BB)
53 ENDDO
54 MM = NCHA/2
55 LK = LVMIN(RIL,MM)
56 C
57 IF (RIL(LK).LE.4.AND.CLUS(M,J,LK+NCHA/2).LT.EMAX)
58 & IGOOD = IGOOD + J
59 C
60 5 CONTINUE
61 C
62 ENDDO
63 ENDDO
64 C
65
66 c 999 RETURN
67 RETURN
68 END

  ViewVC Help
Powered by ViewVC 1.1.23