/[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.40 by mocchiut, Tue Aug 4 14:01:35 2009 UTC revision 1.42 by pam-fi, Tue Nov 10 15:50:16 2009 UTC
# Line 303  c$$$         if(ibest.eq.0)goto 880 !>> Line 303  c$$$         if(ibest.eq.0)goto 880 !>>
303  *     1st) decreasing n.points  *     1st) decreasing n.points
304  *     2nd) increasing chi**2  *     2nd) increasing chi**2
305  *     -------------------------------------------------------  *     -------------------------------------------------------
306           rchi2best=1000000000.           rchi2best=1000000000.
307           ndofbest=0                       ndofbest=0            
308           do i=1,ntracks           do i=1,ntracks
309             ndof=0                           ndof=0              
# Line 706  c         resxPAM = RESXAV Line 706  c         resxPAM = RESXAV
706           stripx  = stripx + corr           stripx  = stripx + corr
707           resxPAM = res           resxPAM = res
708    
709   10   endif   10   continue
710          endif
711            
712  *     -----------------  *     -----------------
713  *     CLUSTER Y  *     CLUSTER Y
# Line 752  c         resyPAM = RESYAV Line 753  c         resyPAM = RESYAV
753           stripy  = stripy + corr           stripy  = stripy + corr
754           resyPAM = res           resyPAM = res
755    
756   20   endif   20   continue
757          endif
758    
759    
760  c===========================================================  c===========================================================
# Line 3882  c$$$         angy    = 180.*atan(tgtemp) Line 3884  c$$$         angy    = 180.*atan(tgtemp)
3884           if(ip.eq.6.and.ssensor.ne.0)ssensor = 3 - ssensor !notazione paolo x align           if(ip.eq.6.and.ssensor.ne.0)ssensor = 3 - ssensor !notazione paolo x align
3885           LS(IP,ntr)      = ssensor+10*sladder           LS(IP,ntr)      = ssensor+10*sladder
3886    
3887           if(id.ne.0)then  c         if(id.ne.0)then
3888    CCCCCC(10 novembre 2009) PATCH X NUCLEI
3889    C     non ho capito perche', ma durante il ritracciamento dei nuclei
3890    C     (quando una traccia viene trovata ma non e' buona) c'e' qualche variabile
3891    C     che non viene reinizializzata correttamente e i cluster esclusi
3892    C     dal fit risultano ancora inclusi...
3893    C
3894             cltrx(ip,ntr)   = 0
3895             cltry(ip,ntr)   = 0
3896             if(
3897         $        xgood_nt(ip,ntr).eq.1.and.ygood_nt(ip,ntr).eq.1
3898         $        .and.
3899         $        id.ne.0)then
3900    
3901  c           >>> is a couple  c           >>> is a couple
3902              cltrx(ip,ntr)   = clx(nplanes-ip+1,icp_cp(id))              cltrx(ip,ntr)   = clx(nplanes-ip+1,icp_cp(id))
3903              cltry(ip,ntr)   = cly(nplanes-ip+1,icp_cp(id))              cltry(ip,ntr)   = cly(nplanes-ip+1,icp_cp(id))
Legend:
Removed from v.1.40  
changed lines
  Added in v.1.42
  ViewVC Help
Powered by ViewVC 1.1.23