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

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

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

revision 1.21 by pam-fi, Tue Aug 7 13:56:29 2007 UTC revision 1.22 by pam-fi, Mon Aug 20 16:07:16 2007 UTC
# Line 27  c$$$      debug = .true. Line 27  c$$$      debug = .true.
27  c$$$      verbose = .true.  c$$$      verbose = .true.
28  c$$$      warning = .true.  c$$$      warning = .true.
29    
30    c$$$      print*,debug,verbose,warning
31    c$$$      debug=1
32    c$$$      verbose=1
33    c$$$      warning=1
34    
35  *     //////////////////////////  *     //////////////////////////
36  *     initialize some parameters  *     initialize some parameters
37  *     //////////////////////////  *     //////////////////////////
# Line 35  c$$$      warning = .true. Line 40  c$$$      warning = .true.
40    
41  c      debug=.true.  c      debug=.true.
42    
43        if(debug)print*,'-- check LEVEL0 status'        if(debug.eq.1)print*,'-- check LEVEL0 status'
44    
45        ievco=-1        ievco=-1
46        mismatch=0        mismatch=0
# Line 60  c               GOOD1(DSPnumber(iv)) = G Line 65  c               GOOD1(DSPnumber(iv)) = G
65                 GOOD1(DSPnumber(iv)) = ior(GOOD1(DSPnumber(iv)),2**1)                 GOOD1(DSPnumber(iv)) = ior(GOOD1(DSPnumber(iv)),2**1)
66   102           format(' * WARNING * Event ',i7,' view',i3   102           format(' * WARNING * Event ',i7,' view',i3
67       $          ,' CRC error')       $          ,' CRC error')
68                 if(debug)write(*,102)eventn(1),DSPnumber(iv)                 if(debug.eq.1)write(*,102)eventn(1),DSPnumber(iv)
69  c               goto 18 !next view  c               goto 18 !next view
70              endif              endif
71  c           ------------------------  c           ------------------------
# Line 81  c               GOOD1(DSPnumber(iv)) = G Line 86  c               GOOD1(DSPnumber(iv)) = G
86                 GOOD1(DSPnumber(iv)) = ior(GOOD1(DSPnumber(iv)),2**2)                 GOOD1(DSPnumber(iv)) = ior(GOOD1(DSPnumber(iv)),2**2)
87   103           format(' * WARNING * Event ',i7,' view',i3   103           format(' * WARNING * Event ',i7,' view',i3
88       $          ,' software alarm')       $          ,' software alarm')
89                 if(debug)write(*,103)eventn(1),DSPnumber(iv)                 if(debug.eq.1)write(*,103)eventn(1),DSPnumber(iv)
90  c               goto 18  c               goto 18
91              endif              endif
92  c           ------------------------  c           ------------------------
# Line 102  c$$$c                  GOOD1(DSPnumber(i Line 107  c$$$c                  GOOD1(DSPnumber(i
107  c$$$                  GOOD1(DSPnumber(iv)) = ior(GOOD1(DSPnumber(iv)),2**3)  c$$$                  GOOD1(DSPnumber(iv)) = ior(GOOD1(DSPnumber(iv)),2**3)
108  c$$$ 104              format(' * WARNING * Event ',i7,' view',i3  c$$$ 104              format(' * WARNING * Event ',i7,' view',i3
109  c$$$     $          ,' counter jump ',i10,i10)  c$$$     $          ,' counter jump ',i10,i10)
110  c$$$                  if(debug)write(*,104)eventn(1),DSPnumber(iv)  c$$$                  if(debug.eq.1)write(*,104)eventn(1),DSPnumber(iv)
111  c$$$     $                 ,eventn_old(iv),eventn(iv))  c$$$     $                 ,eventn_old(iv),eventn(iv))
112  c$$$                  goto 18  c$$$                  goto 18
113  c$$$               endif  c$$$               endif
# Line 126  c           ------------------------ Line 131  c           ------------------------
131    
132  c      print*,'*** ',(eventn(iv),iv=1,12)  c      print*,'*** ',(eventn(iv),iv=1,12)
133                
134        if(mismatch.eq.1.and.debug)        if(mismatch.eq.1.and.debug.eq.1)
135       $     print*,' * WARNING * DSP counter mismatch: '       $     print*,' * WARNING * DSP counter mismatch: '
136       $     ,(eventn(iv),iv=1,12)       $     ,(eventn(iv),iv=1,12)
137    
# Line 149  c     read the variable DATATRACKER from Line 154  c     read the variable DATATRACKER from
154  c     and fill the variable ADC (invertin view 11)  c     and fill the variable ADC (invertin view 11)
155  c--------------------------------------------------  c--------------------------------------------------
156                
157        if(debug)print*,'-- fill ADC vectors'        if(debug.eq.1)print*,'-- fill ADC vectors'
158    
159        call filladc(iflag)        call filladc(iflag)
160        if(iflag.ne.0)then        if(iflag.ne.0)then
# Line 161  c     computes common noise for each VA1 Line 166  c     computes common noise for each VA1
166  c     (excluding strips with signal,  c     (excluding strips with signal,
167  c     tagged with the flag CLSTR)  c     tagged with the flag CLSTR)
168  c--------------------------------------------------  c--------------------------------------------------
169        if(debug)print*,'-- compute CN'        if(debug.eq.1)print*,'-- compute CN'
170    
171        do iv=1,nviews        do iv=1,nviews
172           ima=0           ima=0
# Line 186  c--------------------------------------- Line 191  c---------------------------------------
191              endif              endif
192           enddo           enddo
193   100     format(' * WARNING * Event ',i7,' view',i3,': VK MASK ',24i1)   100     format(' * WARNING * Event ',i7,' view',i3,': VK MASK ',24i1)
194           if(ima.ne.0.and.verbose)write(*,100)eventn(1),iv           if(ima.ne.0.and.verbose.eq.1)write(*,100)eventn(1),iv
195       $        ,(mask_vk_ev(iv,ik),ik=1,nva1_view)       $        ,(mask_vk_ev(iv,ik),ik=1,nva1_view)
196  c         if(ima.ne.0)write(*,100)eventn(1),iv  c         if(ima.ne.0)write(*,100)eventn(1),iv
197  c     $        ,(mask_vk_ev(iv,ik),ik=1,nva1_view)    c     $        ,(mask_vk_ev(iv,ik),ik=1,nva1_view)  
# Line 202  c     sigma informations from histograms Line 207  c     sigma informations from histograms
207  c---------------------------------------------  c---------------------------------------------
208        ind=1                     !clsignal array index        ind=1                     !clsignal array index
209    
210        if(debug)print*,'-- search clusters'        if(debug.eq.1)print*,'-- search clusters'
211        do iv=1,nviews            !loop on views        do iv=1,nviews            !loop on views
212          do is=1,nstrips_view    !loop on strips (1)          do is=1,nstrips_view    !loop on strips (1)
213            if(mod(iv,2).eq.1) then            if(mod(iv,2).eq.1) then
# Line 236  C===  > X view Line 241  C===  > X view
241    
242          if(.not.flag_shower)then          if(.not.flag_shower)then
243             call save_cluster(iv)             call save_cluster(iv)
244             if(debug)print*,'view ',iv,' #clusters ', nclstr_view             if(debug.eq.1)print*,'view ',iv,' #clusters ', nclstr_view
245          else          else
246             fshower(iv) = 1             fshower(iv) = 1
247  c           GOOD1(DSPnumber(iv)) = 11 !AHAHAHAHA ORRORE!!  c           GOOD1(DSPnumber(iv)) = 11 !AHAHAHAHA ORRORE!!
# Line 245  c           GOOD1(iv) = GOOD1(iv) + 2**5 Line 250  c           GOOD1(iv) = GOOD1(iv) + 2**5
250             GOOD1(iv) = ior(GOOD1(iv),2**5)             GOOD1(iv) = ior(GOOD1(iv),2**5)
251   101       format(' * WARNING * Event ',i7,' view',i3   101       format(' * WARNING * Event ',i7,' view',i3
252       $          ,' #clusters > ',i5,' --> MASKED')       $          ,' #clusters > ',i5,' --> MASKED')
253             if(verbose)write(*,101)eventn(1),iv,nclstrmax_view             if(verbose.eq.1)write(*,101)eventn(1),iv,nclstrmax_view
254          endif          endif
255        enddo                     ! end loop on views        enddo                     ! end loop on views
256        do iv=1,nviews        do iv=1,nviews
# Line 265  C--------------------------------------- Line 270  C---------------------------------------
270        do iv = 1,nviews        do iv = 1,nviews
271           ngood = ngood + good1(iv)           ngood = ngood + good1(iv)
272        enddo        enddo
273        if(verbose.and.ngood.ne.0)print*,'* WARNING * Event ',eventn(1)        if(verbose.eq.1.and.ngood.ne.0)
274         $     print*,'* WARNING * Event ',eventn(1)
275       $     ,':LEVEL1 event status: '       $     ,':LEVEL1 event status: '
276       $     ,(good1(i),i=1,nviews)       $     ,(good1(i),i=1,nviews)
277  c------------------------------------------------------------------------  c------------------------------------------------------------------------
Legend:
Removed from v.1.21  
changed lines
  Added in v.1.22
  ViewVC Help
Powered by ViewVC 1.1.23