--- DarthVader/TrackerLevel2/src/F77/functionspfa.f 2007/05/15 16:22:18 1.11 +++ DarthVader/TrackerLevel2/src/F77/functionspfa.f 2007/05/23 15:31:19 1.12 @@ -785,8 +785,8 @@ if((sl2+sl1+sc+sr1).ne.0) $ COG = (sr1-sl1-2*sl2)/(sl2+sl1+sc+sr1) elseif(sl2.le.sr2)then - if((sl2+sl1+sc+sr1).ne.0) - $ COG = (2*sr2+sr1-sl1)/(sl2+sl1+sc+sr1) + if((sr2+sl1+sc+sr1).ne.0) + $ COG = (2*sr2+sr1-sl1)/(sr2+sl1+sc+sr1) endif else print*,'function COG(NCOG,IC) ==> WARNING!! NCOG=',NCOG @@ -818,7 +818,7 @@ COG = COG + ipos*CLSIGNAL(i) SGN = SGN + CLSIGNAL(i) mu = mu + 1 - print*,ipos,CLSIGNAL(i) +c print*,ipos,CLSIGNAL(i) else goto 10 endif @@ -831,14 +831,14 @@ COG = COG + ipos*CLSIGNAL(i) SGN = SGN + CLSIGNAL(i) mu = mu + 1 - print*,ipos,CLSIGNAL(i) +c print*,ipos,CLSIGNAL(i) else goto 20 endif enddo 20 continue if(SGN.le.0)then -c print*,'cog(0,ic) --> ic, dedx ',ic,SGN + print*,'cog(0,ic) --> ic, dedx ',ic,SGN print*,(CLSIGNAL(i)/CLSIGMA(i),i=istart,istop) print*,(CLSIGNAL(i),i=istart,istop) c print*,'cog(0,ic) --> NOT EVALUATED '