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

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

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

revision 1.20 by pam-fi, Fri Apr 27 10:39:58 2007 UTC revision 1.21 by mocchiut, Fri Apr 27 12:11:02 2007 UTC
# Line 331  c$$$         enddo Line 331  c$$$         enddo
331              iimage=0              iimage=0
332           endif           endif
333           if(icand.eq.0)then           if(icand.eq.0)then
334              print*,'HAI FATTO UN CASINO!!!!!! icand = ',icand              if(VERBOSE)then
335       $           ,ibest,iimage                 print*,'HAI FATTO UN CASINO!!!!!! icand = ',icand
336         $              ,ibest,iimage
337                endif
338              return              return
339           endif           endif
340    
# Line 690  c$$$         print*,fbad_cog(4,icx) Line 692  c$$$         print*,fbad_cog(4,icx)
692              resxPAM = resxPAM*fbad_cog(0,icx)              resxPAM = resxPAM*fbad_cog(0,icx)
693    
694           else           else
695              print*,'*** Non valid p.f.a. (x) --> ',PFAx              if(DEBUG) print*,'*** Non valid p.f.a. (x) --> ',PFAx
696           endif           endif
697    
698    
# Line 720  c$$$            print*,icx,' *** ',resxP Line 722  c$$$            print*,icx,' *** ',resxP
722           stripy = float(MAXS(icy))           stripy = float(MAXS(icy))
723    
724           if(icx.ne.0.and.(nply.ne.nplx.or.nldy.ne.nldx))then           if(icx.ne.0.and.(nply.ne.nplx.or.nldy.ne.nldx))then
725              print*,'xyz_PAM   ***ERROR*** invalid cluster couple!!! '              if(DEBUG) then
726       $           ,icx,icy                 print*,'xyz_PAM   ***ERROR*** invalid cluster couple!!! '
727         $              ,icx,icy
728                endif
729              goto 100              goto 100
730           endif           endif
731  *        --------------------------  *        --------------------------
# Line 802  c$$$         print*,fbad_cog(4,icy) Line 806  c$$$         print*,fbad_cog(4,icy)
806              resyPAM = resyPAM*fbad_cog(0,icy)              resyPAM = resyPAM*fbad_cog(0,icy)
807    
808           else           else
809              print*,'*** Non valid p.f.a. (x) --> ',PFAx              if(DEBUG) print*,'*** Non valid p.f.a. (x) --> ',PFAx
810           endif           endif
811    
812    
# Line 832  c     (xi,yi,zi) = mechanical coordinate Line 836  c     (xi,yi,zi) = mechanical coordinate
836  c------------------------------------------------------------------------  c------------------------------------------------------------------------
837           if(((mod(int(stripx+0.5)-1,1024)+1).le.3)           if(((mod(int(stripx+0.5)-1,1024)+1).le.3)
838       $        .or.((mod(int(stripx+0.5)-1,1024)+1).ge.1022)) then !X has 1018 strips...       $        .or.((mod(int(stripx+0.5)-1,1024)+1).ge.1022)) then !X has 1018 strips...
839              print*,'xyz_PAM (couple):',              if(DEBUG) then
840       $          ' WARNING: false X strip: strip ',stripx                 print*,'xyz_PAM (couple):',
841         $              ' WARNING: false X strip: strip ',stripx
842                endif
843           endif           endif
844           xi = acoordsi(stripx,viewx)           xi = acoordsi(stripx,viewx)
845           yi = acoordsi(stripy,viewy)           yi = acoordsi(stripy,viewy)
# Line 925  c            print*,'X-singlet ',icx,npl Line 931  c            print*,'X-singlet ',icx,npl
931  c            if((stripx.le.3).or.(stripx.ge.1022)) then !X has 1018 strips...  c            if((stripx.le.3).or.(stripx.ge.1022)) then !X has 1018 strips...
932              if(((mod(int(stripx+0.5)-1,1024)+1).le.3)              if(((mod(int(stripx+0.5)-1,1024)+1).le.3)
933       $           .or.((mod(int(stripx+0.5)-1,1024)+1).ge.1022)) then !X has 1018 strips...       $           .or.((mod(int(stripx+0.5)-1,1024)+1).ge.1022)) then !X has 1018 strips...
934                 print*,'xyz_PAM (X-singlet):',                 if(DEBUG) then
935       $             ' WARNING: false X strip: strip ',stripx                    print*,'xyz_PAM (X-singlet):',
936         $                 ' WARNING: false X strip: strip ',stripx
937                   endif
938              endif              endif
939              xi   = acoordsi(stripx,viewx)              xi   = acoordsi(stripx,viewx)
940    
# Line 948  c            print*,'X-cl ',icx,stripx,' Line 956  c            print*,'X-cl ',icx,stripx,'
956  c            print*,yi_A,' <--> ',yi_B  c            print*,yi_A,' <--> ',yi_B
957    
958           else           else
959                if(DEBUG) then
960              print *,'routine xyz_PAM ---> not properly used !!!'                 print *,'routine xyz_PAM ---> not properly used !!!'
961              print *,'icx = ',icx                 print *,'icx = ',icx
962              print *,'icy = ',icy                 print *,'icy = ',icy
963                endif
964              goto 100              goto 100
965                            
966           endif           endif
# Line 1016  c--------------------------------------- Line 1025  c---------------------------------------
1025  c         print*,'A-(',xPAM_A,yPAM_A,') B-(',xPAM_B,yPAM_B,')'  c         print*,'A-(',xPAM_A,yPAM_A,') B-(',xPAM_B,yPAM_B,')'
1026    
1027        else        else
1028                       if(DEBUG) then
1029           print *,'routine xyz_PAM ---> not properly used !!!'              print *,'routine xyz_PAM ---> not properly used !!!'
1030           print *,'icx = ',icx              print *,'icx = ',icx
1031           print *,'icy = ',icy              print *,'icy = ',icy
1032                         endif
1033        endif        endif
1034                    
1035    
# Line 1146  c$$$         print*,' resolution ',resxP Line 1155  c$$$         print*,' resolution ',resxP
1155                    
1156        else        else
1157                    
1158           print*  c         print*
1159       $        ,' function distance_to ---> wrong usage!!!'  c     $        ,' function distance_to ---> wrong usage!!!'
1160           print*,' xPAM,yPAM,zPAM ',xPAM,yPAM,zPAM  c         print*,' xPAM,yPAM,zPAM ',xPAM,yPAM,zPAM
1161           print*,' xPAM_A,yPAM_A,zPAM_A,xPAM_b,yPAM_b,zPAM_b '  c         print*,' xPAM_A,yPAM_A,zPAM_A,xPAM_b,yPAM_b,zPAM_b '
1162       $        ,xPAM_A,yPAM_A,zPAM_A,xPAM_b,yPAM_b,zPAM_b  c     $        ,xPAM_A,yPAM_A,zPAM_A,xPAM_b,yPAM_b,zPAM_b
1163        endif          endif  
1164    
1165        distance_to = sngl(distance)        distance_to = sngl(distance)
# Line 1218  c--------------------------------------- Line 1227  c---------------------------------------
1227                 if(((mod(int(stripx+0.5)-1,1024)+1).le.3)                 if(((mod(int(stripx+0.5)-1,1024)+1).le.3)
1228       $              .or.((mod(int(stripx+0.5)-1,1024)+1).ge.1022)) then !X has 1018 strips...       $              .or.((mod(int(stripx+0.5)-1,1024)+1).ge.1022)) then !X has 1018 strips...
1229  c     if((stripx.le.3).or.(stripx.ge.1022)) then !X has 1018 strips...  c     if((stripx.le.3).or.(stripx.ge.1022)) then !X has 1018 strips...
1230                    print*,'whichsensor: ',  c                  print*,'whichsensor: ',
1231       $                ' WARNING: false X strip: strip ',stripx  c     $                ' WARNING: false X strip: strip ',stripx
1232                 endif                 endif
1233                 xi = acoordsi(stripx,viewx)                 xi = acoordsi(stripx,viewx)
1234                 yi = acoordsi(stripy,viewy)                 yi = acoordsi(stripy,viewy)
# Line 1374  c      include 'common_analysis.f' Line 1383  c      include 'common_analysis.f'
1383        is_cp=0        is_cp=0
1384        if(id.lt.0)is_cp=1        if(id.lt.0)is_cp=1
1385        if(id.gt.0)is_cp=2        if(id.gt.0)is_cp=2
1386        if(id.eq.0)print*,'IS_CP ===> wrong couple id !!!'  c      if(id.eq.0)print*,'IS_CP ===> wrong couple id !!!'
1387    
1388        return        return
1389        end        end
Legend:
Removed from v.1.20  
changed lines
  Added in v.1.21
  ViewVC Help
Powered by ViewVC 1.1.23