/[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.2 by pam-fi, Tue May 30 16:30:37 2006 UTC revision 1.3 by pam-fi, Thu Sep 28 14:04:40 2006 UTC
# Line 686  ' Line 686  '
686           if(mod(iv,2).eq.1)incut=incuty           if(mod(iv,2).eq.1)incut=incuty
687           if(mod(iv,2).eq.0)incut=incutx           if(mod(iv,2).eq.0)incut=incutx
688    
689           istart=INDSTART(IC)           istart = INDSTART(IC)
690           istop=TOTCLLENGTH           istop  = TOTCLLENGTH
691           if(ic.lt.NCLSTR1)istop=INDSTART(IC+1)-1           if(ic.lt.NCLSTR1)istop=INDSTART(IC+1)-1
692           COG=0             COG = 0  
693           mu=0           mu  = 0
694           do i=istart,istop           do i = istart,istop
695              ipos=i-INDMAX(ic)              ipos = i-INDMAX(ic)
696              ivk=nvk(MAXS(ic)+ipos)              ivk  = nvk(MAXS(ic)+ipos)
697              is=nst(MAXS(ic)+ipos)              is   = nst(MAXS(ic)+ipos)
698  *            print*,'******************',istart,istop,ipos  *            print*,'******************',istart,istop,ipos
699  *     $           ,MAXS(ic)+ipos,iv,ivk,is  *     $           ,MAXS(ic)+ipos,iv,ivk,is
700              cut=incut*SIGMA(iv,ivk,is)              cut  = incut*SIGMA(iv,ivk,is)
701                if(SIGMA(iv,ivk,is).ne.CLSIGMA(i))
702         $           print*,'cog(0,ic) --> hai fatto qualche cazzata'
703              if(CLSIGNAL(i).ge.cut)then              if(CLSIGNAL(i).ge.cut)then
704                 COG = COG + ipos*CLSIGNAL(i)                 COG = COG + ipos*CLSIGNAL(i)
705                 mu = mu + 1                 mu = mu + 1
706  c               print*,ipos,CLSIGNAL(i),incut,cut  c               print*,ipos,CLSIGNAL(i),incut,cut
707              endif              endif
708           enddo           enddo
709           COG=COG/DEDX(ic)           if(DEDX(ic).le.0)then
710                print*,'cog(0,ic) --> ic, dedx ',ic,DEDX(ic)
711                print*,(CLSIGNAL(i)/CLSIGMA(i),i=istart,istop)
712                print*,(CLSIGNAL(i),i=istart,istop)
713                print*,'cog(0,ic) --> NOT EVALUATED '
714             else
715                COG=COG/DEDX(ic)
716             endif
717  c         if(DEBUG)print*,'COG  (ic ',ic,' m',mu,')'  c         if(DEBUG)print*,'COG  (ic ',ic,' m',mu,')'
718  c     $        ,cog  c     $        ,cog
719                    
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.23