1 |
mocchiut |
1.1 |
SUBROUTINE NOINTER(IGOOD) |
2 |
|
|
C- |
3 |
|
|
C- Created 20-OCT-1994 MIRKO BOEZIO |
4 |
|
|
C- |
5 |
|
|
INCLUDE 'INTEST.TXT' |
6 |
|
|
C |
7 |
mocchiut |
1.2 |
REAL BAR(2,NPLAV) |
8 |
mocchiut |
1.1 |
REAL RIL(NCHA/2) |
9 |
mocchiut |
1.2 |
INTEGER IBAR(2,NPLAV) |
10 |
mocchiut |
1.1 |
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 |
mocchiut |
1.2 |
IF ( NPIANI.gt.NPLA ) NPIANI = NPLA |
23 |
mocchiut |
1.1 |
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 |
mocchiut |
1.3 |
c 999 RETURN |
67 |
|
|
RETURN |
68 |
mocchiut |
1.1 |
END |