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

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

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

revision 1.26 by pam-fi, Thu May 24 16:45:48 2007 UTC revision 1.27 by pam-fi, Fri Aug 17 14:36:05 2007 UTC
# Line 698  c$$$         print*,fbad_cog(4,icx) Line 698  c$$$         print*,fbad_cog(4,icx)
698              stripx  = stripx + pfaeta3(icx,angx)                        stripx  = stripx + pfaeta3(icx,angx)          
699              resxPAM = risxeta3(abs(angx))                                    resxPAM = risxeta3(abs(angx))                      
700              resxPAM = resxPAM*fbad_cog(3,icx)                            resxPAM = resxPAM*fbad_cog(3,icx)              
701              if(DEBUG.and.fbad_cog(3,icx).ne.1)              c            if(DEBUG.and.fbad_cog(3,icx).ne.1)            
702       $           print*,'BAD icx >>> ',viewx,fbad_cog(3,icx)  c     $           print*,'BAD icx >>> ',viewx,fbad_cog(3,icx)
703    
704           elseif(PFAx.eq.'ETA4')then                                   elseif(PFAx.eq.'ETA4')then                        
705    
706              stripx  = stripx + pfaeta4(icx,angx)                          stripx  = stripx + pfaeta4(icx,angx)            
707              resxPAM = risxeta4(abs(angx))                                    resxPAM = risxeta4(abs(angx))                      
708              resxPAM = resxPAM*fbad_cog(4,icx)                            resxPAM = resxPAM*fbad_cog(4,icx)              
709              if(DEBUG.and.fbad_cog(4,icx).ne.1)                c            if(DEBUG.and.fbad_cog(4,icx).ne.1)              
710       $           print*,'BAD icx >>> ',viewx,fbad_cog(4,icx)  c     $           print*,'BAD icx >>> ',viewx,fbad_cog(4,icx)
711    
712           elseif(PFAx.eq.'ETA')then             elseif(PFAx.eq.'ETA')then  
713    
# Line 715  c$$$         print*,fbad_cog(4,icx) Line 715  c$$$         print*,fbad_cog(4,icx)
715  c            resxPAM = riseta(icx,angx)                      c            resxPAM = riseta(icx,angx)                    
716              resxPAM = riseta(viewx,angx)                                  resxPAM = riseta(viewx,angx)                    
717              resxPAM = resxPAM*fbad_eta(icx,angx)                          resxPAM = resxPAM*fbad_eta(icx,angx)            
718              if(DEBUG.and.fbad_cog(2,icx).ne.1)                c            if(DEBUG.and.fbad_cog(2,icx).ne.1)              
719       $           print*,'BAD icx >>> ',viewx,fbad_cog(2,icx)  c     $           print*,'BAD icx >>> ',viewx,fbad_cog(2,icx)
720    
721             elseif(PFAx.eq.'ETAL')then  
722    
723                stripx  = stripx + pfaetal(icx,angx)            
724                resxPAM = riseta(viewx,angx)                    
725                resxPAM = resxPAM*fbad_eta(icx,angx)            
726    c            if(DEBUG.and.fbad_cog(2,icx).ne.1)              
727    c     $           print*,'BAD icx >>> ',viewx,fbad_cog(2,icx)
728    
729           elseif(PFAx.eq.'COG')then                     elseif(PFAx.eq.'COG')then          
730    
# Line 824  c$$$         print*,fbad_cog(4,icy) Line 832  c$$$         print*,fbad_cog(4,icy)
832              stripy  = stripy + pfaeta2(icy,angy)                        stripy  = stripy + pfaeta2(icy,angy)          
833              resyPAM = risyeta2(abs(angy))                            resyPAM = risyeta2(abs(angy))              
834              resyPAM = resyPAM*fbad_cog(2,icy)              resyPAM = resyPAM*fbad_cog(2,icy)
835              if(DEBUG.and.fbad_cog(2,icy).ne.1)  c            if(DEBUG.and.fbad_cog(2,icy).ne.1)
836       $           print*,'BAD icy >>> ',viewy,fbad_cog(2,icy)  c     $           print*,'BAD icy >>> ',viewy,fbad_cog(2,icy)
837    
838           elseif(PFAy.eq.'ETA3')then                                 elseif(PFAy.eq.'ETA3')then                      
839    
840              stripy  = stripy + pfaeta3(icy,angy)              stripy  = stripy + pfaeta3(icy,angy)
841              resyPAM = resyPAM*fbad_cog(3,icy)                resyPAM = resyPAM*fbad_cog(3,icy)  
842              if(DEBUG.and.fbad_cog(3,icy).ne.1)  c            if(DEBUG.and.fbad_cog(3,icy).ne.1)
843       $           print*,'BAD icy >>> ',viewy,fbad_cog(3,icy)  c     $           print*,'BAD icy >>> ',viewy,fbad_cog(3,icy)
844    
845           elseif(PFAy.eq.'ETA4')then             elseif(PFAy.eq.'ETA4')then  
846    
847              stripy  = stripy + pfaeta4(icy,angy)              stripy  = stripy + pfaeta4(icy,angy)
848              resyPAM = resyPAM*fbad_cog(4,icy)              resyPAM = resyPAM*fbad_cog(4,icy)
849              if(DEBUG.and.fbad_cog(4,icy).ne.1)  c            if(DEBUG.and.fbad_cog(4,icy).ne.1)
850       $           print*,'BAD icy >>> ',viewy,fbad_cog(4,icy)  c     $           print*,'BAD icy >>> ',viewy,fbad_cog(4,icy)
851    
852           elseif(PFAy.eq.'ETA')then           elseif(PFAy.eq.'ETA')then
853    
# Line 847  c$$$         print*,fbad_cog(4,icy) Line 855  c$$$         print*,fbad_cog(4,icy)
855  c            resyPAM = riseta(icy,angy)    c            resyPAM = riseta(icy,angy)  
856              resyPAM = riseta(viewy,angy)                resyPAM = riseta(viewy,angy)  
857              resyPAM = resyPAM*fbad_eta(icy,angy)              resyPAM = resyPAM*fbad_eta(icy,angy)
858              if(DEBUG.and.fbad_cog(2,icy).ne.1)  c            if(DEBUG.and.fbad_cog(2,icy).ne.1)
859       $           print*,'BAD icy >>> ',viewy,fbad_cog(2,icy)  c     $           print*,'BAD icy >>> ',viewy,fbad_cog(2,icy)
860    
861             elseif(PFAy.eq.'ETAL')then
862    
863                stripy  = stripy + pfaetal(icy,angy)
864                resyPAM = riseta(viewy,angy)  
865                resyPAM = resyPAM*fbad_eta(icy,angy)
866    c            if(DEBUG.and.fbad_cog(2,icy).ne.1)
867    c     $           print*,'BAD icy >>> ',viewy,fbad_cog(2,icy)
868    
869           elseif(PFAy.eq.'COG')then           elseif(PFAy.eq.'COG')then
870    
Legend:
Removed from v.1.26  
changed lines
  Added in v.1.27
  ViewVC Help
Powered by ViewVC 1.1.23