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

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

Parent Directory | Revision Log | View Patch Patch

-revision 1.18 by pam-fi,
Fri Feb 16 14:56:02 2007 UTC
+revision 1.31 by pam-fi,
Fri Aug 31 14:56:52 2007 UTC
 Line 20
        include 'calib.f'
        include 'level2.f'
- c      include 'momanhough_init.f'
+ c      print*,'======================================================'
+ c$$$      do ic=1,NCLSTR1
+ c$$$         if(.false.
+ c$$$     $        .or.nsatstrips(ic).gt.0
+ c$$$c     $        .or.nbadstrips(0,ic).gt.0
+ c$$$c     $        .or.nbadstrips(4,ic).gt.0
+ c$$$c     $        .or.nbadstrips(3,ic).gt.0
+ c$$$     $        .or..false.)then
+ c$$$            print*,'--- cl-',ic,' ------------------------'
+ c$$$            istart = INDSTART(IC)
+ c$$$            istop  = TOTCLLENGTH
+ c$$$            if(ic.lt.NCLSTR1)istop=INDSTART(IC+1)-1
+ c$$$            print*,'ADC   ',(CLADC(i),i=istart,istop)
+ c$$$            print*,'s/n   ',(CLSIGNAL(i)/CLSIGMA(i),i=istart,istop)
+ c$$$            print*,'sgnl  ',(CLSIGNAL(i),i=istart,istop)
+ c$$$            print*,'strip ',(i-INDMAX(ic),i=istart,istop)
+ c$$$            print*,'view ',VIEW(ic)
+ c$$$            print*,'maxs ',MAXS(ic)
+ c$$$            print*,'COG4 ',cog(4,ic)
+ c$$$            ff = fbad_cog(4,ic)
+ c$$$            print*,'fbad ',ff
+ c$$$            print*,(CLBAD(i),i=istart,istop)
+ c$$$            bb=nbadstrips(0,ic)
+ c$$$            print*,'#BAD (tot)',bb
+ c$$$            bb=nbadstrips(4,ic)
+ c$$$            print*,'#BAD (4)',bb
+ c$$$            bb=nbadstrips(3,ic)
+ c$$$            print*,'#BAD (3)',bb
+ c$$$            ss=nsatstrips(ic)
+ c$$$            print*,'#saturated ',ss
+ c$$$         endif
+ c$$$      enddo
  *-------------------------------------------------------------------------------
  *     STEP 1
-Line 41 
 c      include 'momanhough_init.f'
+Line 74 
 c      include 'momanhough_init.f'
  *-------------------------------------------------------------------------------
  *-------------------------------------------------------------------------------
- c      iflag=0
        call cl_to_couples(iflag)
        if(iflag.eq.1)then        !bad event
           goto 880               !go to next event
-Line 75 
 c      iflag=0
+Line 108 
 c      iflag=0
  *-------------------------------------------------------------------------------
  *-------------------------------------------------------------------------------
- c      iflag=0
        call cp_to_doubtrip(iflag)
        if(iflag.eq.1)then        !bad event
           goto 880               !go to next event
-Line 182 
 c      iflag=0
+Line 215 
 c      iflag=0
  c      include 'momanhough_init.f'
        logical FIMAGE            !
+       real trackimage(NTRACKSMAX)
        real*8 AL_GUESS(5)
  *-------------------------------------------------------------------------------
-Line 239 
 c$$$         if(ibest.eq.0)goto 880 !>>
+Line 273 
 c$$$         if(ibest.eq.0)goto 880 !>>
  *     2nd) increasing chi**2
  *     -------------------------------------------------------
           rchi2best=1000000000.
-          ndofbest=0             !(1)
+          ndofbest=0
           do i=1,ntracks
-            ndof=0               !(1)
+            ndof=0
-            do ii=1,nplanes      !(1)
+            do ii=1,nplanes
-              ndof=ndof          !(1)
+              ndof=ndof
-      $            +int(xgood_store(ii,i)) !(1)
+      $            +int(xgood_store(ii,i))
-      $            +int(ygood_store(ii,i)) !(1)
+      $            +int(ygood_store(ii,i))
-            enddo                !(1)
+            enddo
-            if(ndof.gt.ndofbest)then !(1)
+            if(ndof.gt.ndofbest)then
               ibest=i
               rchi2best=RCHI2_STORE(i)
-              ndofbest=ndof      !(1)
+              ndofbest=ndof
-            elseif(ndof.eq.ndofbest)then !(1)
+            elseif(ndof.eq.ndofbest)then
               if(RCHI2_STORE(i).lt.rchi2best.and.
       $            RCHI2_STORE(i).gt.0)then
                 ibest=i
                 rchi2best=RCHI2_STORE(i)
-                ndofbest=ndof    !(1)
+                ndofbest=ndof
-              endif              !(1)
+              endif
             endif
           enddo
-Line 298 
 c$$$         enddo
+Line 332 
 c$$$         enddo
              iimage=0
           endif
           if(icand.eq.0)then
-             print*,'HAI FATTO UN CASINO!!!!!! icand = ',icand
+             if(VERBOSE.EQ.1)then
-      $           ,ibest,iimage
+                print*,'HAI FATTO UN CASINO!!!!!! icand = ',icand
+      $              ,ibest,iimage
+             endif
              return
           endif
-Line 325 
 c         print*,'## guess: ',al
+Line 361 
 c         print*,'## guess: ',al
           jstep=0                !# minimization steps
           iprint=0
- c         if(DEBUG)iprint=1
+ c         if(DEBUG.EQ.1)iprint=1
-          if(VERBOSE)iprint=1
+          if(VERBOSE.EQ.1)iprint=1
-          if(DEBUG)iprint=2
+          if(DEBUG.EQ.1)iprint=2
           call mini2(jstep,ifail,iprint)
           if(ifail.ne.0) then
-             if(VERBOSE)then
+             if(VERBOSE.EQ.1)then
                 print *,
       $              '*** MINIMIZATION FAILURE *** (after refinement) '
       $              ,iev
-Line 345 
 c$$$               print*,'-------------
+Line 381 
 c$$$               print*,'-------------
  c            chi2=-chi2
           endif
-          if(DEBUG)then
+          if(DEBUG.EQ.1)then
              print*,'----------------------------- improved track coord'
       format(i2,' * ',3f10.4,' --- ',4f10.4,' --- ',2f4.0,2f10.5)
              do ip=1,6
-Line 356 
 c            chi2=-chi2
+Line 392 
 c            chi2=-chi2
           endif
  c         rchi2=chi2/dble(ndof)
-          if(DEBUG)then
+          if(DEBUG.EQ.1)then
              print*,' '
              print*,'****** SELECTED TRACK *************'
              print*,'#         R. chi2        RIG'
-Line 372 
 c         rchi2=chi2/dble(ndof)
+Line 408 
 c         rchi2=chi2/dble(ndof)
  *     ------------- search if the track has an IMAGE -------------
  *     ------------- (also this is stored )           -------------
           if(FIMAGE)goto 122     !>>> jump! (this is already an image)
- *     now search for track-image, by comparing couples IDs
+ *     -----------------------------------------------------
+ *     first check if the track is ambiguous
+ *     -----------------------------------------------------
+ *     (modified on august 2007 by ElenaV)
+          is1=0
+          do ip=1,NPLANES
+             if(SENSOR_STORE(ip,icand).ne.0)then
+                is1=SENSOR_STORE(ip,icand)
+                if(ip.eq.6)is1=3-is1 !last plane inverted
+             endif
+          enddo
+          if(is1.eq.0)then
+             if(WARNING.EQ.1)print*,'** WARNING ** sensor=0'
+             goto 122            !jump
+          endif
+ c         print*,'is1 ',is1
+          do ip=1,NPLANES
+             is2 = SENSOR_STORE(ip,icand) !sensor
+ c            print*,'is2 ',is2,' ip ',ip
+             if(ip.eq.6.and.is2.ne.0)is2=3-is2 !last plane inverted
+             if(
+      $           (is1.ne.is2.and.is2.ne.0)
+      $           )goto 122      !jump (this track cannot have an image)
+          enddo
+          if(DEBUG.eq.1)print*,' >>> ambiguous track! '
+ *     now search for track-image among track candidates
+ c$$$         do i=1,ntracks
+ c$$$            iimage=i
+ c$$$            do ip=1,nplanes
+ c$$$               if(     CP_STORE(nplanes-ip+1,icand).ne.
+ c$$$     $              -1*CP_STORE(nplanes-ip+1,i).and.
+ c$$$     $              CP_STORE(nplanes-ip+1,i).ne.0.and.
+ c$$$     $              CP_STORE(nplanes-ip+1,icand).ne.0 )iimage=0
+ c$$$               print*,' track ',i,' CP ',CP_STORE(nplanes-ip+1,i)
+ c$$$     $              ,CP_STORE(nplanes-ip+1,icand),iimage
+ c$$$            enddo
+ c$$$            if(  iimage.ne.0.and.
+ c$$$c     $           RCHI2_STORE(i).le.CHI2MAX.and.
+ c$$$c     $           RCHI2_STORE(i).gt.0.and.
+ c$$$     $           .true.)then
+ c$$$               if(DEBUG.EQ.1)print*,'Track candidate ',iimage
+ c$$$     $              ,' >>> TRACK IMAGE >>> of'
+ c$$$     $              ,ibest
+ c$$$               goto 122         !image track found
+ c$$$            endif
+ c$$$         enddo
+ *     ---------------------------------------------------------------
+ *     take the candidate with the greatest number of matching couples
+ *     if more than one satisfies the requirement,
+ *     choose the one with more points and lower chi2
+ *     ---------------------------------------------------------------
+ *     count the number of matching couples
+          ncpmax = 0
+          iimage   = 0           !id of candidate with better matching
           do i=1,ntracks
-             iimage=i
+             ncp=0
              do ip=1,nplanes
-                if(     CP_STORE(nplanes-ip+1,icand).ne.
+                if(CP_STORE(nplanes-ip+1,icand).ne.0)then
-      $              -1*CP_STORE(nplanes-ip+1,i) )iimage=0
+                   if(
+      $                 CP_STORE(nplanes-ip+1,i).ne.0
+      $                 .and.
+      $                 CP_STORE(nplanes-ip+1,icand).eq.
+      $                 -1*CP_STORE(nplanes-ip+1,i)
+      $                 )then
+                      ncp=ncp+1  !ok
+                   elseif(
+      $                    CP_STORE(nplanes-ip+1,i).ne.0
+      $                    .and.
+      $                    CP_STORE(nplanes-ip+1,icand).ne.
+      $                    -1*CP_STORE(nplanes-ip+1,i)
+      $                    )then
+                      ncp = 0
+                      goto 117   !it is not an image candidate
+                   else
+                   endif
+                endif
+ c$$$               print*,' track ',i,' CP ',CP_STORE(nplanes-ip+1,i)
+ c$$$     $              ,CP_STORE(nplanes-ip+1,icand),ncp
              enddo
-             if(  iimage.ne.0.and.
+       continue
- c     $           RCHI2_STORE(i).le.CHI2MAX.and.
+             trackimage(i)=ncp   !number of matching couples
- c     $           RCHI2_STORE(i).gt.0.and.
+             if(ncp>ncpmax)then
-      $           .true.)then
+                ncpmax=ncp
-                if(DEBUG)print*,'Track candidate ',iimage
+                iimage=i
-      $              ,' >>> TRACK IMAGE >>> of'
-      $              ,ibest
-                goto 122         !image track found
              endif
           enddo
