/[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.25 by mocchiut, Tue Aug 4 14:01:37 2009 UTC revision 1.26 by mocchiut, Tue Aug 11 12:58:25 2009 UTC
# Line 1006  c      print *,'## ',sl2,sl1,sc,sr1,sr2 Line 1006  c      print *,'## ',sl2,sl1,sc,sr1,sr2
1006  c     ==============================================================  c     ==============================================================
1007           if(ncog.eq.1)then           if(ncog.eq.1)then
1008              COG = 0.              COG = 0.
1009              if(sr1.gt.sc)cog=1.           if(sr1.gt.sc)cog=1.
1010              if(sl1.gt.sc.and.sl1.gt.sr1)cog=-1.           if(sl1.gt.sc.and.sl1.gt.sr1)cog=-1.
1011  c     ==============================================================  c     ==============================================================
1012           elseif(ncog.eq.2)then           elseif(ncog.eq.2)then
1013              COG = 0.              COG = 0.
# Line 1015  c     ================================== Line 1015  c     ==================================
1015                 if((sl1+sc).ne.0)COG = -sl1/(sl1+sc)                         if((sl1+sc).ne.0)COG = -sl1/(sl1+sc)        
1016              elseif(sl1.lt.sr1)then              elseif(sl1.lt.sr1)then
1017                 if((sc+sr1).ne.0)COG = sr1/(sc+sr1)                 if((sc+sr1).ne.0)COG = sr1/(sc+sr1)
1018              elseif( sl1.eq.sr1.and.sl1.ne.-9999.)then           elseif( sl1.eq.sr1.and.sl1.ne.-9999.)then
1019                 if( clsigma(indmax(ic)-1).lt.clsigma(indmax(ic)+1)                 if( clsigma(indmax(ic)-1).lt.clsigma(indmax(ic)+1)
1020       $              .and.(sl1+sc).ne.0 )cog = -sl1/(sl1+sc)       $              .and.(sl1+sc).ne.0 )cog = -sl1/(sl1+sc)
1021                 if( clsigma(indmax(ic)-1).gt.clsigma(indmax(ic)+1)                 if( clsigma(indmax(ic)-1).gt.clsigma(indmax(ic)+1)
1022       $              .and.(sc+sr1).ne.0 )cog = sr1/(sc+sr1)       $              .and.(sc+sr1).ne.0 )cog = sr1/(sc+sr1)
1023              endif           endif
1024  c            if(cog==0)print*,'Strange cluster (2) - @maxs ',MAXS(ic)  c            if(cog==0)print*,'Strange cluster (2) - @maxs ',MAXS(ic)
1025  c     $           ,' : ',sl2,sl1,sc,sr1,sr2  c     $           ,' : ',sl2,sl1,sc,sr1,sr2
1026  c     ==============================================================  c     ==============================================================
Legend:
Removed from v.1.25  
changed lines
  Added in v.1.26
  ViewVC Help
Powered by ViewVC 1.1.23