--- DarthVader/TrackerLevel2/src/F77/cncomp.f 2006/10/26 16:22:37 1.6 +++ DarthVader/TrackerLevel2/src/F77/cncomp.f 2007/08/31 14:56:52 1.9 @@ -33,6 +33,23 @@ strange(i,j,k)=1 !initializes unusually high or low signal enddo ! affected strips flag +c------------------------------------------------------------------------ +c (september 2007) +c remove from CN computation the first and the last 3 channels of +c each X view, becouse they ar not connected to any strip +c------------------------------------------------------------------------ + if(mod(i,2).eq.0)then + if(j.eq.1)then + do k=1,3 + strange(i,j,k)=0 + enddo + elseif(j.eq.nva1_ladder)then + do k=nstrips_va1,nstrips_va1-2,-1 + strange(i,j,k)=0 + enddo + endif + endif + newclstr=1 !flag to warn about new found signal c ! affected strips c------------------------------------------------------------------------ @@ -171,11 +188,11 @@ if(ncn.lt.NSTRIPMIN) then !no signal free strips on this VA1... if(ncn.eq.0)then - if(debug)print*,' WARNING - cnoise: ', + if(debug.eq.1)print*,' WARNING - cnoise: ', $ 'no strips for CN computation on VA1 ',j, $ ', VIEW ',i,' >>> FAILED ' else - if(debug)print*,' WARNING - cnoise: ', + if(debug.eq.1)print*,' WARNING - cnoise: ', $ 'less than ',NSTRIPMIN $ ,' strips for CN computation on VA1 ',j, $ ', VIEW ',i,' >>> FAILED ' @@ -187,7 +204,10 @@ cnrms(i,j)= SQRT( cnrms(i,j)/DBLE(ncn) - cn(i,j)**2 ) cnn(i,j) = ncn gulp=0 - if(debug.and.ABS(cn(i,j)).gt.1000) +c$$$ print*,'Event ',eventn(1) +c$$$ $ ,': cn(',i,',',j,')= ',cn(i,j),' ncn ',ncn + + if(debug.eq.1.and.ABS(cn(i,j)).gt.1000) $ print*,'Event ',eventn(1) $ ,': cn(',i,',',j,')= ',cn(i,j),' ncn ',ncn endif