6 |
* - create LEVEL2 ntuple |
* - create LEVEL2 ntuple |
7 |
* |
* |
8 |
************************************************************************* |
************************************************************************* |
9 |
subroutine analysisflight |
subroutine analysisflight(fin) |
10 |
|
|
11 |
include 'commontracker.f' |
include 'commontracker.f' |
12 |
include 'level1.f' |
include 'level1.f' |
17 |
include 'calib.f' |
include 'calib.f' |
18 |
include 'level2.f' |
include 'level2.f' |
19 |
|
|
20 |
|
* input flag |
21 |
|
* |
22 |
|
* 0 - PFA = ETA |
23 |
|
* 1 - PFA = COG4 |
24 |
|
integer fin |
25 |
|
|
26 |
* flag to chose PFA |
* flag to chose PFA |
27 |
character*10 PFA |
character*10 PFA |
28 |
common/FINALPFA/PFA |
common/FINALPFA/PFA |
43 |
************************************************************ |
************************************************************ |
44 |
************************************************************ |
************************************************************ |
45 |
PFA='ETA' |
PFA='ETA' |
46 |
|
TRACKMODE = 0 |
47 |
|
if(fin.eq.0)PFA='ETA' |
48 |
|
if(fin.eq.2)PFA='ETA2' |
49 |
|
if(fin.eq.3)PFA='ETA3' |
50 |
|
if(fin.eq.4)PFA='ETA4' |
51 |
|
if(fin.eq.10)PFA='COG' |
52 |
|
if(fin.eq.11)PFA='COG1' |
53 |
|
if(fin.eq.12)PFA='COG2' |
54 |
|
if(fin.eq.13)PFA='COG3' |
55 |
|
if(fin.eq.14)PFA='COG4' |
56 |
|
|
57 |
|
c$$$ PFA='ETA' |
58 |
|
c$$$ DEBUG=.true. |
59 |
|
|
60 |
|
if(DEBUG)PRINT*,'P.F.A. --> ',fin,PFA |
61 |
|
|
62 |
if(DEBUG)then |
if(DEBUG)then |
63 |
print*,'----------------------------------' |
print*,'----------------------------------' |
66 |
|
|
67 |
*------------------------------------------------------ |
*------------------------------------------------------ |
68 |
call init_level2 |
call init_level2 |
69 |
c call init_hough |
call init_hough |
70 |
*------------------------------------------------------ |
*------------------------------------------------------ |
71 |
|
|
72 |
*------------------------------------------------------ |
*------------------------------------------------------ |
94 |
goto 880 !fill ntp and go to next event |
goto 880 !fill ntp and go to next event |
95 |
endif |
endif |
96 |
|
|
97 |
c call fill_hough |
call fill_hough |
98 |
|
|
99 |
iflag=0 |
iflag=0 |
100 |
call track_fitting(iflag) |
call track_fitting(iflag) |