inc/F77/common_reduction.f



************************************************************************
*     
*     Common common_reduction.f
*     
*     contains variables used only in the execution of  reduction
*     
*     to be included in:
*     - reduction.f
*     
*     
*************************************************************************

      real adc(nviews,nva1_view,nstrips_va1)
      common/adc_value/adc

c------------------------------------------------------------------------
c     
c     common noise computation
c     
c------------------------------------------------------------------------
      parameter (scut=4.)!(scut=3.)
      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      
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
                                ! common noise computation
      integer strange(nviews,nva1_view,nstrips_va1) !flag matrix to mark 
                                ! strips with unusually high or low signal and
                                ! exclude them from common noise computation

c      common/cn/cn,cnrms,cnn,mask_vk_ev,clstr,strange
      common/cn/cn,cnrms,cnn,clstr,strange

c------------------------------------------------------------------------
c     
c     common noise inclusion cut
c     
c------------------------------------------------------------------------
      parameter (cnincut=2.)    !cut to include strips in cluster definition
                                ! during common noise computation
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)
      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
      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
1
2
3	************************************************************************
4	*
5	* Common common_reduction.f
6	*
7	* contains variables used only in the execution of reduction
8	*
9	* to be included in:
10	* - reduction.f
11	*
12	*
13	*************************************************************************
14
15	real adc(nviews,nva1_view,nstrips_va1)
16	common/adc_value/adc
17
18	c------------------------------------------------------------------------
19	c
20	c common noise computation
21	c
22	c------------------------------------------------------------------------
23	parameter (scut=4.)!(scut=3.)
24	parameter (NSTRIPMIN=10)!(NSTRIPMIN=64)!(NSTRIPMIN=10)!TEMPORANEO
25
26	real cn(nviews,nva1_view) ! common noise values
27	real cnrms(nviews,nva1_view) ! common noise rms
28	integer cnn(nviews,nva1_view) ! >0 = number of strips used for cn computation
29	c 0 = failed (less than NSTRIPMIN strips used) => temporarily masked
30	c -1 = masked a-priori
31	c integer mask_vk_ev(nviews,nva1_view) ! mask of VA1s, local to the event (set to 0 if CN computation fails)
32
33	integer clstr(nviews,nva1_view,nstrips_va1) !flag matrix to mark
34	! strips with signal and exclude them from
35	! common noise computation
36	integer strange(nviews,nva1_view,nstrips_va1) !flag matrix to mark
37	! strips with unusually high or low signal and
38	! exclude them from common noise computation
39
40	c common/cn/cn,cnrms,cnn,mask_vk_ev,clstr,strange
41	common/cn/cn,cnrms,cnn,clstr,strange
42
43	c------------------------------------------------------------------------
44	c
45	c common noise inclusion cut
46	c
47	c------------------------------------------------------------------------
48	parameter (cnincut=2.) !cut to include strips in cluster definition
49	! during common noise computation
50	c------------------------------------------------------------------------
51	c
52	c some common for reductionflight
53	c
54	c------------------------------------------------------------------------
55
56	c parameter (nclstrmax=1000) !maximum number of clusters per event
57
58	parameter (nclstrpmin = 5) !minimum number of stored strips
59	parameter (adc_satx = 2960) ! upper (x) saturation limit
60	parameter (adc_saty = 100) ! lower (y) saturation limit
61
62	integer ind
63	real value(nstrips_view) !per trovare i cluster
64	real clseedcut(nstrips_view)
65	real clinclcut(nstrips_view)
66	integer sat(nstrips_view) !saturation mask
67	common/searchcluster/value,clseedcut,clinclcut,sat,ind
68
69	c added by Elena (08/2006)
70	integer nclstr_view !number of clusters per view
71	integer ladder_view(nclstrmax) !ladder the strip carrying the largest signal
72	! value belongs to
73	integer maxs_view(nclstrmax) !cluster strip carrying the largest signal value
74	integer mult_view(nclstrmax) !cluster multiplicity
75	integer rmax_view(nclstrmax)
76	integer lmax_view(nclstrmax)
77	common/clustersperview/nclstr_view,ladder_view,maxs_view
78	$ ,mult_view,rmax_view,lmax_view
79
80	integer nshowers
81	logical flag_shower
82	common/shower/nshowers,flag_shower
83
84	integer eventn_old(nviews)
85	integer good_old(nviews)
86	common/eventn/eventn_old,good_old