--- DarthVader/TrackerLevel2/src/F77/functionspfa.f	2006/05/30 16:30:37	1.2
+++ DarthVader/TrackerLevel2/src/F77/functionspfa.f	2006/09/28 14:04:40	1.3
@@ -686,25 +686,34 @@
          if(mod(iv,2).eq.1)incut=incuty
          if(mod(iv,2).eq.0)incut=incutx
 
-         istart=INDSTART(IC)
-         istop=TOTCLLENGTH
+         istart = INDSTART(IC)
+         istop  = TOTCLLENGTH
          if(ic.lt.NCLSTR1)istop=INDSTART(IC+1)-1
-         COG=0  
-         mu=0
-         do i=istart,istop
-            ipos=i-INDMAX(ic)
-            ivk=nvk(MAXS(ic)+ipos)
-            is=nst(MAXS(ic)+ipos)
+         COG = 0  
+         mu  = 0
+         do i = istart,istop
+            ipos = i-INDMAX(ic)
+            ivk  = nvk(MAXS(ic)+ipos)
+            is   = nst(MAXS(ic)+ipos)
 *            print*,'******************',istart,istop,ipos
 *     $           ,MAXS(ic)+ipos,iv,ivk,is
-            cut=incut*SIGMA(iv,ivk,is)
+            cut  = incut*SIGMA(iv,ivk,is)
+            if(SIGMA(iv,ivk,is).ne.CLSIGMA(i))
+     $           print*,'cog(0,ic) --> hai fatto qualche cazzata'
             if(CLSIGNAL(i).ge.cut)then
                COG = COG + ipos*CLSIGNAL(i)
                mu = mu + 1
 c               print*,ipos,CLSIGNAL(i),incut,cut
             endif
          enddo
-         COG=COG/DEDX(ic)
+         if(DEDX(ic).le.0)then
+            print*,'cog(0,ic) --> ic, dedx ',ic,DEDX(ic)
+            print*,(CLSIGNAL(i)/CLSIGMA(i),i=istart,istop)
+            print*,(CLSIGNAL(i),i=istart,istop)
+            print*,'cog(0,ic) --> NOT EVALUATED '
+         else
+            COG=COG/DEDX(ic)
+         endif
 c         if(DEBUG)print*,'COG  (ic ',ic,' m',mu,')'
 c     $        ,cog