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

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

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

revision 1.27 by pam-fi, Fri Dec 5 08:26:48 2008 UTC revision 1.29 by mocchiut, Thu Jan 16 15:29:58 2014 UTC
# Line 266  C--------------------------------------- Line 266  C---------------------------------------
266  C     come here if GOOD1=0  C     come here if GOOD1=0
267  C     or the event has too many clusters  C     or the event has too many clusters
268  C---------------------------------------------  C---------------------------------------------
269   200  continue  c 200  continue
270    
271        ngood = 0        ngood = 0
272        do iv = 1,nviews        do iv = 1,nviews
# Line 323  c      good1 = 0 Line 323  c      good1 = 0
323        do id=1,maxlength         !???        do id=1,maxlength         !???
324           clsignal(id) = 0.           clsignal(id) = 0.
325           clsigma(id)  = 0.           clsigma(id)  = 0.
326           cladc(id)    = 0.           cladc(id)    = 0  ! EM GCC4.7 is integer
327           clbad(id)    = 0.           clbad(id)    = 0  ! EM GCC4.7 is integer
328        enddo        enddo
329        do iv=1,nviews        do iv=1,nviews
330  c        crc1(iv)=0  c        crc1(iv)=0
# Line 583  c--------------------------------------- Line 583  c---------------------------------------
583       $                    delta.gt.cut .and.       $                    delta.gt.cut .and.
584       $                    value(iss).gt.clseedcut(iss).and.       $                    value(iss).gt.clseedcut(iss).and.
585       $                    .true.)then       $                    .true.)then
586                          if( value(iss).gt.vmax )then                                                  if( value(iss).gt.vmax )then
587                             imax = iss                             imax = iss
588                             vmax = value(iss)                             vmax = value(iss)
589                          else                          else
# Line 658  c$$$            print*,ind,clsignal(ind) Line 658  c$$$            print*,ind,clsignal(ind)
658              ist=nst(j)              ist=nst(j)
659    
660              clsigma(ind) = sigma(iv,ivk,ist)              clsigma(ind) = sigma(iv,ivk,ist)
661              cladc(ind)   = adc(iv,ivk,ist)              cladc(ind)   = int(adc(iv,ivk,ist))   ! EM GCC4.7 is integer
662              clbad(ind)   = bad(iv,ivk,ist)              clbad(ind)   = bad(iv,ivk,ist)
663  c            clped(ind)   = pedestal(iv,ivk,ist)  c            clped(ind)   = pedestal(iv,ivk,ist)
664    
# Line 740  C===  > Y view Line 740  C===  > Y view
740  c     print*,iv,nvk(is),nst(is),adc(iv,nvk(is),nst(is))  c     print*,iv,nvk(is),nst(is),adc(iv,nvk(is),nst(is))
741  c     $            ,cn(iv,nvk(is))  c     $            ,cn(iv,nvk(is))
742  c     $            ,pedestal(iv,nvk(is),nst(is))  c     $            ,pedestal(iv,nvk(is),nst(is))
743              value(is)= -(DBLE(adc(iv,nvk(is),nst(is)))              value(is)= -(REAL(adc(iv,nvk(is),nst(is)))  ! EM GCC4.7 value(nstrips_view) is real not double
744       $           -cn(iv,nvk(is))-pedestal(iv,nvk(is),nst(is)))       $           -cn(iv,nvk(is))-pedestal(iv,nvk(is),nst(is)))
745       $           *mask(iv,nvk(is),nst(is))       $           *mask(iv,nvk(is),nst(is))
746              clseedcut(is)=clcuty*sigma(iv,nvk(is),nst(is))              clseedcut(is)=clcuty*sigma(iv,nvk(is),nst(is))
# Line 751  c     $            ,pedestal(iv,nvk(is), Line 751  c     $            ,pedestal(iv,nvk(is),
751              if( adc(iv,nvk(is),nst(is)).lt.adc_saty )sat(is)=1              if( adc(iv,nvk(is),nst(is)).lt.adc_saty )sat(is)=1
752           else                       else            
753  C===  > X view  C===  > X view
754              value(is)= (DBLE(adc(iv,nvk(is),nst(is)))              value(is)= (REAL(adc(iv,nvk(is),nst(is)))
755       $           -cn(iv,nvk(is))-pedestal(iv,nvk(is),nst(is)))       $           -cn(iv,nvk(is))-pedestal(iv,nvk(is),nst(is)))
756       $           *mask(iv,nvk(is),nst(is))       $           *mask(iv,nvk(is),nst(is))
757              clseedcut(is)=clcutx*sigma(iv,nvk(is),nst(is))              clseedcut(is)=clcutx*sigma(iv,nvk(is),nst(is))
Legend:
Removed from v.1.27  
changed lines
  Added in v.1.29
  ViewVC Help
Powered by ViewVC 1.1.23