/[PAMELA software]/DarthVader/TrackerLevel2/src/F77/analysisflight.f
ViewVC logotype

Diff of /DarthVader/TrackerLevel2/src/F77/analysisflight.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.1 by mocchiut, Fri May 19 13:15:54 2006 UTC revision 1.10 by pam-fi, Fri Feb 16 14:56:02 2007 UTC
# Line 6  Line 6 
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'
13        include 'common_momanhough.f'        include 'common_momanhough.f'
14        include 'common_mech.f'        include 'common_mech.f'
15        include 'common_xyzPAM.f'        include 'common_xyzPAM.f'
16        include 'common_mini_2.f'        include 'common_mini_2.f'
17        include 'calib.f'        include 'calib.f'
       include 'level1.f'  
18        include 'level2.f'        include 'level2.f'
19    
20  c      include 'momanhough_init.f'  *     input flag
21    *
22  *     flag to set debug mode  *     0 - PFA = ETA
23        logical DEBUG  *     1 - PFA = COG4
24        common/dbg/DEBUG        integer fin
25    
26  *     flag to chose PFA  *     flag to chose PFA
27        character*10 PFA        character*10 PFA
# Line 33  c      parameter (inf=1.e8)      !just a Line 33  c      parameter (inf=1.e8)      !just a
33        external npl        external npl
34        external acoordsi,coordsi,nld,coord,dcoord        external acoordsi,coordsi,nld,coord,dcoord
35    
       DEBUG=.false.!.true.!  
36  ************************************************************  ************************************************************
37  ************************************************************  ************************************************************
38  ************************************************************  ************************************************************
# Line 44  c      parameter (inf=1.e8)      !just a Line 43  c      parameter (inf=1.e8)      !just a
43  ************************************************************  ************************************************************
44  ************************************************************  ************************************************************
45        PFA='ETA'        PFA='ETA'
46          if(fin.eq.0)PFA='ETA'
47          if(fin.eq.2)PFA='ETA2'
48          if(fin.eq.3)PFA='ETA3'
49          if(fin.eq.4)PFA='ETA4'
50          if(fin.eq.10)PFA='COG'
51          if(fin.eq.11)PFA='COG1'
52          if(fin.eq.12)PFA='COG2'
53          if(fin.eq.13)PFA='COG3'
54          if(fin.eq.14)PFA='COG4'
55    
56    c$$$      PFA='ETA'
57    c$$$      DEBUG=.true.
58    
59          if(DEBUG)PRINT*,'P.F.A. --> ',fin,PFA
60    
61        if(DEBUG)then        if(DEBUG)then
62           print*,'----------------------------------'           print*,'----------------------------------'
63           print*,'START... ',good1,nclstr1,nclstrmax_level2  c         print*,'START... ',good1,nclstr1,nclstrmax_level2
64        endif        endif
65    
66  *------------------------------------------------------  *------------------------------------------------------
 *     LEVEL2 N-TUPLE INITIALIZATIONS  
67        call init_level2        call init_level2
68  c      if(.not.good1)then        call init_hough
       if(good1.eq.0)then  
          goto 8800              !fill nt-uple and go to next event  
       endif  
69  *------------------------------------------------------  *------------------------------------------------------
70                
71  *------------------------------------------------------  *------------------------------------------------------
72  *     cut on maximum number of clusters  *     cut on maximum number of clusters
73  *------------------------------------------------------  *------------------------------------------------------
74        if(nclstr1.gt.nclstrmax_level2)then  c$$$      if(nclstr1.gt.nclstrmax_level2)then
75           goto 8800              !fill nt-uple and go to next event  c$$$         goto 8800              !fill nt-uple and go to next event
76        endif  c$$$      endif
77                
78        do i=1,nclstr1        do i=1,nclstr1
79           cl_used(i)=0           !init mask of clusters associated to a track           cl_used(i)=0           !init mask of clusters associated to a track
# Line 84  c      if(.not.good1)then Line 93  c      if(.not.good1)then
93           goto 880               !fill ntp and go to next event               goto 880               !fill ntp and go to next event    
94        endif        endif
95    
96          call fill_hough
97    
98        iflag=0        iflag=0
99        call track_fitting(iflag)        call track_fitting(iflag)
100        if(iflag.eq.1)then        !bad event        if(iflag.eq.1)then        !bad event
# Line 111  c      if(.not.good1)then Line 122  c      if(.not.good1)then
122           print*,''           print*,''
123           print*,''           print*,''
124        endif        endif
125    
126          ngood = 0
127          do iv = 1,nviews
128             ngood = ngood + good1(iv)
129          enddo
130    c$$$      if(ngood.ne.0)print*,'* WARNING * Event '
131    c$$$     $     ,':LEVEL2 event status: '
132    c$$$     $     ,(good2(i),i=1,nviews)
133                
134   8800 continue   8800 continue
135                
# Line 127  c      if(.not.good1)then Line 146  c      if(.not.good1)then
146  ************************************************************  ************************************************************
147    
148    
 c# include "momanhough-subroutines.f"  
   
   
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.23