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

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

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

revision 1.7 by pam-fi, Tue Aug 21 07:11:33 2007 UTC revision 1.11 by mocchiut, Tue Aug 11 12:58:25 2009 UTC
# Line 9  Line 9 
9    
10        c2f_error = 0        c2f_error = 0
11    
12    
13  c      if(C2F_DEBUG.EQ.1.eq.1)print *                    c      if(C2F_DEBUG.EQ.1.eq.1)print *                  
14        if(verbose.eq.1)print *                          if(verbose.eq.1)print *                  
15       $     ,c2f_path(1:c2f_pathlen)//'trk-maskvk.dat'       $     ,c2f_path(1:c2f_pathlen)//'trk-maskvk.dat'
# Line 36  c            if(C2F_DEBUG.EQ.1.eq.1) Line 37  c            if(C2F_DEBUG.EQ.1.eq.1)
37              c2f_error = 1              c2f_error = 1
38              goto 1000              goto 1000
39           endif           endif
40             if(verbose.eq.1)print*,'view ',iv,': ',(mask_vk(iv,i),i=1,24)
41        enddo        enddo
42   1000 close(10)   1000 close(10)
43   1001 continue      c 1001 continue    
44        return        return
45        end        end
46    
# Line 86  c            if(C2F_DEBUG.EQ.1.eq.1)prin Line 88  c            if(C2F_DEBUG.EQ.1.eq.1)prin
88   1000    close(10)   1000    close(10)
89        enddo        enddo
90                
91            return         return
92        end        end
93  *** * * * *** * * * *** * * * *** * * * *** * * * *** * * * ***  *** * * * *** * * * *** * * * *** * * * *** * * * *** * * * ***
94        subroutine readchargeparam        subroutine readchargeparam
# Line 96  c            if(C2F_DEBUG.EQ.1.eq.1)prin Line 98  c            if(C2F_DEBUG.EQ.1.eq.1)prin
98        include 'calib.f'        include 'calib.f'
99        include 'common_c2f.f' !<<<< C2F        include 'common_c2f.f' !<<<< C2F
100        character*60 fname_param        character*60 fname_param
101            c2f_error = 0         c2f_error = 0
102   201  format('charge-l',i1,'.dat')   201  format('charge-l',i1,'.dat')
103        do ilad=1,nladders_view                do ilad=1,nladders_view        
104           write(fname_param,201)ilad           write(fname_param,201)ilad
# Line 120  c         if(C2F_DEBUG.EQ.1.eq.1)print * Line 122  c         if(C2F_DEBUG.EQ.1.eq.1)print *
122       $           )pip,       $           )pip,
123       $            kch(ip,ilad),cch(ip,ilad),sch(ip,ilad)               $            kch(ip,ilad),cch(ip,ilad),sch(ip,ilad)        
124              if(iostat.ne.0)then              if(iostat.ne.0)then
125                             c2f_error = 1                    c2f_error = 1
126                 goto 1000                 goto 1000
127              endif              endif
128           enddo           enddo
# Line 148  c         if(C2F_DEBUG.EQ.1.eq.1)print * Line 150  c         if(C2F_DEBUG.EQ.1.eq.1)print *
150                
151                
152        c2f_error = 0        c2f_error = 0
153                    
154  *     ======================================================  *     ======================================================
155  *     read ANGULAR BINNING info  *     read ANGULAR BINNING info
156  *     ======================================================  *     ======================================================
# Line 171  c         if(C2F_DEBUG.EQ.1.eq.1) Line 173  c         if(C2F_DEBUG.EQ.1.eq.1)
173           c2f_error = 1           c2f_error = 1
174           return           return
175        endif        endif
176  c      print*,'---- ANGULAR BINNING ----'        if(DEBUG.EQ.1)THEN
177  c      print*,'Bin   -   angL   -   angR'           print*,'---- ANGULAR BINNING ----'
178             print*,'Bin   -   angL   -   angR'
179          endif
180   101  format(i2,'       ',f6.2,'     ',f6.2)   101  format(i2,'       ',f6.2,'     ',f6.2)
181        xnn=0        xnn=0
182        do ibin=1,nangmax        do ibin=1,nangmax
# Line 180  c      print*,'Bin   -   angL   -   angR Line 184  c      print*,'Bin   -   angL   -   angR
184       $        ,IOSTAT=iostat       $        ,IOSTAT=iostat
185       $        )xnn,angL(ibin),angR(ibin)       $        )xnn,angL(ibin),angR(ibin)
186           if(iostat.ne.0)goto 1000           if(iostat.ne.0)goto 1000
187  c         write(*,101)int(xnn),angL(ibin),angR(ibin)           if(DEBUG.EQ.1)write(*,101)int(xnn),angL(ibin),angR(ibin)
188        enddo                enddo        
189   1000 nangbin=int(xnn)   1000 nangbin=int(xnn)
190        close(10)        close(10)
191        if(nangbin.eq.0)c2f_error = 1        if(nangbin.eq.0)c2f_error = 1
192  c      print*,'-------------------------'        if(DEBUG.EQ.1)THEN
193  c      print*,nangbin           print*,'-------------------------'
194             print*,nangbin
195          endif
196                
197    
198    
# Line 332  c     $           ii,aal,aar,(fcorr(iv,i Line 338  c     $           ii,aal,aar,(fcorr(iv,i
338  c      character*120 cmd1  c      character*120 cmd1
339  c      character*120 cmd2  c      character*120 cmd2
340    
341            c2f_error = 0         c2f_error = 0
342        call mech_sensor        call mech_sensor
343        do ip=1,nplanes        do ip=1,nplanes
344          fitz(ip)=z_mech_sensor(ip,1,1)*0.1 !cm          fitz(ip)=z_mech_sensor(ip,1,1)*0.1 !cm
# Line 359  c            if(C2F_DEBUG.EQ.1.eq.1)prin Line 365  c            if(C2F_DEBUG.EQ.1.eq.1)prin
365       $           ,STATUS='OLD',IOSTAT=iostat       $           ,STATUS='OLD',IOSTAT=iostat
366       $           )       $           )
367              if(iostat.ne.0)then              if(iostat.ne.0)then
368                           c2f_error = 1                     c2f_error = 1
369                   return                   return
370              endif              endif
371                            
# Line 422  c     ! last plane inversion Line 428  c     ! last plane inversion
428    
429        real xvec(nladders_view),yvec(2),zvec(nplanes)        real xvec(nladders_view),yvec(2),zvec(nplanes)
430    
431            c2f_error = 0         c2f_error = 0
432                    
433  c      if(C2F_DEBUG.EQ.1.eq.1)print *  c      if(C2F_DEBUG.EQ.1.eq.1)print *
434        if(VERBOSE.EQ.1)print *        if(VERBOSE.EQ.1)print *
435       $     ,c2f_path(1:c2f_pathlen)//'mech_pos.dat'       $     ,c2f_path(1:c2f_pathlen)//'mech_pos.dat'
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.23