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

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

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

revision 1.1 by mocchiut, Fri May 19 13:15:52 2006 UTC revision 1.9 by pam-fi, Fri Aug 17 14:36:05 2007 UTC
# Line 7  Line 7 
7  *      *    
8  *************************************************************************  *************************************************************************
9    
10          parameter (clcutx=7.)     !cluster seed cut
11          parameter (clcuty=6.)
12          parameter (incutx=4.)     !cut to include strips in cluster definition
13          parameter (incuty=4.)
14    
15  c------------------------------------------------------------------------  c------------------------------------------------------------------------
16  c      c    
17  c     pedestal, sigma and badstrip variables  c     some parameters to configure data reduction
18  c      c    
19  c------------------------------------------------------------------------  c------------------------------------------------------------------------
20    c$$$      real clcutx               !cluster seed cut
21    c$$$      real clcuty
22    c$$$      real incutx               !cut to include strips in cluster definition
23    c$$$      real incuty
24          integer pfaid             !if of PFA
25          common/sw/
26         $     pfaid
27    c$$$     $     ,clcutx,clcuty
28    c$$$     $     ,incutx,incuty
29    
30        parameter (id_hi_bad=100) !badstrip histo identifier  c------------------------------------------------------------------------
31        parameter (id_hi_ped=200) !pedestal histo identifier  c    
32        parameter (id_hi_sig=300) !sigma histo identifier  c     pedestal, sigma and badstrip variables
33    c    
34    c------------------------------------------------------------------------
35    
36        real pedestal(nviews,nva1_view,nstrips_va1) !pedestal value        real pedestal(nviews,nva1_view,nstrips_va1) !pedestal value
37        real pedestal_t(nviews,nva1_view,nstrips_va1) !pedestal truncated value        real pedestal_t(nviews,nva1_view,nstrips_va1) !pedestal truncated value
# Line 32  c--------------------------------------- Line 48  c---------------------------------------
48        save/pedsigbad/        save/pedsigbad/
49    
50    
       parameter (clcutx=7.)     !cluster seed cut  
       parameter (clcuty=7.)  
       parameter (incutx=4.)     !cut to include strips in cluster definition  
       parameter (incuty=4.)  
51    
52    c------------------------------------------------------------------------
53    c    
54    c     VA1-mask
55    c    
56    c------------------------------------------------------------------------
57    c     VA1 mask, from the DB
58    c      0 = force mask
59    c      1 = ---> run-by-run mask ---> event-by-event mask
60    c     -1 = ---> event-by-event mask
61    c     -------------------------------------
62          integer mask_vk(nviews,nva1_view)    
63    c     -------------------------------------
64    c     VA1 mask run-by-run (set to 0 if <SIG> < SIG_min)
65    c      0 = force mask
66    c      1 = ---> event-by-event mask
67    c     -------------------------------------
68          integer mask_vk_run(nviews,nva1_view)
69    c     -------------------------------------
70    c     VA1 mask event-by-event (set to 0 if CN computation fails)
71    c      0 = force mask
72    c      1 = search clusters
73    c     -------------------------------------
74          integer mask_vk_ev(nviews,nva1_view)
75    
76        parameter (id_hi_mask=400) !mask histo identifier        integer mask(nviews,nva1_view,nstrips_va1) !effective mask
   
       integer mask(nviews,nva1_view,nstrips_va1)  
       integer mask_vk(nviews,nva1_view)  
77    
78        common/mask/mask,mask_vk  *     set from outside F77 routines
79          common/mask/mask_vk,mask_vk_run
80        save/mask/        save/mask/
81    *     set event-by-event
82          common/maskev/mask_vk_ev,mask
83    
84  c------------------------------------------------------------------------  c------------------------------------------------------------------------
85  c      c    
# Line 53  c     Line 88  c    
88  c------------------------------------------------------------------------  c------------------------------------------------------------------------
89  *     angular binning  *     angular binning
90  *     NB the angular binning is defined in common_preanalysis.f  *     NB the angular binning is defined in common_preanalysis.f
91  *     here are defined variables to retrieve calibration info  *     here are defned variables to retrieve calibration info
92  *     (for processing levels next to 1)  *     (for processing levels next to 1)
93        parameter (nangmax=21)    !maximum number of angular bins        parameter (nangmax=21)    !maximum number of angular bins
94        integer nangbin        integer nangbin
95        real angL(nangmax),angR(nangmax) !Left and Right bin limits        real angL(nangmax),angR(nangmax) !Left and Right bin limits
96    
97        parameter (netavalmax=150)   !eta2 points        parameter (netavalmax=500)   !eta2 points
98        integer netaval        integer netaval
99        real eta2(netavalmax,nangmax)        real eta2(netavalmax,nangmax)
100        real feta2(netavalmax,nviews,nladders_view,nangmax)        real feta2(netavalmax,nviews,nladders_view,nangmax)
# Line 67  c--------------------------------------- Line 102  c---------------------------------------
102        real feta3(netavalmax,nviews,nladders_view,nangmax)        real feta3(netavalmax,nviews,nladders_view,nangmax)
103        real eta4(netavalmax,nangmax)        real eta4(netavalmax,nangmax)
104        real feta4(netavalmax,nviews,nladders_view,nangmax)        real feta4(netavalmax,nviews,nladders_view,nangmax)
105          real fcorr(nviews,nladders_view,nangmax)
106    
107        common/angbinning/nangbin,angL,angR  c$$$      common/angbinning/nangbin,angL,angR
108        save/angbinning/  c$$$      save/angbinning/
109        common/pfa/netaval,eta2,feta2,eta3,feta3,eta4,feta4  c$$$      common/pfa/netaval,eta2,feta2,eta3,feta3,eta4,feta4
110    c$$$      save/pfa/
111        
112          common/pfa/nangbin,angL,angR
113         $     ,netaval,eta2,feta2,eta3,feta3,eta4,feta4
114         $     ,fcorr
115        save/pfa/        save/pfa/
116    
117    *     limit of application of pfa
118          parameter (e2fax=0.)
119          parameter (e2tax=10.)
120          parameter (e3fax=10.)
121          parameter (e3tax=15.)
122          parameter (e4fax=90.)
123          parameter (e4tax=90.)
124    
125          parameter (e2fay=0.)
126          parameter (e2tay=20.)
127          parameter (e3fay=90.)
128          parameter (e3tay=90.)
129          parameter (e4fay=90.)
130          parameter (e4tay=90.)
131    
132  c------------------------------------------------------------------------  c------------------------------------------------------------------------
133  c      c    
134  c     CHARGE CORRELATION  c     CHARGE CORRELATION
# Line 79  c     Line 136  c    
136  c------------------------------------------------------------------------  c------------------------------------------------------------------------
137    
138  c      parameter (chcut=3.)      !cut to associate two clusters  c      parameter (chcut=3.)      !cut to associate two clusters
139        parameter (chcut=4.)      !cut to associate two clusters  c      parameter (chcut=4.)      !cut to associate two clusters
140          parameter (chcut=50.)      !cut to associate two clusters
141          parameter (chsatx=1700.)   !saturation limit
142          parameter (chsaty=2500.)   !saturation limit
143          parameter (chmipx=200.)   !mip limit
144          parameter (chmipy=200.)   !mip limit
145    *     above saturation limit correlation conditions are not applied
146        real kch(nplanes,nladders_view) !angular coeff        real kch(nplanes,nladders_view) !angular coeff
147        real cch(nplanes,nladders_view) !const        real cch(nplanes,nladders_view) !const
148        real sch(nplanes,nladders_view) !sigma        real sch(nplanes,nladders_view) !sigma
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.23