--- DarthVader/TrackerLevel2/src/F77/functionspfa.f 2006/12/01 10:43:18 1.6 +++ DarthVader/TrackerLevel2/src/F77/functionspfa.f 2007/01/11 10:20:58 1.7 @@ -682,7 +682,7 @@ endif else print*,'function COG(NCOG,IC) ==> WARNING!! NCOG=',NCOG - print*,' (NCOG must be 0-4)' + $ ,' not implemented' COG = 0. endif @@ -704,18 +704,10 @@ mu = 0 do i = istart,istop ipos = i-INDMAX(ic) -cc ivk = nvk(MAXS(ic)+ipos) -cc is = nst(MAXS(ic)+ipos) -* print*,'******************',istart,istop,ipos -* $ ,MAXS(ic)+ipos,iv,ivk,is -cc cut = incut*SIGMA(iv,ivk,is) cut = incut*CLSIGMA(i) -c if(SIGMA(iv,ivk,is).ne.CLSIGMA(i)) -c $ 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 if(DEDX(ic).le.0)then @@ -726,8 +718,6 @@ else COG=COG/DEDX(ic) endif -c if(DEBUG)print*,'COG (ic ',ic,' m',mu,')' -c $ ,cog else