+ *     check if there are more than one image candidates
+          ngood=0
+          do i=1,ntracks
+             if( ncpmax.ne.0.and.trackimage(i).eq.ncpmax )ngood=ngood+1
+          enddo
+ *     if there are, choose the best one
+          if(ngood.gt.1)then
+ *     -------------------------------------------------------
+ *     order track-candidates according to:
+ *     1st) decreasing n.points
+ *     2nd) increasing chi**2
+ *     -------------------------------------------------------
+             rchi2best=1000000000.
+             ndofbest=0
+             do i=1,ntracks
+                if( trackimage(i).eq.ncpmax )then
+                   ndof=0
+                   do ii=1,nplanes
+                      ndof=ndof
+      $                    +int(xgood_store(ii,i))
+      $                    +int(ygood_store(ii,i))
+                   enddo
+                   if(ndof.gt.ndofbest)then
+                      iimage=i
+                      rchi2best=RCHI2_STORE(i)
+                      ndofbest=ndof
+                   elseif(ndof.eq.ndofbest)then
+                      if(RCHI2_STORE(i).lt.rchi2best.and.
+      $                    RCHI2_STORE(i).gt.0)then
+                         iimage=i
+                         rchi2best=RCHI2_STORE(i)
+                         ndofbest=ndof
+                      endif
+                   endif
+                endif
+             enddo
+          endif
+          if(DEBUG.EQ.1)print*,'Track candidate ',iimage
+      $        ,' >>> TRACK IMAGE >>> of'
+      $        ,ibest
     continue
  *     --- and store the results --------------------------------
           ntrk = ntrk + 1                   !counter of found tracks
           if(.not.FIMAGE
-Line 403 
 c         print*,'++++++++++ iimage,fima
+Line 554 
 c         print*,'++++++++++ iimage,fima
  c     $        ,iimage,fimage,ntrk,image(ntrk)
           if(ntrk.eq.NTRKMAX)then
-             if(verbose)
+             if(verbose.eq.1)
       $           print*,
       $           '** warning ** number of identified '//
       $           'tracks exceeds vector dimension '
-Line 439 
 cc            good2=.false.
+Line 590 
 cc            good2=.false.
       $        rchi2best.le.CHI2MAX.and.
  c     $        rchi2best.lt.15..and.
       $        .true.)then
-             if(DEBUG)then
+             if(DEBUG.EQ.1)then
                 print*,'***** NEW SEARCH ****'
              endif
              goto 11111          !try new search
-Line 523 
 c     $        rchi2best.lt.15..and.
+Line 674 
 c     $        rchi2best.lt.15..and.
        include 'commontracker.f'
        include 'level1.f'
        include 'calib.f'
- c      include 'level1.f'
        include 'common_align.f'
        include 'common_mech.f'
        include 'common_xyzPAM.f'
- c      include 'common_resxy.f'
- c      logical DEBUG
- c      common/dbg/DEBUG
        integer icx,icy           !X-Y cluster ID
        integer sensor
-Line 545 
 c      common/dbg/DEBUG
+Line 691 
 c      common/dbg/DEBUG
        double precision xrt,yrt,zrt
        double precision xrt_A,yrt_A,zrt_A
        double precision xrt_B,yrt_B,zrt_B
- c      double precision xi,yi,zi
- c      double precision xi_A,yi_A,zi_A
- c      double precision xi_B,yi_B,zi_B
        parameter (ndivx=30)
+ c$$$      print*,icx,icy,sensor,PFAx,PFAy,ax,ay,bfx,bfy
        resxPAM = 0
        resyPAM = 0
-Line 565 
 c      double precision xi_B,yi_B,zi_B
+Line 711 
 c      double precision xi_B,yi_B,zi_B
        yPAM_B = 0.
        zPAM_B = 0.
  c      print*,'## xyz_PAM: ',icx,icy,sensor,PFAx,PFAy,angx,angy
+       if(sensor.lt.1.or.sensor.gt.2)then
+          print*,'xyz_PAM   ***ERROR*** wrong input '
+          print*,'sensor ',sensor
+          icx=0
+          icy=0
+       endif
  *     -----------------
  *     CLUSTER X
  *     -----------------
        if(icx.ne.0)then
-          viewx = VIEW(icx)
+          viewx   = VIEW(icx)
-          nldx = nld(MAXS(icx),VIEW(icx))
+          nldx    = nld(MAXS(icx),VIEW(icx))
-          nplx = npl(VIEW(icx))
+          nplx    = npl(VIEW(icx))
-          resxPAM = RESXAV
+ c         resxPAM = RESXAV
-          stripx = float(MAXS(icx))
+          stripx  = float(MAXS(icx))
+          if(
+      $        viewx.lt.1.or.
+      $        viewx.gt.12.or.
+      $        nldx.lt.1.or.
+      $        nldx.gt.3.or.
+      $        stripx.lt.1.or.
+      $        stripx.gt.3072.or.
+      $        .false.)then
+             print*,'xyz_PAM   ***ERROR*** wrong input '
+             print*,'icx ',icx,'view ',viewx,'nld ',nldx,'strip ',stripx
+             icx = 0
+             goto 10
+          endif
  *        --------------------------
  *        magnetic-field corrections
  *        --------------------------
- c$$$         print*,nplx,ax,bfy/10.
+          stripx  = stripx +  fieldcorr(viewx,bfy)
-          angtemp  = ax
+          angx    = effectiveangle(ax,viewx,bfy)
-          bfytemp  = bfy
-          if(nplx.eq.6) angtemp = -1. * ax
+          call applypfa(PFAx,icx,angx,corr,res)
-          if(nplx.eq.6) bfytemp = -1. * bfy
+          stripx  = stripx + corr
-          tgtemp = tan(angtemp*acos(-1.)/180.) + pmuH_h*bfytemp*0.00001
+          resxPAM = res
-          angx = 180.*atan(tgtemp)/acos(-1.)
-          stripx = stripx - 0.5*pmuH_h*bfytemp*0.00001*SiDimZ/pitchX
- c$$$         print*,angx,0.5*pmuH_h*bfytemp*0.00001*SiDimZ/pitchX
- c$$$         print*,'========================'
- *        --------------------------
-          if(PFAx.eq.'COG1')then
-             stripx = stripx
-             resxPAM = resxPAM
-          elseif(PFAx.eq.'COG2')then
-             stripx = stripx + cog(2,icx)
-             resxPAM = resxPAM*fbad_cog(2,icx)
-          elseif(PFAx.eq.'COG3')then
-             stripx = stripx + cog(3,icx)
-             resxPAM = resxPAM*fbad_cog(3,icx)
-          elseif(PFAx.eq.'COG4')then
- c            print*,'COG4'
-             stripx = stripx + cog(4,icx)
-             resxPAM = resxPAM*fbad_cog(4,icx)
-          elseif(PFAx.eq.'ETA2')then
-             stripx = stripx + pfaeta2(icx,angx)
-             resxPAM = risx_eta2(abs(angx))
-             if(DEBUG.and.fbad_cog(2,icx).ne.1)
-      $           print*,'BAD icx >>> ',viewx,fbad_cog(2,icx)
-             resxPAM = resxPAM*fbad_cog(2,icx)
-          elseif(PFAx.eq.'ETA3')then
-             stripx = stripx + pfaeta3(icx,angx)
-             resxPAM = risx_eta3(abs(angx))
-             if(DEBUG.and.fbad_cog(3,icx).ne.1)
-      $           print*,'BAD icx >>> ',viewx,fbad_cog(3,icx)
-             resxPAM = resxPAM*fbad_cog(3,icx)
-          elseif(PFAx.eq.'ETA4')then
-             stripx = stripx + pfaeta4(icx,angx)
-             resxPAM = risx_eta4(abs(angx))
-             if(DEBUG.and.fbad_cog(4,icx).ne.1)
-      $           print*,'BAD icx >>> ',viewx,fbad_cog(4,icx)
-             resxPAM = resxPAM*fbad_cog(4,icx)
-          elseif(PFAx.eq.'ETA')then
- c            print*,'ETA'
-             stripx = stripx + pfaeta(icx,angx)
-             resxPAM = ris_eta(icx,angx)
-             if(DEBUG.and.fbad_cog(2,icx).ne.1)
-      $           print*,'BAD icx >>> ',viewx,fbad_cog(2,icx)
-             resxPAM = resxPAM*fbad_eta(icx,angx)
-          elseif(PFAx.eq.'COG')then
-             stripx = stripx + cog(0,icx)
-             resxPAM = risx_cog(abs(angx))
-             resxPAM = resxPAM*fbad_cog(0,icx)
-          else
-             print*,'*** Non valid p.f.a. (x) --> ',PFAx
-          endif
- c         print*,'%%%%%%%%%%%%'
-       endif
+  endif
  *     -----------------
  *     CLUSTER Y
  *     -----------------
-Line 650 
 c         print*,'%%%%%%%%%%%%'
+Line 765 
 c         print*,'%%%%%%%%%%%%'
           viewy = VIEW(icy)
           nldy = nld(MAXS(icy),VIEW(icy))
           nply = npl(VIEW(icy))
-          resyPAM = RESYAV
+ c         resyPAM = RESYAV
           stripy = float(MAXS(icy))
+          if(
+      $        viewy.lt.1.or.
+      $        viewy.gt.12.or.
+      $        nldy.lt.1.or.
+      $        nldy.gt.3.or.
+      $        stripy.lt.1.or.
+      $        stripy.gt.3072.or.
+      $        .false.)then
+             print*,'xyz_PAM   ***ERROR*** wrong input '
+             print*,'icy ',icy,'view ',viewy,'nld ',nldy,'strip ',stripy
+             icy = 0
+             goto 20
+          endif
           if(icx.ne.0.and.(nply.ne.nplx.or.nldy.ne.nldx))then
-             print*,'xyz_PAM   ***ERROR*** invalid cluster couple!!! '
+             if(DEBUG.EQ.1) then
-      $           ,icx,icy
+                print*,'xyz_PAM   ***ERROR*** invalid cluster couple!!! '
+      $              ,icx,icy
+             endif
              goto 100
           endif
  *        --------------------------
  *        magnetic-field corrections
  *        --------------------------
-          tgtemp = tan(ay*acos(-1.)/180.)+pmuH_e*bfx*0.00001
+          stripy  = stripy +  fieldcorr(viewy,bfx)
-          angy    = 180.*atan(tgtemp)/acos(-1.)
+          angy    = effectiveangle(ay,viewy,bfx)
-          stripy = stripy + 0.5*pmuH_e*bfx*0.00001*SiDimZ/pitchY
- *        --------------------------
-          if(PFAy.eq.'COG1')then !(1)
+          call applypfa(PFAy,icy,angy,corr,res)
-             stripy = stripy     !(1)
+          stripy  = stripy + corr
-             resyPAM = resyPAM   !(1)
+          resyPAM = res
-          elseif(PFAy.eq.'COG2')then
-             stripy = stripy + cog(2,icy)
-             resyPAM = resyPAM*fbad_cog(2,icy)
-          elseif(PFAy.eq.'COG3')then
-             stripy = stripy + cog(3,icy)
-             resyPAM = resyPAM*fbad_cog(3,icy)
-          elseif(PFAy.eq.'COG4')then
-             stripy = stripy + cog(4,icy)
-             resyPAM = resyPAM*fbad_cog(4,icy)
-          elseif(PFAy.eq.'ETA2')then
- c            cog2 = cog(2,icy)
- c            etacorr = pfaeta2(cog2,viewy,nldy,angy)
- c            stripy = stripy + etacorr
-             stripy = stripy + pfaeta2(icy,angy)            !(3)
-             resyPAM = risy_eta2(abs(angy))                       !   (4)
-             resyPAM = resyPAM*fbad_cog(2,icy)
-             if(DEBUG.and.fbad_cog(2,icy).ne.1)
-      $           print*,'BAD icy >>> ',viewy,fbad_cog(2,icy)
-          elseif(PFAy.eq.'ETA3')then                         !(3)
-             stripy = stripy + pfaeta3(icy,angy)            !(3)
-             resyPAM = resyPAM*fbad_cog(3,icy)               !(3)
-             if(DEBUG.and.fbad_cog(3,icy).ne.1)              !(3)
-      $           print*,'BAD icy >>> ',viewy,fbad_cog(3,icy)!(3)
-          elseif(PFAy.eq.'ETA4')then                         !(3)
-             stripy = stripy + pfaeta4(icy,angy)            !(3)
-             resyPAM = resyPAM*fbad_cog(4,icy)               !(3)
-             if(DEBUG.and.fbad_cog(4,icy).ne.1)              !(3)
-      $           print*,'BAD icy >>> ',viewy,fbad_cog(4,icy)!(3)
-          elseif(PFAy.eq.'ETA')then                          !(3)
-             stripy = stripy + pfaeta(icy,angy)             !(3)
-             resyPAM = ris_eta(icy,angy)                     !   (4)
- c            resyPAM = resyPAM*fbad_cog(2,icy)              !(3)TEMPORANEO
-             resyPAM = resyPAM*fbad_eta(icy,angy)            !   (4)
-             if(DEBUG.and.fbad_cog(2,icy).ne.1)              !(3)
-      $           print*,'BAD icy >>> ',viewy,fbad_cog(2,icy)!(3)
-          elseif(PFAy.eq.'COG')then
-             stripy = stripy + cog(0,icy)
-             resyPAM = risy_cog(abs(angy))                        !   (4)
- c            resyPAM = ris_eta(icy,angy)                    !   (4)
-             resyPAM = resyPAM*fbad_cog(0,icy)
-          else
-             print*,'*** Non valid p.f.a. (x) --> ',PFAx
-          endif
-       endif
+  endif
- c      print*,'## stripx,stripy ',stripx,stripy
+ c$$$      print*,'## stripx,stripy ',stripx,stripy
  c===========================================================
  C     COUPLE
