--- DarthVader/TrackerLevel2/inc/F77/common_reduction.f 2006/09/05 12:52:20 1.3 +++ DarthVader/TrackerLevel2/inc/F77/common_reduction.f 2006/11/23 18:51:44 1.6 @@ -24,10 +24,11 @@ parameter (NSTRIPMIN=10)!(NSTRIPMIN=64)!(NSTRIPMIN=10)!TEMPORANEO real cn(nviews,nva1_view) ! common noise values + real cnrms(nviews,nva1_view) ! common noise rms integer cnn(nviews,nva1_view) ! >0 = number of strips used for cn computation c 0 = failed (less than NSTRIPMIN strips used) => temporarily masked -c -1 = masked a-priori - integer mask_vk_ev(nviews,nva1_view) ! mask of VA1s, local to the event (set to 0 if CN computation fails) +c -1 = masked a-priori +c integer mask_vk_ev(nviews,nva1_view) ! mask of VA1s, local to the event (set to 0 if CN computation fails) integer clstr(nviews,nva1_view,nstrips_va1) !flag matrix to mark ! strips with signal and exclude them from @@ -36,7 +37,8 @@ ! strips with unusually high or low signal and ! exclude them from common noise computation - common/cn/cn,cnn,mask_vk_ev,clstr,strange +c common/cn/cn,cnrms,cnn,mask_vk_ev,clstr,strange + common/cn/cn,cnrms,cnn,clstr,strange c------------------------------------------------------------------------ c @@ -53,11 +55,16 @@ c parameter (nclstrmax=1000) !maximum number of clusters per event + parameter (nclstrpmin = 5) !minimum number of stored strips + parameter (adc_satx = 2960) ! upper (x) saturation limit + parameter (adc_saty = 100) ! lower (y) saturation limit + integer ind real value(nstrips_view) !per trovare i cluster real clseedcut(nstrips_view) real clinclcut(nstrips_view) - common/searchcluster/value,clseedcut,clinclcut,ind + integer sat(nstrips_view) !saturation mask + common/searchcluster/value,clseedcut,clinclcut,sat,ind c added by Elena (08/2006) integer nclstr_view !number of clusters per view @@ -73,3 +80,7 @@ integer nshowers logical flag_shower common/shower/nshowers,flag_shower + + integer eventn_old(nviews) + integer good_old(nviews) + common/eventn/eventn_old,good_old