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,NPLA) |
8 |
REAL RIL(NCHA/2) |
9 |
INTEGER IBAR(2,NPLA) |
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.22 ) NPIANI = 22 |
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 |
999 RETURN |
67 |
END |