/[PAMELA software]/DarthVader/TrackerLevel2/inc/F77/common_reduction.f
ViewVC logotype

Diff of /DarthVader/TrackerLevel2/inc/F77/common_reduction.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.2 by pam-fi, Fri Aug 4 08:18:06 2006 UTC revision 1.5 by pam-fi, Thu Oct 26 16:22:37 2006 UTC
# Line 21  c     common noise computation Line 21  c     common noise computation
21  c      c    
22  c------------------------------------------------------------------------  c------------------------------------------------------------------------
23        parameter (scut=4.)!(scut=3.)        parameter (scut=4.)!(scut=3.)
24        parameter (NSTRIPMIN=64)!(NSTRIPMIN=10)!TEMPORANEO        parameter (NSTRIPMIN=10)!(NSTRIPMIN=64)!(NSTRIPMIN=10)!TEMPORANEO
       parameter (NSTRIPWARNING=64)  
25    
26        real    cn(nviews,nva1_view)      ! common noise values        real    cn(nviews,nva1_view)      ! common noise values
27        integer cnflag(nviews,nva1_view)  ! common noise flags:        real    cnrms(nviews,nva1_view)   ! common noise rms
28  c                                         0 = ok        integer cnn(nviews,nva1_view)     ! >0 = number of strips used for cn computation
29  c                                        -1 = failed  (less than NSTRIPMIN strips used)  c                                          0 = failed  (less than NSTRIPMIN strips used) => temporarily masked
30  c                                         1 = warning (less than NSTRIPWARNING strips used)  c                                         -1 = masked a-priori      
31        integer mask_vk_ev(nviews,nva1_view) ! mask of VA1s, local to the event (set to 0 if CN computation fails)        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        integer clstr(nviews,nva1_view,nstrips_va1) !flag matrix to mark
# Line 38  c                                       Line 37  c                                      
37                                  ! strips with unusually high or low signal and                                  ! strips with unusually high or low signal and
38                                  ! exclude them from common noise computation                                  ! exclude them from common noise computation
39    
40        common/cn/cn,cnflag,mask_vk_ev,clstr,strange        common/cn/cn,cnrms,cnn,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.)  
   
41    
42  c------------------------------------------------------------------------  c------------------------------------------------------------------------
43  c      c    
# Line 59  c     Line 46  c    
46  c------------------------------------------------------------------------  c------------------------------------------------------------------------
47        parameter (cnincut=2.)    !cut to include strips in cluster definition        parameter (cnincut=2.)    !cut to include strips in cluster definition
48                                  ! during common noise computation                                  ! 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  
49  c------------------------------------------------------------------------  c------------------------------------------------------------------------
50  c  c
51  c     some common for reductionflight  c     some common for reductionflight
52  c  c
53  c------------------------------------------------------------------------  c------------------------------------------------------------------------
54    
55    c     parameter (nclstrmax=1000) !maximum number of clusters per event
56    
57          parameter (nclstrpmin = 5) !minimum number of stored strips
58          parameter (adc_satx = 2960) ! upper (x) saturation limit
59          parameter (adc_saty = 100)  ! lower (y) saturation limit
60    
61        integer ind        integer ind
62        real value(nstrips_view)  !per trovare i cluster        real value(nstrips_view)  !per trovare i cluster
63        real clseedcut(nstrips_view)        real clseedcut(nstrips_view)
64        real clinclcut(nstrips_view)        real clinclcut(nstrips_view)
65        common/searchcluster/value,clseedcut,clinclcut,ind        integer sat(nstrips_view) !saturation mask
66          common/searchcluster/value,clseedcut,clinclcut,sat,ind
67    
68    c     added by Elena (08/2006)
69          integer nclstr_view            !number of clusters per view
70          integer ladder_view(nclstrmax) !ladder the strip carrying the largest signal
71                                         ! value belongs to
72          integer maxs_view(nclstrmax)   !cluster strip carrying the largest signal value
73          integer mult_view(nclstrmax)   !cluster multiplicity
74          integer rmax_view(nclstrmax)
75          integer lmax_view(nclstrmax)
76          common/clustersperview/nclstr_view,ladder_view,maxs_view
77         $     ,mult_view,rmax_view,lmax_view
78    
79        integer nshowers        integer nshowers
80        logical flag_shower        logical flag_shower
81        common/shower/nshowers,flag_shower        common/shower/nshowers,flag_shower
82    
83          integer eventn_old(nviews)
84          integer good_old(nviews)
85          common/eventn/eventn_old,good_old
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.23