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

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

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.1 - (hide 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 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     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