/[PAMELA software]/DarthVader/TrackerLevel2/src/F77/functionspfa.f
ViewVC logotype

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

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.6 by pam-fi, Fri Dec 1 10:43:18 2006 UTC revision 1.7 by pam-fi, Thu Jan 11 10:20:58 2007 UTC
# Line 682  c         print*,'## ',sl2,sl1,sc,sr1,sr Line 682  c         print*,'## ',sl2,sl1,sc,sr1,sr
682              endif              endif
683           else           else
684              print*,'function COG(NCOG,IC) ==> WARNING!! NCOG=',NCOG              print*,'function COG(NCOG,IC) ==> WARNING!! NCOG=',NCOG
685              print*,'                 (NCOG must be 0-4)'       $           ,' not implemented'
686              COG = 0.              COG = 0.
687           endif           endif
688    
# Line 704  ' Line 704  '
704           mu  = 0           mu  = 0
705           do i = istart,istop           do i = istart,istop
706              ipos = i-INDMAX(ic)              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)  
707              cut  = incut*CLSIGMA(i)              cut  = incut*CLSIGMA(i)
 c            if(SIGMA(iv,ivk,is).ne.CLSIGMA(i))  
 c     $           print*,'cog(0,ic) --> hai fatto qualche cazzata'  
708              if(CLSIGNAL(i).ge.cut)then              if(CLSIGNAL(i).ge.cut)then
709                 COG = COG + ipos*CLSIGNAL(i)                 COG = COG + ipos*CLSIGNAL(i)
710                 mu = mu + 1                 mu = mu + 1
 c               print*,ipos,CLSIGNAL(i),incut,cut  
711              endif              endif
712           enddo           enddo
713           if(DEDX(ic).le.0)then           if(DEDX(ic).le.0)then
# Line 726  c               print*,ipos,CLSIGNAL(i), Line 718  c               print*,ipos,CLSIGNAL(i),
718           else           else
719              COG=COG/DEDX(ic)              COG=COG/DEDX(ic)
720           endif           endif
 c         if(DEBUG)print*,'COG  (ic ',ic,' m',mu,')'  
 c     $        ,cog  
721                    
722        else        else
723                    
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.23