--- DarthVader/TrackerLevel2/src/F77/analysisflight.f 2006/11/07 15:55:11 1.5 +++ DarthVader/TrackerLevel2/src/F77/analysisflight.f 2007/02/16 14:56:02 1.10 @@ -6,7 +6,7 @@ * - create LEVEL2 ntuple * ************************************************************************* - subroutine analysisflight + subroutine analysisflight(fin) include 'commontracker.f' include 'level1.f' @@ -17,6 +17,12 @@ include 'calib.f' include 'level2.f' +* input flag +* +* 0 - PFA = ETA +* 1 - PFA = COG4 + integer fin + * flag to chose PFA character*10 PFA common/FINALPFA/PFA @@ -37,6 +43,20 @@ ************************************************************ ************************************************************ PFA='ETA' + if(fin.eq.0)PFA='ETA' + if(fin.eq.2)PFA='ETA2' + if(fin.eq.3)PFA='ETA3' + if(fin.eq.4)PFA='ETA4' + if(fin.eq.10)PFA='COG' + if(fin.eq.11)PFA='COG1' + if(fin.eq.12)PFA='COG2' + if(fin.eq.13)PFA='COG3' + if(fin.eq.14)PFA='COG4' + +c$$$ PFA='ETA' +c$$$ DEBUG=.true. + + if(DEBUG)PRINT*,'P.F.A. --> ',fin,PFA if(DEBUG)then print*,'----------------------------------' @@ -45,6 +65,7 @@ *------------------------------------------------------ call init_level2 + call init_hough *------------------------------------------------------ *------------------------------------------------------ @@ -72,6 +93,8 @@ goto 880 !fill ntp and go to next event endif + call fill_hough + iflag=0 call track_fitting(iflag) if(iflag.eq.1)then !bad event @@ -104,9 +127,9 @@ do iv = 1,nviews ngood = ngood + good1(iv) enddo - if(ngood.ne.0)print*,'* WARNING * Event ' - $ ,':LEVEL2 event status: ' - $ ,(good2(i),i=1,nviews) +c$$$ if(ngood.ne.0)print*,'* WARNING * Event ' +c$$$ $ ,':LEVEL2 event status: ' +c$$$ $ ,(good2(i),i=1,nviews) 8800 continue