-Line 727 
 c     (xi,yi,zi) = mechanical coordinate
+Line 814 
 c     (xi,yi,zi) = mechanical coordinate
  c------------------------------------------------------------------------
           if(((mod(int(stripx+0.5)-1,1024)+1).le.3)
       $        .or.((mod(int(stripx+0.5)-1,1024)+1).ge.1022)) then !X has 1018 strips...
-             print*,'xyz_PAM (couple):',
+             if(DEBUG.EQ.1) then
-      $          ' WARNING: false X strip: strip ',stripx
+                print*,'xyz_PAM (couple):',
+      $              ' WARNING: false X strip: strip ',stripx
+             endif
           endif
           xi = acoordsi(stripx,viewx)
           yi = acoordsi(stripy,viewy)
           zi = 0.
  c------------------------------------------------------------------------
  c     (xrt,yrt,zrt) = rototranslated coordinates in the silicon sensor frame
  c------------------------------------------------------------------------
-Line 820 
 c            print*,'X-singlet ',icx,npl
+Line 908 
 c            print*,'X-singlet ',icx,npl
  c            if((stripx.le.3).or.(stripx.ge.1022)) then !X has 1018 strips...
              if(((mod(int(stripx+0.5)-1,1024)+1).le.3)
       $           .or.((mod(int(stripx+0.5)-1,1024)+1).ge.1022)) then !X has 1018 strips...
-                print*,'xyz_PAM (X-singlet):',
+                if(DEBUG.EQ.1) then
-      $             ' WARNING: false X strip: strip ',stripx
+                   print*,'xyz_PAM (X-singlet):',
+      $                 ' WARNING: false X strip: strip ',stripx
+                endif
              endif
              xi   = acoordsi(stripx,viewx)
-Line 843 
 c            print*,'X-cl ',icx,stripx,'
+Line 933 
 c            print*,'X-cl ',icx,stripx,'
  c            print*,yi_A,' <--> ',yi_B
           else
+             if(DEBUG.EQ.1) then
-             print *,'routine xyz_PAM ---> not properly used !!!'
+                print *,'routine xyz_PAM ---> not properly used !!!'
-             print *,'icx = ',icx
+                print *,'icx = ',icx
-             print *,'icy = ',icy
+                print *,'icy = ',icy
+             endif
              goto 100
           endif
-Line 911 
 c---------------------------------------
+Line 1002 
 c---------------------------------------
  c         print*,'A-(',xPAM_A,yPAM_A,') B-(',xPAM_B,yPAM_B,')'
        else
+          if(DEBUG.EQ.1) then
-          print *,'routine xyz_PAM ---> not properly used !!!'
+             print *,'routine xyz_PAM ---> not properly used !!!'
-          print *,'icx = ',icx
+             print *,'icx = ',icx
-          print *,'icy = ',icy
+             print *,'icy = ',icy
+          endif
        endif
-Line 926 
 c      print*,'## xPAM_B,yPAM_B,zPAM_B '
+Line 1017 
 c      print*,'## xPAM_B,yPAM_B,zPAM_B '
  continue
        end
+ ************************************************************************
+ *     Call xyz_PAM subroutine with default PFA and fill the mini2 common.
+ *     (done to be called from c/c++)
+ ************************************************************************
+       subroutine xyzpam(ip,icx,icy,lad,sensor,ax,ay,bfx,bfy)
+       include 'commontracker.f'
+       include 'level1.f'
+       include 'common_mini_2.f'
+       include 'common_xyzPAM.f'
+       include 'common_mech.f'
+       include 'calib.f'
+ *     flag to chose PFA
+ c$$$      character*10 PFA
+ c$$$      common/FINALPFA/PFA
+       integer icx,icy           !X-Y cluster ID
+       integer sensor
+       character*4 PFAx,PFAy     !PFA to be used
+       real ax,ay                !X-Y geometric angle
+       real bfx,bfy              !X-Y b-field components
+       ipx=0
+       ipy=0
+ c$$$      PFAx = 'COG4'!PFA
+ c$$$      PFAy = 'COG4'!PFA
+       if(icx.gt.nclstr1.or.icy.gt.nclstr1)then
+             print*,'xyzpam: ***WARNING*** clusters ',icx,icy
+      $           ,' does not exists (nclstr1=',nclstr1,')'
+             icx = -1*icx
+             icy = -1*icy
+             return
+       endif
+       call idtoc(pfaid,PFAx)
+       call idtoc(pfaid,PFAy)
+ cc      print*,PFAx,PFAy
+ c$$$      call xyz_PAM(icx,icy,sensor,PFAx,PFAy,ax,ay,bfx,bfy)
+ c$$$      print*,icx,icy,sensor,PFAx,PFAy,ax,ay,bfx,bfy
+       if(icx.ne.0.and.icy.ne.0)then
+          ipx=npl(VIEW(icx))
+          ipy=npl(VIEW(icy))
+ c$$$         if( (nplanes-ipx+1).ne.ip.or.(nplanes-ipy+1).ne.ip )
+ c$$$     $        print*,'xyzpam: ***WARNING*** clusters ',icx,icy
+ c$$$     $        ,' does not belong to the correct plane: ',ip,ipx,ipy
+          if( (nplanes-ipx+1).ne.ip )then
+             print*,'xyzpam: ***WARNING*** cluster ',icx
+      $           ,' does not belong to plane: ',ip
+             icx = -1*icx
+             return
+          endif
+          if( (nplanes-ipy+1).ne.ip )then
+             print*,'xyzpam: ***WARNING*** cluster ',icy
+      $           ,' does not belong to plane: ',ip
+             icy = -1*icy
+             return
+          endif
+          call xyz_PAM(icx,icy,sensor,PFAx,PFAy,ax,ay,bfx,bfy)
+          xgood(ip) = 1.
+          ygood(ip) = 1.
+          resx(ip)  = resxPAM
+          resy(ip)  = resyPAM
+          xm(ip) = xPAM
+          ym(ip) = yPAM
+          zm(ip) = zPAM
+          xm_A(ip) = 0.
+          ym_A(ip) = 0.
+          xm_B(ip) = 0.
+          ym_B(ip) = 0.
+ c         zv(ip) = zPAM
+       elseif(icx.eq.0.and.icy.ne.0)then
+          ipy=npl(VIEW(icy))
+ c$$$         if((nplanes-ipy+1).ne.ip)
+ c$$$     $        print*,'xyzpam: ***WARNING*** clusters ',icx,icy
+ c$$$     $        ,' does not belong to the correct plane: ',ip,ipx,ipy
+          if( (nplanes-ipy+1).ne.ip )then
+             print*,'xyzpam: ***WARNING*** cluster ',icy
+      $           ,' does not belong to plane: ',ip
+             icy = -1*icy
+             return
+          endif
+          call xyz_PAM(icx,icy,sensor,PFAx,PFAy,ax,ay,bfx,bfy)
+          xgood(ip) = 0.
+          ygood(ip) = 1.
+          resx(ip)  = 1000.
+          resy(ip)  = resyPAM
+          xm(ip) = -100.
+          ym(ip) = -100.
+          zm(ip) = (zPAM_A+zPAM_B)/2.
+          xm_A(ip) = xPAM_A
+          ym_A(ip) = yPAM_A
+          xm_B(ip) = xPAM_B
+          ym_B(ip) = yPAM_B
+ c         zv(ip) = (zPAM_A+zPAM_B)/2.
+       elseif(icx.ne.0.and.icy.eq.0)then
+          ipx=npl(VIEW(icx))
+ c$$$         if((nplanes-ipx+1).ne.ip)
+ c$$$     $        print*,'xyzpam: ***WARNING*** clusters ',icx,icy
+ c$$$     $        ,' does not belong to the correct plane: ',ip,ipx,ipy
+          if( (nplanes-ipx+1).ne.ip )then
+             print*,'xyzpam: ***WARNING*** cluster ',icx
+      $           ,' does not belong to plane: ',ip
+             icx = -1*icx
+             return
+          endif
+          call xyz_PAM(icx,icy,sensor,PFAx,PFAy,ax,ay,bfx,bfy)
+          xgood(ip) = 1.
+          ygood(ip) = 0.
+          resx(ip)  = resxPAM
+          resy(ip)  = 1000.
+          xm(ip) = -100.
+          ym(ip) = -100.
+          zm(ip) = (zPAM_A+zPAM_B)/2.
+          xm_A(ip) = xPAM_A
+          ym_A(ip) = yPAM_A
+          xm_B(ip) = xPAM_B
+          ym_B(ip) = yPAM_B
+ c         zv(ip) = (zPAM_A+zPAM_B)/2.
+       else
+          il = 2
+          if(lad.ne.0)il=lad
+          is = 1
+          if(sensor.ne.0)is=sensor
+ c         print*,nplanes-ip+1,il,is
+          xgood(ip) = 0.
+          ygood(ip) = 0.
+          resx(ip)  = 1000.
+          resy(ip)  = 1000.
+          xm(ip) = -100.
+          ym(ip) = -100.
+          zm(ip) = z_mech_sensor(nplanes-ip+1,il,is)*1000./1.d4
+          xm_A(ip) = 0.
+          ym_A(ip) = 0.
+          xm_B(ip) = 0.
+          ym_B(ip) = 0.
+ c         zv(ip) = z_mech_sensor(nplanes-ip+1,il,is)*1000./1.d4
+       endif
+       if(DEBUG.EQ.1)then
+ c         print*,'----------------------------- track coord'
+   format(i2,' * ',3f10.4,' --- ',4f10.4,' --- ',2f4.0,2f10.5)
+          write(*,22222)ip,zm(ip),xm(ip),ym(ip)
+      $        ,xm_A(ip),ym_A(ip),xm_B(ip),ym_B(ip)
+      $        ,xgood(ip),ygood(ip),resx(ip),resy(ip)
+ c$$$         print*,'-----------------------------'
+       endif
+       end
  ********************************************************************************
  ********************************************************************************
-Line 1001 
 c      print*,'## xPAM_B,yPAM_B,zPAM_B '
+Line 1274 
 c      print*,'## xPAM_B,yPAM_B,zPAM_B '
           endif
           distance=
-      $        ((xmi-XPP)**2+(ymi-YPP)**2)/RE**2
+      $       ((xmi-XPP)**2+(ymi-YPP)**2)!QUIQUI
+ cc     $        ((xmi-XPP)**2+(ymi-YPP)**2)/RE**2
           distance=dsqrt(distance)
  c$$$         print*,xPAM_A,yPAM_A,zPAM_A,xPAM_b,yPAM_b,zPAM_b
-Line 1026 
 c$$$         print*,' resolution ',re
+Line 1300 
 c$$$         print*,' resolution ',re
  *     ----------------------
           distance=
-      $        ((xPAM-XPP)/resxPAM)**2
+      $       ((xPAM-XPP))**2 !QUIQUI
-      $        +
+      $       +
-      $        ((yPAM-YPP)/resyPAM)**2
+      $       ((yPAM-YPP))**2
+ c$$$     $        ((xPAM-XPP)/resxPAM)**2
+ c$$$     $        +
+ c$$$     $        ((yPAM-YPP)/resyPAM)**2
           distance=dsqrt(distance)
  c$$$         print*,xPAM,yPAM,zPAM
-Line 1037 
 c$$$         print*,' resolution ',resxP
+Line 1314 
 c$$$         print*,' resolution ',resxP
        else
-          print*
+ c         print*
-      $        ,' function distance_to ---> wrong usage!!!'
+ c     $        ,' function distance_to ---> wrong usage!!!'
-          print*,' xPAM,yPAM,zPAM ',xPAM,yPAM,zPAM
+ c         print*,' xPAM,yPAM,zPAM ',xPAM,yPAM,zPAM
-          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 '
-      $        ,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
        endif
        distance_to = sngl(distance)
-Line 1109 
 c---------------------------------------
