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

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

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

revision 1.1 by mocchiut, Fri May 19 13:15:55 2006 UTC revision 1.6 by pam-fi, Thu Oct 26 16:22:37 2006 UTC
# Line 7  Line 7 
7  *            *          
8  *************************************************************************  *************************************************************************
9    
10        subroutine cncomp(i,j)    !(view, VA1)        subroutine cncomp(i,j,errflag)    !(view, VA1)
11    
12        include 'commontracker.f'        include 'commontracker.f'
13          include 'level1.f'
14        include 'common_reduction.f'        include 'common_reduction.f'
15        include 'calib.f'        include 'calib.f'
16    
17        integer errflag           !error flag to mark no signal free VA1        integer errflag           !error flag to mark no signal free VA1
   
18        integer clstr_old(nstrips_va1) !flag storage vector        integer clstr_old(nstrips_va1) !flag storage vector
   
19        real signal(nstrips_va1)  !"signal" (=adc-ped) value storage vector        real signal(nstrips_va1)  !"signal" (=adc-ped) value storage vector
   
20        real smean, ssigma        !"signal" mean and sigma        real smean, ssigma        !"signal" mean and sigma
21        real cut                  !"strange" strip exclusion cut        real cut                  !"strange" strip exclusion cut
   
22        integer newclstr          !flag to warn about new found clusters to be        integer newclstr          !flag to warn about new found clusters to be
23                                  ! excluded from common noise computation  c                               ! excluded from common noise computation
   
   
 c     call HBOOK1(20000+100*i+j,' ',30,0.,30.,0.) !???  
24    
25  c------------------------------------------------------------------------  c------------------------------------------------------------------------
26  c      c
27  c     variables initialization  c     variables initialization
28  c      c
29  c------------------------------------------------------------------------  c------------------------------------------------------------------------
30        do k=1,nstrips_va1        !loops on strips        do k=1,nstrips_va1        !loops on strips
31           clstr(i,j,k)=1         !initializes signal affected strips flag           clstr(i,j,k)=1         !initializes signal affected strips flag
# Line 40  c--------------------------------------- Line 34  c---------------------------------------
34        enddo                     ! affected strips flag        enddo                     ! affected strips flag
35    
36        newclstr=1                !flag to warn about new found signal        newclstr=1                !flag to warn about new found signal
37                                  ! affected strips  c                               ! affected strips
   
   
   
38  c------------------------------------------------------------------------  c------------------------------------------------------------------------
39  c      c
40  c     high or low signal affected strips exclusion: computes "signal" (=adc-ped)  c     high or low signal affected strips exclusion: computes "signal" (=adc-ped)
41  c     mean value and sigma, and cuts from common noise computation strips  c     mean value and sigma, and cuts from common noise computation strips
42  c     whose ABS(signal) exceeds scut*sigma  c     whose ABS(signal) exceeds scut*sigma
43  c      c
44  c------------------------------------------------------------------------  c------------------------------------------------------------------------
45        countme=0                 !???        countme=0                 !???
46   666  continue                  !???   666  continue                  !???
# Line 59  c--------------------------------------- Line 50  c---------------------------------------
50        nstr=0        nstr=0
51                
52        do k=1,nstrips_va1        do k=1,nstrips_va1
53           nstr=nstr+strange(i,j,k) !uses only           nstr = nstr + strange(i,j,k) !uses only
54           if(mod(i,2).eq.1) then !odd strip ---> Y view           if(mod(i,2).eq.1) then ! ---> Y view
55              signal(k)= - (DBLE(adc(i,j,k))-pedestal(i,j,k)) !negative signal              signal(k) = - (DBLE(adc(i,j,k)) - pedestal(i,j,k)) !negative signal
56           else                   !even strip ---> X view           else                   ! ---> X view
57              signal(k)= DBLE(adc(i,j,k))-pedestal(i,j,k) !positive signal              signal(k) =    DBLE(adc(i,j,k)) - pedestal(i,j,k) !positive signal
58           endif           endif
59                     smean = smean + signal(k)*strange(i,j,k)
60           smean=smean+signal(k)*strange(i,j,k)           ssigma = ssigma + (signal(k)**2)*strange(i,j,k)
          ssigma=ssigma+(signal(k)**2)*strange(i,j,k)  
           
 c     call HFILL(10000+100*i+j,signal(k),0.,1.) !???  
61        enddo        enddo
62                
63        smean=smean/nstr          !strips value distribution mean        smean=smean/nstr          !strips value distribution mean      
         
