/[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.28 by mocchiut, Tue Aug 4 14:01:39 2009 UTC revision 1.29 by mocchiut, Thu Jan 16 15:29:58 2014 UTC
# 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 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.28  
changed lines
  Added in v.1.29
  ViewVC Help
Powered by ViewVC 1.1.23