/[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.3 by pam-fi, Thu Sep 28 14:04:40 2006 UTC revision 1.16 by pam-fi, Mon Aug 20 16:07:16 2007 UTC
# Line 9  Line 9 
9        subroutine analysisflight        subroutine analysisflight
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    *     input flag
21    *
22    c      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
# Line 36  c      parameter (inf=1.e8)      !just a Line 40  c      parameter (inf=1.e8)      !just a
40  ************************************************************  ************************************************************
41  ************************************************************  ************************************************************
42  ************************************************************  ************************************************************
43        PFA='ETA'  c$$$      TRACKMODE = 0
44    c$$$      FACT = 100.
45    c$$$      ISTEPMIN = 3
46    
47          call idtoc(pfaid,PFA)
48          
49    c$$$      PFA='COG4'
50    c$$$      if(pfaid.eq.0)PFA='ETA'
51    c$$$      if(pfaid.eq.2)PFA='ETA2'
52    c$$$      if(pfaid.eq.3)PFA='ETA3'
53    c$$$      if(pfaid.eq.4)PFA='ETA4'
54    c$$$      if(pfaid.eq.10)PFA='COG'
55    c$$$      if(pfaid.eq.11)PFA='COG1'
56    c$$$      if(pfaid.eq.12)PFA='COG2'
57    c$$$      if(pfaid.eq.13)PFA='COG3'
58    c$$$      if(pfaid.eq.14)PFA='COG4'
59    ***********************************************************
60    
61    c      if(DEBUG.EQ.1)PRINT*,'P.F.A. --> ',fin,PFA
62    
63        if(DEBUG)then        if(DEBUG.EQ.1)then
64           print*,'----------------------------------'           print*,'----------------------------------'
65  c         print*,'START... ',good1,nclstr1,nclstrmax_level2           print*,'Settings: '
66             print*,'PFA                  ',pfaid,' ---> ',PFA
67             print*,'tracking mode        ',trackmode
68             print*,'fit-tolerance factor ',fact
69             print*,'minimum n.step       ',istepmin
70            
71        endif        endif
72    
73  *------------------------------------------------------  *------------------------------------------------------
 *     LEVEL2 N-TUPLE INITIALIZATIONS  
74        call init_level2        call init_level2
75  c      if(.not.good1)then        call init_hough
 c      if(good1.eq.0)then  
 c         goto 8800              !fill nt-uple and go to next event  
 c      endif  
76  *------------------------------------------------------  *------------------------------------------------------
77                
78  *------------------------------------------------------  *------------------------------------------------------
79  *     cut on maximum number of clusters  *     cut on maximum number of clusters
80  *------------------------------------------------------  *------------------------------------------------------
81        if(nclstr1.gt.nclstrmax_level2)then  c$$$      if(nclstr1.gt.nclstrmax_level2)then
82           goto 8800              !fill nt-uple and go to next event  c$$$         goto 8800              !fill nt-uple and go to next event
83        endif  c$$$      endif
84                
85        do i=1,nclstr1        do i=1,nclstr1
86           cl_used(i)=0           !init mask of clusters associated to a track           cl_used(i)=0           !init mask of clusters associated to a track
87        enddo        enddo
88                
89        if(DEBUG)then  c$$$      if(DEBUG.EQ.1)then
90           print*,'----------------------------------'  c$$$         print*,'----------------------------------'
91           print*,iev,'   ** ',nev2  c$$$         print*,iev,'   ** ',nev2
92        endif  c$$$      endif
93    
94  *     ///////////////////////////////////////////////  *     ///////////////////////////////////////////////
95  *     \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\  *     \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\
# Line 77  c      endif Line 100  c      endif
100           goto 880               !fill ntp and go to next event               goto 880               !fill ntp and go to next event    
101        endif        endif
102    
103          call fill_hough
104    
105        iflag=0        iflag=0
106        call track_fitting(iflag)        call track_fitting(iflag)
107        if(iflag.eq.1)then        !bad event        if(iflag.eq.1)then        !bad event
# Line 93  c      endif Line 118  c      endif
118    
119        call fill_level2_siglets        call fill_level2_siglets
120                
121        if(DEBUG)then        if(DEBUG.EQ.1)then
122                            
123           print*,''           print*,''
124           print*,'DONE!'           print*,'DONE!'
# Line 104  c      endif Line 129  c      endif
129           print*,''           print*,''
130           print*,''           print*,''
131        endif        endif
132    
133          ngood = 0
134          do iv = 1,nviews
135             ngood = ngood + good1(iv)
136          enddo
137    c$$$      if(ngood.ne.0)print*,'* WARNING * Event '
138    c$$$     $     ,':LEVEL2 event status: '
139    c$$$     $     ,(good2(i),i=1,nviews)
140                
141   8800 continue   8800 continue
142                
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.16
  ViewVC Help
Powered by ViewVC 1.1.23