+Line 1386 
 c---------------------------------------
                 if(((mod(int(stripx+0.5)-1,1024)+1).le.3)
       $              .or.((mod(int(stripx+0.5)-1,1024)+1).ge.1022)) then !X has 1018 strips...
  c     if((stripx.le.3).or.(stripx.ge.1022)) then !X has 1018 strips...
-                   print*,'whichsensor: ',
+ c                  print*,'whichsensor: ',
-      $                ' WARNING: false X strip: strip ',stripx
+ c     $                ' WARNING: false X strip: strip ',stripx
                 endif
                 xi = acoordsi(stripx,viewx)
                 yi = acoordsi(stripy,viewy)
-Line 1265 
 c      include 'common_analysis.f'
+Line 1542 
 c      include 'common_analysis.f'
        is_cp=0
        if(id.lt.0)is_cp=1
        if(id.gt.0)is_cp=2
-       if(id.eq.0)print*,'IS_CP ===> wrong couple id !!!'
+ c      if(id.eq.0)print*,'IS_CP ===> wrong couple id !!!'
        return
        end
-Line 1364 
 c      include 'level1.f'
+Line 1641 
 c      include 'level1.f'
        integer iflag
        integer badseed,badclx,badcly
+       if(DEBUG.EQ.1)print*,'cl_to_couples:'
  *     init variables
        ncp_tot=0
-Line 1394 
 c      include 'level1.f'
+Line 1673 
 c      include 'level1.f'
  *     mask views with too many clusters
        do iv=1,nviews
           if( ncl_view(iv).gt. nclusterlimit)then
-             mask_view(iv) = 1
+ c            mask_view(iv) = 1
-             if(DEBUG)print*,' * WARNING * cl_to_couple: n.clusters > '
+             mask_view(iv) = mask_view(iv) + 2**0
-      $           ,nclusterlimit,' on view ', iv,' --> masked!'
+             if(DEBUG.EQ.1)
+      $        print*,' * WARNING * cl_to_couple: n.clusters > '
+      $        ,nclusterlimit,' on view ', iv,' --> masked!'
           endif
        enddo
