/[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.14 by pam-fi, Tue Aug 7 13:56:29 2007 UTC revision 1.15 by pam-fi, Fri Aug 17 13:28:02 2007 UTC
# Line 210  c      if(mod(int(VIEW(ic)),2).eq.1)then Line 210  c      if(mod(int(VIEW(ic)),2).eq.1)then
210                            
211        endif        endif
212    
 cc      print*,'---- ',riseta,iview,angle  
213    
214   100  return   100  return
215        end        end
# Line 787  c     ================================== Line 786  c     ==================================
786                 if( clsigma(indmax(ic)-1).gt.clsigma(indmax(ic)+1)                 if( clsigma(indmax(ic)-1).gt.clsigma(indmax(ic)+1)
787       $              .and.(sc+sr1).ne.0 )cog = sr1/(sc+sr1) !NEW       $              .and.(sc+sr1).ne.0 )cog = sr1/(sc+sr1) !NEW
788              endif              endif
789  c$$$            if(cog==0)print*,'Strange cluster (2) - @maxs ',MAXS(ic)  c            if(cog==0)print*,'Strange cluster (2) - @maxs ',MAXS(ic)
790  c$$$     $           ,VIEW(ic),LADDER(ic)  c     $           ,' : ',sl2,sl1,sc,sr1,sr2
 c$$$     $           ,' : ',sl2,sl1,sc,sr1,sr2  
791  c     ==============================================================  c     ==============================================================
792           elseif(ncog.eq.3)then           elseif(ncog.eq.3)then
793               if( (sl1+sc+sr1).ne.0 )COG = (sr1-sl1)/(sl1+sc+sr1)               if( (sl1+sc+sr1).ne.0 )COG = (sr1-sl1)/(sl1+sc+sr1)
794  c$$$             if(cog==0)print*,'Strange cluster (3) - @maxs ',MAXS(ic)  c             if(cog==0)print*,'Strange cluster (3) - @maxs ',MAXS(ic)
795  c$$$     $           ,VIEW(ic),LADDER(ic)  c     $            ,' : ',sl2,sl1,sc,sr1,sr2
 c$$$     $            ,' : ',sl2,sl1,sc,sr1,sr2  
796  c     ==============================================================  c     ==============================================================
797           elseif(ncog.eq.4)then           elseif(ncog.eq.4)then
798              if(sl2.gt.sr2)then              if(sl2.gt.sr2)then
# Line 812  c     ================================== Line 809  c     ==================================
809       $              .and.(sr2+sl1+sc+sr1).ne.0 )       $              .and.(sr2+sl1+sc+sr1).ne.0 )
810       $              cog = (2*sr2+sr1-sl1)/(sr2+sl1+sc+sr1)  !NEW                     $              cog = (2*sr2+sr1-sl1)/(sr2+sl1+sc+sr1)  !NEW              
811              endif              endif
 c$$$            if(cog==0)print*,'Strange cluster (4) - @maxs ',MAXS(ic)  
 c$$$     $           ,VIEW(ic),LADDER(ic)  
 c$$$     $           ,' : ',sl2,sl1,sc,sr1,sr2  
 c     ==============================================================  
812           else           else
813              print*,'function COG(NCOG,IC) ==> WARNING!! NCOG=',NCOG              print*,'function COG(NCOG,IC) ==> WARNING!! NCOG=',NCOG
814       $           ,' not implemented'       $           ,' not implemented'
Legend:
Removed from v.1.14  
changed lines
  Added in v.1.15
  ViewVC Help
Powered by ViewVC 1.1.23