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

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

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

revision 1.2 by pam-fi, Tue Nov 21 17:12:59 2006 UTC revision 1.6 by mocchiut, Fri Jan 17 12:56:52 2014 UTC
# Line 51  c     integer nloops            !number Line 51  c     integer nloops            !number
51        integer ifail             !=-1 if singular matrix error, =0 if not singular        integer ifail             !=-1 if singular matrix error, =0 if not singular
52        integer jfail             !=0 if determinant can be evaluated, =-1 if determinat is probably too small, =+1 if too large        integer jfail             !=0 if determinant can be evaluated, =-1 if determinat is probably too small, =+1 if too large
53    
54        parameter (big=1.e4)      !just a number greater than C(npoints,3)        parameter (ibig=10000)      !just a number greater than C(npoints,3) !EM GCC4.7
55        double precision xxc(big),zzc(big),rrr(big) !centres and radii to be averaged        double precision xxc(ibig),zzc(ibig),rrr(ibig) !centres and radii to be averaged
56    
57        double precision tmp1,tmp2,tmp(npoints)    !temp variables        double precision tmp1,tmp2,tmp(npoints)    !temp variables
58    
# Line 174  c--------------------------------------- Line 174  c---------------------------------------
174              call DFACT(3,d,3,ir,ifail,detd,jfail)              call DFACT(3,d,3,ir,ifail,detd,jfail)
175              if(ifail.eq.-1) then              if(ifail.eq.-1) then
176                if(DEBUG)then                if(DEBUG)then
177                print*,'tricircle: ERROR: singular matrix D:'                   print*,'tricircle: ERROR: singular matrix D:'
178                do i=1,3                   do i=1,3
179                  print*,(d(i,j),j=1,3)                      print*,(d(i,j),j=1,3)
180                enddo                   enddo
181                endif                endif
182                eflag=1                eflag=1
183              endif             endif
184              if(jfail.eq.-1) then             if(jfail.eq.-1) then
185                if(DEBUG)then                if(DEBUG)then
186                print*                   print*
187       $             ,'tricircle: ERROR: matrix D: determinant too small?'       $      ,'tricircle: ERROR: matrix D: determinant too small?'
188                do i=1,3                   do i=1,3
189                  print*,(d(i,j),j=1,3)                      print*,(d(i,j),j=1,3)
190                enddo                   enddo
191                endif                endif
192                eflag=1                eflag=1
193              elseif(jfail.eq.1) then              elseif(jfail.eq.1) then
# Line 212  c--------------------------------------- Line 212  c---------------------------------------
212                enddo                enddo
213                endif                endif
214                eflag=1                eflag=1
215              endif             endif
216              if(jfail.eq.-1) then              if(jfail.eq.-1) then
217                if(DEBUG)then                if(DEBUG)then
218                print*                print*
# Line 296  c--------------------------------------- Line 296  c---------------------------------------
296        zc=0.        zc=0.
297        radius=0.        radius=0.
298        do i=1,k        do i=1,k
299          xc=xc+xxc(i)          xc=xc+REAL(xxc(i)) !EM GCC4.7
300          zc=zc+zzc(i)          zc=zc+REAL(zzc(i)) !EM GCC4.7
301          radius=radius+rrr(i)          radius=radius+REAL(rrr(i)) !EM GCC4.7
302        enddo        enddo
303        xc=xc/k * scale           !back to micrometers        xc=xc/k * scale           !back to micrometers
304        zc=zc/k * scale        zc=zc/k * scale
# Line 335  c     print*,xc,zc,radius !??? Line 335  c     print*,xc,zc,radius !???
335          else          else
336            tmp(i)=-tmp1            tmp(i)=-tmp1
337          endif          endif
338          residual(i)=tmp2 - tmp(i)          residual(i)=REAL(tmp2 - tmp(i)) !EM GCC4.7
339          chi=chi + residual(i)**2.          chi=chi + residual(i)**2.
340  c     print*,dep(i)                 !???  c     print*,dep(i)                 !???
341  c     print*,indep(i)              !???  c     print*,indep(i)              !???
# Line 347  c     it computes the angle between the Line 347  c     it computes the angle between the
347  c     independent variable axis  c     independent variable axis
348  c------------------------------------------------------------------------  c------------------------------------------------------------------------
349        do i=1,npoints        do i=1,npoints
350          angle(i)=(zc-indep(i)) / tmp(i)          angle(i)=REAL((zc-indep(i)) / tmp(i)) !EM GCC4.7
351          angle(i)=ATAN(angle(i)) !-pi/2 <= angle <= pi/2          angle(i)=ATAN(angle(i)) !-pi/2 <= angle <= pi/2
352          angle(i)=angle(i)/pigr*180.          angle(i)=angle(i)/pigr*180.
353        enddo        enddo
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.23