/[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.41 by mocchiut, Tue Aug 11 12:58:24 2009 UTC revision 1.44 by pamelats, Wed Dec 30 14:22:23 2009 UTC
# Line 528  c         rchi2=chi2/dble(ndof) Line 528  c         rchi2=chi2/dble(ndof)
528    
529    
530  *     --- and store the results --------------------------------  *     --- and store the results --------------------------------
          ntrk = ntrk + 1                   !counter of found tracks  
          if(.not.FIMAGE  
      $        .and.iimage.eq.0) image(ntrk)= 0  
          if(.not.FIMAGE  
      $        .and.iimage.ne.0)image(ntrk)=ntrk+1 !this is the image of the next  
          if(FIMAGE)     image(ntrk)=ntrk-1 !this is the image of the previous  
          call fill_level2_tracks(ntrk)     !==> good2=.true.  
   
531           if(ntrk.eq.NTRKMAX)then           if(ntrk.eq.NTRKMAX)then
532              if(verbose.eq.1)              if(verbose.eq.1)
533       $           print*,       $           print*,
# Line 545  c         rchi2=chi2/dble(ndof) Line 537  c         rchi2=chi2/dble(ndof)
537  cc            good2=.false.  cc            good2=.false.
538              goto 880            !fill ntp and go to next event              goto 880            !fill ntp and go to next event
539           endif           endif
540    
541             ntrk = ntrk + 1                   !counter of found tracks
542             if(.not.FIMAGE
543         $        .and.iimage.eq.0) image(ntrk)= 0
544             if(.not.FIMAGE
545         $        .and.iimage.ne.0)image(ntrk)=ntrk+1 !this is the image of the next
546             if(FIMAGE)     image(ntrk)=ntrk-1 !this is the image of the previous
547             call fill_level2_tracks(ntrk)     !==> good2=.true.
548    
549    c$$$         if(ntrk.eq.NTRKMAX)then
550    c$$$            if(verbose.eq.1)
551    c$$$     $           print*,
552    c$$$     $           '** warning ** number of identified '//
553    c$$$     $           'tracks exceeds vector dimension '
554    c$$$     $           ,'( ',NTRKMAX,' )'
555    c$$$cc            good2=.false.
556    c$$$            goto 880            !fill ntp and go to next event
557    c$$$         endif
558           if(iimage.ne.0)then           if(iimage.ne.0)then
559              FIMAGE=.true.       !              FIMAGE=.true.       !
560              goto 1212           !>>> fit image-track              goto 1212           !>>> fit image-track
# Line 3884  c$$$         angy    = 180.*atan(tgtemp) Line 3894  c$$$         angy    = 180.*atan(tgtemp)
3894           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
3895           LS(IP,ntr)      = ssensor+10*sladder           LS(IP,ntr)      = ssensor+10*sladder
3896    
3897           if(id.ne.0)then  c         if(id.ne.0)then
3898    CCCCCC(10 novembre 2009) PATCH X NUCLEI
3899    C     non ho capito perche', ma durante il ritracciamento dei nuclei
3900    C     (quando una traccia viene trovata ma non e' buona) c'e' qualche variabile
3901    C     che non viene reinizializzata correttamente e i cluster esclusi
3902    C     dal fit risultano ancora inclusi...
3903    C
3904             cltrx(ip,ntr)   = 0
3905             cltry(ip,ntr)   = 0
3906    c$$$         if(
3907    c$$$     $        xgood_nt(ip,ntr).eq.1.and.ygood_nt(ip,ntr).eq.1
3908    c$$$     $        .and.
3909    c$$$     $        id.ne.0)then
3910             if(id.ne.0)then        !patch 30/12/09
3911    
3912  c           >>> is a couple  c           >>> is a couple
3913              cltrx(ip,ntr)   = clx(nplanes-ip+1,icp_cp(id))              cltrx(ip,ntr)   = clx(nplanes-ip+1,icp_cp(id))
3914              cltry(ip,ntr)   = cly(nplanes-ip+1,icp_cp(id))              cltry(ip,ntr)   = cly(nplanes-ip+1,icp_cp(id))
# Line 3923  c           >>> is a couple Line 3947  c           >>> is a couple
3947                 ypu(ip,ntr)      = corr                 ypu(ip,ntr)      = corr
3948              endif              endif
3949    
3950           elseif(icl.ne.0)then  c$$$         elseif(icl.ne.0)then
3951             endif                  !patch 30/12/09
3952             if(icl.ne.0)then       !patch 30/12/09
3953    
3954              cl_used(icl) = 1    !tag used clusters                        cl_used(icl) = 1    !tag used clusters          
3955    
Legend:
Removed from v.1.41  
changed lines
  Added in v.1.44
  ViewVC Help
Powered by ViewVC 1.1.23