-Line 1534 
 c                  cut = chcut * sch(npl
+Line 1815 
 c                  cut = chcut * sch(npl
                 endif
                 if(ncp_plane(nplx).gt.ncouplemax)then
-                   if(verbose)print*,
+                   if(verbose.eq.1)print*,
       $                 '** warning ** number of identified '//
       $                 'couples on plane ',nplx,
       $                 'exceeds vector dimention '
       $                 ,'( ',ncouplemax,' ) --> masked!'
-                   mask_view(nviewx(nplx)) = 2
+ c                  mask_view(nviewx(nplx)) = 2
-                   mask_view(nviewy(nply)) = 2
+ c                  mask_view(nviewy(nply)) = 2
+                   mask_view(nviewx(nplx))= mask_view(nviewx(nplx))+ 2**1
+                   mask_view(nviewy(nply))= mask_view(nviewy(nply))+ 2**1
                    goto 10
                 endif
-Line 1570 
 c                  cut = chcut * sch(npl
+Line 1853 
 c                  cut = chcut * sch(npl
        enddo
-       if(DEBUG)then
+       if(DEBUG.EQ.1)then
           print*,'clusters  ',nclstr1
           print*,'good    ',(cl_good(i),i=1,nclstr1)
-          print*,'singles ',(cl_single(i),i=1,nclstr1)
+          print*,'singlets',(cl_single(i),i=1,nclstr1)
           print*,'couples per plane: ',(ncp_plane(ip),ip=1,nplanes)
        endif
-Line 1633 
 c      double precision xm3,ym3,zm3
+Line 1916 
 c      double precision xm3,ym3,zm3
        real xc,zc,radius
  *     -----------------------------
+       if(DEBUG.EQ.1)print*,'cp_to_doubtrip:'
  *     --------------------------------------------
  *     put a limit to the maximum number of couples
-Line 1641 
 c      double precision xm3,ym3,zm3
+Line 1925 
 c      double precision xm3,ym3,zm3
  *     --------------------------------------------
        do ip=1,nplanes
           if(ncp_plane(ip).gt.ncouplelimit)then
-             mask_view(nviewx(ip)) = 8
+ c            mask_view(nviewx(ip)) = 8
-             mask_view(nviewy(ip)) = 8
+ c            mask_view(nviewy(ip)) = 8
+             mask_view(nviewx(ip)) = mask_view(nviewx(ip)) + 2**7
+             mask_view(nviewy(ip)) = mask_view(nviewy(ip)) + 2**7
           endif
        enddo
-Line 1668 
 c     print*,'***',is1,xm1,ym1,zm1
+Line 1954 
 c     print*,'***',is1,xm1,ym1,zm1
                 do ip2=(ip1+1),nplanes !loop on planes - COPPIA 2
                    if(  mask_view(nviewx(ip2)).ne.0 .or.
       $                 mask_view(nviewy(ip2)).ne.0 )goto 20 !skip plane
                    do is2=1,2    !loop on sensors -ndblt COPPIA 2
                       do icp2=1,ncp_plane(ip2) !loop on COPPIA 2
                          icx2=clx(ip2,icp2)
                          icy2=cly(ip2,icp2)
-Line 1688 
 c     $                       (icx2,icy2
+Line 1973 
 c     $                       (icx2,icy2
  *     (2 couples needed)
  *     - - - - - - - - - - - - - - - - - - - - - - - - - - - -
                          if(ndblt.eq.ndblt_max)then
-                            if(verbose)print*,
+                            if(verbose.eq.1)print*,
       $                          '** warning ** number of identified '//
       $                          'doublets exceeds vector dimention '
       $                          ,'( ',ndblt_max,' )'
  c                           good2=.false.
  c                           goto 880 !fill ntp and go to next event
                             do iv=1,12
-                               mask_view(iv) = 3
+ c                              mask_view(iv) = 3
+                               mask_view(iv) = mask_view(iv)+ 2**2
                             enddo
                             iflag=1
                             return
-Line 1751 
 c     $                               (i
+Line 2037 
 c     $                               (i
                                   ym3=yPAM
                                   zm3=zPAM
  *     find the circle passing through the three points
+ c$$$                                 call tricircle(3,xp,zp,angp,resp,chi
+ c$$$     $                                ,xc,zc,radius,iflag)
+                                  iflag = DEBUG
                                   call tricircle(3,xp,zp,angp,resp,chi
       $                                ,xc,zc,radius,iflag)
  c     print*,xc,zc,radius
-Line 1767 
 c     $
+Line 2056 
 c     $
  *     (3 couples needed)
  *     - - - - - - - - - - - - - - - - - - - - - - - - - - - -
                                   if(ntrpt.eq.ntrpt_max)then
-                                     if(verbose)print*,
+                                     if(verbose.eq.1)print*,
       $                     '** warning ** number of identified '//
       $                     'triplets exceeds vector dimention '
       $                    ,'( ',ntrpt_max,' )'
  c                                    good2=.false.
  c                                    goto 880 !fill ntp and go to next event
                                      do iv=1,nviews
-                                        mask_view(iv) = 4
+ c                                       mask_view(iv) = 4
+                                        mask_view(iv)=mask_view(iv)+ 2**3
                                      enddo
                                      iflag=1
                                      return
-Line 1828 
 c     print*,alfaxz1(ntrpt),alfaxz2(ntrp
+Line 2118 
 c     print*,alfaxz1(ntrpt),alfaxz2(ntrp
   continue
        enddo                     !end loop on planes  - COPPIA 1
-       if(DEBUG)then
+       if(DEBUG.EQ.1)then
           print*,'--- doublets ',ndblt
           print*,'--- triplets ',ntrpt
        endif
-Line 1875 
 c      include 'momanhough_init.f'
+Line 2165 
 c      include 'momanhough_init.f'
        integer cp_useds1(ncouplemaxtot) ! sensor 1
        integer cp_useds2(ncouplemaxtot) ! sensor 2
+       if(DEBUG.EQ.1)print*,'doub_to_YZcloud:'
  *~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  *     classification of DOUBLETS
-Line 2001 
 c         if(ncpused.lt.ncpyz_min)goto 2
+Line 2292 
 c         if(ncpused.lt.ncpyz_min)goto 2
  *     >>> NEW CLOUD <<<
           if(nclouds_yz.ge.ncloyz_max)then
-             if(verbose)print*,
+             if(verbose.eq.1)print*,
       $           '** warning ** number of identified '//
       $           'YZ clouds exceeds vector dimention '
       $           ,'( ',ncloyz_max,' )'
  c               good2=.false.
  c     goto 880         !fill ntp and go to next event
              do iv=1,nviews
-                mask_view(iv) = 5
+ c               mask_view(iv) = 5
+                mask_view(iv) = mask_view(iv) + 2**4
              enddo
              iflag=1
              return
-Line 2028 
 c     ptcloud_yz_nt(nclouds_yz)=npt
+Line 2320 
 c     ptcloud_yz_nt(nclouds_yz)=npt
  c     print*,'>> ',ipt,db_all(ipt)
           enddo
           npt_tot=npt_tot+npt
-          if(DEBUG)then
+          if(DEBUG.EQ.1)then
              print*,'-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~'
              print*,'>>>> cloud ',nclouds_yz,' --- ',npt,' points'
-             print*,'- alfayz1 ',alfayz1_av(nclouds_yz)
+             print*,'- alfayz1  ',alfayz1_av(nclouds_yz)
-             print*,'- alfayz2 ',alfayz2_av(nclouds_yz)
+             print*,'- alfayz2  ',alfayz2_av(nclouds_yz)
-             print*,'cp_useds1 ',(cp_useds1(icp),icp=1,ncp_tot)
+             print*,'cp_useds1  ',(cp_useds1(icp),icp=1,ncp_tot)
-             print*,'cp_useds2 ',(cp_useds2(icp),icp=1,ncp_tot)
+             print*,'cp_useds2  ',(cp_useds2(icp),icp=1,ncp_tot)
-             print*,'hit_plane ',(hit_plane(ip),ip=1,nplanes)
+             print*,'cpcloud_yz '
+      $           ,(cpcloud_yz(nclouds_yz,icp),icp=1,ncp_tot)
+             print*,'hit_plane  ',(hit_plane(ip),ip=1,nplanes)
  c$$$            print*,'nt-uple: ptcloud_yz(',nclouds_yz,') = '
  c$$$     $           ,ptcloud_yz(nclouds_yz)
  c$$$            print*,'nt-uple: db_cloud(...) = '
-Line 2053 
 c$$$     $           ,(db_cloud(iii),iii
+Line 2347 
 c$$$     $           ,(db_cloud(iii),iii
          goto 90
        endif
-       if(DEBUG)then
+       if(DEBUG.EQ.1)then
           print*,'---------------------- '
           print*,'Y-Z total clouds ',nclouds_yz
           print*,' '
-Line 2102 
 c      include 'momanhough_init.f'
+Line 2396 
 c      include 'momanhough_init.f'
        integer cp_useds1(ncouplemaxtot) ! sensor 1
        integer cp_useds2(ncouplemaxtot) ! sensor 2
+       if(DEBUG.EQ.1)print*,'trip_to_XZcloud:'
  *~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  *     classification of TRIPLETS
  *     according to distance in parameter space
-Line 2159 
 c         tr_temp(1)=itr1
+Line 2455 
 c         tr_temp(1)=itr1
       $              +((alfaxz2(itrref)-alfaxz2(itr2))/Dalfaxz2)**2
                 distance = sqrt(distance)
-                if(distance.lt.cutdistxz)then
+ *     ------------------------------------------------------------------------
+ *     FORCE INCLUSION OF TRIPLETS COMPOSED BY SAME COUPLES, IGNORING THE IMAGE
+ *     ------------------------------------------------------------------------
+ *     (added in august 2007)
+                istrimage=0
+                if(
+      $              abs(cpxz1(itrref)).eq.abs(cpxz1(itr2)).and.
+      $              abs(cpxz2(itrref)).eq.abs(cpxz2(itr2)).and.
+      $              abs(cpxz3(itrref)).eq.abs(cpxz3(itr2)).and.
+      $              .true.)istrimage=1
+                if(distance.lt.cutdistxz.or.istrimage.eq.1)then
  c     print*,idb1,idb2,distance,' cloud ',nclouds_yz
                    if(cpxz1(itr2).gt.0)cp_useds2(cpxz1(itr2))=1
                    if(cpxz1(itr2).lt.0)cp_useds1(-cpxz1(itr2))=1
-Line 2218 
 c     print*,'check cp_used'
+Line 2525 
 c     print*,'check cp_used'
           enddo
  c         if(ncpused.lt.ncpxz_min)goto 22288 !next triplet
           if(npt.lt.nptxz_min)goto 22288     !next triplet
-          if(nplused.lt.nplxz_min)goto 22288 !next doublet
+          if(nplused.lt.nplxz_min)goto 22288 !next triplet
  *     ~~~~~~~~~~~~~~~~~
  *     >>> NEW CLOUD <<<
           if(nclouds_xz.ge.ncloxz_max)then
-             if(verbose)print*,
+             if(verbose.eq.1)print*,
       $           '** warning ** number of identified '//
       $           'XZ clouds exceeds vector dimention '
       $           ,'( ',ncloxz_max,' )'
  c     good2=.false.
  c     goto 880         !fill ntp and go to next event
              do iv=1,nviews
-                mask_view(iv) = 6
+ c               mask_view(iv) = 6
+                mask_view(iv) =  mask_view(iv) + 2**5
              enddo
              iflag=1
              return
-Line 2249 
 c     goto 880         !fill ntp and go
+Line 2557 
 c     goto 880         !fill ntp and go
           enddo
           npt_tot=npt_tot+npt
-          if(DEBUG)then
+          if(DEBUG.EQ.1)then
              print*,'-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~'
              print*,'>>>> cloud ',nclouds_xz,' --- ',npt,' points'
-             print*,'- alfaxz1 ',alfaxz1_av(nclouds_xz)
+             print*,'- alfaxz1  ',alfaxz1_av(nclouds_xz)
-             print*,'- alfaxz2 ',alfaxz2_av(nclouds_xz)
+             print*,'- alfaxz2  ',alfaxz2_av(nclouds_xz)
-             print*,'- alfaxz3 ',alfaxz3_av(nclouds_xz)
+             print*,'- alfaxz3  ',alfaxz3_av(nclouds_xz)
-             print*,'cp_useds1 ',(cp_useds1(icp),icp=1,ncp_tot)
+             print*,'cp_useds1  ',(cp_useds1(icp),icp=1,ncp_tot)
-             print*,'cp_useds2 ',(cp_useds2(icp),icp=1,ncp_tot)
+             print*,'cp_useds2  ',(cp_useds2(icp),icp=1,ncp_tot)
+             print*,'cpcloud_xz '
+      $           ,(cpcloud_xz(nclouds_xz,icp),icp=1,ncp_tot)
              print*,'hit_plane ',(hit_plane(ip),ip=1,nplanes)
  c$$$            print*,'nt-uple: ptcloud_xz(',nclouds_xz,') = '
  c$$$     $           ,ptcloud_xz(nclouds_xz)
-Line 2274 
 c$$$     $           ,(tr_cloud(iii),iii
+Line 2584 
 c$$$     $           ,(tr_cloud(iii),iii
           goto 91
         endif
-       if(DEBUG)then
+       if(DEBUG.EQ.1)then
           print*,'---------------------- '
           print*,'X-Z total clouds ',nclouds_xz
           print*,' '
-Line 2295 
 c$$$     $           ,(tr_cloud(iii),iii
+Line 2605 
 c$$$     $           ,(tr_cloud(iii),iii
  **************************************************
        subroutine clouds_to_ctrack(iflag)
- c*****************************************************
- c     02/02/2006 modified by Elena Vannuccini --> (1)
- c*****************************************************
        include 'commontracker.f'
        include 'level1.f'
-Line 2305 
 c***************************************
+Line 2612 
 c***************************************
        include 'common_xyzPAM.f'
        include 'common_mini_2.f'
        include 'common_mech.f'
- c      include 'momanhough_init.f'
  *     output flag
-Line 2329 
 c      include 'momanhough_init.f'
+Line 2636 
 c      include 'momanhough_init.f'
  *     -----------------------------------------------------------
  *     variables for track fitting
        double precision AL_INI(5)
- c      double precision tath
  *     -----------------------------------------------------------
- c      real fitz(nplanes)        !z coordinates of the planes in cm
+       if(DEBUG.EQ.1)print*,'clouds_to_ctrack:'
        ntracks=0                 !counter of track candidates
-Line 2354 
 c      real fitz(nplanes)        !z coor
+Line 2660 
 c      real fitz(nplanes)        !z coor
                 enddo
              enddo
              ncp_ok=0
-             do icp=1,ncp_tot    !loop on couples
+             do icp=1,ncp_tot    !loop over couples
  *     get info on
                 cpintersec(icp)=min(
       $              cpcloud_yz(iyz,icp),
-Line 2363 
 c      real fitz(nplanes)        !z coor
+Line 2669 
 c      real fitz(nplanes)        !z coor
       $    (cpcloud_yz(iyz,icp).eq.1.and.cpcloud_xz(ixz,icp).eq.2).or.
       $    (cpcloud_yz(iyz,icp).eq.2.and.cpcloud_xz(ixz,icp).eq.1).or.
       $              .false.)cpintersec(icp)=0
+ *     cpintersec is >0 if yz and xz clouds contain the same image of couple icp
                 if(cpintersec(icp).ne.0)then
                    ncp_ok=ncp_ok+1
-Line 2395 
 c      real fitz(nplanes)        !z coor
+Line 2702 
 c      real fitz(nplanes)        !z coor
                 nplused=nplused+ hit_plane(ip)
              enddo
- c            if(nplused.lt.nplxz_min)goto 888 !next doublet
-             if(nplused.lt.nplyz_min)goto 888 !next doublet
-             if(ncp_ok.lt.ncpok)goto 888 !next cloud
-             if(DEBUG)then
+             if(DEBUG.EQ.1)then
                 print*,'Combination ',iyz,ixz
       $              ,' db ',ptcloud_yz(iyz)
       $              ,' tr ',ptcloud_xz(ixz)
       $              ,'  -----> # matching couples ',ncp_ok
              endif
+ c            if(nplused.lt.nplxz_min)goto 888 !next combination
+             if(nplused.lt.nplyz_min)goto 888 !next combination
+             if(ncp_ok.lt.ncpok)goto 888 !next combination
  c$$$  print*,'~~~~~~~~~~~~~~~~~~~~~~~~~'
  c$$$  print*,'Configurazione cluster XZ'
  c$$$  print*,'1 -- ',(clx(1,i),i=1,ncp_plane(1))
-Line 2450 
 c$$$            AL_INI(5) = (1.e2*alfaxz
+Line 2759 
 c$$$            AL_INI(5) = (1.e2*alfaxz
  c$$$
  c$$$            if(AL_INI(5).gt.defmax)goto 888 !next cloud
-             if(DEBUG)then
+             if(DEBUG.EQ.1)then
+                print*,'track candidate', ntracks+1
                 print*,'1 >>> ',(cp_match(6,i),i=1,ncp_match(6))
                 print*,'2 >>> ',(cp_match(5,i),i=1,ncp_match(5))
                 print*,'3 >>> ',(cp_match(4,i),i=1,ncp_match(4))
-Line 2483 
 c$$$            if(AL_INI(5).gt.defmax)g
+Line 2793 
 c$$$            if(AL_INI(5).gt.defmax)g
                                hit_plane(6)=icp6
                                if(ncp_match(6).eq.0)hit_plane(6)=0 !-icp6
+ *                             ---------------------------------------
+ *                             check if this group of couples has been
+ *                             already fitted
+ *                             ---------------------------------------
+                               do ica=1,ntracks
+                                  isthesame=1
+                                  do ip=1,NPLANES
+                                     if(hit_plane(ip).ne.0)then
+                                        if(  CP_STORE(nplanes-ip+1,ica)
+      $                                      .ne.
+      $                                      cp_match(ip,hit_plane(ip)) )
+      $                                      isthesame=0
+                                     else
+                                        if(  CP_STORE(nplanes-ip+1,ica)
+      $                                      .ne.
+      $                                      0 )
+      $                                      isthesame=0
+                                     endif
+                                  enddo
+                                  if(isthesame.eq.1)then
+                                     if(DEBUG.eq.1)
+      $                                   print*,'(already fitted)'
+                                     goto 666 !jump to next combination
+                                  endif
+                               enddo
                                call track_init !init TRACK common
-                               do ip=1,nplanes !loop on planes
+                               do ip=1,nplanes !loop on planes (bottom to top)
                                   if(hit_plane(ip).ne.0)then
                                      id=cp_match(ip,hit_plane(ip))
                                      is=is_cp(id)
-Line 2530 
 c$$$                              enddo
+Line 2865 
 c$$$                              enddo
                                ifail=0 !error flag in chi^2 computation
                                jstep=0 !number of  minimization steps
                                iprint=0
- c                              if(DEBUG)iprint=1
+ c                              if(DEBUG.EQ.1)iprint=1
-                               if(DEBUG)iprint=2
+                               if(DEBUG.EQ.1)iprint=2
                                call mini2(jstep,ifail,iprint)
                                if(ifail.ne.0) then
-                                  if(DEBUG)then
+                                  if(DEBUG.EQ.1)then
                                      print *,
       $                              '*** MINIMIZATION FAILURE *** '
       $                              //'(clouds_to_ctrack)'
-Line 2559 
 c                                 chi2=-
+Line 2894 
 c                                 chi2=-
  *     --------------------------
                                if(ntracks.eq.NTRACKSMAX)then
-                                  if(verbose)print*,
+                                  if(verbose.eq.1)print*,
       $                 '** warning ** number of candidate tracks '//
       $                 ' exceeds vector dimension '
       $                ,'( ',NTRACKSMAX,' )'
  c                                 good2=.false.
  c                                 goto 880 !fill ntp and go to next event
                                   do iv=1,nviews
-                                     mask_view(iv) = 7
+ c                                    mask_view(iv) = 7
+                                     mask_view(iv) = mask_view(iv) + 2**6
                                   enddo
                                   iflag=1
                                   return
-Line 2574 
 c                                 goto 8
+Line 2910 
 c                                 goto 8
                                ntracks = ntracks + 1
- c$$$                              ndof=0
+                               do ip=1,nplanes !top to bottom
-                               do ip=1,nplanes
- c$$$                                 ndof=ndof
- c$$$     $                                +int(xgood(ip))
- c$$$     $                                +int(ygood(ip))
                                   XV_STORE(ip,ntracks)=sngl(xv(ip))
                                   YV_STORE(ip,ntracks)=sngl(yv(ip))
                                   ZV_STORE(ip,ntracks)=sngl(zv(ip))
-Line 2594 
 c$$$     $
+Line 2927 
 c$$$     $
                                   AYV_STORE(ip,ntracks)=sngl(ayv(ip))
                                   XGOOD_STORE(ip,ntracks)=sngl(xgood(ip))
                                   YGOOD_STORE(ip,ntracks)=sngl(ygood(ip))
+ *                                NB! hit_plane is defined from bottom to top
                                   if(hit_plane(ip).ne.0)then
                                      CP_STORE(nplanes-ip+1,ntracks)=
       $                                   cp_match(ip,hit_plane(ip))
+                                     SENSOR_STORE(nplanes-ip+1,ntracks)
+      $                              = is_cp(cp_match(ip,hit_plane(ip)))
+                                     LADDER_STORE(nplanes-ip+1,ntracks)
+      $                                   = LADDER(
+      $                                   clx(ip,icp_cp(
+      $                                   cp_match(ip,hit_plane(ip)
+      $                                   ))));
                                   else
                                      CP_STORE(nplanes-ip+1,ntracks)=0
+                                     SENSOR_STORE(nplanes-ip+1,ntracks)=0
+                                     LADDER_STORE(nplanes-ip+1,ntracks)=0
                                   endif
-                                  CLS_STORE(nplanes-ip+1,ntracks)=0
+                                  BX_STORE(ip,ntracks)=0!I dont need it now
+                                  BY_STORE(ip,ntracks)=0!I dont need it now
+                                  CLS_STORE(ip,ntracks)=0
                                   do i=1,5
                                      AL_STORE(i,ntracks)=sngl(AL(i))
                                   enddo
                                enddo
- c$$$  *                             Number of Degree Of Freedom
- c$$$  ndof=ndof-5
- c$$$  *                             reduced chi^2
- c$$$  rchi2=chi2/dble(ndof)
                                RCHI2_STORE(ntracks)=chi2
  *     --------------------------------
-Line 2633 
 c$$$  rchi2=chi2/dble(ndof)
+Line 2974 
 c$$$  rchi2=chi2/dble(ndof)
           return
        endif
-       if(DEBUG)then
+ c$$$      if(DEBUG.EQ.1)then
-          print*,'****** TRACK CANDIDATES ***********'
+ c$$$         print*,'****** TRACK CANDIDATES ***********'
-          print*,'#         R. chi2        RIG'
+ c$$$         print*,'#         R. chi2        RIG'
-          do i=1,ntracks
+ c$$$         do i=1,ntracks
-             print*,i,' --- ',rchi2_store(i),' --- '
+ c$$$            print*,i,' --- ',rchi2_store(i),' --- '
-      $           ,1./abs(AL_STORE(5,i))
+ c$$$     $           ,1./abs(AL_STORE(5,i))
-          enddo
+ c$$$         enddo
-          print*,'***********************************'
+ c$$$         print*,'***********************************'
+ c$$$      endif
+       if(DEBUG.EQ.1)then
+         print*,'****** TRACK CANDIDATES *****************'
+         print*,'#         R. chi2        RIG         ndof'
+         do i=1,ntracks
+           ndof=0                !(1)
+           do ii=1,nplanes       !(1)
+             ndof=ndof           !(1)
+      $           +int(xgood_store(ii,i)) !(1)
+      $           +int(ygood_store(ii,i)) !(1)
+           enddo                 !(1)
+           print*,i,' --- ',rchi2_store(i),' --- '
+      $         ,1./abs(AL_STORE(5,i)),' --- ',ndof
+         enddo
+         print*,'*****************************************'
        endif
-Line 2659 
 c$$$  rchi2=chi2/dble(ndof)
+Line 3015 
 c$$$  rchi2=chi2/dble(ndof)
        subroutine refine_track(ibest)
- c******************************************************
- cccccc 06/10/2005 modified by elena vannuccini ---> (1)
- cccccc 31/01/2006 modified by elena vannuccini ---> (2)
- cccccc 12/08/2006 modified by elena vannucicni ---> (3)
- c******************************************************
        include 'commontracker.f'
        include 'level1.f'
-Line 2671 
 c***************************************
+Line 3022 
 c***************************************
        include 'common_xyzPAM.f'
        include 'common_mini_2.f'
        include 'common_mech.f'
- c      include 'momanhough_init.f'
- c      include 'level1.f'
        include 'calib.f'
  *     flag to chose PFA
        character*10 PFA
        common/FINALPFA/PFA
+       real k(6)
+       DATA k/1.099730,0.418900,0.220939,0.220907,0.418771,1.100674/
        real xp,yp,zp
        real xyzp(3),bxyz(3)
        equivalence (xp,xyzp(1)),(yp,xyzp(2)),(zp,xyzp(3))
+       if(DEBUG.EQ.1)print*,'refine_track:'
  *     =================================================
  *     new estimate of positions using ETA algorithm
  *                          and
-Line 2695 
 c      include 'level1.f'
+Line 3048 
 c      include 'level1.f'
           yP=YV_STORE(nplanes-ip+1,ibest)
           zP=ZV_STORE(nplanes-ip+1,ibest)
           call gufld(xyzp,bxyz)
- c$$$         bxyz(1)=0
+          BX_STORE(nplanes-ip+1,ibest)=bxyz(1)
+          BY_STORE(nplanes-ip+1,ibest)=bxyz(2)
+ c$$$  bxyz(1)=0
  c$$$         bxyz(2)=0
  c$$$         bxyz(3)=0
  *     |||||||||||||||||||||||||||||||||||||||||||||||||
-Line 2729 
 c     $           AYV_STORE(nplanes-ip+1
+Line 3084 
 c     $           AYV_STORE(nplanes-ip+1
       $           bxyz(1),
       $           bxyz(2)
       $           )
- c$$$  call xyz_PAM(icx,icy,is,
- c$$$  $              'COG2','COG2',
- c$$$  $              0.,
- c$$$  $              0.)
              xm(nplanes-ip+1) = xPAM
              ym(nplanes-ip+1) = yPAM
              zm(nplanes-ip+1) = zPAM
-Line 2741 
 c$$$  $              0.)
+Line 3093 
 c$$$  $              0.)
              resx(nplanes-ip+1) = resxPAM
              resy(nplanes-ip+1) = resyPAM
- c            dedxtrk(nplanes-ip+1) = (sgnl(icx)+sgnl(icy))/2. !(1)
+             dedxtrk_x(nplanes-ip+1)=sgnl(icx)/mip(VIEW(icx),LADDER(icx))
-             dedxtrk_x(nplanes-ip+1)=sgnl(icx)/mip(VIEW(icx),LADDER(icx)) !(1)(2)
+             dedxtrk_y(nplanes-ip+1)=sgnl(icy)/mip(VIEW(icy),LADDER(icy))
-             dedxtrk_y(nplanes-ip+1)=sgnl(icy)/mip(VIEW(icy),LADDER(icy)) !(1)(2)
  *     |||||||||||||||||||||||||||||||||||||||||||||||||
  *     -------------------------------------------------
-Line 2758 
 c            dedxtrk(nplanes-ip+1) = (sg
+Line 3109 
 c            dedxtrk(nplanes-ip+1) = (sg
  *     --------------------------------------------------------------
  *     determine which ladder and sensor are intersected by the track
- c$$$            xP=XV_STORE(nplanes-ip+1,ibest)
- c$$$            yP=YV_STORE(nplanes-ip+1,ibest)
- c$$$            zP=ZV_STORE(nplanes-ip+1,ibest)
              call whichsensor(ip,xP,yP,nldt,ist)
  *     if the track hit the plane in a dead area, go to the next plane
              if(nldt.eq.0.or.ist.eq.0)goto 133
+             SENSOR_STORE(nplanes-ip+1,IBEST)=ist
+             LADDER_STORE(nplanes-ip+1,IBEST)=nldt
  *     --------------------------------------------------------------
-             if(DEBUG)then
+             if(DEBUG.EQ.1)then
                 print*,
       $              '------ Plane ',ip,' intersected on LADDER ',nldt
       $              ,' SENSOR ',ist
-Line 2777 
 c$$$            zP=ZV_STORE(nplanes-ip+1
+Line 3128 
 c$$$            zP=ZV_STORE(nplanes-ip+1
  *     ===========================================
  *     STEP 1 >>>>>>>  try to include a new couple
  *     ===========================================
- c            if(DEBUG)print*,'>>>> try to include a new couple'
+ c            if(DEBUG.EQ.1)print*,'>>>> try to include a new couple'
              distmin=1000000.
              xmm = 0.
              ymm = 0.
-Line 2798 
 c     $              cl_used(icy).eq.1.o
+Line 3149 
 c     $              cl_used(icy).eq.1.o
       $              cl_used(icy).ne.0.or. !or the Y cluster is already used !(3)
       $              .false.)goto 1188 !then jump to next couple.
  *
- c               call xyz_PAM(icx,icy,ist,
- c     $              PFA,PFA,
- c     $              AXV_STORE(nplanes-ip+1,ibest),
- c     $              AYV_STORE(nplanes-ip+1,ibest))
                 call xyz_PAM(icx,icy,ist,
       $              PFA,PFA,
       $              AXV_STORE(nplanes-ip+1,ibest),
-Line 2811 
 c     $              AYV_STORE(nplanes-i
+Line 3158 
 c     $              AYV_STORE(nplanes-i
       $              )
                 distance = distance_to(XP,YP)
-                distance = distance / RCHI2_STORE(ibest)!<<< MS
+ c               distance = distance / RCHI2_STORE(ibest)!<<< MS !QUIQUI
                 id=id_cp(ip,icp,ist)
-                if(DEBUG)print*,'( couple ',id
+                if(DEBUG.EQ.1)print*,'( couple ',id
-      $              ,' ) normalized distance ',distance
+      $              ,' ) distance ',distance
                 if(distance.lt.distmin)then
                    xmm = xPAM
                    ymm = yPAM
-Line 2823 
 c     $              AYV_STORE(nplanes-i
+Line 3170 
 c     $              AYV_STORE(nplanes-i
                    rymm = resyPAM
                    distmin = distance
                    idm = id
- c                 dedxmm = (sgnl(icx)+sgnl(icy))/2. !(1)
                    dedxmmx = sgnl(icx)/mip(VIEW(icx),LADDER(icx)) !(1)(2)
                    dedxmmy = sgnl(icy)/mip(VIEW(icy),LADDER(icy)) !(1)(2)
+ c     QUIQUI --> non devo moltiplicare per clinc?!?!?!
+                   clincnewc=10*sqrt(rymm**2+rxmm**2 !QUIQUI
+      $                 +RCHI2_STORE(ibest)*k(ip)*(cov(1,1)+cov(2,2))) !QUIQUI
                 endif
          continue
              enddo               !end loop on couples on plane icp
-             if(distmin.le.clinc)then
+ c            if(distmin.le.clinc)then     !QUIQUI
+             if(distmin.le.clincnewc)then     !QUIQUI
  *              -----------------------------------
-                xm(nplanes-ip+1) = xmm         !<<<
+                xm(nplanes-ip+1) = xmm !<<<
-                ym(nplanes-ip+1) = ymm         !<<<
+                ym(nplanes-ip+1) = ymm !<<<
-                zm(nplanes-ip+1) = zmm         !<<<
+                zm(nplanes-ip+1) = zmm !<<<
-                xgood(nplanes-ip+1) = 1        !<<<
+                xgood(nplanes-ip+1) = 1 !<<<
-                ygood(nplanes-ip+1) = 1        !<<<
+                ygood(nplanes-ip+1) = 1 !<<<
-                resx(nplanes-ip+1)=rxmm        !<<<
+                resx(nplanes-ip+1)=rxmm !<<<
-                resy(nplanes-ip+1)=rymm        !<<<
+                resy(nplanes-ip+1)=rymm !<<<
- c              dedxtrk(nplanes-ip+1) = dedxmm !<<<  !(1)
+                dedxtrk_x(nplanes-ip+1) = dedxmmx !<<<
-                dedxtrk_x(nplanes-ip+1) = dedxmmx    !(1)
+                dedxtrk_y(nplanes-ip+1) = dedxmmy !<<<
-                dedxtrk_y(nplanes-ip+1) = dedxmmy    !(1)
  *              -----------------------------------
                 CP_STORE(nplanes-ip+1,ibest)=idm
-                if(DEBUG)print*,'%%%% included couple ',idm
+                if(DEBUG.EQ.1)print*,'%%%% included couple ',idm
-      $              ,' (norm.dist.= ',distmin,', cut ',clinc,' )'
+      $              ,' (dist.= ',distmin,', cut ',clinc,' )'
                 goto 133         !next plane
              endif
  *     ================================================
  *     STEP 2 >>>>>>>  try to include a single cluster
  *                     either from a couple or single
  *     ================================================
- c            if(DEBUG)print*,'>>>> try to include a new cluster'
+ c            if(DEBUG.EQ.1)print*,'>>>> try to include a new cluster'
              distmin=1000000.
              xmm_A = 0.          !---------------------------
              ymm_A = 0.          ! init variables that
-Line 2887 
 c     $              AXV_STORE(nplanes-i
+Line 3236 
 c     $              AXV_STORE(nplanes-i
       $              bxyz(2)
       $              )
                 distance = distance_to(XP,YP)
-                distance = distance / RCHI2_STORE(ibest)!<<< MS
+ c               distance = distance / RCHI2_STORE(ibest)!<<< MS !QUIQUI
-                if(DEBUG)print*,'( cl-X ',icx
+                if(DEBUG.EQ.1)print*,'( cl-X ',icx
-      $              ,' in cp ',id,' ) normalized distance ',distance
+      $              ,' in cp ',id,' ) distance ',distance
                 if(distance.lt.distmin)then
                    xmm_A = xPAM_A
                    ymm_A = yPAM_A
-Line 2920 
 c     $              0.,AYV_STORE(nplane
+Line 3269 
 c     $              0.,AYV_STORE(nplane
       $              bxyz(2)
       $              )
                 distance = distance_to(XP,YP)
-                distance = distance / RCHI2_STORE(ibest)!<<< MS
+ c               distance = distance / RCHI2_STORE(ibest)!<<< MS !QUIQUI
-                if(DEBUG)print*,'( cl-Y ',icy
+                if(DEBUG.EQ.1)print*,'( cl-Y ',icy
-      $              ,' in cp ',id,' ) normalized distance ',distance
+      $              ,' in cp ',id,' ) distance ',distance
                 if(distance.lt.distmin)then
                    xmm_A = xPAM_A
                    ymm_A = yPAM_A
-Line 2944 
 c                 dedxmm = sgnl(icy)  !(
+Line 3293 
 c                 dedxmm = sgnl(icy)  !(
  c            print*,'## ncls(',ip,') ',ncls(ip)
              do ic=1,ncls(ip)    !loop on single clusters
                 icl=cls(ip,ic)
- c              print*,'## ic ',ic,' ist ',ist
  c               if(cl_used(icl).eq.1.or.     !if the cluster is already used
                 if(cl_used(icl).ne.0.or.     !if the cluster is already used !(3)
       $              LADDER(icl).ne.nldt.or. !or the ladder number does not match
       $              .false.)goto 18882      !jump to the next singlet
                 if(mod(VIEW(icl),2).eq.0)then!<---- X view
- c                  call xyz_PAM(icl,0,ist,
- c     $                 PFA,PFA,
- c     $                 AXV_STORE(nplanes-ip+1,ibest),0.)
                    call xyz_PAM(icl,0,ist,
       $                 PFA,PFA,
       $                 AXV_STORE(nplanes-ip+1,ibest),0.,
-Line 2960 
 c     $                 AXV_STORE(nplane
+Line 3305 
 c     $                 AXV_STORE(nplane
       $                 bxyz(2)
       $                 )
                 else                         !<---- Y view
- c                  call xyz_PAM(0,icl,ist,
- c     $                 PFA,PFA,
- c     $                 0.,AYV_STORE(nplanes-ip+1,ibest))
                    call xyz_PAM(0,icl,ist,
       $                 PFA,PFA,
       $                 0.,AYV_STORE(nplanes-ip+1,ibest),
-Line 2972 
 c     $                 0.,AYV_STORE(npl
+Line 3314 
 c     $                 0.,AYV_STORE(npl
                 endif
                 distance = distance_to(XP,YP)
-                distance = distance / RCHI2_STORE(ibest)!<<< MS
+ c               distance = distance / RCHI2_STORE(ibest)!<<< MS !QUIQUI
-                if(DEBUG)print*,'( cl-s ',icl
+                if(DEBUG.EQ.1)print*,'( cl-s ',icl
-      $              ,' ) normalized distance ',distance,'<',distmin,' ?'
+      $              ,' ) distance ',distance
                 if(distance.lt.distmin)then
-                   if(DEBUG)print*,'YES'
+ c                  if(DEBUG.EQ.1)print*,'YES'
                    xmm_A = xPAM_A
                    ymm_A = yPAM_A
                    zmm_A = zPAM_A
-Line 2987 
 c     $                 0.,AYV_STORE(npl
+Line 3329 
 c     $                 0.,AYV_STORE(npl
                    rymm = resyPAM
                    distmin = distance
                    iclm = icl
- c                  dedxmm = sgnl(icl)                   !(1)
                    if(mod(VIEW(icl),2).eq.0)then !<---- X view
-                      dedxmmx = sgnl(icl)/mip(VIEW(icl),LADDER(icl)) !(1)(2)
+                      dedxmmx = sgnl(icl)/mip(VIEW(icl),LADDER(icl))
-                      dedxmmy = 0.                       !(1)
+                      dedxmmy = 0.
                    else          !<---- Y view
-                      dedxmmx = 0.                       !(1)
+                      dedxmmx = 0.
-                      dedxmmy = sgnl(icl)/mip(VIEW(icl),LADDER(icl)) !(1)(2)
+                      dedxmmy = sgnl(icl)/mip(VIEW(icl),LADDER(icl))
                    endif
                 endif
          continue
              enddo               !end loop on single clusters
  c            print*,'## distmin ', distmin,' clinc ',clinc
-             if(distmin.le.clinc)then
+ c     QUIQUI------------
-                CLS_STORE(nplanes-ip+1,ibest)=iclm !<<<<
+ c     anche qui: non ci vuole clinc???
- *              ----------------------------
+             if(iclm.ne.0)then
- c               print*,'~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~'
                 if(mod(VIEW(iclm),2).eq.0)then
-                   XGOOD(nplanes-ip+1)=1.
+                   clincnew=
-                   resx(nplanes-ip+1)=rxmm
+      $                 20*
-                   if(DEBUG)print*,'%%%% included X-cl ',iclm
+      $                 sqrt(rxmm**2+RCHI2_STORE(ibest)*k(ip)*cov(1,1))
- c                  if(.true.)print*,'%%%% included X-cl ',iclm
+                else if(mod(VIEW(iclm),2).ne.0)then
-      $                 ,'( chi^2, ',RCHI2_STORE(ibest)
+                   clincnew=
-      $                 ,', norm.dist.= ',distmin
+      $                 10*
-      $                 ,', cut ',clinc,' )'
+      $                 sqrt(rymm**2+RCHI2_STORE(ibest)*k(ip)*cov(2,2))
-                else
+                endif
-                   YGOOD(nplanes-ip+1)=1.
+ c     QUIQUI------------
-                   resy(nplanes-ip+1)=rymm
-                   if(DEBUG)print*,'%%%% included Y-cl ',iclm
+                if(distmin.le.clincnew)then   !QUIQUI
- c                  if(.true.)print*,'%%%% included Y-cl ',iclm
+ c     if(distmin.le.clinc)then          !QUIQUI
-      $                 ,'( chi^2, ',RCHI2_STORE(ibest)
-      $                 ,', norm.dist.= ', distmin
+                   CLS_STORE(nplanes-ip+1,ibest)=iclm !<<<<
-      $                 ,', cut ',clinc,' )'
+ *     ----------------------------
+ c     print*,'~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~'
+                   if(mod(VIEW(iclm),2).eq.0)then
+                      XGOOD(nplanes-ip+1)=1.
+                      resx(nplanes-ip+1)=rxmm
+                      if(DEBUG.EQ.1)print*,'%%%% included X-cl ',iclm
+      $                    ,'( chi^2, ',RCHI2_STORE(ibest)
+      $                    ,', dist.= ',distmin
+      $                    ,', cut ',clinc,' )'
+                   else
+                      YGOOD(nplanes-ip+1)=1.
+                      resy(nplanes-ip+1)=rymm
+                      if(DEBUG.EQ.1)print*,'%%%% included Y-cl ',iclm
+      $                    ,'( chi^2, ',RCHI2_STORE(ibest)
+      $                    ,', dist.= ', distmin
+      $                    ,', cut ',clinc,' )'
+                   endif
+ c     print*,'~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~'
+ *     ----------------------------
+                   xm_A(nplanes-ip+1) = xmm_A
+                   ym_A(nplanes-ip+1) = ymm_A
+                   xm_B(nplanes-ip+1) = xmm_B
+                   ym_B(nplanes-ip+1) = ymm_B
+                   zm(nplanes-ip+1) = (zmm_A+zmm_B)/2.
+                   dedxtrk_x(nplanes-ip+1) = dedxmmx !<<<
+                   dedxtrk_y(nplanes-ip+1) = dedxmmy !<<<
+ *     ----------------------------
                 endif
- c               print*,'~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~'
- *              ----------------------------
-                xm_A(nplanes-ip+1) = xmm_A
-                ym_A(nplanes-ip+1) = ymm_A
-                xm_B(nplanes-ip+1) = xmm_B
-                ym_B(nplanes-ip+1) = ymm_B
-                zm(nplanes-ip+1) = (zmm_A+zmm_B)/2.
- c              dedxtrk(nplanes-ip+1) = dedxmm !<<<    !(1)
-                dedxtrk_x(nplanes-ip+1) = dedxmmx !<<< !(1)
-                dedxtrk_y(nplanes-ip+1) = dedxmmy !<<< !(1)
- *              ----------------------------
              endif
           endif
     continue
-Line 3050 
 c              dedxtrk(nplanes-ip+1) = d
+Line 3405 
 c              dedxtrk(nplanes-ip+1) = d
  *                                                 *
  *                                                 *
  **************************************************
- cccccc 12/08/2006 modified by elena ---> (1)
  *
        subroutine clean_XYclouds(ibest,iflag)
        include 'commontracker.f'
        include 'level1.f'
        include 'common_momanhough.f'
- c      include 'momanhough_init.f'
+       include 'level2.f'
-       include 'level2.f'        !(1)
- c      include 'calib.f'
- c      include 'level1.f'
+       if(DEBUG.EQ.1)print*,'clean_XYclouds:'
        do ip=1,nplanes           !loop on planes
-Line 3072 
 c      include 'level1.f'
+Line 3423 
 c      include 'level1.f'
              if(id.ne.0)then
                 iclx=clx(ip,icp_cp(id))
                 icly=cly(ip,icp_cp(id))
- c               cl_used(iclx)=1  !tag used clusters
+ c$$$               cl_used(iclx)=ntrk  !tag used clusters
- c               cl_used(icly)=1  !tag used clusters
+ c$$$               cl_used(icly)=ntrk  !tag used clusters
-                cl_used(iclx)=ntrk  !tag used clusters !(1)
-                cl_used(icly)=ntrk  !tag used clusters !(1)
              elseif(icl.ne.0)then
- c               cl_used(icl)=1   !tag used clusters
+ c$$$               cl_used(icl)=ntrk   !tag used clusters
-                cl_used(icl)=ntrk   !tag used clusters !1)
              endif
- c               if(DEBUG)then
- c                  print*,ip,' <<< ',id
- c               endif
  *     -----------------------------
  *     remove the couple from clouds
  *     remove also vitual couples containing the
-Line 3100 
 c               endif
+Line 3445 
 c               endif
       $              cly(ip,icp).eq.icl
       $              )then
                    id=id_cp(ip,icp,1)
-                   if(DEBUG)then
+                   if(DEBUG.EQ.1)then
                       print*,ip,' <<< cp ',id
       $                    ,' ( cl-x '
       $                    ,clx(ip,icp)
-Line 3144 
 c               endif
+Line 3489 
 c               endif
        include 'level1.f'
        include 'common_momanhough.f'
        include 'level2.f'
- c      include 'level1.f'
+ *     ---------------------------------
+ *     variables initialized from level1
+ *     ---------------------------------
        do i=1,nviews
           good2(i)=good1(i)
+          do j=1,nva1_view
+             vkflag(i,j)=1
+             if(cnnev(i,j).le.0)then
+                vkflag(i,j)=cnnev(i,j)
+             endif
+          enddo
        enddo
+ *     ----------------
+ *     level2 variables
+ *     ----------------
        NTRK = 0
        do it=1,NTRKMAX
           IMAGE(IT)=0
-Line 3161 
 c      include 'level1.f'
+Line 3515 
 c      include 'level1.f'
              ZM_nt(IP,IT) = 0
              RESX_nt(IP,IT) = 0
              RESY_nt(IP,IT) = 0
+             TAILX_nt(IP,IT) = 0
+             TAILY_nt(IP,IT) = 0
+             XBAD(IP,IT) = 0
+             YBAD(IP,IT) = 0
              XGOOD_nt(IP,IT) = 0
              YGOOD_nt(IP,IT) = 0
+             LS(IP,IT) = 0
              DEDX_X(IP,IT) = 0
              DEDX_Y(IP,IT) = 0
              CLTRX(IP,IT) = 0
              CLTRY(IP,IT) = 0
+             multmaxx(ip,it) = 0
+             seedx(ip,it)    = 0
+             xpu(ip,it)      = 0
+             multmaxy(ip,it) = 0
+             seedy(ip,it)    = 0
+             ypu(ip,it)      = 0
           enddo
           do ipa=1,5
              AL_nt(IPA,IT) = 0
-Line 3295 
 c      include 'level1.f'
+Line 3660 
 c      include 'level1.f'
        include 'commontracker.f'
- c      include 'level1.f'
        include 'level1.f'
        include 'common_momanhough.f'
        include 'level2.f'
        include 'common_mini_2.f'
-       real sinth,phi,pig
+       include 'calib.f'
+       character*10 PFA
+       common/FINALPFA/PFA
+       real sinth,phi,pig
+       integer ssensor,sladder
        pig=acos(-1.)
+ *     -------------------------------------
        chi2_nt(ntr)        = sngl(chi2)
        nstep_nt(ntr)       = nstep
+ *     -------------------------------------
        phi   = al(4)
        sinth = al(3)
        if(sinth.lt.0)then
-Line 3318 
 c      include 'level1.f'
+Line 3689 
 c      include 'level1.f'
       $     phi = phi + 2*pig
        al(4) = phi
        al(3) = sinth
        do i=1,5
           al_nt(i,ntr)     = sngl(al(i))
           do j=1,5
              coval(i,j,ntr) = sngl(cov(i,j))
           enddo
        enddo
+ *     -------------------------------------
        do ip=1,nplanes           ! loop on planes
           xgood_nt(ip,ntr) = int(xgood(ip))
           ygood_nt(ip,ntr) = int(ygood(ip))
-Line 3334 
 c      include 'level1.f'
+Line 3704 
 c      include 'level1.f'
           zm_nt(ip,ntr)    = sngl(zm(ip))
           RESX_nt(IP,ntr)  = sngl(resx(ip))
           RESY_nt(IP,ntr)  = sngl(resy(ip))
+          TAILX_nt(IP,ntr) = 0.
+          TAILY_nt(IP,ntr) = 0.
           xv_nt(ip,ntr)    = sngl(xv(ip))
           yv_nt(ip,ntr)    = sngl(yv(ip))
           zv_nt(ip,ntr)    = sngl(zv(ip))
           axv_nt(ip,ntr)   = sngl(axv(ip))
           ayv_nt(ip,ntr)   = sngl(ayv(ip))
- c     l'avevo dimenticato!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
           factor = sqrt(
-      $        sin( acos(-1.) * sngl(axv(ip)) /180. )**2 +
+      $        tan( acos(-1.) * sngl(axv(ip)) /180. )**2 +
-      $        sin( acos(-1.) * sngl(ayv(ip)) /180. )**2 +
+      $        tan( acos(-1.) * sngl(ayv(ip)) /180. )**2 +
       $        1. )
- c     !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
           dedx_x(ip,ntr)   = sngl(dedxtrk_x(ip)/factor)
           dedx_y(ip,ntr)   = sngl(dedxtrk_y(ip)/factor)
-          id  = CP_STORE(ip,IDCAND)
+          ax   = axv_nt(ip,ntr)
+          ay   = ayv_nt(ip,ntr)
+          bfx  = BX_STORE(ip,IDCAND)
+          bfy  = BY_STORE(ip,IDCAND)
+ c$$$         if(ip.eq.6) ax = -1. * axv_nt(ip,ntr)
+ c$$$         if(ip.eq.6) bfy = -1. * BY_STORE(ip,IDCAND)
+ c$$$         tgtemp   = tan(ax*acos(-1.)/180.) + pmuH_h*bfy*0.00001
+ c$$$         angx     = 180.*atan(tgtemp)/acos(-1.)
+ c$$$         tgtemp = tan(ay*acos(-1.)/180.)+pmuH_e*bfx*0.00001
+ c$$$         angy    = 180.*atan(tgtemp)/acos(-1.)
+          angx = effectiveangle(ax,2*ip,bfy)
+          angy = effectiveangle(ay,2*ip-1,bfx)
+ c         print*,'* ',ip,bfx,bfy,angx,angy
+          id  = CP_STORE(ip,IDCAND) ! couple id
           icl = CLS_STORE(ip,IDCAND)
+          ssensor = -1
+          sladder = -1
+          ssensor = SENSOR_STORE(ip,IDCAND)
+          sladder = LADDER_STORE(ip,IDCAND)
+          if(ip.eq.6.and.ssensor.ne.0)ssensor = 3 - ssensor !notazione paolo x align
+          LS(IP,ntr)      = ssensor+10*sladder
           if(id.ne.0)then
+ c           >>> is a couple
              cltrx(ip,ntr)   = clx(nplanes-ip+1,icp_cp(id))
              cltry(ip,ntr)   = cly(nplanes-ip+1,icp_cp(id))
- c            print*,ip,' ',cltrx(ip,ntr),cltry(ip,ntr)
+             cl_used(cltrx(ip,ntr)) = 1     !tag used clusters
+             cl_used(cltry(ip,ntr)) = 1     !tag used clusters
+             xbad(ip,ntr)= nbadstrips(4,clx(nplanes-ip+1,icp_cp(id)))
+             ybad(ip,ntr)= nbadstrips(4,cly(nplanes-ip+1,icp_cp(id)))
+             if(nsatstrips(clx(nplanes-ip+1,icp_cp(id))).gt.0)
+      $           dedx_x(ip,ntr)=-dedx_x(ip,ntr)
+             if(nsatstrips(cly(nplanes-ip+1,icp_cp(id))).gt.0)
+      $           dedx_y(ip,ntr)=-dedx_y(ip,ntr)
+             multmaxx(ip,ntr) = maxs(cltrx(ip,ntr))
+      $                         +10000*mult(cltrx(ip,ntr))
+             seedx(ip,ntr)    = clsignal(indmax(cltrx(ip,ntr)))
+      $           /clsigma(indmax(cltrx(ip,ntr)))
+             call applypfa(PFA,cltrx(ip,ntr),angx,corr,res)
+             xpu(ip,ntr)      = corr
+             multmaxy(ip,ntr) = maxs(cltry(ip,ntr))
+      $                         +10000*mult(cltry(ip,ntr))
+             seedy(ip,ntr)    = clsignal(indmax(cltry(ip,ntr)))
+      $           /clsigma(indmax(cltry(ip,ntr)))
+             call applypfa(PFA,cltry(ip,ntr),angy,corr,res)
+             ypu(ip,ntr)      = corr
           elseif(icl.ne.0)then
-             if(mod(VIEW(icl),2).eq.0)cltrx(ip,ntr)=icl
-             if(mod(VIEW(icl),2).eq.1)cltry(ip,ntr)=icl
+             cl_used(icl) = 1    !tag used clusters
- c            print*,ip,' ',cltrx(ip,ntr),cltry(ip,ntr)
+             if(mod(VIEW(icl),2).eq.0)then
+                cltrx(ip,ntr)=icl
+                xbad(ip,ntr) = nbadstrips(4,icl)
+                if(nsatstrips(icl).gt.0)dedx_x(ip,ntr)=-dedx_x(ip,ntr)
+                multmaxx(ip,ntr) = maxs(cltrx(ip,ntr))
+      $                         +10000*mult(cltrx(ip,ntr))
+                seedx(ip,ntr)    = clsignal(indmax(cltrx(ip,ntr)))
+      $           /clsigma(indmax(cltrx(ip,ntr)))
+                call applypfa(PFA,cltrx(ip,ntr),angx,corr,res)
+                xpu(ip,ntr)      = corr
+             elseif(mod(VIEW(icl),2).eq.1)then
+                cltry(ip,ntr)=icl
+                ybad(ip,ntr) = nbadstrips(4,icl)
+                if(nsatstrips(icl).gt.0)dedx_y(ip,ntr)=-dedx_y(ip,ntr)
+                multmaxy(ip,ntr) = maxs(cltry(ip,ntr))
+      $                         +10000*mult(cltry(ip,ntr))
+                seedy(ip,ntr)    = clsignal(indmax(cltry(ip,ntr)))
+      $           /clsigma(indmax(cltry(ip,ntr)))
+                call applypfa(PFA,cltry(ip,ntr),angy,corr,res)
+                ypu(ip,ntr)      = corr
+             endif
           endif
        enddo
+       if(DEBUG.eq.1)then
+          print*,'> STORING TRACK ',ntr
+          print*,'clusters: '
+          do ip=1,6
+             print*,'> ',ip,' -- ',cltrx(ip,ntr),cltry(ip,ntr)
+          enddo
+       endif
+ c$$$      print*,(xgood(i),i=1,6)
+ c$$$      print*,(ygood(i),i=1,6)
+ c$$$      print*,(ls(i,ntr),i=1,6)
+ c$$$      print*,(dedx_x(i,ntr),i=1,6)
+ c$$$      print*,(dedx_y(i,ntr),i=1,6)
+ c$$$      print*,'-----------------------'
        end
-Line 3374 
 c            print*,ip,' ',cltrx(ip,ntr)
+Line 3841 
 c            print*,ip,' ',cltrx(ip,ntr)
  *     -------------------------------------------------------
        include 'commontracker.f'
- c      include 'level1.f'
        include 'calib.f'
        include 'level1.f'
        include 'common_momanhough.f'
-Line 3382 
 c      include 'level1.f'
+Line 3848 
 c      include 'level1.f'
        include 'common_xyzPAM.f'
  *     count #cluster per plane not associated to any track
- c      good2=1!.true.
        nclsx = 0
        nclsy = 0
        do iv = 1,nviews
-          if( mask_view(iv).ne.0 )good2(iv) = 20+mask_view(iv)
+ c         if( mask_view(iv).ne.0 )good2(iv) = 20+mask_view(iv)
+          good2(iv) = good2(iv) + mask_view(iv)*2**8
        enddo
+       if(DEBUG.eq.1)then
+          print*,'> STORING SINGLETS '
+       endif
        do icl=1,nclstr1
+          ip=nplanes-npl(VIEW(icl))+1
           if(cl_used(icl).eq.0)then !cluster not included in any track
-             ip=nplanes-npl(VIEW(icl))+1
              if(mod(VIEW(icl),2).eq.0)then !=== X views
                 nclsx = nclsx + 1
                 planex(nclsx) = ip
-                sgnlxs(nclsx) = sgnl(icl)/mip(VIEW(icl),LADDER(icl))!(2)
+                sgnlxs(nclsx) = sgnl(icl)/mip(VIEW(icl),LADDER(icl))
+                if(nsatstrips(icl).gt.0)sgnlxs(nclsx)=-sgnlxs(nclsx)
                 clsx(nclsx)   = icl
                 do is=1,2
  c                  call xyz_PAM(icl,0,is,'COG1',' ',0.,0.)
-Line 3412 
 c$$$               print*,'xs(2,nclsx)
+Line 3885 
 c$$$               print*,'xs(2,nclsx)
              else                          !=== Y views
                 nclsy = nclsy + 1
                 planey(nclsy) = ip
-                sgnlys(nclsy) = sgnl(icl)/mip(VIEW(icl),LADDER(icl))!(2)
+                sgnlys(nclsy) = sgnl(icl)/mip(VIEW(icl),LADDER(icl))
+                if(nsatstrips(icl).gt.0)sgnlys(nclsy)=-sgnlys(nclsy)
                 clsy(nclsy)   = icl
                 do is=1,2
  c                  call xyz_PAM(0,icl,is,' ','COG1',0.,0.)
-Line 3427 
 c$$$               print*,'ys(1,nclsy)
+Line 3901 
 c$$$               print*,'ys(1,nclsy)
  c$$$               print*,'ys(2,nclsy)   ',ys(2,nclsy)
              endif
           endif
- c      print*,icl,cl_used(icl),cl_good(icl),ip,VIEW(icl)!nclsx(ip),nclsy(ip)
  ***** LO METTO QUI PERCHE` NON SO DOVE METTERLO
           whichtrack(icl) = cl_used(icl)
+ *     --------------------------------------------------
+ *     per non perdere la testa...
+ *     whichtrack(icl) e` una variabile del common level1
+ *     che serve solo per sapere quali cluster sono stati
+ *     associati ad una traccia, e permettere di salvare
+ *     solo questi nell'albero di uscita
+ *     --------------------------------------------------
+ c$$$         print*,' cl ',icl,' --> ',cl_used(icl)
+ c$$$
+ c$$$         if( cl_used(icl).ne.0 )then
+ c$$$            if(
+ c$$$     $           mod(VIEW(icl),2).eq.0.and.
+ c$$$     $           cltrx(ip,whichtrack(icl)).ne.icl )
+ c$$$     $           print*,'**WARNING** cltrx(',ip,',',whichtrack(icl)
+ c$$$     $           ,')=',cltrx(ip,whichtrack(icl)),'.ne.',icl
+ c$$$            if(
+ c$$$     $           mod(VIEW(icl),2).eq.1.and.
+ c$$$     $           cltry(ip,whichtrack(icl)).ne.icl )
+ c$$$     $           print*,'**WARNING** cltry(',ip,',',whichtrack(icl)
+ c$$$     $           ,')=',cltry(ip,whichtrack(icl)),'.ne.',icl
+ c$$$         endif
        enddo
        end

 Legend:



Removed from v.1.18
 


changed lines


 
Added in v.1.31
 Legend:



Removed from v.1.18
 


changed lines


 
Added in v.1.31
-Removed from v.1.18
+Added in v.1.31

	ViewVC Help
Powered by ViewVC 1.1.23