--- DarthVader/TrackerLevel2/inc/F77/common_reduction.f	2006/08/04 08:18:06	1.2
+++ DarthVader/TrackerLevel2/inc/F77/common_reduction.f	2008/10/22 15:17:35	1.7
@@ -21,15 +21,14 @@
 c     
 c------------------------------------------------------------------------
       parameter (scut=4.)!(scut=3.)
-      parameter (NSTRIPMIN=64)!(NSTRIPMIN=10)!TEMPORANEO
-      parameter (NSTRIPWARNING=64)
+      parameter (NSTRIPMIN=10)!(NSTRIPMIN=64)!(NSTRIPMIN=10)!TEMPORANEO
 
       real    cn(nviews,nva1_view)      ! common noise values
-      integer cnflag(nviews,nva1_view)  ! common noise flags: 
-c                                         0 = ok 
-c                                        -1 = failed  (less than NSTRIPMIN strips used)
-c                                         1 = warning (less than NSTRIPWARNING strips used)
-      integer mask_vk_ev(nviews,nva1_view) ! mask of VA1s, local to the event (set to 0 if CN computation fails)
+      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      
+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
@@ -38,19 +37,8 @@
                                 ! strips with unusually high or low signal and
                                 ! exclude them from common noise computation
 
-      common/cn/cn,cnflag,mask_vk_ev,clstr,strange
-
-
-c$$$c------------------------------------------------------------------------
-c$$$c     
-c$$$c     cluster seed cuts
-c$$$c     
-c$$$c------------------------------------------------------------------------
-c$$$      parameter (clcutx=7.)     !cluster seed cut (to esclude particles in CN
-c$$$                                ! computation and to find clusters in data
-c$$$                                ! reduction)
-c$$$      parameter (clcuty=7.)
-
+c      common/cn/cn,cnrms,cnn,mask_vk_ev,clstr,strange
+      common/cn/cn,cnrms,cnn,clstr,strange
 
 c------------------------------------------------------------------------
 c     
@@ -59,43 +47,42 @@
 c------------------------------------------------------------------------
       parameter (cnincut=2.)    !cut to include strips in cluster definition
                                 ! during common noise computation
-
-
-c$$$c------------------------------------------------------------------------
-c$$$c     
-c$$$c     cluster finding inclusion cut
-c$$$c     
-c$$$c------------------------------------------------------------------------
-c$$$      parameter (incutx=4.)     !cut to include strips in cluster definition
-c$$$                                ! during cluster finding procedure
-c$$$      parameter (incuty=4.)
-
-c$$$c------------------------------------------------------------------------
-c$$$c
-c$$$c     variables passed from TRACKEREXE.cra
-c$$$c
-c$$$c------------------------------------------------------------------------
-c$$$      
-c$$$      integer trkpaw
-c$$$      integer procerr
-c$$$      character*20 data_file
-c$$$      character*50 data_file_level1
-c$$$      character*50 data_file_level2
-c$$$      character*74 data_file_calib 
-c$$$
-c$$$      common/c2red/trkpaw,procerr,data_file,data_file_level1
-c$$$     $     ,data_file_level2,data_file_calib
 c------------------------------------------------------------------------
 c
 c     some common for reductionflight
 c
 c------------------------------------------------------------------------
+
+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/calibratedsignal/value
+      save/calibratedsignal/
+      common/searchclustr/clseedcut,clinclcut,sat,ind 
+
+c     added by Elena (08/2006)
+      integer nclstr_view            !number of clusters per view
+      integer ladder_view(nclstrmax) !ladder the strip carrying the largest signal 
+                                     ! value belongs to
+      integer maxs_view(nclstrmax)   !cluster strip carrying the largest signal value
+      integer mult_view(nclstrmax)   !cluster multiplicity
+      integer rmax_view(nclstrmax)
+      integer lmax_view(nclstrmax)
+      common/clustersperview/nclstr_view,ladder_view,maxs_view
+     $     ,mult_view,rmax_view,lmax_view
 
       integer nshowers
       logical flag_shower
       common/shower/nshowers,flag_shower
+
+      integer eventn_old(nviews)
+      integer good_old(nviews)
+      common/eventn/eventn_old,good_old