| 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 |
|
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' |
PFA='COG4' |
| 44 |
|
TRACKMODE = 0 |
| 45 |
|
FACT = 100. |
| 46 |
|
ISTEPMIN = 3 |
| 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='COG' |
| 55 |
|
c$$$ PFA='ETA' |
| 56 |
|
c$$$ DEBUG=.true. |
| 57 |
|
c$$$ DEBUG=.false. |
| 58 |
|
|
| 59 |
|
if(DEBUG)PRINT*,'P.F.A. --> ',fin,PFA |
| 60 |
|
|
| 61 |
if(DEBUG)then |
if(DEBUG)then |
| 62 |
print*,'----------------------------------' |
print*,'----------------------------------' |
| 65 |
|
|
| 66 |
*------------------------------------------------------ |
*------------------------------------------------------ |
| 67 |
call init_level2 |
call init_level2 |
| 68 |
c call init_hough |
call init_hough |
| 69 |
*------------------------------------------------------ |
*------------------------------------------------------ |
| 70 |
|
|
| 71 |
*------------------------------------------------------ |
*------------------------------------------------------ |
| 93 |
goto 880 !fill ntp and go to next event |
goto 880 !fill ntp and go to next event |
| 94 |
endif |
endif |
| 95 |
|
|
| 96 |
c call fill_hough |
call fill_hough |
| 97 |
|
|
| 98 |
iflag=0 |
iflag=0 |
| 99 |
call track_fitting(iflag) |
call track_fitting(iflag) |