/[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.32 by bongi, Tue Sep 4 15:44:49 2007 UTC revision 1.34 by pam-fi, Wed Mar 5 17:00:20 2008 UTC
# Line 713  c$$$      print*,icx,icy,sensor,PFAx,PFA Line 713  c$$$      print*,icx,icy,sensor,PFAx,PFA
713        xPAM_B = 0.D0        xPAM_B = 0.D0
714        yPAM_B = 0.D0        yPAM_B = 0.D0
715        zPAM_B = 0.D0        zPAM_B = 0.D0
716  c      print*,'## xyz_PAM: ',icx,icy,sensor,PFAx,PFAy,angx,angy  cc      print*,'## xyz_PAM: ',icx,icy,sensor,PFAx,PFAy,angx,angy
717    
718        if(sensor.lt.1.or.sensor.gt.2)then        if(sensor.lt.1.or.sensor.gt.2)then
719           print*,'xyz_PAM   ***ERROR*** wrong input '           print*,'xyz_PAM   ***ERROR*** wrong input '
# Line 805  c         resyPAM = RESYAV Line 805  c         resyPAM = RESYAV
805    
806   20   endif   20   endif
807    
808  c$$$      print*,'## stripx,stripy ',stripx,stripy  cc      print*,'## stripx,stripy ',stripx,stripy
809    
810  c===========================================================  c===========================================================
811  C     COUPLE  C     COUPLE
# Line 1053  c$$$      PFAy = 'COG4'!PFA Line 1053  c$$$      PFAy = 'COG4'!PFA
1053    
1054        if(icx.gt.nclstr1.or.icy.gt.nclstr1)then        if(icx.gt.nclstr1.or.icy.gt.nclstr1)then
1055              print*,'xyzpam: ***WARNING*** clusters ',icx,icy              print*,'xyzpam: ***WARNING*** clusters ',icx,icy
1056       $           ,' does not exists (nclstr1=',nclstr1,')'       $           ,' do not exists (n.clusters=',nclstr1,')'
1057              icx = -1*icx              icx = -1*icx
1058              icy = -1*icy              icy = -1*icy
1059              return              return
# Line 3052  c$$$      endif Line 3052  c$$$      endif
3052        call track_init        call track_init
3053        do ip=1,nplanes           !loop on planes        do ip=1,nplanes           !loop on planes
3054    
3055             if(DEBUG.EQ.1)print*,' ........... plane ',ip,' ........... '
3056    
3057           xP=XV_STORE(nplanes-ip+1,ibest)           xP=XV_STORE(nplanes-ip+1,ibest)
3058           yP=YV_STORE(nplanes-ip+1,ibest)           yP=YV_STORE(nplanes-ip+1,ibest)
3059           zP=ZV_STORE(nplanes-ip+1,ibest)           zP=ZV_STORE(nplanes-ip+1,ibest)
# Line 3153  c            if(DEBUG.EQ.1)print*,'>>>> Line 3155  c            if(DEBUG.EQ.1)print*,'>>>>
3155                 if(LADDER(icx).ne.nldt.or. !If the ladder number does not match                 if(LADDER(icx).ne.nldt.or. !If the ladder number does not match
3156  c     $              cl_used(icx).eq.1.or. !or the X cluster is already used  c     $              cl_used(icx).eq.1.or. !or the X cluster is already used
3157  c     $              cl_used(icy).eq.1.or. !or the Y cluster is already used  c     $              cl_used(icy).eq.1.or. !or the Y cluster is already used
3158       $              cl_used(icx).ne.0.or. !or the X cluster is already used !(3)       $              cl_used(icx).ne.0.or. !or the X cluster is already used
3159       $              cl_used(icy).ne.0.or. !or the Y cluster is already used !(3)       $              cl_used(icy).ne.0.or. !or the Y cluster is already used
3160       $              .false.)goto 1188 !then jump to next couple.       $              .false.)goto 1188 !then jump to next couple.
3161  *            *          
3162                 call xyz_PAM(icx,icy,ist,                 call xyz_PAM(icx,icy,ist,
# Line 3300  c                 dedxmm = sgnl(icy)  !( Line 3302  c                 dedxmm = sgnl(icy)  !(
3302  *----- single clusters -----------------------------------------------    *----- single clusters -----------------------------------------------  
3303  c            print*,'## ncls(',ip,') ',ncls(ip)  c            print*,'## ncls(',ip,') ',ncls(ip)
3304              do ic=1,ncls(ip)    !loop on single clusters              do ic=1,ncls(ip)    !loop on single clusters
3305    c               print*,'-',ic,'-'
3306                 icl=cls(ip,ic)                 icl=cls(ip,ic)
3307  c               if(cl_used(icl).eq.1.or.     !if the cluster is already used  c               if(cl_used(icl).eq.1.or.     !if the cluster is already used
3308                 if(cl_used(icl).ne.0.or.     !if the cluster is already used !(3)                 if(cl_used(icl).ne.0.or.     !if the cluster is already used !(3)
# Line 3559  c$$$               cl_used(icl)=ntrk   ! Line 3562  c$$$               cl_used(icl)=ntrk   !
3562          ys(1,ip)=0          ys(1,ip)=0
3563          ys(2,ip)=0          ys(2,ip)=0
3564          sgnlys(ip)=0          sgnlys(ip)=0
3565            sxbad(ip)=0
3566            sybad(ip)=0
3567            multmaxsx(ip)=0
3568            multmaxsy(ip)=0
3569        enddo        enddo
3570        end        end
3571    
# Line 3787  c           >>> is a couple Line 3794  c           >>> is a couple
3794    
3795           elseif(icl.ne.0)then           elseif(icl.ne.0)then
3796    
3797    
3798              cl_used(icl) = 1    !tag used clusters                        cl_used(icl) = 1    !tag used clusters          
3799    
3800              if(mod(VIEW(icl),2).eq.0)then              if(mod(VIEW(icl),2).eq.0)then
3801                 cltrx(ip,ntr)=icl                 cltrx(ip,ntr)=icl
   
3802                 xbad(ip,ntr) = nbadstrips(4,icl)                 xbad(ip,ntr) = nbadstrips(4,icl)
3803    
3804                 if(nsatstrips(icl).gt.0)dedx_x(ip,ntr)=-dedx_x(ip,ntr)                 if(nsatstrips(icl).gt.0)dedx_x(ip,ntr)=-dedx_x(ip,ntr)
# Line 3805  c           >>> is a couple Line 3812  c           >>> is a couple
3812    
3813              elseif(mod(VIEW(icl),2).eq.1)then              elseif(mod(VIEW(icl),2).eq.1)then
3814                 cltry(ip,ntr)=icl                 cltry(ip,ntr)=icl
   
3815                 ybad(ip,ntr) = nbadstrips(4,icl)                 ybad(ip,ntr) = nbadstrips(4,icl)
3816    
3817                 if(nsatstrips(icl).gt.0)dedx_y(ip,ntr)=-dedx_y(ip,ntr)                 if(nsatstrips(icl).gt.0)dedx_y(ip,ntr)=-dedx_y(ip,ntr)
# Line 3873  c         if( mask_view(iv).ne.0 )good2( Line 3879  c         if( mask_view(iv).ne.0 )good2(
3879           ip=nplanes-npl(VIEW(icl))+1                       ip=nplanes-npl(VIEW(icl))+1            
3880                    
3881           if(cl_used(icl).eq.0)then !cluster not included in any track           if(cl_used(icl).eq.0)then !cluster not included in any track
3882    
3883              if(mod(VIEW(icl),2).eq.0)then !=== X views              if(mod(VIEW(icl),2).eq.0)then !=== X views
3884    
3885                 nclsx = nclsx + 1                 nclsx = nclsx + 1
3886                 planex(nclsx) = ip                 planex(nclsx) = ip
3887                 sgnlxs(nclsx) = sgnl(icl)/mip(VIEW(icl),LADDER(icl))                 sgnlxs(nclsx) = sgnl(icl)/mip(VIEW(icl),LADDER(icl))
3888                 if(nsatstrips(icl).gt.0)sgnlxs(nclsx)=-sgnlxs(nclsx)                 if(nsatstrips(icl).gt.0)sgnlxs(nclsx)=-sgnlxs(nclsx)
3889                 clsx(nclsx)   = icl                 clsx(nclsx)   = icl
3890                   sxbad(nclsx)  = nbadstrips(1,icl)
3891                   multmaxsx(nclsx) = maxs(icl)+10000*mult(icl)
3892                  
3893    cc               print*,icl,' >>>> ',sxbad(nclsx)
3894    
3895                 do is=1,2                 do is=1,2
3896  c                  call xyz_PAM(icl,0,is,'COG1',' ',0.,0.)  c                  call xyz_PAM(icl,0,is,'COG1',' ',0.,0.)
3897  c                  call xyz_PAM(icl,0,is,PFAdef,' ',0.,0.)  c                  call xyz_PAM(icl,0,is,PFAdef,' ',0.,0.)
# Line 3896  c$$$               print*,'xs(2,nclsx)   Line 3909  c$$$               print*,'xs(2,nclsx)  
3909                 sgnlys(nclsy) = sgnl(icl)/mip(VIEW(icl),LADDER(icl))                 sgnlys(nclsy) = sgnl(icl)/mip(VIEW(icl),LADDER(icl))
3910                 if(nsatstrips(icl).gt.0)sgnlys(nclsy)=-sgnlys(nclsy)                 if(nsatstrips(icl).gt.0)sgnlys(nclsy)=-sgnlys(nclsy)
3911                 clsy(nclsy)   = icl                 clsy(nclsy)   = icl
3912                   sybad(nclsy)  = nbadstrips(1,icl)
3913                   multmaxsy(nclsy) = maxs(icl)+10000*mult(icl)
3914    
3915    cc               print*,icl,' >>>> ',sybad(nclsy)
3916    
3917                 do is=1,2                 do is=1,2
3918  c                  call xyz_PAM(0,icl,is,' ','COG1',0.,0.)  c                  call xyz_PAM(0,icl,is,' ','COG1',0.,0.)
3919  c                  call xyz_PAM(0,icl,is,' ',PFAdef,0.,0.)  c                  call xyz_PAM(0,icl,is,' ',PFAdef,0.,0.)
# Line 3919  c$$$               print*,'ys(2,nclsy)   Line 3937  c$$$               print*,'ys(2,nclsy)  
3937  *     associati ad una traccia, e permettere di salvare  *     associati ad una traccia, e permettere di salvare
3938  *     solo questi nell'albero di uscita  *     solo questi nell'albero di uscita
3939  *     --------------------------------------------------  *     --------------------------------------------------
3940                            
   
 c$$$         print*,' cl ',icl,' --> ',cl_used(icl)  
 c$$$  
 c$$$         if( cl_used(icl).ne.0 )then  
 c$$$            if(  
 c$$$     $           mod(VIEW(icl),2).eq.0.and.  
 c$$$     $           cltrx(ip,whichtrack(icl)).ne.icl )  
 c$$$     $           print*,'**WARNING** cltrx(',ip,',',whichtrack(icl)  
 c$$$     $           ,')=',cltrx(ip,whichtrack(icl)),'.ne.',icl  
 c$$$            if(  
 c$$$     $           mod(VIEW(icl),2).eq.1.and.  
 c$$$     $           cltry(ip,whichtrack(icl)).ne.icl )  
 c$$$     $           print*,'**WARNING** cltry(',ip,',',whichtrack(icl)  
 c$$$     $           ,')=',cltry(ip,whichtrack(icl)),'.ne.',icl  
 c$$$         endif  
           
   
3941        enddo        enddo
3942        end        end
3943    
Legend:
Removed from v.1.32  
changed lines
  Added in v.1.34
  ViewVC Help
Powered by ViewVC 1.1.23