64        ssigma=SQRT((ssigma/nstr)-smean**2) !strips value distribution sigma        ssigma=SQRT((ssigma/nstr)-smean**2) !strips value distribution sigma
65                
66        cut=scut*ssigma           !exclusion cut        cut=scut*ssigma           !exclusion cut
67                
68          nco=0
69          nbo=0
70        do k=1,nstrips_va1        do k=1,nstrips_va1
 c     call HFILL(20000+100*i+j,ABS(signal(k)-smean)/ssigma,0.,1.) !???  
71           if(ABS(signal(k)-smean).gt.cut) then           if(ABS(signal(k)-smean).gt.cut) then
 c     print*,i,j,k,signal(k),abs(signal(k)),cut,strange(i,j,k) !???  
72              strange(i,j,k)=0    !marks strips exceeding cut              strange(i,j,k)=0    !marks strips exceeding cut
73    c            print*,i,j,k,signal(k),smean
74           endif           endif
75             nco=nco+strange(i,j,k)
76             nbo=nbo+bad(i,j,k)
77        enddo                     ! in order not to use them in CN computation        enddo                     ! in order not to use them in CN computation
78    
79    c$$$      if(i.eq.12.and.(j.eq.2.or.j.eq.3))then
80    c$$$         print*,'view ',i,' vk ',j
81    c$$$         print*,'ADC (1-51-128) = ',adc(i,j,1),adc(i,j,52),adc(i,j,128)
82    c$$$         print*,'<ADC-PED> = ',smean
83    c$$$         print*,'s         = ',ssigma
84    c$$$         print*,'nstrange  = ',128-nco
85    c$$$         print*,'nbad      = ',128-nbo
86    c$$$      endif
87    
88        countme=countme+1         !???        countme = countme + 1         !???
89        if (countme.le.3) goto 666 !???        if (countme.le.3) goto 666 !???
90    
   
91  c------------------------------------------------------------------------  c------------------------------------------------------------------------
92  c      c    
93  c     common noise computation  c     common noise computation
94  c      c    
95  c------------------------------------------------------------------------  c-----------------------------------------------------------------------
96        do while(newclstr.eq.1)   !loops on this VA1 till no new signal  *     loops on this VA1 till no new signal affected strips are found
97                                  ! affected strips are found        do while(newclstr.eq.1)  
98    
99           newclstr=0             !to exit from loop if no new cluster is           newclstr=0             !to exit from loop if no new cluster is found
                                 ! found  
100                    
101           errflag=0           errflag=0
102             call cnoise(i,j,errflag) !(view, VA1, error flag) computes cn
103           call cnoise(i,j,errflag) !(view, VA1, error flag) computes common           if(errflag.eq.1) goto 10 !goes to next VA1: this one has no signal-free strips...
                                 ! noise  
   
 c     print*,cn(i,j)         !???  
104                    
105           if(errflag.eq.1) goto 10 !goes to next VA1: this one has no signal           call cutcn(i,j)        !(view, VA1) excludes clusters from cn computation
                                 ! free strips...  
           
          call cutcn(i,j)        !(view, VA1) excludes clusters from  
                                 ! common noise calculation  
106                    
107           ncs=0                  !initializes number of strips not excluded by cncut           ncs=0                  !initializes number of strips not excluded by cncut
   
108           do k=1,nstrips_va1     !loops on strips           do k=1,nstrips_va1     !loops on strips
109              if(clstr(i,j,k).ne.clstr_old(k)) then !checks if there are  *           checks if there are new found clusters, and if so sets
110                                  ! new found clusters, and if so sets              if(clstr(i,j,k).ne.clstr_old(k)) then
111                 newclstr=1       ! newclstr flag = 1                 newclstr=1                      
112                                 clstr_old(k)=clstr(i,j,k)  
113                 clstr_old(k)=clstr(i,j,k) !stores cluster flags in              endif              
             endif               ! clstr_old variable  
   
114              iok=strange(i,j,k)*bad(i,j,k)*clstr(i,j,k)              iok=strange(i,j,k)*bad(i,j,k)*clstr(i,j,k)
   
115              ncs=ncs+iok         !counts number of good strips for cn computation              ncs=ncs+iok         !counts number of good strips for cn computation
   
