/[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.1 - (show annotations) (download)
Fri May 19 13:15:51 2006 UTC (18 years, 6 months ago) by mocchiut
Branch: MAIN
Branch point for: DarthVader
Initial revision

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

  ViewVC Help
Powered by ViewVC 1.1.23