--- DarthVader/TrackerLevel2/src/F77/analysisflight.f 2007/02/16 14:56:02 1.10 +++ DarthVader/TrackerLevel2/src/F77/analysisflight.f 2007/08/20 16:07:16 1.16 @@ -6,7 +6,7 @@ * - create LEVEL2 ntuple * ************************************************************************* - subroutine analysisflight(fin) + subroutine analysisflight include 'commontracker.f' include 'level1.f' @@ -19,9 +19,7 @@ * input flag * -* 0 - PFA = ETA -* 1 - PFA = COG4 - integer fin +c integer fin * flag to chose PFA character*10 PFA @@ -42,25 +40,34 @@ ************************************************************ ************************************************************ ************************************************************ - 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$$$ TRACKMODE = 0 +c$$$ FACT = 100. +c$$$ ISTEPMIN = 3 + + call idtoc(pfaid,PFA) + +c$$$ PFA='COG4' +c$$$ if(pfaid.eq.0)PFA='ETA' +c$$$ if(pfaid.eq.2)PFA='ETA2' +c$$$ if(pfaid.eq.3)PFA='ETA3' +c$$$ if(pfaid.eq.4)PFA='ETA4' +c$$$ if(pfaid.eq.10)PFA='COG' +c$$$ if(pfaid.eq.11)PFA='COG1' +c$$$ if(pfaid.eq.12)PFA='COG2' +c$$$ if(pfaid.eq.13)PFA='COG3' +c$$$ if(pfaid.eq.14)PFA='COG4' +*********************************************************** -c$$$ PFA='ETA' -c$$$ DEBUG=.true. +c if(DEBUG.EQ.1)PRINT*,'P.F.A. --> ',fin,PFA - if(DEBUG)PRINT*,'P.F.A. --> ',fin,PFA - - if(DEBUG)then + if(DEBUG.EQ.1)then print*,'----------------------------------' -c print*,'START... ',good1,nclstr1,nclstrmax_level2 + print*,'Settings: ' + print*,'PFA ',pfaid,' ---> ',PFA + print*,'tracking mode ',trackmode + print*,'fit-tolerance factor ',fact + print*,'minimum n.step ',istepmin + endif *------------------------------------------------------ @@ -79,10 +86,10 @@ cl_used(i)=0 !init mask of clusters associated to a track enddo - if(DEBUG)then - print*,'----------------------------------' - print*,iev,' ** ',nev2 - endif +c$$$ if(DEBUG.EQ.1)then +c$$$ print*,'----------------------------------' +c$$$ print*,iev,' ** ',nev2 +c$$$ endif * /////////////////////////////////////////////// * \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\ @@ -111,7 +118,7 @@ call fill_level2_siglets - if(DEBUG)then + if(DEBUG.EQ.1)then print*,'' print*,'DONE!'