116           enddo           enddo
117    
118        enddo                     !ends do while loop when there are no new        enddo                     !ends do while
                                 ! clusters  
   
 c      call HFILL(666,FLOAT(ncs),0.,1.) !???  
   
   
 c$$$      if(ncs.lt.20) then        !warns if too many strips have been excluded from CN  
 c$$$                                ! computation  
 c$$$         print*,'cncomp: WARNING, LESS THAN 20 STRIPS PASSED CN CUT'  
 c$$$     $        //' ON VA1 ',j,', VIEW ',i !NB questo errore e' "un po'" in conflitto  
 c$$$                                ! con quello che setta errflag (vedi cnoise.f)...  
 c$$$  
 c$$$      endif  
119    
120   10   continue   10   continue
121    
# Line 174  c$$$      endif Line 146  c$$$      endif
146        subroutine cnoise(i,j,gulp) !(view, VA1)        subroutine cnoise(i,j,gulp) !(view, VA1)
147    
148        include 'commontracker.f'        include 'commontracker.f'
149          include 'level0.f'
150          include 'level1.f'
151        include 'common_reduction.f'        include 'common_reduction.f'
152        include 'calib.f'        include 'calib.f'
153                
# Line 182  c$$$      endif Line 156  c$$$      endif
156                
157        ncn=0                     !number of strips in cn computation        ncn=0                     !number of strips in cn computation
158        cn(i,j)=0                 !initializes cn variable        cn(i,j)=0                 !initializes cn variable
159          cnrms(i,j)=0              !initializes cn rms
160          cnn(i,j)=0                !initialize cn flag
161    
162        do k=1,nstrips_va1        !loops on strips        do k=1,nstrips_va1        !loops on strips
163           iok=strange(i,j,k)*bad(i,j,k)*clstr(i,j,k) !flag to mark strange, bad  *        tags strange, bad or signal-affected strips
164                                  ! or signal affected strips           iok = strange(i,j,k)*bad(i,j,k)*clstr(i,j,k)
165  ccc         print*,i,j,k,strange(i,j,k),bad(i,j,k),clstr(i,j,k),iok !???           cn(i,j) = cn(i,j) + (DBLE(adc(i,j,k)) - pedestal(i,j,k))*iok
166             cnrms(i,j) = cnrms(i,j)
167           cn(i,j)=cn(i,j) + (DBLE(adc(i,j,k)) - pedestal(i,j,k))*iok !sums ADC-PED       $        + (DBLE(adc(i,j,k)) - pedestal(i,j,k))
168                                  ! values to compute common noise       $        *(DBLE(adc(i,j,k)) - pedestal(i,j,k))*iok
169           ncn = ncn + iok            !counts number of strips in cn computation           ncn = ncn + iok            !counts number of strips in cn computation
170        enddo        enddo
171          
172  ccc      print*,'ncn= ',ncn        if(ncn.lt.NSTRIPMIN) then         !no signal free strips on this VA1...
173        if(ncn.eq.0) then         !no signal free strips on this VA1...           if(ncn.eq.0)then
174           print*,'cnoise: WARNING, NO SIGNAL FREE STRIPS ON VA1 ',j,              if(debug)print*,' WARNING - cnoise: ',
175       $        ', VIEW ',i       $        'no strips for CN computation on VA1 ',j,
176         $        ', VIEW ',i,'  >>> FAILED '
177             else
178                if(debug)print*,' WARNING - cnoise: ',
179         $        'less than ',NSTRIPMIN
180         $           ,' strips for CN computation on VA1 ',j,
181         $        ', VIEW ',i,'  >>> FAILED '
182             endif
183           gulp=1           gulp=1
184             cnn(i,j) = 0
185        else        else
186           cn(i,j)=cn(i,j)/DBLE(ncn) !computes common noise           cn(i,j)=cn(i,j)/DBLE(ncn) !<<<< computes common noise
187           gulp=0                 !resets error flag           cnrms(i,j)= SQRT( cnrms(i,j)/DBLE(ncn) - cn(i,j)**2 )
188             cnn(i,j) = ncn
189             gulp=0                
190             if(debug.and.ABS(cn(i,j)).gt.1000)
191         $        print*,'Event ',eventn(1)
192         $        ,': cn(',i,',',j,')= ',cn(i,j),' ncn ',ncn
193        endif        endif
194    
195        return        return
# Line 227  ccc      print*,'ncn= ',ncn Line 216  ccc      print*,'ncn= ',ncn
216        subroutine cutcn(i,j)     !(view, VA1)        subroutine cutcn(i,j)     !(view, VA1)
217    
218        include 'commontracker.f'        include 'commontracker.f'
219          include 'level1.f'
220        include 'common_reduction.f'        include 'common_reduction.f'
221        include 'calib.f'        include 'calib.f'
222    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.23