/[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.4 by pam-fi, Thu Sep 28 14:04:40 2006 UTC revision 1.5 by pam-fi, Fri Sep 29 08:13:04 2006 UTC
# Line 311  c         print*,'++++++++++ iimage,fima Line 311  c         print*,'++++++++++ iimage,fima
311  c     $        ,iimage,fimage,ntrk,image(ntrk)  c     $        ,iimage,fimage,ntrk,image(ntrk)
312    
313           if(ntrk.eq.NTRKMAX)then           if(ntrk.eq.NTRKMAX)then
314              if(DEBUG)              if(verbose)
315       $           print*,       $           print*,
316       $           '** warning ** number of identified '//       $           '** warning ** number of identified '//
317       $           'tracks exceeds vector dimension '       $           'tracks exceeds vector dimension '
# Line 1807  c$$$               endif Line 1807  c$$$               endif
1807                                
1808  *     ------------------> COUPLE <------------------  *     ------------------> COUPLE <------------------
1809  *     check to do not overflow vector dimentions  *     check to do not overflow vector dimentions
1810                 if(ncp_plane(nplx).gt.ncouplemax)then  c$$$               if(ncp_plane(nplx).gt.ncouplemax)then
                   if(DEBUG)print*,  
      $                    ' ** warning ** number of identified'//  
      $                    ' couples on plane ',nplx,  
      $                    ' exceeds vector dimention'//  
      $                    ' ( ',ncouplemax,' )'  
 c     good2=.false.  
 c     goto 880   !fill ntp and go to next event  
                   iflag=1  
                   return  
                endif  
                 
 c$$$               if(ncp_plane(nplx).eq.ncouplemax)then  
1811  c$$$                  if(DEBUG)print*,  c$$$                  if(DEBUG)print*,
1812  c$$$     $                 '** warning ** number of identified '//  c$$$     $                    ' ** warning ** number of identified'//
1813  c$$$     $                 'couples on plane ',nplx,  c$$$     $                    ' couples on plane ',nplx,
1814  c$$$     $                 'exceeds vector dimention '  c$$$     $                    ' exceeds vector dimention'//
1815  c$$$     $                 ,'( ',ncouplemax,' )'  c$$$     $                    ' ( ',ncouplemax,' )'
1816  c$$$c     good2=.false.  c$$$c     good2=.false.
1817  c$$$c     goto 880   !fill ntp and go to next event                      c$$$c     goto 880   !fill ntp and go to next event
1818  c$$$                  iflag=1  c$$$                  iflag=1
1819  c$$$                  return  c$$$                  return
1820  c$$$               endif  c$$$               endif
1821                                
1822                   if(ncp_plane(nplx).eq.ncouplemax)then
1823                      if(verbose)print*,
1824         $                 '** warning ** number of identified '//
1825         $                 'couples on plane ',nplx,
1826         $                 'exceeds vector dimention '
1827         $                 ,'( ',ncouplemax,' )'
1828    c     good2=.false.
1829    c     goto 880   !fill ntp and go to next event                    
1830                      iflag=1
1831                      return
1832                   endif
1833                  
1834                 ncp_plane(nplx) = ncp_plane(nplx) + 1                 ncp_plane(nplx) = ncp_plane(nplx) + 1
1835                 clx(nplx,ncp_plane(nplx))=icx                 clx(nplx,ncp_plane(nplx))=icx
1836                 cly(nply,ncp_plane(nplx))=icy                 cly(nply,ncp_plane(nplx))=icy
# Line 1868  c$$$               endif Line 1868  c$$$               endif
1868  c     if(ncp_tot.gt.ncp_max)goto 100!next event (TEMPORANEO!!!)  c     if(ncp_tot.gt.ncp_max)goto 100!next event (TEMPORANEO!!!)
1869                
1870        if(ncp_tot.gt.ncp_max)then        if(ncp_tot.gt.ncp_max)then
1871           if(DEBUG)print*,           if(verbose)print*,
1872       $           '** warning ** number of identified '//       $           '** warning ** number of identified '//
1873       $           'couples exceeds upper limit for Hough tr. '       $           'couples exceeds upper limit for Hough tr. '
1874       $           ,'( ',ncp_max,' )'                   $           ,'( ',ncp_max,' )'            
# Line 2027  c$$$               if(abs(ddd).gt.cut)go Line 2027  c$$$               if(abs(ddd).gt.cut)go
2027                                
2028  *     ------------------> COUPLE <------------------  *     ------------------> COUPLE <------------------
2029  *     check to do not overflow vector dimentions  *     check to do not overflow vector dimentions
2030                 if(ncp_plane(nplx).gt.ncouplemax)then  c$$$               if(ncp_plane(nplx).gt.ncouplemax)then
                   if(DEBUG)print*,  
      $                    ' ** warning ** number of identified'//  
      $                    ' couples on plane ',nplx,  
      $                    ' exceeds vector dimention'//  
      $                    ' ( ',ncouplemax,' )'  
 c     good2=.false.  
 c     goto 880   !fill ntp and go to next event  
                   iflag=1  
                   return  
                endif  
                 
 c$$$               if(ncp_plane(nplx).eq.ncouplemax)then  
2031  c$$$                  if(DEBUG)print*,  c$$$                  if(DEBUG)print*,
2032  c$$$     $                 '** warning ** number of identified '//  c$$$     $                    ' ** warning ** number of identified'//
2033  c$$$     $                 'couples on plane ',nplx,  c$$$     $                    ' couples on plane ',nplx,
2034  c$$$     $                 'exceeds vector dimention '  c$$$     $                    ' exceeds vector dimention'//
2035  c$$$     $                 ,'( ',ncouplemax,' )'  c$$$     $                    ' ( ',ncouplemax,' )'
2036  c$$$c     good2=.false.  c$$$c     good2=.false.
2037  c$$$c     goto 880   !fill ntp and go to next event                      c$$$c     goto 880   !fill ntp and go to next event
2038  c$$$                  iflag=1  c$$$                  iflag=1
2039  c$$$                  return  c$$$                  return
2040  c$$$               endif  c$$$               endif
2041                                
2042                   if(ncp_plane(nplx).eq.ncouplemax)then
2043                      if(verbose)print*,
2044         $                 '** warning ** number of identified '//
2045         $                 'couples on plane ',nplx,
2046         $                 'exceeds vector dimention '
2047         $                 ,'( ',ncouplemax,' )'
2048    c     good2=.false.
2049    c     goto 880   !fill ntp and go to next event                    
2050                      iflag=1
2051                      return
2052                   endif
2053                  
2054                 ncp_plane(nplx) = ncp_plane(nplx) + 1                 ncp_plane(nplx) = ncp_plane(nplx) + 1
2055                 clx(nplx,ncp_plane(nplx))=icx                 clx(nplx,ncp_plane(nplx))=icx
2056                 cly(nply,ncp_plane(nplx))=icy                 cly(nply,ncp_plane(nplx))=icy
# Line 2088  c$$$               endif Line 2088  c$$$               endif
2088  c     if(ncp_tot.gt.ncp_max)goto 100!next event (TEMPORANEO!!!)  c     if(ncp_tot.gt.ncp_max)goto 100!next event (TEMPORANEO!!!)
2089                
2090        if(ncp_tot.gt.ncp_max)then        if(ncp_tot.gt.ncp_max)then
2091           if(DEBUG)print*,           if(verbose)print*,
2092       $           '** warning ** number of identified '//       $           '** warning ** number of identified '//
2093       $           'couples exceeds upper limit for Hough tr. '       $           'couples exceeds upper limit for Hough tr. '
2094       $           ,'( ',ncp_max,' )'                   $           ,'( ',ncp_max,' )'            
# Line 2402  c     $                       (icx2,icy2 Line 2402  c     $                       (icx2,icy2
2402  *     (2 couples needed)  *     (2 couples needed)
2403  *     - - - - - - - - - - - - - - - - - - - - - - - - - - - -  *     - - - - - - - - - - - - - - - - - - - - - - - - - - - -
2404                          if(ndblt.eq.ndblt_max)then                          if(ndblt.eq.ndblt_max)then
2405                             if(DEBUG)print*,                             if(verbose)print*,
2406       $                          '** warning ** number of identified '//       $                          '** warning ** number of identified '//
2407       $                          'doublets exceeds vector dimention '       $                          'doublets exceeds vector dimention '
2408       $                          ,'( ',ndblt_max,' )'       $                          ,'( ',ndblt_max,' )'
# Line 2472  c     $                                 Line 2472  c     $                                
2472  *     (3 couples needed)  *     (3 couples needed)
2473  *     - - - - - - - - - - - - - - - - - - - - - - - - - - - -  *     - - - - - - - - - - - - - - - - - - - - - - - - - - - -
2474                                   if(ntrpt.eq.ntrpt_max)then                                   if(ntrpt.eq.ntrpt_max)then
2475                                      if(DEBUG)print*,                                      if(verbose)print*,
2476       $                     '** warning ** number of identified '//       $                     '** warning ** number of identified '//
2477       $                     'triplets exceeds vector dimention '       $                     'triplets exceeds vector dimention '
2478       $                    ,'( ',ntrpt_max,' )'       $                    ,'( ',ntrpt_max,' )'
# Line 2699  c     print*,'>>>> ',ncpused,npt,nplused Line 2699  c     print*,'>>>> ',ncpused,npt,nplused
2699  *     >>> NEW CLOUD <<<  *     >>> NEW CLOUD <<<
2700    
2701           if(nclouds_yz.ge.ncloyz_max)then           if(nclouds_yz.ge.ncloyz_max)then
2702              if(DEBUG)print*,              if(verbose)print*,
2703       $           '** warning ** number of identified '//       $           '** warning ** number of identified '//
2704       $           'YZ clouds exceeds vector dimention '       $           'YZ clouds exceeds vector dimention '
2705       $           ,'( ',ncloyz_max,' )'       $           ,'( ',ncloyz_max,' )'
# Line 2911  c     print*,'check cp_used' Line 2911  c     print*,'check cp_used'
2911  *     ~~~~~~~~~~~~~~~~~  *     ~~~~~~~~~~~~~~~~~
2912  *     >>> NEW CLOUD <<<  *     >>> NEW CLOUD <<<
2913           if(nclouds_xz.ge.ncloxz_max)then           if(nclouds_xz.ge.ncloxz_max)then
2914              if(DEBUG)print*,              if(verbose)print*,
2915       $           '** warning ** number of identified '//       $           '** warning ** number of identified '//
2916       $           'XZ clouds exceeds vector dimention '       $           'XZ clouds exceeds vector dimention '
2917       $           ,'( ',ncloxz_max,' )'       $           ,'( ',ncloxz_max,' )'
# Line 3211  c     $                                 Line 3211  c     $                                
3211  *     --------------------------  *     --------------------------
3212                                if(ntracks.eq.NTRACKSMAX)then                                if(ntracks.eq.NTRACKSMAX)then
3213                                                                    
3214                                   if(DEBUG)print*,                                   if(verbose)print*,
3215       $                 '** warning ** number of candidate tracks '//       $                 '** warning ** number of candidate tracks '//
3216       $                 ' exceeds vector dimension '       $                 ' exceeds vector dimension '
3217       $                ,'( ',NTRACKSMAX,' )'       $                ,'( ',NTRACKSMAX,' )'
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.23