17 |
include 'calib.f' |
include 'calib.f' |
18 |
include 'level2.f' |
include 'level2.f' |
19 |
|
|
20 |
|
* input flag |
21 |
|
* |
22 |
|
c integer fin |
23 |
|
|
24 |
* flag to chose PFA |
* flag to chose PFA |
25 |
character*10 PFA |
character*10 PFA |
26 |
common/FINALPFA/PFA |
common/FINALPFA/PFA |
40 |
************************************************************ |
************************************************************ |
41 |
************************************************************ |
************************************************************ |
42 |
************************************************************ |
************************************************************ |
43 |
PFA='ETA' |
c$$$ TRACKMODE = 0 |
44 |
|
c$$$ FACT = 100. |
45 |
|
c$$$ ISTEPMIN = 3 |
46 |
|
|
47 |
|
call idtoc(pfaid,PFA) |
48 |
|
|
49 |
|
c$$$ PFA='COG4' |
50 |
|
c$$$ if(pfaid.eq.0)PFA='ETA' |
51 |
|
c$$$ if(pfaid.eq.2)PFA='ETA2' |
52 |
|
c$$$ if(pfaid.eq.3)PFA='ETA3' |
53 |
|
c$$$ if(pfaid.eq.4)PFA='ETA4' |
54 |
|
c$$$ if(pfaid.eq.10)PFA='COG' |
55 |
|
c$$$ if(pfaid.eq.11)PFA='COG1' |
56 |
|
c$$$ if(pfaid.eq.12)PFA='COG2' |
57 |
|
c$$$ if(pfaid.eq.13)PFA='COG3' |
58 |
|
c$$$ if(pfaid.eq.14)PFA='COG4' |
59 |
|
*********************************************************** |
60 |
|
|
61 |
|
c if(DEBUG)PRINT*,'P.F.A. --> ',fin,PFA |
62 |
|
|
63 |
if(DEBUG)then |
if(DEBUG)then |
64 |
print*,'----------------------------------' |
print*,'----------------------------------' |
65 |
c print*,'START... ',good1,nclstr1,nclstrmax_level2 |
print*,'Settings: ' |
66 |
|
print*,'PFA ',pfaid,PFA |
67 |
|
print*,'tracking mode ',trackmode |
68 |
|
print*,'fit-tolerance factor ',fact |
69 |
|
print*,'minimum n.step ',istepmin |
70 |
endif |
endif |
71 |
|
|
72 |
*------------------------------------------------------ |
*------------------------------------------------------ |
73 |
call init_level2 |
call init_level2 |
74 |
|
call init_hough |
75 |
*------------------------------------------------------ |
*------------------------------------------------------ |
76 |
|
|
77 |
*------------------------------------------------------ |
*------------------------------------------------------ |
99 |
goto 880 !fill ntp and go to next event |
goto 880 !fill ntp and go to next event |
100 |
endif |
endif |
101 |
|
|
102 |
|
call fill_hough |
103 |
|
|
104 |
iflag=0 |
iflag=0 |
105 |
call track_fitting(iflag) |
call track_fitting(iflag) |
106 |
if(iflag.eq.1)then !bad event |
if(iflag.eq.1)then !bad event |
128 |
print*,'' |
print*,'' |
129 |
print*,'' |
print*,'' |
130 |
endif |
endif |
131 |
|
|
132 |
|
ngood = 0 |
133 |
|
do iv = 1,nviews |
134 |
|
ngood = ngood + good1(iv) |
135 |
|
enddo |
136 |
|
c$$$ if(ngood.ne.0)print*,'* WARNING * Event ' |
137 |
|
c$$$ $ ,':LEVEL2 event status: ' |
138 |
|
c$$$ $ ,(good2(i),i=1,nviews) |
139 |
|
|
140 |
8800 continue |
8800 continue |
141 |
|
|